Conservation matrices chemical reactions

The results presented above were discussed in terms of the special case of elementary reactions. However, if we relax the condition that the coefficients vfai and uTai must be integers, (5.1) is applicable to nearly all chemical reactions occurring in practical applications. In this general case, the element conservation constraints are no longer applicable. Nevertheless, all of the results presented thus far can be expressed in terms of the reaction coefficient matrix T, defined as before by... 

In Section 5.1, we have seen (Fig. 5.2) that the molar concentration vector c can be transformed using the SVD of the reaction coefficient matrix T into a vector c that has Nr reacting components cr and N conserved components cc.35 In the limit of equilibrium chemistry, the behavior of the Nr reacting scalars will be dominated by the transformed chemical source term S. 36 On the other hand, the behavior of the N conserved scalars will depend on the turbulent flow field and the inlet and initial conditions for the flow domain. However, they will be independent of the chemical reactions, which greatly simplifies the mathematical description. 

George and Ross34 set out to derive symmetry rules for chemical reactions as a set of selection rules on elements of the transition matrix. Each element of this matrix describes the probability of transition from a specified state of the reactants to a specified state of the products. One selection rule on such a matrix is the approximate conservation of total electron spin by making the Born-Oppenheimer approxima-... 

This set of chemical reactions is not unique for example, the reference reaction can be written with H2P04. Additional reactions are involved if Mg2 + or other cations are bound reversibly by these species. The conservation matrix for this... 

Systems of biochemical reactions like glycolysis, the citric acid cycle, and larger and smaller sequential and cyclic sets of enzyme-catalyzed reactions present challenges to make calculations and to obtain an overview. The calculations of equilibrium compositions for these systems of reactions are different from equilibrium calculations on chemical reactions because additional constraints, which arise from the enzyme mechanisms, must be taken into account. These additional constraints are taken into account when the stoichiometric number matrix is used in the equilibrium calculation via the program equcalcrx, but they must be explicitly written out when the conservation matrix is used with the program equcalcc. The stoichiometric number matrix for a system of reactions can also be used to calculate net reactions and pathways. 

Glycolysis involves 10 biochemical reactions and 16 reactants. Water is not counted as a reactant in writing the stoichiometric number matrix or the conservation matrix for reasons described in Section 6.3. Thus there are six components because C = N — R = 16 — 10 = 6. From a chemical standpoint this is a surprise because the reactants involve only C, H, O, N, and P. Since H and O are not conserved at specified pH in dilute aqueous solution, there are only three conservation equations based on elements. Thus three additional conservation relations arise from the mechanisms of the enzyme-catalyzed reactions in glycolysis. Some of these conservation relations are discussed in Alberty (1992a). At specified pH in dilute aqueous solutions the reactions in glycolysis are... 

The principle of mass conservation, on a macroscopic basis, using equations representing the conservation of total mass, the conservation of each chemical element and the mass balance for each of the molecular species was applied. Chemical reactions were taken into account in these equations. The stoichiometric coefficients of individual reactants were evaluated using Microsoft Excel software by expressing the equations in matrix notation. 

Because the vector m is constrained by the mass conservation requirement = const, the space of possible m values has — 1 dimensions. If the number of independent chemical reactions, R, is less than - 1, then some vectors m are not accessible at some assigned M this, as will be seen, has important consequences in the consideration of heterogeneous chemical equilibria. Now consider the special case where R = N — 1, so that indeed all admissible m s are accessible. Because the kernel of a contains only the zero vector, there exists an M X N matrix A such that... 

One way to recognize the significance of this equation is to remember that the ultimate objective of chemical thermodynamics is to calculate the equilibrium composition of a system of reactions. A chemical reaction system has R independent equilibrium constant expressions and C conservation equations, and this is just enough information to calculate the equilibrium concentrations of N species. Equation 7.1-9 is useful because it makes it possible to calculate a conservation matrix from a stoichiometric number matrix. In doing this with the operation NullSpace we will see again that it yields a basis for the conservation matrix. 

It is important to notice that row reduction changes the components from atoms of elements to combinations of atoms. The first three species are selected as components if they contain all the different atoms. Thus a set of conservation equations can be written to conserve the first three species, rather than the atoms of C, H, and O. When row reduction yields a matrix of this form, the chemical reaction can be read from the last column. This shows how H2 O is made up from the three components H2 O = (1/2)CH4 -I- O2 - (1/2)C02. This can be rearranged to give equation 7.1-1. 

The central concept involved in coupling is the identification of components, which are the things that are conserved in a reaction system. When chemical reactions are studied, atoms of elements are conserved, but some of these conservation equations may not be independent. Redundant conservation equations are not counted as components C. When the pH is specified, the conservation equation for hydrogen atoms is omitted, and so the number of components for a given system is reduced by one C = C - 1. A test of the conservation matrix A is that the equation A v = 0 must yield a suitable basis for the stoichiometric number matrix v. When it is necessary to recognize that oxygen atoms are available from h2o. A must be used, and C " = C - 1. A test of the conservation matrix A" is that the equation A V " = 0 must yield a suitable basis for the stoichiometric number matrix v. ... 

Equilibrium compositions of systems of chemical reactions or systems of enzyme-catalyzed reactions can only be calculated by iterative methods, like the Newton-Raphson method, and so computer programs are required. These computer programs involve matrix operations for going back and forth between conservation matrices and stoichiometric number matrices. A more global view of biochemical equilibria can be obtained by specifying steady-state concentrations of coenzymes. These are referred to as calculations at the third level to distinguish them from the first level (chemical thermodynamic calculations in terms of species) and the second level (biochemical thermodynamic calculations at specified pH in terms of reactants). 

The conservation of chemical mass in the subsurface is mainly controlled by the transport of solutes through the water-saturated porous medium, the chemical reactions in the groundwater and the reactions between the groundwater plus solutes and the solid matrix. 

In this section, we first introduce the standard form of the chemical source term for both elementary and non-elementary reactions. We then show how to transform the composition vector into reacting and conserved vectors based on the form of the reaction coefficient matrix. We conclude by looking at how the chemical source term is affected by Reynolds averaging, and define the chemical time scales based on the Jacobian of the chemical source term. 

For elementary reactions (Hill 1977), the values of the stoichiometric coefficients are constrained by the fact that all chemical elements must be conserved in (5.1). Mathematically, this can be expressed in terms of an E x K element matrix A where E is the total number of chemical elements present in the reacting flow. Each column of A thus corresponds to a particular chemical species, and each row to a particular chemical element. As an example, consider a system containing E = 2 elements O and H, and K = 3 species H2, O2 and H20. The 2x3 element matrix for this system is... 

In a long series of papers on the master equation, Pritchard and his coworkers elucidated for the first time the effects of rotational and vibrational disequilibrium on the dissociation and recombination of a dilute diatomic gas. Ultrasonic dispersion in a diatomic gas was analyzed by similar computational experiments, and the first example of the breakdown of the linear mixture rule in chemical kinetics was demonstrated. A major difficulty in these calculations is that the eigenvalue of the reaction matrix (corresponding to the rate constant) differs from the zero eigenvalue (required by species conservation) by less than... 

We normally use a full set of NS material balances, even though the number of independent balances is limited to the rank NK of the stoichiometric matrix u for the given reaction scheme as shown in Section 2.1. cmd in Aris (1969). Exceptions must be made, however, when one or more constraints are imposed, such as quasi-equilibrium for some reactions or pseudo-steady state (better called quasi-conservation) for some chemical species then each active constraint will replace a mass balance. By these procedures, we avoid catastrophic cancellations that might occur in subtractions performed to reduce the number of species variables from NS to NK. 

---