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
solvers:enthalpy [2014/06/20 13:08]
gag [General Description]
solvers:enthalpy [2017/04/21 11:28] (current)
tzwinger [Examples]
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, defined as <m>C_p = AT+Bp</m>+  * <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 120: Line 125:
  
 ==== Examples ==== ==== Examples ====
-An example solving for the enthalpy within the Tete Rousse glacier assuming an elevation dependant enthalpy at the upper surface can be found in ''[ELMER_TRUNK]/elmerice/Tests/Enthalpy''.+An example solving for the enthalpy within the Tete Rousse glacier assuming an elevation dependent enthalpy at the upper surface can be found in ''[ELMER_TRUNK]/elmerice/Tests/Enthalpy''.
  
 ==== 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.1403269709.txt.gz · Last modified: 2014/06/20 13:08 by gag
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