An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation

Martin, Roland; Komatitsch, Dimitri

doi:10.1111/j.1365-246x.2009.04278.x

Cited by 155 publications

(86 citation statements)

References 44 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…A very suitable version of the CPML ABC is given in [18]; this is because its implementation in a FDTD code without a PML ABC is simply obtained by replacing the spatial derivative ∂x, with 1/k∂x + Ψ, where Ψ is the so-called convolutional term of the CPML ABC-whose time evolution is the same as the time evolution of the other variables-and k is the stretching function always greater than 1 [5] and given by (13). According to [18], Ψ can be implemented as a recursive sum in time, so we can update the memory variable of the F (H or E) field in the i = x, y direction for each time step n as…”

Section: Cpml Abc Formulationmentioning

confidence: 99%

Analysis of Electromagnetic Propagation from MHz to THz with a Memory-Optimised CPML-FDTD Algorithm

Rodríguez-Sánchez

Couder-Castañeda

Hernández-Gómez

et al. 2018

International Journal of Antennas and Propagation

View full text Add to dashboard Cite

FDTD method opened a fertile research area on the numerical analysis of electromagnetic phenomena under a wide range of media and propagation conditions, providing an extensive analysis of electromagnetic behaviour like propagation, reflection, refraction, and multitrajectory phenomena. In this paper, we present an optimised FDTD-CPML algorithm, focused in saving memory while increasing the performance of the algorithm. We particularly implement FDTD-CPML method at high frequency bands, used in several telecommunications applications as well as in nanoelectromagnetism. We show an analysis of the performance of the algorithm in single and double precision, as well as a stability of the algorithm analysis, from where we conclude that the implemented CPML ABC constitutes a robust choice in terms of precision and accuracy for the high frequencies herein considered. It is important to recall that the CPML ABC parameters provided in this paper are fixed for the tested range of frequencies, that is, from MHz to THz.

show abstract

Section: Cpml Abc Formulationmentioning

confidence: 99%

Analysis of Electromagnetic Propagation from MHz to THz with a Memory-Optimised CPML-FDTD Algorithm

Rodríguez-Sánchez

Couder-Castañeda

Hernández-Gómez

et al. 2018

International Journal of Antennas and Propagation

View full text Add to dashboard Cite

show abstract

“…η ν are the PML decay factors, α ν are point-wise coefficients related to the CFS-PML (Roden & Gedney 2000). The essence of the real parameter α ν is to absorb the evanescent waves and to improve the absorption performance at grazing angles (Bérenger 2002;Komatitsch & Martin 2007;Martin et al 2009). In this study, the parameter κ ν is set as 1, since the effect of these parameters for the absorption is not significant (Bérenger 2002).…”

Section: Implementation Of Cfs-pmlmentioning

confidence: 99%

“…We use an optimal maximum attenuation (η max ) formulation similar as Komatitsch & Martin (2007), Martin et al (2009) and Pan et al (2012) …”

Section: Implementation Of Cfs-pmlmentioning

confidence: 99%

Application of the perfectly matched layer in 3-D marine controlled-source electromagnetic modelling

Han

et al. 2017

Geophysical Journal International

View full text Add to dashboard Cite

SUMMARYIn this study, the complex frequency-shifted perfectly matched layer (CFS-PML) in stretching Cartesian coordinates, is successfully applied to three-dimensional (3D) frequency-domain marine controlled-source electromagnetic (CSEM) field modelling.The Dirichlet boundary, which is usually used within the traditional framework of EM modeling algorithms, assumes the electric or magnetic field values are zero at the boundaries. This requires the boundaries be sufficiently far away from the sources in the area of interest. To mitigate the boundary artifacts, a large modelling area may be necessary even though cell sizes are allowed to grow toward the boundaries due to the diffusion of the electromagnetic wave propagation. Compared with the conventional Dirichlet boundary, the PML boundary is preferred as the modelling area of interest could be restricted to the target region and only a few absorbing layers surrounding can effectively depress the artificial boundary effect without losing the numerical accuracy. Furthermore, for joint inversion of seismic and marine CSEM data, if we used the PML for CSEM field simulation instead of the conventional Dirichlet, the modeling area for these two different geophysical data collected from the same survey area could the same, which is convenient for joint inversion grid Downloaded from https://academic.oup.com/gji/article-abstract/doi/10.1093/gji/ggx382/4157794/Application-of-the-perfectly-matched-layer-in-3D by GEOMAR Bibliothek Helmholtz-Zentrum fuer Ozeanforschung user on 22 September 2017 2 G. Li and B. Li matching. We apply the CFS-PML boundary to 3D marine CSEM modelling by using the staggered finite-difference (SFD) discretization. Numerical test indicates that the modeling algorithm using the CFS-PML also shows good accuracy compared to the Dirichlet. Furthermore, the modeling algorithm using the CFS-PML shows advantages in computational time and memory saving than that using the Dirichlet boundary. For the 3D example in this study, the memory saving using the PML is nearly 42 % and the time saving is around 48% compared to using the Dirichlet.

show abstract

“…In this study, we use an accurate implementation of the standard staggered-grid (SSG) finitedifference scheme (Virieux, 1986;Levander, 1988;Zhang, 2010) with the convolution perfectly matched layer absorbing boundary condition (Komatitsch and Martin, 2007;Martin and Komatitsch, 2009;Zhang and Shen, 2010) to fully model elastic waves in the presence of heterogeneity. The internal interfaces are represented by the so-called effective medium parameters to avoid spurious numerical diffractions caused by sharp material discontinuity due to the spatial grid.…”

Section: Introductionmentioning

confidence: 99%

Numerical modeling of elastic wave scattering by near-surface heterogeneities

Almuhaidib¹,

Toksöz²

2013

SEG Technical Program Expanded Abstracts 2013

View full text Add to dashboard Cite

In land seismic data, scattering from surface and near-surface heterogeneities adds complexity to the recorded signal and masks weak primary reflections. To understand the effects of near-surface heterogeneities on seismic reflections, we simulated seismic-wave scattering from arbitrary-shaped, shallow, subsurface heterogeneities through the use of a perturbation method for elastic waves and finite-difference forward modeling. The near-surface scattered wavefield was modeled by looking at the difference between the calculated incident (i.e., in the absence of scatterers) and the total wavefields. Wave propagation was simulated for several earth models with different nearsurface characteristics to isolate and quantify the influence of scattering on the quality of the seismic signal. The results indicated that the direct surface waves and the upgoing reflections were scattered by the near-surface heterogeneities. The scattering took place from body waves to surface waves and from surface waves to body waves. The scattered waves consisted mostly of body waves scattered to surface waves and were, generally, as large as, or larger than, the reflections. They often obscured weak primary reflections and could severely degrade the image quality. The results indicated that the scattered energy depended strongly on the properties of the shallow scatterers and increased with increasing impedance contrast, increasing size of the scatterers relative to the incident wavelength, decreasing depth of the scatterers, and increasing attenuation factor of the background medium. Also, sources deployed at depth generated weak surface waves, whereas deep receivers recorded weak surface and scattered body-to-surface waves. The analysis and quantified results helped in the understanding of the scattering mechanisms and, therefore, could lead to developing new acquisition and processing techniques to reduce the scattered surface wave and enhance the quality of the seismic image.

show abstract

An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation

Cited by 155 publications

References 44 publications

Analysis of Electromagnetic Propagation from MHz to THz with a Memory-Optimised CPML-FDTD Algorithm

Analysis of Electromagnetic Propagation from MHz to THz with a Memory-Optimised CPML-FDTD Algorithm

Application of the perfectly matched layer in 3-D marine controlled-source electromagnetic modelling

Numerical modeling of elastic wave scattering by near-surface heterogeneities

Contact Info

Product

Resources

About