Studies on gas turbulence and particle fluctuation in dense gas-particle flows

Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carried out first by using RANS modeling with a Reynolds stress equation turbulence model and sufficiently fine grids, and then by using LES. The turbulence enhancement by the particle wake effect is studied under various particle sizes and relative gas velocities, and the turbulence enhancement is found proportional to the particle diameter and the square of velocity. Based on the above results, a turbulence enhancement model for the particle-wake effect is proposed and is incorporated as a sub-model into a comprehensive two-phase flow model, which is then used to simulate dilute gas-particle flows in a horizontal channel. The simulation results show that the predicted gas turbulence by using the present model accounting for the particle wake effect is obviously in better agreement with the experimental results than the prediction given by the model not accounting for the wake effect. Finally, the proposed model is incorporated into another two-phase flow model to simulate dense gasparticle flows in a downer. The results show that the particle wake effect not only enhances the gas turbulence, but also amplifies the particle fluctuation.     

湍流扩散多相流

湍流扩散多相流在许多工程和环境科学实际应用中十分常见.流体相中的湍流和扩散相的随机本质使得湍流扩散多相流中遇到的问题远比单一流体中遇到的湍流现象复杂得多.首先评述了湍流扩散多相流研究中实验技术和数值计算方面的现状、各自的优点与局限性,以及该领域研究中未来面临的挑战.主要关注湍流扩散多相流研究中以下3方面重点内容:颗粒、液滴和气泡的选择性聚积;湍流对流体相和扩散相耦合作用的影响以及颗粒物和气泡的存在对流体相中湍流的调节机制.     

Evaluation of a Modified Reynolds Stress Model for Turbulent Dispersed Two-Phase Flows Including Two-Way Coupling

A modified Reynolds stress turbulence model for the pressure rate of strain can be derived for dispersed two-phase flows taking into account gas-particle interaction. The transport equations for the Reynolds stresses as well as the equation for the fluctuating pressure can be derived starting from the multiphase Navier–Stokes equations. The unknown pressure rate of strain correlation in the Reynolds stress equations is then modelled by considering the multiphase equation for the fluctuating pressure. This leads to a multiphase pressure rate of strain model. The extra particle interaction source terms occurring in the model for the pressure rate of strain can be constructed easily, with no noticeable extra computational cost. Eulerian–Lagrangian simulation results of a turbulent dispersed two-phase jet are presented to show the differences in results with and without the new two-way coupling terms.     

基于颗粒湍流和颗粒碰撞相互作用规律的湍流模型的研究

颗粒湍流和颗粒碰撞的相互作用规律是两相流动中的核心问题。用颗粒湍流模型和颗粒碰撞的动力论模型叠加的方法在研究两相湍流流动方面取得了一定的成效,但是还有待改进。本文基于颗粒湍流形成大尺度脉动和颗粒间碰撞引起小尺度脉动的概念,从双流体模型出发,建立了两相流动的双尺度kp-pε两相湍流模型。利用该模型对下行床和突扩室内的气固...     

Eulerian prediction of the fluid/particle correlated motion in turbulent two-phase flows

An approximate equation governing the turbulent fluid velocity encountered along discrete particle path is used to derive the fluid/particle turbulent moments required for dispersed two-phase flows modelling. Then, closure model predictions are compared with results obtained from large-eddy simulation of particle fluctuating motion in forced isotropic fluid turbulence.     

A scaling analysis of added-mass and history forces and their coupling in dispersed multiphase flows

Accurate momentum coupling model is vital to simulation of dispersed multiphase flows. The overall force exerted on a particle is divided into four physically meaningful contributions, i.e., quasi-steady, stress-gradient, added-mass, and viscous-unsteady (history) forces. Time scale analysis on the turbulent multiphase flow and the viscous-unsteady kernel shows that the integral representation of the viscous-unsteady force is required except for a narrow range of particle size around the Kolmogorov length scale when particle-to-fluid density ratio is large. Conventionally, the particle-to-fluid density ratio is used to evaluate the relative importance of the unsteady forces (stress-gradient, added-mass, and history forces) in the momentum coupling. However, it is shown from our analysis that when particle-to-fluid density ratio is large, the importance of the unsteady forces depends on the particle-to-fluid length scale ratio and not on the density ratio. Provided the particle size is comparable to the smallest fluid length scale (i.e., Kolmogorov length scale for turbulence or shock thickness for shock-particle interaction) or larger, unsteady forces are important in evaluating the particle motion. Furthermore, the particle mass loading is often used to estimate the importance of the back effect of particles on the fluid. An improved estimate of backward coupling for each force contribution is established through a scaling argument. The back effects of stress-gradient and added-mass forces depend on particle volume fraction. For large particle-to-fluid density ratio, the importance of the quasi-steady force in backward coupling depends on the particle mass fraction; while that of the viscous-unsteady force is related to both particle mass and volume fractions.     

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the single-scale second-order moment two-phase turbulence model superposed with a particle collision model (USM-θ model) in most regions. The project supported by the Special Funds for Major State Basic Research, China (G-1999-0222-08), and the Postdoctoral Science Foundation (2004036239) The English text was polished by Keren Wang     

A nonlinear kp-εp particle two-scale turbulence model and its application

A particle nonlinear two-scale turbulence model is proposed for simulating the anisotropic turbulent two-phase flow. The particle kinetic energy equation for two-scale fluctuation, particle energy transfer rate equation for large-scale fluctuation, and particle turbulent kinetic energy dissipation rate equation for small-scale fluctuation are derived and closed. This model is used to simulate gas–particle flows in a sudden-expansion chamber. The simulation is compared with the experiment and with those obtained by using another two kinds of tow-phase turbulence model, such as the single-scale two-phase turbulence model and the particle two-scale second-order moment (USM) two-phase turbulence model. It is shown that the present model gives simulation in much better agreement with the experiment than the single-scale two-phase turbulence model does and is almost as good as the particle two-scale USM turbulence model. The project supported by China Postdoctoral Science Foundation (2004036239).     

视密度加权的时平均统一二阶矩两相湍流模型

在气粒两相湍流的双流体模型中,颗粒相的视（表观）密度是有脉动的,在时平均的统一二阶矩（USM）模型中出现了和颗粒数密度或视密度脉动有关的项和方程,使模型方程比较复杂。实际上,用LDV或PDPA测量的流体（用小颗粒代表）和颗粒速度都是颗粒数加权平均的结果。因此,在视密度加权平均基础上推导两相湍流模型更为合理。通过推导和封闭了视密度加权平均的统一二阶矩模型（MUSM）方程组,改进了两相速度脉动关联的封闭,并引入了颗粒遇到的气体脉动速度及其输运方程。MUSM模型可以减少所用方程数,节省计算量。视密度加权平均的统一二阶矩两相湍流模型是一种对原有时间平均的统一二阶矩模型和改进和发展。     

Possibilities and Limitations of Computer Simulations of Industrial Turbulent Dispersed Multiphase Flows

We give an overview on the usage of computer simulations in industrial turbulent dispersed multiphase flows. We present a few examples of industrial flows: bubble columns and bubbly pipe flows, stirred tanks, cyclones, and a fluid catalytic cracking unit. The fluid catalytic cracking unit is used to illustrate the complexity of the physical phenomena involved, and the possibilities and limitations of the different approaches used: Eulerian–Lagrangian (particle-tracking) and Eulerian–Eulerian (two-fluid). In the first approach, the continuous phase is solved using either RANS simulations (Reynolds-Averaged Navier–Stokes simulations) or DNS/LES (Direct Numerical Simulations/Large-Eddy Simulations), and the individual particles are tracked. In the second approach, the dispersed phase is averaged, leading to two sets equations, which are quite similar to the RANS equations of single-phase flows. The Eulerian–Eulerian approach is the most commonly used in industrial applications, however, it requires a significant amount of modelling. Eulerian–Lagrangian RANS can be simpler to use; in particular in situations involving complex boundary conditions, polydisperse flows and agglomeration/breakup. The key issue for the success of the simulations is to have good models for the complex physics involved. A major weakness is the lack of good models for: the turbulence modification promoted by the particles, the inter-particle interactions, and the near-wall effects. Eulerian–Lagrangian DNS/LES can play an important role as a research tool, in order to get a better physical understanding, and to improve the models used in the RANS simulations (either Eulerian–Eulerian or Eulerian–Lagrangian).     

On basic equations of turbulent swirling gas-solid flows and their application in cyclones

The basic equations of turbulent gas-solid flows are derived by using the pseudo-fluid model of particle phase with a refined two-phase turbulence model. These equations are then applied to swirling gas-particle flows for analyzing the collection efficiency in cyclone separators.     

Interactions of particles and microbubbles with turbulence

Dilute, dispersed two-phase flows arise in many contexts ranging from solid particles or droplets in gas flows to bubbles in liquids. Many of the flows of interest are turbulent, which presents a complex problem to analyze or to determine the dominant physical processes contributing to the observed phenomena. Advances in experimental techniques have made it possible to measure directly turbulent and particle velocity fluctuations in dilute systems. This has provided a counterpart to advances in computational and analytical models and a basis on which to test these models. Three specific areas are considered: the fluctuating forces on an individual particle in an unsteady flow, the response of a solid particle to a turbulent air flow, and the corresponding response of a small bubble in turbulent liquid flows. Results from direct numerical simulations are presented for each of these, including the nonuniform distribution of particles generated by local instantaneous features of the flow. The issue of turbulence modulation at low to moderate void fractions is discussed.     

离散型湍流多相流动的研究进展和需求

离散型多相流动,指气体-颗粒(气-固)、液体-颗粒(液-固)、液体-气泡、气体-液雾以及气泡-液体-颗粒等两相或三相流动.这种类型的多相流动广泛存在于能源, 航天和航空, 化工和冶金,交通运输, 水利, 核能等领域.本文阐述了离散型多相流动的国内外基础研究,包括颗粒/液滴/气泡在流场中受流体动力作用力的研究, 颗粒-颗粒,液滴-液滴,气泡-气泡之间以及颗粒/液滴和壁面之间碰撞和聚集规律的研究,颗粒-气体和气泡-液体湍流相互作用的研究, 和数值模拟的研究,包括多相流动的雷诺平均模拟、大涡模拟和直接数值模拟的研究进展.最后, 归纳了目前尚待研究的需求.      

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

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

Two-phase turbulence models for simulating dense gas–particle flows

The two-fluid model is widely adopted in simulations of dense gas–particle flows in engineering facilities. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas–particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas–particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-Θ and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle–particle collision using the transport equation model for the two-phase velocity correlation.     

EFFECT OF EMPIRICAL COEFFICIENTS ON SIMULATION IN TWO-SCALE SECOND-ORDER MOMENT PARTICLE-PHASE TURBULENCE MODEL

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.     

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.     

Two-phase turbulence models for simulating dense gas-particle flows

The two-fluid model is widely adopted in simulations of dense gas-particle flows in engineering facili- ties. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas-particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas-particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-O and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle-particle collision using the transport equation model for the two-phase velocity correlation.     

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.