

Home Astronomy research Software Infrastructure: MESA FLASH STARLIB My codes White dwarf supernova: Remnant metallicities Colliding white dwarfs Merging white dwarfs Ignition conditions Metallicity effects Central density effects Detonation density effects Tracer particle burning Subsonic burning fronts Supersonic burning fronts W7 profiles Massive star supernova: Rotating progenitors 3D evolution ^{26}Al & ^{60}Fe ^{44}Ti, ^{60}Co & ^{56}Ni Yields of radionuclides Effects of ^{12}C +^{12}C SN 1987A light curve Constraints on Ni/Fe ratios An rprocess Compact object IMF Stars: Pulsating white dwarfs Pop III with JWST Monte Carlo massive stars Neutrinos from preSN PreSN variations Monte Carlo white dwarfs SAGB stars Classical novae He shell convection Presolar grains He burn on neutron stars BBFH at 40 years Chemical Evolution: Hypatia catalog Zone models H to Zn Mixing ejecta γrays within 100 Mpc Thermodynamics & Networks Stellar EOS 12C(α,γ)16O Rate Protonrich NSE Reaction networks Bayesian reaction rates Verification Problems: Validating an astro code SuOlson Cog8 Mader RMTV Sedov Noh Software instruments Presentations Illustrations Videos Bicycle adventures AAS Journals 2019 JINA Rprocess Workshop 2019 MESA Marketplace 2019 MESA Summer School 2019 AST111 Earned Admission Teaching materials Education and Public Outreach Contact: F.X.Timmes my one page vitae, full vitae, research statement, and teaching statement. 
Integration of Nuclear Reaction Networks for Stellar Hydrodynamics
In this paper, methods for solving the stiff system of ordinary differential equations that constitute nuclear reaction networks are surveyed. Three semiimplicit time integration algorithms are examined; a traditional firstorderaccurate Euler method, a fourthorderaccurate KapsRentrop method, and a variableorder BaderDeuflhard method. These three integration methods are coupled to eight different linear algebra packages. Four of the linear algebra packages operate on dense matrices (LAPACK, LUDCMP, LEQS, GIFT), three of them are designed for the direct solution of sparse matrices (MA28, UMFPACK, Y12M), and one uses an iterative method for sparse matrices (BiCG). The scaling properties and behavior of the 24 combinations (3 time integration methods times 8 linear algebra packages) are analyzed by running each combination on seven different nuclear reaction networks.
An Inexpensive Nuclear Energy Generation Network for Stellar Hydrodynamics In this paper, we compare the nuclear energy generation rate and abundance levels given by an αchain nuclear reaction network that contains only7 isotopes with a standard 13 isotope αchain reaction network. The energy generation rate of these two small networks are also compared to the energy generation rate given by a 489 isotope reaction network with weak reactions turned on and off. The comparison between the seven isotope and αchain reaction networks indicate the extent to which one can be replaced by the other, and the comparison with the 489 isotope reaction network roughly indicates under what physical conditions it is safe to use the seven isotope and αchain reaction networks.
Nuclear Reaction Network Codes Open source codes are avaliable from this link. 


