Kinetic algorithm

The Beeman integration scheme uses a more accurate expression for the velocity. As a consequence it often gives better energy conservation, because the kinetic energy is calculated directly from the velocities. However, the expressions used are more complex than those of the Verlet algorithm and so it is computationally more expensive. 

Computer simulation of the reactor kinetic hydrodynamic and transport characteristics reduces dependence on phenomenological representations and idealized models and provides visual representations of reactor performance. Modem quantitative representations of laminar and turbulent flows are combined with finite difference algorithms and other advanced mathematical methods to solve coupled nonlinear differential equations. The speed and reduced cost of computation, and the increased cost of laboratory experimentation, make the former increasingly usehil. 

Mechanism. The thermal cracking of hydrocarbons proceeds via a free-radical mechanism (20). Siace that discovery, many reaction schemes have been proposed for various hydrocarbon feeds (21—24). Siace radicals are neutral species with a short life, their concentrations under reaction conditions are extremely small. Therefore, the iategration of continuity equations involving radical and molecular species requires special iategration algorithms (25). An approximate method known as pseudo steady-state approximation has been used ia chemical kinetics for many years (26,27). The errors associated with various approximations ia predicting the product distribution have been given (28). 

When the kinetics are unknown, still-useful information can be obtained by finding equilibrium compositions at fixed temperature or adiabatically, or at some specified approach to the adiabatic temperature, say within 25°C (45°F) of it. Such calculations require only an input of the components of the feed and produc ts and their thermodynamic properties, not their stoichiometric relations, and are based on Gibbs energy minimization. Computer programs appear, for instance, in Smith and Missen Chemical Reaction Equilibrium Analysis Theory and Algorithms, Wiley, 1982), but the problem often is laborious enough to warrant use of one of the several available commercial services and their data banks. Several simpler cases with specified stoichiometries are solved by Walas Phase Equilibiia in Chemical Engineering, Butterworths, 1985). 

The Monte Carlo method as described so far is useful to evaluate equilibrium properties but says nothing about the time evolution of the system. However, it is in some cases possible to construct a Monte Carlo algorithm that allows the simulated system to evolve like a physical system. This is the case when the dynamics can be described as thermally activated processes, such as adsorption, desorption, and diffusion. Since these processes are particularly well defined in the case of lattice models, these are particularly well suited for this approach. The foundations of dynamical Monte Carlo (DMC) or kinetic Monte Carlo (KMC) simulations have been discussed by Eichthom and Weinberg (1991) in terms of the theory of Poisson processes. The main idea is that the rate of each process that may eventually occur on the surface can be described by an equation of the Arrhenius type ... 

A simple algorithm [17] makes it possible to find the probability of any fragment of macromolecules of Gordonian polymers. Comparison of these probabilities with the data obtained by NMR spectroscopy provides the possibility to evaluate the adequacy of a chosen kinetic model of a synthesis process of a particular polymer specimen. The above-mentioned probabilities are also involved in the expressions for the glass transition temperature and some structure-additive properties of branched polymers [18,19]. 

Timm, Gilbert, Ko, and Simmons O) presented a dynamic model for an isothermal, continuous, well-mixed polystyrene reactor. This model was in turn based upon the kinetic model developed by Timm and co-workers (2-4) based on steady state data. The process was simulated using the model and a simple steady state optimization and decoupling algorithm was tested. The results showed that steady state decoupling was adequate for molecular weight control, but not for the control of production rate. In the latter case the transient fluctuations were excessive. 

---