TY - JOUR

T1 - A general analysis of exact nonlinear lumping in chemical kinetics

AU - Li, Genyuan

AU - Rabitz, Herschel

AU - Tóth, János

N1 - Funding Information:
Acknowledgments-The authors acknowledge support from the Denartment of Enerev and Exxon Cornoration. One of the authors (J.T.) has Go been support& by the Soros Foundation and by the National Scientific Research Fund, Hungary, No. 328s.

PY - 1994

Y1 - 1994

N2 - A general analysis of exact nonlinear lumping is presented. This analysis can be applied to the kinetics of any reaction system with n species described by a set of first-order ordinary differential equations dy/dt = f(y), where y is an n-dimensional vector and f(y) is an arbitrary n-dimensional function vector. We consider lumping by means of n̂ (n̂ ≤ n)-dimensional arbitrary transformation ŷ = h(y). The lumped differential equation system is d y ̂ Dt = y(h(ŷ))f(h(ŷ)), where hy(y) is teh Jacobian matrix of h(y), h is a generalized inverse transformation of h satisfying the relation h(h) = In̂. Three necessary and sufficient conditions of the existence of exact nonlinear lumping schemes have been determined. The geometric and algebraic interpretations of these conditions are discussed. It is found that a system is exactly lumpable by h only if h(y) = 0 is its invariant manifold. A linear partial differential operator A = Σni=1 fi(y)θ{symbol}/θ{symbol}yi corresponding to dy dt = f(y) is defined. Using the eigenfunctions and the generalized eigenfunctions of A, the operator can be transformed to Jordan or diagonal canonical forms which give the lumped differential equation systems without determination of h. These approaches are illustrated by a simple example. The results of this analysis serve as a theoretical basis for the development of approaches for approximate nonlinear lumping.

AB - A general analysis of exact nonlinear lumping is presented. This analysis can be applied to the kinetics of any reaction system with n species described by a set of first-order ordinary differential equations dy/dt = f(y), where y is an n-dimensional vector and f(y) is an arbitrary n-dimensional function vector. We consider lumping by means of n̂ (n̂ ≤ n)-dimensional arbitrary transformation ŷ = h(y). The lumped differential equation system is d y ̂ Dt = y(h(ŷ))f(h(ŷ)), where hy(y) is teh Jacobian matrix of h(y), h is a generalized inverse transformation of h satisfying the relation h(h) = In̂. Three necessary and sufficient conditions of the existence of exact nonlinear lumping schemes have been determined. The geometric and algebraic interpretations of these conditions are discussed. It is found that a system is exactly lumpable by h only if h(y) = 0 is its invariant manifold. A linear partial differential operator A = Σni=1 fi(y)θ{symbol}/θ{symbol}yi corresponding to dy dt = f(y) is defined. Using the eigenfunctions and the generalized eigenfunctions of A, the operator can be transformed to Jordan or diagonal canonical forms which give the lumped differential equation systems without determination of h. These approaches are illustrated by a simple example. The results of this analysis serve as a theoretical basis for the development of approaches for approximate nonlinear lumping.

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

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

U2 - 10.1016/0009-2509(94)87006-3

DO - 10.1016/0009-2509(94)87006-3

M3 - Article

AN - SCOPUS:0028378405

VL - 49

SP - 343

EP - 361

JO - Chemical Engineering Science

JF - Chemical Engineering Science

SN - 0009-2509

IS - 3

ER -