2012/05/25 by Grady I. Lemoine, Lemoine, Grady I., Matthew Ou +4 · 1 citation
Computer Science · Earth and Planetary Sciences · Engineering · Mathematics · #65M08 (Primary) 74S10 #74F10 #74J10 #74L05 #74L15 #86-08 (Secondary) #Advanced Numerical Methods in Computational Mathematics #FOS: Mathematics #Hydraulic Fracturing and Reservoir Analysis #Numerical Analysis (math.NA) #Seismic Imaging and Inversion Techniques #cs.NA #math.NA #msc:65M08 #msc:74F10 #msc:74J10 #msc:74L05 #msc:74L15 #msc:74S10 #msc:86-08
paper · pdf · doi:10.48550/arxiv.1205.5824
26 pages, 7 figures, to appear in SIAM Journal on Scientific Computing Section 2, code provided at https://bitbucket.org/grady_lemoine/poro-2d-cartesian-archive
openalex publication_date 2012/05/25 · arxiv created 2012/11/27 · arxiv updated 2012/11/29 · openalex created_date 2025/10/24 · openalex updated_date 2026/07/28
Poroelasticity theory models the dynamics of porous, fluid-saturated media. It was pioneered by Maurice Biot in the 1930s through 1960s, and has applications in several fields, including geophysics and modeling of in vivo bone. A wide variety of methods have been used to model poroelasticity, including finite difference, finite element, pseudospectral, and discontinuous Galerkin methods. In this work we use a Cartesian-grid high-resolution finite volume method to numerically solve Biot's equations in the time domain for orthotropic materials, with the stiff relaxation source term in the equations incorporated using operator splitting. This class of finite volume method has several useful properties, including the ability to use wave limiters to reduce numerical artifacts in the solution, ease of incorporating material inhomogeneities, low memory overhead, and an explicit time-stepping approach. To the authors' knowledge, this is the first use of high-resolution finite volume methods to model poroelasticity. The solution code uses the CLAWPACK finite volume method software, which also includes block-structured adaptive mesh refinement in its AMRCLAW variant. We present convergence results for known analytic plane wave solutions, achieving second-order convergence rates outside of the stiff regime of the system. Our convergence rates are degraded in the stiff regime, but we still achieve similar levels of error on the finest grids examined. We also demonstrate good agreement against other numerical results from the literature. To aid in reproducibility, we provide all of the code used to generate the results of this paper, at https://bitbucket.org/gradylemoine/poro-2d-cartesian-archive .