Second-Order Adjoint Sensitivity Analysis Methodology for Computing Exactly Response Sensitivities to Uncertain Parameters and Boundaries of Linear Systems: Mathematical Framework
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)” for the efficient and exact computation of 1st- and 2nd-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (i.e., model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2nd-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1st-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2nd-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (N2/2 + 3 N/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2nd-LASS requires very little additional effort beyond the construction of the 1st-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1st-LASS and 2nd-LASS requires the same computational solvers as needed for solving (i.e., “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1st-LASS and the 2nd-LASS. Since neither the 1st-LASS nor the 2nd-LASS involves any differentials of the operators underlying the original system, the 1st-LASS is designated as a “first-level” (as opposed to a “first-order”) adjoint sensitivity system, while the 2nd-LASS is designated as a “second-level” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2nd-LASS that involve the imprecisely known boundary parameters. Notably, the 2nd-LASS encompasses an automatic, inherent, and independent “solution verification”
References
[1]
Wigner, E.P. (1945) Effect of Small Perturbations on Pile Period. Chicago Report CP-G-3048.
[2]
Dan Cacuci, G. (1981) Sensitivity Theory for Nonlinear Systems: I. Nonlinear Functional Analysis Approach. Journal of Mathematical Physics, 22, 2794-2802.
https://doi.org/10.1063/1.525186
[3]
Dan Cacuci, G. (1981) Sensitivity Theory for Nonlinear Systems: II. Extensions to Additional Classes of Responses. Journal of Mathematical Physics, 22, 2803-2812.
https://doi.org/10.1063/1.524870
[4]
Dan Cacuci, G. (2003) Sensitivity and Uncertainty Analysis: Theory, Volume 1. Chapman & Hall/CRC, Boca Raton.
[5]
Cacuci, D.G., Ionescu-Bujor, M. and Navon, M.I. (2005) Sensitivity and Uncertainty Analysis: Applications to Large Scale Systems, Volume 2. Chapman & Hall/CRC, Boca Raton. https://doi.org/10.1201/9780203483572
[6]
Dan Cacuci, G. (2015) 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. Journal of Computational Physics, 284, 687-699. https://doi.org/10.1016/j.jcp.2014.12.042
[7]
Dan Cacuci, G. (2016) Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) for Large-Scale Nonlinear Systems: I. Theory. Nuclear Science and Engineering, 184, 16-30. https://doi.org/10.13182/NSE16-16
[8]
Dan Cacuci, G. (2018) The Second-Order Adjoint Sensitivity Analysis Methodology. CRC Press, Taylor & Francis Group, Boca Raton.
https://doi.org/10.1201/9781315120270
[9]
Komata, M. (1977) A Generalized Perturbation Theory Applicable to Reactor Boundary Changes. Nuclear Science and Engineering, 64, 811.
https://doi.org/10.13182/NSE77-A14496
[10]
Larsen, E.W. and Pomraning, G.C. (1981) Boundary Perturbation Theory. Nuclear Science and Engineering, 77, 415. https://doi.org/10.13182/NSE81-A18954
[11]
Rahnema, F. and Pomraning, G.C. (1983) Boundary Perturbation Theory for Inhomogeneous Transport Equations. Nuclear Science and Engineering, 84, 313-319.
https://doi.org/10.13182/NSE83-A15451
[12]
Rahnema, F. (1996) Boundary Condition Perturbations in Transport Theory. Nuclear Science and Engineering, 124, 320. https://doi.org/10.13182/NSE96-A28581
[13]
McKinley, M.S. and Rahnema, F. (2002) High-Order Boundary Condition Perturbation Theory for the Neutron Transport Equation. Nuclear Science and Engineering, 140, 285-294. https://doi.org/10.13182/NSE02-A2261
[14]
Favorite, J.A. (2006) Variational Estimates of Internal Interface Perturbations and a New Variational Functional for Inhomogeneous Transport Problems. Nuclear Science and Engineering, 152, 180. https://doi.org/10.13182/NSE06-A2574
[15]
Favorite, J.A. and Gonzalez, E. (2017) Revisiting Boundary Perturbation Theory for Inhomogeneous Transport Problems. Nuclear Science and Engineering, 185, 445-459.
https://doi.org/10.1080/00295639.2016.1277108
[16]
Dan Cacuci, G. (2019) Second-Order Sensitivities of a General Functional of the Forward and Adjoint Fluxes in a Multiplying Nuclear System with Source. Nuclear Engineering and Design, 344, 83-106.
https://doi.org/10.1016/j.nucengdes.2019.01.007
[17]
Dan Cacuci, G. (2019) Second Order Adjoint Sensitivity Analysis of a General Ratio of Functionals of the Forward and Adjoint Fluxes in a Multiplying Nuclear System with Source. Annals of Nuclear Energy, 135, Article ID: 106956.
https://doi.org/10.1016/j.anucene.2019.106956
[18]
Cacuci, D.G. (2020) Illustrative Application of the 2nd-Order Adjoint Sensitivity Analysis Methodology to a Paradigm Linear Evolution/Transmission Model: Point-Detector Response. AJCM, 10, 355-381.
[19]
Cacuci, D.G. (2020) Illustrative Application of the 2nd-Order Adjoint Sensitivity Analysis Methodology to a Paradigm Linear Evolution/Transmission Model: Reaction-Rate Detector Response. AJCM, 10, 382-397.
