3.2.1. Second-Order Sensitivities Corresponding to
The application of the 4th-CASAM-N to compute efficiently and exactly second-order sensitivities will be illustrated in this subsection by considering the 1st-order sensitivity , the expression of which has been obtained in Equation (82). The sensitivity is representative of the procedure involved when applying the 4th-CASAM-N to a 1st-order sensitivity that involves a single state function, which in this case is the 1st-level adjoint function .
The second-order sensitivities corresponding to
are obtained by determining the G-differential of Equation (83), which has, by definition, the following expression:
The direct-effect term
is defined to depend only on parameter variations and, in the case of Equation (98), stems only from the derivative with respect to
of the upper limit of integration, i.e.,
As indicated in Equation (99), the direct-effect term vanishes because of the boundary condition satisfied by the 1st-level adjoint function
provided in Equation (79). Therefore, the G-differential
comprises only the indirect-effect term
, which depends only on the variation
and is defined as follows:
The indirect-effect term defined by Equation (100) can be evaluated only after having determined the variational function
, which is the solution of the system obtained by taking the G-differential of the 1st-LASS defined by Equations (78)–(80), and which has the following expression:
The system comprising Equations (101)–(103) is actually the 2nd-level variational sensitivity system (2nd-LVSS) for the function
. Its solution,
, could be used to determine the G-differential
shown in Equation (98). As
depends on parameter variations, however, solving repeatedly the 2nd-LVSS for all parameter variations is avoided by constructing a corresponding 2nd-LASS, which would need to be solved only once, as its solution would be independent of parameter variations and would be used to obtain the G-differential
. As the 2nd-LVSS comprising Equations (101)–(103) has the same structure as the 1st-LVSS comprising Equations (26)–(28), it follows that the 2nd-LASS that corresponds to the 2nd-LVSS comprising Equations (101)–(103) will be constructed by applying the same principles that were applied to construct the 1st-LASS in
Section 3.1 when determining the 1st-order sensitivities of the response
. Thus, the inner product defined in Equation (30) will be used to construct the inner product of Equation (101) with a yet undefined one-component function
—which will ultimately become the sought-after 2nd-level adjoint function—to obtain the following relation:
The superscript “2” in the notation
indicates “2nd-level” while the subscript “
Q” indicates that that this 2nd-level adjoint function corresponds to the parameter “Q”. Integrating the left side of Equation (104) twice by parts yields the following relation:
Using the relation provided in Equation (104) and inserting into Equation (105) the boundary conditions provided in Equations (102) and (103) makes it possible to recast Equation (105) into the following form:
The left side of Equation (106) is now required to represent the G-differential defined in Equation (100). Furthermore, the definition of the 2nd-level adjoint function
is completed by requiring that it satisfy boundary conditions, which would eliminate the unknown boundary terms from Equation (106). Imposing these requirements yields the following 2nd-LASS for
:
Recalling Equation (100) and implementing the relations represented by the 2nd-LASS into Equation (106) transforms the latter relation into the following form:
Identifying the quantities that multiply the various parameter variations in Equation (110) yields the following expressions for the corresponding 2nd-order sensitivities:
The expressions on the right sides of Equations (112)–(115) are to be evaluated at the nominal values of the respective parameters and state functions but the notation has been omitted for simplicity.
The expressions for the 2nd-order sensitivities obtained in Equations (112)–(115) can be evaluated inexpensively, using quadrature formulas, after having solved the 2nd-LASS once to obtain the 2nd-level adjoint function. Thus, solving the 2nd-LASS is the sole “large-scale” computation needed in order to compute the respective sensitivities. In contradistinction, using forward methods would have required at least as many “large-scale” computations as there are model parameters.
Solving Equations (107)–(109) yields the following explicit expression for the 2nd-level adjoint function
:
Using the expressions for
and
in Equations (112)–(115) yields the following closed-form expressions for the respective 2nd-order sensitivities:
Of course, the above closed-form expressions are for verification purposes; in practice, the values of the respective sensitivities are computed numerically using Equations (112)–(115), as has been mentioned in the foregoing.
3.2.2. Second-Order Sensitivities Corresponding to
The application of the 4th-CASAM-N to compute efficiently and exactly second-order sensitivities will be illustrated in this subsection by considering the 1st-order sensitivity , the expression of which has been obtained in Equation (83). The sensitivity is representative of the procedure involved when applying the 4th-CASAM-N to a 1st-order sensitivity that involves both the original forward function and the 1st-level adjoint function .
The second-order sensitivities corresponding to
are obtained by determining the G-differential of Equation (83), which has, by definition, the following expression:
where the direct-effect term
depends only on parameter variations and is defined as follows:
and where the indirect-effect term
depends only on variations
and
, and is defined as follows:
The direct-effect term defined by Equation (122) can be evaluated immediately, as the state functions
and
are available already; in practice, this direct-effect term is computed numerically. However, the indirect-effect term defined by Equation (123) can be evaluated only after having determined the functions
and
. The function
is the solution of the 1st-LVSS represented by Equations (26)–(28). The function
is the solution of the system obtained by taking the G-differential of the 1st-LASS defined by Equations (78)–(80), which has the following expression:
Concatenating the 1st-LVSS and the G-differentiated 1st-LASS yields the following 2nd-level variational sensitivity system (2nd-LVSS), comprising Equations (127) and (128) below, for obtaining the 2nd-level variational function
, which is also defined below:
where
with
To distinguish vectors from matrices, two capital bold letters were used (and will henceforth be used) to denote matrices, as in the case of “the second-level variational matrix” . The “2nd-level” is indicated by the superscript “(2)”. Subsequently in this work, levels higher than the second will also be indicated by a corresponding superscript attached to the appropriate (block-) vectors and/or (block-) matrices. The argument “”, which appears in the list of arguments of , indicates that this matrix is a -dimensional matrix.
The need for solving the 2nd-LVSS is circumvented by deriving an alternative expression for the indirect-effect term defined in Equation (123), in which the function
is replaced by a 2nd-level adjoint function that is independent of variations in the model parameter and state functions. This 2nd-level adjoint function will satisfy a 2nd-level adjoint sensitivity system (2nd-LASS) which will be constructed by using the 2nd-LVSS as the starting point and following the same principles as outlined in
Section 3.1. The 2nd-LASS will be constructed in a Hilbert space, which will be denoted as
and which is comprised of element vectors of the same form as
; a generic vector in
is a two-component column vector of the form
. The inner product of two vectors
and
in the Hilbert space
will be denoted as
and defined as follows:
The inner product defined in Equation (133) is continuous in
in a neighborhood of
. Using the definition of the inner product defined in Equation (133), construct the inner product of Equation (127) with a vector
to obtain the following relation:
The notation chosen for the vector indicates the following characteristics of this vector: (i) the bold letter indicates that this quantity is a vector; (ii) the superscript “2” indicates that this quantity is a “second-level” function; (iii) the subscript “k” indicates that this quantity corresponds to the 1st-order sensitivity of the response with respect to the parameter ; (iv) the digit “2” in the argument of indicates that this vector has two components; (v) the letter “x” in the argument of indicates that the two components of this vector are functions of the independent variable “x”.
In component form, Equation (134) has the following expression:
The inner product on the left side of Equation (134) is now further transformed by using the definition of the adjoint operator to obtain, after integrating twice by parts over the independent variable
, the following relation:
In Equation (136), the adjoint (matrix-valued) operator
has the following form:
while the quantity
denotes the corresponding bilinear concomitant on the domain’s boundary, evaluated at the nominal values for the parameters and respective state functions, and having the following expression:
Inserting into Equation (138) the 2nd-LVSS’s boundary conditions provided in Equation (128) reduces the bilinear concomitant to the following expression:
The first term on the right side of Equation (136) is now required to represent the indirect-effect term defined in Equation (123), which is accomplished by imposing the following relation:
where
The definition domain of the adjoint (matrix-valued) operator
is now specified by requiring the function
to satisfy adjoint boundary/initial conditions which fulfill the following conditions: (a) they must be independent of unknown values of
; (b) they must cause all terms containing unknown values of
in the expression of
to vanish. The above requirements are satisfied by imposing the following boundary conditions:
where:
The system of equations comprising Equations (140) and (142) constitutes the 2nd-LASS for the 2nd-level adjoint function .
Using the equations underlying the 2nd-LASS together with Equations (134) and (139) in Equation (136) yields the following expression for the indirect-effect term defined in Equation (123):
Replacing in Equation (147) the source terms and with their expressions provided in Equation (130) and adding the resulting equation to the expression for the direct-effect term provided in Equation (122) yields the following expression for the G-differential, defined in Equation (121):
Identifying in Equation (148) the quantities that multiply the various parameter variations yields the following expressions for the respective partial second-order sensitivities:
The expressions presented in Equations (149)–(153) are to be evaluated at the nominal values of the model parameters and state functions but the indicator
has been omitted for simplicity. As the expressions and values for the functions
and
are already available, the expressions provided in Equations (149)–(153) can be evaluated by using quadrature formulas after having determined the 2nd-level adjoint function
by solving the 2nd-LASS. Solving Equation (140) subject to the boundary conditions provided in Equation (142) yields the following expressions:
Inserting the expressions for
,
and
in Equations (149)–(153) and performing the operations indicated on the right sides of these equations yields the following closed-form expressions for the respective second-order sensitivities:
It is important to note that the symmetry relationship
provides the following equality between the expressions provided in Equations (112) and (149):
The identity provided in Equation (161) enables an intrinsic mutual verification of the accuracy of the numerical computation of the respective adjoint functions.
As this paradigm illustrative heat conduction model comprises five uncertain model parameters, there will be a total of 25 second-order sensitivities of the response
with respect to the model parameters. Of these 25 second-order sensitivities, 15 will be distinct, while the others are not because of the intrinsic symmetries of the 2nd-order mixed derivatives. Of the 25 second-order sensitivities, 10 have been obtained in
Section 3.2. As the computation of the remaining 15 second-order sensitivities of the response
with respect to the model parameters, are to be performed by applying the same 4th-CASAM-N principles, as illustrated thus far in
Section 3.2, the respective derivations will be omitted for brevity. It is important to note that if all of these computations are performed, the mixed 2nd-order sensitivities of
will have been computed twice, due to their symmetries, in the same manner as indicated in Equation (161). These symmetry features enable the intrinsic mutual verification of the accuracy of the numerical computation of all of the 1st-level and 2nd-level adjoint functions involved in the respective computations.