Systematic construction of upwind constrained transport schemes for MHD

JOURNAL OF COMPUTATIONAL PHYSICS

卷:424

Article

Mignone, A.¹  Del Zanna, L.^2,3,4 

[1] Univ Torino, Dipartimento Fis, Via Pietro Giuria1, I-10125 Turin, Italy

[2] Univ Firenze, Dipartimento Fis & Astron, Via G Sansone1, I-50019 Sesto Fiorentino, Italy

[3] INFN, Sez Firenze, Via G Sansone 1, I-50019 Sesto Fiorentino, Italy

[4] INAF, Osservatorio Astrofis Arcetri, Largo E Fermi 5, I-50125 Florence, Italy

关键词: Magnetohydrodynamics (MHD);    Methods: numerical;    Constrained transport;    Finite volume;    Riemann solvers;   

DOI  :  10.1016/j.jcp.2020.109748

来源: Elsevier

【 摘 要 】

The constrained transport (CT) method reflects the state of the art numerical technique for preserving the divergence-free condition of magnetic field to machine accuracy in multidimensional MHD simulations performed with Godunov-type, or upwind, conservative codes. The evolution of the different magnetic field components, located at zone interfaces using a staggered representation, is achieved by calculating the electric field components at cell edges, in a way that has to be consistent with the Riemann solver used for the update of cell-centered fluid quantities at interfaces. Albeit several approaches have been undertaken, the purpose of this work is, on the one hand, to compare existing methods in terms of robustness and accuracy and, on the other, to extend the upwind constrained transport (UCT) method by Londrillo & Del Zanna (2004) [22] and Del Zanna et al. (2007) [23] for the systematic construction of new averaging schemes. In particular, we propose a general formula for the upwind fluxes of the induction equation which simply involves the information available from the base Riemann solver employed for the fluid part, provided it does not require full spectral decomposition, and 1D reconstructions of velocity and magnetic field components from nearby intercell faces to cell edges. Our results are presented here in the context of second-order schemes for classical MHD, but they can be easily generalized to higher than second order schemes, either based on finite volumes or finite differences, and to other physical systems retaining the same structure of the equations, such as that of relativistic or general relativistic MHD. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free