A perfectly matched layer for the absorption of electromagnetic waves
A new technique of free-space simulation has been developed for solving unbounded electromagnetic problems with the finite-difference time-domain method. Referred to as PML, the new technique is based on the use of an absorbing layer especially designed to absorb without reflection the electromagnetic waves. The first part of the paper presents the theory of the PML technique. The second part is devoted to numerical experiments and to numerical comparisons with the previously used techniques of free-space simulation. These comparisons show that the PML technique works better than the others in all cases; using it allows us to obtain a higher accuracy in some problems and a release of computational requirements in some others.
References (12)
- J. Blaschak et al.
J. Comput. Phys.
(1988) - A. Taflove
Wave Motion
(1988) - K.S. Yee
IEEE Trans. Antennas Propag.
(1966) - D. Merewether
IEEE Trans. Electromagn. Compat.
(1971) - R. Holland
IEEE Trans. Nucl. Sci.
(1977) - J.P. Berenger, Note Technique DGA/ETCA/DET/390, 1977...
Cited by (9301)
Exact domain truncation for the Morse-ingard equations
2024, Computers and Mathematics with ApplicationsMorse and Ingard [1] give a coupled system of time-harmonic equations for the temperature and pressure of an excited gas. These equations form a critical aspect of modeling trace gas sensors. Like other wave propagation problems, the computational problem must be closed with suitable far-field boundary conditions. Working in a scattered-field formulation, we adapt a nonlocal boundary condition proposed in [2] for the Helmholtz equation to this coupled system. This boundary condition uses a Green's formula for the true solution on the boundary, giving rise to a nonlocal perturbation of standard transmission boundary conditions. However, the boundary condition is exact and so Galerkin discretization of the resulting problem converges to the restriction of the exact solution to the computational domain. Numerical results demonstrate that accuracy can be obtained on relatively coarse meshes on small computational domains, and the resulting algebraic systems may be solved by GMRES using the local part of the operator as an effective preconditioner. These numerical results taken together combine several advanced techniques, including higher-order finite elements, geometric multigrid in curvilinear geometry, native use of complex arithmetic, and incorporation of nonlocal operators. These are tied together in a high-level simulation using the Firedrake library.
On the representation of reinforced concrete structures in FDTD-based simulations to compute lightning-induced magnetic fields
2024, Electric Power Systems ResearchReinforced concrete structures are generally well protected against the adverse effects of a direct lightning strike. Nevertheless, when it occurs, the electromagnetic field generated by the lighting current flowing through their reinforcement can cause a malfunction of sensitive electronic devices. The electromagnetic field can be estimated using full-wave methods, although the accuracy of the estimation depends on the representatives of the model, which is why experimental validation is essential. This paper presents the numerical model of a reinforced concrete structure where a transient current was injected to emulate a direct lightning strike. Full-wave simulations are conducted using TEMSI-FD and the results are compared to the measurements made at the testing facility of EDF TEGG. The simulation results were in good agreement with the measurements, especially when the soil was used as the return path for the current.
Analytic pulse technique for computational electromagnetics
2024, Computer Physics CommunicationsNumerical modeling of electromagnetic waves is an important tool for understanding the interaction of light and matter, and lies at the core of computational electromagnetics. Traditional approaches to injecting and evolving electromagnetic waves, however, can be prohibitively expensive and complex for emerging problems of interest and can restrict the comparisons that can be made between simulation and theory. As an alternative, we demonstrate that electromagnetic waves can be incorporated analytically by decomposing the physics equations into analytic and computational parts. In particle-in-cell simulation of laser–plasma interaction, for example, treating the laser pulse analytically enables direct examination of the validity of approximate solutions to Maxwell's equations including Laguerre–Gaussian beams, allows lower-dimensional simulations to capture 3-D–like focusing, and facilitates the modeling of novel space–time structured laser pulses such as the flying focus. The flexibility and new routes to computational savings introduced by this analytic pulse technique are expected to enable new simulation directions and significantly reduce computational cost in existing areas.
Semi-analytical solution for sound propagation from a moving directional source in a shallow-water waveguide
2024, Journal of Sound and VibrationUnderstanding the waveguide fields produced by moving directional sources is of crucial importance in underwater acoustic detection. In practice, moving underwater targets often exhibit directionality, especially at high frequencies. Herein, we present a semi-analytical method to calculate the Doppler-shifted field (DSF) from a directional source moving horizontally in a range-independent, shallow-water waveguide. We develop a novel normal-mode model using a modal-projection approach with a perfectly matched layer. This comprehensively accounts for changes in modal shape, eigenvalue shifts, and the branch-integral contribution. We derive an analytical expression for the DSF based on Huygens’ principle, and the solution is presented as a sum over propagating eigenmodes, with the modal excitation determined by the source directivity depending upon the grazing angle of each pair of modal plane waves. This expression is valid for any type of radiator moving in a range-independent environment with a smoothly varying sound-speed profile. We then present a theoretical analysis of the Doppler beam shift produced by a piston-like radiator. This reveals that the observed angular offset of modal plane waves can be attributed to the Doppler beam shift and confirms that this is a new mechanism causing considerable differences between the fields generated by a directional moving source and by the same source at rest. This has potential value in various applications, particularly for underwater acoustic detection of moving targets.
TLM-CFSPML for 3D dynamic responses of a layered transversely isotropic half-space
2024, Computers and GeotechnicsThe perfectly matched layer (PML) has proven to be efficient for absorbing elastic waves in the unbounded domain. This strategy is frequently coupled with the semi-analytical approach, such as the thin layer method (TLM), to simulate wave propagation in a layered half-space subjected to dynamic sources. However, beyond spurious reflections due to discretization and zero-frequency singularity, traditional PML may cause exponentially increasing backpropagation waves attributed to material anisotropy when simulating transversely isotropic (TI) media, resulting in significant computational errors. This paper explores the coupling of the complex frequency shifted perfectly matched layer (CFSPML) with the TLM to provide a highly efficient and accurate tool for modelling wave propagation in a layered TI medium. The 3D fundamental solution for a dynamic point source acting on a multi-layered TI half-space is therefore derived. Numerical cases are conducted to evaluate the performance of the TLM-CFSPML in comparison to the TLM-PML. Various factors such as the material anisotropy, stiffness, load and observation positions, and loading frequency are systematically considered. The numerical results demonstrate that the CFSPML can address the computational instability issues related to the anisotropy of the medium. The TLM-CFSPML exhibits computational stability with all the considered scenarios while maintaining high computational efficiency.
Time domain modeling of elastic waves using a stress-based unsplit-field perfectly matched layer with enhanced numerical stability
2024, Applied Mathematical ModellingA stress-based unsplit-field perfectly matched layer (PML) method with enhanced numerical stability is proposed to simulate elastic waves in heterogeneous unbounded domains. The proposed PML method attenuates the waves by forcing the satisfaction of mixed-field equations constructed using a complex coordinate stretching in the frequency domain. Unlike other unsplit-field PMLs, which introduce coordinate stretching before combining the constitutive and compatibility equations, the proposed method first combines the two equations and then stretches the resulting equation. The stretched equation is converted back to the time domain using the inverse Fourier transform to yield transient elastic wave equations in the PML-truncated domain. Semi-discrete equations of motion for displacement and stresses were obtained using standard Lagrange-family finite elements. In a series of numerical examples involving semi-infinite elastic media, the accuracy of the wave solutions calculated using the proposed PML method was evaluated using a displacement -norm, a normalized -error, and an energy norm. The examples further validated the wave absorption performance of the PML and evaluated the long-term stability of solutions utilizing the Ricker pulse excitation with various central frequencies. The numerical experiments indicate that the solutions obtained using the proposed PML are significantly more stable than those obtained using the existing strain-based unsplit-field PMLs. The proposed method is particularly useful for applications requiring explicit time integration of semi-discrete equations of motion because the mass matrix can be diagonalized. The PML method and its associated solution strategies can be applied to various engineering problems, such as soil-structure interaction analysis, monitoring of structures, and geophysical probing.