Explicit propagation steps

The action of the catalysts in the wet-oxidation has not been explicitly clarified. In the usual liquid-phase oxidation of hydrocarbons, the most important function of the catalyst (metal ions) is to decompose hydroperoxide (ROOH) formed during the reaction (known as Haber-Weiss mechanism) and accelerate chain propagation step by producing active radicals (ROO- and RO-) according to eq. 14.13 and 14.14. 

In this sequence R represents the active center, which normally is an atom or radical in homogeneous reactions and a catalytic site in heterogeneous reactions. Other molecules in each step are not shown explicitly, and we let the product of their concentrations with the appropriate rate constants (requiring unit stoichiometric coefficient for R) be designated by r for the initiation step, tp for the propagation step, and for the termination step. The rate equation for the active center R is... 

When solvent effects on the propagation step occnr in free-radical copolymerization reactions, they result not only in deviations from the expected overall propagation rate, but also in deviations from the ejqiected copolymer composition and microstracture. This may be trae even in bulk copolymerization, if either of the monomers exerts a direct effect or if strong cosolvency behavior causes preferential solvation. A number of models have been proposed to describe the effect of solvents on the composition, microstmcture and propagation rate of copolymerization. In deriving each of these models, an appropriate base model for copolymerization kinetics is selected (such as the terminal model or the implicit or explicit penultimate models), and a mechanism by which the solvent influences the propagation step is assumed. The main mechanisms by which the solvent (which may be one or both of the comonomers) can affect the propagation kinetics of free-radical copolymerization reactions are as follows ... 

A one-dimensional mesh through time (temporal mesh) is constructed as the calculation proceeds. The new time step is calculated from the solution at the end of the old time step. The size of the time step is governed by both accuracy and stability. Imprecisely speaking, the time step in an explicit code must be smaller than the minimum time it takes for a disturbance to travel across any element in the calculation by physical processes, such as shock propagation, material motion, or radiation transport [18], [19]. Additional limits based on accuracy may be added. For example, many codes limit the volume change of an element to prevent over-compressions or over-expansions. 

The RWP method also has features in common with several other accurate, iterative approaches to quantum dynamics, most notably Mandelshtam and Taylor s damped Chebyshev expansion of the time-independent Green s operator [4], Kouri and co-workers time-independent wave packet method [5], and Chen and Guo s Chebyshev propagator [6]. Kroes and Neuhauser also implemented damped Chebyshev iterations in the time-independent wave packet context for a challenging surface scattering calculation [7]. The main strength of the RWP method is that it is derived explicitly within the framework of time-dependent quantum mechanics and allows one to make connections or interpretations that might not be as evident with the other approaches. For example, as will be shown in Section IIB, it is possible to relate the basic iteration step to an actual physical time step. 

Architectural models explicitly specify the di.stribution of roots in space. An alternative approach, which is also useful for rhizosphere studies, is the continuum approach where only the amount of roots per unit soil volume is specified. Rules are defined that specify how roots propagate in the vertical and horizontal dimensions, and root propagation is u.sually viewed as a diffusive phenomenon (i.e., root proliferation favors unexploited soil). This defines the exploitation intensity per unit volume of soil and, under the assumption of even di.stribution, provides the necessary information for the integration step above. Acock and Pachepsky (68) provide an excellent review of the different assumptions made in the various continuum models formulated and show how such models can explain root distribution data relating to chrysanthemum. 

Hybrid MPC-MD schemes may be constructed where the mesoscopic dynamics of the bath is coupled to the molecular dynamics of solute species without introducing explicit solute-bath intermolecular forces. In such a hybrid scheme, between multiparticle collision events at times x, solute particles propagate by Newton s equations of motion in the absence of solvent forces. In order to couple solute and bath particles, the solute particles are included in the multiparticle collision step [40]. The above equations describe the dynamics provided the interaction potential is replaced by Vj(rJVs) and interactions between solute and bath particles are neglected. This type of hybrid MD-MPC dynamics also satisfies the conservation laws and preserves phase space volumes. Since bath particles can penetrate solute particles, specific structural solute-bath effects cannot be treated by this rule. However, simulations may be more efficient since the solute-solvent forces do not have to be computed. 

Despite the evident complexities of handling additional initiation or propagation reactions, the algebra of these steps is always first-order in radical concentration and the stationary-state equations, although cumbersome, can always be solved explicitly. However, as soon as we introduce more than one termination reaction involving bimolecular participation of two radicals, the equations become nonlinear and explicit solutions arc not always possible. Suppose for example that we were to include in the simple scheme Eq. (XIII.10.5), the additional termination reactions... 

A variety of explicit (Dufort-Frankel, Lax-Wendroff, Runge-Kutta) and implicit (approximate factorization, LU-SGS) or hybrid schemes have been employed for integration in time. Because of the complexity of the incompressible Navier-Stokes equations, stability analyses to determine critical time steps are difficult. As a general rule, the allowable time step for an explicit method is proportional to the ratio of the smallest grid size to the largest convective velocity (or the wave propagation speed for an artificial compressibility method). 

The convergence of the normalization algorithm may be proved if one is able to estimate the perturbation left after a given number of normalization step. This may be done with analytical estimates. Thus, we can apply a formal statement of KAM theorem to a Hamiltonian with a remainder dramatically reduced, thanks to both explicit calculation of the expansion and recursive estimates. A fully rigorous result may be achieved by performing all the calculations using interval arithmetics, so that we have full control on the propagation of roundoff errors. 

---