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/compiling-seissol.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.