Trajectory computation particle trajectories

If we wish to know the number of (VpV)-collisions that actually take place in this small time interval, we need to know exactly where each particle is located and then follow the motion of all the particles from time tto time t+ bt. In fact, this is what is done in computer simulated molecular dynamics. We wish to avoid this exact specification of the particle trajectories, and instead carry out a plausible argument for the computation of r To do this, Boltzmann made the following assumption, called the Stosszahlansatz, which we encountered already in the calculation of the mean free path ... 

Computational issues that are pertinent in MD simulations are time complexity of the force calculations and the accuracy of the particle trajectories including other necessary quantitative measures. These two issues overwhelm computational scientists in several ways. MD simulations are done for long time periods and since numerical integration techniques involve discretization errors and stability restrictions which when not put in check, may corrupt the numerical solutions in such a way that they do not have any meaning and therefore, no useful inferences can be drawn from them. Different strategies such as globally stable numerical integrators and multiple time steps implementations have been used in this respect (see [27, 31]). 

Simulations. In addition to analytical approaches to describe ion—soHd interactions two different types of computer simulations are used Monte Cado (MC) and molecular dynamics (MD). The Monte Cado method rehes on a binary coUision model and molecular dynamics solves the many-body problem of Newtonian mechanics for many interacting particles. As the name Monte Cado suggests, the results require averaging over many simulated particle trajectories. A review of the computer simulation of ion—soUd interactions has been provided (43). 

By integrating Eq. (13.35) step by step in time, the particle trajectory of the particle may be obtained. In the integration, the interaction between the particle and the wall may be approximated as being fully elastic however, when the particle hits the sidewall of the cyclone, the particle may be treated as being collected and the computation for the particle may terminated in order to save the computational time that may be required to track the particle to the bottom of the cyclone. If the particle trajectories for a range of particle diameters at different rates of fluid flow through the cyclone are determined, then the particle efficiency curve and the cut-off particle diameter of the cyclone may be obtained. 

Particle trajectories can be calculated by utilizing the modern CFD (computational fluid dynamics) methods. In these calculations, the flow field is determined with numerical means, and particle motion is modeled by combining a deterministic component with a stochastic component caused by the air turbulence. This technique is probably an effective means for solving particle collection in complicated cleaning systems. Computers and computational techniques are being developed at a fast pace, and one can expect that practical computer programs for solving particle collection in electrostatic precipitators will become available in the future. 

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. 

To compute 4>j, x, f jc, 0), particle trajectories must be obtained. In the numerical approaches, these particle trajectories are generated by following fluid particles in the numerical flow field. In the Monte Carlo method, an algorithm is formulated to produce the particle velocity and position as a function of time. Perhaps the simplest such algorithm is the following ... 

In the first case (Figure 8a), the side walls are adiabatic, and the reactor height (2 cm) is low enough to make natural convection unimportant. The fluid-particle trajectories are not perturbed, except for the gas expansion at the beginning of the reactor that is caused by the thermal expansion of the cold gas upon approaching the hot susceptor. On the basis of the mean temperature, the effective Rayleigh number, Rat, is 596, which is less than the Rayleigh number of 1844 necessary for the existence of a two-dimensional, stable, steady-state solution with flow in the transverse direction that was computed for equivalent Boussinesq conditions (188). 

The integration of forces between all molecules over several thousand time-steps produces particle trajectories from which time-averaged macroscopic properties can be computed. In MD the simulation is limited by the computer storage capacity and speed, so that short-lived phenomena (100-1000 ps) are generally calculated. 

Evident progress in studies of liquids has been achieved up to now with the use of computer simulations and of the models based on analytical theory. These methods provide different information and are mutually complementary. The first method employs rather rigorous potential functions and yields usually a chaotic picture of the multiple-particle trajectories but has not been able to give, as far as we know, a satisfactory description of the wideband spectra. The analytical theory is based on a phenomenological consideration (which possibly gives more regular trajectories of the particles than arise in reality ) in terms of a potential well. It can be tractable only if the profile of such a well is rather... 

P 62] A Lagrangian particle tracking technique, i.e. the computation of trajectories of massless tracer particles, which allows the computation of interfacial stretching factors, was coupled to CFD simulation [47]. Some calculations concerning the residence time distribution were also performed. A constant, uniform velocity and pressure were applied at the inlet and outlet, respectively. The existence of a fully developed flow without any noticeable effect of the inlet and outlet boundaries was assured by inspection of the computed flow fields obtained in the third mixer segment for all Reynolds numbers under study. 

The GCE model is more detailed than the TCE and VFD chemistry models. However, only through the availability of detailed information about the dependence of reaction cross sections on the precollision states of the colliding particles can the parameters g, P, and 7 be determined. In the results section, such information for the exchange reaction involving formation of nitric oxide is available from detailed trajectory computations and is used to determine suitable parameters for the GCE model for this reaction. 

With a Eulerian-Lagrangian approach, processes occurring at the particle surface can be modeled when simulating particle trajectories (for example, the process of dissolution or evaporation can be simulated). However, as the volume fraction of dispersed phase increases, the Eulerian-Lagrangian approach becomes increasingly computation intensive. A Eulerian-Eulerian approach more efficiently simulate such dispersed multiphase flows. 

It must be noted here that most industrial fluidized bed reactors operate in a turbulent flow regime. Trajectory simulations of individual particles in a turbulent field may become quite complicated and time consuming. Details of models used to account for the influence of turbulence on particle trajectories are discussed in Chapter 4. These complications and constraints on available computational resources may restrict the number of particles considered in DPM simulations. Eulerian-Eulerian approaches based on the kinetic theory of granular flows may be more suitable to model such cases. Application of this approach to simulations of fluidized beds is discussed below. 

MD allows the study of the time evolution of an V-body system of interacting particles. The approach is based on a deterministic model of nature, and the behavior of a system can be computed if we know the initial conditions and the forces of interaction. For a detail description see Refs. [14,15]. One first constructs a model for the interaction of the particles in the system, then computes the trajectories of those particles and finally analyzes those trajectories to obtain observable quantities. A very simple method to implement, in principle, its foundations reside on a number of branches of physics classical nonlinear dynamics, statistical mechanics, sampling theory, conservation principles, and solid state physics. 

Figure 8.1. An example of particle trajectory crossing in a 2D domain. The fully developed particle number density moo is shown together with the mean particle velocity vector (mio/moo, moi/moo). This example was computed using QBMM with a KBFVM (Yuan Fox, 2011) and realizable second-order fluxes (Vikas et al, 2011a). Note that the mean particle velocity is not defined in regions where the density is exactly zero. Owing to a smali amount of numerical diffusion, the density is very small, but nonzero, outside of the two jets.

Particle trajectories are then computed for both the small and the large-scale models in order to computing the flow resistances. For this purpose, the programa Sdtransp was used (Hendricks, 2001). Figure 5 shows the particle trajectories in a small-scale realization and in the associated upscaled model. The shape of the trajectories is very similar in both models the strong conductivity contrast between the two formations produces the observed refraction the lower formation acts almost as an impermeable boundary to flow from/to the upper formation. 

Two-way coupling allows for interaction between both phases by including the effects of the particulate phase on the fluid phase. In order to simplify the model and computation, the particles are usually assumed to be fully dispersed, i.e., they are not interacting with each other. The particle trajectory is updated in fixed intervals (so-called length scales) along the particle path. Additionally, the particle trajectory is updated each time the particle enters a neighboring cell. 

A list of typical dependent and independent variables for a furnace simulation is shown in Table VI. Coal particles typically are treated in a procedure that couples the continuum (Eulerian or fixed reference fame) gas-phase grid and the discrete (Lagrangian or moving reference frame) particles. Numerical solutions are repeated until the gas flow field is converged for the computed particle source terms, radiative fluxes, and gaseous reactions. Lagrangian particle trajectories are then calculated. After solving all particle-class trajectories, the new source terms... 

---