Sensitivity analysis is a powerful first-order mathematical technique that identifies the major reactions, the effects of errors in rate constants, pressure effects, and temperature effects in a plasma chemistry type model. Several methods of sensitivity analysis have been discussed in the literature. This paper describes the successful implementation of the state transition matrix method to the sensitivity analysis of a Xe| plasma chemistry model using both e-beams and ions for excitation. The ion excited plasma study is important to the field of nuclear-pumped lasers. This model examines 8 species and 37 reactions.
IntroductionDynamic modeling of plasma kinetics is very crucial to the understanding of system performance. However, due to the large number of possible reactions in typical plasma chemistry models, and the uncertainties associated with the rate constants for the reactions, the brute force approach of analyzing each reaction separately to seek out important reactions over different regimes of operation (e.g. pressure, and temperature) can be computationally tedious and prohibitive.Methods of first-order sensitivity analysis can prove to be a powerful mathematical tool which identifies the major reactions, pressure effects, temperature dependence and rate constant uncertainties. Several methods have been examined: the Fourier Amplitude Sensitivity Test (FAST) (Cukier 1973;; the Direct Method (DM) (Atherton 1975); the Green's Function Method (GFM) (Hwang 1978; Dougherty 1979;Demiralp 1981); and the State Transition Matrix Method (STMM). The FAST and the DM focus upon the evaluation of sensitivity coefficients (i.e. the partial derivatives of the system solutions with respect to the system's parameters). In the FAST approach, each parameter is varied as a periodic function of a "search variable", s. Hence the system solutions are periodic in s and can be Fourier analysed. Despite the success of the FAST approach, it still has several weaknesses: it is very expensive to use; and it yields time averaged sensitivity coefficients (impossible to explain transient effects). The DM solves the original set of system equations together with a large set of sensitivity equations. This approach has an advantage in that the sensitivity coefficients are local and therefore useful in transient analysis. However, the stiffness ratio of the system is very large (i.e. a total of mn + n equations must be solved simultaneously where m = number of parameters and n = number of system equations) (Atherton 1973).The Green's Function Method transforms the set of sensitivity equations into integral equations. For well-defined physical systems, the integrands are smooth. This reduces not only the stiffness ratio but also the computational cost. The STMM approach follows the GFM except that the state transition matrix is solvable by methods developed by linear systems engineering. The latter method was chosen to