Kinetic equation discretized

An exhaustive statistical description of living copolymers is provided in the literature [25]. There, proceeding from kinetic equations of the ideal model, the type of stochastic process which describes the probability measure on the set of macromolecules has been rigorously established. To the state Sa(x) of this process monomeric unit Ma corresponds formed at the instant r by addition of monomer Ma to the macroradical. To the statistical ensemble of macromolecules marked by the label x there corresponds a Markovian stochastic process with discrete time but with the set of transient states Sa(x) constituting continuum. Here the fundamental distinction from the Markov chain (where the number of states is discrete) is quite evident. The role of the probability transition matrix in characterizing this chain is now played by the integral operator kernel ... 

The fast stage of relaxation of a complex reaction network could be described as mass transfer from nodes to correspondent attractors of auxiliary dynamical system and mass distribution in the attractors. After that, a slower process of mass redistribution between attractors should play a more important role. To study the next stage of relaxation, we should glue cycles of the first auxiliary system (each cycle transforms into a point), define constants of the first derivative network on this new set of nodes, construct for this new network an (first) auxiliary discrete dynamical system, etc. The process terminates when we get a discrete dynamical system with one attractor. Then the inverse process of cycle restoration and cutting starts. As a result, we create an explicit description of the relaxation process in the reaction network, find estimates of eigenvalues and eigenvectors for the kinetic equation, and provide full analysis of steady states for systems with well-separated constants. 

An auxiliary reaction network is associated with the auxiliary discrete dynamical system. This is the set of reactions A, A q with kinetic constants k,. The correspondent kinetic equation is... 

In the simplest case, the auxiliary discrete dynamical system for the reaction network W is acyclic and has only one attractor, a fixed point. Let this point be A (n is the number of vertices). The correspondent eigenvectors for zero eigenvalue are r = S j and Z = 1. For such a system, it is easy to find explicit analytic solution of kinetic equation (32). 

Until now we considered the limiting situation in which the plasma consists only of atoms. The kinetic equation was given by Eq. (4.62). However, we are also interested in describing the partially ionized plasma, and especially the ionization and recombination reactions. Formally, this corresponds to taking into account scattering states in Eq. (4.44). This means that the quantum numbers a may be discrete numbers and. may also run over the continuous spectrum. Taking into account scattering states, we are faced with difficulties. These are connected essentially with the application of Eq. (4.28) in the equation of motion. 

If we replace the quantities a, px by discrete numbers n[, we get from (4.76) a kinetic equation that describes the inelastic scattering of (4.74). 

In the second case, if the species mobilities differ greatly, the dimensionality of the system of kinetic equations decreases [103], Let all the components be divided into two groups of species a slow (5) and a rapid (r) one. This yields three types of pair functions. For the rapid species the condition of the equilibrium distribution can be considered as satisfied. Then, for the pair functions of types sr and rr instead of the kinetic Eqs. (32) algebraic relations in Appendix A apply, whose dimensionality can be lowered using the method of substitution variables according to Appendix B. In this case the kinetic Eqs (31) for the local concentrations and Eq. (32) for the pair functions type ss do not change. A similar situation remains in passing to the one-dimensional discrete and point-like models. 

Fundamentals of a method for developing models of mass transfer of low-molecular substances in non-reconstructing microheterogeneous membranes were formulated. The local properties of membranes differ in sorbability with respect to species and in the probability for a species to jump from one sorption site to another. Because of this, the permeability of a membrane depends on the amounts of different-type sites, their mutual arrangement, mutual influence of adjacent molecules, and the probabilities of jumps between different sites. The probabilities of occupation of different sorption sites are described by kinetic equations, which take into account the interactions between species. The atomic-molecular discrete and continuous models of mass transfer for thin and thick films are constructed. 

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). 

Consider first a discrete monomolecular system in which every species present can transform to any other one by a first-order reaction. In such a system, the kinetic equation at all possible compositions is given by Eq. (70), which is repeated here ... 

Equation (98) should be discussed in some detail. When a approaches < (which is the case where only reactant x = 1 is present in the mixture), it correctly predicts the single-reactant result dCldt = —C. This is worded by saying that the alias satisfies the single-component identity (SCI) When the initial concentration distribution approaches a delta function, one recovers the result for a single component. It is clear that the SCI requirement must be satisfied for every kinetic equation, not only the linear one. In fact, one can generalize the requirement to that of the discrete component identity—when the initial concentration distribution approaches the sum of N distinct delta functions, one must recover the corresponding discrete description for A components (Aris, 1991a). 

When the intrinsic kinetics are nonlinear, some interesting problems arise that are best discussed first with a discrete description. A possible assumption for the kinetics is that they are independent, that is, that the rate at which component 7 disappears depends only on the concentration of component I itself (this assumption is obviously correct in the first-order case). The difficulties associated with the assumption of independence are best illustrated by considering the case of parallel th order reactions, which has been analyzed by Luss and Hutchinson (1971), who write the kinetic equation for component / as -dcj/dt = kjc", 1=1, 2,. . ., N, where Cj is the (dimensional) concentration of component / at time t. The total initial concentration C(0) is , and this is certainly finite. Now consider the following special, but perfectly legitimate case. The value of C(0) is fixed, and the initial concentrations of all reactants are equal, so that C/(0) = C(0)/N. Furthermore, all the k/ s are equal to each other, k/ = k. One now obtains, for the initial rate of decrease of the overall concentration ... 

Quantitative approaches to describing reactions in micelles differ markedly from treatments of reactions in homogeneous solution primarily because discrete statistical distributions of reactants among the micelles must be used in place of conventional concentrations [74], Further, the kinetic approach for bimolecular reactions will depend on how the reactants partition between micelles and bulk solution, and where they are located within the microphase region. Distinct microphase environments have been sensed by NMR spectrometry for hydrophobic molecules such as pyrene, cyclohexane and isopropylbenzene, which are thought to lie within a hydrophobic core , and less hydrophobic molecules such as nitrobenzene and N,N-dimethylaniline, which are preferentially located at the micelle-water interface [75]. Despite these complexities, relatively simple kinetic equations for electron-transfer reactions can be derived for cases where both donors and acceptors are uniformly distributed inside the micelle or on its surface. 

Thus the Debye equation [Eq. (1)] may be satisfactorily explained in terms of the thermal fluctuations of an assembly of dipoles embedded in a heat bath giving rise to rotational Brownian motion described by the Fokker-Planck or Langevin equations. The advantage of a formulation in terms of the Brownian motion is that the kinetic equations of that theory may be used to extend the Debye calculation to more complicated situations [8] involving the inertial effects of the molecules and interactions between the molecules. Moreover, the microscopic mechanisms underlying the Debye behavior may be clearly understood in terms of the diffusion limit of a discrete time random walk on the surface of the unit sphere. 

We shall now almost exclusively concentrate on the fractal time random walk excluding inertial effects and the discrete orientation model of dielectric relaxation. We shall demonstrate how in the diffusion limit this walk will yield a fractional generalization of the Debye-Frohlich model. Just as in the conventional Debye relaxation, a fractional generalization of the Debye-Frohlich model may be derived from a number of very different models of the relaxation process (compare the approach of Refs. 22, 23, 28 and 34—36). The advantage of using an approach based on a kinetic equation such as the fractional Fokker-Planck equation (FFPE) however is that such a method may easily be extended to include the effects of the inertia of the dipoles, external potentials, and so on. Moreover, the FFPE (by use of a theorem of operational calculus generalized to fractional exponents and continued fraction methods) clearly indicates how many existing results of the classical theory of the Brownian motion may be extended to include fractional dynamics. 

The answer to this question is mainly driven by the computational cost of solving the kinetic equation due to the large number of independent variables. In the simplest example of a 3D velocity-distribution function n t, x, v) the number of independent variables is 1 + 3 + 3 = 1. However, for polydisperse multiphase flows the number of mesoscale variables can be much larger than three. In comparison, the moment-transport equations involve four independent variables (physical space and time). Furthermore, the form of the moment-transport equations is such that they can be easily integrated into standard computational-fluid-dynamics (CFD) codes. Direct solvers for the kinetic equation are much more difficult to construct and require specialized numerical methods if accurate results are to be obtained (Filbet Russo, 2003). For example, with a direct solver it is necessary to discretize all of phase space since a priori the location of nonzero values of n is unknown, which can be very costly when phase space is not bounded. 

Equation (7.54) holds at the equilibrium and when the initial NDF is isotropic. In order to illustrate how the discretization for the kinetic equation is typically carried out, some additional simplifications must be introduced. It is convenient to write the collision integral as... 

At elastic scatter collisions and when a discrete energy level is excited in an (n, n ) reaction, the standard kinetic equations may be used. The following equations are used to calculate the scatter angle (j) in the laboratory system from the scatter angle 6 in the center of mass system, and the emergent energy E in terms of the incident energy . 

LBM has been developed to simulate flows in microchannels based on kinetics equations and statistical physics. In LBM, the motion of the fluid is modeled by a lattice-Boltzmann equation for the distribution function of the fluid molecules. The discrete velocity Boltzmann equation corresponding to the Navier-Stokes equations can be written as... 

Action of which mass The naming law of mass action is incorrect for several reasons. First, from the viewpoint of chemistry, which in modern literature (although not always) prefers the naming law of reaction equilibrium, because the reference to the mass of a chemical substance into kinetic equations or into equilibrium equations is abandoned for long to the benefit of the number of moles (which is a measure of the quantity of substance based on a counting of discretized entities). 

---