Numerical analysis of swirling turbulent droplet-laden flow and heat transfer in a sudden pipe expansion

The effect of swirling intensity on the structure and heat transfer of a turbulent gas–droplet flow after a sudden pipe expansion has been numerically simulated. Air is used as the carrier phase, and water, ethanol, and acetone are used as the dispersed phase. The Eulerian approach is applied to simulate the dynamics and heat transfer in the dispersed phase. The gas phase is described by a system of Reynolds-averaged Navier-Stokes (RANS) equations, taking into account the effect of droplets on mean transport and turbulent characteristics in the carrier phase. Gas phase turbulence is predicted using the second-moment closure. A swirling droplet-laden flow is characterized by an increase in the number of small particles on the pipe axis due to their accumulation in the zone of flow recirculation and the action of the turbulent migration (turbophoresis) force. A rapid dispersion of fine droplets over the pipe cross-section is observed without swirling. With an increase in swirling intensity, a significant reduction in the length of the separation region occurs. The swirling of a two-phase flow with liquid droplets leads to an increase in the level of turbulence for all three types of liquid droplets investigated in this work due to their intensive evaporation. It is shown that the addition of droplets leads to a significant increase in heat transfer in comparison with a single-phase swirling flow. The greatest effect of flow swirling on heat transfer intensification in a two-phase gas-droplet flow is obtained for the droplets of ethanol and water and smallest effect is for the acetone droplets.     

Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.     

Numerical investigation of flow and heat transfer in a gas-droplet separated turbulent stream using the Reynolds stress transfer model

The results of the numerical modeling of flow structure, turbulence, and heat transfer in a gas-droplet stream after sudden tube expansion on the basis of the Eulerian approach are presented. The gas phase turbulence was modeled using the Reynolds stress transfer model modified to allow for the presence of particles. The results are compared with those obtained using the two-equation k-ε model. The latter results overestimate the heat transfer in the separation flow as compared with the Reynolds stress transfer model. The heat transfer is shown to considerably increase, when evaporating droplets are incorporated in the separation flow (by a factor of more than 1.5 compared with the case of a single-phase flow at a small mass concentration of the droplets M _L1 ≤ 0.05). The addition of the disperse phase in the turbulent gas flow leads a slight increase in the recirculation zone length. Good agreement with the experimental data indicates the adequacy of the numerical model developed.     

Unsteady turbulent flow of a gas suspension in a channel under conditions of injection and forced pressure oscillations

A turbulent flow of a suspension of solid particles in a gas is considered. The suspension is located in a channel with permeable walls (the pressure at the left end face of the channel follows a sinusoidal law). The flow considered here reflects the principal features of the flow in the combustion chamber of a solid-propellant rocket motor. The unsteady flow of the gas suspension is described by using the Eulerian-Lagrangian approach. A stochastic variant of the discrete-trajectory approach is used for modeling the particle motion. The influence of the condensed phase on the turbulence characteristics and acoustic oscillations of the parameters of the working medium in the channel in the case of injection is discussed. The calculated results are compared with data obtained in a physical experiment.     

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.     

Validation of a numerical model for the simulation of an electrostatic powder coating process

Numerical modeling of a complete powder coating process is carried out to understand the gas-particle two-phase flow field inside a powder coating booth and results of the numerical simulations are compared with experimental data to validate the numerical results. The flow inside the coating booth is modeled as a three-dimensional turbulent continuous gas flow with solid powder particles as a discrete phase. The continuous gas flow is predicted by solving Navier–Stokes equations using a standard k–ε turbulence model with non-equilibrium wall functions. The discrete phase is modeled based on a Lagrangian approach. In the calculation of particle propagation, a particle size distribution obtained through experiments is applied. The electrostatic field, including the effect of space charge due to free ions, is calculated with the use of the user defined scalar transport equations and user defined scalar functions in the software package, FLUENT, for the electrostatic potential and charge density.     

Prediction of the particle-laden jet with a two-equation turbulence model

A two-equation turbulence model for two-phase flows has recently been proposed by Elghobashi & Abou-Arab (1983). They derived the exact equations of the kinetic energy of turbulence and the rate of dissipation of that energy, and modeled the turbulent correlations, resulting from time-averaging, up to third order. In order to validate the proposed model, a turbulent axisymmetric gaseous jet laden with spherical uniform-size solid particles is studied here. The predictions of the mean flow properties of the two-phases and the turbulence kinetic energy and shear stress of the carrier phase show good agreement with the experimental data.     

A two-equation turbulence model for jet flows laden with vaporizing droplets

A two-equation turbulence model for steady incompressible two-phase flows including phase change has been recently developed by Mostafa & Elghobashi (1984). This model is tested for the flow of a turbulent axisymmetric gaseous jet laden with evaporating liquid droplets. To avoid the problem of density fluctuations of the carrier phase at this stage, only isothermal flow is considered and vaporization is assumed to be due to the vapor concentration gradient. The continuous size distribution of the droplets is approximated by finite size groups. Each group is considered as a continuous phase interpenetrating and interacting with the carrier phase. Two test cases have been predicted by the model. The first is for a Freon-11 spray issuing from a round nozzle, where experimental data are available at distances equal to or greater than 170 nozzle diameters. Good agreement between the data and the predictions was achieved. The second is for a methanol spray where no experiments are available yet and the predictions consider the flow region close to the nozzle ($\text{z/D < 40}$). The results of the methanol spray include distributions of the mean velocity, volume fractions of the different phases, concentration of the evaporated material in the carrier phase, turbulence intensity and shear stress of the carrier phase, droplet diameter distribution, and the jet spreading rate. In this case the results are analyzed based on a qualitative comparison with the corresponding single phase jet flow.     

Turbulence modulation in dilute particle-laden flow

A new particle source term to account for the effect of particles on the turbulence equations based on the Euler/Lagrange approach is introduced and compared with existing models and experimental data. Three different sizes of particles are considered to cover the range of large particles, where augmentation of the carrier phase turbulence is expected, and small particles, for which attenuation is expected. The new model is derived directly from the balance equations of fluid flow and represents a combination of the so-called standard and consistent approaches. The performance of the new model surpasses that of the standard and consistent models and it is able to predict both the suppression and enhancement of fluid turbulence for small and large particles.     

Particles-induced turbulence: A critical review of physical concepts,numerical modelings and experimental investigations

The presence of solid particles or water droplets in continuous fluid flow can either induce turbulence attenuation or amplification. The modification of the state of the turbulence depends on the characteristics of the particles, such as volume fraction, mean diameter, mass density, or carrier phase flow properties. In this brief review, the main physical concepts related to the most important physical aspects of turbulence modulation are summarized. Different criteria used to distinguish the enhancement or the attenuation effects of the particles on the carrier phase flows are recalled. For the interest of large-scale industrial applications, several theoretical,experimental and empirical approaches are discussed, which provides an interesting framework for the study of the effect of particles on turbulence behavior modification.     

Evaporation of polydispersed droplets in a highly turbulent channel flow

A model experiment for the study of evaporating turbulent two-phase flows is presented here. The study focuses on a situation where pre-atomized and dispersed droplets vaporize and mix in a heated turbulent flow. The test bench consists in a channel flow with characteristics of homogeneous and isotropic turbulence where fluctuations levels reach very high values (25% in the established zone). An ultrasonic atomizer allows the injection of a mist of small droplets of acetone in the carrier flow. The large range diameters ensure that every kind of droplet behavior with regards to turbulence is possible. Instantaneous concentration fields of the vaporized phase are extracted from fluorescent images (PLIF) of the two phase flow. The evolution of the mixing of the acetone vapor is analyzed for two different liquid mass loadings. Despite the high turbulence levels, concentration fluctuations remain significant, indicating that air and acetone vapor are not fully mixed far from the injector.     

稠密气固两相湍流流动的实验和数值模拟

基于气固两相流动模型计算循环流化床内稠密气固两相流湍流动,颗粒动理学方法模拟颗粒相湍动能,ＳＧＳ模型模拟气相湍流,采用γ－射线密度计和非等速取样管测量局部颗粒浓度和流率,利用ＦＦＴ方法计算颗粒浓度功率谱密度。模拟计算得到上升管内气相和固相速度和浓度分布等。同时数值模拟与Ｔｓｕｊｉ等和Ｋｎｏｗｌｔｏｎ等试验结果进行了比较,结果表明数值模拟计算与实验结果相吻合。     

Particle statistics in a gas–solid coaxial strongly swirling flow: A direct numerical simulation

A direct numerical simulation of a strongly coaxial swirling particle-laden flow is conducted with reference to a previous experiment. The carrier phase is simulated as a coaxial swirling flow through a short nozzle injecting into a large container. The particle phase is carried by the primary jet, and simulated in the Lagrangian approach. The drag force, slip-shear force and slip-rotation force experienced by particles are calculated. A partial validation of the results is followed. The results are analyzed in Eulerian approach focusing on the statistical behavior of particle motion. The relative importance of the drag, slip-shear and slip-rotation forces under different Stokes numbers is indicated quantitatively. The particle velocity profiles, fluctuations, Reynolds stress, and turbulence intensity are demonstrated and analyzed respectively. An important “choke” behavior for large particles within the mainstream is found and interpreted. Additionally, the patterns of particle distribution and the helical structures of particle motion under different Stokes numbers are demonstrated qualitatively and analyzed quantitatively.     

Structure of the nonisothermal swirling gas-droplet flow behind an abrupt tube expansion

Flow structure and heat and mass transfer in a swirling two-phase stream is numerically modeled using the Reynolds stress transport model. The gas phase is described by the 3DRANS system of equations with account for the inverse influence of particles on the transport processes in the gas. The gas phase turbulence is calculated using the Reynolds stress transport model with account for the presence of disperse particles. The two-phase nonswirling flow behind an abrupt tube expansion contains a secondary corner vortex which is absent from the swirling flow. The disperse phase is redistributed over the tube cross-section. Large particles are concentrated in the wall region of the channel under the action of the centrifugal forces, while the smaller particles are in the central zone of the chamber.     

Investigation of turbulence modification in a non-reactive two-phase flow

In a two-phase flow the influence of a dispersed phase on the turbulence properties of a continuous phase, known as turbulence modification, is investigated. An experimental approach is discussed that is suitable for studying the decay of grid-generated turbulence in a vertically orientated wind tunnel with a cross-section large enough to avoid influences from walls. Phase Doppler anemometry is used to characterize both single and two-phase flow by measuring mean axial and radial velocity components, velocity fluctuations, turbulent kinetic energy, and integral time scales. By direct comparison of results from single- and two-phase flows, the feedback of the dispersed phase on the continuous phase can be isolated. The data is used to deduce a source term for particle-induced turbulence production appropriate for a numerical simulation of the flow, based on the Reynolds-averaged Navier-Stokes equations. Although of special importance for a detailed understanding of turbulent two-phase combustion, additional complexity introduced by evaporation and chemical reactions is avoided by using glass beads as dispersed phase.     

Prediction of mono-disperse gas–liquid turbulent flow in a vertical pipe

A two-fluid model in the Eulerian–Eulerian framework has been implemented for the prediction of gas volume fraction, mean phasic velocities, and the liquid phase turbulence properties for gas–liquid upward flow in a vertical pipe. The governing two-fluid transport equations are discretized using the finite volume method and a low Reynolds number $k-\varepsilon$ model is used to predict the turbulence field for the continuous liquid phase. In the present analysis, a fully developed one-dimensional flow is considered where the gas volume fraction profile is predicted using the radial force balance for the bubble phase. The current study investigates: (1) the turbulence modulation terms which represent the effect of bubbles on the liquid phase turbulence in the k–ε transport equations; (2) the role of the bubble induced turbulent viscosity compared to turbulence generated by shear; and (3) the effect of bubble size on the radial forces which results in either a center-peak or a wall-peak in the gas volume fraction profiles. The results obtained from the current simulation are generally in good agreement with the experimental data, and somewhat improved over the predictions of some previous numerical studies.     

Characterization of particle-laden, confined swirling flows by phase-doppler anemometry and numerical calculation

The particle dispersion characteristics in a confined swirling flow with a swirl number of approx. 0.5 were studied in detail by performing measurements using phase-Doppler anemometry (PDA) and numerical predictions. A mixture of gas and particles was injected without swirl into the test section, while the swirling airstream was provided through a co-flowing annular inlet. Two cases with different primary jet exit velocities were considered. For these flow conditions, a closed central recirculation bubble was established just downstream of the inlet.

The PDA measurements allowed the correlation between particle size and velocity to be obtained and also the spatial change in the particle size distribution throughout the flow field. For these results, the behaviour of different size classes in the entire particle size spectrum, ranging from about 15 to 80 μm, could be studied, and the response of the particles to the mean flow and the gas turbulence could be characterized. Due to the response characteristics of particles with different diameters to the mean flow and the flow turbulence, a considerable separation of the particles was observed which resulted in a streamwise increase in the particle mean number diameter in the core region of the central recirculation bubble. For the lower particle inlet velocity (i.e. low primary jet exit velocity), this effect is more pronounced, since here the particles have more time to respond to the flow reversal and the swirl velocity component. This also gave a higher mass of recirculating particle material.

The numerical predictions of the gas flow were performed by solving the time-averaged Navier-Stokes equations in connection with the well known kε turbulence model. Although this turbulence model is based on the assumption of isotropic turbulence, the agreement of the calculated mean velocity profiles compared to the measured gas velocities is very good. The gas-phase turbulent kinetic energy, however, is considerably underpredicted in the initial mixing region. The particle dispersion characteristics were calculated by using the Lagrangian approach, where the influence of the particulate phase on the gas flow could be neglected, since only very low mass loadings were considered. The calculated results for the particle mean velocity and the mass flux are also in good agreement with the experiments. Furthermore, the change in the particle mean diameter throughout the flow field was predicted approximately, which shows that the applied simple stochastic dispersion model also gives good results for such very complex flows. The variation of the gas and particle velocity in the primary inlet had a considerable impact on the particle dispersion behaviour in the swirling flow and the particle residence time in the central recirculation bubble, which could be determined from the numerical calculations. For the lower particle inlet velocity, the maximum particle size-dependence residence time within the recirculation region was considerably shifted towards larger particles.     

液固两相槽道流中湍流调制的数值研究

结合雷诺应力模型 (Reynolds Stress Model, RSM) 和混合模型 (Mixture Model) 对槽道湍流向下流动中的颗粒调制湍流问题进行了研究．该模型考虑了颗粒流的动能理论和颗粒对湍流的反馈作用．着重分析了颗粒对湍流的调制作用,以及颗粒参数变量（如颗粒密度和质量载荷）对湍流调制的影响．结果表明：(1)在颗粒抑制湍流的范围内,当颗粒密度小于载流体密度时,湍流强度的改变量与颗粒密度成反比;当颗粒密度大于载流体密度时,湍流强度的改变量与颗粒密度成正比;(2)在一定范围内,颗粒抑制湍流的能力随颗粒质量载荷增加而变强．     

Gas turbulence modulation in the pneumatic conveying of massive particles in vertical tubes

The present work is concerned with the interaction between large particles and gas phase turbulence. Gas turbulence modulation in these systems is considered to be dominated by a generation mechanism which arises due to the presence of wakes behind particles. Following a recent proposal, a closure for gas turbulence modulation accounting for the effect of wakes is employed within the context of a mathematical model for particle-laden, turbulent flows. The model accounts for particle particle and particle-wall interactions associated with larger particles based on concepts from gas kinetic theory. It is shown that due to the significant flattening of the mean gas velocity profile with the addition of particles, and the corresponding decrease in turbulent energy production, a generation mechanism must be present in order to produce gas velocity fluctuation predictions which are consistent with the experimental measurements, even in the case where the experimental results indicate a net suppression of gas phase turbulence in the presence of particles.