of chemical kinetics from the standpoint of the method of invariant manifold (MIM). MIM is based on a formulation of the condition of invariance as an equation, ...
Available online at www.sciencedirect.com
Physica A 333 (2004) 106 – 154
www.elsevier.com/locate/physa
Invariant grids for reaction kinetics Alexander N. Gorbana; b; c , Iliya V. Karlina; b;∗ , Andrei Yu. Zinovyevb; c a Department
of Materials, Institute of Polymers, ETH-Zentrum, Sonneggstr. 3, ML J19, CH-8092 Z%urich, Switzerland b Institute of Computational Modeling SB RAS, Akademgorodok, Krasnoyarsk 660036, Russia c Institut des Hautes Etudes Scienti1ques, Le Bois-Marie, 35, route de Chartres, F-91440, Bures-sur-Yvette, France Received 12 June 2003; received in revised form 9 September 2003
Abstract In this paper, we construct low-dimensional manifolds of reduced description for equations of chemical kinetics from the standpoint of the method of invariant manifold (MIM). MIM is based on a formulation of the condition of invariance as an equation, and its solution by Newton iterations. A grid-based version of MIM is developed (the method of invariant grids). We describe the Newton method and the relaxation method for the invariant grids construction. The problem of the grid correction is fully decomposed into the problems of the grid’s nodes correction. The edges between the nodes appear only in the calculation of the tangent spaces. This fact determines high computational e9ciency of the method of invariant grids. The method is illustrated by two examples: the simplest catalytic reaction (Michaelis–Menten mechanism), and the hydrogen oxidation. The algorithm of analytical continuation of the approximate invariant manifold from the discrete grid is proposed. Generalizations to open systems are suggested. The set of methods covered makes it possible to e
