Towards Overcoming the Curse of Dimensionality: The Third-Order Adjoint Method for Sensitivity Analysis of Response-Coupled Linear Forward/Adjoint Systems, with Applications to Uncertainty Quantification and Predictive Modeling
Abstract
:1. Introduction
2. Mathematical Description of the Physical System
- denotes a -dimensional column vector whose components are the physical system’s imprecisely known parameters, which are subject to uncertainties; , where denotes a subset of the -dimensional real vector space . The symbol “” will be used to denote “is defined as” or “is by definition equal to.” The vector is considered to include any imprecisely known model parameters that may enter into defining the system’s boundary in the phase space of independent variables.
- denotes the -dimensional phase-space position vector, defined on a phase-space domain denoted as .
- denotes a -dimensional column vector whose components represent the system’s dependent variables (also called “state functions”). In virtually all of the physical systems represented by Equation (1), the components are square-integrable functions and , where is a Hilbert space endowed with an inner product that will be denoted as , , and which is defined as follows:
- denotes a -dimensional column vector. The components of are operators acting (in general nonlinearly) on and ;
- denotes a -component column vector. The components of are operators acting linearly on and nonlinearly on . When contains differential operators, a set of boundary and/or initial conditions which define the domain of must also be given. Since is considered to act linearly on , the accompanying boundary and/or initial conditions must also be linear in . Such linear boundary and/or initial conditions are represented in the following operator form:
3. The Third-Order Adjoint Sensitivity Analysis Methodology (3rd-ASAM) for Coupled Linear Forward and Adjoint Systems: Another Step Towards Overcoming the Curse of Dimensionality in the Exact Computation of High-Order Response Sensitivities
3.1. The First-Level Adjoint System (1st-LASS) for Computing Exactly and Efficiently the First-Order Model Response Sensitivities to Parameters
- (i)
- Consider two vector-valued functions and , each having two -dimensional vector-components defined as follows: and . The components of these vectors are assumed to be square-integrable functions.
- (ii)
- Introduce a Hilbert space, denoted as , endowed with the following inner product, denoted as , , , between the two functions defined in item (i), above:
- (iii)
- In the Hilbert space , form the inner product of Equation (16) with a yet undefined vector-valued function to obtain the following relation, evaluated at , in which the superscript “zero” is omitted to simplify the notation:
- (iv)
- Use the definition of the adjoint operator in the Hilbert space to recast the left-side of Equation (20) as follows:
- (v)
- Identify the term on the left-side of Equation (21) with the indirect effect term defined in Equation (15), i.e., require that
- (vi)
- The boundary conditions given in Equation (17) are now implemented in Equation (21), thereby reducing by half the number of unknown boundary-values in the bilinear concomitant . The boundary conditions for the adjoint functions and are chosen next so as to eliminate the remaining unknown boundary-values of the functions and while ensuring that Equation (22) is well posed. The boundary conditions thus chosen for the adjoint functions and can be represented in operator form as follows:
- (vii)
- In most cases, the above choice of boundary conditions for the 1st-level adjoint function will cause the bilinear concomitant in Equation (21) to vanish. When the boundary conditions for the original system are non-homogeneous, however, the bilinear concomitant may not vanish. Even when it does not vanish, however, this bilinear concomitant will be reduced to a quantity, denoted here as , which will contain only known values of its arguments.
- (viii)
- Use the 1st-LASS defined by Equations (22) and (23) together with Equations (20) and (21) to obtain the following expression for the indirect-effect term defined in Equation (15), in terms of the adjoint functions and :
3.2. The Second-Level Adjoint System (2nd-LASS) for Computing Exactly and Efficiently the Second-Order Model Response Sensitivities to Parameters
- (ii)
- Use the definition of the adjoint operator in the Hilbert space to recast the left-side of Equation (33) as follows:
- (iii)
- Identify the first term on the right-side of Equation (34) with the indirect-effect term defined in Equation (27) by requiring that the following system of equations be satisfied for :
- (iv)
- The boundary conditions given in Equations (17) and (29) are now implemented in Equation (34), thereby reducing by half the number of unknown boundary-values in the bilinear concomitant . The boundary conditions for the 2nd-level adjoint functions are now chosen so as to eliminate the remaining unknown boundary-values of the functions , , and while ensuring that Equation (35) is well posed. The boundary conditions thus chosen for the adjoint functions and can be represented in operator form as follows:
3.3. The Third-Level Adjoint System (3rd-LASS) for Computing Exactly and Efficiently the Third-Order Model Response Sensitivities to Parameters
- (ii)
- Use the definition of the adjoint operator in the Hilbert space to recast the left-side of Equation (64) as follows:
- (iii)
- Identify the first term on the right-side of Equation (65) with the indirect-effect term defined in Equation (41) by requiring that:
- (iv)
- The boundary conditions given for the 2nd-LFSS and those given in Equation (62) are now implemented in Equation (65), thereby reducing by half the number of unknown boundary-values in the bilinear concomitant . The boundary conditions for the 3rd-level adjoint functions are chosen next so as to eliminate the remaining unknown boundary-values of the functions while ensuring that Equation (66) is well posed. The boundary conditions thus chosen for the adjoint functions can be represented in operator form as follows:
- (vi)
- Replace Equation (68) in Equation (40) to obtain the following expression for the total 2nd-order response sensitivity to model parameters:
- (vii)
- Note that the 2nd-LASS is independent of parameter variations . Thus, the exact computation of all of the partial third-order sensitivities, , requires at most large-scale (adjoint) computations using the 3rd-LASS, rather than large-scale computations as would be required by forward methods. In order to implement the practical computation of the 3rd-level adjoint functions, it is important to note that, in component form, Equation (66) has the following structure, for each :
4. Third-Order Expressions for the Cumulants of the Response Distribution in Parameter Space
5. 2nd/3rd-Order Best-Estimated Results with Reduced Uncertainties Predictive Modeling (2nd/3rd-BERRU-PM) in the Joint Phase-Space of Responses and Parameters
5.1. 2nd/3rd-BERRU-PM: A Priori Information
5.1.1. Expected Values and Covariances of Measured Responses
5.1.2. Expectations and Covariances of Computed Responses Including Second- and Third-Order Sensitivities to Model Parameters
5.2. 2nd/3rd-BERRU-PM: Analytical Expressions for Best-Estimate Results with Reduced Uncertainties for Responses and Parameters in the Joint Phase-Space of Responses and Parameters
5.2.1. Predicted Best-Estimate Expected Values for the Responses and Parameters in the Joint Phase-Space of Responses and Parameters
5.2.2. Predicted Best-Estimate Covariances for the Responses and Parameters in the Joint Phase-Space of Responses and Parameters
5.2.3. Data Consistency Indicator
6. Conclusions
Funding
Conflicts of Interest
Appendix A
- (i)
- The vector of mean values (first-order moments), denoted as , of , and defined as ;
- (ii)
- The second-order moment or covariance, , of two parameters, and , defined as . The covariances constitute the elements of a symmetric, positive-definite parameter covariance matrix of dimension , denoted as
References
- 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]
- 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, 4100. [Google Scholar] [CrossRef]
- 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]
- 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]
- Cacuci, D.G. BERRU Predictive Modeling: Best Estimate Results with Reduced Uncertainties; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2018. [Google Scholar]
- Lewis, J.M.; Lakshmivarahan, S.; Dhall, S.K. Dynamic Data Assimilation: A Least Square Approach; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
- Lahoz, W.; Khattatov, B.; Ménard, R. Data Assimilation: Making Sense of Observations; Springer: New York, NY, USA, 2010. [Google Scholar]
- Goldstein, M.; Woof, D. Bayes Linear Statistics: Theory and Methods; John Wiley & Sons: Chichester, UK, 2007. [Google Scholar]
- Jaynes, E.T. Information Theory and Statistical Mechanics. Phys. Rev. 1957, 106, 620–630. [Google Scholar] [CrossRef]
- Granziol, D.; Ru, B.; Zohren, S.; Dong, X.; Osborne, M.; Roberts, S. MEMe: An Accurate Maximum Entropy Method for Efficient Approximations in Large-Scale Machine Learning. Entropy 2019, 21, 551. [Google Scholar] [CrossRef]
- Rowlands, J.L.; Dean, C.J.; MacDougall, J.D.; Smith, R.W. The Production and Performance of the Adjusted Cross-Section Set FGL5. In Proceedings of the International Symposium “Physics of Fast Reactors”, Tokio, Japan, 16–23 October 1973. [Google Scholar]
- Gandini, A.; Petilli, M. AMARA: A Code Using the Lagrange Multipliers Method for Nuclear Data Adjustment; CNEN-RI/FI(73)39; Comitato Nazionale Energia Nucleare: Rome, Italy, 1973.
- Kuroi, H.; Mitani, H. Adjustment to Cross-Section Data to Fit Integral Experiments by Least Squares Method. J. Nucl. Sci. Technol. 1975, 12, 663. [Google Scholar] [CrossRef]
- Dragt, J.B.; Dekker, J.W.M.; Gruppelaar, H.; Janssen, A.J. Methods of Adjustment and Error Evaluation of Neutron Capture Cross Sections. Nucl. Sci. Eng. 1977, 62, 11. [Google Scholar] [CrossRef]
- Weisbin, C.R.; Oblow, E.M.; Marable, J.H.; Peelle, R.W.; Lucius, J.L. Application of Sensitivity and Uncertainty Methodology to Fast Reactor Integral Experiment Analysis. Nucl. Sci. Eng. 1978, 66, 307. [Google Scholar] [CrossRef]
- Barhen, J.; Cacuci, D.G.; Wagschal, J.J.; Bjerke, M.A.; Mullins, C.B. Uncertainty Analysis of Time-Dependent Nonlinear Systems: Theory and Application to Transient Thermal Hydraulics. Nucl. Sci. Eng. 1982, 81, 23–44. [Google Scholar] [CrossRef]
- Cacuci, D.G.; Ionescu-Bujor, M. Sensitivity and Uncertainty Analysis, Data Assimilation and Predictive Best-Estimate Model Calibration. In Handbook of Nuclear Engineering; Cacuci, D.G., Ed.; Springer: New York, NY, USA; Berlin/Heidelberg, Germnay, 2010; Volume 3, Chapter 17; pp. 1913–2051. ISBN 978-0-387-98150-5. [Google Scholar]
- Badea, M.C.; Cacuci, D.G.; Badea, A.F. Best-Estimate Predictions and Model Calibration for Reactor Thermal-Hydraulics. Nucl. Sci. Eng. 2012, 172, 1–19. [Google Scholar] [CrossRef]
- Arslan, E.; Cacuci, D.G. Predictive Modeling of Liquid-Sodium Thermal-Hydraulics Experiments and Computations. Ann. Nucl. Energy 2014, 63, 355–370. [Google Scholar] [CrossRef]
- Cacuci, D.G.; Arslan, E. Reducing Uncertainties via Predictive Modeling: FLICA4 Calibration Using BFBT Benchmarks. Nucl. Sci. Eng. 2014, 176, 339–349. [Google Scholar] [CrossRef]
- Peltz, J.J.; Cacuci, D.G.; Badea, A.F.; Badea, M.C. Predictive Modeling Applied to a Spent Fuel Dissolver Model: II. Uncertainty Quantification and Reduction. Nucl. Sci. Eng. 2016, 183, 332–346. [Google Scholar] [CrossRef]
- Badea, A.F.; Cacuci, D.G. Predictive Uncertainty Reduction in Coupled Neutron-Kinetics/Thermal Hydraulics Modeling of the BWR-TT2 Benchmark. Nucl. Eng. Des. 2017, 313, 330–344. [Google Scholar] [CrossRef]
- Peltz, J.J.; Cacuci, D.G. Inverse Predictive Modeling of a Spent Fuel Dissolver Model. Nucl. Sci. Eng. 2016, 184, 1–15. [Google Scholar] [CrossRef]
- 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]
© 2019 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 (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Cacuci, D.G. Towards Overcoming the Curse of Dimensionality: The Third-Order Adjoint Method for Sensitivity Analysis of Response-Coupled Linear Forward/Adjoint Systems, with Applications to Uncertainty Quantification and Predictive Modeling. Energies 2019, 12, 4216. https://doi.org/10.3390/en12214216
Cacuci DG. Towards Overcoming the Curse of Dimensionality: The Third-Order Adjoint Method for Sensitivity Analysis of Response-Coupled Linear Forward/Adjoint Systems, with Applications to Uncertainty Quantification and Predictive Modeling. Energies. 2019; 12(21):4216. https://doi.org/10.3390/en12214216
Chicago/Turabian StyleCacuci, Dan Gabriel. 2019. "Towards Overcoming the Curse of Dimensionality: The Third-Order Adjoint Method for Sensitivity Analysis of Response-Coupled Linear Forward/Adjoint Systems, with Applications to Uncertainty Quantification and Predictive Modeling" Energies 12, no. 21: 4216. https://doi.org/10.3390/en12214216
APA StyleCacuci, D. G. (2019). Towards Overcoming the Curse of Dimensionality: The Third-Order Adjoint Method for Sensitivity Analysis of Response-Coupled Linear Forward/Adjoint Systems, with Applications to Uncertainty Quantification and Predictive Modeling. Energies, 12(21), 4216. https://doi.org/10.3390/en12214216