Dissipativeness of the hyperbolic quadrature method of moments for kinetic equations

Ruixi Zhang,Yihong Chen,Qian Huang,Wen-An Yong

2024-06-20

Abstract:This paper presents a dissipativeness analysis of a quadrature method of moments (called HyQMOM) for the one-dimensional BGK equation. The method has exhibited its good performance in numerous applications. However, its mathematical foundation has not been clarified. Here we present an analytical proof of the strict hyperbolicity of the HyQMOM-induced moment closure systems by introducing a polynomial-based closure technique. As a byproduct, a class of numerical schemes for the HyQMOM system is shown to be realizability preserving under CFL-type conditions. We also show that the system preserves the dissipative properties of the kinetic equation by verifying a certain structural stability condition. The proof uses a newly introduced affine invariance and the homogeneity of the HyQMOM and heavily relies on the theory of orthogonal polynomials associated with realizable moments, in particular, the moments of the standard normal distribution.

Numerical Analysis

What problem does this paper attempt to address?

The problems that this paper attempts to solve mainly focus on two aspects: one is to prove the strict hyperbolicity of the HyQMOM (Hyperbolic Quadrature Method of Moments) system based on the moment method, and the other is to verify whether the system preserves the dissipative properties of the original kinetic equations.

### 1. Proof of strict hyperbolicity

The paper mentions that although HyQMOM has shown good performance in many applications, its mathematical foundation has not been clear. For this reason, by introducing a polynomial - based closure technique, the authors give an analytical proof of the strict hyperbolicity of the moment - closure system induced by HyQMOM. Specifically, they prove that the characteristic polynomial can be decomposed as:

\[ F(X;M) = Q_n \cdot \left[ (X - a_n)Q_{n - 2} - \frac{2n+\gamma}{n} b_n Q_{n - 1} \right] := Q_n R_{n + 1} \]

where \( Q_n \) and \( R_{n + 1} \) are orthogonal polynomials. By analyzing the roots of these polynomials, it can be concluded that the characteristic polynomial has \( 2n + 1 \) distinct real roots, thus proving the strict hyperbolicity of the system.

### 2. Verification of dissipative properties

To ensure that the moment system induced by HyQMOM correctly inherits the dissipative characteristics of the original kinetic equations, the authors verify the structural stability condition. This condition is a suitable substitute for the H - theorem and is used to characterize the dissipative properties of the moment system. The specific steps include:

- **Condition (I)**: For any equilibrium state \( M\in E \), there exist invertible matrices \( P \) and \( \hat{T} \) such that \( P S_M(M)P^{-1} \) has a specific form.
- **Condition (III)**: The coupling between the coefficient matrix and the source term satisfies a certain inequality.

Through these verifications, the authors prove that the moment system induced by HyQMOM satisfies the structural stability condition, thus ensuring the dissipative properties of the system.

### Conclusion

The main contribution of the paper lies in providing theoretical proofs of the strict hyperbolicity and dissipative properties of the HyQMOM system, filling the gaps in existing research. This not only deepens the understanding of HyQMOM but also provides a solid theoretical basis for the application of this method in fields such as rarefied gas dynamics, multiphase flow, and plasma simulation.