Research reports

Years: 2024 2023 2022 2021 2020 2019 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2004 2003 2002 2001 2000 1999 1998 1997 1996 1995 1994 1993 1992 1991

Splitting-Based Structure Preserving Discretizations for Magnetohydrodynamics

by R. Hiptmair and C. Pagliantini

(Report number 2017-18)

Abstract
We start from the splitting of the equations of single-fluid magnetohydrodynamics (MHD) into a magnetic induction part and a fluid part. We design novel numerical methods for the MHD system based on the coupling of Galerkin schemes for the electromagnetic fields via finite element exterior calculus with finite volume methods for the conservation laws of fluid mechanics. Using a vector potential based formulation, the magnetic induction problem is viewed as an instance of a generalized transient advection problem of differential forms. For the latter, we rely on an Eulerian method of lines with explicit Runge-Kutta timestepping and on structure preserving spatial upwind discretizations of the Lie derivative based on the duality between the contraction of differential forms and the extrusion of chains. The balance laws for the fluid constitute a system of conservation laws with the magnetic induction field as a space and time dependent coefficient, supplied at every time step by the structure preserving discretization of the magnetic induction problem. We describe finite volume schemes based on approximate Riemann solvers adapted to accommodate the electromagnetic contributions to the momentum and energy conservation. A set of benchmark tests for the two- dimensional planar ideal MHD equations provide numerical evidence that the resulting lowest order coupled scheme has excellent conservation properties, is first order accurate for smooth solutions, conservative and stable.

Keywords: Magnetohydrodynamics, hyperbolic conservation laws, Lie derivative, structure preserving discretizations, discrete differential forms, finite volume methods.

@Techreport{HP17_714,
  author = {R. Hiptmair and C. Pagliantini},
  title = {Splitting-Based Structure Preserving Discretizations for
Magnetohydrodynamics},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2017-18},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2017/2017-18.pdf },
  year = {2017}
}

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).