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: 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]

This is an old revision of the document!


User Function IceProperties

General Informations

  • USF Fortran File: USF_IceProperties.f90
  • USF Name: IceConductivity, IceCapacity and IcePressureMeltingPoint
  • Required Input Variable(s): Temperature (IceConductivity, IceCapacity), Pressure (IcePressureMeltingPoint)

General Description

The aim of these user functions is to provide a Fortran version of the else as 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 therm-mechanically coupled problems.

IceConductivity

The heat conductivity of ice as a function of temperature (T) is defined (in SI units) as:

\kappa_{ice} = 9.828 . exp(-5.7^{-3}.T) [W m^{−1} K^{−1}]

IceCapacity

The capacity of ice as a function of temperature (T) is defined (in SI units) as:

c = 146.3 + (7.253 . T) [J kg^{−1} K^{−1}]

GetIcePressureMeltingPoint

The pressure melting point of ice as a function of pressure (p) is defined as (in Kelvin):

T_{pmp} = 273.15 - C_{cc} . max(p, 0) [K]

where C_{cc} is the Clausius Clapeyron constant. In case of negative ice pressures, the function uses zero value.

SIF contents

The required keywords in the SIF file for these user functions are:

Constants
   Clausius Clapeyron = Real 9.8e-08
End

Material 1
  Name = "ice"

  ! Heat transfer stuff (converted to MPa-m-a system)
  Temp Heat Capacity = Variable Temp
    Real Procedure "ElmerIceUSF" "IceCapacity"
  Heat Capacity Scaling Factor = Real $(secondsperyear)^(2.0)

  Temp Heat Conductivity = Variable Temp
    Real Procedure "ElmerIceUSF" "IceConductivity"  
  Heat Conductivity Scaling Factor = Real $(secondsperyear)*1.0E-06 
  
  Temp Upper Limit = Variable HydroPressure
    Real Procedure "ElmerIceUSF" "IcePressureMeltingPoint"
  Pressure Scaling Factor = Real 1.0E06 ! from MPa to Pa
End

Examples

An example demonstrating the use of the thermal properties of ice can be found in [ELMER_TRUNK]/elmerice/Tests/TemperateIceTestFct.

References

Ritz, C. 1987. Time dependent boundary conditions for calculation of temperature fields in ice sheets. In: E. D. Waddington and J. S. Walder (Eds.), The Physical Basis of Ice Sheet Modelling, IAHS Publication No. 170, pp. 207–216. IAHS Press, Wallingford, UK.

userfunctions/iceproperties.1540736578.txt.gz · Last modified: 2018/10/28 14:22 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