WOMBAT - Astrophysical Fluid Code
WOMBAT (sWift Objects for Mhd BAsed on Tvd) is a new implementation of the non-relativistic MHD TVD scheme described in Ryu & Jones, 1995[1], Ryu, Jones, & Frank, 1995[2], and Ryu, Miniati, Jones, & Frank, 1998[3]. An extension for relativistic MHD has recently been added. The code operates on 1, 2, and 3D Eulerian meshes (cartesian and cylindrical coordinates) with magnetic field divergence restriction controlled by a constrained transport (CT) scheme. A lot of effort has gone into creating hooks for the user to tune code performance to a given processor based on chip cache sizes. Proper settings yield significant speed-ups due to efficient cache reuse. The primary developer is Peter Mendygral with an enormous amount of wonderful parallelization technique ideas provided by David Porter and contributions/testing by Paul Edmon and Jennifer Delgado.
Scaling
WOMBAT has a hybrid MPI + OpenMP framework, and it exhibits near perfect weak and strong scaling up to at least 8192 cores by relying heavily on overlapped computation and communication.
 |
 |
WOMBAT is organized into psuedo objects that operate on a set of physics. For example, one object produces the MHD eigenvectors while another utilizes the results to produce fluxes at cell boundaries. This organization makes it quite straight forward to add additional physics objects. WOMBAT currently has the following additional physics:
- 2nd order, fully conservative gravity (static or evolving)
- Test particle cosmic rays (CRs) subjected to first-order Fermi acceleration at shocks, radiative losses, and adiabatic changes[4]
- CRs subject to full diffusive shock acceleration (DSA) with feedback on the ambient fluid[5] (see Multidimensional DSA)
- Directionally un-split hydrodynamic shock detection algorithm (e.g.[6])
- Radiative cooling based on the exact solution to a piecewise power-law cooling curve[7]
- Relativistic MHD[8-10]
Visualization
Output from WOMBAT comes in three flavors; restarts, compressed dumps, and visualization bricks-of-bytes (BOBs). The formats are all essentially raw files that are stacked into bricks corresponding to each MPI rank. Each rank is a member of an output team that sends output to an asynchronous I/O (AIO) server. The AIO servers were written by David Porter. This system allows all output to be written to disk while computation continues. The BOBs can be read into the LCSE visualization pipeline for rendering on the Powerwall. Interactive computing with WOMBAT is therefore a native operation.
Obtaining WOMBAT
Please visit the code request page for information on obtaining WOMBAT.
[1]Ryu, D., & Jones, T. W., 1995, ApJ, 442, 228 (ADS link)
[2]Ryu, D., Jones, T. W., & Frank, A., 1995, ApJ, 452, 785 (ADS link)
[3]Ryu, D., Miniati, F., Jones, T. W., & Frank, A., 1998, ApJ, 509, 244 (ADS link)
[4]Jones, T. W., & Kang, H., 2005, AstroPart, 24, 75-91 (ADS link)
[5]Edmon, P., & Jones, T. W., 2010, in progress
[6]Skillman, S. W., O'Shea, B. W., Hallman, E. J., Burns, J. O., & Norman, M. L., 2008, ApJ, 689, 1063-1077 (ADS Link)
[7]Townsend, R.H.D., 2009, ApJ, 181, 391-397 (ADS Link)
[8]Balsara, D., 2001, ApJS, 132, 83-101 (ADS Link)
[9]Anton, L., Miralles, J. A., Marti, J. M., Ibanez, J. M., Aloy, M. A., & Mimica, P. 2010, ApJS, 188, 1-31 (ADS Link)
[10]Beckwith, K. & Stone, J. M. 2011, ApJS, 193, 6 (ADS Link)