- Mathematical Modeling
- Numerical Analysis
- Scientific Computing
Computational simulation is an essential research tool in plasma and accelerator physics. Physical experiments in these fields are often prohibitively costly or complex, and sometimes even impractical or impossible. (For example, the extreme temperature of close proximity to the sun would destroy any space probe sent to collect data on solar plasma near its surface.) Hence researchers conduct virtual experiments via computer simulations as an alternative. Additionally, practice runs of experiments are often simulated before the true physical experiment actually takes place, thereby serving as a tool to plan, design, and optimize the experiment. Computational plasma physics has in fact been a driving force behind the modern field of scientific computing.
At a fundamental descriptive level plasma is simply a fluid consisting of disassociated electrons and ions. Thus a natural approach to describe plasma is to model it as a collection of charged particles. A key problem in these particle simulation models is solving for charged particle dynamics in which their motion is primarily driven by electromagnetic fields, a task called the particle pushing problem. In the case of strongly magnetized plasma this problem is computationally challenging due to multi-scale phenomenon. On the one hand, charged particles gyrate about magnetic field lines on fast time scales in a phenomenon known as gyromotion. On the other hand, however, computer simulations typically evolve over significantly longer time scales in order to determine long term macroscopic behavior of the system. This wide disparity in time scales between fast-scale microscopic gyromotion and slow-scale macroscopic behavior makes the particle pushing problem numerically stiff. This difficulty is further compounded by the fact that any reasonably accurate computational model of plasma requires simulation of a large number of particles. Consequently, the particle pushing problem is a computationally intensive task.
This research project is a collaborative effort with scientists at the Lawrence Livermore National Laboratory directed at finding more computationally efficient algorithms to solve the particle pushing problem. We are investigating novel approaches to particle pushing using a technique called exponential integration. Exponential integrators are a lesser-known class of numerical time integrators that solve linear problems exactly (to within computational limits) and are A-stable. With these favorable accuracy and stability properties exponential integrators offer a promising alternative to conventional numerical methods.
ACKNOWLEDGEMENTS: This work is supported by the National Science Foundation under Award No. 1840265 administrated through the UC Merced Department of Applied Mathematics Data-Intentsive Research And Computing (DIRAC) Research Training Group (RTG). In addition, part of this work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52 07NA27344.
DISCLAIMER: Any opinions, findings, and conclusions or recommendations expressed here are those of the author and do not necessarily reflect the views of the the UC Merced Department of Applied Mathematics, UC Merced, the National Science Foundation, Lawrence Livermore National Laboratory, or the U.S. Department of Energy.