Fourth-Order Comprehensive Adjoint Sensitivity Analysis (4th-CASAM) of Response-Coupled Linear Forward/Adjoint Systems: I. Theoretical Framework
Abstract
:1. Introduction
2. Background: Mathematical Description of the Physical System
3. The Fourth-Order Comprehensive Adjoint Sensitivity Analysis Methodology (4th-CASAM) for Linear Systems
3.1. Computation of the First-Order Sensitivities of
3.1.1. Finite-Difference Approximation Using Re-Computations with User-Modified Parameters
3.1.2. Forward Sensitivity Analysis Methodology (FSAM)
3.1.3. First-Order Comprehensive Adjoint Sensitivity Analysis Methodology (1st-CASAM)
- 1
- Introduce a Hilbert space, denoted as , comprising square-integrable functions vector-valued elements of the form , with and endowed with an inner product between two elements, , , of this Hilbert space, which will be denoted as and defined as follows:
- 2
- In the Hilbert , form the inner product of Equation (36) with a yet undefined vector-valued function to obtain the following relation:
- 3
- Using the definition of the adjoint operator in the Hilbert space , recast the left-side of Equation (43) as follows:
- 4
- Require the first term on right-side of Equation (44) to represent the indirect-effect term defined in Equation (35), to obtain the following relation:
- 5
- Implement the boundary conditions given in Equation (37) into Equation (44) and eliminate the remaining unknown boundary-values of the functions and from the expression of the bilinear concomitant by selecting appropriate boundary conditions for the function , to ensure that Equation (46) is well-posed while being independent of unknown values of , , and . The boundary conditions thus chosen for the function can be represented in operator form as follow
- 6
- The selection of the boundary conditions for the adjoint function represented by Equation (48) eliminates the appearance of the unknown values of in and reduces this bilinear concomitant to a residual quantity that contains boundary terms involving only known values of , , , and . This residual quantity will be denoted as . In general, this residual quantity does not automatically vanish, although it may do so occasionally.
- 7
- The system of equations comprising Equation (46) together with the boundary conditions represented by Equation (48) constitute the 1st-Level Adjoint Sensitivity System (1st-LASS). the solution of the 1st-LASS will be called the 1st-level adjoint function. The 1st-LASS is called “first-level” (as opposed to “first-order”) because it does not contain any differential or functional-derivatives, but its solution will be used below to compute the first-order sensitivities of the response with respect to the model parameters. This terminology will be also used in the sequel, when deriving the expressions for the 2nd- and 3rd-order sensitivities.
- 8
- It follows from Equations (43) and (44) that the following relation holds:
- 9
- Recalling that the first term on the right-side of Equation (49) is, in view of Equation (46), the indirect-effect term , it follows from Equation (49) that the indirect-effect term can be expressed in terms of the 1st-level adjoint function as follows:
3.1.4. Comparison of Computational Requirements for Computing the First-Order Response Sensitivities with Respect to the Model Parameters
- The 1st-Order Comprehensive Adjoint Sensitivity Analysis Methodology (1st-CASAM) requires large-scale computations: one large-scale computation for computing the first-level adjoint state function and one large-scale computation for computing the first-level adjoint state function .
- The Forward Sensitivity Analysis Method (FSAM) requires large-scale computations: large-scale computations for computing the first-order derivatives of the forward state function and large-scale computations for computing the first-order derivatives of the adjoint state function .
- The finite-difference (FD) method requires large-scale computations: large-scale computations for computing the forward state functions and large-scale computations for computing the adjoint state functions .
3.2. Computation of the Second-Order Sensitivities of
3.2.1. Finite-Difference Approximation Using Re-Computations with User-Modified Parameters
3.2.2. Forward Sensitivity Analysis Methodology (FSAM)
3.2.3. The 2nd-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)
- (i)
- The boundary conditions Equation (87) together with the operator Equation (85) constitute a well posed problem for the functions .
- (ii)
- The implementation in Equation (83) of the boundary conditions provided in Equation (77) together with those provided in Equation (87) eliminates all of the unknown values of the functions and in the expression of the bilinear concomitant . The bilinear concomitant may vanish after these boundary conditions are implemented, but if it does not, it will be reduced to a residual quantity which will be denoted as and which will comprise only known values of , , and . In principle, the bilinear concomitant can always be “forced” to vanish by introducing delta-function terms in the definition of the adjoint operator, but such a practice could lead to severe (and usually unnecessary) difficulties when attempting to solve the equations that involve such extended adjoint operators.
3.2.4. Comparison of Computational Requirements for Computing the Second-Order Response Sensitivities with Respect to the Model Parameters
3.3. Computation of the Third-Order Sensitivities of
3.3.1. Finite-Difference Approximation Using Re-Computations with User-Modified Parameters
3.3.2. Forward Sensitivity Analysis Methodology
3.3.3. The 3rd-Order Comprehensive Adjoint Sensitivity Analysis Methodology (3rd-CASAM)
- 1
- Define a Hilbert space, denoted as , comprising vector-valued elements of the form , with The inner product between two elements, and , of this Hilbert space, will be denoted as and is defined as follows:
- 2
- In the Hilbert , form the inner product of Equation (119) with a set of yet undefined vector-valued functions ,, to obtain the following relation:
- 3
- Using the definition of the adjoint operator in the Hilbert space , recast the left-side of Equation (135) as follows:
- 4
- The first term on right-side of Equation (136) is now required to represent the indirect-effect term defined in Equation (109). This requirement is satisfied by imposing the following relation on each element , :
- 5
- The definition of the set of vectors will now be completed by selecting boundary conditions for this set of vectors, which will be represented in operator form as follows:
- (i)
- The boundary conditions Equation (148) together with the operator Equation (138) constitute a well posed problem for the functions .
- (ii)
- Implementation in Equation (136) of the boundary conditions provided in Equation (120) together with those provided in Equation (148) eliminates all of the unknown values of the functions and in the expression of the bilinear concomitant . This bilinear concomitant may vanish after implementing the boundary conditions represented by Equation (148), but if it does not, it will be reduced to a residual quantity which will be denoted as and which will comprise only known values of , , and .
3.3.4. Comparison of Computational Requirements for Computing the Third-Order Response Sensitivities with Respect to Model Parameters
3.4. Computation of the Fourth-Order Sensitivities of
3.4.1. Finite-Difference Approximation Using Re-Computations with User-Modified Parameters
3.4.2. Forward Sensitivity Analysis Methodology
3.4.3. The Fourth-Order Adjoint Sensitivity Analysis Methodology (4th-CASAM)
- (i)
- The boundary conditions Equation (192) together with the operator Equation (189) constitute a well posed problem for the functions .
- (ii)
- Implementation in Equation (187) of the boundary conditions (for the 3rd-LVSSS) provided in Equation (148) together with those provided in Equation (192) eliminates all of the unknown values of the functions and in the expression of the bilinear concomitant . This bilinear concomitant may vanish after implementing the boundary conditions represented by Equation (148), but if it does not, it will be reduced to a residual quantity which will be denoted as and which will comprise only known values of , , and .
4. Discussion and Conclusions
Funding
Conflicts of Interest
References
- Lewins, J. IMPORTANCE: The Adjoint Function; Pergamon Press Ltd.: Oxford, UK, 1965. [Google Scholar]
- Stacey, W.M. Variational Methods in Nuclear Reactor Physics; Academic Press, Inc.: New York, NY, USA, 1974. [Google Scholar]
- Bellman, R.E. Dynamic Programming; Rand Corporation: Santa Monica, CA, USA; Princeton University Press: Princeton, NJ, USA, 1957; ISBN 978-0-691-07951-6. [Google Scholar]
- Cacuci, D.G. Sensitivity theory for nonlinear systems: I. Nonlinear functional analysis approach. J. Math. Phys. 1981, 22, 2794–2802. [Google Scholar] [CrossRef]
- Cacuci, D.G. Sensitivity theory for nonlinear systems: II. Extensions to additional classes of responses. J. Math. Phys. 1981, 22, 2803–2812. [Google Scholar] [CrossRef]
- Cacuci, D.G. Sensitivity and Uncertainty Analysis: Theory, Volume 1; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
- Cacuci, D.G. Second-order adjoint sensitivity analysis methodology (2nd-ASAM) for computing exactly and efficiently first- and second-order sensitivities in large-scale linear systems: I. Computational methodology. J. Comput. Phys. 2015, 284, 687–699. [Google Scholar] [CrossRef] [Green Version]
- Cacuci, D.G. Second-order adjoint sensitivity analysis methodology (2nd-ASAM) for large-scale nonlinear systems: I. Theory. Nucl. Sci. Eng. 2016, 184, 16–30. [Google Scholar] [CrossRef]
- Cacuci, D.G. The Second-Order Adjoint Sensitivity Analysis Methodology; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2018. [Google Scholar]
- Cacuci, D.G.; Fang, R.; Favorite, J.A. Comprehensive second-order adjoint sensitivity analysis methodology (2nd-ASAM) applied to a subcritical experimental reactor physics benchmark: I. Effects of imprecisely known microscopic total and capture cross sections. Energies 2019, 12, 4219. [Google Scholar] [CrossRef] [Green Version]
- Fang, R.; Cacuci, D.G. Comprehensive second-order adjoint sensitivity analysis methodology (2nd-ASAM) applied to a subcritical experimental reactor physics benchmark: II. Effects of imprecisely known microscopic scattering cross sections. Energies 2019, 12, 4114. [Google Scholar] [CrossRef] [Green Version]
- Cacuci, D.G.; Fang, R.; Favorite, J.A.; Badea, M.C.; Di Rocco, F. Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark: III. Effects of Imprecisely Known Microscopic Fission Cross Sections and Average Number of Neutrons per Fission. Energies 2019, 12, 4100. [Google Scholar] [CrossRef] [Green Version]
- Fang, R.; Cacuci, D.G. Comprehensive second-order adjoint sensitivity analysis methodology (2nd-ASAM) applied to a subcritical experimental reactor physics benchmark: IV. Effects of imprecisely known source parameters. Energies 2020, 13, 1431. [Google Scholar] [CrossRef] [Green Version]
- Fang, R.; Cacuci, D.G. Comprehensive second-order adjoint sensitivity analysis methodology (2nd-ASAM) applied to a subcritical experimental reactor physics benchmark: V. Computation of mixed 2nd-order sensitivities involving isotopic number densities. Energies 2020, 13, 2580. [Google Scholar] [CrossRef]
- Cacuci, D.G.; Fang, R.; Favorite, J.A. Comprehensive second-order adjoint sensitivity analysis methodology (2nd-ASAM) applied to a subcritical experimental reactor physics benchmark: VI. Overall impact of 1st-and 2nd-order sensitivities on response uncertainties. Energies 2020, 13, 1674. [Google Scholar] [CrossRef] [Green Version]
- Valentine, T.E. Polyethylene-Reflected Plutonium Metal Sphere Subcritical Noise Measurements, SUB-PU-METMIXED-001. In International Handbook of Evaluated Criticality Safety Benchmark Experiments, NEA/NSC/DOC(95)03/I-IX; Organization for Economic Co-Operation and Development: Paris, France; Nuclear Energy Agency: Paris, France, 2006. [Google Scholar]
- Alcouffe, R.E.; Baker, R.S.; Dahl, J.A.; Turner, S.A.; Ward, R. PARTISN: A Time-Dependent, Parallel Neutral Particle Transport Code System; Los Alamos National Laboratory: New Mexico, NM, USA, 2008; LA-UR-08-07258 (Revised November 2008). [Google Scholar]
- Cacuci, D.G. Third order adjoint sensitivity analysis of reaction rate responses in a multiplying nuclear system with source. Ann. Nucl. Energy 2021, 151, 107924. [Google Scholar] [CrossRef]
- Cacuci, D.G.; Fang, R. Third-order adjoint sensitivity analysis of an OECD/NEA reactor physics benchmark: I. Mathematical framework. Am. J. Comput. Math. (AJCM) 2020, 10, 503–528. [Google Scholar] [CrossRef]
- Fang, R.; Cacuci, D.G. Third-Order Adjoint Sensitivity Analysis of an OECD/NEA Reactor Physics Benchmark: II Computed Sensitivities. Am. J. Comput. Math. (AJCM) 2020, 10, 529–558. [Google Scholar] [CrossRef]
- Fang, R.; Cacuci, D.G. Third-order adjoint sensitivity analysis of an OECD/NEA reactor physics benchmark: III Response moments. Am. J. Comput. Math. (AJCM) 2020, 10, 559–570. [Google Scholar] [CrossRef]
- Cacuci, D.G. BERRU Predictive Modeling: Best Estimate Results with Reduced Uncertainties; Springer: Heidelberg, Germany; New York, NY, USA, 2018. [Google Scholar]
- Cacuci, D.G. Inverse predictive modeling of radiation transport through optically thick media in the presence of counting uncertainties. Nucl. Sci. Eng. 2017, 186, 199–223. [Google Scholar] [CrossRef]
- Cacuci, D.G.; Fang, R.; Badea, M.C. MULTI-PRED: A software module for predictive modeling of coupled multi-physics systems. Nucl. Sci. Eng. 2018, 191, 187–202. [Google Scholar] [CrossRef]
- Cacuci, D.G. Second-order sensitivities of a general functional of the forward and adjoint fluxes in a multiplying nuclear system with source. Nucl. Eng. Des. 2019, 344, 83–106. [Google Scholar] [CrossRef]
- Cacuci, D.G. Application of the second-order comprehensive adjoint sensitivity analysis methodology to compute 1st-and 2nd-order sensitivities of flux functionals in a multiplying system with source. Nucl. Sci. Eng. 2019, 193, 555–600. [Google Scholar] [CrossRef] [Green Version]
- Cacuci, D.G. The Roussopoulos and Schwinger functionals for nuclear systems involving imprecisely known fluxes and/or parameters: Distinctions and equivalences. Nucl. Sci. Eng. 2019, 193, 681–721. [Google Scholar] [CrossRef] [Green Version]
- Cacuci, D.G. Nuclear Thermal-Hydraulics Applications Illustrating the Key Roles of Adjoint-Computed Sensitivities for Overcoming the Curse of Dimensionality in Sensitivity Analysis, Uncertainty Quantification and Predictive Modeling. Nucl. Eng. Des. 2019, 351, 20–32. [Google Scholar] [CrossRef]
- Cacuci, D.G. Second-order adjoint sensitivity analysis of ratios of functionals of the forward and adjoint fluxes in a multiplying nuclear system with source. Ann. Nucl. Energy 2020, 135, 106956. [Google Scholar] [CrossRef]
- Cacuci, D.G. The first-order comprehensive sensitivity analysis methodology (1st-CASAM) for scalar-valued responses: I. Theory. Am. J.Comput. Math. (AJCM) 2020, 10, 275–289. [Google Scholar] [CrossRef]
- Cacuci, D.G. The first-order comprehensive sensitivity analysis methodology (1st-CASAM) for scalar-valued responses: II. Illustrative application to a heat transport benchmark model. Am. J. Comput. Math. (AJCM) 2020, 10, 290–310. [Google Scholar] [CrossRef]
- Cacuci, D.G. Second-order adjoint sensitivity analysis methodology for computing exactly response sensitivities to uncertain parameters and boundaries of linear systems: Mathematical framework. Am. J. Comput. Math. (AJCM) 2020, 10, 329–354. [Google Scholar] [CrossRef]
- Cacuci, D.G. Illustrative application of the 2nd-order adjoint sensitivity analysis methodology to a paradigm linear evolution/transmission model: Point-detector response. Am. J. Comput. Math. (AJCM) 2020, 10, 355–381. [Google Scholar] [CrossRef]
- Cacuci, D.G. Illustrative Application of the 2nd-Order Adjoint Sensitivity Analysis Methodology to a Paradigm Linear Evolution/Transmission Model: Reaction-Rate Detector Response. Am. J. Comput. Math. (AJCM) 2020, 10, 382–397. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Cacuci, D.G. Fourth-Order Comprehensive Adjoint Sensitivity Analysis (4th-CASAM) of Response-Coupled Linear Forward/Adjoint Systems: I. Theoretical Framework. Energies 2021, 14, 3335. https://doi.org/10.3390/en14113335
Cacuci DG. Fourth-Order Comprehensive Adjoint Sensitivity Analysis (4th-CASAM) of Response-Coupled Linear Forward/Adjoint Systems: I. Theoretical Framework. Energies. 2021; 14(11):3335. https://doi.org/10.3390/en14113335
Chicago/Turabian StyleCacuci, Dan Gabriel. 2021. "Fourth-Order Comprehensive Adjoint Sensitivity Analysis (4th-CASAM) of Response-Coupled Linear Forward/Adjoint Systems: I. Theoretical Framework" Energies 14, no. 11: 3335. https://doi.org/10.3390/en14113335
APA StyleCacuci, D. G. (2021). Fourth-Order Comprehensive Adjoint Sensitivity Analysis (4th-CASAM) of Response-Coupled Linear Forward/Adjoint Systems: I. Theoretical Framework. Energies, 14(11), 3335. https://doi.org/10.3390/en14113335