Stress Balance Solution

Physical basis

Conservation of linear momentum

The conservation of momentum reads:

where:

• is the ice density
• is the velocity vector
• is the Cauchy stress tensor

• is a body force Now if we assume that:

• The ice motion is a Stokes flow (acceleration negligible)
• The only body force is due to gravity (Coriolis effect negligible) The equation of momentum conservation is reduced to:

Conservation of angular momentum

For a non-polar material body, the balance of angular momentum imposes the stress tensor to be symmetrical:

Ice constitutive equations

Ice is treated as a purely viscous incompressible material [Cuffey2010]. Its constitutive equation therefore only involves the deviatoric stress and the strain rate tensor:

where:

• is the deviatoric stress tensor ()
• is the ice effective viscosity
• is the strain rate tensor Ice is a non-Newtonian fluid, its viscosity follows the generalized Glen’s flow law [Glen1955]:

where:

• is the ice hardness or rigidity
• is Glen’s flow law exponent, generally taken as equal to 3
• is the effective strain rate The effective strain rate is defined as:

where is the Frobenius norm.

Full-Stokes (FS) field equations

Without any further approximation, the previous system of equations are called the Full-Stokes model.

Higher-Order (HO) field equations

We make two assumptions:

1. Bridging effects are neglected
2. Horizontal gradient of vertical velocities are neglected compared to vertical gradients of horizontal velocities With these two assumptions, the Full-Stokes equations are reduced to a system of 2 equations with 2 unknowns [Blatter1995, Pattyn2003]:

with:

Shelfy-Stream Approximation (SSA) field equations

We make the following assumption:

1. Vertical shear is negligible With this assumption, we have a system of 2 equations with 2 unknowns in the horizontal plane [Morland1987a, MacAyeal1989]:

with:

where:

• is the depth-averaged viscosity
• is the ice thickness
• is the basal friction coefficient

Boundary conditions

At the surface of the ice sheet, , we assume a stress-free boundary condition. A viscous friction law is applied at the base of the ice sheet, , and water pressure is applied at the ice/water interface . For FS, these boundary conditions are:

where

• is the outward-pointing unit normal vector
• is the water density
• is the vertical coordinate with respect to sea level

For HO, these boundary conditions become:

where .

For SSA, these boundary conditions are:

Model parameters

The parameters relevant to the stress balance solution can be displayed by typing:

>> md.stressbalance

• md.stressbalance.restol: mechanical equilibrium residue convergence criterion
• md.stressbalance.reltol: velocity relative convergence criterion, (NaN if not applied)
• md.stressbalance.abstol: velocity absolute convergence criterion, (NaN if not applied)
• md.stressbalance.maxiter: maximum number of nonlinear iterations (default is 100)
• md.stressbalance.spcvx: x-axis velocity constraint (NaN means no constraint)
• md.stressbalance.spcvy: y-axis velocity constraint (NaN means no constraint)
• md.stressbalance.spcvz: z-axis velocity constraint (NaN means no constraint)
• md.stressbalance.rift_penalty_threshold: threshold for instability of mechanical constraints
• md.stressbalance.rift_penalty_lock: number of iterations before rift penalties are locked
• md.stressbalance.penalty_factor: offset used by penalties:

• md.stressbalance.vertex_pairing: pairs of vertices that are penalized
• md.stressbalance.shelf_dampening: use dampening for floating ice? Only for Stokes model
• md.stressbalance.referential: local referential
• md.stressbalance.requested_outputs: additional outputs requested

The solution will also use the following model fields:

• md.flowequations: FS, HO or SSA
• md.materials: material parameters
• md.initialization.vx: x component of velocity (used as an initial guess)
• md.initialization.vy: y component of velocity (used as an initial guess)
• md.initialization.vz: y component of velocity (used as an initial guess)

Running a simulation

To run a simulation, use the following command:

>> md = solve(md, 'Stressbalance');

The first argument is the model, the second is the nature of the simulation one wants to run.

References

• H. Blatter. Velocity And Stress-Fields In Grounded Glaciers: A Simple Algorithm For Including Deviatoric Stress Gradients. J. Glaciol., 41(138):333-344, 1995.

• K. M. Cuffey and W. S. B. Paterson. The Physics of Glaciers, 4th Edition. Elsevier, Oxford, 2010.

• J. W. Glen. The creep of polycrystalline ice. Proc. R. Soc. A, 228(1175):519-538, 1955.

• L. W. Morland. Unconfined ice shelf flow. Proceedings of Workshop on the Dynamics of the West Antarctic Ice Sheet, University of Utrecht, May 1985. Published by Reidel, ed. C.J. van der Veen and J. Oerlemans:99-116, 1987.