This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:temperature [2012/11/14 23:05] gag [General Informations] |
solvers:temperature [2018/10/28 14:28] tzwinger [Examples] |
||
---|---|---|---|
Line 4: | Line 4: | ||
==== General Informations ==== | ==== General Informations ==== | ||
- | * **Solver Fortran File: | + | * **Solver Fortran File: |
- | * **Solver Name: | + | * **Solver Name: |
- | * **Required Output Variable(s): | + | * **Required Output Variable(s): |
- | * **Required Input Variable(s): | + | * **Required Input Variable(s): |
- | * **Optional Output Variable(s): | + | * **Optional Output Variable(s): |
- | * **Optional Input Variable(s): | + | * **Optional Input Variable(s): |
==== General Description ==== | ==== General Description ==== | ||
- | This solver treats the heat transfer problem with respect to an upper limit of the temperature (usually with ice the pressure-melting point, T< T_pm). Optionally, such a limit (and furthermore also a lower limit, e.g., T > 0 K) is introduced by solving the consequent | + | This solver treats the heat transfer problem with respect to an upper limit of the temperature (usually with ice the pressure-melting point, T< T_pm). Optionally, such a limit (and furthermore also a lower limit, e.g., T > 0 K) is introduced by solving the consequent |
+ | |||
+ | 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 20: | 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 42: | Line 38: | ||
Variable DOFs = 1 | Variable DOFs = 1 | ||
- | procedure = "./ | + | procedure = "ElmerIceSolvers" " |
Linear System Solver = direct | Linear System Solver = direct | ||
Line 50: | Line 46: | ||
Solver 4 | Solver 4 | ||
Equation = String " | Equation = String " | ||
- | Procedure = File "TemperateIce" " | + | Procedure = File "ElmerIceSolvers" " |
! Comment next line in parallel, as EliminateDirichlet does | ! Comment next line in parallel, as EliminateDirichlet does | ||
! not work in parallel | ! not work in parallel | ||
Line 92: | 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 110: | Line 106: | ||
Boundary Condition 1 | Boundary Condition 1 | ||
... | ... | ||
- | Temp = Real -10.0 | + | Temp = Real 263.15 |
End | End | ||
Line 120: | 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. |