January 25, 2007
Benchmark problem description. Formerly known as benchmark 5b.
Viscoelastic (Maxwell) relaxation of stresses from a single, finite, reverse-slip earthquake in 3-D without gravity. Evaluate results with imposed displacement boundary conditions on a cube with sides of length 24 km. The displacements imposed are the analytic elastic solutions. Symmetry boundary conditions are imposed at y = 0, so the solution is equivalent to that for a domain with a 48 km length in the y direction.
Bottom and side displacements set to analytic solution. (Note: the side at y = 0 km has zero y- displacements because of symmetry.) Top of the model is a free surface.
The model should be discretized with nominal spatial resolutions of 1000 m, 500 m, and 250 m. If possible, also run the models with a nominal spatial resolution of 125 m. Optionally, use meshes with variable (optimal) spatial resolution with the same number of nodes as the uniform resolution meshes.
Linear and/or quadratic and tetrahedral and/or hexahedral.
Displacements at all nodes at times of 0, 1, 5, and 10 years as well as the mesh topology (i.e., element connectivity arrays and coordinates of vertices) and basis functions.
Okada routines are available to generate an elastic solution. The ‘best’ viscoelastic answer will be derived via mesh refinement. Analytical solutions to the viscoelastic solution are being sought if anyone has information.