Compressible Astrophysics: CASTRO
CASTRO Contributors We would like to gratefully acknowledge the direct contributions to the CASTRO code by:
You can find the latest version of the CASTRO repository at https://github.com/BoxLib-Codes/Castro
The CASTRO Users' Guide is available here.
CASTRO Email Support List We have set up a mailing list for CASTRO users to communicate with one another and with the developers. We will be using this list to announce new features in CASTRO, as well as any changes in CASTRO or BoxLib which may affect people running the CASTRO code.
To use the list, you subscribe here:
You can e-mail everyone on the list by sending to
Any e-mails from the list will have a subject starting with "[Castro-help]" so you can easily filter those.
To see past conversations, the archives are here:
Please encourage anyone who is interested in CASTRO to use this list, as it will make supporting the code easier for everyone.
A. S. Almgren, V.E. Beckner, J.B. Bell, M.S. Day, L.H. Howell, C.C. Joggerst, M.J. Lijewski, A. Nonaka, M. Singer, M. Zingale, "CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity", Astrophysical Journal, 715, 1221-1238, June 2010. [pdf] [IOP] [arxiv]
W. Zhang, L. Howell, A. Almgren, A. Burrows, and J. Bell "CASTRO: A New Compressible Astrophysical Solver. II. Gray Radiation Hydrodynamics", Astrophysical Journal Supplement Series, 196, 20, (2011) [pdf]
Weiqun Zhang, L. Howell, A. Almgren, A. Burrows, J. Dolence, J. Bell, "CASTRO: A New Compressible Astrophysical Solver. III. Multigroup Radiation Hydrodynamics", Astrophysical Journal Supplement Series, 204, 7, 2013. [pdf]
Software FrameworkThe CASTRO software is based on the hybrid C++/Fortran BoxLib software framework developed by CCSE.
Additional Publications Using CASTRO
J.C. Dolence, A. Burrows, and W. Zhang, "Two-Dimensional Core-Collapse Supernova Models with Multi-Dimensional Transport," submitted to ApJ [arxiv].
C. M. Malone, A. Nonaka, S. E. Woosley, A. S. Almgren, J. B. Bell, S. Dong, and M. Zingale, "The Deflagration Stage of Chandrasekhar Mass Models for Type Ia Supernovae: I. Early Evolution", Astrophysical Journal, 782, 11, 2014. [iop]
Ke-Jung Chen, Alexander Heger, Stan Woosley, Ann Almgren, and Daniel Whalen, "Pair Instability Supernovae of Very Massive Population III Stars" Astrophysical Journal, 792, 44, 2014. [arxiv].
Ke-Jung Chen, Alexander Heger, S.E. Woosley, Ann S. Almgren, and Daniel J. Whalen, "Two-Dimensional Simulations of Pulsational Pair-Instability Supernova", Astrophysical Journal, 792, 28, 2014. [arxiv].
Ke-Jung Chen, Alexander Heger, S.E. Woosley, Ann Almgren, and Daniel J. Whalen, and Jarrett L. Johnson, "The General Relativitistic Instability Supernova of a Supermassive Population III Star", Astrophysical Journal, 790, 162, 2014.
Ke-Jung Chen, Alexander Heger, and Ann S. Almgren, "Numerical Approaches for Multidimensional Simulations of Stellar Explosions", Astronomy and Computing, 3-4, pp. 70-78, Nov.-Dec. 2013. [doi]
K. Chen, A. Heger, S. Woosley, A. Almgren, and W. Zhang, "The Most Powerful Stellar Explosion", Bulletin of the American Physical Society, 2013, vol. 58, no. 4.
Haitao Ma, Stan Woosley, Chris Malone, Ann Almgren, and J.B. Bell, "Carbon Deflagration in Type Ia Supernovae: I. Centrally Ignited Models", Astrophysical Journal, to appear.
Adam Burrows, Joshua C. Dolence, Jeremiah W. Murphy, "An Investigation into the Character of Pre-Explosion Core-Collapse Supernova Shock Motion," [arxiv].
Emmanouela Rantsiou, A. Burrows, J. Nordhaus, Ann Almgren, "Induced Rotation in Three-Dimensional Simulations of Core-Collapse Supernovae: Implications for Pulsar Spins," Astrophysical Journal, 732, 57, 2011. [arxiv].
A. Almgren, J. Bell, D. Kasen, M. Lijewski, A. Nonaka, P. Nugent, C. Rendleman, R. Thomas, M. Zingale, "MAESTRO, CASTRO and SEDONA -- Petascale Codes for Astrophysical Applications," SciDAC 2010, J. of Physics: Conference Series, Chattanooga, Tennessee, July 2010. [arxiv]
C. C. Joggerst, A. Almgren, and S. E. Woosley, "Three-dimensional simulations of Rayleigh-Taylor mixing in core collapse supernovae with CASTRO", Astrophysical Journal, 723, 353, October 2010.
C. C. Joggerst, A. Almgren, J. Bell, Alexander Heger, Daniel Whalen, and S. E. Woosley, "Primordial Core-Collapse Supernovae and the Chemical Abundances of Metal-Poor Stars," Astrophysical Journal, 709, 11-26, January 2010 [arxiv].
Ke-Jung Chen, Alexander Heger, and Ann Almgren, "Two Dimensional Simulations of Pair-Instability Supernovae", Conference on The First Stars and Galaxies: Challenges for the Next Decade, Austin, Texas, March 8-11, 2010, Conference Proceedings published by the American Institute of Physics. [AIP]
Posters Using CASTRO
A. Almgren, J. Bell, M. Day, L. Howell, C. Joggerst, E. Myra, A. Nonaka, J. Nordhaus, M. Singer, and M. Zingale, "CASTRO: A New AMR Radiation-Hydrodynamics Code for Compressible Astrophyics, presented at the American Astronomical Society (AAS) Meeting in Washington D.C., January 3--7, 2010. [pdf]
OverviewAs part of the SciDAC Computational Astrophysics Consortium, CCSE, in collaboration with Louis Howell and Mike Singer at LLNL, have developed CASTRO, a new multi-dimensional Eulerian AMR radiation-hydrodynamics code that includes stellar equations of state, nuclear reaction networks, and self-gravity. Initial target applications for CASTRO include Type Ia and Type II supernovae.
Coordinate SystemsCASTRO supports calculations in 1-d, 2-d and 3-d Cartesian coordinates, as well as 1-d spherical and 2-d cylindrical (r-z) coordinate systems.
HydrodynamicsTime integration of the hydrodynamics equations is based on an unsplit version of the the piecewise parabolic method (PPM) with new limiters that avoid reducing the accuracy of the scheme at smooth extrema.
Equation of State
CASTRO can follow an arbitrary number of isotopes or elements. The atomic weights and amounts of these elements are used to calculate the mean molecular weight of the gas required by the equation of state.
Nuclear Reaction Networks
CASTRO has a number of nuclear reaction networks available with the code. The reactions are included in the time integration scheme in a second-order accurate operator-split formulation (Strang splitting).
Radiation in CASTRO
CASTRO supports several different approaches to solving for self-gravity. The most general is a full Poisson solve for the gravitational potential. In this case, the user controls the choice of boundary conditions for the gravitational potentialat ''outflow'' boundaries, either homogeneous Dirichlet, or calculated from a monopole approximation at the coarsest level.
CASTRO also supports a monopole approximation for gravity. For the monopole approximation, a radial average of the density is computed from the full grid to create a one-dimensional density profile. This profile is then used to compute a one-dimensional gravitational potential which is then mapped back onto the full grid.
A constant gravity option is also available, where the user specifies the magnitude of the gravitational vector and the direction is assumed to align with the y-axis in 2-d, and the z-axis in 3-d.
AMR in CASTRO
Our approach to adaptive refinement in CASTRO uses a nested hierarchy of logically-rectangular grids with simultaneous refinement of the grids in both space and time. The integration algorithm on the grid hierarchy is a recursive procedure in which coarse grids are advanced in time, fine grids are advanced multiple steps to reach the same time as the coarse grids and the data at different levels are then synchronized. The AMR methodology was introduced by Berger and Oliger (1984) for hyperbolic problems; it has been demonstrated to be highly successful for gas dynamics by Berger and Colella (1989) in two dimensions and by Bell et al. (1994) in three dimensions. During the regridding step, increasingly finer grids are recursively embedded in coarse grids until the solution is sufficiently resolved. An error estimation procedure based on user-specified criteria evaluates where additional refinement is needed and grid generation procedures dynamically create or remove rectangular fine grid patches as resolution requirements change.
Software FrameworkThe CASTRO software is written in C++ and Fortran, and is based on the BoxLib software framework developed by CCSE. If you are interested in becoming a CASTRO user, please contact Ann Almgren.
Questions?Contact Ann Almgren.
The work at LBNL was supported by the Office of High Energy Physics and the Office of Mathematics, Information, and Computational Sciences as part of the SciDAC Program under the U.S. Department of Energy under contract No. DE-AC02-05CH11231.
The work performed at LLNL was under the auspices of the U.S. Department of Energy under contract No. DE-AC52-07NA27344.
Mike Zingale was supported by Lawrence Livermore National Lab under contracts B568673, B574691, and B582735.
We thank Haitao Ma, Jason Nordhaus and Ken Chen for being patient early users of CASTRO.