High order boundary element methods

The thesis is a major contribution to the state of research in the field of boundary element methods in general and to the electromagnetic engineering sciences in particular. The first contribution is the development and implementation of high order boundary element methods. It is explained how to set up a high order boundary element implementation which provides solvers for elliptic, Maxwell and mixed problems. The second contribution is a complete mathematical theory for a stabilized boundary variational formulation describing scattering problems. The stabilized formulation presented in the thesis is equivalent to the classical formulations, however, it does not suffer from the low-frequency break-down. The high order methods have been applied to the classical and stabilized formulations describing electromagnetic scattering and the theoretical results on asymptotic convergence and stability have been verified by the numerical computations.Die Dissertation befaßt sich mit der Entwicklung von Randelementmethoden höherer Ordnung. Die Implementierung, die im Zuge der Arbeit entstand, unterscheidet sich von herkömmlicher Software dadurch, dass sie zur numerischen Lösung von elliptischen Problemen, von Maxwell Problemen und insbesondere zur numerischen Lösung von gemischten Formulierungen benutzt werden kann. In der Arbeit werden die theoretisch bewiesenen Konvergenzraten für all diese Problemklassen verifiziert. Um die Leistungsfähigkeit der Software zu demonstrieren, wird insbesondere eine gemischte Formulierung zur Lösung eines elektromagnetischen Streuproblems entwickelt. Diese sogenannte stabilisierte Formulierung ist gleichsam das zweite Forschungsergebnis dieser Arbeit. Im Unterschied zu der klassischen Formulierung gewährleistet die stabilierte Formulierung numerische Stabilität im Grenzfall quasi-elektrostatischer Prozesse. Die theoretischen und numerischen Resultate, die diese Aussage rechtfertigen, werden in der Arbeit geliefert

Formulation and Solution of Electromagnetic Integral Equations Using Constraint-Based Helmholtz Decompositions

This dissertation develops surface integral equations using constraint-based Helmholtz decompositions for electromagnetic modeling. This new approach is applied to the electric field integral equation (EFIE), and it incorporates a Helmholtz decomposition (HD) of the current. For this reason, the new formulation is referred to as the EFIE-hd. The HD of the current is accomplished herein via appropriate surface integral constraints, and leads to a stable linear system. This strategy provides accurate solutions for the electric and magnetic fields at both high and low frequencies, it allows for the use of a locally corrected Nyström (LCN) discretization method for the resulting formulation, it is compatible with the local global solution framework, and it can be used with non-conformal meshes. To address large-scale and complex electromagnetic problems, an overlapped localizing local-global (OL-LOGOS) factorization is used to factorize the system matrix obtained from an LCN discretization of the augmented EFIE (AEFIE). The OL-LOGOS algorithm provides good asymptotic performance and error control when used with the AEFIE. This application is used to demonstrate the importance of using a well-conditioned formulation to obtain efficient performance from the factorization algorithm

Numerical methods for electromagnetic wave propagation and scattering in complex media

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2004.Vita.Includes bibliographical references (p. 227-242).Numerical methods are developed to study various applications in electromagnetic wave propagation and scattering. Analytical methods are used where possible to enhance the efficiency, accuracy, and applicability of the numerical methods. Electromagnetic induction (EMI) sensing is a popular technique to detect and discriminate buried unexploded ordnance (UXO). Time domain EMI sensing uses a transient primary magnetic field to induce currents within the UXO. These currents induce a secondary field that is measured and used to determine characteristics of the UXO. It is shown that the EMI response is difficult to calculate in early time when the skin depth is small. A new numerical method is developed to obtain an accurate and fast solution of the early time EMI response. The method is combined with the finite element method to provide the entire time domain response. The results are compared with analytical solutions and experimental data, and excellent agreement is obtained. A fast Method of Moments is presented to calculate electromagnetic wave scattering from layered one dimensional rough surfaces. To facilitate the solution, the Forward Backward method with Spectral Acceleration is applied. As an example, a dielectric layer on a perfect electric conductor surface is studied. First, the numerical results are compared with the analytical solution for layered flat surfaces to partly validate the formulation. Second, the accuracy, efficiency, and convergence of the method are studied for various rough surfaces and layer permittivities. The Finite Difference Time Domain (FDTD) method is used to study metamaterials exhibiting both negative permittivity and permeability in certain frequency bands.(cont.) The structure under study is the well-known periodic arrangement of rods and split-ring resonators, previously used in experimental setups. For the first time, the numerical results of this work show that fields propagating inside the metamaterial with a forward power direction exhibit a backward phase velocity and negative index of refraction. A new metamaterial design is presented that is less lossy than previous designs. The effects of numerical dispersion in the FDTD method are investigated for layered, anisotropic media. The numerical dispersion relation is derived for diagonally anisotropic media. The analysis is applied to minimize the numerical dispersion error of Huygens' plane wave sources in layered, uniaxial media. For usual discretization sizes, a typical reduction of the scattered field error on the order of 30 dB is demonstrated. The new FDTD method is then used to study the Angular Correlation Function (ACF) of the scattered fields from continuous random media with and without a target object present. The ACF is shown to be as much as 10 dB greater when a target object is present for situations where the target is undetectable by examination of the radar cross section only.by Christopher D. Moss.Ph.D

Annual Review of Progress in Applied Computational Electromagnetics

Approved for public release; distribution is unlimited

Compact support wavelet representations for solution of quantum and electromagnetic equations: Eigenvalues and dynamics

Wavelet-based algorithms are developed for solution of quantum and electromagnetic differential equations. Wavelets offer orthonormal localized bases with built-in multiscale properties for the representation of functions, differential operators, and multiplicative operators. The work described here is part of a series of tools for use in the ultimate goal of general, efficient, accurate and automated wavelet-based algorithms for solution of differential equations. The most recent work, and the focus here, is the elimination of operator matrices in wavelet bases. For molecular quantum eigenvalue and dynamics calculations in multiple dimensions, it is the coupled potential energy matrices that generally dominate storage requirements. A Coefficient Product Approximation (CPA) for the potential operator and wave function wavelet expansions dispenses with the matrix, reducing storage and coding complexity. New developments are required, however. It is determined that the CPA is most accurate for specific choices of wavelet families, and these are given here. They have relatively low approximation order (number of vanishing wavelet function moments), which would ordinarily be thought to compromise both wavelet reconstruction and differentiation accuracy. Higher-order convolutional coefficient filters are determined that overcome both apparent problems. The result is a practical wavelet method where the effect of applying the Hamiltonian matrix to a coefficient vector can be calculated accurately without constructing the matrix. The long-familiar Lanczos propagation algorithm, wherein one constructs and diagonalizes a symmetric tridiagonal matrix, uses both eigenvalues and eigenvectors. We show here that time-reversal-invariance for Hermitian Hamiltonians allows a new algorithm that avoids the usual need to keep a number Lanczos vectors around. The resulting Conjugate Symmetric Lanczos (CSL) method, which will apply for wavelets or other choices of basis or grid discretization, is simultaneously low-operation-count and low-storage. A modified CSL algorithm is used for solution of Maxwell's time-domain equations in Hamiltonian form for non-lossy media. The matrix-free algorithm is expected to complement previous work and to decrease both storage and computational overhead. It is expected- that near-field electromagnetic solutions around nanoparticles will benefit from these wavelet-based tools. Such systems are of importance in plasmon-enhanced spectroscopies

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described