Implementation of CHyQMOM in OpenFOAM for the simulation of non-equilibrium gas-particle flows under one-way and two-way coupling

The modeling of dilute gas-particle flows is challenging due to particle-trajectory crossing (PTC). Lagrangian particle tracking can be used, but requires a large number of parcels resulting in high computational costs. A less costly method is the Eulerian number density function (NDF) approach, based on the Boltzmann equation, often solved in terms of lower-order moments of the NDF. In this context the conditional hyperbolic quadrature method of moments (CHyQMOM) was developed and is here implemented for the first time in the OpenFOAM-7, together with high-order advection schemes and a new operator splitting procedure. The resulting solver is used to simulate different test cases: phase segregation problems, collision-less and weakly-collisional PTC flows, asymmetric and symmetric Taylor-Green vortex flow and a dilute gas-particle riser. Results, validated against analytical solutions and predictions obtained with Lagrangian particle tracking, show that the implemented CHyQMOM can be used to handle highly non-equilibrium flows. (c) 2021 Elsevier B.V. All rights reserved.