TY - JOUR

T1 - An improved computational method for sensitivity analysis

T2 - Green's function method with 'AIM'

AU - Kramer, M. A.

AU - Calo, J. M.

AU - Rabitz, H.

N1 - Funding Information: The authors acknowledge the Office of Naval Research, the Department of Energy, and the National Science Foundation for support of this work.

PY - 1981/12

Y1 - 1981/12

N2 - Sensitivity analysis of initial value models via the calculation of linear sensitivity coefficients is quite important for model evaluation and validation. Direct solution of the sensitivity equations for n-dimensional, m-parameter systems of ordinary differential equations requires the solution of m × n differential equations, which can become quite expensive for large-scale models. When m > n (the usual case for chemical kinetic systems, for example), the Green's function method (GFM), which requires solutions of n2 differential equations with m × n subsequent numerical quadratures, is the most efficient computational technique for determining linear sensitivity coefficients. Even so, associated computing costs can still become quite large. In the current work, an algorithm, known as the analytical integrated Magnus (AIM) modification of the GFM, is presented which dramatically reduces the computational effort required to determine linear sensitivity coefficients. The technique employs the piecewise Magnus method for more efficient calculation of Green's function kernels, and treats the sensitivity integrals analytically. An application of this technique to a chemical kinetics system is presented in which the computational effort is reduced by an order of magnitude in comparison to the unmodified GFM.

AB - Sensitivity analysis of initial value models via the calculation of linear sensitivity coefficients is quite important for model evaluation and validation. Direct solution of the sensitivity equations for n-dimensional, m-parameter systems of ordinary differential equations requires the solution of m × n differential equations, which can become quite expensive for large-scale models. When m > n (the usual case for chemical kinetic systems, for example), the Green's function method (GFM), which requires solutions of n2 differential equations with m × n subsequent numerical quadratures, is the most efficient computational technique for determining linear sensitivity coefficients. Even so, associated computing costs can still become quite large. In the current work, an algorithm, known as the analytical integrated Magnus (AIM) modification of the GFM, is presented which dramatically reduces the computational effort required to determine linear sensitivity coefficients. The technique employs the piecewise Magnus method for more efficient calculation of Green's function kernels, and treats the sensitivity integrals analytically. An application of this technique to a chemical kinetics system is presented in which the computational effort is reduced by an order of magnitude in comparison to the unmodified GFM.

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

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

U2 - 10.1016/S0307-904X(81)80027-3

DO - 10.1016/S0307-904X(81)80027-3

M3 - Article

SN - 0307-904X

VL - 5

SP - 432

EP - 441

JO - Applied Mathematical Modelling

JF - Applied Mathematical Modelling

IS - 6

ER -