Simulations kinetic theory

A further paper [167] explains the lamellar thickness selection in the row model. The minimum thickness lmin is derived from the similation and found to be consistent with equilibrium results. The thickness deviation 81 = l — lmin is approximately constant with /. It is established that the model fulfills the criteria of a kinetic theory Firstly, a driving force term (proportional to 81) and a barrier term (proportional to /) are indentified. Secondly, the competition between the two terms leads to a maximum in growth rate (see Fig. 2.4) which is located at the average thickness l obtained by simulation. Further, the role of fluctuations becomes apparent when the dependence on the interaction energy e is investigated. Whereas downwards (i.e. decreasing l) fluctuations are approximately independent... 

The theory was very similar to that described earlier, but was simplified in view of the complexity of the problem. A number of reaction intermediates were considered explicitly, and the corresponding signals were calculated by molecular dynamics simulation. Kinetic equations governing the reaction sequence were established and were solved numerically. The main simplification of the theory is that, when calculating A5[r, r], the lower limit of the Fourier integral was shifted from 0 to a small value q. The authors wrote [59]... 

This relation is consistent with previous results observed experimentally [49]. Although the kinetic theory of Lauritzen and Hoffman predicts the same law as Eq. 16, it predicts a divergence in L at lower undercoolings. The simulations do not show any evidence for such a catastrophe. 

Fig. 29. Snapshots of particle volume fraction fields obtained while solving a kinetic theory-based TFM. 75 pm fluid catalytic particles in ambient air. Simulations were done over a 16 x 32 cm periodic domain. The average particle volume fraction in the domain is 0.05. Dark (light) color indicates regions of high (low) particle volume fractions. (See Refs. Agrawal et al., 2001 Andrews et al., 2005) for other parameter values.) Source Andrews and Sundaresan (2005).

Fig. 32. Filtered particle phase pressure (in CGS units) extracted from simulations over 16 x 16 cm domain using 128 x 128 cells. Source Andrews and Sundaresan (2005). The filtered particle-phase pressure includes the Reynolds stress-like fluctuations and the kinetic theory pressure.

The theory was very similar to that described earlier but was simplified in view of the complexity of the problem. A number of reaction intermediates were considered explicitly, and the corresponding signals were calculated by molecular dynamics simulation. Kinetic equations governing the reaction... 

We would like to conclude this introductory Chapter by the following general comment. Most of the papers dealing with the fluctuation-controlled reactions, focus their attention on the simplest bimolecular A + B —> B and A + B —> 0 reactions. To our mind, main results in this field are already obtained and the situation is quite clear. In the nearest future the most prospective direction of kinetic theory seems to be many-stage catalytic processes the first results are discussed in Chapters 8 and 9. Their study (stimulated also by the technological importance) should be continued using in parallel both refined mathematical formalisms of the fluctuation-controlled kinetics and full-scale computer simulations. 

In the formulations developed from the renormalized kinetic theory approach, these self-consistencies were avoided either by using values obtained from computer simulation and experiments or by using some exactly known limiting values for the transport coefficient. For example, in the treatment of Mazenko [5-7], and of Mehaffey and Cukier s [8] the transport coefficients are replaced by their Enskog values. In the theory developed by Sjogren and Sjolander [9], the velocity autocorrelation function is required as an input that was obtained from the computer simulated values. This limits the validity of the theories only to certain regimes and for certain systems where the experimental or computer-simulated results are available. 

The above-mentioned computer simulation and experimental studies have addressed various aspects of mass dependence, but they all show that the selfdiffusion coefficient of a tagged molecule exhibits a weak mass dependence, especially for solutes with size comparable to or larger than the size of the solvent molecules. Sometimes this mass dependence can be fitted to a power law, with a small exponent less than 0.1 [99]. This weak mass dependence has often been considered as supportive of the hydrodynamic picture. In hydrodynamics the diffusion of a solute is conventionally described by the well-known Stokes-Einstein (SE) relation, which predicts that the diffusion is totally independent of the mass of the solute. Kinetic theory, on the other... 

According to the SE relation, the product Drj should remain constant for systems having particles of same size and studied at the same temperature. Recent studies [102, 104, 105] have found that the SE relation does not hold when the mass of the particles are changed. Walser et al. [104] have performed MD simulations of water molecules with different mass and different molecular mass distribution. They have shown that although the viscosity increases and the diffusion decreases with mass, the product of the two does not remain constant. They have found that the product Dr) is not correlated with the molecular mass, but it is correlated for those systems with the same mass distribution. Thus, while the mass dependence is not as strong as predicted by the kinetic theory, it is also not totally negligible. 

Ray Kapral came to Toronto from the United States in 1969. His research interests center on theories of rate processes both in systems close to equilibrium, where the goal is the development of a microscopic theory of condensed phase reaction rates,89 and in systems far from chemical equilibrium, where descriptions of the complex spatial and temporal reactive dynamics that these systems exhibit have been developed.90 He and his collaborators have carried out research on the dynamics of phase transitions and critical phenomena, the dynamics of colloidal suspensions, the kinetic theory of chemical reactions in liquids, nonequilibrium statistical mechanics of liquids and mode coupling theory, mechanisms for the onset of chaos in nonlinear dynamical systems, the stochastic theory of chemical rate processes, studies of pattern formation in chemically reacting systems, and the development of molecular dynamics simulation methods for activated chemical rate processes. His recent research activities center on the theory of quantum and classical rate processes in the condensed phase91 and in clusters, and studies of chemical waves and patterns in reacting systems at both the macroscopic and mesoscopic levels. 

The simulator used was a DISMOL, described previously by Batistella and Maciel (2). All explanations of the equations used, the solution methods, and the routine of solution are described in Batistella and Maciel (5). DISMOL is a simulator that permits changes in feed composition, feed temperaturethe evaporation rate, as well as feed flow rate. The effective rate of surface evaporation is obtained from the kinetic theory of gases. The liquid film thickness is obtained by mass balance and geometry of the evaporator. The temperature in the liquid obeys the Fourier-Kirchhoff equation. The solution of the velocity profile requires knowledge of the viscosity and the liquid film thickness over the evaporator. The solution for the temperature and the concentration profiles requires knowledge of the velocity profiles, which determine the convective heat and mass fluxes. 

Correlation Between Monte-Carlo Simulations and Kinetic Theory. 435... 

Except the kinetic equations, now various numerical techniques are used to study the dynamics of surfaces and gas-solid interface processes. The cellular automata and MC techniques are briefly discussed. Both techniques can be directly connected with the lattice-gas model, as they operate with discrete distribution of the molecules. Using the distribution functions in a kinetic theory a priori assumes the existence of the total distribution function for molecules of the whole system, while all numerical methods have to generate this function during computations. A success of such generation defines an accuracy of simulations. Also, the well-known molecular dynamics technique is used for interface study. Nevertheless this topic is omitted from our consideration as it requires an analysis of a physical background for construction of the transition probabilities. This analysis is connected with an oscillation dynamics of all species in the system that is absent in the discussed kinetic equations (Section 3). 

In the review information only about the first steps of MC simulation is given as today this method is dominant by comparison with the kinetic theory. The calculations based on the dynamic MC methods for the lattice-gas model are carried out using the master equation (24). The calculation results depend appreciably on the way of assigning the probabilities of transitions Wa. This was repeatedly pointed out in applying both the cluster methods (Section 3) and the MC method (see, e.g. Ref. [269]). Nevertheless, practically in all the papers of Section 7 the expressions (29) and (30) do not take into account the interaction between AC and its neighbors (i.e., the collision model was used). It means s (r) = 0, whereas analysis of the cluster simulations demonstrated important influence of the parameter s (r) (that restricts obtained MC results). 

We proceed now to the problems (Problem 2) and (Problem 3). At least two levels of description are involved in direct molecular simulations. The first one is the level of the np-particle kinetic theory and the second is the level of fluid mechanics on which the external forces and the final results that we seek are formulated. We shall use the multiscale formulation developed above and combine the two levels. The two levels that we consider in this section are... 

William Russel May I follow up on that and sharpen the issue a bit In the complex fluids that we have talked about, three types of nonequilibrium phenomena are important. First, phase transitions may have dynamics on the time scale of the process, as mentioned by Matt Tirrell. Second, a fluid may be at equilibrium at rest but is displaced from equilibrium by flow, which is the origin of non-Newtonian behavior in polymeric and colloidal fluids. And third, the resting state itself may be far from equilibrium, as for a glass or a gel. At present, computer simulations can address all three, but only partially. Statistical mechanical or kinetic theories have something to say about the first two, but the dynamics and the structure and transport properties of the nonequilibrium states remain poorly understood, except for the polymeric fluids. 

Professor Wakeham is interested in the relationship between the bulk thermophysical properties of fluids and the intermolecular forces between the molecules that comprise them. Thus, at one extreme, he is involved in the determination of intermolecular forces from measurements of macroscopic properties and the development and application of the statistical mechanics and kinetic theory that interrelate them. He is also actively involved in the measurement of the thermophysical properties of fluids under a very wide variety of thermodynamic states. The same thermophysical properties find application in the process industries within the design of a plant. A part of Professor Wakeham s activities are therefore concerned with the representation and extension of a body of accurate information on thermophysical properties in a fashion that allows their use with software packages for process simulation. 

---