Warning: Undefined array key 1 in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 78

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 431

Warning: Undefined array key "DOKU_PREFS" in /home/np29546/public_html/elmerice/wiki/inc/common.php on line 2082

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/common.php on line 2090

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/actions.php on line 38
userfunctions:iceproperties [Elmer/Ice Wiki]

Warning: Undefined array key -1 in /home/np29546/public_html/elmerice/wiki/inc/html.php on line 1458

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
userfunctions:iceproperties [2018/10/28 14:24]
tzwinger [SIF contents]
userfunctions:iceproperties [2018/11/03 12:14]
tzwinger [General Description]
Line 6: Line 6:
  
 ==== General Description ==== ==== General Description ====
-The aim of these user functions is to provide a Fortran version of the else as [[tips:thermoprop|MATC functions]] prescribed material parameters for ice (except for viscosity, which is handled by Glen's flow law). Fortran functions are way faster in execution time, which, in a run that repeatedly calls those parameters, can lead to tremendous speed-ups. Hence, if computing thermo-mechanically coupled problems.+The aim of these user functions is to provide a Fortran version of the else as [[tips:thermoprop|MATC functions]] prescribed material parameters for ice (except for viscosity, which is handled by Glen's flow law). Fortran functions are way faster in execution time, which, in a run that repeatedly calls those parameters, can lead to tremendous speed-ups. Hence, if computing thermo-mechanically coupled problems, rather stick to those. 
 + 
 +All input is expected to be in SI units (Temperature in Kelvin). All outputs by default are in SI units, except the user provides scaling factors (see below).
  
 ===IceConductivity=== ===IceConductivity===
Line 18: Line 20:
  <m>c = 146.3 + (7.253 . T) [J kg^{−1} K^{−1}]</m>\\  <m>c = 146.3 + (7.253 . T) [J kg^{−1} K^{−1}]</m>\\
  
-===GetIcePressureMeltingPoint=== +===IcePressureMeltingPoint=== 
 The pressure melting point of ice as a function of pressure (<m>p</m>) is defined as (in Kelvin):\\ The pressure melting point of ice as a function of pressure (<m>p</m>) is defined as (in Kelvin):\\
  
  <m>T_{pmp} = 273.15 - C_{cc} . max(p, 0) [K]</m>\\  <m>T_{pmp} = 273.15 - C_{cc} . max(p, 0) [K]</m>\\
  
-where <m>C_{cc}</m> is the Clausius Clapeyron constant. In case of negative ice pressures, the function uses zero value.  +where <m>C_{cc}</m> is the Clausius Clapeyron constant. In case of negative ice pressures (actually, any below atmospheric pressures), the function uses the reference value at atmospheric pressure.  
  
  
userfunctions/iceproperties.txt · Last modified: 2018/11/03 12:14 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