Abstract:The numerical method of dynamical low-rank approximation (DLRA) has recently been applied to various kinetic equations showing a significant reduction of the computational effort. In this paper, we apply this concept to the linear Boltzmann-Bhatnagar-Gross-Krook (Boltzmann-BGK) equation which due its high dimensionality is challenging to solve. Inspired by the special structure of the non-linear Boltzmann-BGK problem, we consider a multiplicative splitting of the distribution function. We propose a rank-adaptive DLRA scheme making use of the basis update & Galerkin integrator and combine it with an additional basis augmentation to ensure numerical stability, for which an analytical proof is given and a classical hyperbolic Courant-Friedrichs-Lewy (CFL) condition is derived. This allows for a further acceleration of computational times and a better understanding of the underlying problem in finding a suitable discretization of the system. Numerical results of a series of different test examples confirm the accuracy and efficiency of the proposed method compared to the numerical solution of the full system.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to provide a stable dynamic low - rank discretization method for the linear Boltzmann - BGK equation. Specifically, the author focuses on how to effectively and stably solve the high - dimensional linear Boltzmann - BGK equation through the dynamic low - rank approximation (DLRA) technique, thereby significantly reducing the computational complexity and memory requirements.
### Main problem decomposition
1. **Computational complexity of high - dimensional problems**:
The linear Boltzmann - BGK equation requires huge computational resources and memory when numerically solved due to its high - dimensional characteristics (including time, space, and velocity variables). Traditional full - order methods are inefficient in dealing with such high - dimensional problems and are difficult to meet the needs of practical applications.
2. **Numerical stability**:
When applying dynamic low - rank approximation, ensuring the stability of the numerical scheme is a key issue. An inappropriate discretization method may lead to numerical instability, which in turn affects the accuracy and reliability of the results. Therefore, it is crucial to design a numerical scheme that is both efficient and stable.
3. **Multiscale problems**:
The Boltzmann - BGK equation involves multiple physical scales. Especially in the case of frequent collisions, the solution is close to the Maxwellian equilibrium distribution, which brings additional challenges to numerical simulation. How to effectively handle this multiscale characteristic is one of the research focuses.
### Solutions
To address the above challenges, the author proposes a dynamic low - rank discretization method based on multiplicative splitting. Specific measures include:
- **Multiplicative splitting**: Decompose the distribution function \( f(t, x, v) \) into \( f = M g \), where \( M \) is the Maxwellian equilibrium distribution and \( g \) is the part deviating from equilibrium. This method makes \( g \) have a lower rank, thus facilitating the application of dynamic low - rank approximation.
- **Stable numerical scheme**: Ensure the stability of the numerical scheme by introducing the classical hyperbolic CFL condition. The paper gives a strict mathematical proof and derives the corresponding CFL condition, ensuring that the choice of time step can effectively reduce the amount of computation.
- **Adaptive rank update**: Adopt the rank - adaptive basis update & Galerkin (BUG) integrator to dynamically adjust the rank size, avoiding the limitations of fixing the rank in advance and improving the flexibility and efficiency of the algorithm.
### Summary
By introducing multiplicative splitting and dynamic low - rank approximation techniques, combined with strict stability analysis and adaptive rank update strategies, this paper successfully solves the problem of efficiently and stably solving the linear Boltzmann - BGK equation in high - dimensional cases. This method not only reduces the computational complexity and memory requirements but also provides new ideas and tools for dealing with similar high - dimensional kinetic equations.