An efficient solution algorithm for force-driven continuum and rarefied flows

TL;DR

This paper introduces a novel, highly efficient numerical algorithm based on the Boltzmann-BGK equation that accurately models multi-scale, force-driven gas flows across continuum and rarefied regimes, significantly improving computational speed.

Contribution

A new solution algorithm combining Hermite spectral methods with a multi-scale model for force-driven flows, achieving up to 50 times faster computation with maintained accuracy.

Findings

Significantly improved computational efficiency (up to 50x)

Accurate modeling of multi-scale, force-driven flows

Preservation of key gas-dynamic features in simulations

Abstract

Gaseous flows under an external force are intrinsically defined by their multi-scale nature due to the large variation of densities along the forcing direction. Devising a numerical method capable of accurately and efficiently solving force-driven cross-scale flow dynamics, encompassing both continuum and rarefied regimes, continues to pose a formidable and enduring challenge. In this work, a novel solution algorithm for multi-scale and non-equilibrium flow transport under an external force is developed based on the Boltzmann-BGK equation. The core innovation lies in the fusion of the Hermite spectral method (employed to characterize non-equilibrium particle distributions) with a multi-scale evolution model (sourced from the unified gas-kinetic scheme), achieving a seamless connection between computational methods and physical models. To accommodate the properties of the spectral-collocation method, a series of collocation points and weights are adapted based on the Gauss-Hermite quadrature. As a result, the computational efficiency of the solution algorithm is significantly improved (up to 50 times) while maintaining comparable accuracy as the classical discrete velocity method. It is demonstrated that the solution algorithm effectively preserves the key structural features of gas-dynamic systems subjected to an external force, e.g., the well-balanced property. Extensive numerical experiments have been performed to verify the accuracy and efficiency of the proposed method, including the one-dimensional hydrostatic equilibrium problem, the Sod shock tube, the Fourier flow, the Poiseuille flow, and the Rayleigh-Taylor instability problem. The proposed methodology can provide substantive theoretical insights into a wide range of engineering challenges involving force-driven multi-scale flows.