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 22:25] gag [SIF contents] |
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]]. |
< | < | ||
- | Solver 1 | + | ! Units : MPa - m - yr |
+ | $yearinsec = 365.25*24*60*60 | ||
+ | |||
+ | !Compute the heat generated by ice deformation | ||
+ | Solver | ||
+ | Equation = DeformationalHeat | ||
+ | Variable = W | ||
+ | Variable DOFs = 1 | ||
+ | |||
+ | procedure = " | ||
+ | |||
+ | Linear System Solver = direct | ||
+ | Linear System direct Method = umfpack | ||
+ | End | ||
+ | |||
+ | 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 58: | Line 80: | ||
Body Force 1 | Body Force 1 | ||
... | ... | ||
- | Temp Volume Source = Real 1.0E00 | + | Temp Volume Source = Equals W ! The volumetric heat source |
End | End | ||
Line 66: | 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 81: | Line 103: | ||
End | End | ||
+ | !Upper surface | ||
Boundary Condition 1 | Boundary Condition 1 | ||
... | ... | ||
- | | + | Temp = Real 263.15 |
- | !--------------------------------------------------------- | + | |
- | | + | |
- | Real MATC " | + | |
End | End | ||
+ | !Bedrock | ||
Boundary Condition 2 | Boundary Condition 2 | ||
... | ... | ||
Line 95: | 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. |