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.

The wave equation in generalized coordinates

Geophysics ◽  
1994 ◽  
Vol 59 (12) ◽  
pp. 1911-1919 ◽  
Author(s):  
José M. Carcione
Keyword(s):  
Wave Equation ◽  
Pseudospectral Method ◽  
Special Treatment ◽  
Lamb’S Problem ◽  
Generalized Coordinates ◽  
Topographic Features ◽  
Boundary Treatment ◽  
Collocation Scheme ◽  
Spatial Derivatives ◽  
Derivatives Of

This work introduces a spectral collocation scheme for the viscoelastic wave equation transformed from Cartesian to generalized coordinates. Both the spatial derivatives of field variables and the metrics of the transformation are calculated by the Chebychev pseudospectral method. The technique requires a special treatment of the boundary conditions, which is based on 1-D characteristics normal to the boundaries. The numerical solution of Lamb’s problem requires two 1-D stretching transformations for each Cartesian direction. The results show excellent agreement between the elastic numerical and analytical solutions, demonstrating the effectiveness of the differential operator and boundary treatment. Another example, requiring 1-D transformations, tests the propagation of a Rayleigh wave around a corner of the numerical mesh. Two‐dimensional transformations adapt the grid to topographic features: a syncline, and an anticlinal structure formed with fine layers.

Efficient modeling of wave propagation in a vertical transversely isotropic attenuative medium based on fractional Laplacian

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. T121-T131 ◽  
Author(s):  
Tieyuan Zhu ◽  
Tong Bai
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Fractional Laplacian ◽  
Numerical Solutions ◽  
Pseudospectral Method ◽  
Transversely Isotropic ◽  
Seismic Attenuation ◽  
Fourier Pseudospectral Method ◽  
Frequency Independent ◽  
Theoretical Predictions

To efficiently simulate wave propagation in a vertical transversely isotropic (VTI) attenuative medium, we have developed a viscoelastic VTI wave equation based on fractional Laplacian operators under the assumption of weak attenuation ([Formula: see text]), where the frequency-independent [Formula: see text] model is used to mathematically represent seismic attenuation. These operators that are nonlocal in space can be efficiently computed using the Fourier pseudospectral method. We evaluated the accuracy of numerical solutions in a homogeneous transversely isotropic medium by comparing with theoretical predictions and numerical solutions by an existing viscoelastic-anisotropic wave equation based on fractional time derivatives. To accurately handle heterogeneous [Formula: see text], we select several [Formula: see text] values to compute their corresponding fractional Laplacians in the wavenumber domain and interpolate other fractional Laplacians in space. We determined its feasibility by modeling wave propagation in a 2D heterogeneous attenuative VTI medium. We concluded that the new wave equation is able to improve the efficiency of wave simulation in viscoelastic-VTI media by several orders and still maintain accuracy.

Staggered mesh for the anisotropic and viscoelastic wave equation

Geophysics ◽  
1999 ◽  
Vol 64 (6) ◽  
pp. 1863-1866 ◽  
Author(s):  
José M. Carcione
Keyword(s):  
Wave Equation ◽  
Differential Operators ◽  
Dynamic Range ◽  
Pseudospectral Method ◽  
Transversely Isotropic ◽  
Staggered Grid ◽  
Physical Situation ◽  
Fourier Pseudospectral Method ◽  
Regular Grids ◽  
Grid Solution

Computation of the spatial derivatives with nonlocal differential operators, such as the Fourier pseudospectral method, may cause strong numerical artifacts in the form of noncausal ringing. This situation happens when regular grids are used. The problem is attacked by using a staggered pseudospectral technique, with a different scheme for each rheological relation. The nature and causes of acausal ringing in regular grid methods and the reasons why staggered‐grid methods eliminate this problem are explained in papers by Fornberg (1990) and Özdenvar and McMechan (1996). Thus, the objective here is not to propose a new method but to develop the algorithm for the viscoelastic and transversely isotropic (VTI) wave equation, for which the technique can be implemented without interpolation. The algorithm is illustrated for one physical situation that requires very high accuracy, such as a fluid‐solid interface, where very large contrasts in material properties occur. The staggered‐grid solution is noise free in the dynamic range where regular grids generate artifacts that may have amplitudes similar to those of physical arrivals.

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.

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.

Sign in / Sign up

Export Citation Format

Share Document