Simplified unified wave-particle method for diatomic gases based on Rykov model

TL;DR

This paper introduces a simplified unified wave-particle method for diatomic gases that effectively models multi-scale flows by integrating continuum and rarefied flow behaviors using the Rykov model.

Contribution

It develops a novel QMC mechanism and a unified wave-particle method that combines N-S and DSMC approaches for diatomic gases, addressing multi-scale flow challenges.

Findings

Validated across various flow regimes and geometries.

Accurately captures non-equilibrium phenomena.

Demonstrates effectiveness in hypersonic and micro-electro-mechanical systems flows.

Abstract

During the past decades, the numerical methods based on Navier-Stokes (N-S) equations and direct simulation Monte Carlo (DSMC) methods have been proved effective in simulating flows in the continuum and rarefied regimes, respectively. However, as single-scale methods, they face challenges in addressing common multi-scale problems, which are essential to simulate hypersonic flows around near-space vehicles and the flows in the micro-electro-mechanical systems. Hence, there is an urgent need for a method to predict multi-scale flows. In this work, a quantified model-competition (QMC) mechanism for diatomic multi-scale flows is derived from the integral solution of the Rykov model equations. This mechanism encapsulates both continuum and rarefied behaviors in a cell, weighted according to its local physical scale. By building upon the QMC mechanism, the N-S solver and DSMC solver are directly integrated within a cell to devise a simplified unified wave-particle (SUWP) method for diatomic gases. Specifically, the two-temperature equations considering the rotational energy are introduced into the kinetic inviscid flux (KIF) scheme and the N-S solver. As to the particle part, the collisionless DSMC solver is utilized to describe the non-equilibrium phenomenon. The proposed SUWP method for diatomic gases undergoes validation across a series of cases, including zero-dimensional homogeneous gas relaxation, one-dimensional normal shock structure, two-dimensional flow around the flat and the cylinder, and three-dimensional flows past the sphere and the blunt cone. Additionally, the implementation details of multi-scale wave-particle methods analysis and discussion are also undertaken in this work.