Step 3: Multiple Earthquake Ruptures and Velocity Boundary Conditions

Features

• Triangular cells

• pylith.meshio.MeshIOPetsc

• pylith.problems.TimeDependent

• pylith.materials.Elasticity

• pylith.materials.IsotropicLinearElasticity

• pylith.faults.FaultCohesiveKin

• pylith.faults.KinSrcStep

• field split preconditioner

• Schur complement preconditioner

• pylith.bc.DirichletTimeDependent

• spatialdata.spatialdb.UniformDB

• pylith.meshio.OutputSolnBoundary

• pylith.meshio.DataWriterHDF5

• Quasi-static simulation

• spatialdata.spatialdb.SimpleDB

Simulation parameters

This example involves a quasistatic simulation that solves for the deformation from velocity boundary conditions and multiple earthquake ruptures on the fault. The velocity boundary conditions match those in Step 2. We prescribe the first earthquake rupture to occur at 100 years with 1 meter of right-lateral slip and the second earthquake rupture to occur at 200 years with 3 meters of right-lateral slip. Fig. 63 shows the boundary conditions on the domain. The parameters specific to this example are in step03_multislip_velbc.cfg.

Fig. 63 Boundary conditions for quasi-static simulation with velocity boundary conditions and coseismic slip. We set the x displacement to zero on the +x and -x boundaries. We set the y velocity to -1 cm/yr on the +x boundary and +1 cm/yr on the -x boundary. We prescribe 1 meter of right-lateral slip to occur at 100 years and 3 meters of right-lateral slip to occur at 200 years.

Listing 83 Prescribed slip rupture parameters for Step 3 with two earthquake ruptures on the strike-slip fault.

[pylithapp.problem.interfaces.fault]
eq_ruptures = [one, two]

[pylithapp.problem.interfaces.fault.eq_ruptures.one]
db_auxiliary_field = spatialdata.spatialdb.UniformDB
db_auxiliary_field.description = Fault rupture one
db_auxiliary_field.values = [initiation_time, final_slip_left_lateral, final_slip_opening]
db_auxiliary_field.data = [100.0*year, -1.0*m, 0.0*m]

[pylithapp.problem.interfaces.fault.eq_ruptures.two]
db_auxiliary_field = spatialdata.spatialdb.UniformDB
db_auxiliary_field.description = Fault rupture two
db_auxiliary_field.values = [initiation_time, final_slip_left_lateral, final_slip_opening]
db_auxiliary_field.data = [200.0*year, -3.0*m, 0.0*m]

Running the simulation

Listing 84 Run Step 3 simulation

$ pylith step03_multislip_velbc.cfg

# The output should look something like the following.
 >> /software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/apps/PyLithApp.py:77:main
 -- pylithapp(info)
 -- Running on 1 process(es).
 >> /software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/meshio/MeshIOObj.py:38:read
 -- meshiopetsc(info)
 -- Reading finite-element mesh
 >> /src/cig/pylith/libsrc/pylith/meshio/MeshIO.cc:85:void pylith::meshio::MeshIO::read(pylith::topology::Mesh *, const bool)
 -- meshiopetsc(info)
 -- Component 'reader': Domain bounding box:
    (-50000, 50000)
    (-75000, 75000)
# -- many lines omitted --

25 TS dt 0.1 time 2.4
    0 SNES Function norm 1.078084164657e-03
    Linear solve converged due to CONVERGED_ATOL iterations 3
    1 SNES Function norm 1.269009043820e-12
  Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 1
26 TS dt 0.1 time 2.5
 >> /software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/problems/Problem.py:199:finalize
 -- timedependent(info)
 -- Finalizing problem.

The beginning of the output written to the terminal is identical to that from Step 1. At the end of the output, we see that the simulation advanced the solution 26 time steps. Remember that the PETSc TS monitor shows the nondimensionalized time and time step values.

Visualizing the results

In Fig. 64 we use the pylith_viz utility to visualize the y displacement field. You can move the slider or use the p and n keys to change the increment or decrement time.

Listing 85 Visualize PyLith output using pylith_viz.

pylith_viz --filename=output/step03_multislip_velbc-domain.h5 warp_grid --component=y

Fig. 64 Solution for Step 3 at t=190 yr. The colors of the shaded surface indicate the y displacement, and the deformation is exaggerated by a factor of 1000. The undeformed configuration is shown by the gray wireframe.