GEOS
|
Solving the ODE for seismicity rate from Dieterich, 1994. More...
#include <SeismicityRateKernels.hpp>
Solving the ODE for seismicity rate from Dieterich, 1994.
This solver finds a solution R(x, t) - the seismicity rate - to the ordinary differential equation (ODE) formulated by Dieterich, 1994 given a certain stressing history. The stressing history can consist of mechanical stresses and pore pressure. The solver class includes a member variable pointing to the stress solver that is specified in the XML file. SolverStep for the stress solver is then called in the SolverStep function for the seismicity rate, to take the updated stress history as the input.
Solving the ODE is currently implemented by computing the closed-form integral solution to the ODE which involves numerical calculation of an integral of a stress functional. We initially solve for the log of the seismicity rate in order to avoid overflow that typically occurs in the exponential of the stress history.