Document Actions

ZEUS-MP 1.0

by streeter — last modified 2007-03-30 03:36

Summary

ZEUS-MP is a computational fluid dynamics code developed at the Laboratory for Computational Astrophysics (NCSA, SDSC, University of Illinois at Urbana-Champaign, UC San Diego) for the simulation of astrophysical phenomena. ZEUS-MP is developed based on ZEUS-3D and parallelized using the MPI message-passing library. ZEUS-MP solves problems in three spatial dimensions with a wide variety of boundary conditions. A source-code preprocessor allows the user to customize the ZEUS-MP algorithm for the desired physics, geometry, and output.

Physics Included

ZEUS-MP solves the equations of ideal (non-resistive), non-relativistic, hydrodynamics and magnetohydrodynamics, including externally applied gravitational fields and self-gravity. The gas can be adiabatic or isothermal, and the thermal pressure is isotropic. Boundary conditions may be specified as reflecting, periodic, inflow, or outflow.

Geometries

ZEUS-MP solves problems in Cartesian (x, y, z), cylindrical-polar (z, r, phi), and spherical-polar (r, theta, phi) coordinate systems.

Algorithm

ZEUS-MP is a finite difference code with a fixed or moving orthogonal Eulerian mesh.
Time integration is fully explicit, with the time step controlled by the Courant condition.
The mesh is staggered, with zone-centered scalar quantities (density, energy) and face-centered vector components (velocity, magnetic field) to increase accuracy.
The algorithm is written in a "covariant" form to minimize the effect of the choice of coordinate systemson the structure of the code.
von-Neumann-Richtmyer artificial viscosity is used to spread shocks over several mesh zones.
The upstream weighted, monotonicity-preserving advection schemes are exactly conservative.
The magnetic field components are evolved in a manner that ensures that the divergence vanishes and that all possible MHD modes propagate accurately and stably (Evans and Hawley Constrained Transport scheme modified by Stone and Norman using the Method of Characteristics).
A Fast Fourier Transform algorithm (FFTW) is included for solving Poisson's equation for periodic boudary conditions. For more general boundary conditions, a multigrid algorithm (MGMPI) can be used to solve Poisson's equation for the gravitational potential.

Data I/O

Users can output data in ASCII format or in Hierarchical Data Format (HDF) developed and supported by NCSA.

Test Suite

There are 5 simple test problems accompanying the ZEUS-MP code in this release. The test problems are for testing the ZEUS-MP code and serve as demonstrations of the wide variety of physical problems that ZEUS-MP can handle. The following is a list of the test problems:

Supersonic Jet (in Cartesian coordinates)
Sedov Blastwave (in spherical coordinates)
1D Shock Tube (in cylindrical coordinates)
MHD Blast Problem (in Cartesian, spherical, and cylindrical coordinates)
Gravitational Collapse of gas sphere (in Cartesian, spherical, and cylindrical coordinates)

User Guide

Information about compilation and running test suite of ZEUS-MP are described in a brief online user guide. The user guide also contains some information about the postprocessing of data, data output format and the Poisson solvers.

The Code

Here is the tarred file of ZEUS-MP source code and related utility files (version 1.0b).

Benchmarks

For More Information

NCSA/LCA Technical Report "OPTIMIZATION AND SCALING OF SHARED-MEMORY AND MESSAGE-PASSING IMPLEMENTATIONS OF THE ZEUS HYDRODYNAMICS ALGORITHM" -- Comparison of the ZEUS-MP and ZEUS-3D implementations and performance on the SGI POWER Challenge Array and HP/Convex Exemplar SPP-1200.
See the series of 3 papers on ZEUS-2D in the Astrophysical Journal Supplement Series, vol. 80, 1992:

LCA Portal

Sections

Personal tools

Document Actions

ZEUS-MP 1.0