Locally refined discrete velocity grids for stationary rarefied flow simulations

Céline Baranger,Jean Claudel,Nicolas Hérouard,Luc Mieussens
DOI: https://doi.org/10.1016/j.jcp.2013.10.014
2013-10-05
Abstract:Most of deterministic solvers for rarefied gas dynamics use discrete velocity (or discrete ordinate) approximations of the distribution function on a Cartesian grid. This grid must be sufficiently large and fine to describe the distribution functions at every space position in the computational domain. For 3-dimensional hypersonic flows, like in re-entry problems, this induces much too dense velocity grids that cannot be practically used, for memory storage requirements. In this article, we present an approach to generate automatically a locally refined velocity grid adapted to a given simulation. This grid contains much less points than a standard Cartesian grid and allows us to make realistic 3-dimensional simulations at a reduced cost, with a comparable accuracy.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is how to effectively reduce the number of grid points required for velocity - space discretization in simulating steady supersonic flows in rarefied gas dynamics, thereby reducing the computational cost (including CPU time and memory storage requirements), while maintaining an accuracy comparable to that of the traditional Cartesian grid. Specifically: 1. **The problem of velocity - space discretization for supersonic flows**: For three - dimensional supersonic flows (such as re - entry problems), the traditional Cartesian velocity grid is too dense, resulting in excessive computational resource requirements and being impractical for actual applications. 2. **The proposal of locally refined velocity grids**: To address the above challenges, the author proposes a method for automatically generating locally refined velocity grids. This grid is refined where needed (for example, in the core region of the distribution function) and coarsened in other places (for example, in the tail region of the distribution function). This can significantly reduce the number of grid points, thereby reducing the computational cost. 3. **Adaptive grid generation algorithm**: The author estimates the macroscopic parameters of the flow field (such as velocity and temperature) by using compressible Navier - Stokes pre - simulation and generates an adaptive velocity grid based on these parameters. This algorithm can automatically adjust the grid density according to the characteristics of the flow field, ensuring sufficient resolution in key areas. 4. **Numerical verification**: Through numerical experiments, the author demonstrates the effectiveness of this method. For example, in a re - entry flow simulation with a Mach number of 25, using the new method can reduce CPU time and memory requirements by as much as 24 times while maintaining an accuracy comparable to that of the standard Cartesian grid. In summary, this paper aims to solve the computational efficiency problem in supersonic flow simulation in rarefied gas dynamics by introducing locally refined velocity grids, making complex three - dimensional simulations possible.