Compressible Astrophysics: CASTRO
|This is an animation of nucleosynthesis in a Type Ia supernovae from a simulation run by Haitao Ma at UC Santa Cruz. This movie was created using the VisIt visualization software.|
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 at http://boxlib-codes.github.io/Castro.
CASTRO Email Support List
If you are interested in using CASTRO, please join our CASTRO mailing list to
receive any updates and see questions other users are asking:
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
Ke-Jung Chen, Alexander Heger, Stan Woosley, Ann Almgren and Daniel J. Whalen, "Pair-Instability Supernovae of Non-Zero Metallicity Stars," Numerical Modeling of Space Plasma Flows, ASTRONUM-2014, 498, 2015. [arxiv].
Ke-Jung Chen, "Pair-instability supernovae of fast rotating stars," Modern Physics Letters A, 30, Issue 2, 2015.
Max P. Katz, Michale Zingale, Alan C. Calder, F. Douglas Swesty, Ann S. Almgren, Weiqun Zhang, "White Dwarf Mergers on Adaptive Meshes. I. Methodology and Code Verification", submitted for publication, 2015.
Joshua C. Dolence, Adam Burrows, and Weiqun Zhang, "Two-Dimensional Core-Collapse Supernova Models with Multi-Dimensional Transport," Astrophysical Journal, 800, No. 1, 2015. [arxiv].
Ke-Jung Chen, "Supernovae at the Cosmic Dawn," International Jounrnal of Modern Physics D, 23, No. 5, 2014.
Joseph Smidt, Daniel J. Whalen, Brandon K. Wiggins, Wesley Even, Jarrett L. Johnson, and Chris L. Fryer, "Population III Hypernovae," Astrophysical Journal, 797, No. 2, 2014.
Daniel J. Whalen, Joseph Smidt, Alexander Heger, Raphael Hirschi, Norhasliza Yusof, Wesley Even, Chris L. Fryer, Massimo Stiavelli, Ke-Jung Chen, and Candace C. Joggerst, "Pair-Instability Supernovae in the Local Universe," Astrophysical Journal, 797, No. 1, 2014.
Candace C. Joggerst, Anthony Nelson, Paul Woodward, Catherine Lovekin, Thomas Masser, Chris L. Fryer, P. Ramaprabhu, Marianne Francois, Gabriel Rockefeller, "Cross-code Comparisons of mixing during the implosion of dense cylindrical and spherical shells", J. Comp. Phys., 275, 154-173, 2013.
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 Supernovae", 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.
Cody Raskin, Daniel Kasen, Rainer Moll, Josiah Schwab and S.E. Woosley, "Type Ia Supernovae from Merging White Dwarfs. II. Post-Merger Detonations," Astrophysical Journal, 788, No. 1, 2014.
R. Moll, C. Raskin, D. Kasen and S.E. Woosley, "Type Ia Supernovae from Merging White Dwarfs. I. Prompt Detonations," Astrophysical Journal, 785, No. 2, 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]
Lucille H. Frey, Wesley Even, Daniel J. Whalen, Chris L. Fryer, Aimee L. Hungerford, Christopher J. Fontes, James Colgan, "The Los Alamos Supernova Light-Curve Project: Computational Methods," Astrophysical Journal Supp. Series, 204, No. 2, 2013.
Daniel J. Whalen, Wesley Even, Joseph Smidt, Alexander Heger, K.-J. Chen, Chris L. Fryer, Massimo Stiavelli, Hao Xu, Candace C. Joggerst, "Supermassive Population III Supernovae and the Birth of the First Quasars," Astrophysical Journal, 778, No. 1, 2013.
Jarrett L. Johnson, Daniel J. Whalen, Wesley Even, Chris L. Fryer, Alex Heger, Joseph Smidt, Ke-Jung Chen, "The Biggest Explosions in the Universe," Astrophysical Journal, 775, No. 2, 2013.
R. Moll and S.E. Woosley, "Multi-Dimensional Models for Double Detonation in Sub-Chandrasekhar Mass White Dwarfs," Astrophysical Journal, 774, No. 2, 2013.
Elizabeth Lovegrove and S.E. Woosley, "Very Low Energy Supernovae From Neutrino Mass Loss," Astrophysical Journal, 769, No. 2, 2013.
Daniel J. Whalen, Wesley Even, C.C. Lovekin, Chris L. Fryer, Massimo Stiavelli, P.W.A. Roming, Jeff Cooke, T.A. Pritchard, Daniel E. Holz, Cynthia Knight, "Illuminating the Primeval Universe with Type IIn Supernovae," Astrophysical Journal, 768, No. 2, 2013.
Daniel J. Whalen, Candace C. Joggerst, Chris L. Fryer, Massimo Stiavelli, Alexander Heger, Daniel E. Holz, "Finding the First Cosmic Explosions. II. Core-Collapse Supernovae," Astrophysical Journal, 768, No. 1, 2013.
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, 771, 58, 2013.
Adam Burrows, Joshua C. Dolence, Jeremiah W. Murphy, "An Investigation into the Character of Pre-Explosion Core-Collapse Supernova Shock Motion," [arxiv].
J. Nordhaus, T.D. Brandt, A. Burrows, A. Almgren, "The Hydrodynamic Origin of Neutron Star Kicks", Monthly Notices of the Royal Astronomical Society, 423, Issue 2, June 2012.
Ke-Jung Chen, Alexander Heger and Ann Almgren, "Conservative Initial Mapping for Multidimensional Simulations of Stellar Explosions", Journal of Physics: Conference Series, 402, Conf. 1, 2012.
R. Sisneros, C. Malone, A. Nonaka, and S. Woosley, "Investigation of Turbulence in the Early Stages of a High Resolution Supernova Simulation", Proceedings of the Supercomputing 2012 Conference.
A. J. Aspden, J. B. Bell, S. Dong, and S. E. Woosley, "Burning Thermals in Type Ia Supernovae", Astrophysical Journal, 738, 94-107, (2011). [ApJ]
Ke-Jung Chen, Alexander Heger, and Ann Almgren, "Multidimensional Simulations of Pair-Instability Supernovae", Computer Physics Communications, 182:1, 254-256, January 2011. [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].
Adam Burrows, Jason Nordhaus, Ann Almgren, John Bell, "The potential role of spatial dimension in the neutrino-driving mechanism of core-collapse supernova explosions,", Comp. Phys. Comm., 182, Issue 9, pp. 1764-1766, 2011.
Ke-Jung Chen, Alexander Heger, Ann Almgren, and Daniel Whalen, "Pair Instability Supernovae of Very Massive Population III Stars" Astrophysical Journal, 792, 44, 2014. [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.