Numerical solution of time domain integral equations for arbitrarily shaped conductor/dielectric composite bodies

In this work, we present a numerical solution of the coupled time domain integral equations to obtain induced currents and scattered far fields on a three-dimensional, arbitrary shaped conducting/dielectric composite body illuminated by a Gaussian electromagnetic plane wave pulse. The coupled integral equations are derived utilizing the equivalence principle. The solution method is based on the method of moments and involves the triangular patch modeling of the composite body, in conjunction with the patch basis functions. Detailed mathematical steps along with several numerical results are presented to illustrate the efficacy of this approach.     

Volume integral equations for analysis of dielectric branchingwaveguides

Novel forms of volume integral equations are developed for the exact treatment of wave propagation in two-dimensional dielectric branching waveguides. The integral equations can be obtained by considering the condition at a point far away from the junction section. An approximate solution by the Born approximation and a numerical solution by the moment method establish the validity of the new volume integral equations. The numerical results are discussed from the viewpoint of energy conservation and reciprocity. The solution is exact if sufficiently large computer memory and computational time are used. The method can be extended to problems of a more general nature (i.e. the incident TM mode), and complex configurations of branching waveguides. The basic idea is also applicable to techniques using boundary (surface) integral equations which are applicable to three-dimensional problems     

透层天线电流积分方程的数值解

在已有的三层有耗媒质中的透层天线（ＢＰＡ）电流积分方程的基础上,给出了电流积分方程中的索末菲积分的低频近似解析公式并给出了其物理解释。利用三项式全域电流基函数,求得了透层天线的输入阻抗。采用静态场近似方法,很方便的求出了透层天线在地面的电场分布。     

Single integral equation for electromagnetic scattering bythree-dimensional homogeneous dielectric objects

A single integral equation formulation for electromagnetic scattering by three-dimensional (3-D) homogeneous dielectric objects is developed. In this formulation, a single effective electric current on the surface S of a dielectric object is used to generate the scattered fields in the interior region. The equivalent electric and magnetic currents for the exterior region are obtained by enforcing the continuity of the tangential fields across S. A single integral equation for the effective electric current is obtained by enforcing the vanishing of the total field due to the exterior equivalent currents inside S. The single integral equation is solved by the method of moments. Numerical results for a dielectric sphere obtained with this method are in good agreement with the exact results. Furthermore, the convergence speed of the iterative solution of the matrix equation in this formulation is significantly greater than that of the coupled integral equations formulation     

Iterative integral equation solution for cylindrical dielectric resonator

An exact integral equation describing the substrate-mounted open cylindrical dielectric resonator is derived, and a first approximate solution for the fundamental resonant mode is calculated. Results are presented and are compared with both experiment and another theoretical approach.     

Volterra integral equations for nonstationary electromagnetic processes in time-varying dielectric waveguides

Time-domain Volterra integral equations are obtained for the electromagnetic fields in longitudinal uniform dielectric waveguides with time-varying media in their core. Resolvent operators are constructed for the case when the core permittivity changes abruptly in time and new features are observed in the consequent field transients in comparison with those occurring in a simple unbounded space.     

Fast solution of mixed-potential time-domain integral equations for half-space environments

A fast Fourier transform-accelerated integral-equation based algorithm to efficiently analyze transient scattering from planar perfect electrically conducting objects residing above or inside a potentially lossy dielectric half-space is presented. The algorithm requires O(N/sub t/N/sub s/(logN/sub s/+log/sup 2/N/sub t/)) CPU and O(N/sub t/N/sub s/) memory resources when analyzing electromagnetic wave interactions with uniformly meshed planar structures. Here, N/sub t/ and N/sub s/ are the numbers of simulation time steps and spatial unknowns, respectively. The proposed scheme is therefore far more efficient than classical time-marching solvers, the CPU and memory requirements of which scale as O(N/sub t//sup 2/N/sub s//sup 2/) and O(N/sub t/N/sub s//sup 2/). In the proposed scheme, all pertinent time-domain half-space Green functions are (pre) computed from their frequency-domain counterparts via inverse discrete Fourier transformation. In this process, in-band aliasing is avoided through the application of a smooth and interpolatory window. Numerical results demonstrate the accuracy and efficiency of the proposed algorithm.     

Improved quadrature formulas for boundary integral equations with conducting or dielectric edge singularities

In this paper we derive new two-dimensional (2-D) quadrature formulas for the discretization of boundary integral equations in the presence of conducting or dielectric edges. The proposed formulas allow us to exactly integrate polynomials of degree less than or equal to five, multiplied by an algebraic singular factor that diverges along one side of the triangular integration domain. This is the kind of singularity that occurs when physical edges are present in both conducting and dielectric bodies. Numerical tests are performed on the presented formulas, in order to validate the achieved improvement in accuracy, and examples are given of their application to the determination of radar cross-section of 3-D metallic objects.     

A fast combined field volume integral equation solution to EM scattering by 3-D dielectric objects of arbitrary permittivity and permeability

A fast solution to the combined field volume integral equation (CFVIE) for electromagnetic scattering by large three-dimensional dielectric bodies of arbitrary permittivity and permeability is presented. The CFVIE is formulated in the region of the scatterers by expressing the total fields as the sum of the incident wave and the radiated wave due to both the electric and magnetic polarization currents. The resultant integral equation is solved using the method of moments (MoM). Then the precorrected fast Fourier transform (P-FFT) method is applied to reduce the memory requirement and accelerate the matrix-vector multiplication in the MoM solution. In the implementation of the P-FFT method, two sets of projection operators are constructed respectively for the projections of the electric sources and magnetic sources. In addition, two sets of interpolation operators are also applied respectively for the computation of the vector/scalar potentials and the curl of the vector potentials in the support of the testing functions. The resultant method has a memory requirement of O(N) and a computational complexity of O(NlogN) respectively, where N denotes the number of unknowns.     

A technique for solution of Maxwell's equations in a moving dielectric medium

A simple technique is presented for converting a known solution for the electric and magnetic vector fields in a dielectric medium at rest into the corresponding fields in a moving dielectric medium. The technique combines methods presented by Tai [1] with a scaling procedure developed by Clemmow [2]. Tai's work reduces the moving medium problem to the solution of Maxwell's equations in a uniaxial medium, and Clemmow's procedure enables one to convert a known solution in an isotropic medium to the corresponding solution in a uniaxial medium. Thus by first solving for the fields in the medium at rest, then following Clemmow's procedure to obtain the fields in Tai's uniaxial medium, and finally applying Tai's reasoning, one may easily obtain the solution of Maxwell's equations in the moving medium.     

Electromagnetic scattering from 3-D curved dielectric bodies usingtime-domain integral equations

A boundary integral equation (BIE) approach is developed to calculate transient scattering from dielectric bodies. The treatment is directly in terms of the E and H fields rather than magnetic and electric currents. It employs curvilinear (quadratic) modeling, which allows accurate representation of arbitrarily shaped curved bodies. The treatment is isoparametric with the same quadratic representation of the spatial field variation and with the temporal variation modeled by similar quadratic elements. Integration employs high-order Gaussian quadrature with special treatment of the singular and hypersingular integrals that arise. The treatment is implicit, requiring the solution of a sparse matrix equation at each timestep. This adds only trivially to the cost at each timestep and, by freeing the timestep from the constraint that it be smaller than the smallest nodal spacing, can greatly reduce the number of timesteps that must be employed. Additionally, it produces stable results without resort to the averaging processes proposed elsewhere. Example calculations of scattering from a sphere, a cube, and an almond are presented and compared with earlier published transient results and with results from a frequency domain treatment. Good agreement and improved accuracy is found     

Fast convergent integral equation solution of strip gratings on dielectric substrate

A fast convergent integral equation solution to the scattering problem of a transverse electric/transverse magnetic (TE/TM) plane wave by a one-dimensional periodic array of thin metal strips on a dielectric substrate is described. The formulation of the integral equation is similar to that derived by Montgomery for a two-dimensional periodic array of thin conductors on a dielectric substrate. However, the basis functions which satisfy the appropriate edge conditions are incorporated here for the unknown current expansion on the strips. Following the standard Galerkin's procedure, one may readily determine the induced currents on the strips and thus the reflected and transmitted fields. Sample numerical results are given and good agreement with previously published data is obtained. It is found that the convergence rate of this method is improved by an order of magnitude. Also it is shown that the dielectric substrate has a strong effect on the scattering from the large spacing strip grating.     

Fast solution of electromagnetic integral equations using adaptivewavelet packet transform

The adaptive wavelet packet transform is applied to sparsify the moment matrices for the fast solution of electromagnetic integral equations. In the algorithm, a cost function is employed to adaptively select the optimal wavelet packet expansion/testing functions to achieve the maximum sparsity possible in the resulting transformed system. The search for the best wavelet packet basis and the moment matrix transformation are implemented by repeated two-channel filtering of the original moment matrix with a pair of quadrature filters. It is found that the sparsified matrix has above-threshold elements that grow only as O(N^1.4) for typical scatterers. Consequently the operations to solve the transformed moment equation using the conjugate gradient method scales as O(N^1.4). The additional computational cost for carrying out the adaptive wavelet packet transform is evaluated and discussed     

On the determination of resonant modes of dielectric objects using surface Integral equations

Although surface integral equations have been extensively used for solving the scattering problem of arbitrarily shaped dielectric objects, when applied to the resonance problem, there are still some issues not fully addressed by the literature. In this paper, the method of moments with Rao-Wilton-Glisson basis functions is applied to the electric field integral equation (EFIE) for solving the resonance problem of dielectric objects. The resonant frequency is obtained by searching for the minimum of the reciprocal of the condition number of the impedance matrix in the complex frequency plane, and the modal field distribution is obtained through singular value decomposition (SVD). The determinant of the impedance matrix is not used since it is difficult to find its roots. For the exterior EFIE, the original basis functions are used as testing functions; for the interior EFIE, the basis functions rotated by 90/spl deg/ are used as testing functions. To obtain an accurate modal field solution, the impedance matrix needs to be reduced by half before SVD is applied to it. Numerical results are given and compared with those obtained by using the volume integral equation.     

Extended boundary condition integral equations for perfectly conducting and dielectric bodies: Formulation and uniqueness

The equivalence theorem is used to derive novel generalized boundary condition (GBC) integral equations for the tangential components of the electric and magnetic fields on the interfaces of a finite number of dielectric or conducting scatterers. Closed surface, plane, and line extended boundary conditions (EBC) equivalent to the GBC are introduced. The GBC integral equations can now be replaced by any of these EBC integral equations whose solutions are unique and easy to obtain numerically using the moment method. A perfectly conducting sphere and a dielectric sphere in the electrostatic field of two equal and opposite point charges are presented as simple examples of the general procedure.     

An iterative solution of fredholm integral equations of the first kind

A formal solution is obtained for the Fredholm integral equation of the first kind. The unknown function is represented by a completely general Fourier series, and the Fourier coefficients are obtained by an iterative process. The formulation also yields an easily obtained approximate solution, as well as the estimate of its accuracy.     

The closed-form solution of vector finite element integral equations for three dimensional electromagnetic analysis

In this paper, a set of closed-form formulas for vector Finite Element Method (FEM) to analyze three dimensional electromagnetic problems is presented on the basis of Gaussian quadrature integration scheme. By analyzing the open region problems, the first-order Absorbing Boundary Condition (ABC) is considered as the truncation boundary condition and the equation is carried out in a closed-form. Based on the formulas, the hybrid Expanded Cholesky Method (ECM) and MultiFrontal algorithm (MF) is applied to solve finite element equations. Using the closed-form solution, the electromagnetic field of three dimensional targets can be studied sententiously and accurately. Simulation results show that the presented formulas are successfully and concise, which can be easily used to analyze three dimensional electromagnetic problems.     

Resonant scattering for characterization of axisymmetric dielectric objects

Resonances in the electromagnetic scattering by dielectric objects are investigated both numerically and experimentally with particular emphasis on their direct application to object characterization. Low ka calculations for spheres and an infinite circular cylinder illustrate the modal nature of the resonances and the dependence of the resonance spectra on refractive index and size. Experiments show that the specific wavelengths at which resonances occur in the scattering intensity from a glass fiber can be used to determine its diameter to a high accuracy. The extended boundary condition method is used to calculate the low ka resonances of a prolate spheroid and a finite circulate cylinder with the anticipation that the scattering resonances may be particularly suitable for characterizing randomly oriented nonspherical objects.     

New boundary integral equations for CAD of waveguide circuits:guided-mode extracted integral equations

Novel boundary integral equations which are applicable to the analysis of many kinds of waveguide circuits are presented. The new integral equations can treat the waveguide discontinuity problems like the scattering by the isolated finite-sized metallic objects or cavity problems and do not use normal-mode expansion techniques. They are suitable for the basic theory of CAD software for various waveguide circuits. The two-port and H-plane waveguide discontinuity problems which satisfy the single-mode and two-mode conditions are treated. The case of waveguide corner bend is considered as an example. Numerical examples are shown in order to confirm the validity of the new integral equations