Assessment of numerical methods for fully resolved simulations of particle-laden turbulent flows

J.C. Brändle de Motta (1, 2) , P. Costa (3, 4) , J.J. Derksen (5) , C. Peng (6) , L.P. Wang (6, 7) , W.P. Breugem (3) , Jean-Luc Estivalezes (1, 8, 9) , S. Vincent (10) , Éric Climent (1) , Pascal Fede (1) , P. Barbaresco (9) , N. Renon (9)

1 IMFT - Institut de mécanique des fluides de Toulouse
2 CORIA - Complexe de recherche interprofessionnel en aérothermochimie
3 TU Delft - Delft University of Technology
4 KTH - Department of Mechanics
5 University of Aberdeen
6 University of Delaware [Newark]
7 SUSTech - Southern University of Science and Technology
8 DMPE, ONERA, Université de Toulouse [Toulouse]
9 CALMIP - Calcul en Midi-Pyrénée
10 MSME - Laboratoire de Modélisation et Simulation Multi Echelle

During the last decade, many approaches for resolved-particle simulation (RPS) have been developed for numerical studies of finite size particle-laden turbulent flows. In this paper, three RPS approaches are compared for a particle-laden decaying turbulence case. These methods are, the Volume-of-Fluid Lagrangian method, based on the viscosity penalty method (VoF-Lag); a direct forcing Immersed Boundary Method, based on a regularized delta-function approach for the fluid/solid coupling (IBM); and the Bounce Back scheme developed for Lattice Boltzmann method (LBM-BB). The physics and the numerical performances of the methods are analyzed. Modulation of turbulence is observed for all the methods, with a faster decay of turbulent kinetic energy compared to the single-phase case. Lagrangian particle statistics, such as the velocity probability density function and the velocity autocorrelation function, show minor differences among the three methods. However, major differences between the codes are observed in the evolution of the particle kinetic energy. These differences are related to the treatment of the initial condition when the particles are inserted in an initially single-phase turbulence. The averaged particle/fluid slip velocity is also analyzed, showing similar behavior as compared to the results referred in the literature. The computational performances of the different methods differ significantly. The VoF-Lag method appears to be computationally most expensive. Indeed, this method is not adapted to turbulent cases. The IBM and LBM-BB implementations show very good scaling.

Mots clés: PARTICLE-LADEN FLOW; FINITE-SIZE PARTICLE; DIRECT NUMERICAL SIMULATION

Mots clés: ECOULEMENT DE PARTICULE CHARGE; TURBULENCE; SIMULATION NUMERIQUE DIRECTE

Domaines: Milieux fluides et réactifs; Mécanique des fluides [physics.class-ph]

J.C. Brändle de Motta, P. Costa, J.J. Derksen, C. Peng, L.P. Wang, et al.. Assessment of numerical methods for fully resolved simulations of particle-laden turbulent flows. Computers and Fluids, 2018, 179, pp.1-14. ⟨10.1016/j.compfluid.2018.10.016⟩. ⟨hal-02006358⟩