Skip to content
forked from wqseis/drdg3d

Discontinuous Galerkin Finite Element Code for 3D Dynamic Rupture Modeling

Notifications You must be signed in to change notification settings

bjmarfito/drdg3d

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

37 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

DRDG3D DOI

For 3D dynamic rupture (DR) modelling, using a nodal discontinuous Galerkin (DG) finite element method [1].

Installation

Here we show how to install the software in MacOS. For the installation on Linux system, please refer to manual_drdg3d.pdf.

Before the compilation of the software, make sure you have make, gcc, openmpi, netcdf-fortran, gmsh, meshio and metis installed.

> brew install make
> brew install gcc
> brew install open-mpi
> brew install netcdf
> brew install netcdf-fortran
> pip install gmsh
> pip install meshio
> pip install metis

Compilation

After installing requirements, the program is installed as follows.

> git clone https://github.com/wqseis/drdg3d.git
> cd drdg3d
> mkdir -p bin obj
> make

Three executable files will be generated in the bin directory: exe_solver, exe_get_neigh and exe_part_mesh.

Before running the examples, set the path first:

> bash scripts/setenv.sh

Usage

We take tpv5 as an example:

> cd examples/tpv5
  1. Generate the mesh:
> bash genmesh.sh
  1. Run the preprocess in MATLAB:

    run_preprocess.m

  2. Run the program

> bash go.sh

Plot the results

Use draw_fault_snap.m to draw the snapshots of fault slip rate (or slip, stresses, etc.). There are other scripts such as draw_fault_seismo.m, draw_grdsurf_snap.m, draw_grdsurf_seismo.m, please refer to manual_drdg3d.pdf for more details.

Reference

[1] Wenqiang Zhang, Yajing Liu, Xiaofei Chen. A Mixed-Flux-Based Nodal Discontinuous Galerkin Method for 3D Dynamic Rupture Modeling. ESS Open Archive. November 01, 2022. DOI: 10.1002/essoar.10512657.1

[2] Zhang, W., Liu, Y., & Chen, X. (2023). A mixed-flux-based nodal discontinuous Galerkin method for 3D dynamic rupture modeling. Journal of Geophysical Research: Solid Earth, 128, e2022JB025817. https://doi.org/10.1029/2022JB025817

About

Discontinuous Galerkin Finite Element Code for 3D Dynamic Rupture Modeling

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • Fortran 84.8%
  • MATLAB 14.4%
  • Other 0.8%