Modeling Electromagnetic Radiation in the Time Domain using the Finite Element Method
Most of my before this involved the solution of Maxwell's Equations in the frequency domain (which is generally simpler). Modeling the open region, where waves can radiate outwards, never to return, is one of the “Holy Grails” of computational electromagnetics. In this work, I start with Maxwell's equations (in two dimensions now, for simplicity)
The first expression is the electric field curl equation, the second and third equations form the components of the magnetic field curl equation. The last expression is the divergence of the electric field (which vanishes because we are assuming a chargefree region). Note that these equations are equivalent to solving a corresponding second order scalar wave equation for h_{z} . The z direction is assumed to be out of the page (the xy plane is in the plane of the page). This corresponds to the TM propagation modes. TE modes can be modeled by replacing the components of H with E and E with H and the wave impedance Z_{0} with wave conductance G_{0}. The first curl equation will will become inhomgeneous, with source term J_{z} while the second and third equations' source terms vanish.
By discretising these equations in time (using a CrankNicolson or implicit Euler scheme) and constructing a highorder finiteelement discretisation in space, one can generate accurate representation of large bandwidth electromagnetic pulses in the time domain. By assuming that the waves obey
on the exit surface (where p_{0}, p_{1} are constants which govern the reflection from our wave from the surface and is the curvature of the boundary, the subscripts n and s indicate directions normal and tangential to the boundary surface), we can simulate, with varying degrees of success, the disappearence of wave energy into the exterior, freespace region.
Material available for download
File 
Description 
Last Modified on 

Paper that appears in the Proceedings of the 2004 WSEAS conference on Applied Mathematics, Corfu, Greece. 
07/10/04 

A draft (unpublished) paper which covers aspects of constructing the absorbing boundary condition using the leastsquares finite element method. 
25/09/2005 

Talk presented to the 2004 WSEAS conference on Applied Mathematics. 
06/06/04 
Last updated 20050925, Bill Slade