I was doing some research in the temperature of the antarctic ice sheet,I add the temperature as the website said
Code: Select all
http://elmerice.elmerfem.org/wiki/doku.php?id=solvers:temperature
Every thing goes well ,but when I want to obtain the temeprature,the problem occurs.the temperature of the bottom is 274 but the pressure melting point I set was
In my opinion ,the temperature must be under the pressure melting point,the results I obtain as follows:!! pressuremeltingpoint (in SI units)
$ function pressuremeltingpoint2(D) {\
P = 910*D*9.81;\
if (P<0.0) P=0.0;\
beta=9.8E-08;\
_pressuremeltingpoint2=273.15-(beta*P);\
}
[img]
[/img]
By the way ,I also have some questions about the parameter of the glen's law .in my sif file I set the parameters as follows:
I do not know where I can obtain these parameters or what parameters should I set .! Prefactor from Cuffey and Paterson (2010) in MPa^{-3} a^{-1}
$A1 = 1.258e13
$A2 = 6.046e28
$Q1 = 78.80e3
$Q2 = 139.0e3
I am a newer,any help is welcome
looking forward to everyone of you