A discontinuous Galerkin method for wave propagation in orthotropic poroelastic media with memory terms

In this paper, we investigate wave propagation in orthotropic poroelastic media by studying the time-domain poroelastic equations. Both the low frequency Biot's (LF-Biot) equations and the Biot-Johnson-Koplik-Dashen (Biot-JKD) models are considered. In LF-Biot equations, the dissipation terms are proportional to the relative velocity between the fluid and the solid by a constant. Contrast to this, the dissipation terms in the Biot-JKD model are in the form of time convolution (memory) as a result of the frequency-dependence of fluid-solid interaction at the underlying microscopic scale in the frequency domain. The dynamic tortuosity and permeability described by Darcy's law are two crucial factors in this problem, and highly linked to the viscous force. In the Biot model, the key difficulty is to handle the viscous term when the pore fluid is viscous flow. In the Biot-JKD dynamic permeability model, the convolution operator involves order $1/2$ shifted fractional derivatives in the time domain, which is challenging to discretize. In this work, we utilize the multipoint Pad$\acute{e}$ (or Rational) approximation for Stieltjes function to approximate the dynamic tortuosity and then obtain an augmented system of equations which avoids storing the solutions of the past time. The Runge-Kutta discontinuous Galerkin (RKDG) method is used to compute the numerical solution, and numerical examples are presented to demonstrate the high order accuracy and stability of the method

A CDG-FE method for the two-dimensional Green-Naghdi model with the enhanced dispersive property

In this work, we investigate numerical solutions of the two-dimensional shallow water wave using a fully nonlinear Green-Naghdi model with an improved dispersive effect. For the purpose of numerics, the Green-Naghdi model is rewritten into a formulation coupling a pseudo-conservative system and a set of pseudo-elliptic equations. Since the pseudo-conservative system is no longer hyperbolic and its Riemann problem can only be approximately solved, we consider the utilization of the central discontinuous Galerkin method which possesses an important feature of needlessness of Riemann solvers. Meanwhile, the stationary elliptic part will be solved using the finite element method. Both the well-balanced and the positivity-preserving features which are highly desirable in the simulation of the shallow water wave will be embedded into the proposed numerical scheme. The accuracy and efficiency of the numerical model and method will be illustrated through numerical tests.Comment: 27 pages, 35 figure