Overview

This documentation is a collection of useful dynamic simulation examples to help users build models from scratch with little effort. Each example is demonstrated carefully with geometry building, parameter setup, and result visualization. Users are suggested to repeat each example in order to get a comprehensive idea of how to set up dynamic simulation models with SeisSol.

SeisSol is a part of SCEC dynamic code validation project (Harris et al. 2018) (http://strike.scec.org/cvws/). Here we show several SCEC benchmarks for beginners to quickly catch up with SeisSol workflow. Each benchmark example is composed of a short problem description, a section of geometry, initial setups (stress, nucleation, friction, etc.), and simulation results.

Please note that the examples used here are only for demonstration purpose. For detailed benchmark tests please refer to SCEC benchmark center.

Description

Source

Faulting
mechanism 1

Friction
law 2

Number
of faults

Further details

TPV5

dynamic

SS

LSW

1

3 stress asperities, see SCEC TPV5

TPV6

dynamic

SS

LSW

1

Bi-material fault, heterogeneous initial stress, see SCEC TPV6

TPV12

dynamic

N

LSW

1

depth-dependent initial stress conditions, see SCEC TPV12

TPV13

dynamic

N

LSW

1

Same as TPV12 with non-associative Drucker-Prager plastic
with yielding in shear, see SCEC TPV13

TPV16

dynamic

SS

LSW

1

Randomly-generated heterogeneous initial stress conditions,
see SCEC TPV16/17

TPV24

dynamic

SS

LSW

2

Rightward branch forming a 30 degree angle, see SCEC TPV24

TPV29

dynamic

SS

LSW

1

Stochastic roughness, see SCEC TPV29

TPV34

dynamic

SS

LSW

1

Imperial Fault model with 3D velocity structure, see SCEC TPV34

TPV104

dynamic

SS

fvw-RS

1

see SCEC TPV104

LOH.1

point

n/a

n/a

n/a

point-source benchmark, see SISMOWINE WP2_LOH1

Northridge

kinematic

R

n/a

1

see Kinematic source example - 1994 Northridge earthquake

1

SS: strike-slip, N: normal, R: reverse, O: oblique

2

LSW: linear slip-weakening friction, fvw-RS: fast-velocity weakening rate-and-state friction

Prerequisites

Before you begin any of the examples, you will need to install the latest SeisSol from (https://github.com/SeisSol/SeisSol). The instruction can be found at https://seissol.readthedocs.io/en/latest/compilation.html. All geometry and tetrahedral meshes are generated using free software Gmsh (http://gmsh.info/). If you do not wish to create your own mesh at this time, the meshes are also provided as part of the example. The ParaView visualization package (https://www.paraview.org/) may be used to view simulation results. You may use other visualization software, but some adaptions from what is described here will be necessary. Furthermore, you can complete a subset of the example using files provided (as described below), skipping the steps for which you do not have the proper software packages installed.

Input file resources

The files needed to work through the examples are provided. All files necessary to set up the cookbook examples can be downloaded at https://github.com/SeisSol/Examples

References

Harris, R. A., Michael Barall, B. T. Aagaard, S. Ma, and K. B. O. Daniel Roten, Benchun Duan, Dunyu Liu, Bin Luo, Kangchen Bai, Jean-Paul Ampuero, Yoshihiro Kaneko, Alice-Agnes Gabriel, Kenneth Duru, Thomas Ulrich, Stephanie Wollherr, Zheqiang Shi, Eric Dunham, Sam Bydlon, Zhenguo Zhang, Xiaofei Chen, Surendra N. Somala, Christian Pelties, Josue Tago, Victor Manuel Cruz-Atienza, Jeremy Kozdon, Eric Daub, Khurram Aslam, Yuko Kase, Kyle Withers (2018), A Suite of Exercises for Verifying Dynamic Earthquake Rupture Codes, Seismol. Res. Lett., 89(3), 1146-1162, doi:10.1785/0220170222.