A posteriori error estimator and AMR for discrete ordinates nodal transport methods

Jose I. Duo, Yousry Y. Azmy, Ludmil Tomov Zikatanov

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

In the development of high fidelity transport solvers, optimization of the use of available computational resources and access to a tool for assessing quality of the solution are key to the success of large-scale nuclear systems' simulation. In this regard, error control provides the analyst with a confidence level in the numerical solution and enables for optimization of resources through Adaptive Mesh Refinement (AMR). In this paper, we derive an a posteriori error estimator based on the nodal solution of the Arbitrarily High Order Transport Method of the Nodal type (AHOT-N). Furthermore, by making assumptions on the regularity of the solution, we represent the error estimator as a function of computable volume and element-edges residuals. The global L2 error norm is proved to be bound by the estimator. To lighten the computational load, we present a numerical approximation to the aforementioned residuals and split the global norm error estimator into local error indicators. These indicators are used to drive an AMR strategy for the spatial discretization. However, the indicators based on forward solution residuals alone do not bound the cell-wise error. The estimator and AMR strategy are tested in two problems featuring strong heterogeneity and highly transport streaming regime with strong flux gradients. The results show that the error estimator indeed bounds the global error norms and that the error indicator follows the cell-error's spatial distribution pattern closely. The AMR strategy proves beneficial to optimize resources, primarily by reducing the number of unknowns solved for to achieve prescribed solution accuracy in global L2 error norm. Likewise, AMR achieves higher accuracy compared to uniform refinement when resolving sharp flux gradients, for the same number of unknowns.

Original languageEnglish (US)
Pages (from-to)268-273
Number of pages6
JournalAnnals of Nuclear Energy
Volume36
Issue number3
DOIs
StatePublished - Apr 1 2009

Fingerprint

Fluxes
Spatial distribution

All Science Journal Classification (ASJC) codes

  • Nuclear Energy and Engineering

Cite this

@article{80cf089bfeac4af3970eb1b2e818cdc4,
title = "A posteriori error estimator and AMR for discrete ordinates nodal transport methods",
abstract = "In the development of high fidelity transport solvers, optimization of the use of available computational resources and access to a tool for assessing quality of the solution are key to the success of large-scale nuclear systems' simulation. In this regard, error control provides the analyst with a confidence level in the numerical solution and enables for optimization of resources through Adaptive Mesh Refinement (AMR). In this paper, we derive an a posteriori error estimator based on the nodal solution of the Arbitrarily High Order Transport Method of the Nodal type (AHOT-N). Furthermore, by making assumptions on the regularity of the solution, we represent the error estimator as a function of computable volume and element-edges residuals. The global L2 error norm is proved to be bound by the estimator. To lighten the computational load, we present a numerical approximation to the aforementioned residuals and split the global norm error estimator into local error indicators. These indicators are used to drive an AMR strategy for the spatial discretization. However, the indicators based on forward solution residuals alone do not bound the cell-wise error. The estimator and AMR strategy are tested in two problems featuring strong heterogeneity and highly transport streaming regime with strong flux gradients. The results show that the error estimator indeed bounds the global error norms and that the error indicator follows the cell-error's spatial distribution pattern closely. The AMR strategy proves beneficial to optimize resources, primarily by reducing the number of unknowns solved for to achieve prescribed solution accuracy in global L2 error norm. Likewise, AMR achieves higher accuracy compared to uniform refinement when resolving sharp flux gradients, for the same number of unknowns.",
author = "Duo, {Jose I.} and Azmy, {Yousry Y.} and Zikatanov, {Ludmil Tomov}",
year = "2009",
month = "4",
day = "1",
doi = "10.1016/j.anucene.2008.12.008",
language = "English (US)",
volume = "36",
pages = "268--273",
journal = "Annals of Nuclear Energy",
issn = "0306-4549",
publisher = "Elsevier Limited",
number = "3",

}

A posteriori error estimator and AMR for discrete ordinates nodal transport methods. / Duo, Jose I.; Azmy, Yousry Y.; Zikatanov, Ludmil Tomov.

In: Annals of Nuclear Energy, Vol. 36, No. 3, 01.04.2009, p. 268-273.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A posteriori error estimator and AMR for discrete ordinates nodal transport methods

AU - Duo, Jose I.

AU - Azmy, Yousry Y.

AU - Zikatanov, Ludmil Tomov

PY - 2009/4/1

Y1 - 2009/4/1

N2 - In the development of high fidelity transport solvers, optimization of the use of available computational resources and access to a tool for assessing quality of the solution are key to the success of large-scale nuclear systems' simulation. In this regard, error control provides the analyst with a confidence level in the numerical solution and enables for optimization of resources through Adaptive Mesh Refinement (AMR). In this paper, we derive an a posteriori error estimator based on the nodal solution of the Arbitrarily High Order Transport Method of the Nodal type (AHOT-N). Furthermore, by making assumptions on the regularity of the solution, we represent the error estimator as a function of computable volume and element-edges residuals. The global L2 error norm is proved to be bound by the estimator. To lighten the computational load, we present a numerical approximation to the aforementioned residuals and split the global norm error estimator into local error indicators. These indicators are used to drive an AMR strategy for the spatial discretization. However, the indicators based on forward solution residuals alone do not bound the cell-wise error. The estimator and AMR strategy are tested in two problems featuring strong heterogeneity and highly transport streaming regime with strong flux gradients. The results show that the error estimator indeed bounds the global error norms and that the error indicator follows the cell-error's spatial distribution pattern closely. The AMR strategy proves beneficial to optimize resources, primarily by reducing the number of unknowns solved for to achieve prescribed solution accuracy in global L2 error norm. Likewise, AMR achieves higher accuracy compared to uniform refinement when resolving sharp flux gradients, for the same number of unknowns.

AB - In the development of high fidelity transport solvers, optimization of the use of available computational resources and access to a tool for assessing quality of the solution are key to the success of large-scale nuclear systems' simulation. In this regard, error control provides the analyst with a confidence level in the numerical solution and enables for optimization of resources through Adaptive Mesh Refinement (AMR). In this paper, we derive an a posteriori error estimator based on the nodal solution of the Arbitrarily High Order Transport Method of the Nodal type (AHOT-N). Furthermore, by making assumptions on the regularity of the solution, we represent the error estimator as a function of computable volume and element-edges residuals. The global L2 error norm is proved to be bound by the estimator. To lighten the computational load, we present a numerical approximation to the aforementioned residuals and split the global norm error estimator into local error indicators. These indicators are used to drive an AMR strategy for the spatial discretization. However, the indicators based on forward solution residuals alone do not bound the cell-wise error. The estimator and AMR strategy are tested in two problems featuring strong heterogeneity and highly transport streaming regime with strong flux gradients. The results show that the error estimator indeed bounds the global error norms and that the error indicator follows the cell-error's spatial distribution pattern closely. The AMR strategy proves beneficial to optimize resources, primarily by reducing the number of unknowns solved for to achieve prescribed solution accuracy in global L2 error norm. Likewise, AMR achieves higher accuracy compared to uniform refinement when resolving sharp flux gradients, for the same number of unknowns.

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

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

U2 - 10.1016/j.anucene.2008.12.008

DO - 10.1016/j.anucene.2008.12.008

M3 - Article

AN - SCOPUS:61349157276

VL - 36

SP - 268

EP - 273

JO - Annals of Nuclear Energy

JF - Annals of Nuclear Energy

SN - 0306-4549

IS - 3

ER -