Input Parameters
In the example input files for benchmark problem 1 (Thrase.jl/examples/bp1-qd.dat) and benchmark problem 6 (Thrase.jl/examples/bp6.dat) there are multiple parameters that can be updated to change the problem. They are desribed in the table below.
Sometimes the benchmark problem 1 description paper, and benchmark problem 6 description paper contain slightly different names than the input file, we try to indicate in the description in parantheses if this is the case.
| Input File Parameter | Units | Data Type | Description |
|---|---|---|---|
| pth | - | string | path to the directory that stores output data |
| generate_plots | - | boolean | option to generate example plots |
| stride_space | - | Real <: Number | write-out every "stride_space" grid points along fault |
| stride_time | - | Real <: Number | write-out every "stride_time" time steps |
| Lx1 | $km$ | Real <: Number | x-domain is (Lx1, Lx2) |
| Lx2 | $km$ | Real <: Number | x-domain is (Lx1, Lx2) |
| Lz1 | $km$ | Real <: Number | z-domain is (Lz1, Lz2) |
| Lz2 | $km$ | Real <: Number | z-domain is (Lz1, Lz2) |
| Nx (Nr) | - | Real <: Number | number of grid points in each dimension |
| Nz (Ns) | - | Real <: Number | number of grid points in each dimension |
| Hx | - | Real <: Number | dimensions in which constant grid-spacing is applied |
| Hz | - | Real <: Number | dimensions in which constant grid-spacing is applied |
| dx | $km$ | Real <: Number | desired grid spacing in [0, Hx] in physical domain |
| dz | $km$ | Real <: Number | desired grid spacing in[-Hz, Hz] in physical domain |
| el_r | - | Real <: Number | length scales over which grid stretching occurs |
| el_s | - | Real <: Number | length scales over which grid stretching occurs |
| sim_years | $years$ | Real <: Number | how many years to simulate ($t_f$) |
| Vp | $m/s$ | Real <: Number | plate rate |
| ρ | $kg/m^3$ | Real <: Number | density |
| cs | $km/s$ | Real <: Number | shear wave speed |
| σn (σn_0) | $MPa$ | Real <: Number | initial effective normal stress on fault |
| RSamin | - | Real <: Number | rate-and-state parameter ($a_0$) |
| RSamax | - | Real <: Number | rate-and-state parameter ($a_max$) |
| RSb | - | Real <: Number | rate-and-state parameter ($b_0$) |
| RSDc | $m$ | Real <: Number | critical slip distance ($D_c$) |
| RSf0 | $m/s$ | Real <: Number | reference friction coefficient (f) |
| RSV0 | $m/s$ | Real <: Number | reference slip rate initial slip rate ($V_0$) |
| RSVinit | $m/s$ | Real <: Number | initial slip rate ($V_{init}$) |
| RSH1 | $km$ | Real <: Number | depth extent of uniform VS region (H) |
| RSH2 | $km$ | Real <: Number | width of VS-VS transition zone (h) |
| RSWf | $km$ | Real <: Number | width of rate-and-state fault ($W_f$) |
| RSD_RS | $mm$ | Real <: Number | state evolution distance |
| RSLf | $km$ | Real <: Number | half-length of rate-and-state fault |
| τ_init | $MPa$ | Real <: Number | initial shear stress |
| q_0 | $m/s$ | Real <: Number | fluid injection rate |
| t_off | $days$ | Real <: Number | injection turn-off time |
| α | $m^2/s$ | Real <: Number | hydraulic diffusivity |
| β | $Pa^{-1}$ | Real <: Number | pore and fluid compressibility |
| φ | - | Real <: Number | porosity |
| k | $m^2$ | Real <: Number | permeability |
| η_visc | $Pa\cdot s$ | Real <: Number | fluid viscosity |
| SBPp | - | Real <: Number | SBP interior spatial order of accuracy |