On the Anisotropic Gaussian velocity closure for inertial-particle laden flows

Aymeric Vié (1, 2, 3) , François Doisneau (1, 3, 4) , Marc Massot (1, 3, 5)

1 EM2C - Laboratoire d'Énergétique Moléculaire et Macroscopique, Combustion
2 CTR - Center for Turbulence Research [Stanford]
3 Ecole Centrale Paris
4 ONERA - The French Aerospace Lab [Palaiseau]
5 FR3487 - Fédération de Mathématiques de l'Ecole Centrale Paris

The accurate simulation of disperse two-phase flows, where a discrete partic- ulate condensed phase is transported by a carrier gas, is crucial for many applications; Eulerian approaches are well suited for high performance computations of such flows. However when the particles from the disperse phase have a significant inertia compared to the time scales of the flow, particle trajectory crossing (PTC) occurs i.e. the particle ve- locity distribution at a given location can become multi-valued. To properly account for such a phenomenon many Eulerian moment methods have been recently proposed in the literature. The resulting models hardly comply with a full set of desired criteria involving: 1- ability to reproduce the physics of PTC, at least for a given range of particle inertia, 2- well-posedness of the resulting set of PDEs on the chosen moments as well as guaran- teed realizability, 3- capability of the model to be associated with a high order realizable numerical scheme for the accurate resolution of particle segregation in turbulent flows. The purpose of the present contribution is to introduce a multi-variate Anisotropic Gaus- sian closure for such particulate flows, in the spirit of the closure that has been suggested for out-of-equilibrium gas dynamics and which satisfies the three criteria. The novelty of the contribution is three-fold. First we derive the related moment systems of conservation laws with source terms, and justify the use of such a model in the context of high Knudsen numbers, where collision operators play no role. We exhibit the main features and advan- tages in terms of mathematical structure and realizability. Then a second order accurate and realizable MUSCL/HLL scheme is proposed and validated. Finally the behavior of the method for the description of PTC is thoroughly investigated and its ability to account accurately for inertial particulate flow dynamics in typical configurations is assessed.

Mots clés: disperse phase flows; numerical methods; particle trajectory crossing; entropy maximization

Domaines: Mécanique des fluides [physics.class-ph]; Mécanique des fluides [physics.class-ph]; Milieux fluides et réactifs; Analyse numérique [math.NA]

Aymeric Vié, François Doisneau, Marc Massot. On the Anisotropic Gaussian velocity closure for inertial-particle laden flows. Communications in Computational Physics, 2015, 17 (1), pp.1-46. ⟨10.4208/cicp.021213.140514a⟩. ⟨hal-00912319v2⟩