On the performance of exponential integrators for problems in magnetohydrodynamics
- Univ. Innsbruck (Austria)
- Univ. of California, Merced, CA (United States)
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Exponential integrators have been introduced as an efficient alternative to explicit and implicit methods for integrating large stiff systems of differential equations. Over the past decades these methods have been studied theoretically and their performance was evaluated using a range of test problems. While the results of these investigations showed that exponential integrators can provide significant computational savings, the research on validating this hypothesis for large scale systems and understanding what classes of problems can particularly benefit from the use of the new techniques is in its initial stages. Resistive magnetohydrodynamic (MHD) modeling is widely used in studying large scale behavior of laboratory and astrophysical plasmas. In many problems numerical solution of MHD equations is a challenging task due to the temporal stiffness of this system in the parameter regimes of interest. In this paper we evaluate the performance of exponential integrators on large MHD problems and compare them to a state-of-the-art implicit time integrator. Both the variable and constant time step exponential methods of EPIRK-type are used to simulate magnetic reconnection and the Kevin–Helmholtz instability in plasma. Performance of these methods, which are part of the EPIC software package, is compared to the variable time step variable order BDF scheme included in the CVODE (part of SUNDIALS) library. We study performance of the methods on parallel architectures and with respect to magnitudes of important parameters such as Reynolds, Lundquist, and Prandtl numbers. We find that the exponential integrators provide superior or equal performance in most circumstances and conclude that further development of exponential methods for MHD problems is warranted and can lead to significant computational advantages for large scale stiff systems of differential equations such as MHD.
- Research Organization:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA); National Science Foundation (NSF)
- Grant/Contract Number:
- AC52-07NA27344; 1115978; P25346
- OSTI ID:
- 1860902
- Alternate ID(s):
- OSTI ID: 1398578
- Report Number(s):
- LLNL-JRNL-832030; 1049472; TRN: US2305993
- Journal Information:
- Journal of Computational Physics, Vol. 330; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Efficient adaptive step size control for exponential integrators
|
journal | October 2022 |
| Exponential collocation methods for conservative or dissipative systems | preprint | January 2017 |
| Exponential methods for solving hyperbolic problems with application to kinetic equations | text | January 2019 |
| Robust Quantum Optimal Control with Trajectory Optimization | preprint | January 2021 |
Similar Records
Partitioned exponential methods for coupled multiphysics systems
User documentation for PVODE, an ODE solver for parallel computers