Research reports

Entropy stable numerical schemes for two-fluid MHD equations

by H. Kumar and S. Mishra

(Report number 2011-22)

Abstract
Two-fluid ideal magnetohydrodynamics (MHD) equations are a generalized form of the ideal MHD equations in which electrons and ions are considered as separate species. The design of efficient numerical schemes for the these equations is complicated on account of their non-linear nature and the presence of stiff source terms, especially for high charge to mass ratios. In this article, we design entropy stable finite difference schemes for the two-fluid equations by combining entropy conservative fluxes and suitable numerical diffusion operators. Furthermore, to overcome the time step restrictions imposed by the stiff source terms, we devise time-stepping routines based on implicit-explicit (IMEX)-Runge Kutta (RK) schemes. The special structure of the two-fluid MHD equations is exploited by us to design IMEX schemes in which only local (in each cell) linear equations need to be solved at each time step. Benchmark numerical experiments are presented to illustrate the robustness and accuracy of these schemes.

Keywords: Convection-dffusion problem, discrete differential forms, discrete Lie derivative, semi-Lagrangian methods

BibTeX
@Techreport{KM11_98,
  author = {H. Kumar and S. Mishra},
  title = {Entropy stable numerical schemes for two-fluid MHD equations},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2011-22},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2011/2011-22.pdf },
  year = {2011}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser