skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Constrained-transport magnetohydrodynamics with adaptive mesh refinement in CHARM

Abstract

We present the implementation of a three-dimensional, second-order accurate Godunov-type algorithm for magnetohydrodynamics (MHD) in the adaptive-mesh-refinement (AMR) cosmological code CHARM. The algorithm is based on the full 12-solve spatially unsplit corner-transport-upwind (CTU) scheme. The fluid quantities are cell-centered and are updated using the piecewise-parabolic method (PPM), while the magnetic field variables are face-centered and are evolved through application of the Stokes theorem on cell edges via a constrained-transport (CT) method. The so-called multidimensional MHD source terms required in the predictor step for high-order accuracy are applied in a simplified form which reduces their complexity in three dimensions without loss of accuracy or robustness. The algorithm is implemented on an AMR framework which requires specific synchronization steps across refinement levels. These include face-centered restriction and prolongation operations and a reflux-curl operation, which maintains a solenoidal magnetic field across refinement boundaries. The code is tested against a large suite of test problems, including convergence tests in smooth flows, shock-tube tests, classical two- and three-dimensional MHD tests, a three-dimensional shock-cloud interaction problem, and the formation of a cluster of galaxies in a fully cosmological context. The magnetic field divergence is shown to remain negligible throughout.

Authors:
 [1];  [2]
  1. ETH, Zurich (Switzerland)
  2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21), Scientific Discovery through Advanced Computing (SciDAC) (SC21.1 )
OSTI Identifier:
1511331
Grant/Contract Number:  
AC02-05CH11231
Resource Type:
Accepted Manuscript
Journal Name:
Astrophysical Journal, Supplement Series
Additional Journal Information:
Journal Volume: 195; Journal Issue: 1; Journal ID: ISSN 0067-0049
Publisher:
American Astronomical Society
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS

Citation Formats

Miniati, Francesco, and Martin, Daniel F. Constrained-transport magnetohydrodynamics with adaptive mesh refinement in CHARM. United States: N. p., 2011. Web. doi:10.1088/0067-0049/195/1/5.
Miniati, Francesco, & Martin, Daniel F. Constrained-transport magnetohydrodynamics with adaptive mesh refinement in CHARM. United States. doi:10.1088/0067-0049/195/1/5.
Miniati, Francesco, and Martin, Daniel F. Thu . "Constrained-transport magnetohydrodynamics with adaptive mesh refinement in CHARM". United States. doi:10.1088/0067-0049/195/1/5. https://www.osti.gov/servlets/purl/1511331.
@article{osti_1511331,
title = {Constrained-transport magnetohydrodynamics with adaptive mesh refinement in CHARM},
author = {Miniati, Francesco and Martin, Daniel F.},
abstractNote = {We present the implementation of a three-dimensional, second-order accurate Godunov-type algorithm for magnetohydrodynamics (MHD) in the adaptive-mesh-refinement (AMR) cosmological code CHARM. The algorithm is based on the full 12-solve spatially unsplit corner-transport-upwind (CTU) scheme. The fluid quantities are cell-centered and are updated using the piecewise-parabolic method (PPM), while the magnetic field variables are face-centered and are evolved through application of the Stokes theorem on cell edges via a constrained-transport (CT) method. The so-called multidimensional MHD source terms required in the predictor step for high-order accuracy are applied in a simplified form which reduces their complexity in three dimensions without loss of accuracy or robustness. The algorithm is implemented on an AMR framework which requires specific synchronization steps across refinement levels. These include face-centered restriction and prolongation operations and a reflux-curl operation, which maintains a solenoidal magnetic field across refinement boundaries. The code is tested against a large suite of test problems, including convergence tests in smooth flows, shock-tube tests, classical two- and three-dimensional MHD tests, a three-dimensional shock-cloud interaction problem, and the formation of a cluster of galaxies in a fully cosmological context. The magnetic field divergence is shown to remain negligible throughout.},
doi = {10.1088/0067-0049/195/1/5},
journal = {Astrophysical Journal, Supplement Series},
number = 1,
volume = 195,
place = {United States},
year = {2011},
month = {6}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 24 works
Citation information provided by
Web of Science

Save / Share: