The Friction Changing with Temperature

aisundong

I found a bottom friction parameterization scheme in the model：
md.friction=frictiontemp();
It's roughly described like this:
disp(sprintf('Basal shear stress parameters: tau_b = coefficient² * Neff ^r * |u_b|^s-1 * u_b * 1/f(T)\n(effective stress Neff=rho_icegthickness+rho_watergbed, r=q/p and s=1/p)'));
Does this mean that using this method can consider the impact of temperature changes on bottom friction?Because the bottom melting caused by temperature raising sometimes has a significant impact on the friction coefficient, I want to consider this process in the simulation.

mathieumorlighem

Hello aisundong

yes, this is one possible parameterization of basal friction that accounts for basal temperature. It scales the friction by 1/f(T) where f(T) = exp((T-Tpmp)/gamma) (see source code here: src/c/classes/Loads/Friction.cpp function Friction::GetAlpha2Temp).
Cheers
Mathieu

aisundong

mathieumorlighem
Thank you very much for your answer!
I would like to ask if it is possible to provide a reference for determining f (T), or how did you conclude that f (T) should be an exponential expression like this?

mathieumorlighem

Hi @aisundong

this expression is definitely made up. If you would like to try out another relationship, feel free to change the code and test it out. Let us know if you need help
Mathieu