This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:temperature [2013/02/28 23:29] gag [SIF contents] |
solvers:temperature [2018/10/28 14:28] tzwinger [Examples] |
||
---|---|---|---|
Line 16: | Line 16: | ||
The volumetric heat source term can be estimated from the ice flow deformational heat using the [[: | The volumetric heat source term can be estimated from the ice flow deformational heat using the [[: | ||
+ | |||
+ | **Looping option, added 6th June 2013 (version 6205)**: It is possible in some cases that the nonlinear convergence tolerance can be reached before all nodes have been properly constrained. | ||
+ | < | ||
+ | Loop While Unconstrained Nodes = Logical True | ||
+ | </ | ||
==== Known bugs ==== | ==== Known bugs ==== | ||
Line 22: | Line 27: | ||
==== SIF contents ==== | ==== SIF contents ==== | ||
- | The required keywords in the SIF file for this solver are: | + | The required keywords in the SIF file for this solver are given below. The MATC functions used here are explain and given on this [[tips:thermoprop|page]]. |
< | < | ||
- | ! Some useful MATC functions | + | ! Units : MPa - m - yr |
- | !! conductivity | + | $yearinsec |
- | $ function conductivity(T) | + | |
- | !! capacity | + | |
- | $ function capacity(T) { _capacity=146.3+(7.253*T)} | + | |
- | + | ||
- | !! pressuremeltingpoint | + | |
- | $ function pressuremeltingpoint(PIN) {\ | + | |
- | P = PIN;\ | + | |
- | if (P<0.0) P=0.0;\ | + | |
- | beta=9.8E-08*1.0E06;\ | + | |
- | _pressuremeltingpoint=273.15-(beta*P);\ | + | |
- | } | + | |
!Compute the heat generated by ice deformation | !Compute the heat generated by ice deformation | ||
Line 94: | Line 88: | ||
| | ||
Temp Heat Capacity = Variable Temp | Temp Heat Capacity = Variable Temp | ||
- | Real MATC " | + | Real MATC " |
! the heat conductivity as a MATC function of temperature itself | ! the heat conductivity as a MATC function of temperature itself | ||
| | ||
Temp Heat Conductivity = Variable Temp | Temp Heat Conductivity = Variable Temp | ||
- | Real MATC " | + | Real MATC " |
! Upper limit - pressure melting point | ! Upper limit - pressure melting point | ||
| | ||
Line 122: | Line 116: | ||
!-------------------- | !-------------------- | ||
Temp Flux BC = Logical True | Temp Flux BC = Logical True | ||
- | Temp Heat Flux = Real 56.05E-03 | + | Temp Heat Flux = Real $56.05E-03*yearinsec*1.0E-6 |
+ | !------------------- | ||
+ | ! frictional heat | ||
+ | !-------------------- | ||
+ | Temp Load = Variable Velocity 1 | ||
+ | Real Procedure | ||
End | End | ||
</ | </ | ||
+ | |||
+ | See also [[: | ||
==== Examples ==== | ==== Examples ==== | ||
- | Download an example | + | An example |
==== Reference ==== | ==== Reference ==== | ||
When used this solver can be cited using the following reference: | When used this solver can be cited using the following reference: | ||
Zwinger T. , R. Greve, O. Gagliardini , T. Shiraiwa and M. Lyly, 2007. A full Stokes-flow thermo-mechanical model for firn and ice applied to the Gorshkov crater glacier, Kamchatka. Annals of Glaciol., 45, p. 29-37. | Zwinger T. , R. Greve, O. Gagliardini , T. Shiraiwa and M. Lyly, 2007. A full Stokes-flow thermo-mechanical model for firn and ice applied to the Gorshkov crater glacier, Kamchatka. Annals of Glaciol., 45, p. 29-37. |