*
Cococubed.com


On Integrating Nuclear Reaction Networks




Home

Teaching materials
Astronomy research
Astronomy codes
... Stellar equation of states
... EOS with ionization
... EOS for supernovae
... Chemical potentials
... Stellar atmospheres
... Voigt Function
... Polytropic stars
... Cold white dwarfs
... Hotter white dwarfs
... Cold neutron stars
... Stellar opacities
... Neutrino energy loss rates
... Ephemeris routines
... Fermi-Dirac functions
... Galactic chemical evolution

... Nuclear reaction networks
... Nuclear statistical equilibrium
... Laminar deflagrations
... CJ detonations
... ZND detonations
... Fitting to conic sections
... Unusual linear algebra
... Derivatives on uneven grids
... Pentadiagonal solver
... Quadratics, Cubics, Quartics
... Supernova light curves
... Exact Riemann solutions
... 1D PPM Hydrodynamics
... Verification problems

... EZ stellar evolution
... FLASH code
... Mesa code
Astronomy talks
Astronomy images
Outreach
Family album
Bicycle adventures
Artwork

Contact us:
J.D. Maldonado
F.X.Timmes, my vitae
This journal article discusses these aspects of integrating nuclear reaction networks, while "An Inexpensive Nuclear Energy Generation Network For Stellar Hydrodynamics" discusses the relative correctness of the answers generated. Code for various thermonuclear reaction networks is given here.


image
13 isotopes alpha-chain 		Jacobian matrix for a 13 isotope reaction network at a time shortly after the onset of hydrostatic helium burning. image
19 isotopes alpha-chain 		Jacobian matrix for the 19 isotope reaction network at a time shortly after the onset of helium burning.
image
47 isotopes 		Jacobian matrix for the 47 isotope reaction network at a time shortly after the onset of helium burning. image
76 isotopes 		Jacobian matrix for the 76 isotope reaction network at a time shortly after the onset of a hydrostatic helium burning.
image
127 isotopes 		Jacobian matrix for the 127 isotope reaction network at a time shortly after the onset of hydrostatic helium burning. image
200 isotopes 		Jacobian matrix for the 200 isotope reaction network at a time shortly after the onset of helium burning.
image
489 isotopes 		Jacobian matrix for the 489 isotope reaction network at a time shortly after the onset of helium burning. image
CPU of Euler method 		CPU time of the linear algebra for a single time step with the first-order accurate Euler method.
image
CPU of Kaps-Rentrop 		CPU time of the linear algebra for a single, successful time step with the fourth-order accurate Kaps-Rentrop method. image
CPU of Bader-Deuflhard 		Minimum CPU time of the linear algebra for a single successful time step with the variable order Bader-Deuflhard method.
image
Cost per decomposition 		CPU time required by the linear algebra for a single, successful time step across the three time integration methods. image
Cost per decomposition 		CPU time required by the linear algebra for a single, successful step across the three time integration methods.
image
Hot helium burn 		Evolution of the abundance levels during hydrostatic helium burning. image
Kaps-Rentrop helium burn 		Total CPU time with the fourth-order Kaps-Rentrop method for the hydrostatic helium burning test case.
image
Bader-Deuflhard burn 		Total CPU time consumed by the variable order Bader-Deuflhard for the hydrostatic helium burning test case. image
Relative efficiency 		Total CPU time consumed by the fourth-order Kaps-Rentrop method relative to the variable order Bader-Deuflhard.
image
CPU for hot helium burn 		Total CPU time consumed by the variable order Bader-Deuflhard method for the hydrostatic helium burning test case.