A numerical algorithm for MHD of free surface flows at low magnetic Reynolds numbers

  • Authors:

  • Roman Samulyak;Jian Du;James Glimm;Zhiliang Xu

  • Affiliations:

  • Computational Science Center, Brookhaven National Laboratory, Upton, NY 11973, USA;Department of Applied Mathematics and Statistics, SUNY at Stony Brook, Stony Brook, NY 11794, USA;Computational Science Center, Brookhaven National Laboratory, Upton, NY 11973, USA and Department of Applied Mathematics and Statistics, SUNY at Stony Brook, Stony Brook, NY 11794, USA;Computational Science Center, Brookhaven National Laboratory, Upton, NY 11973, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.45

Visualization

Abstract

We have developed a numerical algorithm and computational software for the study of magnetohydrodynamics (MHD) of free surface flows at low magnetic Reynolds numbers. The governing system of equations is a coupled hyperbolic-elliptic system in moving and geometrically complex domains. The numerical algorithm employs the method of front tracking and the Riemann problem for material interfaces, second order Godunov-type hyperbolic solvers, and the embedded boundary method for the elliptic problem in complex domains. The numerical algorithm has been implemented as an MHD extension of FronTier, a hydrodynamic code with free interface support. The code is applicable for numerical simulations of free surface flows of conductive liquids or weakly ionized plasmas. The code has been validated through the comparison of numerical simulations of a liquid metal jet in a non-uniform magnetic field with experiments and theory. Simulations of the Muon Collider/Neutrino Factory target have also been discussed.