Chemical reactors equilibrium constant method

A useful empirical approach to the design of heterogeneous chemical reactors often consists of selecting a suitable equation, such as one in Table 3.3 which, with numerical values substituted for the kinetic and equilibrium constants, represents the chemical reaction in the absence of mass transfer effects. Graphical methods are often employed to aid the selection of an appropriate equation140 and the constants determined by a least squares approach<40). It is important to stress, however, that while the equation selected may well represent the experimental data, it does not... 

The concentration of any of these species depends on the total concentration of dissolved aluminum and on the pH, and this makes the system complex from the mathematical point of view and consequently, difficult to solve. To simplify the calculations, mass balances were applied only to a unique aluminum species (the total dissolved aluminum, TDA, instead of the several species considered) and to hydroxyl and protons. For each time step (of the differential equations-solving method), the different aluminum species and the resulting proton and hydroxyl concentration in each zone were recalculated using a pseudoequilibrium approach. To do this, the equilibrium equations (4.64)-(4.71), and the charge (4.72), the aluminum (4.73), and inorganic carbon (IC) balances (4.74) were considered in each zone (anodic, cathodic, and chemical), and a nonlinear iterative procedure (based on an optimization method) was applied to satisfy simultaneously all the equilibrium constants. In these equations (4.64)-(4.74), subindex z stands for the three zones in which the electrochemical reactor is divided (anodic, cathodic, and chemical). 

---