## Compressible Astrophysics: CASTRO
## CASTRO ContributorsWe would like to gratefully acknowledge the direct contributions to the CASTRO code by:
Mike Zingale, Stonybrook University
Max Katz, Stonybrook University -- multipole boundary conditions
## CASTRO Download
You can find the latest version of the
The ## CASTRO Email Support ListWe 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: http://bender.astro.sunysb.edu/mailman/listinfo/castro-help You can e-mail everyone on the list by sending to castro-help@bender.astro.sunysb.edu 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: http://bender.astro.sunysb.edu/pipermail/castro-help Please encourage anyone who is interested in CASTRO to use this list, as it will make supporting the code easier for everyone. ## CASTRO Papers
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,
W. Zhang, L. Howell, A. Almgren, A. Burrows, and J. Bell
Weiqun Zhang, L. Howell, A. Almgren, A. Burrows, J. Dolence, J. Bell,
## 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,
C. M. Malone, A. Nonaka, S. E. Woosley, A. S. Almgren, J. B. Bell, S. Dong, and M. Zingale,
Ke-Jung Chen, Alexander Heger, Stan Woosley, Ann Almgren, and Daniel Whalen,
Ke-Jung Chen, Alexander Heger, S.E. Woosley, Ann S. Almgren, and Daniel J. Whalen,
Ke-Jung Chen, Alexander Heger, S.E. Woosley, Ann Almgren, and Daniel J. Whalen, and Jarrett L. Johnson,
Ke-Jung Chen, Alexander Heger, and Ann S. Almgren,
K. Chen, A. Heger, S. Woosley, A. Almgren, and W. Zhang,
Haitao Ma, Stan Woosley, Chris Malone, Ann Almgren, and J.B. Bell,
Adam Burrows, Joshua C. Dolence, Jeremiah W. Murphy,
Emmanouela Rantsiou,
A. Burrows,
J. Nordhaus,
Ann Almgren,
C. C. Joggerst, Daniel Whalen,
A. Almgren, J. Bell, D. Kasen, M. Lijewski, A. Nonaka, P. Nugent, C. Rendleman, R. Thomas, M. Zingale,
J. Nordhaus,
A. Burrows,
A. Almgren,
and
J. Bell,
C. C. Joggerst,
A. Almgren, and S. E. Woosley,
C. C. Joggerst,
A. Almgren, J. Bell, Alexander Heger, Daniel Whalen, and S. E. Woosley,
Ke-Jung Chen,
Alexander Heger,
and
Ann Almgren,
Ke-Jung Chen,
Alexander Heger,
and
Ann Almgren,
## Posters Using CASTRO
Ke-Jung Chen,
Alexander Heger, Ann Almgren, and Shuxia Zhang,
A. Almgren, J. Bell, M. Day, L. Howell, C. Joggerst, E. Myra, A. Nonaka,
J. Nordhaus, M. Singer, and M. Zingale,
## 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 StateCASTRO 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 NetworksCASTRO 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 CASTROFor more information about the radiation component of CASTRO, see Paper II and Paper III that describe the implementation and performance of single-group and multigroup radiation in CASTRO. ## Self-gravityCASTRO 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 CASTROOur 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. ## Visualization## 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.## AcknowledgementsThe 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. |