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.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas–solid flows

In gas–solid flows, particle–particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian–Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.     

Measurement on particle rotation speed in gas–solid flow based on identification of particle rotation axis

When a two-dimensional (2D) imaging system is used to visualize particle motion in a 3D gas–solid flow, the particle rotation speed was found extremely difficult to be accurately measured due to the fact that the direction of rotation axis was usually random and hard to be distinguished. The paper presents a method to calculate the particle rotation speed from particle images based on the identification of its rotation axis using two or more characteristic points on its surface. The idea was analyzed and realized in a mathematical way and based on which a calculation program was given. The measurement method was verified with an experiment using a small sphere with known rotation axis and rotation speed. The effects of several factors, including the direction of the particle rotation axis, the particle image resolution, the types and positions of characteristic points, etc., on the measurement error are discussed. The error is found to be acceptable for most cases. The measurement method was finally applied to those small glass beads in a real 3D gas–solid flow inside a cold circulating fluidized bed (CFB) riser, which indicates that the problems of 2D imaging system applying to 3D particulate system could be solved by using this mathematical method.     

Numerical simulations of gas–liquid–solid flows in a hydrocyclone separator

The flow behavior in hydrocyclones is quite complex. In this study, the computational fluid dynamics (CFD) method was used to simulate the flow fields inside a hydrocyclone in order to investigate its separation efficiency. In the computational fluid dynamics study of hydrocyclones, the air-core dimension is a key to predicting the mass split between the underflow and overflow. In turn, the mass split influences the prediction of the size classification curve. Three models, the model, the Reynolds stress model (RSM) without considering the air-core, and the Reynolds stress turbulence model with the volume of fluid (VOF) multiphase model for simulating the air-core, were compared in terms of their predictions of velocity, axial and tangential velocity distributions, and separation proportion. The RSM with air-core simulation model, since it reproduces some detailed features of the turbulence and multiphase, clearly predicted the experimental data more closely than did the other two models.     

Direct numerical simulations of gas–solid–liquid interactions in dilute fluids

The dynamic interactions between gas bubbles, rigid particles and liquid can lead to profound nonlinearities in the aggregate behavior of a multiphase fluid. Predicting the nonlinear dynamics of the multiphase mixture hence requires understanding how the phases interact at the scale of individual interfaces, but these interactions are notoriously difficult to resolve in models. The goal of this paper is to develop and validate a computational method capable of capturing the complex flow interactions between gas bubbles and rigid particles immersed in a Newtonian liquid. We focus on multiphase systems that are dilute enough for the solid and gas components to move through and be moved by the ambient liquid. We use level sets with a topology-preserving advection scheme to track the gas interfaces. To include the motion of the rigid particles, we couple distributed Lagrange multipliers to an immersed-boundary method. The high viscosity contrast between the liquid and the gas requires both time splitting and approximate factorization to efficiently solve the governing equations consisting of the conservation of mass, momentum and energy. To resolve interactions between interfaces that vary drastically in size, we refine our mesh adaptively in the vicinity of the boundary.     

Separation efficiency of liquid–solid undergoing vibration based on breakage of liquid bridge

Vibrating separation is a significant method for liquid–solid separation. A typical example is the vibrating screen to dewater wet granular matter. The properties of granular matter and the vibrating parameters significantly affect the separation efficiency. This study investigates the effect of vibration parameters in separation based on the breakage of large-scale liquid bridge numerically by using a calibrated simulation model. Through analysing the simulation results, the liquid bridge shape and the volume between two sphere particles for various particle sizes and particle distances were studied in the static condition under the effect of gravity. The results show a general reducing trend of liquid bridge volume when the radius ratio of two particles increases, particularly when the ratio increases to 5. Additionally, a set of vibrating motion was applied to the liquid bridge in the simulation model. A group of experiments were also performed to validate the simulation model with vibration. Then, the effect of vibrating peak acceleration, distance between spheres and radius on the separation efficiency which was reflected by the residual water were investigated. It is found that separation efficiency increased obviously with the peak acceleration and the increase slowed down after the peak acceleration over 1 m/s².     

Study of fluid cell coarsening for CFD-DEM simulations of polydisperse gas–solid flows

Particle polydispersity is ubiquitous in industrial fluidized beds, which possesses a significant impact on hydrodynamics of gas–solid flow. Computational fluid dynamics-discrete element method (CFD-DEM) is promising to adequately simulate gas–solid flows with continuous particle size distribution (PSD) while it still suffers from high computational cost. Corresponding coarsening models are thereby desired. This work extends the coarse-grid model to polydisperse systems. Well-resolved simulations with different PSDs are processed through a filtering procedure to modify the gas–particle drag force in coarse-grid simulations. We reveal that the drag correction of individual particle exhibits a dependence on filtered solid volume fraction and filtered slip velocity for both monodisperse and polydisperse systems. Subsequently, the effect of particle size and surrounding PSD is quantified by the ratio of particle size to Sauter mean diameter. Drag correction models for systems with monodisperse and continuous PSD are developed. A priori analysis demonstrates that the developed models exhibit reliable prediction accuracy.     

Order and disorder in particle–liquid flows in a pipe

The spherical expanded polystyrene particle–oil two-phase flow in a vertical pipe was used to simulate the dispersed phase distribution in laminar bubbly flows. A three-dimensional particle image tracking technique was used to track the particles in the flow to study the ordered structure of dispersed phase distribution and its transition to disorder. The ordered structures behaved as particle strings aligned in the flow direction as induced by the flow shear. The structures were quite durable in high liquid velocity flows and dispersed gradually as the liquid velocity decreased. In lower velocity flows, the particles tended to form clusters in the horizontal direction, as predicted by potential theory for spherical bubbles rising in a quiescent inviscid liquid and as observed in experiments on non-shear bubbly water flows.     

Theoretical and numerical studies on the flow multiplicity phenomenon for gas–solids two-phase flows in CFB risers

The dependence of the fully-developed flow profiles on the inlet flow conditions for gas–solids two-phase flows, i.e. the flow multiplicity phenomenon, in circulating fluidized bed (CFB) risers was proposed and discussed in this article. The flow multiplicity phenomenon for gas–solids two-phase flows was first proved mathematically based on the conservation equations of mass and momentum. Then the CFD model using Eulerian–Eulerian approach with k–ε turbulence model for each phase was further adopted to analyze the details of this flow multiplicity phenomenon. It is theoretically and numerically revealed that for gas–solids two-phase flows, the flow profiles in the fully-developed region are always dominated by the flow profiles at the inlet. The solids concentration profile is closely coupled with the velocity profile, and the inlet solids concentration and velocity profiles can largely influence the fully-developed concentration and velocity profiles.     

Direct simulation of liquid–gas–solid flow with a free surface lattice Boltzmann method

Direct numerical simulation of liquid–gas–solid flows is uncommon due to the considerable computational cost. As the grid spacing is determined by the smallest involved length scale, large grid sizes become necessary – in particular, if the bubble–particle aspect ratio is on the order of 10 or larger. Hence, it arises the question of both feasibility and reasonability. In this paper, we present a fully parallel, scalable method for direct numerical simulation of bubble–particle interaction at a size ratio of 1–2 orders of magnitude that makes simulations feasible on currently available super-computing resources. With the presented approach, simulations of bubbles in suspension columns consisting of more than 100,000 fully resolved particles become possible. Furthermore, we demonstrate the significance of particle-resolved simulations by comparison to previous unresolved solutions. The results indicate that fully resolved direct numerical simulation is indeed necessary to predict the flow structure of bubble–particle interaction problems correctly.     

Drag models for simulating gas–solid flow in the turbulent fluidization of FCC particles

This paper examines the suitability of various drag models for predicting the hydrodynamics of the turbulent fluidization of FCC particles on the Fluent V6.2 platform.The drag models included those of Syamlal–O’Brien,Gidaspow,modified Syamlal–O’Brien,and McKeen.Comparison between experimental data and simulated results showed that the Syamlal–O’Brien,Gidaspow,and modified Syamlal–O’Brien drag models highly overestimated gas–solid momentum exchange and could not predict the formation of dense phase in the fl...     

A multi-scale architecture for multi-scale simulation and its application to gas–solid flows

A multi-scale hardware and software architecture implementing the EMMS (energy-minimization multi-scale) paradigm is proven to be effective in the simulation of a two-dimensional gas–solid suspension. General purpose CPUs are employed for macro-scale control and optimization, and many integrated cores (MICs) operating in multiple-instruction multiple-data mode are used for a molecular dynamics simulation of the solid particles at the meso-scale. Many cores operating in single-instruction multiple-data mode, such as general purpose graphics processing units (GPGPUs), are employed for direct numerical simulation of the fluid flow at the micro-scale using the lattice Boltzmann method. This architecture is also expected to be efficient for the multi-scale simulation of other complex systems.     

Closure proposals for the tracking of turbulence-agitated gas–liquid interfaces in stratified flows

The present work focuses on the formulation of new modelling approaches to ensemble-averaged equation for the interface tracking in stratified two-phase flow regime with separated gas and liquid layers. The model aims to represent the influence of the unresolved part of the surface. The modelling constants have been estimated based on the k–? model and the empirical observations of Brocchini and Peregrine (2001a) for the air–water flows. The new modelling proposal was implemented numerically for the 1D case. Results of the simulations are presented in the paper.     

Technique for quantitative mapping of three-dimensional liquid–gas phase boundaries in microchannel flows

A diagnostic technique capable of characterizing interfaces between transparent, immiscible fluids is developed and demonstrated by investigating the morphology of liquid–gas interfaces in an adiabatic two-phase flow through a microchannel of 500 μm × 500 μm square cross section. Water seeded with 0.5 μm-diameter fluorescent polystyrene particles is pumped through the channel, and the desired adiabatic two-phase flow regime is achieved through controlled air injection. The diagnostic technique relies on obtaining particle position data through epifluorescent imaging of the flow at excitation and emission wavelengths of 532 nm and 620 nm, respectively. The particle position data are then used to resolve interface locations to within ±1 μm in the focal plane. By mapping the interface within individual focal planes at various depths within the channel, it is possible to determine the complete liquid–gas interface geometry across the channel cross section in a dynamic flow environment. Utilizing this approach, the liquid–gas phase boundaries of annular flows within a microchannel have been successfully characterized.     

Friction factor correlations for gas–liquid/liquid–liquid flows with drag-reducing polymers in horizontal pipes

In the present paper, two correlations have been developed to predict the effect of DRP on friction factor of two-phase flow for any pipe diameter. The correlations have been verified using published experimental results of DRP added to air–liquid annular flow and for oil–water flows with any flow pattern at the asymptotic state (maximum drag reduction). Such correlations are not available in literature and considered to be very useful in predicting the drag reduction using DRP and in understanding the most significant parameters that affecting the mechanism of drag reduction by polymers in multiphase flow.     

Experimental study and numerical simulation of the characteristics of a percussive gas–solid separator

The presence of solid particles in the flow of hypersonic wind tunnels damages the appearance of the experiment models in the wind tunnel and influences the accuracy of experimental results. The design of a highly efficient gas–solid separator was therefore undertaken. Particle trajectory imaging methods were used to measure trajectories under different conditions. The flow field and particle movement characteristics for different head angles (HAs) and separation tooth angles (STAs), inlet velocities, and the exhaust gas outlet pressures in the separator, were calculated using simulations based on the discrete phase model. The particle separation efficiency, pressure loss, and flow loss resulting from different structural parameters were also studied. In line with experimental observations, the characteristic angle of particle movements in the separator and the separation efficiency of the separator were found to increase with decreasing HA and with increasing STA. Separation efficiency improves with increasing inlet velocity and with increasing negative pressure of the exhaust gas outlet; however, the corresponding pressure loss and the flow rate of the waste gas also increased.     

A compressible two-layer model for transient gas–liquid flows in pipes

This work is dedicated to the modeling of gas–liquid flows in pipes. As a first step, a new two-layer model is proposed to deal with the stratified regime. The starting point is the isentropic Euler set of equations for each phase where the classical hydrostatic assumption is made for the liquid. The main difference with the models issued from the classical literature is that the liquid as well as the gas is assumed compressible. In that framework, an averaging process results in a five-equation system where the hydrostatic constraint has been used to define the interfacial pressure. Closure laws for the interfacial velocity and source terms such as mass and momentum transfer are provided following an entropy inequality. The resulting model is hyperbolic with non-conservative terms. Therefore, regarding the homogeneous part of the system, the definition and uniqueness of jump conditions is studied carefully and acquired. The nature of characteristic fields and the corresponding Riemann invariants are also detailed. Thus, one may build analytical solutions for the Riemann problem. In addition, positivity is obtained for heights and densities. The overall derivation deals with gas–liquid flows through rectangular channels, circular pipes with variable cross section and includes vapor–liquid flows.     

Modeling of multiple noncatalytic gas–solid reactions in a moving bed of porous pellets based on finite volume method

A mathematical model is presented to simulate the multiple heterogeneous reactions with complex set of physicochemical and thermal phenomena in a moving bed of porous pellets. This model is based on both heat and mass transfer phenomena of gaseous species in a porous medium including chemical reactions at interfaces whose areas vary during the conversion. This model accounts for both the exothermic and endothermic reactions which can be equimolar or nonequimolar. Furthermore it considers simultaneously the reactions in the nonisothermal transient condition. A powerful technique based upon finite volume fully implicit approach has been implemented to solve the complicated governing equations numerically. The model has been validated by comparing with various experimental and analytical results in two cases: the single pellet scale as well as the counter current moving bed reactor.

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.