Integration of nuclear reaction networks for stellar hydrodynamics

Research output: Contribution to journalArticle

123 Citations (Scopus)

Abstract

Methods for solving the stiff system of ordinary differential equations that constitute nuclear reaction networks are surveyed. Three semi-implicit time integration algorithms are examined; a traditional first-order-accurate Euler method, a fourth-order-accurate Kaps-Rentrop method, and a variable-order Bader-Deuflhard 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. These reaction networks range from a hardwired 13 isotope a-chain and heavy-ion reaction network, which is suitable for most multidimensional simulations of stellar phenomena, to a 489 isotope reaction network, which is suitable for determining the yields of isotopes lighter than technetium in spherically symmetric models of Type II supernovae. Each of the time integration methods and linear algebra packages are capable of generating accurate results, but the efficiency of the various methods - evaluated across several different machine architectures and compiler options - differ dramatically. If the execution speed of reaction networks that contain less than about 50 isotopes is an overriding concern, then the variable-order Bader-Deuflhard time integration method coupled with routines generated from the GIFT matrix package or LAPACK with vendor-optimized BLAS routines is a good choice. If the amount of storage needed for any reaction network is a concern, then any of the sparse matrix packages will reduce the storage costs by 70%-90%. When a balance between accuracy, overall efficiency, and ease of use is desirable, then the variable-order Bader-Deuflhard time integration method coupled with the MA28 sparse matrix package is a strong choice.

Original languageEnglish (US)
Pages (from-to)241-263
Number of pages23
JournalAstrophysical Journal, Supplement Series
Volume124
Issue number1
DOIs
StatePublished - Sep 1999
Externally publishedYes

Fingerprint

nuclear reactions
hydrodynamics
algebra
isotopes
matrix
isotope
compilers
technetium
matrices
method
nuclear reaction
supernovae
heavy ions
differential equations
costs
scaling
simulation
ion

Keywords

  • Hydrodynamics
  • Methods: numerical
  • Nuclear reactions, nucleosynthesis, abundances
  • Stars: interiors

ASJC Scopus subject areas

  • Space and Planetary Science

Cite this

Integration of nuclear reaction networks for stellar hydrodynamics. / Timmes, Francis.

In: Astrophysical Journal, Supplement Series, Vol. 124, No. 1, 09.1999, p. 241-263.

Research output: Contribution to journalArticle

@article{da905f13972b49019e02ab45a6ab9789,
title = "Integration of nuclear reaction networks for stellar hydrodynamics",
abstract = "Methods for solving the stiff system of ordinary differential equations that constitute nuclear reaction networks are surveyed. Three semi-implicit time integration algorithms are examined; a traditional first-order-accurate Euler method, a fourth-order-accurate Kaps-Rentrop method, and a variable-order Bader-Deuflhard 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. These reaction networks range from a hardwired 13 isotope a-chain and heavy-ion reaction network, which is suitable for most multidimensional simulations of stellar phenomena, to a 489 isotope reaction network, which is suitable for determining the yields of isotopes lighter than technetium in spherically symmetric models of Type II supernovae. Each of the time integration methods and linear algebra packages are capable of generating accurate results, but the efficiency of the various methods - evaluated across several different machine architectures and compiler options - differ dramatically. If the execution speed of reaction networks that contain less than about 50 isotopes is an overriding concern, then the variable-order Bader-Deuflhard time integration method coupled with routines generated from the GIFT matrix package or LAPACK with vendor-optimized BLAS routines is a good choice. If the amount of storage needed for any reaction network is a concern, then any of the sparse matrix packages will reduce the storage costs by 70{\%}-90{\%}. When a balance between accuracy, overall efficiency, and ease of use is desirable, then the variable-order Bader-Deuflhard time integration method coupled with the MA28 sparse matrix package is a strong choice.",
keywords = "Hydrodynamics, Methods: numerical, Nuclear reactions, nucleosynthesis, abundances, Stars: interiors",
author = "Francis Timmes",
year = "1999",
month = "9",
doi = "10.1086/313257",
language = "English (US)",
volume = "124",
pages = "241--263",
journal = "Astrophysical Journal, Supplement Series",
issn = "0067-0049",
publisher = "IOP Publishing Ltd.",
number = "1",

}

TY - JOUR

T1 - Integration of nuclear reaction networks for stellar hydrodynamics

AU - Timmes, Francis

PY - 1999/9

Y1 - 1999/9

N2 - Methods for solving the stiff system of ordinary differential equations that constitute nuclear reaction networks are surveyed. Three semi-implicit time integration algorithms are examined; a traditional first-order-accurate Euler method, a fourth-order-accurate Kaps-Rentrop method, and a variable-order Bader-Deuflhard 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. These reaction networks range from a hardwired 13 isotope a-chain and heavy-ion reaction network, which is suitable for most multidimensional simulations of stellar phenomena, to a 489 isotope reaction network, which is suitable for determining the yields of isotopes lighter than technetium in spherically symmetric models of Type II supernovae. Each of the time integration methods and linear algebra packages are capable of generating accurate results, but the efficiency of the various methods - evaluated across several different machine architectures and compiler options - differ dramatically. If the execution speed of reaction networks that contain less than about 50 isotopes is an overriding concern, then the variable-order Bader-Deuflhard time integration method coupled with routines generated from the GIFT matrix package or LAPACK with vendor-optimized BLAS routines is a good choice. If the amount of storage needed for any reaction network is a concern, then any of the sparse matrix packages will reduce the storage costs by 70%-90%. When a balance between accuracy, overall efficiency, and ease of use is desirable, then the variable-order Bader-Deuflhard time integration method coupled with the MA28 sparse matrix package is a strong choice.

AB - Methods for solving the stiff system of ordinary differential equations that constitute nuclear reaction networks are surveyed. Three semi-implicit time integration algorithms are examined; a traditional first-order-accurate Euler method, a fourth-order-accurate Kaps-Rentrop method, and a variable-order Bader-Deuflhard 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. These reaction networks range from a hardwired 13 isotope a-chain and heavy-ion reaction network, which is suitable for most multidimensional simulations of stellar phenomena, to a 489 isotope reaction network, which is suitable for determining the yields of isotopes lighter than technetium in spherically symmetric models of Type II supernovae. Each of the time integration methods and linear algebra packages are capable of generating accurate results, but the efficiency of the various methods - evaluated across several different machine architectures and compiler options - differ dramatically. If the execution speed of reaction networks that contain less than about 50 isotopes is an overriding concern, then the variable-order Bader-Deuflhard time integration method coupled with routines generated from the GIFT matrix package or LAPACK with vendor-optimized BLAS routines is a good choice. If the amount of storage needed for any reaction network is a concern, then any of the sparse matrix packages will reduce the storage costs by 70%-90%. When a balance between accuracy, overall efficiency, and ease of use is desirable, then the variable-order Bader-Deuflhard time integration method coupled with the MA28 sparse matrix package is a strong choice.

KW - Hydrodynamics

KW - Methods: numerical

KW - Nuclear reactions, nucleosynthesis, abundances

KW - Stars: interiors

UR - http://www.scopus.com/inward/record.url?scp=0002035985&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0002035985&partnerID=8YFLogxK

U2 - 10.1086/313257

DO - 10.1086/313257

M3 - Article

VL - 124

SP - 241

EP - 263

JO - Astrophysical Journal, Supplement Series

JF - Astrophysical Journal, Supplement Series

SN - 0067-0049

IS - 1

ER -