Assessment of numerical methods for fully resolved simulations of particle-laden turbulent flows. (30th January 2019)

• Record Type: Journal Article
• Authors: Brändle de Motta, J.C.; Costa, P.; Derksen, J.J.; Peng, C.; Wang, L.-P.; Breugem, W.-P.; Estivalezes, J.L.; Vincent, S.; Climent, E.; Fede, P.; Barbaresco, P.; Renon, N.
• Is Part Of: Computers & fluids. Volume 179(2019)
• Page Start: 1
• Page End: 14
• Publication Date: 2019-01-30
• Subjects: Particle-laden flows -- Finite-size particles -- Turbulence -- Direct numerical simulations; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2018.10.016
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: Three different approaches for resolved-particle simulation are compared for a particle-laden decaying turbulence case. Turbulence modulation, particle statistics and averaged particle/fluid slip velocity are analyzed and show similar behavior among these approaches. The computational performances of the different approaches are also computed and differ significantly. Abstract: 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.