This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:enthalpy [2014/06/18 12:20] gag [General Description] |
solvers:enthalpy [2017/04/21 11:28] (current) tzwinger [Examples] |
||
---|---|---|---|
Line 15: | Line 15: | ||
Solves the enthalpy equation: | Solves the enthalpy equation: | ||
- | < | + | <m> |
+ | rho {{\partial H}/ | ||
+ | </m> | ||
- | where | + | where |
* < | * < | ||
Line 26: | Line 28: | ||
* <m>tr (sigma epsilon)</ | * <m>tr (sigma epsilon)</ | ||
* < | * < | ||
+ | |||
+ | Enthalpy is defined as a function of the water content < | ||
+ | |||
+ | <m> | ||
+ | If H < H_f, then H(T, omega) = int_T0^T C_p (T) dT | ||
+ | </m> | ||
+ | |||
+ | <m> | ||
+ | If H > H_f, then H(T, omega) = int_T0^Tm C_p (T) dT + omega L | ||
+ | </m> | ||
+ | |||
+ | where | ||
+ | |||
+ | * < | ||
+ | * < | ||
+ | * < | ||
+ | |||
+ | |||
+ | For the boundary conditions, a flux ('' | ||
+ | |||
+ | |||
==== SIF contents ==== | ==== SIF contents ==== | ||
+ | In this example, ice velocity are in m/s and pressure en MPa. | ||
< | < | ||
Solver 2 | Solver 2 | ||
Line 48: | Line 72: | ||
Stabilize = True | Stabilize = True | ||
- | Exported Variable 1 = String "Phase Change Enthalpy" | + | Exported Variable 1 = String "Phase Change Enthalpy" |
Exported Variable 1 DOFs = 1 | Exported Variable 1 DOFs = 1 | ||
- | Exported Variable 2 = String "Water Content" | + | Exported Variable 2 = String "Water Content" |
Exported Variable 2 DOFs = 1 | Exported Variable 2 DOFs = 1 | ||
- | Exported Variable 3 = String " | + | Exported Variable 3 = String " |
Exported Variable 3 DOFs = 1 | Exported Variable 3 DOFs = 1 | ||
End | End | ||
Constants | Constants | ||
- | | + | |
- | | + | |
+ | ! Cp(T) = A*T + B | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
End | End | ||
Line 68: | Line 98: | ||
Material 1 | Material 1 | ||
- | Enthalpy Density = real 917.0 | + | Enthalpy Density = real 917.0 !(kg m-3) |
- | Enthalpy Heat Diffusivity = Real $2.1/2050.0 ! = k / Cp | + | Enthalpy Heat Diffusivity = Real $2.1/2050.0 ! = k / Cp (kg m-1 s-1) |
- | Enthalpy Water Diffusivity = real 1.045e-4! | + | Enthalpy Water Diffusivity = real 1.045e-4 ! (kg m-1 s-1) |
End | End | ||
Line 81: | Line 111: | ||
Enthalpy Heat Flux BC = logical True | Enthalpy Heat Flux BC = logical True | ||
- | Enthalpy Heat Flux = real 0.02 | + | Enthalpy Heat Flux = real 0.02 !(W m-2) |
End | End | ||
Line 88: | Line 118: | ||
Target Boundaries = 2 | Target Boundaries = 2 | ||
Enthalpy_h = variable coordinate 3 | Enthalpy_h = variable coordinate 3 | ||
- | real MATC " | + | real MATC " |
+ | End | ||
End | End | ||
</ | </ | ||
Line 94: | Line 125: | ||
==== Examples ==== | ==== Examples ==== | ||
- | An example solving for the enthalpy within the Tete Rousse glacier assuming an elevation | + | An example solving for the enthalpy within the Tete Rousse glacier assuming an elevation |
==== References ==== | ==== References ==== | ||
+ | Gilbert, A., O. Gagliardini, |