Velocity probability distribution

This method has been devised as an effective numerical teclmique of computational fluid dynamics. The basic variables are the time-dependent probability distributions f x, f) of a velocity class a on a lattice site x. This probability distribution is then updated in discrete time steps using a detenninistic local rule. A carefiil choice of the lattice and the set of velocity vectors minimizes the effects of lattice anisotropy. This scheme has recently been applied to study the fomiation of lamellar phases in amphiphilic systems [92, 93]. 

If Restart is not checked then the velocities are randomly assigned in a way that leads to a Maxwell-Boltzmann distribution of velocities. That is, a random number generator assigns velocities according to a Gaussian probability distribution. The velocities are then scaled so that the total kinetic energy is exactly 12 kT where T is the specified starting temperature. After a short period of simulation the velocities evolve into a Maxwell-Boltzmann distribution. 

In its more advanced aspects, kinetic theory is based upon a description of the gas in terms of the probability of a particle having certain values of coordinates and velocity, at a given time. Particle interactions are developed by the ordinary laws of mechanics, and the results of these are averaged over the probability distribution. The probability distribution function that is used for a given macroscopic physical situation is determined by means of an equation, the Boltzmann transport equation, which describes the space, velocity, and time changes of the distribution function in terms of collisions between particles. This equation is usually solved to give the distribution function in terms of certain macroscopic functions thus, the macroscopic conditions imposed upon the gas are taken into account in the probability function description of the microscopic situation. 

This is the probability of finding particle 1 with coordinate rx and velocity vx (within drx and dVj), particle 2 with coordinate r2 and velocity v2 (within phase space with velocity rather than momentum for convenience since only one type of particle is being considered, this causes no difficulties in Liouville s equation.) The -particle probability distribution function ( < N) is... 

The velocity gradient leads to an altered distribution of configuration. This distortion is in opposition to the thermal motions of the segments, which cause the configuration of the coil to drift towards the most probable distribution, i.e. the equilibrium s configurational distribution. Rouse derivations confirm that the motions of the macromolecule can be divided into (N-l) different modes, each associated with a characteristic relaxation time, iR p. In this case, a generalised Maxwell model is obtained with a discrete relaxation time distribution. 

Fig. 2.8.13 Velocity probability distributions at different positions across the gap in a 5 mm-9 mm Couette cell at a shear of 0.101 s-1 and following long pre-shearing at high shear rate.

The RTD quantifies the number of fluid particles which spend different durations in a reactor and is dependent upon the distribution of axial velocities and the reactor length [3]. The impact of advection field structures such as vortices on the molecular transit time in a reactor are manifest in the RTD [6, 33], MRM measurement of the propagator of the motion provides the velocity probability distribution over the experimental observation time A. The residence time is a primary means of characterizing the mixing in reactor flow systems and is provided directly by the propagator if the velocity distribution is invariant with respect to the observation time. In this case an exact relationship between the propagator and the RTD, N(t), exists... 

Two template examples based on a capillary geometry are the plug flow ideal reactor and the non-ideal Poiseuille flow reactor [3]. Because in the plug flow reactor there is a single velocity, v0, with a velocity probability distribution P(v) = v0 16 (v - Vo) the residence time distribution for capillary of length L is the normalized delta function RTD(t) = T 1S(t-1), where x = I/v0. The non-ideal reactor with the para-... 

Detailed modeling study of practical sprays has a fairly short history due to the complexity of the physical processes involved. As reviewed by O Rourke and Amsden, 3l() two primary approaches have been developed and applied to modeling of physical phenomena in sprays (a) spray equation approach and (b) stochastic particle approach. The first step toward modeling sprays was taken when a statistical formulation was proposed for spray analysis. 541 Even with this simplification, however, the mathematical problem was formidable and could be analyzed only when very restrictive assumptions were made. This is because the statistical formulation required the solution of the spray equation determining the evolution of the probability distribution function of droplet locations, sizes, velocities, and temperatures. The spray equation resembles the Boltzmann equation of gas dynamics[542] but has more independent variables and more complex terms on its right-hand side representing the effects of nucleations, collisions, and breakups of droplets. 

In the second method, i.e., th particle method 546H5471 a spray is discretized into computational particles that follow droplet characteristic paths. Each particle represents a number of droplets of identical size, velocity, and temperature. Trajectories of individual droplets are calculated assuming that the droplets have no influence on surrounding gas. A later method, 5481 that is restricted to steady-state sprays, includes complete coupling between droplets and gas. This method also discretizes the assumed droplet probability distribution function at the upstream boundary, which is determined by the atomization process, by subdividing the domain of coordinates into computational cells. Then, one parcel is injected for each cell. 

We start by considering an arbitrary measurable10 one-point11 scalar function of the random fields U and 0 Q U, 0). Note that, based on this definition, Q is also a random field parameterized by x and t. For each realization of a turbulent flow, Q will be different, and we can define its expected value using the probability distribution for the ensemble of realizations.12 Nevertheless, the expected value of the convected derivative of Q can be expressed in terms of partial derivatives of the one-point joint velocity, composition PDF 13... 

Before we present our idea, we give a short summary of the mean-field theory. The core concept there is the velocity (v)-position (r) probability distribution of stars,/(r, v, t), a positive and integrable function that is a priori time-dependent. The number density p(r, t) is the integral over velocities of/(r, v, t) ... 

Figure 3.1 Energy levels and wave functions of harmonic oscillator. Heavy line bounding potential (3.2). Light solid lines quantum-mechanic probability density distributions for various quantum vibrational numbers see section 1.16.1). Dashed lines classical probability distribution maximum classical probability is observed in the zone of inversion of motion where velocity is zero. From McMillan (1985). Reprinted with permission of The Mineralogical Society of America.

In this section, we begin the description of Brownian motion in terms of stochastic process. Here, we establish the link between stochastic processes and diffusion equations by giving expressions for the drift velocity and diffusivity of a stochastic process whose probability distribution obeys a desired diffusion equation. The drift velocity vector and diffusivity tensor are defined here as statistical properties of a stochastic process, which are proportional to the first and second moments of random changes in coordinates over a short time period, respectively. In Section VILA, we describe Brownian motion as a random walk of the soft generalized coordinates, and in Section VII.B as a constrained random walk of the Cartesian bead positions. 

An equivalent dehnition of the drift velocity E" may be obtained by using the diffusion equation alone to calculate the average flux velocity in a statishcal ensemble characterized by a probability distribution... 

Here, is the mobility tensor in the chosen system of coordinates, which is a constrained mobility for a constrained system and an unconstrained mobility for an unconstrained system. As discussed in Section VII, in the case of a constrained system, Eq. (2.344) may be applied either to the drift velocities for the / soft coordinates, for which is a nonsingular / x / matrix, or to the drift velocities for a set of 3N unconstrained generalized or Cartesian coordinates, for a probability distribution (X) that is dynamically constrained to the constraint surface, for which is a singular 3N x 3N matrix. The equilibrium distribution is. (X) oc for unconstrained systems and... 

Let us now consider a NESS where the trap is moved at constant velocity, x t) = vt. It is not possible to solve the Fokker-Planck equation to find the probability distribution in the steady state for arbitrary potentials. Only for harmonic potentials, U x) = ioc /2, can the Fokker-Planck equation be solved exactly. The result is... 

Most probable settling velocity from sedimentation data Particle-size determination from sedimentation equation Sedimentation in an ultracentrifuge Solvation and ellipticity from sedimentation data Diffusion and Gaussian distribution Temperature-dependence of diffusion coefficients... 

For compounds with Kii/Vl larger than about 10 2 the overall air-water transfer velocity is approximately equal to the water-phase exchange velocity viw The latter is related to wind speed uw by a nonlinear relation (Table 20.2, Eq. 20-16). The annual mean of viw calculated from Eq. 20-16 with the annual mean wind speed ul0 would underestimate the real mean air-water exchange velocity. Thus, we need information not only on the average wind speed, but also on the wind-speed probability distribution. 

The velocity probability distribution function of Eq. 10.20 is the well-known Maxwell-Boltzmann distribution of velocities. Integrating over vx = —cc — oo shows that P(vx) is normalized. It is also easy to calculate the expectation value for the one-dimensional translational energy of a mole of gas as... 

Exercise. Let X stand for the three components of the velocity of a molecule in a gas. Give its range and probability distribution. 

Exercise. A cannon projects a ball with initial velocity v at angle 9 with the horizontal. Both v and 9 are subject to uncertainties that can be described by a Gaussian distribution for each. The distributions are centered at v0 and 60, respectively, and so sharp that non-physical values, such as negative v or 9, may be ignored. What is the probability distribution of the distance covered by the cannon ball ... 

Suppose a series of observations of the same Brownian particle gives a sequence of positions Xl9X2,-..- Each displacement Xk + 1 — Xk is subject to chance, but its probability distribution does not depend on the previous history, i.e., it is independent of Xk-l9 Xk 2.Hence, not only the velocity... 

For any r-component stochastic process one may ignore a number of components and the remaining s components again constitute a stochastic process. But, if the r-component process is Markovian, the process formed by the sfirst example above each velocity component is itself Markovian in chemical reactions, however, the future probability distribution of the amount of each chemical component is determined by the present amounts of all components. 

We shall call this a quasilinear Fokker-Planck equation, to indicate that it has the form (1.1) with constant B but nonlinear It is clear that this equation can only be correct if F(X) varies so slowly that it is practically constant over a distance in which the velocity is damped. On the other hand, the Rayleigh equation (4.6) involves only the velocity and cannot accommodate a spatial inhomogeneity. It is therefore necessary, if F does not vary sufficiently slowly for (7.1) to hold, to describe the particle by the joint probability distribution P(X, V, t). We construct the bivariate Fokker-Planck equation for it. 

The first complication, already pointed out in IX.7, is that the process X(t) is not uniquely determined by an initial value X(t0 but one needs to know V(t0) as well, either its value or its probability distribution. The second complication is that the first passage through R may occur with different values of V and it may happen, depending on the physical context, that these different events cannot be counted indiscriminately. A simple example is a particle diffusing in the presence of a potential wall, over which it escapes when it reaches it with a sufficient velocity. 

So far we only considered transport of particles by diffusion. As mentioned in 1 the continuous description was not strictly necessary, because diffusion can be described as jumps between cells and therefore incorporated in the multivariate master equation. Now consider particles that move freely and should therefore be described by their velocity v as well as by their position r. The cells A are six-dimensional cells in the one-particle phase space. As long as no reaction occurs v is constant but r changes continuously. As a result the probability distribution varies in a way which cannot be described as a succession of jumps but only in terms of a differential operator. Hence the continuous description is indispensable, but the method of compounding moments can again be used. 

---