Validation of Lagrangian Two-Way Coupled Point-Particle Models in Large-Eddy Simulations

The accuracy of point-particle models with two-way coupling for particles of Kolmogorov-length-scale size is assessed. Turbulent kinetic energy budgets are analyzed in physical and in spectral space. It is shown that the force projection of the two-way coupling consistently models the direct transfer of kinetic energy on the particle surfaces and the enhanced viscous dissipation in the vicinity of the particles. Direct and large-eddy simulations of particle-laden flows in isotropic decaying turbulence are conducted and compared with direct-particle fluid simulations, where the particle-fluid interaction is fully resolved. An analysis in spectral space shows that turbulence modulation by particles mainly occurs at larger scales, although the momentum transfer takes place at the smallest scales. Therefore, the turbulent kinetic energy cascade of the single phase dominates in particle-laden flows. It is shown that point-particle models do not interfere with subgrid scale models, which usually act on the smallest scale. Consequently, point-particle models predict sufficiently accurate the turbulence modulation in direct numerical simulations and even when a subgrid scale model is used. The resolution of the LES does not affect the accuracy of the point-particle model, when the subgrid kinetic energy is negligible.     

Characteristics of turbulence transport for momentum and heat in particle-laden turbulent vertical channel flows

The dynamic and thermal performance of particle-laden turbulent flow is investigated via direction numerical simulation combined with the Lagrangian point-particle tracking under the condition of two-way coupling, with a focus on the contributions of particle feedback effect to momentum and heat transfer of turbulence. We take into account the effects of particles on flow drag and Nusselt number and explore the possibility of drag reduction in con-junction with heat transfer enhancement in particle-laden turbulent flows.The effects of particles on momentum and heat transfer are analyzed,and the possibility of drag reduc-tion in conjunction with heat transfer enhancement for the prototypical case of particle-laden turbulent channel flows is addressed.We present results of turbulence modification and heat transfer in turbulent particle-laden channel flow,which shows the heat transfer reduction when large inertial parti-cles with low specific heat capacity are added to the flow. However,we also found an enhancement of the heat transfer and a small reduction of the flow drag when particles with high specific heat capacity are involved.The present results show that particles,which are active agents,interact not only with the velocity field,but also the temperature field and can cause a dissimilarity in momentum and heat transport.This demonstrates that the possibility to increase heat transfer and suppress friction drag can be achieved with addition of par-ticles with different thermal properties.     

Point-Particle DNS and LES of Particle-Laden Turbulent flow - a state-of-the-art review

Particle-laden or droplet-laden turbulent flows occur in many industrial applications and in natural phenomena. Knowledge about the properties of these flows can help to improve the design of unit operations in industry and to predict for instance the occurrence of rain showers. This knowledge can be obtained from experimental research and from numerical simulations. In this paper a review is given of numerical simulation methods for particle-laden flows. There are various simulation methods possible. They range from methods in which all details, including the flow around each particle, are resolved, via point-particle methods, in which for each particle an equation of motion is solved, to Eulerian methods in which equations for particle concentration and velocity are solved. This review puts the emphasis on the intermediate class of methods, the Euler-Lagrange methods in which the continuous phase is described by an Eulerian approach and the dispersed phase in a Lagrangian way with equations of motion for each individual particle.     

Direct numerical simulation of a particle-laden low Reynolds number turbulent round jet

Direct numerical simulation method is used for the investigating of particle-laden turbulent flows in a spatially evolution of low Reynolds number axisymmetric jet, and the Eulerian–Lagrangian point-particle approach is employed in the simulation. The simulation uses an explicit coupling scheme between particles and the fluid, which considers two-way coupling between the particle and the fluid. The DNS results are compared well with experimental data with equal Reynolds number (R_e = 1700). Our objects are: (i) to investigate the correlation between the particle number density and the fluctuating of fluid streamwise velocity; (ii) to examine whether the three-dimensional vortex structures in the particle-laden jet are the same as that in the free-air jet and how the particles modulate the thee-dimensional vortex structures and turbulence properties with different Stokes number particles; (iii) to discover the particle circumferential dispersion with different Stokes number particles. Our findings: (i) all the particles, regardless of their particle size, tend to preferentially accumulate in the region with large-than-mean fluid streamwise velocity; (ii) the small Stokes number particles take an important part in the modulation of three-dimensional vortex structures, but for the intermediate and larger sized particles, this modulation effect seems not so apparent; (iii) the particle circumferential dispersion is more effective for the smaller and intermediate sized particles, especially for the intermediate sized particles.     

Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion

A numerical study based on the Eulerian–Lagrangian formulation is performed for dispersed phase motion in a turbulent flow. The effect of spatial filtering, commonly employed in large-eddy simulations, and the role of the subgrid scale turbulence on the statistics of heavy particles, including preferential concentration, are studied through a priori analysis of DNS of particle-laden forced isotropic turbulence. In simulations where the subgrid scale kinetic energy attains 30–35% of the total we observe the impact of residual fluid motions on particles of a smaller inertia. It is shown that neglecting the influence of subgrid scale fluctuations has a significant effect on the preferential concentration of those particles. A stochastic Langevin model is proposed to reconstruct the residual (or subgrid scale) fluid velocity along particle trajectories. The computation results for a selection of particle inertia parameters are performed to appraise the model through comparisons of particle turbulent kinetic energy and the statistics of preferential concentrations.     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 1. Isothermal flows

This paper presents a review of authors' collective works in the field of two-phase flow modeling done in the past few decades. The paper is aimed at the construction of mathematical models for simulation of particle-laden turbulent flows. A kinetic equation was obtained for the probability density function (PDF) of the particle velocity distribution in turbulent flows. The proposed kinetic equation describes both the interaction of particles with turbulent eddies of the carrier phase and particle-particle collisions. This PDF equation is used for the derivation of different schemes describing turbulent momentum transfer in the dispersed particle phase. The turbulent characteristics of the gaseous phase are calculated on the basis of the k - turbulence model with a modulation effect of particles on the turbulence.

The constructed models have been applied to the calculation of various two-phase gas-particle turbulent flows in jets and channels as well as particle deposition in tubes and separators. For validating the theoretical and numerical results, a wide range of comparisons with experimental data from Russian and foreign sources has been done.     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 2. Flows with heat transfer and combustion

The objective of this part of the paper is to summarize the information concerning the authors' works in the field of simulation of two-phase gas-particle turbulent flows with heat transfer and combustion. A kinetic equation had been derived for the probability density function (PDF) of the particle velocity, temperature, and mass distributions in turbulent flows. This PDF equation is used for the construction of the governing conservation equations of mass, momentum, and heat transfer in the dispersed particle phase.The numerical scheme incorporates two-phase fluid dynamics, convective and radiative heat transfer, and combustion. The proposed models have been applied to the calculation of various particle-laden turbulent flows in jets, combustion and gasification chambers, and furnaces.     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

A sharp interface Cartesian grid method for viscous simulation of shocked particle-laden flows

A Cartesian grid-based sharp interface method is presented for viscous simulations of shocked particle-laden flows. The moving solid–fluid interfaces are represented using level sets. A moving least-squares reconstruction is developed to apply the no-slip boundary condition at solid–fluid interfaces and to supply viscous stresses to the fluid. The algorithms developed in this paper are benchmarked against similarity solutions for the boundary layer over a fixed flat plate and against numerical solutions for moving interface problems such as shock-induced lift-off of a cylinder in a channel. The framework is extended to 3D and applied to calculate low Reynolds number steady supersonic flow over a sphere. Viscous simulation of the interaction of a particle cloud with an incident planar shock is demonstrated; the average drag on the particles and the vorticity field in the cloud are compared to the inviscid case to elucidate the effects of viscosity on momentum transfer between the particle and fluid phases. The methods developed will be useful for obtaining accurate momentum and heat transfer closure models for macro-scale shocked particulate flow applications such as blast waves and dust explosions.     

Modeling anisotropic Reynolds-stress dissipation in particle- or droplet-laden flows

Particles or droplets dispersed in turbulent flows at sufficiently high volume loadings lead to modifications of turbulence characteristics. More specifically, in a detailed experimental investigation by Poelma and coworkers, where the particle phase moves with a non-zero mean velocity relative to the fluid phase, it was found that anisotropic Reynolds-stress dissipation is induced. Recently, we have proposed a model that can account for this effect in RANS-based and PDF method simulations. In our previous work, however, no simulation results of the RANS model formulation were presented. In the present work, a new compact tensorial RANS formulation is presented and the new formulation is validated against the experimental data of Poelma and coworkers.     

Suspended load and bed-load transport of particle-laden gravity currents: the role of particle–bed interaction

The development of particle-enriched regions (bed-load) at the base of particle-laden gravity currents has been widely observed, yet the controls and relative partitioning of material into the bed-load is poorly understood. We examine particle-laden gravity currents whose initial mixture (particle and fluid) density is greater than the ambient fluid, but whose interstitial fluid density is less than the ambient fluid (such as occurs in pyroclastic flows produced during volcanic eruptions or when sediment-enriched river discharge enters the ocean, generating hyperpycnal turbidity currents). A multifluid numerical approach is employed to assess suspended load and bed-load transport in particle-laden gravity currents under varying boundary conditions. Particle-laden flows that traverse denser fluid (such as pyroclastic flows crossing water) have leaky boundaries that provide the conceptual framework to study suspended load in isolation from bed-load transport. We develop leaky and saltation boundary conditions to study the influence of flow substrate on the development of bed-load. Flows with saltating boundaries develop particle–enriched basal layers (bed-load) where momentum transfer is primarily a result of particle–particle collisions. The grain size distribution is more homogeneous in the bed-load and the saltation boundaries increase the run-out distance and residence time of particles in the flow by as much as 25% over leaky boundary conditions. Transport over a leaky substrate removes particles that reach the bottom boundary and only the suspended load remains. Particle transport to the boundary is proportional to the settling velocity of particles, and flow dilution results in shear and buoyancy instabilities at the upper interface of these flows. These instabilities entrain ambient fluid, and the continued dilution ultimately results in these currents becoming less dense than the ambient fluid. A unifying concept is energy dissipation due to particle–boundary interaction: leaky boundaries dissipate energy more efficiently at the boundary than their saltating counterparts and have smaller run-out distance.

---

     

Efficient simulation of particle-laden turbulent flows with high mass loadings using LES

The paper is concerned with the simulation of particle-laden two-phase flows based on the Euler–Lagrange approach. The methodology developed is driven by two major requirements: (i) the necessity to tackle complex turbulent flows by eddy-resolving schemes such as large-eddy simulation; (ii) the demand to predict dispersed multiphase flows at high mass loadings. First, a highly efficient particle tracking algorithm was developed working on curvilinear, block-structured grids. Second, to allow the prediction of dense two-phase flows, the fluid–particle interaction (two-way coupling) as well as particle–particle collisions (four-way coupling) had to be taken into account. For the latter instead of a stochastic collision model, in the present study a deterministic collision model is considered. Nevertheless, the computational burden is minor owing to the concept of virtual cells, where only adjacent particles are taken into account in the search for potential collision partners. The methodology is applied to different test cases (plane channel flow, combustion chamber flow). The computational results are compared with experimental measurements and good agreement is found.     

非球形颗粒两相流的数值模拟研究进展

非球形颗粒两相流是多相流的重要研究方向之一, 常见于自然界及工业生产过程中. 不同于球形颗粒, 由于非球形颗粒形状的各向异性, 除了颗粒平动行为, 还需要考虑颗粒的转动与取向行为, 颗粒的取向与转动行为会影响颗粒所受的力和力矩. 为了准确模拟非球形颗粒的运动行为, 目前非球形颗粒两相流的数值模拟研究主要基于欧拉?拉格朗日的求解框架展开, 常见的非球形颗粒两相流数值模拟方法主要包括点颗粒法与全分辨颗粒法. 本文将对这两类方法进行介绍, 同时会全面介绍非球形颗粒两相流研究的基础理论模型, 并系统总结非球形颗粒在简单基本流和复杂湍流中的研究进展, 包括对于非球形颗粒在湍流中的取向与转动行为机理, 以及颗粒对湍流减阻调制作用的研究. 最后, 本文提出了非球形颗粒两相流研究存在的问题及未来研究方向.      

Steady-state particle distributions in shear flows of colloids and finely dispersed suspensions

On the basis of an analysis of the pseudoturbulent motion of both the suspended particles and the carrier fluid, the normal stress components in the dispersed phase are obtained for the problem of inclined confined flows of finely dispersed suspensions and colloids. These hydrodynamic pulsations are due to the shear and the work done by the average relative flow of the fluid phase on random concentration fluctuations of the disperse system because of the substantial slip of the phases of the suspension under gravity. The momentum conservation equations for the particles are obtained with allowance for the angle of inclination of the flow to the vertical and on the basis of these equations the suspension capacity of the flow as a function of the angle of inclination, particle size, Galileo number and other parameters is illustrated.Ekaterinburg. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 78–84, January–February, 1996.     

Effect of inertial particles with different specific heat capacities on heat transfer in particle-laden turbulent flow

The effect of inertial particles with different specific heat on heat transfer in particle-laden turbulent channel flows is studied using the direct numerical simulation(DNS) and the Lagrangian particle tracking method. The simulation uses a two-way coupling model to consider the momentum and thermal interactions between the particles and turbulence. The study shows that the temperature fields display differences between the particle-laden flow with different specific heat particles and the particle-free flow,indicating that the particle specific heat is an important factor that affects the heat transfer process in a particle-laden flow. It is found that the heat transfer capacity of the particle-laden flow gradually increases with the increase of the particle specific heat. This is due to the positive contribution of the particle increase to the heat transfer. In addition,the Nusselt number of a particle-laden flow is compared with that of a particle-free flow.It is found that particles with a large specific heat strengthen heat transfer of turbulent flow, while those with small specific heat weaken heat transfer of turbulent flow.     

Large eddy simulation of gas-particle turbulent channel flow with momentum exchange between the phases

This paper presents results of a large eddy simulation (LES) combined with Lagrangian particle tracking and a point-force approximation for the feedback effect of particles on the downward turbulent gaseous flow in a vertical channel. The LES predictions are compared with the results obtained by direct numerical simulation (DNS) of a finer computational mesh. A parametric study is conducted for particles with two response times in simulations with and without streamwise gravitational settling and elastic, binary interparticle collisions. It is shown that the classical and the dynamic Smagorinsky turbulence models adequately predict the particle-induced changes in the mean streamwise velocity and the Reynolds stresses of the carrier phase for the range of parameters studied. However, the largest discrepancies between the LES and DNS results are found in the cases of particle-laden flows. Conditional sampling of the instantaneous resolved flow fields indicates that the mechanisms by which particles directly oppose the production of momentum and vorticity of the organized fluid motions are also observed in the LES results. However, the geometric features of the near-wall quasistreamwise vortices are overestimated by the use of both turbulence models compared to the DNS predictions.     

An analytical approach for the closure equations of gas–solid flows with inter-particle collisions

This work examines the effect of inter-particle collisions on the motion of solid particles in two-phase turbulent pipe and channel flows. Two mechanisms for the particle–particle collisions are considered, with and without friction sliding. Based on these collision mechanisms, the correlations of the various velocity components of colliding particles are obtained analytically by using an averaging procedure. This takes into account three collision coordinates, two angles and the distance between the centers of colliding particles. The various stress tensor components are obtained and then introduced in the mass, linear momentum and angular momentum equations of the dispersed phase. The current approach applies to particle–particle collisions that result from both the average velocity difference and the turbulent velocity fluctuations. In order to close the governing equations of the dispersed phase, the pseudo-viscosity coefficients are defined and determined by the time of duration of the inter-particle collision process. The model is general enough to apply to both polydisperse and monodisperse particulate systems and has been validated by comparisons with experimental data.