Nearly perfectly matched layer absorber for viscoelastic wave equations

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. T335-T345
Author(s):  
Enjiang Wang ◽  
José M. Carcione ◽  
Jing Ba ◽  
Mamdoh Alajmi ◽  
Ayman N. Qadrouh
Keyword(s):  
Wave Equations ◽  
Pseudospectral Method ◽  
Perfectly Matched Layer ◽  
Zener Model ◽  
Fourier Pseudospectral Method ◽  
Time Stepping ◽  
First Case ◽  
Stress Strain Relation ◽  
Viscoelastic Wave ◽  
Spatial Derivatives

We have applied the nearly perfectly matched layer (N-PML) absorber to the viscoelastic wave equation based on the Kelvin-Voigt and Zener constitutive equations. In the first case, the stress-strain relation has the advantage of not requiring additional physical field (memory) variables, whereas the Zener model is more adapted to describe the behavior of rocks subject to wave propagation in the whole frequency range. In both cases, eight N-PML artificial memory variables are required in the absorbing strips. The modeling simulates 2D waves by using two different approaches to compute the spatial derivatives, generating different artifacts from the boundaries, namely, 16th-order finite differences, where reflections from the boundaries are expected, and the staggered Fourier pseudospectral method, where wraparound occurs. The time stepping in both cases is a staggered second-order finite-difference scheme. Numerical experiments demonstrate that the N-PML has a similar performance as in the lossless case. Comparisons with other approaches (S-PML and C-PML) are carried out for several models, which indicate the advantages and drawbacks of the N-PML absorber in the anelastic case.

Fractional Laplacians viscoacoustic wavefield modeling with k-space-based time-stepping error compensating scheme

Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. T1-T13 ◽  
Author(s):  
Ning Wang ◽  
Tieyuan Zhu ◽  
Hui Zhou ◽  
Hanming Chen ◽  
Xuebin Zhao ◽  
...  
Keyword(s):  
Wave Equation ◽  
Wave Equations ◽  
Time Derivative ◽  
Time Stepping ◽  
Long Time ◽  
Sampling Intervals ◽  
Spatial Derivatives ◽  
The Stability ◽  
Wavefield Extrapolation ◽  
Space Approach

The spatial derivatives in decoupled fractional Laplacian (DFL) viscoacoustic and viscoelastic wave equations are the mixed-domain Laplacian operators. Using the approximation of the mixed-domain operators, the spatial derivatives can be calculated by using the Fourier pseudospectral (PS) method with barely spatial numerical dispersions, whereas the time derivative is often computed with the finite-difference (FD) method in second-order accuracy (referred to as the FD-PS scheme). The time-stepping errors caused by the FD discretization inevitably introduce the accumulative temporal dispersion during the wavefield extrapolation, especially for a long-time simulation. To eliminate the time-stepping errors, here, we adopted the [Formula: see text]-space concept in the numerical discretization of the DFL viscoacoustic wave equation. Different from existing [Formula: see text]-space methods, our [Formula: see text]-space method for DFL viscoacoustic wave equation contains two correction terms, which were designed to compensate for the time-stepping errors in the dispersion-dominated operator and loss-dominated operator, respectively. Using theoretical analyses and numerical experiments, we determine that our [Formula: see text]-space approach is superior to the traditional FD-PS scheme mainly in three aspects. First, our approach can effectively compensate for the time-stepping errors. Second, the stability condition is more relaxed, which makes the selection of sampling intervals more flexible. Finally, the [Formula: see text]-space approach allows us to conduct high-accuracy wavefield extrapolation with larger time steps. These features make our scheme suitable for seismic modeling and imaging problems.

scholarly journals Differential form and numerical implementation of Biot’s poroelasticity equations with squirt dissipation

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. N55-N64 ◽  
Author(s):  
José M. Carcione ◽  
Boris Gurevich
Keyword(s):  
Differential Equations ◽  
Wave Attenuation ◽  
Pseudospectral Method ◽  
Confining Pressure ◽  
Fourier Pseudospectral Method ◽  
High Frequencies ◽  
Splitting Technique ◽  
Stress Strain Relation ◽  
Attenuation Mechanism ◽  
The Time Domain

The squirt-flow wave attenuation mechanism is implemented in Biot’s theory of poroelasticity in the form of differential equations. All the stiffnesses involved in the stress-strain relation become complex and frequency dependent, which can exactly be expressed in terms of kernels based on the Zener mechanical model. In the time domain, this approach implies time convolutions, which are circumvented by introducing memory variables. The differential equations are consistent with Gassmann’s and Mavko-Jizba equations at low and high frequencies, respectively. All the coefficients in the poro-viscoelastic differential equations have a clear physical meaning and can be obtained or estimated from independent measurements. The key additional parameters are the dry-rock bulk modulus at a confining pressure where all the compliant pores are closed, i.e., a hypothetical rock without the soft porosity, the grain-contact aspect ratio and the compliant porosity. We recasted the wave equation in the particle-velocity/stress formulation and solved it by using a time-splitting technique and the Fourier pseudospectral method to compute the spatial derivatives. The algorithm can be used to obtain synthetic wave fields in inhomogeneous media.

A generalization of the Fourier pseudospectral method

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. A53-A56 ◽  
Author(s):  
José M. Carcione
Keyword(s):  
Wave Equation ◽  
Fractional Laplacian ◽  
Differential Operators ◽  
Density Wave ◽  
Pseudospectral Method ◽  
Uniform Density ◽  
Fourier Pseudospectral Method ◽  
Spatial Derivatives ◽  
Fourier Pseudospectral ◽  
Derivatives Of

The Fourier pseudospectral (PS) method is generalized to the case of derivatives of nonnatural order (fractional derivatives) and irrational powers of the differential operators. The generalization is straightforward because the calculation of the spatial derivatives with the fast Fourier transform is performed in the wavenumber domain, where the operator is an irrational power of the wavenumber. Modeling constant-[Formula: see text] propagation with this approach is highly efficient because it does not require memory variables or additional spatial derivatives. The classical acoustic wave equation is modified by including those with a space fractional Laplacian, which describes wave propagation with attenuation and velocity dispersion. In particular, the example considers three versions of the uniform-density wave equation, based on fractional powers of the Laplacian and fractional spatial derivatives.

scholarly journals Prediction of Dynamic Stability Using Mapped Chebyshev Pseudospectral Method

2018 ◽  
Vol 2018 ◽  
pp. 1-14 ◽  
Author(s):  
Jae-Young Choi ◽  
Dong Kyun Im ◽  
Jangho Park ◽  
Seongim Choi
Keyword(s):  
Unsteady Flow ◽  
Three Dimensional ◽  
Pseudospectral Method ◽  
Accurate Method ◽  
Spectral Approach ◽  
Fourier Pseudospectral Method ◽  
Chebyshev Pseudospectral Method ◽  
Fourier Pseudospectral ◽  
Chebyshev Pseudospectral ◽  
Good Agreement

A mapped Chebyshev pseudospectral method is extended to solve three-dimensional unsteady flow problems. As the classical Chebyshev spectral approach can lead to numerical instabilities due to ill conditioning of the spectral matrix, the Chebyshev points are evenly redistributed over the domain by an inverse sine mapping function. The mapped Chebyshev pseudospectral method can be used as an alternative time-spectral approach that uses a Chebyshev collocation operator to approximate the time derivative terms in the unsteady flow governing equations, and the method can make general applications to both nonperiodic and periodic problems. In this study, the mapped Chebyshev pseudospectral method is employed to solve three-dimensional periodic problem to verify the spectral accuracy and computational efficiency with those of the Fourier pseudospectral method and the time-accurate method. The results show a good agreement with both of the Fourier pseudospectral method and the time-accurate method. The flow solutions also demonstrate a good agreement with the experimental data. Similar to the Fourier pseudospectral method, the mapped Chebyshev pseudospectral method approximates the unsteady flow solutions with a precise accuracy at a considerably effective computational cost compared to the conventional time-accurate method.

scholarly journals Convergent Analysis of Energy Conservative Algorithm for the Nonlinear Schrödinger Equation

2015 ◽  
Vol 2015 ◽  
pp. 1-11
Author(s):  
Lv Zhong-Quan ◽  
Gong Yue-Zheng ◽  
Wang Yu-Shun
Keyword(s):  
Schrödinger Equation ◽  
Nonlinear Schrödinger Equation ◽  
Schrodinger Equation ◽  
Pseudospectral Method ◽  
Field Method ◽  
Nonlinear Schrodinger Equation ◽  
Fourier Pseudospectral Method ◽  
Nonlinear Schrödinger ◽  
Vector Field Method ◽  
Energy Conservation Law

Using average vector field method in time and Fourier pseudospectral method in space, we obtain an energy-preserving scheme for the nonlinear Schrödinger equation. We prove that the proposed method conserves the discrete global energy exactly. A deduction argument is used to prove that the numerical solution is convergent to the exact solution in discreteL2norm. Some numerical results are reported to illustrate the efficiency of the numerical scheme in preserving the energy conservation law.

Chebyshev multidomain pseudospectral method to solve cardiac wave equations with rotational anisotropy

2018 ◽  
Vol 09 (04) ◽  
pp. 1850025 ◽  
Author(s):  
Jairo Rodríguez-Padilla ◽  
Daniel Olmos-Liceaga
Keyword(s):  
Wave Propagation ◽  
Numerical Methods ◽  
Wave Equations ◽  
Numerical Solutions ◽  
Three Dimensional ◽  
Pseudospectral Method ◽  
Finite Difference Schemes ◽  
Cardiac Models ◽  
Computational Resources ◽  
Realistic Geometries

The implementation of numerical methods to solve and study equations for cardiac wave propagation in realistic geometries is very costly, in terms of computational resources. The aim of this work is to show the improvement that can be obtained with Chebyshev polynomials-based methods over the classical finite difference schemes to obtain numerical solutions of cardiac models. To this end, we present a Chebyshev multidomain (CMD) Pseudospectral method to solve a simple two variable cardiac models on three-dimensional anisotropic media and we show the usefulness of the method over the traditional finite differences scheme widely used in the literature.

Sign in / Sign up

Export Citation Format

Share Document