skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Final Report for Time Domain Boundary Element and Hybrid Finite Element Simulation for Maxwell's Equations

Abstract

This report summarizes the work performed for Lawrence Livermore National Laboratory (LLNL) at the University of Washington between September 2004 and May 2006. This project studied fast solvers and stability for time domain integral equations (TDIE), especially as applied to radiating boundary for a massively parallel FEM solver.

Authors:
;
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
902353
Report Number(s):
UCRL-SR-228619
TRN: US200717%%552
DOE Contract Number:
W-7405-ENG-48
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; INTEGRAL EQUATIONS; LAWRENCE LIVERMORE NATIONAL LABORATORY; SIMULATION; STABILITY

Citation Formats

Pingenot, J, and Jandhyala, V. Final Report for Time Domain Boundary Element and Hybrid Finite Element Simulation for Maxwell's Equations. United States: N. p., 2007. Web. doi:10.2172/902353.
Pingenot, J, & Jandhyala, V. Final Report for Time Domain Boundary Element and Hybrid Finite Element Simulation for Maxwell's Equations. United States. doi:10.2172/902353.
Pingenot, J, and Jandhyala, V. Thu . "Final Report for Time Domain Boundary Element and Hybrid Finite Element Simulation for Maxwell's Equations". United States. doi:10.2172/902353. https://www.osti.gov/servlets/purl/902353.
@article{osti_902353,
title = {Final Report for Time Domain Boundary Element and Hybrid Finite Element Simulation for Maxwell's Equations},
author = {Pingenot, J and Jandhyala, V},
abstractNote = {This report summarizes the work performed for Lawrence Livermore National Laboratory (LLNL) at the University of Washington between September 2004 and May 2006. This project studied fast solvers and stability for time domain integral equations (TDIE), especially as applied to radiating boundary for a massively parallel FEM solver.},
doi = {10.2172/902353},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Mar 01 00:00:00 EST 2007},
month = {Thu Mar 01 00:00:00 EST 2007}
}

Technical Report:

Save / Share:
  • Hybrid finite element (FE)--boundary integral (BI) analysis of infinite periodic arrays is extended to include planar multilayered Green's functions. In this manner, a portion of the volumetric dielectric region can be modeled via the finite element method whereas uniform multilayered regions can be modeled using a multilayered Green's function. As such, thick uniform substrates can be modeled without loss of efficiency and accuracy. The multilayered Green's function is analytically computed in the spectral domain and the resulting BI matrix-vector products are evaluated via the fast spectral domain algorithm (FSDA). As a result, the computational cost of the matrix-vector products ismore » kept at O(N). Furthermore, the number of Floquet modes in the expansion are kept very few by placing the BI surfaces within the computational unit cell. Examples of frequency selective surface (FSS) arrays are analyzed with this method to demonstrate the accuracy and capability of the approach. One example involves complicated multilayered substrates above and below an inhomogeneous filter element and the other is an optical ring-slot array on a substrate several hundred wavelengths in thickness. Comparisons with measurements are included.« less
  • In this paper, the authors present a two-dimensional Galerkin finite element formulation of the Maxwell's equations in the time domain. The Galerkin element integrals are computed analytically and an explicit forward-backward time integration scheme is employed for advancing the resulting set of ordinary differential equations in time. The authors will compare two formulations of the simple 4-node bilinear element: (i) an equal-order interpolation element in which both the electric field and the magnetic field are approximated by bilinear basis functions; and (ii) a mixed-interpolation element in which the electric field is approximated as piecewise constant and the magnetic field asmore » piecewise bilinear functions. The mixed formulation may be viewed as a finite element analog to certain staggered finite difference representations. Numerical examples will be presented to evaluate the accuracy of these two elements.« less
  • In this paper, we present a two-dimensional Galerkin finite element formulation of the Maxwell's equations in the time domain. The Galerkin element integrals are computed analytically and an explicit forward-backward time integration scheme is employed for advancing the resulting set of ordinary differential equations in time. We will compare two formulations of the simple 4-node bilinear element: an equal-order interpolation element in which both the electric field and the magnetic field are approximated by bilinear basis functions; and a mixed-interpolation element in which the electric field is approximated as piecewise constant and the magnetic field as piecewise bilinear functions. Themore » mixed formulation may be viewed as a finite element analog to certain staggered finite difference representations. Numerical examples will be presented to evaluate the accuracy of these two elements. 5 refs., 3 figs.« less
  • MAX3D is a three-dimensional computer program for solving Maxwell's electromagnetic field equations. The equations solved by this code neglect the displacement current and charge density. This report documents Maxwell's electromagnetic field equations, the simplifying assumptions, and the boundary conditions. The field equations are then written in terms of a scalar potential and a vector potential. A variational principle, which is presented in terms of these potentials, is equivalent to Maxwell's theory and to the boundary conditions. Finite-element approximations applied to the variational principle give a set of linear algebraic equations that can be solved to approximate the electric and magneticmore » fields.« less
  • A new class of hybrid finite element methods for the numerical analysis of second-order elliptic boundary-value problems is presented. The methods are characterized by the use of particular solutions of the differential equation being solved, in contrast to conventional hybrid methods, in which polynomial approximations are used. As a model problem the Dirichlet problem for Laplace's equation is studied. A priori error estimates are derived, and the results of numerical experiments are presented. 3 figures, 1 table.