Dynamics of fibres in a turbulent flow field – A particle-level simulation technique

A particle-level simulation technique has been developed for modelling the flow of fibres in a turbulent flow field. A single fibre is conceived here as a chain of segments, thus enabling the model fibre to have all the degrees of freedom (translation, rotation, bending and twisting) needed to realistically reproduce the dynamics of real fibres. Equations of motion are solved for each segment, accounting for the interaction forces with the fluid, the contact forces with other fibres and the forces that maintain integrity of the fibre.The motion of the fluid is resolved as a combination of 3D mean flow velocities obtained from a CFD code and fluctuating turbulent velocities derived from the Langevin equation. A case of homogeneous turbulence is treated in this paper.The results obtained show that fibre flocs in air-fibre flows can be created even when attractive forces are not present. In such a case, contacts between fibres, properties of an individual fibre (such as flexibility and equilibrium shapes) and properties of the flow of the carrying fluid are shown to govern the physics behind formation and breaking up of fibre flocs. Highly irregular fibre shapes and stiff fibres lead to strong flocculation.The modelling framework applied in this work aims at making possible a numerical model applicable for designing processes involving transport of fibres by air at industrial scale.     

Liquid/liquid dispersion in a chaotic advection flow

Mixing by chaotic advection in a twisted-pipe flow is used here to investigate the efficiency of this flow in the liquid/liquid dispersion process. This study focuses on water/oil dispersions produced by continuous water injection into a main oil flow, for small Dean numbers. The drop sizes obtained with the chaotic-advection twisted-pipe flow are compared with those in a straight pipe and a helically coiled flow for the same conditions. It is found that the resulting dispersions are finer and more mono-dispersed in the chaotic advection flow. These results are compared with the theoretical maximum diameter _dmax$d_{\text{max}}$ determined by the Grace theory in which the viscous stress controls the breakup phenomena. For this purpose, the kinematic field is computed from the theoretical formulae for Dean flow. The strain rate fields in the pipe cross-section are then analytically computed and used to predict the maximum drop diameter. The theoretical values are identical for the three configurations (straight, helically coiled, and twisted pipe) up to a critical Dean number, where the secondary flow becomes significant. Beyond this value, the shear stress is enhanced in the twisted-pipe flow compared with the straight-pipe flow, and the predicted drop diameters are smaller. An interpretation of the higher dispersive performance of the chaotic flow is provided by the Lagrangian trajectories of the particles.     

A closed model of fluctuating particle motion in a turbulent flow

Random particle motion in a turbulent and molecular velocity fluctuation field is considered. Using a spectral representation of the carrier-phase Eulerian velocity fluctuation correlations, a closed system of integral equations for calculating the carrier-phase velocity correlation along the particle trajectory and the particle Lagrangian velocity fluctuation correlation is obtained. Based on this system, the fluctuations of the particle parameters are analyzed. In the limiting case of a passive admixture, an estimate is found for the ratio of the integral Lagrangian and Eulerian time scales and the Kolmogorov constant for the Lagrangian structure function of the carrier-phase velocity fluctuations.     

A dissipative particle dynamics study of flow in periodically grooved nanochannels

The effect of periodic rectangular wall roughness on planar nanochannel flow is investigated by dissipative particle dynamics simulation. The wall protrusion length is varied, and its effect on the flow is examined. Analysis of particle trajectories and average residence time reveals temporary trapping of fluid particles inside the rectangular cavities for a considerable amount of time. This trapping affects the density, velocity, pressure, and temperature distribution inside and close to the cavities. Inside the cavities, low‐velocity regions and regions of high density related to high pressure and high temperature are observed. When compared with that of the channel with flat walls case, lower flow velocities, temperatures, and pressures are observed for grooved channels. The reduction of the above quantities is more pronounced as the protrusion length, that is, the roughness characteristic length, decreases. Finally, the relation of friction factor, f, with the flow Reynolds number is discussed. The model predicts = constant in the range . The results of this work are of direct relevance to the design of nanofluidic devices. Copyright © 2011 John Wiley & Sons, Ltd.     

Consistency issues of Lagrangian particle tracking applied to a spray jet in crossflow

Numerical simulations are performed for multiphase jets in crossflow. The flow solver uses an Eulerian/Lagrangian approach. Turbulence in the gas phase is modeled in the framework of large eddy simulation. The dispersed phase is handled using Lagrangian particle tracking. The model assumptions of solvers for Lagrangian particle tracking are critically assessed for typical flow conditions of spray jets in crossflow. The droplets are assumed to be spherical and isolated. It is shown that several model assumptions are apparently inconsistent in larger portions of the flow field. Firstly, average Weber numbers can be so large that the model assumption to regard droplets as spherical is questionable, not only near the nozzle, but also in the far-field. Secondly, the average droplet spacing can be so low that droplets directly interact with each other, again also in the far-field. Thirdly, the average Stokes numbers in the jet region can be so large that the phase coupling between the dispersed and continuous phase is weak. Some remedies to these deficiencies are proposed.     

Dynamic simulation of arbitrarily shaped particles in shear flow

A new simulation framework was created for modeling the dynamics of arbitrarily shaped particles dispersed in Newtonian fluid. Theoretical complexity usually restricts suspension simulations to those for spheroids. This new simulation is loosely based on the Stokesian Dynamics method including long range hydrodynamic interaction and uses spheres as building components for greater particulates of arbitrary shape. This approach is capable of accurately reproducing the dynamics of an isolated arbitrarily shaped particle. As verification, the simulated results are compared against known results for a rod-like particle. An elongated rod-shaped structure made from linked spheres is shown to reproduce the well-known elongated ellipsoidal particle dynamics described by [Jeffrey Proc R Soc Lond A 102:161–179, 1923]. The predicted orbital period and spin rates for a fiber in shear are reproduced and compare well with theoretical prediction over a wide aspect ratio range. Predicted particle dynamics for other shaped particles are then demonstrated.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Forces on a spherical particle in shear flow at intermediate Reynolds numbers

When particles are submerged in a shear flow, there are lateral (lift) forces on the particles, and these lateral forces affect the dispersion of the particles very much. Recent literature survey indicates that there are large discrepancies among the results from the previous numerical investigations on this subject. A small computational domain ranging between 20–30 sphere radii was used in all the previous numerical investigations. However, the result from the present study reveals that the value of lift coefficient strongly depends on the size of computational domain. To provide correct numerical data and physical interpretation for the forces on a spherical particle in linear shear flow, accurate numerical computations were performed for 5≤Re≤200 using a computational domain of 101 sphere radii.     

Physics and modelling of turbulent particle deposition and entrainment: Review of a systematic study

Deposition and entrainment of particles in turbulent flows are crucial in a number of technological applications and environmental processes. We present a review of recent results from our previous works, which led to physical insights on these phenomena. These results were obtained from a systematic numerical study based on the accurate resolution – Direct Numerical Simulation via a pseudo-spectral approach – of the turbulent flow field, and on Lagrangian tracking of particles under different modelling assumptions. We underline the multiscale aspect of wall turbulence, which has challenged scientists to devise simple theoretical models adequate to fit experimental data, and we show that a sound rendering of wall turbulence mechanisms is required to produce a physical understanding of particle deposition and re-entrainment. This physical understanding can be implemented in more applied simulation techniques, such as Large-Eddy Simulation. Our arguments are based also on the phenomenology of coherent structures and on the examination of flow topology in connection with particle preferential distribution. Starting from these concepts, reasons why theoretical predictions may fail are examined together with the requirements which must be fulfilled by suitable predictive models.     

Influence of gravity and lift on particle velocity statistics and transfer rates in turbulent vertical channel flow

The present study reports detailed statistics for velocity and transfer rates of heavy particles dispersed in turbulent boundary layers. Statistics have been extracted from a homogeneous source of data covering a large target parameter space and are used here to analyze the effects of gravity and lift on particle dispersion and deposition in a systematic way. Datasets were obtained performing Direct Numerical Simulation (DNS) of particle-laden turbulent upward/downward flow in a vertical channel. Six values for the particle timescale (the particle Stokes number, St) ranging three orders of magnitude were considered to analyze the deposition process as the controlling mechanism was shifting from turbulent diffusion to inertia-moderated crossing trajectories. For the particle timescales examined, gravity and lift do not influence the qualitative behavior of particles even though velocity profiles and deposition coefficients are modified in a non-monotonic fashion, reaching an optimum for St ? 15. Physical mechanisms for the different behavior are discussed. Raw data and statistics obtained from the present DNS are made available at http://cfd.cineca.it (mirror site: http://158.110.32.35/download/database) and can be used to test simple models and closure equations for multiphase RANS and Large Eddy simulations.     

考虑颗粒间碰撞的气固两相流拉格朗日模拟

在均匀,稳定的各向同性气固两相紊流场颗粒弥散的拉格朗日模拟计算方法基础上,进一步考虑了流场中颗粒之间的碰撞对于模拟计算结果的影响。与Lavieville用大涡模拟所做的计算结果进行了对比,以对本方法进行验证,并考察了颗粒间的碰撞分别对流体相和颗粒相的影响。     

Transmission fluctuation method for particle analysis in multiphase flow

The transmission fluctuations can be used for measurements on particle size and particle concentration in multiphase flows. In Gregory’s method, the beam diameter is required to be much larger than the particle diameter, which limits the application seriously. In this work, we introduce a new method of the transmission fluctuation measurement, in which the beam diameter can be less, equal to and larger than the particle diameter. The theoretical analysis proves that, when the beam-to-particle diameter ratio is within the range of 0.05 and 10, the new method is able to achieve satisfying measurement results. However, Gregory’s method is only suitable in the range of 5 and 10. Therefore, this new method enhances the transmission fluctuation measurements greatly.     

Chaotic advection in the near-wake of a low-rise building

Properties of the mean-streamlines of a turbulent wake flow behind a low-rise building have been investigated in numerical simulations. For some near wind directions, the solutions to the equations for the mean-streamlines in the recirculating region of flow are found to be chaotic. This is shown to be analogous to the occurrence of chaotic advection in laminar flows and consistent with the current understanding of chaotic dynamics in Hamiltonian systems. The transitions between regular and chaotic solutions for the mean-streamlines are found to occur abruptly over incremental changes, Δθ<2°, in the mean wind direction, θ. These transitions are shown to be accompanied by abrupt changes in the characteristics of plumes formed from tracer-particles undergoing turbulent dispersion. Small fluctuations in the mean wind direction at the transition between regular and chaotic solutions for the mean-streamlines are predicted to result in strongly intermittent turbulent dispersion.     

Numerical study of short-circuiting flow and particles in a gas cyclone

Short-circuiting flow is an important secondary flow in gas cyclones, which has a negative impact on the separation performance. To improve the understanding of the short-circuiting flow and guide the optimization of gas cyclones, this paper presents a numerical study of a cyclone using computational fluid dynamics. Based on the steady flow field, three methods were adopted to investigate the formation mechanism and characteristics of the short-circuiting flow and particles. The temporal variation of the tracer species concentration distribution reveals that the formation mechanism of the short-circuiting flow is the squeeze between the airflows entering the annular space of the gas cyclone at different times. The short-circuiting flow region, distinguished through the spatial distribution of the moments of age, is characterized by a small mean age and a large coefficient of variation. The proportion of the short-circuiting particles increases with the increase of the inlet velocity only for small particles. But with the increase of particle size, the proportion of the short-circuiting particles decreases faster at higher inlet velocities, resulting in significant differences in collection efficiency curves.     

Particle and bubble dynamics in a creeping flow

Dynamics of a solid particle and non-deformable gaseous bubble in viscous fluid are studied analytically and numerically within the framework of creeping flow regime (flow at vanishingly small Reynolds numbers). Equations of motion for the particle and bubble include the consideration of the buoyancy force, Stokes drag force and memory-integral drag force. Exact analytical solutions are obtained and categorised in terms of inclusion (particle or bubble) density with respect to the density of a surrounding fluid. Through the analytical and numerical solutions, the dynamics of solid particle and air bubble in water have been found to behave differently especially at the early stages of motion, whereas some qualitative similarities exist in the long-term asymptotic.     

Averaging of particle data from phase Doppler anemometry in unsteady two-phase flow: Validation by numerical simulation

A novel post-processing algorithm is proposed to correct statistical bias observed in the treatment of time series obtained by a phase Doppler anemometer (PDA) at flow locations with variable particle velocity and concentration. Extensive properties of each validated particle are weighted with their inverse measuring (validation) volume to account for the procedure of particle sampling and fluctuations in the particle concentration. To compensate for the short characteristic length of the validation volume, the properties of particles are expressed by properties of fields of particle groups, using a local averaging time. A window shift and a decorrelation scheme are applied on the fields to increase their frequency resolution. This algorithm has been tested on numerical time series, provided by an Eulerian/Lagrangian code representing a gas/solids flow past a bluff body. Moments and spectral estimates of concentration and velocity of particle groups were successfully validated by the numerical simulation using the PDA data algorithm and control volume averaging. The control volume was much larger than the PDA validation volume, but the centre positions of the two volumes were identical.     

Estimation of particle dynamics in 2-D fluidized beds using particle tracking velocimetry

The experimental characterization of particle dynamics in fluidized beds is of great importance in fostering an understanding of solid phase motion and its effect on particle properties in granulation processes. Commonly used techniques such as particle image velocimetry rely on the cross-correlation of illumination intensity and averaging procedures. It is not possible to obtain single particle velocities with such techniques. Moreover, the estimated velocities may not accurately represent the local particle velocities in regions with high velocity gradients. Consequently, there is a need for devices and methods that are capable of acquiring individual particle velocities. This paper describes how particle tracking velocimetry can be adapted to dense particulate flows. The approach presented in this paper couples high-speed imaging with an innovative segmentation algorithm for particle detection, and employs the Voronoi method to solve the assignment problem usually encountered in densely seeded flows. Lagrangian particle tracks are obtained as primary information, and these serve as the basis for calculating sophisticated quantities such as the solid-phase flow field, granular temperature, and solid volume fraction. We show that the consistency of individual trajectories is sufficient to recognize collision events.     

Effect of heat-source geometry on distribution and deposition of particulates in a ventilated chamber

To investigate the effect of near-wall heat-source shape on particulate motion, the particulate distribution and deposition in a ventilated chamber with different heat-source configurations were numerically modeled. Using the discrete random walk model of the Lagrangian method, the trajectories of 3200 mono-disperse particulates ranging from 1 to 10 μm with a density of 1400 kg/m³ were tracked. Airflow pattern, temperature fields, the distribution of particulate concentrations, and deposition patterns are calculated and presented. The results show that the shape of a near-wall heat source has an influence on the airflow as well as the temperature field in the chamber and hence affects the particulate distribution and deposition.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured.The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured. The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.