Reaction front prediction

In a reactive transport model, the domain of interest is divided into nodal blocks, as shown in Figure 2.11. Fluid enters the domain across one boundary, reacts with the medium, and discharges at another boundary. In many cases, reaction occurs along fronts that migrate through the medium until they either traverse it or assume a steady-state position (Lichtner, 1988). As noted by Lichtner (1988), models of this nature predict that reactions occur in the same sequence in space and time as they do in simple reaction path models. The reactive transport models, however, predict how the positions of reaction fronts migrate through time, provided that reliable input is available about flow rates, the permeability and dispersivity of the medium, and reaction rate constants. 

Empirical Methods. The grcphical deactivation plot is a very useful empirical method for prediction of the catalyst performance and for estimation of catalyst lifetime (18,19). The deactivation plot shows the length of the reaction front as a function of time. This illustrates the movement of the temperature profile caused by the progressive deactivation of the catalyst. The method is illustrated in Figure 3. The temperature increase over the catalyst bed is calculated as AT = Texii - Twet and a certain percentage hereof, e.g. 90% (AT90) is calculated. The axial distance in the... 

Co and Ci are the continuous phase solute concentrations in teed and inside the mixer, respectively, and are the feed rates of the continuous and emulsion phases, respectively, Cir is the internal reagent concentration (based on volume of emulsion), R is the emulsion globule radius, Kr is the emulsion phase holdup volume in the mixer, a is the partition coefficient for solute between external phase and emulsion, is the effective solute diffusivity in the emulsion, and % is the dimensionless reaction front position. Hatton and Wardius [51] also extended their analysis to develop simple graphical and numerical procedures for the prediction of multistage extraction performance of mixer-settler trains operating either cocurrendy or countercurrently without any external recycle over individual stages. For a typical stage i in a cocurrent mixer-settler, they defined the parameter 6 as... 

During the oxidation of CO, CH4 and C3H8, the ignited state is characterized by a reaction front stabilized in a thin portion of the bed near the reactor inlet. This condition, corresponding to a diffusion-controlled reaction, is predicted by the known models of exothermic catalytic reactions [4], The chemical factors determining this dynamics are the heat of reaction and the activation energy. For all of the reactants considered in this study, a similar behaviour in the ignited state is observed. 

The same chemical factors that determine the mobility of the reaction front, seem also to play a role in the mechanism of oscillations. These results can help define general criteria for predicting unstable behaviour in processes carried out in the ignited state. 

Reversibility Effects. Figure 2 demonstrates that differences between advancing front and reversible reaction model predictions are significant when oj Is less than 1 or when oj Is greater than 10. When o5 Is small, as measured by the deviation of (1+o )/o from 1, the solute concentration Is too small to force the reagent to react completely. Reaction reversibility causes the globule extraction capacity to depend on the outlet solute concentration. 

For oj equal to 10 or larger, the advancing front and reversible reaction models predict identical curves which depend only on the ratio of Oj/oJ. In this situation, extraction curves like those in Figure 2 can be developed by modifying the pseudosteady-state advancing front calculation slightly ... 

Figures 6-8 compare adjusted and nonadjusted advancing front predictions with the small oj approximation and the numerically-solved reversible reaction model for o,/oJ of 200 and oj values of 1.0, 0.5 and 0.1. The nonadjusted advancing front curve is the same in all three plots, because the advancing front solution depends only on the ratio o,/oS, which is fixed. When ojis 0.1, differences between the small o and reversible reaction curves are small and within the error of the numerical solution. As expected, the small o approximation is poor when oj is 1 and the condition that (1 + oS) 1 no longer applies. In this situation, the adjusted advancing front approach can be used to estimate reasonable design parameters. For oj of 0.5, the adjusted advancing front and small oj predictions bracket the reversible reaction solution.

It is interesting to note that calculations of turbulent flows during fast chemical reactions, predicted that the chemical reaction rate constant influences the effective diffusion coefficient and accelerates micromixing, due to an increase of the local reactant concentration gradients [13]. The dependence of the lower boundaries of the reaction front macrostructure formation, in particular, the plane and the torch front, which characterise different scales of liquid flow mixing, on the values of the chemical reaction constants is experimental evidence of the correlation between the kinetic and diffusive parameters of the process. At the same time, one can suppose that the formation of the characteristic reaction front macrostructures is defined by the mixing at the macro- and microlevels. 

In the absence of a planar reaction front, we cannot, in general, predict the rate of selective oxidation of alloys. The criterion (9.72) is nevertheless useful, because it allows us to understand qualitatively the influence of different parameters on the oxidation behavior of binary alloys. 

Sorbed pesticides are not available for transport, but if water having lower pesticide concentration moves through the soil layer, pesticide is desorbed from the soil surface until a new equiUbrium is reached. Thus, the kinetics of sorption and desorption relative to the water conductivity rates determine the actual rate of pesticide transport. At high rates of water flow, chances are greater that sorption and desorption reactions may not reach equihbrium (64). NonequiUbrium models may describe sorption and desorption better under these circumstances. The prediction of herbicide concentration in the soil solution is further compHcated by hysteresis in the sorption—desorption isotherms. Both sorption and dispersion contribute to the substantial retention of herbicide found behind the initial front in typical breakthrough curves and to the depth distribution of residues. 

This equation has been derived as a model amplitude equation in several contexts, from the flow of thin fluid films down an inclined plane to the development of instabilities on flame fronts and pattern formation in reaction-diffusion systems we will not discuss here the validity of the K-S as a model of the above physicochemical processes (see (5) and references therein). Extensive theoretical and numerical work on several versions of the K-S has been performed by many researchers (2). One of the main reasons is the rich patterns of dynamic behavior and transitions that this model exhibits even in one spatial dimension. This makes it a testing ground for methods and algorithms for the study and analysis of complex dynamics. Another reason is the recent theory of Inertial Manifolds, through which it can be shown that the K-S is strictly equivalent to a low dimensional dynamical system (a set of Ordinary Differentia Equations) (6). The dimension of this set of course varies as the parameter a varies. This implies that the various bifurcations of the solutions of the K-S as well as the chaotic dynamics associated with them can be predicted by low-dimensional sets of ODEs. It is interesting that the Inertial Manifold Theory provides an algorithmic approach for the construction of this set of ODEs. 

The conditions in the reaction zone determine the release rate of the N-precursors HCN and NH3 [11], Among these conditions are the properties of the fuel (e.g., N-content and particle size), parameters related to the combustion front (temperature and propagation velocity) and the gas composition in and directly above the combustion front. As the prediction of the mass fractions of the N-precursors is important for the final goal of this research, i.e., the prediction of NO formation of the complete furnace, a model is needed in which all these conditions are represented. 

---