Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userWENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
67
0
0.0
(
0
)
Authors
J. M.F. Donnert
Ask ChatGPT about the research
No Arabic abstract
Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine Hazel Hen, with a sustained performance of about 5 percent of peak at scale.
rate research
Read More
In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on small scales via Faraday tomography. This leap forward demands new numerical models for the amplification of magnetic fields in cosmic structure formation - the cosmological magnetic dynamo. Here we present a novel numerical approach to astrophyiscal MHD simulations aimed to resolve this small-scale dynamo in future cosmological simulations. As a first step, we implement a fifth order WENO scheme in the new code WOMBAT. We show that this scheme doubles the effective resolution of the simulation and is thus less expensive than common second order schemes. WOMBAT uses a novel approach to parallelization and load balancing developed in collaboration with performance engineers at Cray Inc. This will allow us scale simulation to the exaflop regime and achieve kpc resolution in future cosmological simulations of galaxy clusters. Here we demonstrate the excellent scaling properties of the code and argue that resolved simulations of the cosmological small scale dynamo within the whole virial radius are possible in the next years.
We investigate the applicability of curvilinear grids in the context of astrophysical simulations and WENO schemes. With the non-smooth mapping functions from Calhoun et al. (2008), we can tackle many astrophysical problems which were out of scope with the standard grids in numerical astrophysics. We describe the difficulties occurring when implementing curvilinear coordinates into our WENO code, and how we overcome them. We illustrate the theoretical results with numerical data. The WENO finite difference scheme works only for high Mach number flows and smooth mapping functions whereas the finite volume scheme gives accurate results even for low Mach number flows and on non-smooth grids.
We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium.
We demonstrate that, for the case of quasi-equipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics alpha-model (LAMHD) reproduces well both the large-scale and small-scale properties of turbulent flows; in particular, it displays no increased (super-filter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the sub-filter-scales. This is in contrast to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes alpha-model is somewhat limited in its applications because of the formation of spatial regions with no internal degrees of freedom and subsequent contamination of super-filter-scale spectral properties. No such regions are found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of 200 when compared to a direct numerical simulation on a large grid of 1536^3 points at the same Reynolds number.
We describe the magnetohydrodynamics (MHD) code CRONOS, which has been used in astrophysics and space physics studies in recent years. CRONOS has been designed to be easily adaptable to the problem at hand, where the user can expand or exchange core modules or add new functionality to the code. This modularity comes about through its implementation using a C++ class structure. The core components of the code include solvers for both hydrodynamical (HD) and MHD problems. These problems are solved on different rectangular grids, which currently support Cartesian, spherical, and cylindrical coordinates. CRONOS uses a finite-volume description with different approximate Riemann solvers that can be chosen at runtime. Here, we describe the implementation of the code with a view toward its ongoing development. We illustrate the codes potential through several (M)HD test problems and some astrophysical applications.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
