Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

Physics and modelling of turbulent particle deposition and entrainment: Review of a systematic study

Deposition and entrainment of particles in turbulent flows are crucial in a number of technological applications and environmental processes. We present a review of recent results from our previous works, which led to physical insights on these phenomena. These results were obtained from a systematic numerical study based on the accurate resolution – Direct Numerical Simulation via a pseudo-spectral approach – of the turbulent flow field, and on Lagrangian tracking of particles under different modelling assumptions. We underline the multiscale aspect of wall turbulence, which has challenged scientists to devise simple theoretical models adequate to fit experimental data, and we show that a sound rendering of wall turbulence mechanisms is required to produce a physical understanding of particle deposition and re-entrainment. This physical understanding can be implemented in more applied simulation techniques, such as Large-Eddy Simulation. Our arguments are based also on the phenomenology of coherent structures and on the examination of flow topology in connection with particle preferential distribution. Starting from these concepts, reasons why theoretical predictions may fail are examined together with the requirements which must be fulfilled by suitable predictive models.     

Averaging of particle data from phase Doppler anemometry in unsteady two-phase flow: Validation by numerical simulation

A novel post-processing algorithm is proposed to correct statistical bias observed in the treatment of time series obtained by a phase Doppler anemometer (PDA) at flow locations with variable particle velocity and concentration. Extensive properties of each validated particle are weighted with their inverse measuring (validation) volume to account for the procedure of particle sampling and fluctuations in the particle concentration. To compensate for the short characteristic length of the validation volume, the properties of particles are expressed by properties of fields of particle groups, using a local averaging time. A window shift and a decorrelation scheme are applied on the fields to increase their frequency resolution. This algorithm has been tested on numerical time series, provided by an Eulerian/Lagrangian code representing a gas/solids flow past a bluff body. Moments and spectral estimates of concentration and velocity of particle groups were successfully validated by the numerical simulation using the PDA data algorithm and control volume averaging. The control volume was much larger than the PDA validation volume, but the centre positions of the two volumes were identical.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.