Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Last revision Both sides next revision
solvers:enthalpy [2014/06/20 13:07]
gag [References]
solvers:enthalpy [2015/04/13 07:33]
gag [References]
Line 42: Line 42:
  
   * <m>H_f</m> is the enthalpy of fusion, defined from the fusion temperature according to the pressure dependent Clausius-Clapeyron relationship.     * <m>H_f</m> is the enthalpy of fusion, defined from the fusion temperature according to the pressure dependent Clausius-Clapeyron relationship.  
-  * <m>C_p</m> is the temperature dependant heat capacity+  * <m>C_p</m> is the temperature dependant heat capacity, defined as <m>C_p = AT+B</m>
   * <m>L</m> is the latent heat of fusion   * <m>L</m> is the latent heat of fusion
  
Line 51: Line 51:
  
 ==== SIF contents ==== ==== SIF contents ====
 +In this example, ice velocity are in m/s and pressure en MPa. 
 <code> <code>
 Solver 2 Solver 2
Line 71: Line 72:
   Stabilize = True   Stabilize = True
  
-  Exported Variable 1 = String "Phase Change Enthalpy"+  Exported Variable 1 = String "Phase Change Enthalpy" ! (J kg-1)
   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 "temperature"+  Exported Variable 3 = String "temperature" ! (°C)
   Exported Variable 3 DOFs = 1   Exported Variable 3 DOFs = 1
 End End
  
 Constants Constants
- T_ref_enthalpy = real 200.0 + T_ref_enthalpy = real 200.0 !(J kg-1) 
- L_heat = real 334000.0+ L_heat = real 334000.0 !(J kg-1)
  ! Cp(T) = A*T + B  ! Cp(T) = A*T + B
- Enthalpy Heat Capacity A = real 7.253  + Enthalpy Heat Capacity A = real 7.253 !(J kg-1 K-2) 
- Enthalpy Heat Capacity B = real 146.3+ Enthalpy Heat Capacity B = real 146.3 !(J kg-1 K-1) 
 + P_triple = real 0.061173 !Triple point pressure for water (MPa) 
 + P_surf = real 0.1013 ! Surface atmospheric pressure(MPa) 
 + beta_clapeyron = real 0.0974 ! clausus clapeyron relationship (K MPa-1)
 End End
  
Line 94: 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! Material+  Enthalpy Water Diffusivity = real 1.045e-4 ! (kg m-1 s-1)
 End End
  
Line 107: 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 114: Line 118:
   Target Boundaries = 2   Target Boundaries = 2
   Enthalpy_h = variable coordinate 3   Enthalpy_h = variable coordinate 3
-    real MATC "25000.0/150.0*(tx-3250)+140000.0"+    real MATC "25000.0/150.0*(tx-3250)+140000.0" ! (J kg-1) 
 +End
 End End
 </code> </code>
Line 124: Line 129:
 ==== References ==== ==== References ====
  
-Gilbert A., O. Gagliardini, C. Vincent and P. Wagnon, accepted JGR. A 3-D thermal regime model suitable for cold accumulation zones of polythermal mountain glaciers. +GilbertA., O. Gagliardini, C. Vincentand P. Wagnon, 2014. A 3-D thermal regime model suitable for cold accumulation zones of polythermal mountain glaciers, J. Geophys. Res. Earth Surf., 119, doi:10.1002/2014JF003199
solvers/enthalpy.txt · Last modified: 2017/04/21 11:28 by tzwinger
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0