Bounds-Preserving Lax-Wendroff Discontinuous Galerkin Schemes for Quadrature-Based Moment-Closure Approximations of Kinetic Models

by   Erica R. Johnson, et al.

The quadrature-based method of moments (QMOM) offers a promising class of approximation techniques for reducing kinetic equations to fluid equations that are valid beyond thermodynamic equilibrium. A major challenge with these and other closures is that whenever the flux function must be evaluated (e.g., in a numerical update), a moment-inversion problem must be solved that computes the flux from the known input moments. In this work we study a particular five-moment variant of QMOM known as HyQMOM and establish that this system is moment-invertible over a convex region in solution space. We then develop a high-order Lax-Wendroff discontinuous Galerkin scheme for solving the resulting fluid system. The scheme is based on a predictor-corrector approach, where the prediction step is a localized space-time discontinuous Galerkin scheme. The nonlinear algebraic system that arises in this prediction step is solved using a Picard iteration. The correction step is a straightforward explicit update using the predicted solution in order to evaluate space-time flux integrals. In the absence of additional limiters, the proposed high-order scheme does not in general guarantee that the numerical solution remains in the convex set over which HyQMOM is moment-invertible. To overcome this challenge, we introduce novel limiters that rigorously guarantee that the computed solution does not leave the convex set over which moment-invertible and hyperbolicity of the fluid system is guaranteed. We develop positivity-preserving limiters in both the prediction and correction steps, as well as an oscillation-limiter that damps unphysical oscillations near shocks and rarefactions. Finally, we perform convergence tests to verify the order of accuracy of the scheme, as well as test the scheme on Riemann data to demonstrate the shock-capturing and robustness of the method.



There are no comments yet.


page 1

page 2

page 3

page 4


Entropy stable discontinuous Galerkin methods for ten-moment Gaussian closure equations

In this article, we propose high order discontinuous Galerkin entropy st...

High order asymptotic preserving discontinuous Galerkin methods for gray radiative transfer equations

In this paper, we will develop a class of high order asymptotic preservi...

Energy stable Runge-Kutta discontinuous Galerkin schemes for fourth order gradient flows

We present unconditionally energy stable Runge-Kutta (RK) discontinuous ...

A new coordinate-transformed discretization method for minimum entropy moment approximations of linear kinetic equations

In this contribution we derive and analyze a new numerical method for ki...

A Discontinuous Galerkin method for Shock Capturing using a mixed high-order and sub-grid low-order approximation space

This article considers a new discretization scheme for conservation laws...

On energy preserving high-order discretizations for nonlinear acoustics

This paper addresses the numerical solution of the Westervelt equation, ...

An Efficient ADER-DG Local Time Stepping Scheme for 3D HPC Simulation of Seismic Waves in Poroelastic Media

Many applications from geosciences require simulations of seismic waves ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.