Viscoelasticity in Underworld2

Rebecca Farrington

Jun 24, 2016

A short introduction to using viscoelasticity

Viscoelastic materials exhibit the properies of both solids and liquids, with deformation rates dependent on both the viscous stress and elastic stress rate. In a Maxwell viscoelastic material the strain rate is proportional to the sum of the stress and the stress rate ,

\begin{equation} D_{ij} = D^v_{ij} + D^e_{ij} = \frac{\tau_{ij}}{2\eta} + \frac{\dot\tau_{ij}}{2\mu} \end{equation}

with viscosity and shear modulus. In Underworld this rheology can be implemented using an effective viscosity and an additional force term within the momentum equation,

\begin{equation} (2\eta_{eff}D_{ij})_{,j} <!–

with pressure, the effective viscosity, the Maxwell relaxation time, the elastic time step and the stress history term. See Farrington et al (2014) for full details.

Viscoelasticity in simple shear

A uniform viscoelastic material undergoes simple shear at a constant rate until . The shearing velocity is then taken to zero with the stress decaying with time.

Viscoelastic stress history term for differing relaxation time

The graph above shows the shear component of the stress history term from equation (2) for a range of Maxwell relaxation times (alpha). See the viscoelastic example for details, my_uw/doc/2_11_ViscoelasticityInSimpleShear.ipynb

An objective stress rate tensor

The stress history term (stress rate tensor) is carried on particles that can be advected and rotated through time by the velocity field. To maintain objectivity it must be both advected and rotated for the current reference frame. This can be acheived using the Jaumann stress rate, defined by with the vorticity tensor. If significant rotation occurs the use of the Jaumann stress rate is encouranged. There are other terms which may be added to this definition accounting for other changes in reference frame including stretching. However, including these additional terms couples stress and pressure, adding a considerable complication that is not required for mathematical consistency.

For full details regarding the background, numerics and use in geodynamic modelling please see Farrington et al (2014) and the references within.

References

Farrington, R. J., L.-N. Moresi, and F. A. Capitanio (2014), The role of viscoelasticity in subducting plates, Geochem. Geophys. Geosyst., 15, 4291–4304