Search

In this work, we investigate the propagation of electromagnetic waves in the Cole-Cole dispersive medium by using the discontinuous Galerkin (DG) method to solve the coupled time-domain Maxwell's equations and polarization equation. We define a new and sharpened total energy function for the Cole-Cole model, which better describes the behaviors of the energy than what is available in the current literature. A major theme in the time-domain numerical modeling of this problem has been tackling the difficulty of handling the nonlocal term involved in the time-domain polarization equation. Based on the diffusive representation and the quadrature formula, we derive an approximate system, where the convolution kernel is replaced by a finite number of auxiliary variables that satisfy local-in-time ordinary differential equations. To ensure the resulted approximate system is stable, a nonlinear constrained optimization numerical scheme is established to determine the quadrature coefficients. By a special choice of the numerical fluxes and projections, we obtain {for the constant coefficient case } an optimal-order convergence result for the semi-discrete DG scheme. The temporal discretization is achieved by the standard two-step backward difference formula and a fast algorithm with linear complexity is constructed. Numerical examples are provided for demonstrating the efficiency of the proposed algorithm, validating the theoretical results and illustrating the behaviors of the energy., Comment: 24 pages, 6 figures, 3 tables

This paper concerns the Time-Domain Full Waveform Inversion (FWI) for dispersive and dissipative poroelastic materials. The forward problem is an initial boundary value problem (IBVP) of the poroelastic equations with a memory term; the FWI is formulated as a minimization problem of a least-square misfit function with the (IBVP) as the constraint. In this paper, we derive the adjoint problem of this minimization problem, whose solution can be applied to computed the direction of steepest descent in the iterative process for minimization. The adjoint problem has a similar numerical structure as the forward problem and hence can be solved by the same numerical solver. Because the tracking of the energy evolution plays an important role in the FWI for dissipative and dispersive equations, the energy analysis of the forward system is also carried out in this paper., Comment: 14 pages

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.