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 (9307)
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.
Machine learning-based model inference for spectral response of photonic crystals
2024, Micro and NanostructuresPhotonic Crystals (PhCs) are materials with a periodic arrangement of dielectric or metallic components that can manipulate the flow of photons. Conventional techniques for computing the spectrum response of these artificial structures are time-consuming, laborious, and susceptible to human errors. This paper presents a novel approach incorporating Machine Learning (ML) in forming PhC-based periodic structures. The structures are designed from the Two-Dimensional (2D) PhCs having air holes in a dielectric material. These crystalline structures work on the Guided Mode Resonance (GMR) principle and find their use in numerous applications, including optical filters in Near Infrared Range (NIR). However, the conventional methods for analyzing the output spectra of PhCs need to be revised due to the complexities in designing compound optical filters. Therefore, an automated process employing ML is required. As a result, in our work, we form 2D crystalline structures using the Finite-Difference Time Domain (FDTD) method and predict various spectral responses of PhCs via ML-based linear regression. The proposed mathematical models provide efficient results and considerably less simulation efforts and time compared to the traditional manual methods.
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.