Domain Electric Field Integral Equation Research Articles

Finite difference delay modelling for analysis of shielding effectiveness of perforated conductive enclosure

AbstractThis paper introduces finite difference delay modelling (FDDM) for computing the shielding effectiveness (SE) of a conductive enclosure with an aperture against electromagnetic pulses. FDDM offers optimal accuracy and stability for analysing complex structures. The time domain electric field integral equation (TD‐EFIE) for the conductive enclosure is derived by imposing boundary conditions on the perfect electrical conductor (PEC) surface in the Laplace domain. Time discretization is based on finite differences, and the Laplace‐to‐Z transform mapping is utilized. Fast Fourier Transform (FFT) is employed to expedite the FDDM solution process. Both frequency domain shielding effectiveness (FD‐SE) and time domain shielding effectiveness (TD‐SE) are evaluated using this approach. Finally, to validate the accuracy of the proposed method, results are compared with simulations using CST‐MWS software and the frequency domain moment method.

Applying Loop-Star Functions With Multibranch Rao-Wilton-Glisson Basis Functions in the Time Domain Electric Field Integral Equation

In this paper, novel loop-star basis functions with the multibranch Rao-Wilton-Glisson (MB-RWG) basis functions are applied to obtain the low-frequency stable solution of the time domain electric field integral equation (TD-EFIE). With the help of the MB-RWG basis functions, it is possible to realize the quasi-Helmholtz decomposition in time domain between different regions that are discretized separately with nonconformal meshes, by modifying slightly the generation of the loop and star bases at the interfaces. Especially, the temporally differentiated form of the TD-EFIE is tested by the star basis, while the undifferentiated form of the TD-EFIE is tested by the loop basis. By means of the filtering technique, the components of the solution which have the potential of growing exponentially have been treated. Numerical results illustrate the accuracy and robustness of the proposed method in analyzing the transient scattering problems of perfectly electrical conducting (PEC) objects.

A novel space hybrid meshing technique combined with Laguerre polynomial to solve TD-EFIE

In this paper, we have developed a system of Time Domain Electric Field Integral Equation (TD-EFIE) for the resolution of a scattering electromagnetic problem by various shapes illuminated by an incident transient wave. The solvability of the TD-EFIE system is based on the Time Domain Method of Moments (TD-MoM) using Galerkin's method in space and time domains. In time domain, a linear combination of Laguerre functions is used as temporal basis functions. In order to approximate the unknown space variation of the transient current, we develop a novel space hybrid basis functions derived from the Rao-Wilton-Glisson function. The elaborated hybrid mesh technique using hexagonal and triangular mesh is applied for various analyzed structures: regular grid for planar structures and irregular grid for curved structures. The validity, the efficiency and the precision of developed meshing techniques are presented and discussed.

Large Time Step and DC Stable TD-EFIE Discretized With Implicit Runge–Kutta Methods

The Time Domain-Electric Field Integral Equation (TD-EFIE) and its differentiated version are widely used to simulate the transient scattering of a time dependent electromagnetic field by a Perfect Electrical Conductor (PEC). The time discretization of the TD-EFIE can be achieved by a space-time Galerkin approach or, as it is considered in this contribution, by a convolution quadrature using Implicit Runge-Kutta methods. The solution is then computed using the Marching-On-in-Time (MOT) algorithm. The differentiated TD-EFIE has two problems: (i) the system matrix suffers from ill-conditioning when the time step increases (low frequency breakdown) and (ii) it suffers from the DC instability, i.e. the formulation allows for the existence of spurious solenoidal currents that grow slowly in the solution. In this work, we show that (i) and (ii) can be alleviated by leveraging quasi-Helmholtz projectors to separate the Helmholtz components of the induced current and rescale them independently. The efficacy of the approach is demonstrated by numerical examples including benchmarks and real life applications.

A leader‐return‐stroke consistent macroscopic model for calculations of return stroke current and its optical and electromagnetic emissions

AbstractA downward lightning flash usually starts with a downward leader and an upward connecting leader followed by an upward return stroke. It is the preceding leader that governs the following return stroke property. Besides, the return stroke property evolves with height and time. These two aspects, however, are not well addressed in most existing return stroke models. In this paper, we present a leader‐return stroke consistent model based on the time domain electric field integral equation, which is a growth and modification of Kumar's macroscopic model. The model is further extended to simulate the optical and electromagnetic emissions of a return stroke by introducing a set of equations relating the return stroke current and conductance to the optical and electromagnetic emissions. With a presumed leader initiation potential, the model can then simulate the temporal and spatial evolution of the current, charge transfer, channel size, and conductance of the return stroke, furthermore the optical and electromagnetic emissions. The model is tested with different leader initiation potentials ranging from −10 to −140 MV, resulting in different return stroke current peaks ranging from 2.6 to 209 kA with different return stroke speed peaks ranging from 0.2 to 0.8 speed of light and different optical power peaks ranging from 4.76 to 248 MW/m. The larger of the leader initiation potential, the larger of the return stroke current and speed. Both the return stroke current and speed attenuate exponentially as it propagates upward. All these results are qualitatively consistent with those reported in the literature.

New mapping properties of the time domain electric field integral equation

We show some improved mapping properties of the Time Domain Electric Field Integral Equation and of its Galerkin semidiscretization in space. We relate the weak distributional framework with a stronger class of solutions using a group of strongly continuous operators. The stability and error estimates we derive are sharper than those in the literature.

A DC-Stable, Well-Balanced, Calderón Preconditioned Time Domain Electric Field Integral Equation

The marching-on-in-time (MOT) solution of the time domain electric field integral equation (TD-EFIE) has traditionally suffered from a number of problems, including: 1) instability; 2) spurious static contributions plaguing the solution; 3) low-frequency breakdown; and 4) dense discretization breakdown. The first issue can be resolved by employing proper space-time Galerkin discretization schemes and accurate quadrature methods. The second and the third issue have been resolved by the quasi-Helmholtz Projected TD-EFIE (qHP-TDEFIE). This contribution introduces a multiplicative preconditioner which can be applied to the qHP-TDEFIE, without further modifying the original scheme. This preconditioner is based on Calderon techniques and guarantees that the MOT system can be solved efficiently using iterative methods, not only for large time step sizes but also for dense spatial discretizations, and for both simply and multiply connected geometries.

Marching-on in degree method for electromagnetic coupling analysis of carbon nanotubes (CNT) dipoles array

Transient electromagnetic analysis of carbon nanotubes (CNT) dipoles is presented. The thin wire approximation and the scattering theory are employed to formulate the problem. Therefore, an accurate system of coupled time domain electric field integral equations (TD-EFIE) is developed in order to predict the responses of the CNT dipoles excited by transient electromagnetic waves. To obtain stable and accurate solutions based on proposed system, the marching-on in degree method is applied with piecewise functions as space basis and weighted Laguerre Polynomials as temporal basis. Numerical results involving transient currents on each CNT dipole are presented. The electromagnetic coupling analysis between the CNT dipoles is conducted depending on the separate distances and the number of the CNT dipoles. Hence, the physical coupling effects are proved.

Time domain simulations of preliminary breakdown pulses in natural lightning.

Lightning discharge is a complicated process with relevant physical scales spanning many orders of magnitude. In an effort to understand the electrodynamics of lightning and connect physical properties of the channel to observed behavior, we construct a simulation of charge and current flow on a narrow conducting channel embedded in three-dimensional space with the time domain electric field integral equation, the method of moments, and the thin-wire approximation. The method includes approximate treatment of resistance evolution due to lightning channel heating and the corona sheath of charge surrounding the lightning channel. Focusing our attention on preliminary breakdown in natural lightning by simulating stepwise channel extension with a simplified geometry, our simulation reproduces the broad features observed in data collected with the Huntsville Alabama Marx Meter Array. Some deviations in pulse shape details are evident, suggesting future work focusing on the detailed properties of the stepping mechanism.Key PointsPreliminary breakdown pulses can be reproduced by simulated channel extension Channel heating and corona sheath formation are crucial to proper pulse shape Extension processes and channel orientation significantly affect observations

A Higher Order Space-Time Galerkin Scheme for Time Domain Integral Equations

Stability of time domain integral equation (TDIE) solvers has remained an elusive goal for many years. Advancement of this research has largely progressed on four fronts: (1) Exact integration, (2) Lubich quadrature, (3) smooth temporal basis functions, and (4) Space-time separation of convolutions with the retarded potential. The latter method was explored in [Pray et al. IEEE TAP 2012]. This method's efficacy in stabilizing solutions to the time domain electric field integral equation (TD-EFIE) was demonstrated on first order surface descriptions (flat elements) in tandem with 0th order functions as the temporal basis. In this work, we develop the methodology necessary to extend to higher order surface descriptions as well as to enable its use with higher order temporal basis functions. These higher order temporal basis functions are used in a Galerkin framework. A number of results that demonstrate convergence, stability, and applicability are presented.

Time Domain Analysis of Electromagnetic Scattering Problems by Using Integral Equation Method With Space-Delayed Temporal Basis Functions

This communication presents a new numerical method for the efficient solution of marching on in time (MOT)-based time domain integral equations (TDIE) for electromagnetic scattering from perfect electrically conducting (PEC) object with large smooth surfaces. A set of space-delayed temporal basis functions is proposed for the temporal discretization, and the curvilinear Rao-Wilson-Glisson (CRWG) basis functions for the spatial discretization. The space-delayed temporal basis functions are firstly extended to the time domain electric field integral equations (TDEFIE), then easily to the time domain combined field integral equations (TDCFIE). Numerical results demonstrate that the new method is stable and effective for the transient scattering problems from PEC objects with large smooth surfaces and a single plane wave incidence. A significant reduction of the number of spatial unknowns is achieved by using the space-delayed temporal basis functions.

Corrections to “A stable solution of time domain electric field integral equation using weighted laguerre polynomials”

ABSTRACTCorrections to the final equations in (Chung et al., Microwave Opt Technol Lett 49 (2007), 2789–2793) are presented. © 2014 Wiley Periodicals, Inc. Microwave Opt Technol Lett 56:1263–1264, 2014

TRANSIENT ANALYSIS OF A RECTANGULAR CAVITY CONTAINING AN INTERIOR SCATTERER USING TD-EFIE WITH WEIGHTED LAGUERRE POLYNOMIALS AS TEMPORAL BASIS FUNCTIONS

Novel 2-D Time Domain Electric Field Integral Equations (TD-EFIE) are established in order to predict transient response of a wire enclosed within a rectangular cavity. The wire and cavity are excited by an external incident transient electromagnetic wave through a slot in the cavity wall. The formulation of the TD-EFIE is based on equivalence principle and boundary conditions taking account the effect of reflection from cavity walls. The equations are efficiently solved by Method of Moments. The transient unknown coefficients of the electric current at the wire and magnetic current at the slot are approximated using a set of orthonormal temporal basis functions derived from Laguerre Polynomials. The analysis demonstration is presented to prove that the novel TD-EFIE combined to MoM is able to solve this critical problem. No late-time instability is encountered.

On the Internal Resonant Modes in Marching-on-in-Time Solution of the Time Domain Electric Field Integral Equation

Internal resonant modes are always observed in the marching-on-in-time (MOT) solution of the time domain electric field integral equation (EFIE), although “relaxed initial conditions,” which are enforced at the beginning of time marching, should in theory prevent these spurious modes from appearing. It has been conjectured that, numerical errors built up during time marching establish the necessary initial conditions and induce the internal resonant modes. However, this conjecture has never been proved by systematic numerical experiments. Our numerical results in this communication demonstrate that, the internal resonant modes' amplitudes are indeed dictated by the numerical errors. Additionally, it is shown that in a few cases, the internal resonant modes can be made “invisible” by significantly suppressing the numerical errors. These tests prove the conjecture that the internal resonant modes are induced by numerical errors when the time domain EFIE is solved by the MOT method.

High-Order Calderón Preconditioned Time Domain Integral Equation Solvers

Two high-order accurate Calderon preconditioned time domain electric field integral equation (TDEFIE) solvers are presented. In contrast to existing Calderon preconditioned time domain solvers, the proposed preconditioner allows for high-order surface representations and current expansions by using a novel set of fully-localized high-order div- and quasi curl-conforming (DQCC) basis functions. Numerical results demonstrate that the linear systems of equations obtained using the proposed basis functions converge rapidly, regardless of the mesh density and of the order of the current expansion.

A Space-Time Mixed Galerkin Marching-on-in-Time Scheme for the Time-Domain Combined Field Integral Equation

The time domain combined field integral equation (TD-CFIE), which is constructed from a weighted sum of the time domain electric and magnetic field integral equations (TD-EFIE and TD-MFIE) for analyzing transient scattering from closed perfect electrically conducting bodies, is free from spurious resonances. The standard marching-on-in-time technique for discretizing the TD-CFIE uses Galerkin and collocation schemes in space and time, respectively. Unfortunately, the standard scheme is theoretically not well understood: stability and convergence have been proven for only one class of space-time Galerkin discretizations. Moreover, existing discretization schemes are nonconforming, i.e., the TD-MFIE contribution is tested with divergence conforming functions instead of curl conforming functions. We therefore introduce a novel space-time mixed Galerkin discretization for the TD-CFIE. A family of temporal basis and testing functions with arbitrary order is introduced. It is explained how the corresponding interactions can be computed efficiently by existing collocation-in-time codes. The spatial mixed discretization is made fully conforming and consistent by leveraging both Rao-Wilton-Glisson and Buffa-Christiansen basis functions and by applying the appropriate bi-orthogonalization procedures. The combination of both techniques is essential when high accuracy over a broad frequency band is required.

Stability Properties of the Time Domain Electric Field Integral Equation Using a Separable Approximation for the Convolution With the Retarded Potential

The state of art of time domain integral equation (TDIE) solvers has grown by leaps and bounds over the past decade. During this time, advances have been made in (i) the development of accelerators that can be retrofitted with these solvers and (ii) understanding the stability properties of the electric field integral equation. As is well known, time domain electric field integral equation solvers have been notoriously difficult to stabilize. Research into methods for understanding and prescribing remedies have been on the uptick. The most recent of these efforts are (i) Lubich quadrature and (ii) exact integration. In this paper, we re-examine the solution to this equation using (i) the undifferentiated form of the TD-EFIE and (ii) a separable approximation to the spatio-temporal convolution. The proposed scheme can be constructed such that the spatial integrand over the source and observer domains is smooth and integrable. As several numerical results will demonstrate, the proposed scheme yields stable results for long simulation times and a variety of targets, both of which have proven extremely challenging in the past.

Analysis of Convolution Quadrature Applied to the Time-Domain Electric Field Integral Equation

AbstractWe show how to apply convolution quadrature (CQ) to approximate the time domain electric field integral equation (EFIE) for electromagnetic scattering. By a suitable choice of CQ, we prove that the method is unconditionally stable and has the optimal order of convergence. Surprisingly, the resulting semi discrete EFIE is dispersive and dissipative, and we analyze this phenomena. Finally, we present numerical results supporting and extending our convergence analysis.

An Improved Marching-on-in-Degree Method Using a New Temporal Basis

The marching-on-in-degree (MOD) method has been presented earlier for solving time domain electric field integral equations in a stable fashion. This is accomplished by expanding the transient responses by a complete set of orthogonal entire domain associated Laguerre functions, which helps one to analytically integrate out the time variable from the final computations in a Galerkin methodology. So, the final computations are carried out using only the spatial variables. However, the existing MOD method suffers from low computational efficiency over a marching-on-in-time (MOT) method. The two main causes of the computational inefficiency in the previous MOD method are now addressed using a new form of the temporal basis functions and through a different computational arrangement for the Green's function. In this paper, it is shown that incorporating these two new concepts can speed up the computational process and make it comparable to a MOT algorithm. Sample numerical results are presented to illustrate the validity of these claims in solution of large problems using the MOD method.

Time Domain CalderÓn Identities and Their Application to the Integral Equation Analysis of Scattering by PEC Objects Part I: Preconditioning

Time domain electric field integral equations often are used to analyze transient scattering from perfect electrically conducting objects. When discretized using marching-on-in-time recipes they give rise to linear systems of equations that can be solved for the induced currents for all time steps. Unfortunately, when the scatterer is approximated by increasingly dense meshes, the condition number of these systems grows rapidly, slowing down the convergence of iterative solvers. Here, time domain Calderon identities are derived and subsequently used to construct a Calderon-preconditioned time domain electric field integral equation that can be discretized even with dense meshes using Buffa-Christiansen basis functions. Numerical results that demonstrate the effectiveness and accuracy of the proposed method are presented.