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 |
Friction |
Number |
Further details |
---|---|---|---|---|---|
dynamic |
SS |
LSW |
1 |
3 stress asperities, see SCEC TPV5 |
|
dynamic |
SS |
LSW |
1 |
Bi-material fault, heterogeneous initial stress, see SCEC TPV6 |
|
dynamic |
N |
LSW |
1 |
depth-dependent initial stress conditions, see SCEC TPV12 |
|
dynamic |
N |
LSW |
1 |
Same as TPV12 with non-associative Drucker-Prager plastic |
|
dynamic |
SS |
LSW |
1 |
Randomly-generated heterogeneous initial stress conditions, |
|
dynamic |
SS |
LSW |
2 |
Rightward branch forming a 30 degree angle, see SCEC TPV24 |
|
dynamic |
SS |
LSW |
1 |
Stochastic roughness, see SCEC TPV29 |
|
dynamic |
SS |
LSW |
1 |
Imperial Fault model with 3D velocity structure, see SCEC TPV34 |
|
dynamic |
SS |
fvw-RS |
1 |
see SCEC TPV104 |
|
point |
n/a |
n/a |
n/a |
point-source benchmark, see SISMOWINE WP2_LOH1 |
|
kinematic |
R |
n/a |
1 |
- 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.