Improving the accuracy of the Chebyshev rational approximation method using substeps

Aarno E. Isotalo, Maria Pusa

Research output: Contribution to journalArticleScientificpeer-review

6 Citations (Scopus)

Abstract

The Chebyshev rational approximation method (CRAM) for solving the decay and depletion of nuclides is shown to have a remarkable decrease in error when advancing the system with the same time step and microscopic reaction rates as the previous step. This property is exploited here to achieve high accuracy in any end-of-step solution by dividing a step into equidistant substeps. The computational cost of identical substeps can be reduced significantly below that of an equal number of regular steps, as the lower-upper decompositions for the linear solutions required in CRAM need to be formed only on the first substep. The improved accuracy provided by substeps is most relevant in decay calculations, where there have previously been concerns about the accuracy and generality of CRAM. With substeps, CRAM can solve any decay or depletion problem with constant microscopic reaction rates to an extremely high accuracy for all nuclides with concentrations above an arbitrary limit.
Original languageEnglish
Pages (from-to)65-77
JournalNuclear Science and Engineering
Volume183
Issue number1
DOIs
Publication statusPublished - 2016
MoE publication typeA1 Journal article-refereed

Fingerprint

Reaction rates
Isotopes
Decomposition
Costs

Keywords

  • CRAM
  • decay calculations
  • substeps

Cite this

Isotalo, Aarno E. ; Pusa, Maria. / Improving the accuracy of the Chebyshev rational approximation method using substeps. In: Nuclear Science and Engineering. 2016 ; Vol. 183, No. 1. pp. 65-77.
@article{6e22f256d9e540d282dbae926a13255b,
title = "Improving the accuracy of the Chebyshev rational approximation method using substeps",
abstract = "The Chebyshev rational approximation method (CRAM) for solving the decay and depletion of nuclides is shown to have a remarkable decrease in error when advancing the system with the same time step and microscopic reaction rates as the previous step. This property is exploited here to achieve high accuracy in any end-of-step solution by dividing a step into equidistant substeps. The computational cost of identical substeps can be reduced significantly below that of an equal number of regular steps, as the lower-upper decompositions for the linear solutions required in CRAM need to be formed only on the first substep. The improved accuracy provided by substeps is most relevant in decay calculations, where there have previously been concerns about the accuracy and generality of CRAM. With substeps, CRAM can solve any decay or depletion problem with constant microscopic reaction rates to an extremely high accuracy for all nuclides with concentrations above an arbitrary limit.",
keywords = "CRAM, decay calculations, substeps",
author = "Isotalo, {Aarno E.} and Maria Pusa",
year = "2016",
doi = "10.13182/NSE15-67",
language = "English",
volume = "183",
pages = "65--77",
journal = "Nuclear Science and Engineering",
issn = "0029-5639",
publisher = "American Nuclear Society ANS",
number = "1",

}

Improving the accuracy of the Chebyshev rational approximation method using substeps. / Isotalo, Aarno E.; Pusa, Maria.

In: Nuclear Science and Engineering, Vol. 183, No. 1, 2016, p. 65-77.

Research output: Contribution to journalArticleScientificpeer-review

TY - JOUR

T1 - Improving the accuracy of the Chebyshev rational approximation method using substeps

AU - Isotalo, Aarno E.

AU - Pusa, Maria

PY - 2016

Y1 - 2016

N2 - The Chebyshev rational approximation method (CRAM) for solving the decay and depletion of nuclides is shown to have a remarkable decrease in error when advancing the system with the same time step and microscopic reaction rates as the previous step. This property is exploited here to achieve high accuracy in any end-of-step solution by dividing a step into equidistant substeps. The computational cost of identical substeps can be reduced significantly below that of an equal number of regular steps, as the lower-upper decompositions for the linear solutions required in CRAM need to be formed only on the first substep. The improved accuracy provided by substeps is most relevant in decay calculations, where there have previously been concerns about the accuracy and generality of CRAM. With substeps, CRAM can solve any decay or depletion problem with constant microscopic reaction rates to an extremely high accuracy for all nuclides with concentrations above an arbitrary limit.

AB - The Chebyshev rational approximation method (CRAM) for solving the decay and depletion of nuclides is shown to have a remarkable decrease in error when advancing the system with the same time step and microscopic reaction rates as the previous step. This property is exploited here to achieve high accuracy in any end-of-step solution by dividing a step into equidistant substeps. The computational cost of identical substeps can be reduced significantly below that of an equal number of regular steps, as the lower-upper decompositions for the linear solutions required in CRAM need to be formed only on the first substep. The improved accuracy provided by substeps is most relevant in decay calculations, where there have previously been concerns about the accuracy and generality of CRAM. With substeps, CRAM can solve any decay or depletion problem with constant microscopic reaction rates to an extremely high accuracy for all nuclides with concentrations above an arbitrary limit.

KW - CRAM

KW - decay calculations

KW - substeps

U2 - 10.13182/NSE15-67

DO - 10.13182/NSE15-67

M3 - Article

VL - 183

SP - 65

EP - 77

JO - Nuclear Science and Engineering

JF - Nuclear Science and Engineering

SN - 0029-5639

IS - 1

ER -