Generalised higher-order FDTD-PML algorithm for enhanced analysis of 3-D waveguiding EMC structures in curvilinear coordinates

Kantartzis, Nikolaos V.

doi:10.1049/ip-map:20030269

Cited by 11 publications

(1 citation statement)

References 27 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…2009). Recently, many researchers have extended the PML from Cartesian coordinates to cylindrical and spherical coordinates, or even more general curvilinear coordinates (Collino and Monk ; Teixeira and Chew ; Zheng and Huang ; Kantartzis ). However, no free surface boundary condition (FSBC) was involved in their studies.…”

Section: Introductionmentioning

confidence: 99%

Application of a perfectly matched layer in seismic wavefield simulation with an irregular free surface

Lan

Chen

Zhang

et al. 2015

Geophysical Prospecting

View full text Add to dashboard Cite

A B S T R A C TRecently, an effective and powerful approach for simulating seismic wave propagation in elastic media with an irregular free surface was proposed. However, in previous studies, researchers used the periodic condition and/or sponge boundary condition to attenuate artificial reflections at boundaries of a computational domain. As demonstrated in many literatures, either the periodic condition or sponge boundary condition is simple but much less effective than the well-known perfectly matched layer boundary condition. In view of this, we intend to introduce a perfectly matched layer to simulate seismic wavefields in unbounded models with an irregular free surface. We first incorporate a perfectly matched layer into wave equations formulated in a frequency domain in Cartesian coordinates. We then transform them back into a time domain through inverse Fourier transformation. Afterwards, we use a boundary-conforming grid and map a rectangular grid onto a curved one, which allows us to transform the equations and free surface boundary conditions from Cartesian coordinates to curvilinear coordinates. As numerical examples show, if free surface boundary conditions are imposed at the top border of a model, then it should also be incorporated into the perfectly matched layer imposed at the top-left and top-right corners of a 2D model where the free surface boundary conditions and perfectly matched layer encounter; otherwise, reflections will occur at the intersections of the free surface and the perfectly matched layer, which is confirmed in this paper. So, by replacing normal second derivatives in wave equations in curvilinear coordinates with free surface boundary conditions, we successfully implement the free surface boundary conditions into the perfectly matched layer at the top-left and top-right corners of a 2D model at the surface. A number of numerical examples show that the perfectly matched layer constructed in this study is effective in simulating wave propagation in unbounded media and the algorithm for implementation of the perfectly matched layer and free surface boundary conditions is stable for long-time wavefield simulation on models with an irregular free surface.

show abstract

Section: Introductionmentioning

confidence: 99%