Arrhenius law
Ver:1
DESCRIPTION
S = 2 * eta * D
where S is the extra-stress tensor, D is the rate of deformation tensor, and eta is the viscosity.
The temperature dependence
of the viscosity is denoted by a function H(T), where T is the temperature.
H(T) comes as a factor
in the viscosity law:
eta(g,T) = F(g)*H(T)
For the Arrhenius law, the function H(T) is written as follows :
H(T) = exp [ alpha / ( T - T0 ) - alpha / ( Talpha - T0 ) ]
where alpha is the
ratio of the activation energy to the Blotzman constant, Talpha is a reference
temperature,
and T0 is a scaling
temperature when non-absolute temperature scales are used.
By default, alpha =
Talpha = T0 = 0.
| -1 | Upper level menu |
| 1 | Modify alfa = 0.0000000E+00 |
| 2 | Modify talfa = 0.0000000E+00 |
| 3 | Modify t0 = 0.0000000E+00 |
A high value of the
parameter alpha leads to strong non-linearities. A Picard iteration is
recommended
instead of a full
Newton iteration. In case the program diverges, an evolution scheme (on
alpha) must be used.
The units for alpha,
Talpha and T0 are :
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|