Time-Domain Solvers for Complex-Media Electrodynamics and Plasma Physics
In this dissertation, several extensions for two popular electromagnetic simulation methods: finite-different time-domain (FDTD) and finite-element time-domain (FETD), are presented. These extensions aim to increase the geometrical flexibility and modeling capabilities in simulation of Maxwell's equations. Since straight-forward extensions to these methods produce numerical artifacts that pollute the results and reduce accuracy, alternative strategies have been sought. Various methodologies are explored here to address these issues: A E-B mixed-vector FETD implementation based on first order Maxwell equations is introduced. In this method a mix of electric and magnetic field variables are used, where an edge element expansion is used for the electric field and face element expansion is used for the magnetic field. Compared to the standard FETD methods, it can produce several advantages without any significant computational drawback: (i) it eliminates the spurious linear growth in time that may exist in the standard schemes; (ii) it produces energy-conserving schemes under appropriate time-discretization; (iii) it can be easily extended to frequency-dispersive media; (iv) it provides a natural path for hybridization with FDTD. Exploiting item (iii), E-B mixed-vector FETD is extended to inhomogeneous doubly-dispersive media. A conformal-PML implementation is also proposed for efficient simulation of open-domain boundaries by significantly reducing the buffer regions in the computational domain. Despite the advantages, finite-element simulations require solution of a linear system of equations which, in most practical problems, highly computationally intense. A hybrid FDTD-FETD method is introduced to alleviate this issue. The hybridization can result in significant optimization in computational cost by assigning detailed portion of the computational domain to FETD, while assigning the remaining to FDTD. Finally, a subgridding by domain-overriding (SGDO) methodology for full electromagnetic particle-in-cell (PIC) simulations is presented. Combined with relaxation methods, SGDO can produce significant improvements in PIC simulation accuracy.
School:The Ohio State University
School Location:USA - Ohio
Source Type:Master's Thesis
Keywords:fdtd fetd finite difference element complex media plasma physics pic particle in cell
Date of Publication:01/01/2008