EULAG, a computational model for multiscale flows: An MHD extension

EULAG is an established high-performance computational model for simulating fluid flows across a wide range of scales and physical scenarios [Prusa et al., Comput. Fluids 37 (2008) 1193]. Historically driven by interests in simulating weather and climate processes, the numerics of EULAG are unique,...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2013-03, Vol.236, p.608-623
Main Authors: Smolarkiewicz, Piotr K., Charbonneau, Paul
Format: Article
Language:English
Subjects:
Anelastic equations
Computation
Computational fluid dynamics
Computer simulation
EULAG
Fluid flow
Fluids
Magnetohydrodynamics
Magnetohydrodynamics (MHDs)
Mathematical models
MHD
Multidimensional positive definite advection transport algorithm
Nonoscillatory forward-in-time schemes
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:EULAG is an established high-performance computational model for simulating fluid flows across a wide range of scales and physical scenarios [Prusa et al., Comput. Fluids 37 (2008) 1193]. Historically driven by interests in simulating weather and climate processes, the numerics of EULAG are unique, owing to a synergistic blend of non-oscillatory forward-in-time MPDATA methods, robust elliptic solver, and generalized coordinate formulation enabling grid adaptivity. In this paper the numerical apparatus of an ideal magnetohydrodynamic (MHD) extension of the EULAG model is discussed, the robust workings of which have been recently revealed in global large-eddy simulations of solar magneto-convection producing solar-like magnetic cycles and dynamo action [Ghizaru et al., ApJL 715 (2010) L133; Racine et al., ApJ 735 (2011) 46]. Here, a specialized nonoscillatory forward-in-time scheme for integrating ideal anelastic MHD equations is presented in detail, and illustrated with an abstract example of magnetized three-dimensional flow in time-dependent geometry for a weak, moderate and strong magnetic field. An analysis of the model performance reveals that multiple solutions of elliptic problems do not have to imply proportionally larger computational expense.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2012.11.008