BISICLES AMR ice sheet model
0.9
|
Arrhenius Rate Factor. More...
#include <ConstitutiveRelation.H>
Public Member Functions | |
ArrheniusRateFactor (Real a_seconds_per_unit_time) | |
void | setDefaultParameters (Real a_seconds_per_unit_time) |
void | setParameters (Real a_n, Real a_enhance, Real a_B0, Real a_theta_r, Real a_K, Real a_C, Real a_R, Real a_Q) |
void | computeA (FArrayBox &a_A, const FArrayBox &a_thetaStar, const FArrayBox &a_pressure, const Box &a_box) const |
RateFactor * | getNewRateFactor () const |
![]() | |
virtual | ~RateFactor () |
Arrhenius Rate Factor.
The Arrhenius rate factor is used widely in ice sheet modelling, e.g (Pattyn, 2003)
ArrheniusRateFactor::ArrheniusRateFactor | ( | Real | a_seconds_per_unit_time | ) |
|
virtual |
Implements RateFactor.
|
virtual |
Implements RateFactor.
void ArrheniusRateFactor::setDefaultParameters | ( | Real | a_seconds_per_unit_time | ) |
power law exponent
Pattyns enhancement factor
flow rate factor (2.207 Pa a^(1/n) )
limit temperature in flow-rate factor (273.39 K)
flow rate exponent (1.17)
flow rate factor 0.16612 K^m_K
universal gas constant 8.31 J/(mol K) )
activation energy for creep (7.88e4 J/mol
void ArrheniusRateFactor::setParameters | ( | Real | a_n, |
Real | a_enhance, | ||
Real | a_B0, | ||
Real | a_theta_r, | ||
Real | a_K, | ||
Real | a_C, | ||
Real | a_R, | ||
Real | a_Q | ||
) |
power law exponent
Pattyn's "enhancement factor" = 1.0
flow rate factor (2.207 Pa a^(1/n) )
limit temperature in flow-rate factor (273.39 K)
flow rate exponent (1.17)
flow rate factor 0.16612 K^m_K
universal gas constant 8.31 J/(mol K) )
activation energy for creep (7.88e4 J/mol