skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: On the performance of exponential integrators for problems in magnetohydrodynamics

; ;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 330; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-07 03:24:44; Journal ID: ISSN 0021-9991
Country of Publication:
United States

Citation Formats

Einkemmer, Lukas, Tokman, Mayya, and Loffeld, John. On the performance of exponential integrators for problems in magnetohydrodynamics. United States: N. p., 2017. Web. doi:10.1016/
Einkemmer, Lukas, Tokman, Mayya, & Loffeld, John. On the performance of exponential integrators for problems in magnetohydrodynamics. United States. doi:10.1016/
Einkemmer, Lukas, Tokman, Mayya, and Loffeld, John. Wed . "On the performance of exponential integrators for problems in magnetohydrodynamics". United States. doi:10.1016/
title = {On the performance of exponential integrators for problems in magnetohydrodynamics},
author = {Einkemmer, Lukas and Tokman, Mayya and Loffeld, John},
abstractNote = {},
doi = {10.1016/},
journal = {Journal of Computational Physics},
number = C,
volume = 330,
place = {United States},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved functionals. We propose a new discretisation strategy for these equations based on a discrete variational principle applied to a formal Lagrangian. The resulting integrator preserves important quantities like the total energy, magnetic helicity and cross helicity exactly (up to machine precision). As the integrator is free of numerical resistivity, spurious reconnection along current sheets is absent in the ideal case. If effects of electron inertia are added, reconnection of magnetic field linesmore » is allowed, although the resulting model still possesses a noncanonical Hamiltonian structure. After reviewing the conservation laws of the model equations, the adopted variational principle with the related conservation laws is described both at the continuous and discrete level. We verify the favourable properties of the variational integrator in particular with respect to the preservation of the invariants of the models under consideration and compare with results from the literature and those of a pseudo-spectral code.« less
  • The effect of non-exponential and multi-exponential decay or relaxation behavior on the performance of the direct exponential curve resolution algorithm (DECRA) is investigated through a series of numerical simulations. Three different combinations of decay or relaxation behavior were investigated through DECRA analysis of simulated pulse gradient spin echo (PGSE) NMR diffusion spectra that contained the combination of two individual components. The diffusion decay behavior of one component was described by a single-exponential decay, while the second component was described by either (1) a multi-exponential decay, (2) a decay behavior described by the empirical Kohlrausch-Williams-Watts (KWW) relation or (3) a multi-exponentialmore » decay behavior correlated with variations in the NMR spectral line shape. The magnitudes and types of errors produced during the DECRA analysis of spectral data with deviations from a pure single-exponential decay behavior are presented. It is demonstrated that the deviation from single-exponential decay impacts the resulting calculated line shapes, the calculated relative concentrations and the quantitative estimation of the decay or relaxation time constants of both components present in the NMR spectra.« less
  • We introduce new projective versions of second-order accurate Runge-Kutta and Adams-Bashforth methods, and demonstrate their use as outer integrators in solving stiff differential systems. An important outcome is that the new outer integrators, when combined with an inner telescopic projective integrator, can result in fully explicit methods with adaptive outer step size selection and solution accuracy comparable to those obtained by implicit integrators. If the stiff differential equations are not directly available, our formulations and stability analysis are general enough to allow the combined outer-inner projective integrators to be applied to black-box legacy codes or perform a coarse-grained time integrationmore » of microscopic systems to evolve macroscopic behavior, for example.« less
  • Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative interactions. In this Letter, we develop the “slimplectic” integrator, a new type of numerical integrator that shares many of the benefits of traditional symplectic integrators yet is applicable to general nonconservative systems. We utilize a fixed-time-step variational integrator formalism applied to the principle of stationary nonconservative action developed in Galley et al. As a result, the generalized momenta and energy (Noether current) evolutions are well-tracked. Wemore » discuss several example systems, including damped harmonic oscillators, Poynting–Robertson drag, and gravitational radiation reaction, by utilizing our new publicly available code to demonstrate the slimplectic integrator algorithm. Slimplectic integrators are well-suited for integrations of systems where nonconservative effects play an important role in the long-term dynamical evolution. As such they are particularly appropriate for cosmological or celestial N-body dynamics problems where nonconservative interactions, e.g., gas interactions or dissipative tides, can play an important role.« less