Development of high-order realizable finite-volume schemes for quadrature-based moment method

Kinetic equations containing terms for spatial transport, gravity, fluid drag and particle-particle collisions can be used to model dilute gas-particle flows. However, the enormity of independent variables makes direct numerical simulation of these equations almost impossible for practical problems. A viable alternative is to reformulate the problem in terms of moments of velocity distribution. Recently, a quadrature-based moment method was derived by Fox for approximating solutions to kinetic equation for arbitrary Knudsen number. Fox also described 1st- and 2nd-order finite-volume schemes for solving the equations. The success of the new method is based on a moment-inversion algorithm that is used to calculate non-negative weights and abscissas from moments. The moment-inversion algorithm does not work if the moments are non-realizable, meaning they do not correspond to a distribution function. Not all the finite-volume schemes lead to realizable moments. Desjardins et al. showed that realizability is guaranteed with the 1 st-order finite-volume scheme, but at the expense of excess numerical diffusion. In the present work, the nonrealizability of the standard 2 nd-order finite-volume scheme is demonstrated and a generalized idea for the development of high-order realizable finite-volume schemes for quadrature-based moment methods is presented. This marks a significant improvement in the accuracy of solutions using the quadrature-based moment method as the use of 1st-order scheme to guarantee realizability is no longer a limitation

A realizability-preserving high-order kinetic scheme using WENO reconstruction for entropy-based moment closures of linear kinetic equations in slab geometry

We develop a high-order kinetic scheme for entropy-based moment models of a one-dimensional linear kinetic equation in slab geometry. High-order spatial reconstructions are achieved using the weighted essentially non-oscillatory (WENO) method, and for time integration we use multi-step Runge-Kutta methods which are strong stability preserving and whose stages and steps can be written as convex combinations of forward Euler steps. We show that the moment vectors stay in the realizable set using these time integrators along with a maximum principle-based kinetic-level limiter, which simultaneously dampens spurious oscillations in the numerical solutions. We present numerical results both on a manufactured solution, where we perform convergence tests showing our scheme converges of the expected order up to the numerical noise from the numerical optimization, as well as on two standard benchmark problems, where we show some of the advantages of high-order solutions and the role of the key parameter in the limiter

ON THE IMPLEMENTATION OF MOMENT TRANSPORT EQUATIONS IN OPENFOAM TO PRESERVE CONSERVATION, BOUNDEDNESS AND REALIZABILITY

Different industrial scale multiphase systems can be successfully described by considering their polydispersity (e.g. particle/droplet/bubble size and velocity distributions) and phase coupling issues are properly overcome only by considering the evolution in space and time of such distributions, dictated by the so-called Generalized Population Balance Equation (GPBE). A computationally efficient approach for solving the GPBE is represented by the quadrature-based moment methods, where the evolution of the entire particle/droplet/bubble population is recovered by tracking some specific moments of the distribution and the quadrature approximation is used to solve the "closure problem" typical of moment-based methods. In this contribution some crucial computational and numerical details concerning the implementation of these methods into the opensource Computational Fluid Dynamics (CFD) code OpenFOAM are discussed. These aspects are in fact very often overlooked, resulting in implementations that do not satisfy the properties of conservation, realizability and boundedness. These constraints have to be satisfied in a consistent way, with respect to what done with the other conserved transported variables (e.g. volume fraction of the disperse phase) also when higher-order discretization schemes are used. These issues are illustrated on examples taken on our work on the simulation of fluid-fluid multiphase system

Strongly coupled fluid-particle flows in vertical channels. II. Turbulence modeling

In Part I, simulations of strongly coupled fluid-particle flow in a vertical channel were performed with the purpose of understanding, in general, the fundamental physics of wall-bounded multiphase turbulence and, in particular, the roles of the spatially correlated and uncorrelated components of the particle velocity.The exact Reynolds-averaged (RA) equations for high-mass-loading suspensions were presented, and the unclosed terms that are retained in the context of fully developed channel flow were evaluated in an Eulerian–Lagrangian (EL) framework. Here, data from the EL simulations are used to validate a multiphase Reynolds-stress model (RSM) that predicts the wall-normal distribution of the two-phase, one-point turbulence statistics up to second order. It is shown that the anisotropy of the Reynolds stresses both near the wall and far away is a crucial component for predicting the distribution of the RA particle-phase volume fraction. Moreover, the decomposition of the phase-average (PA) particle-phase fluctuating energy into the spatially correlated and uncorrelated components is necessary to account for the boundary conditions at the wall. When these factors are properly accounted for in the RSM, the agreement with the EL turbulence statistics is satisfactory at first order (e.g., PA velocities) but less so at second order (e.g., PA turbulent kinetic energy). Finally, an algebraic stress model for the PA particle-phase pressure tensor and the Reynolds stresses is derived from the RSM using the weak-equilibrium assumption