{"title":"A kinetic energy and entropy preserving scheme for compressible magnetohydrodynamics","authors":"Takashi Shiroto, Tomo-Hiko Watanabe","doi":"10.1016/j.jcp.2025.114212","DOIUrl":null,"url":null,"abstract":"<div><div>This paper extends the kinetic energy and entropy preserving (KEEP) scheme, which was originally designed for compressible fluid dynamics, to include compressible magnetohydrodynamics (MHD). In MHD, the magnetic energy included in the total energy equation is consistent with that derived from the induction equation. Discretizing the spatial derivative in the induction equation using the central finite difference preserves the solenoidal constraint of the magnetic field. Numerical fluxes for the momentum and the total energy related to the magnetic field can be derived for the momentum and total energy from the Lorentz force and the solenoidal constraint, even in discrete form. Numerical results confirm that the KEEP scheme enables non-dissipative and stable MHD simulations, unlike the conventional central difference method.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"539 ","pages":"Article 114212"},"PeriodicalIF":3.8000,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125004954","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
This paper extends the kinetic energy and entropy preserving (KEEP) scheme, which was originally designed for compressible fluid dynamics, to include compressible magnetohydrodynamics (MHD). In MHD, the magnetic energy included in the total energy equation is consistent with that derived from the induction equation. Discretizing the spatial derivative in the induction equation using the central finite difference preserves the solenoidal constraint of the magnetic field. Numerical fluxes for the momentum and the total energy related to the magnetic field can be derived for the momentum and total energy from the Lorentz force and the solenoidal constraint, even in discrete form. Numerical results confirm that the KEEP scheme enables non-dissipative and stable MHD simulations, unlike the conventional central difference method.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.