Heat transfer effects on particle motion under rarefied conditions

Gas–solid flows occurring on very small spatial scales (of the order of micro and nanometres) are of great relevance in a number of industrial applications. It is currently not well established how particle motion and filtration are affected by non-isothermal conditions at such scales. Furthermore, when the particle size is comparable to the mean free path of the gas, rarefaction effects become important. In the present work we investigate the effects of heat transfer and non-isothermal conditions on the motion of small particles in rarefied flow. For that purpose, a suitable framework is developed here as a generic multiphase DNS method for rarefied flows. The resulting model is valid for low particle Reynolds number flows, irrespective of the Biot number, and for particle Knudsen numbers up to unity in unbounded flow. Using this model, we show that there is different settling behaviour of particles with an internal heat source in rarefied and continuum cases of the carrier gas respectively. It is shown that the chances for thermal levitation and/or lifting up of a particle due to buoyancy effects are significantly reduced under rarefied conditions.     

A scalable parallel algorithm for the direct numerical simulation of three-dimensional incompressible particulate flow

Particulate flow is of great importance from both the scientific and engineering points of view. Owing to the complexity of particle-flow interactions, direct numerical simulations (DNS) of inertial particulate flow with finite-size particles have been limited to a very small number of particles, while the industrial applications involve larger numbers with many orders of magnitude. This article presents a parallel implementation of a fictitious domain method for the DNS of particulate flows. The method is thoroughly tested and its parallel performance on distributed memory clusters is evaluated on a large-scale problem. Finally, we present the results for the separation of 21,336 particles of two different densities in a viscous fluid. Although there is still a significant gap between DNS and the industrial applications, the present algorithm allows to simulate significantly large number of particles so that a meaningful statistical analysis can be performed. This will help in the development of new closure relations for the averaged models of multiphase flows.     

A multiphase DNS approach for handling solid particles motion with heat transfer

In the current work we propose a multiphase DNS method capable of resolving the motion of solid particles coupled with heat transfer effects. The method is based on solving a shared set of momentum and energy balance equations for the carrier phase and the particulate phase. Individual particles are tracked using a number of volume fraction advection equations. The proposed method is in very good agreement with the available data in the literature for the following cases: isothermal particle motion (in the presence of walls and other particles), natural convection around a stationary particle and solid particles motion accompanied with heat transfer effects. In addition, we show that the method is inherently capable of handling deformable particles (i.e. droplets and bubbles) co-existing with solid particles. The method is thus well suited to deal with challenging multiphase systems, such as diesel spray combustion with soot formation, spray drying with particle nucleation, and biological treatment of waste water.     

多相流局部混合型质点网格法

提出模拟多相流的局部混合型质点网格法,该方法能稳定地模拟高浓度流体颗粒两相 流. 在每一个颗粒团尺度的欧拉网格下(本文称之为欧拉微元),基于Lagrangian 追踪原理,可直接估计体积内颗粒的总量,从而准确求出欧拉微元和控制容积内颗粒的 浓度(即颗粒在容积内的体积含量). 同时,假设在新的时间步下,颗粒在欧拉微元里充分 混合,形成新颗粒团. 作者对竖直和倾斜容器中单粒径颗粒沉降和竖直容器中双粒径颗粒双 峰悬浮液沉降过程进行了计算,结果与实测数据相符.     

湍流横掠固定颗粒的全尺度直接数值模拟

气固两相流模拟中,当固相尺度接近或大于Kolmogorov尺度时,普通的点源模型将不再适用,固体相的体积效应和表面效应将对流体相产生显著的影响。通过采用直接数值模拟方法,结合内嵌边界方法对湍流中不同湍流强度流体横掠大于Kolmogorov尺度的固相颗粒进行了全尺度模拟,讨论分析了在两种湍流度下方形颗粒对湍流的调制影响以及颗粒的受力情况。     

Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

多相流动的直接数值模拟进展

多相流动的直接数值模拟的讨论就是把流场中颗粒周围计算网格缩小到颗粒尺寸以下进行流动的计算,颗粒的受力不是通过模型计算,而是通过积分表面的黏性力与压力获得.直接模拟方法的出现标志着人们对多相流动的认识从宏观扩展到微观层次.主要介绍了几种先进的直接模拟的方法:基于体适应的非结构化移动网格方法;基于固定网格方法;其它方法.      

Effect of Particle Clusters on Carrier Flow Turbulence: A Direct Numerical Simulation Study

Experiments indicate that particle clusters that form in fluidized–bed risers can enhance gas-phase velocity fluctuations. Direct numerical simulations (DNS) of turbulent flow past uniform and clustered configurations of fixed particle assemblies at the same solid volume fraction are performed to gain insight into particle clustering effects on gas-phase turbulence, and to guide model development. The DNS approach is based on a discrete-time, direct-forcing immersed boundary method (IBM) that imposes no-slip and no-penetration boundary conditions on each particle’s surface. Results are reported for mean flow Reynolds number Re _p ?=?50 and the ratio of the particle diameter d _p to Kolmogorov scale is 5.5. The DNS confirm experimental observations that the clustered configurations enhance the level of fluid-phase turbulent kinetic energy (TKE) more than the uniform configurations, and this increase is found to arise from a lower dissipation rate in the clustered particle configuration. The simulations also reveal that the particle-fluid interaction results in significantly anisotropic fluid-phase turbulence, the source of which is traced to the anisotropic nature of the interphase TKE transfer and dissipation tensors. This study indicates that when particles are larger than the Kolmogorov scale (d _p ?>?η), modeling the fluid-phase TKE alone may not be adequate to capture the underlying physics in multiphase turbulence because the Reynolds stress is anisotropic. It also shows that multiphase turbulence models should consider the effect of particle clustering in the dissipation model.     

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.     

Comparative study of two corrective gradient models in the simulation of multiphase flows using moving particle semi-implicit method

In this study, two corrective gradient models (CGMs) are compared in the simulation of multiphase flows. Linear consistency of the gradient model of moving particle semi-implicit (MPS) method has been recovered by introducing corrective matrix. However, it is found that particles tend to disperse along the streamline while using the CGM proposed in a previous study. Particle shifting (PS) schemes are necessary to reduce the irregularity of particle distribution to stabilize the calculation. To enhance the accuracy and stability, another CGM with dummy particle (CGMD) was proposed in our previous study. This enhanced CGM is featured by linear consistency and purely repulsive pressure gradient force. In this study, this enhanced CGM is modified and applied to multiphase flow simulation. Comparative study suggests that the modified CGM with PS scheme is capable of calculating various multiphase flows and predicting the interface evolution both clearly and accurately.     

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.     

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.     

A scaling analysis for point–particle approaches to turbulent multiphase flows

Simple dimensional arguments are used in establishing three different regimes of particle time scale, where explicit expression for particle Reynolds number and Stokes number are obtained as a function of nondimensional particle size (d/η)$(d/\eta )$ and density ratio. From a comparative analysis of the different computational approaches available for turbulent multiphase flows it is argued that the point–particle approach is uniquely suited to address turbulent multiphase flows where the Stokes number, defined as the ratio of particle time scale to Kolmogorov time scale (_τp/_τk)$({\tau }_p/{\tau }_k)$, is greater than 1. The Stokes number estimate has been used to establish parameter range where point–particle approach is ideally suited. The point–particle approach can be extended to handle “finite-sized” particles whose diameter approach that of the smallest resolved eddies. However, new challenges arise in the implementation of Lagrangian–Eulerian coupling between the particles and the carrier phase. An approach where the inter-phase momentum and energy coupling can be separated into a deterministic and a stochastic contribution has been suggested.     

An improved SPH model for multiphase flows with large density ratios

This paper presents a new smoothed particle hydrodynamics (SPH) model for simulating multiphase fluid flows with large density ratios. The new SPH model consists of an improved discretization scheme, an enhanced multiphase interface treatment algorithm, and a coupled dynamic boundary treatment technique. The presented SPH discretization scheme is developed from Taylor series analysis with kernel normalization and kernel gradient correction and is then used to discretize the Navier‐Stokes equation to obtain improved SPH equations of motion for multiphase fluid flows. The multiphase interface treatment algorithm involves treating neighboring particles from different phases as virtual particles with specially updated density to maintain pressure consistency and a repulsive interface force between neighboring interface particles into the pressure gradient to keep sharp interface. The coupled dynamic boundary treatment technique includes a soft repulsive force between approaching fluid and solid particles while the information of virtual particles are approximated using the improved SPH discretization scheme. The presented SPH model is applied to 3 typical multiphase flow problems including dam breaking, Rayleigh‐Taylor instability, and air bubble rising in water. It is demonstrated that inherent multiphase flow physics can be well captured while the dynamic evolution of the complex multiphase interfaces is sharp with consistent pressure across the interfaces.     

Turbulent particle dispersion in arbitrary wall-bounded geometries: A coupled CFD-Langevin-equation based approach

A Lagrangian continuous random walk (CRW) model is developed to predict turbulent particle dispersion in arbitrary wall-bounded flows with prevailing anisotropic, inhomogeneous turbulence. The particle tracking model uses 3D mean flow data obtained from the Fluent CFD code, as well as Eulerian statistics of instantaneous quantities computed from DNS databases. The turbulent fluid velocities at the current time step are related to those of the previous time step through a Markov chain based on the normalized Langevin equation which takes into account turbulence inhomogeneities. The model includes a drift velocity correction that considerably reduces unphysical features common in random walk models. It is shown that the model satisfies the well-mixed criterion such that tracer particles retain approximately uniform concentrations when introduced uniformly in the domain, while their deposition velocity is vanishingly small, as it should be. To handle arbitrary geometries, it is assumed that the velocity rms values in the boundary layer can locally be approximated by the DNS data of fully developed channel flows. Benchmarks of the model are performed against particle deposition data in turbulent pipe flows, 90° bends, as well as more complex 3D flows inside a mouth-throat geometry. Good agreement with the data is obtained across the range of particle inertia.     

Stochastic simulations of particle-laden isotropic turbulent flow

Numerical simulations are performed of dispersion and polydispersity of particles in isotropic incompressible turbulence. The mass loading of the particles is assumed to be small; thus the effects of particles on turbulence is neglected (one-way coupling). A stochastic model is employed to simulate the carrier phase. The results of the simulations are compared with direct numerical simulation (DNS) data and theoretical results. The stochastic model predicts most of the trends as portrayed by DNS and theory. However, the continuity effect associated with the crossing trajectories effect is not captured. Also, the peaking in the variation of the particle asymptotic diffusivity coefficient with the particle time constant is not observed. For evaporating particles, the stochastic model predicts thinner probability density functions (pdfs) for the particle diameter as compared with DNS generated pdfs. The model is implemented to investigate the effects of gravity on evaporation. It is shown that the depletion rate increases with increase of the drift velocity at short and intermediate times, but an opposite trend is observed at long times. The standard deviation and skewness of the particle diameter indicate peak values in their variations with the drift velocity. Dispersion of evaporating particles decreases with respect to that of non-evaporating particles at small drift velocities; an opposite trend is observed at large drift velocities. The effects of the initial evaporation rate and the particle Schmidt number on the evaporation in the gravity environment are also studied.     

An incompressible smoothed particle hydrodynamics scheme for Newtonian/non-Newtonian multiphase flows including semi-analytical solutions for two-phase inelastic Poiseuille flows

An incompressible smoothed particle hydrodynamics (ISPH) method is developed for the modeling of multiphase Newtonian and inelastic non-Newtonian flows at low density ratios. This new method is the multiphase extension of Xenakis et al, J. Non-Newtonian Fluid Mech., 218, 1-15, which has been shown to be stable and accurate, with a virtually noise-free pressure field for single-phase non-Newtonian flows. For the validation of the method a semi-analytical solution of a two-phase Newtonian/non-Newtonian (inelastic) Poiseuille flow is derived. The developed method is also compared with the benchmark multiphase case of the Rayleigh Taylor instability and a submarine landslide, thereby demonstrating capability in both Newtonian/Newtonian and Newtonian/non-Newtonian two-phase applications. Comparisons with analytical solutions, experimental and previously published results are conducted and show that the proposed methodology can accurately predict the free-surface and interface profiles of complex incompressible multi-phase flows at low-density ratios relevant, for example, to geophysical environmental applications.     

柔性微粒介电泳分离过程的多尺度模拟

介电泳分离是一种高效的微细颗粒分离技术,利用非均匀电场极化并操纵分离微流道中的颗粒. 柔性微粒在介电泳分离过程中同时受多种物理场、多相流和微粒变形等复杂因素的影响,仅用单一的计算方法对其进行模拟存在一定的难度,本文采用有限单元——格子玻尔兹曼耦合计算的方法处理这一难题.介观尺度的格子玻尔兹曼方法将流体看成由大量微小粒子组成,在离散格子上求解玻尔兹曼输运方程,易于处理多相流及大变形问题,特别适合模拟柔性颗粒在介电泳分离过程中的变形情况.另一方面,介电泳分离过程的模拟需求解流体、电场和微粒运动方程,计算量相当庞大,通过有限单元法求解介电泳力,提高计算效率.利用这种多尺度耦合计算方法,对一款现有的介电泳芯片分离过程进行了模拟.分析了微粒在电场作用下产生的介电泳力,揭示了介电泳力与电场变化率等因素之间的关系.对微粒运动轨迹及其变形的情况进行了研究,发现微粒的变形主要与流体剪切作用有关.这种多尺度耦合计算方法,为复杂微流体的计算提供了一种有效的解决方案.      

Grid-averaged Lagrangian LES model for multiphase turbulent flow

A grid-averaged Lagrangian (GAL) model for dispersed particle motion in multiphase turbulent flow is presented to provide a large eddy simulation (LES) model for multiphase turbulent flow in which a quite large number of particles are involved. The GAL model is based on an averaging operation for a Lagrangian-type equation of motion of a particle over a computational grid volume and a procedure of reallocation of a dispersed particle cloud with its centroid movement to each grid. The model is therefore a mixed Eulerian–Lagrangian model which can effectively reduce computational time compared with existing Lagrangian-type models, without losing the advantage of Lagrangian-type models that they can properly describe the dynamical evolution of particles. Since the GAL model adopts the grid-volume averaging operation it can easily provide an effective SGS model for LES modeling of multiphase turbulent flow. The validity of the multiphase LES model developed, which is named the GAL-LES model, is confirmed through its application to a particle plume, in which the present model is found to simulate large-eddy motion usually observed in a jet and plume, and to give good agreements with experimental data.