A Hybrid Gas-Kinetic Scheme and Discrete Velocity Method for Continuum and Rarefied Flows

TL;DR

This paper introduces a hybrid gas-kinetic and discrete velocity method that effectively simulates both continuum and rarefied gas flows, improving accuracy and efficiency across regimes.

Contribution

The work develops a novel hybrid approach combining GKS and DVM, with adaptive strategies for efficient multi-regime flow simulation.

Findings

Achieves asymptotic preserving solutions in continuum limit

Captures free molecular flows in rarefied limit

Demonstrates high accuracy and efficiency in diverse test cases

Abstract

The gas-kinetic scheme (GKS) provides high computational efficiency and accuracy for continuum flow simulations but is unable to reliably capture rarefaction effects. In contrast, although the discrete velocity method (DVM) is better suited for rarefied flows, it exhibits reduced accuracy and slow convergence when applied to continuum regimes. To overcome these limitations, this work proposes a hybrid GKS-DVM method that integrates the strengths of both approaches. The hybrid approach balances the equilibrium distribution function in GKS with the upwind-reconstructed non-equilibrium distribution function in DVM through a numerical collision time. This balancing strategy ensures to recover Navier-Stokes solutions in the continuum limit (asymptotic preserving), while naturally capturing free molecular flows in the rarefied limit. Moreover, the introduction of a numerical collision time significantly enhances robustness in shock capturing for continuum flow applications. To further reduce computational cost of the hybrid approach, several adaptive strategies based on the local Knudsen number and Mach number have been proposed. The effectiveness and accuracy of the proposed hybrid method are systematically assessed through four representative test cases: a flat-plate boundary layer, a lid-driven cavity flow, shock structures, and flow past a semi-cylinder. The first case is subjected to continuum conditions, while the latter two span a broad range of Knudsen numbers. The results demonstrate that the proposed method achieves high solution accuracy and computational efficiency across both continuum and rarefied flow regimes.