Topic/Type: 1. Plasma Simulation, Poster
A.J. Christlieb, B.W. Ong
Michigan State University, Department of Mathematics - USA
Basic plasma science plays an increasingly significant role in applications of importance to the United States Air Force. Many of these applications require a fully kinetic description in at least part of the domain. The most common approach is to use a fully Lagrangian framework, where the model is reduced to tracking the evolution of test particles in phase space. Of the many varieties, the most accepted approach is Particle-In-Cell. The PI, with his collaborators, is developing an alternative approach, the Boundary Integral Treecode (BIT). BIT is based on fast summation algorithms and boasts arbitrary accuracy.
In this talk, we demonstrate how the existing BIT framework can be extended to electrodynamic problems. The key idea is to consider the transpose of the standard method of lines methodology, i.e., we choose to discretize in time and directly solve the resulting Helmholtz equation using an integral formulation. Observe that Maxwell\'s equations can be cast as
where is the charge density, is the current density, (, ) are the permittivity and permeability and is the speed of light. To illustrate the implicit treecode time stepping methodology for and , it suffices to consider the wave equation, . Using a centered difference approximation to and evaluating at time level gives,
The integral solution for is
where , the free space Green\'s function for the Helmholtz operator , is in . Here, and is the normalization. The volumetric and boundary integrals are approximated at the midpoints, and the resulting sums can be computed using fast summation algorithms which are time in evaluating the sum. We have applied this approach to the wave equation in 1D, 2D and 3D with Direcltet boundary conditions. For the wave equation, the method is able to take time steps much larger than the imposed CFL restrictions for an explicit integrator. Additionally, this implicit time stepping methodology has been extended to higher order approximations in time. We working on extned these ideas to Maxwell's equations.