This paper discusses an accurate numerical approach based on complex variables for the computation of the Jacobian matrix of complex chemical reaction mechanisms. The Jacobian matrix is required in the calculation of low dimensional manifolds during kinetic chemical mechanism reduction. The approach is suitable for numerical computations of large-scale problems and is more accurate than the finite difference approach of computing Jacobians. The method is demonstrated via a nonlinear reaction mechanism for the synthesis of Bromide acid and a H2/Air mechanism using a modified CHEMKIN package. The Bromide mechanism consisted of five species participating in six elementary chemical reactions and the H2/Air mechanism consisted of 11 species and 23 reactions. In both cases it is shown that the method is superior to the finite difference approach of computing derivatives with an arbitrary computational step size h.

Issue Section:
Technical Briefs
Keywords:
chemical reactions, combustion, simulation, air, hydrogen
Topics:
Chemical reactions, Combustion, Computation, Errors, Hydrogen, Jacobian matrices, Manifolds, Simulation, Differential equations, Computer software
1.
Maas
,
U.
, and
Pope
,
S. B.
,
1992
, “
Simplifying Chemical Kinetics: Intrinsic Low-Dimensional Manifolds in Composition Space
,”
Combust. Flame
,
88
, pp.
239
–
264
.
2.
Butuk, N. K., and Pemba, J.-P., 2001, “Computing Low Dimensional Invariant Manifolds of Reaction Mechanism Using Genetic Algorithms,” AIAA Paper 2001-1076.
3.
Tomlin, A. S., Turanyi, T., and Pilling, M. J., 1997, “Mathematical Tools for the Construction, Investigation and Reduction of Combustion Mechanisms,” Low Temperature Combustion and Autoignition, M. J. Pilling and G. Hancock, eds., Elsevier, New York, pp. 293–437.
4.
Rabitz
,
H.
,
Kramer
,
M.
, and
Dacol
,
D.
,
1983
, “
Sensitivity Analysis in Chemical Kinetics
,”
Annu. Rev. Phys. Chem.
,
34
, pp.
419
–
461
.
5.
Hwang
,
J.-T.
,
Dougherty
,
E. P.
,
Rabitz
,
S.
, and
Rabitz
,
H.
,
1987
, “
The Green’s Function Method of Sensitivity Analysis in Chemical Kinetics
,”
J. Chem. Phys.
,
69
(
11
), pp.
5180
–
5191
.
6.
Carmichael
,
G. R.
,
Sandu
,
A.
, and
Potra
,
F. A.
,
1997
, “
Sensitivity Analysis for Atmospheric Chemistry Models via Automatic Differentiation
,”
Atmos. Environ.
,
31
(
3
), pp.
475
–
489
.
7.
Shen
,
J.
,
1999
, “
A Direct Method of Calculation Sensitivity Coefficients of Chemical Kinetics
,”
J. Chem. Phys.
,
111
(
6
), pp.
7209
–
7214
.
8.
Newman, III, J. C., Anderson, K. W., and Whitfield, D. L., 1998, “Multi-disciplinary Sensitivity Derivatives Using Complex Variables,” Computational Fluid Dynamics Laboratory, Mississippi State University, MS, MSSU-COE-ERC-98-08.
9.
Squire
,
W.
, and
Trapp
,
G.
,
1998
, “
Using Complex Variables to Estimate Derivatives of Real Functions
,”
SIAM Rev.
,
10
(
1
), pp.
110
–
112
.
10.
Anderson, K., Nielson, E., Whitefield, D., and Newman, J., 1999, “Sensitivity Analysis for the Navier-Stokes Equations on Unstructured Meshes Using Complex Variables,” AIAA Paper 99-3294.
11.
Ye, F., Butuk, N. K., Huque, Z., and Biney, P., 1999, “CFD Optimization of an Integrated Multi-Tube Pulse Detonation Engine Inlet Concept Using Complex Variables,” AIAA Paper 99-3229.
12.
Pemba, J.-P., 2001, “The Complex Method in Numerical Differentiation,” preprint, submitted for publication.
13.
Martins, J. R. R. A., Sturdza, P., and Alonso, J. J., 2001, “The Connection Between the Complex-Step Derivative Approximation and Algorithmic Differentiation,” AIAA Paper 2001-0921.
14.
Kee, R. J., Rupley, F. M., Meeks, E., and Miller, J. A., 1996, Chemkin III User Manual, Reaction Design, San Diego, CA
15.
Deco, G., and Obradovic, D., 1996, An Information-Theoretic Approach to Neural Computing, Springer-Verlag, New York, pp. 152–156.
Copyright © 2003
 by ASME
You do not currently have access to this content.