A Lagrangian–Eulerian model of particle dispersion in a turbulent plane mixing layer

A Lagrangian–Eulerian model for the dispersion of solid particles in a two‐dimensional, incompressible, turbulent flow is reported and validated. Prediction of the continuous phase is done by solving an Eulerian model using a control‐volume finite element method (CVFEM). A Lagrangian model is also applied, using a Runge–Kutta method to obtain the particle trajectories. The effect of fluid turbulence upon particle dispersion is taken into consideration through a simple stochastic approach. Validation tests are performed by comparing predictions for both phases in a particle‐laden, plane mixing layer airflow with corresponding measurements formerly reported by other authors. Even though some limitations are detected in the calculation of particle dispersion, on the whole the validation results are rather successful. Copyright © 2002 John Wiley & Sons, Ltd.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Evaluation of erosion in a two‐way coupled fluid–particle system

In this study, the effects of flow turbulence intensity, temperature, particle sizes and impinging velocity on erosion by particle impact are demonstrated numerically. Underlying turbulent flow on an Eulerian frame is described by the compressible Reynolds averaged Navier–Stokes equations with a RNG k–ε turbulence model. The particle trajectories and particle–wall interactions are evaluated by a Eulerian–Lagrangian approach in a two‐way coupling system. An erosion model considering material weight removal from surfaces is used to predict erosive wear. Computational validation against measured data is demonstrated satisfactorily. The analysis of erosion shows that the prevention of erosion is enhanced by increasing the effects of flow temperature and turbulence intensity and reducing particle inertial momentum. Copyright © 2001 John Wiley & Sons, Ltd.     

The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method

We perform direct numerical simulation of three‐dimensional turbulent flows in a rectangular channel, with a lattice Boltzmann method, efficiently implemented on heavily parallel general purpose graphical processor units. After validating the method for a single fluid, for standard boundary layer problems, we study changes in mean and turbulent properties of particle‐laden flows, as a function of particle size and concentration. The problem of physical interest for this application is the effect of water droplets on the turbulent properties of a high‐speed air flow, near a solid surface. To do so, we use a Lagrangian tracking approach for a large number of rigid spherical point particles, whose motion is forced by drag forces caused by the fluid flow; particle effects on the latter are in turn represented by distributed volume forces in the lattice Boltzmann method. Results suggest that, while mean flow properties are only slightly affected, unless a very large concentration of particles is used, the turbulent vortices present near the boundary are significantly damped and broken down by the turbulent motion of the heavy particles, and both turbulent Reynolds stresses and the production of turbulent kinetic energy are decreased because of the particle effects. We also find that the streamwise component of turbulent velocity fluctuations is increased, while the spanwise and wall‐normal components are decreased, as compared with the single fluid channel case. Additionally, the streamwise velocity of the carrier (air) phase is slightly reduced in the logarithmic boundary layer near the solid walls. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

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.     

A numerical continuous model for the hydrodynamics of fluid particle systems

In order to understand the hydrodynamic interactions that can appear in a fluid particle motion, an original method based on the equations governing the motion of two immiscible fluids has been developed. These momentum equations are solved for both the fluid and solid phases. The solid phase is assumed to be a fluid phase with physical properties, such as its behaviour can be assimilated to that of pseudo‐rigid particles. The only unknowns are the velocity and the pressure defined in both phases. The unsteady two‐dimensional momentum equations are solved by using a staggered finite volume formulation and a projection method. The transport of each particle is solved by using a second‐order explicit scheme. The physical model and the numerical method are presented, and the method is validated through experimental measurements and numerical results concerning the flow around a circular cylinder. Good agreement is observed in most cases. The method is then applied to study the trajectory of one settling particle initially off‐centred between two parallel walls and the corresponding wake effects. Different particle trajectories related to particulate Reynolds numbers are presented and commented. A two‐body interaction problem is investigated too. This method allows the simulation of the transport of particles in a dilute suspension in reasonable time. One of the important features of this method is the computational cost that scales linearly with the number of particles. Copyright © 1999 John Wiley & Sons, Ltd.     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

An implicit three‐dimensional numerical model to simulate transport processes in coastal water bodies

A three‐dimensional baroclinic numerical model has been developed to compute water levels and water particle velocity distributions in coastal waters. The numerical model consists of hydrodynamic, transport and turbulence model components. In the hydrodynamic model component, the Navier–Stokes equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. The transport model component consists of the pollutant transport model and the water temperature and salinity transport models. In this component, the three‐dimensional convective diffusion equations are solved for each of the three quantities. In the turbulence model, a two‐equation k–ϵ formulation is solved to calculate the kinetic energy of the turbulence and its rate of dissipation, which provides the variable vertical turbulent eddy viscosity. Horizontal eddy viscosities can be simulated by the Smagorinsky algebraic sub grid scale turbulence model. The solution method is a composite finite difference–finite element method. In the horizontal plane, finite difference approximations, and in the vertical plane, finite element shape functions are used. The governing equations are solved implicitly in the Cartesian co‐ordinate system. The horizontal mesh sizes can be variable. To increase the vertical resolution, grid clustering can be applied. In the treatment of coastal land boundaries, the flooding and drying processes can be considered. The developed numerical model predictions are compared with the analytical solutions of the steady wind driven circulatory flow in a closed basin and of the uni‐nodal standing oscillation. Furthermore, model predictions are verified by the experiments performed on the wind driven turbulent flow of an homogeneous fluid and by the hydraulic model studies conducted on the forced flushing of marinas in enclosed seas. Copyright © 2000 John Wiley & Sons, Ltd.     

Simulation of gas–solid particle flows over a wide range of concentration

A two-fluid model of gas–solid particle flows that is valid for a wide range of the solid-phase volume concentration (dilute to dense) is presented. The governing equations of the fluid phase are obtained by volume averaging the Navier–Stokes equations for an incompressible fluid. The solid-phase macroscopic equations are derived using an approach that is based on the kinetic theory of dense gases. This approach accounts for particle–particle collisions. The model is implemented in a control-volume finite element method for simulations of the flows of interest in two-dimensional, planar or axisymmetric, domains. The chosen mathematical model and the proposed numerical method are applied to three test problems and one demonstration problem. © 1998 John Wiley & Sons, Ltd.     

Three‐dimensional remeshed smoothed particle hydrodynamics for the simulation of isotropic turbulence

We present a remeshed particle‐mesh method for the simulation of three‐dimensional compressible turbulent flow. The method is related to the meshfree smoothed particle hydrodynamics method, but the present method introduces a mesh for efficient calculation of the pressure gradient, and laminar and turbulent diffusion. In addition, the mesh is used to remesh (reorganise uniformly) the particles to ensure a regular particle distribution and convergence of the method. The accuracy of the presented methodology is tested for a number of benchmark problems involving two‐ and three‐dimensional Taylor‐Green flow, thin double shear layer, and three‐dimensional isotropic turbulence. Two models were implemented, direct numerical simulations, and Smagorinsky model. Taking advantage of the Lagrangian advection, and the finite difference efficiency, the method is capable of providing quality simulations while maintaining its robustness and versatility.     

Large eddy simulation and a simple wall model for turbulent flow calculation by a particle method

A turbulent channel flow and the flow around a cubic obstacle are calculated by the moving particle semi‐implicit method with the subparticle‐scale turbulent model and a wall model, which is based on the zero equation RANS (Reynolds Averaged Navier‐Stokes). The wall model is useful in practical problems that often involve high Reynolds numbers and wall turbulence, because it is difficult to keep high resolution in the near‐wall region in particle simulation. A turbulent channel flow is calculated by the present method to validate our wall model. The mean velocity distribution agrees with the log‐law velocity profile near the wall. Statistical values are also the same order and tendency as experimental results with emulating viscous layer by the wall model. We also investigated the influence of numerical oscillations on turbulence analysis in using the moving particle semi‐implicit method. Finally, the turbulent flow around a cubic obstacle is calculated by the present method to demonstrate capability of calculating practical turbulent flows. Three characteristic eddies appear in front of, over, and in the back of the cube both in our calculation and the experimental result that was obtained by Martinuzzi and Tropea. Mean velocity and turbulent intensity profiles are predicted in the same order and have similar tendency as the experimental result. Copyright © 2012 John Wiley & Sons, Ltd.     

Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Prediction of turbulent gas–solids flow in curved ducts using the Eulerian–Lagrangian method

The flow of particulate two‐phase flow mixtures occur in several components of solid fuel combustion systems, such as the pressurised fluidised bed combustors (PFBC) and suspension‐fired coal boilers. A detailed understanding of the mixture characteristics in the conveying component can aid in refining and optimising its design. In this study, the flow of an isothermal, dilute two‐phase particulate mixture has been examined in a high curvature duct, which can be representative of that transporting the gas–solid mixture from the hot clean‐up section to the gas turbine combustor in a PFBC plant. The numerical study has been approached by utilising the Eulerian–Lagrangian methodology for describing the characteristics of the fluid and particulate phases. By assuming that the mixture is dilute and the particles are spherical, the governing particle momentum equations have been solved with appropriately prescribed boundary conditions. Turbulence effects on the particle dispersion were represented by a statistical model that accounts for both the turbulent eddy lifetime and the particle transit time scales. For the turbulent flow condition examined it was observed that mixtures with small particle diameters had low interphase slip velocities and low impaction probability with the pipe walls. Increasing the particle diameters (>50 μm) resulted in higher interphase slip velocities and, as expected, their impaction probability with the pipe walls was significantly increased. The particle dispersion is significant for the smaller sizes, whereas the larger particles are relatively insensitive to the gas turbulence. The main particle impaction region, and locations most prone to erosion damage, is estimated to be within an outer duct length of two to six times the duct diameter, when the duct radius of curvature to the duct diameter ratio is equal to unity. Copyright © 1999 John Wiley & Sons, Ltd.     

An embedded boundary framework for compressible turbulent flow and fluid–structure computations on structured and unstructured grids

The finite volume method with exact two‐phase Riemann problems (FIVER) is a two‐faceted computational method for compressible multi‐material (fluid–fluid, fluid–structure, and multi‐fluid–structure) problems characterized by large density jumps, and/or highly nonlinear structural motions and deformations. For compressible multi‐phase flow problems, FIVER is a Godunov‐type discretization scheme characterized by the construction and solution at the material interfaces of local, exact, two‐phase Riemann problems. For compressible fluid–structure interaction (FSI) problems, it is an embedded boundary method for computational fluid dynamics (CFD) capable of handling large structural deformations and topological changes. Originally developed for inviscid multi‐material computations on nonbody‐fitted structured and unstructured grids, FIVER is extended in this paper to laminar and turbulent viscous flow and FSI problems. To this effect, it is equipped with carefully designed extrapolation schemes for populating the ghost fluid values needed for the construction, in the vicinity of the fluid–structure interface, of second‐order spatial approximations of the viscous fluxes and source terms associated with Reynolds averaged Navier–Stokes (RANS)‐based turbulence models and large eddy simulation (LES). Two support algorithms, which pertain to the application of any embedded boundary method for CFD to the robust, accurate, and fast solution of FSI problems, are also presented in this paper. The first one focuses on the fast computation of the time‐dependent distance to the wall because it is required by many RANS‐based turbulence models. The second algorithm addresses the robust and accurate computation of the flow‐induced forces and moments on embedded discrete surfaces, and their finite element representations when these surfaces are flexible. Equipped with these two auxiliary algorithms, the extension of FIVER to viscous flow and FSI problems is first verified with the LES of a turbulent flow past an immobile prolate spheroid, and the computation of a series of unsteady laminar flows past two counter‐rotating cylinders. Then, its potential for the solution of complex, turbulent, and flexible FSI problems is also demonstrated with the simulation, using the Spalart–Allmaras turbulence model, of the vertical tail buffeting of an F/A‐18 aircraft configuration and the comparison of the obtained numerical results with flight test data. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

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.