SULEC

SULEC is a finite element code that solves the incompressible Navier-Stokes equations for slow creeping flows. Its current characteristics are:

• 2D/3D
• ALE (Arbitrary Lagrangian Eulerian)
• Tracers track materials
• 4-node linear velocity, constant pressure elements
• 9-node quadratic velocity, linear pressure elements
• Thermomechanical
• Uzawa pressure iterations
• Powerlaw creep viscosity
• Elasticity
• Drucker-Prager type plasticity
• Simple implementations of sedimentation and erosion

Tests passed:

• Reproduce lithostatic pressure for gravitational body force model without external forces
• Viscous pure shear (test of free surface)
• Viscous Couette flow (velocity and pressure)
• Poiseuille flow (velocity and pressure)
• Stokes flow (sinking cylinder)
• Cavity flow
• Viscous relaxation of a sinusoidal perturbation
• Pressure field around a rigid inclusion (Schmid and Podladchikov, GJI, 2003)
• Viscoelastic pure shear
• Viscoelastic simple shear
• Recovery of original shape of deformed elastic block
• Mohr-Coulomb shear zones in compression (see also Lemiale et al, PEPI, 2008)
• Halfspace cooling
• Steady-state temperature gradient with basal temperature or basal heatflux
• Schmeling et al subduction benchmark (PEPI 171, 2008)

Developers: Susan Ellis (GNS Sciences, NZ) and Susanne Buiter (NGU)

(c) Susanne Buiter, 25 June 2010