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.     

Particle/Turbulence Interactions, Mass Transfer and Gas/Solid Chemistry in a CFBC Riser

Gas/solid chemistry in the upper, dilute region of a circulating fluidised bed combustor (CFBC) riser is addressed. The limitations of turbulent mixing are illustrated by the example of the heterogeneous NO/CO/char reaction, relevant in CFB combustion of nitrogen-containing solid fuels. The mass transfer of the gaseous reactants to the char surface is determined, and how the conversion is influenced by the degree of mixing of the multiphase system by turbulent dispersion. Particle/turbulence interactions are taken into account by a (Lagrangian) frequency spectrum of velocity fluctuations, which determine the dispersion of the char particles described here with the Tchen–Hinze model. Chars from solid fuels characterised by fuel ratio (FR) ranging from 0.1 (wood) via 0.5 (peat) and 1 (coal A) to 2 (coal B) were considered. The effective rate of the NO/CO/char reaction is determined as a function of the size and type of the char particle, temperature, particle concentration, reactor dimensions and fluidization velocity, at atmospheric pressure. It was found that for this case the effective gas/solid conversion rate in the upper, dilute region of the riser is much lower than the gas/solid chemistry, mainly due to mass transfer limitations for char particles with sizes of typically 300 μm. The concentration of NO at the char particle surface is only a few % of that in the bulk gas phase. Strong influences were found for particle size and temperature, whilst the fluidisation velocity and the reactor size have only a small influence. It is concluded that for a typical CFB riser, for particles larger than approx. 20 μm, mass transfer has a stronger influence on the heterogeneous NO/CO/char reduction mechanism than the ``unmixedness' due to particle eddy dispersion limitations. It is recommended that this or a similar approach to turbulent dispersive mixing is implemented in CFD codes when these are used for boiler and furnace calculations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Investigation of particle-laden turbulent pipe flow at high-Reynolds-number using particle image/tracking velocimetry (PIV/PTV)

An experimental investigation of a high Reynolds number flow (Re = 320 000) of a dilute liquid-solid mixture (<1% by volume) was conducted. The turbulent motion of both the liquid phase (water) and particles (0.5, 1, and 2 mm glass beads) was evaluated in an upward pipe flow using a particle image/tracking velocimetry (PIV/PTV) technique. Results show that the Eulerian mean axial velocity of the glass beads is lower than that of the liquid phase in the central region but higher in the near-wall region. Moreover, the presence of the coarse particles has a negligible effect on the turbulence intensity of the liquid phase. Particles show higher streamwise and radial fluctuations than the liquid-phase at the tested conditions. The profiles of particle concentration across the pipe radius show almost constant concentration in the core of the pipe with a decrease towards the near wall region for 0.5 and 1 mm particles. For the 2 mm particles, a nearly linear concentration gradient from centre to the pipe wall is observed. The results presented here provide new information concerning the effect of a dispersed particulate phase on the turbulence modulation of the liquid carrier phase, especially at high Reynolds numbers. The present study also demonstrates how correlations developed to determine if particles cause turbulence attenuation/augmentation are not applicable for solid-liquid flows at high Reynolds numbers. Finally, the importance of particle-fluid slip velocity on fluid phase turbulence modulation is illustrated.     

Experimental investigation of turbulence modulation in particle-laden coaxial jets by Phase Doppler Anemometry

The effect of solid particles on the flow characteristics of axisymmetric turbulent coaxial jets for two flow conditions was studied. Simultaneous measurements of size and velocity distributions of continuous and dispersed phases in a two-phase flow are presented using a Phase Doppler Anemometry (PDA) technique. Spherical glass particles with a particle diameter range from 102 to 212 μm were used in this two-phase flow, the experimental results indicate a significant influence of the solid particles and the Re on the flow characteristics. The data show that the gas phase has lower mean velocity in the near-injector region and a higher mean velocity at the developed region. Near the injector at low Reynolds number (Re = 2839) the presence of the particles dampens the gas-phase turbulence, while at higher Reynolds number (Re = 11 893) the gas-phase turbulence and the velocity fluctuation of particle-laden jets are increased. The particle velocity at higher Reynolds number (Re = 11 893) and is lower at lower Reynolds number (Re = 2839). The slip velocity between particles and gas phase existed over the flow domain was examined. More importantly, the present experiment results suggest that, consideration of the gas characteristic length scales is insufficient to predict gas-phase turbulence modulation in gas-particle flows.     

Hydrodynamic and thermal fields analysis in gas-solid two-phase flow

The present work aims to investigate numerically the flowfield and heat transfer process in gas-solid suspension in a vertical pneumatic conveying pipe. The Eulerian-Lagrangian model is used to simulate the flow of the two-phases. The gas phase is simulated based on Reynolds Average Navier-Stokes equations (RANS) with low Reynolds number k-ε model, while particle tracking procedure is used for the solid phase. An anisotropic model is used to calculate the Reynolds stresses and the turbulent Prandtl number is calculated as a function of the turbulent viscosity. The model takes into account the lift and drag forces and the effect of particle rotation as well as the particles dispersion by turbulence effect. The effects of inter-particles collisions and turbulence modulation by the solid particles, i.e. four-way coupling, are also included in the model. Comparisons between different models for turbulence modulation with experimental data are carried out to select the best model. The model is validated against published experimental data for velocities of the two phases, turbulence intensity, solids concentration, pressure drop, heat transfer rates and Nusselt number distribution. The comparisons indicate that the present model is able to predict the complex interaction between the two phases in non-isothermal gas-solid flow in the tested range. The results indicate that the particle-particle collision, turbulence dispersion and lift force play a key role in the concentration distribution. In addition, the heat transfer rate increases as the mass loading ratio increases and Nusselt number increases as the pipe diameter increases.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.     

Calculation of a turbulent two-phase isobaric jet

A numerical calculation is carried out by the finite-difference method based on proposed equations for a turbulent submerged jet containing an admixture of solid particles. The relative longitudinal particle velocity and the influence of particles on the turbulence intensity are taken into account. The calculated results adequately agree with available experimental data. A turbulent two-phase jet is examined in [1] on the basis of the theory for a variable density jet, assuming equal mean velocities for the gas and particles and not considering the influence of particles on the turbulence intensity. Particles are analogously taken into account by a noninertial gas mixture in [2, 3], and a particle Schmidt number of 1.1 is assumed in [4]. A model is proposed in [5] which takes into account the influence of particles on the turbulence intensity of the gas phase. Problems concerning the initial and main sections of a submerged jet were solved in [6] by the integral method on the basis of this model and the assumed equality of the mean velocities of the gas and particles. Turbulent mixing of homogeneous two-phase flows with allowance made for dynamic nonequilibrium of the phases is considered in [7]. However, the neglect of turbulent transfer of particle mass and momentum led to a physically unrealistic solution for the particle concentration in the far field of the mixture. A two-phase jet is considered in the present work on the basis of the theory of a two-velocity continuous medium [8, 9] with allowance made for turbulent transfer of particle mass and momentum. The influence of particles on the turbulence intensity of the gas phase is taken into account with the model of [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 57–63, September–October, 1976.The author acknowledges useful comments and discussion.of the work by G. N. Abramovich and participants of his seminar. The author sincerely thanks I. N. Murzinov for scientific supervision of the work.     

Experimental studies on particle behaviour and turbulence modification in horizontal channel flow with different wall roughness

Detailed measurements in a developed particle-laden horizontal channel flow (length 6 m, height 35 mm, the length is about 170 channel heights) are presented using phase-Doppler anemometry for simultaneous determination of air and particle velocity. The particles were spherical glass beads with mean diameters in the range of 60 µm-1 mm. The conveying velocity could be varied between about 10 m/s and 25 m/s, and the particle mass loading could reach values of about 2 (the mass loading is defined as the ratio of particle to gas phase mass flow rates), depending on particle size. For the first time, the degree of wall roughness could be modified by exchanging the wall plates. The influence of these parameters and the effect of inter-particle collisions on the profiles of particle mean and fluctuating velocities and the normalised concentration in the developed flow were examined. It was shown that wall roughness decreases the particle mean velocity and enhances fluctuating velocities due to irregular wall bouncing and an increase in wall collision frequency, i.e. reduction in mean free path. Thereby, the larger particles are mainly more uniformly distributed across the channel, and gravitational settling is reduced. Both components of the particle velocity fluctuation were reduced with increasing mass loading due to inter-particle collisions and the momentum loss involved. Moreover, the effect of the particles on the air flow and the turbulent fluctuations was studied on the basis of profiles in the developed flow and turbulence spectra determined for the streamwise velocity component. In addition to the effect of particle size and mass loading on turbulence modulation, the influence of wall roughness was analysed. It was clearly shown that increasing wall roughness also results in a stronger turbulence dissipation due to two-way coupling.     

Modulation of homogeneous turbulence seeded with finite size bubbles or particles

The dynamics of homogeneous, isotropic turbulence seeded with finite sized particles or bubbles is investigated in a series of numerical simulations, using the force-coupling method for the particle phase and low wavenumber forcing of the flow to sustain the turbulence. Results are given on the modulation of the turbulence due to massless bubbles, neutrally buoyant particles and inertial particles of specific density 1.4 at volumetric concentrations of 6%. Buoyancy forces due to gravity are excluded to emphasize finite size and inertial effects for the bubbles or particles and their interactions with the turbulence. Besides observing the classical entrapment of bubbles and the expulsion of inertial particles by vortex structures, we analyze the Lagrangian statistics for the velocity and acceleration of the dispersed phase. The turbulent fluctuations are damped at mid-range wavenumbers by the bubbles or particles while the small-scale kinetic energy is significantly enhanced. Unexpectedly, the modulation of turbulence depends only slightly on the dispersion characteristics (bubble entrapment in vortices or inertial sweeping of the solid particles) but is closely related to the stresslet component (finite size effect) of the flow disturbances. The pivoting wavenumber characterizing the transition from damped to enhanced energy content is shown to vary with the size of the bubbles or particles. The spectrum for the energy transfer by the particle phase is examined and the possibility of representing this, at large scales, through an additional effective viscosity is discussed.     

Turbulent thermal diffusion in a multi-fan turbulence generator with imposed mean temperature gradient

We studied experimentally the effect of turbulent thermal diffusion in a multi-fan turbulence generator which produces a nearly homogeneous and isotropic flow with a small mean velocity. Using particle image velocimetry and image processing techniques, we showed that in a turbulent flow with an imposed mean vertical temperature gradient (stably stratified flow) particles accumulate in the regions with the mean temperature minimum. These experiments detected the effect of turbulent thermal diffusion in a multi-fan turbulence generator for relatively high Reynolds numbers. The experimental results are in compliance with the results of the previous experimental studies of turbulent thermal diffusion in oscillating grid turbulence (Buchholz et al. 2004; Eidelman et al. 2004). We demonstrated that the turbulent thermal diffusion is an universal phenomenon. It occurs independently of the method of turbulence generation, and the qualitative behavior of particle spatial distribution in these very different turbulent flows is similar. Competition between turbulent fluxes caused by turbulent thermal diffusion and turbulent diffusion determines the formation of particle inhomogeneities.     

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.     

Non-equilibrium turbulent phenomena in the ?ow over a backward-facing ramp

Non-equilibrium turbulence phenomena have raised great interests in recent years. Significant efforts have been devoted to non-equilibrium turbulence properties in canonical flows, e.g., grid turbulence, turbulent wakes, and homogeneous isotropic turbulence(HIT). The non-equilibrium turbulence in non-canonical flows, however, has rarely been studied due to the complexity of the flows. In the present contribution, a directnumerical simulation(DNS) database of a turbulent flow is analyzed over a backwardfacing ramp, the flow near the boundary is demonstrated, and the non-equilibrium turbulent properties of the flow in the wake of the ramp are presented by using the characteristic parameters such as the dissipation coefficient C and the skewness of longitudinal velocity gradient Sk, but with opposite underlying turbulent energy transfer properties. The equation of Lagrangian velocity gradient correlation is examined, and the results show that non-equilibrium turbulence is the result of phase de-coherence phenomena, which is not taken into account in the modeling of non-equilibrium turbulence. These findings are expected to inspire deeper investigation of different non-equilibrium turbulence phenomena in different flow conditions and the improvement of turbulence modeling.     

Modeling fluid-particle interaction in dilute-phase turbulent liquid-particle flow simulation

The present work examines the predictive capability of a two-fluid CFD model that is based on the kinetic theory of granular flow in simulating dilute-phase turbulent liquid-particle pipe flows in which the inter-stitial fluid effect on the particle fluctuating motion is significant.The impacts of employing different drag correlations and turbulence closure models to describe the fluid-particle interactions(i.e.drag force and long-range interaction)are examined at both the mean and fluctuating velocity levels.The model pre-dictions are validated using experimental data of turbulent liquid-particle flows in a vertical pipe at different particle Reynolds numbers(ReP > 400 and ReP < 400),which characterize the importance of the vortex shedding phenomenon in the fluid-phase turbulence modulation.The results indicate that(1)the fluctuating velocity level predictions at different ReP are highly sensitive to the drag correlation selec-tion and(2)different turbulence closure models must be employed to accurately describe the long-range fluid-particle interaction in each phase.In general,good agreement is found between the model predic-tions and the experimental data at both the mean and fluctuating velocity levels provided that appropriate combinations of the drag correlation and the turbulence closure model are selected depending on Rep.     

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.     

Measurement of turbulent spectra of particle-laden flows in a horizontal channel

Frequency spectra of air turbulence of particle-laden flows were investigated by use of a laser-Doppler velocimeter to discover the eddy-length scales that are influenced by the transported particles. The influence of glass and steel particles of 100–1,000 μm diameter was measured in a horizontal channel and a horizontal pipe for the streamwise and transverse components of the velocity vector. Particles that were small compared with the integral length scale of the particle-laden flow decrease the turbulent power density of the greatest eddies in varying degrees, depending on mass loading and distance from the wall. All fractions create turbulence in their wakes, the size of which depends on loading and slip velocity. These results support the hypothesis that the particles consume energy by following the large eddies that are much greater than the particle diameters, and in so doing, turbulence is created by this energy. Received: 28 September 2000/Accepted: 9 April 2001     

Statistical models for predicting the effect of bidisperse particle collisions on particle velocities and stresses in homogeneous anisotropic turbulent flows

The purpose of this paper is to present and compare two statistical models for predicting the effect of collisions on particle velocities and stresses in bidisperse turbulent flows. These models start from a kinetic equation for the probability density function (PDF) of the particle velocity distribution in a homogeneous anisotropic turbulent flow. The kinetic equation describes simultaneously particle–turbulence and particle–particle interactions. The paper is focused on deriving the collision terms in the governing equations of the PDF moments. One of the collision models is based on a Grad-like expansion for the PDF of the velocity distributions of two particles. The other model stems from a Grad-like expansion for the joint fluid–particle PDF. The validity of these models is explored by comparing with Lagrangian simulations of particle tracking in uniformly sheared and isotropic turbulent flows generated by LES. Notwithstanding the fact that the fluid turbulence may be isotropic, the particle velocity fluctuations are anisotropic due to the impact of gravitational settling. Comparisons of the model predictions and the numerical simulations show encouraging agreement.     

Turbulent transport dissimilarity with modulated turbulence structure in channel flow of viscoelastic fluid

This experimental study investigated the turbulent transport dissimilarity with a modulated turbulence structure in a channel flow of a viscoelastic fluid using simultaneous particle image velocimetry and planar laser-induced fluorescence measurements. An instantaneous dye concentration field with fluctuating velocity vectors showed that mass was transferred by hierarchically large-scale wavy motions with inclination. A co-spectral analysis showed that the spatial phase modulation of the streamwise velocity and dye concentration fluctuations for the wall-normal velocity fluctuation corresponded to the relaxation time. The occurrence of intense dye concentration fluctuation and small streamwise velocity fluctuation in a thin boundary layer caused dissimilar turbulent transport because of the non-zero negative correlation of the streamwise velocity and dye concentration fluctuations for the wall-normal velocity fluctuation only on large scales. This explains the turbulent transport dissimilarity which leads to the zero averaged Reynolds shear stress and non-zero wall-normal turbulent mass flux.     

Experimental and computational studies on flow behavior of gas-solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally.Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695--0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

A laser-fluorescence technique for turbulent two-phase flow measurements

A non-intrusive measurement technique has been developed for accurate determination of gas and particle velocities in a turbulent two-phase flow field. The principle of the technique is based on the discrimination between the scattered light from particles and the fluorescence emission from particles coated with a fluorescent dye. A high-powered, argon-ion based, single-channel, on-axis backscatter laser-Doppler velocimetry system was used. The fluorescent dye was Rhodamine 6G. A study of the gas-solid two-phase flow behaviour in the freeboard of a cold gas-fluidized bed was undertaken. The solid phase contained two particle groups: bed material (sand) and fuel particles (wood). Measurements of the axial velocity and turbulence intensity distributions of the gas phase and both particle groups within the solid phase were made along the column centre and across the freeboard. Excellent discrimination of velocities from the two phases and from the two particle groups within the solid phase was achieved.