Benchmark Description

Benchmark problem description.

Summary

Viscoelastic (Maxwell) relaxation of stresses from the 1992 M7.3 Landers earthquake, focusing on the deformation in the area of the 1999 M7.1 Hector Mine earthquake.

Problem Specification

Model size: [NEED SPECS FOR CARL'S MESH]

Material properties

Elastic

The material properties are a simplified 1-D version of the 3-D SCEC Community Velocity Model. The 1-D model contains 11 layers with uniform material properties within each layer and a minimum layer thickness of 2 km. The elastic properties are given in an ASCII file.

Viscoelastic

Maxwell linear viscoelasticity (based on values in Pollitz, EPSL, 2003)

Upper crust (-19 km ≤ z): η = 1.0e+25 Pa-s (essentially elastic)
Lower crust (-30 km ≤ z < -19 km): η = 32.2e+18 Pa-s
Mantle (z < -30 km): η = 4.6e+18 Pa-s

Fault geometry and slip distribution

The Landers and Hector Mine fault geometries and slip distribution for Landers are incorporated into the LaGriT mesh.

Boundary conditions

Bottom and side displacements are pinned. Top of the model is a free surface.

Discretization

[GET SPECS FROM CARL'S MESH]

Element types

Linear and/or quadratic tetrahedral elements.

Requested Output

Solution

Displacements at all nodes at times of 0, 0.5, 1, 2, 4, and 7 years as well as the mesh topology (i.e., element connectivity arrays and coordinates of vertices) and basis functions. Also compute the traction vector computed at the quadrature points of the faces making up the Hector Mine faults.

June 30, 2006: Use ASCII output for now. In the future we will switch to using HDF5 files.

Sections

Personal tools