基于联合流体粒子方法的高速非平衡流模拟

安家琦; 刘晨宇; 赵文文; 樊国超; 杨俊沅; 陈伟芳

doi:10.7638/kqdlxxb-2026.0071

欢迎来到能源电力期刊网！

English Version 期刊登录

您当前的位置：

首页 >

文章列表页 >

基于联合流体粒子方法的高速非平衡流模拟

研究专栏——稀薄气体动力学理论与高效计算方法Ⅱ | 更新时间：2026-10-14

- 基于联合流体粒子方法的高速非平衡流模拟
- Simulation of high-speed non-equilibrium flows via the joint hydrodynamics-particle method
- 空气动力学学报 2026年44卷第8期页码：116-126
- 作者机构：
  
  1. 浙江大学航天航空学院,杭州,310027
  2. 北京宇航系统工程研究所,北京,100076
  3. 杭州电子科技大学机械工程学院,杭州,310018
- 作者简介：
  
  赵文文*，研究员，研究方向：稀薄气体动力学. E-mail：wwzhao@zju.edu.cn
- 基金信息：
  
  国家自然科学基金（92271204，U2570216，92371206）
- DOI：10.7638/kqdlxxb-2026.0071
  中图分类号： V211.3;O356
- 纸质出版：2026
- 稿件说明：
移动端阅览
安家琦, 刘晨宇, 赵文文, 等. 基于联合流体粒子方法的高速非平衡流模拟[J]. 空气动力学学报, 2026, 44(8): 116−126.

An Jiaqi, Liu Chenyu, Zhao Wenwen, et al. Simulation of high-speed non-equilibrium flows via the joint hydrodynamics-particle method[J]. 2026, 44(8): 116-126.
安家琦, 刘晨宇, 赵文文, 等. 基于联合流体粒子方法的高速非平衡流模拟[J]. 空气动力学学报, 2026, 44(8): 116−126. DOI： 10.7638/kqdlxxb-2026.0071

An Jiaqi, Liu Chenyu, Zhao Wenwen, et al. Simulation of high-speed non-equilibrium flows via the joint hydrodynamics-particle method[J]. 2026, 44(8): 116-126. DOI： 10.7638/kqdlxxb-2026.0071.

摘要

针对高速流动中的显著非平衡特性，本文基于联合流体粒子（joint hydrodynamics-particle

JHP）方法，对典型高速流动开展数值模拟研究。该方法采用随机粒子描述无碰撞迁移过程，并通过积分解框架结合单元平均碰撞时间和全局推进步长引入竞争机制，实现宏观Navier-Stokes（N-S）方程与介观粒子输运的耦合，能够准确刻画局部强非平衡流动特征。为验证该方法对高速非平衡流动的模拟能力，选取马赫数20条件下的一维激波结构和二维高速圆柱绕流，以及马赫数5条件下的钝楔绕流作为典型算例。通过与统一气体动理学格式（UGKS）及统一气体动理学波粒法（UGKWP）结果对比，验证了JHP方法在温度、速度等宏观量分布及流场结构预测方面的准确性。进一步与N-S方程结果对比，并引入局部Knudsen数对流场中非平衡区域的空间分布特征进行分析，发现相较于传统宏观方程，JHP能够更准确地捕捉激波层、近壁区及尾迹区等局部强非平衡区域的流动特征。研究表明，高速流场中的非平衡效应具有显著的区域性特征，并对流场演化产生重要影响。在钝楔绕流算例中，JHP方法的计算时间约为参考方法（IUGKS）的37.9%，内存消耗约为2.5%，表明该方法可为高速非平衡流的高效高精度模拟提供有效手段。

Abstract

Targeting the pronounced non-equilibrium characteristics of high-speed flows

this study conducted numerical simulations of typical high-speed flows based on the joint hydrodynamics–particle (JHP) method. The method employed stochastic particles to describe collisionless transport

and a competition mechanism was introduced within an integral solution framework through the coupling of the cell-averaged collision time and the global time step. This enabled a unified coupling between the macroscopic Navier–Stokes (N-S) equations and mesoscopic particle transport

allowing accurate resolution of locally strong non-equilibrium features. To validate the applicability of the method for high-speed non-equilibrium flows

three benchmark cases were selected: a one-dimensional shock structure and a high-speed flow over a circular cylinder at Mach 20

as well as a high-speed blunt-wedge flow at Mach 5. Comparisons with the unified gas-kinetic scheme (UGKS) and the unified gas-kinetic wave–particle (UGKWP) method demonstrated that the JHP method achieves high accuracy in predicting macroscopic quantities

including temperature and velocity distributions

as well as overall flow structures. Further comparisons with N-S solutions

combined with local Knudsen number analysis

revealed the spatial distribution of non-equilibrium regions within the flow field. Compared with conventional macroscopic approaches

the JHP method more accurately captures flow features in regions with strong non-equilibrium effects

such as the shock layer

near-wall region

and wake. These results indicate that non-equilibrium effects in high-speed flows exhibit pronounced spatial heterogeneity and significantly influence flow evolution. In the blunt wedge flow case

the computational time of the JHP method is approximately 37.9% of that of the reference method (IUGKS)

and its memory consumption is about 2.5%

indicating that the method provides an efficient and high-fidelity approach for the simulation of hypersonic non-equilibrium flows.

关键词

Keywords

references

Xu K. A unified computational fluid dynamics framework from rarefied to continuum regimes[M]. Cambridge University Press, 2021.

Wang M R, Li Z X. Simulations for gas flows in microgeometries using the direct simulation Monte Carlo method[J]. International Journal of Heat and Fluid Flow, 2004, 25(6): 975−985. doi: 10.1016/j.ijheatfluidflow.2004.02.024

Zhao W W, Chen W F, Agarwal R K. Formulation of a new set of Simplified Conventional Burnett equations for computation of rarefied hypersonic flows[J]. Aerospace Science and Technology, 2014, 38: 64−75. doi: 10.1016/j.ast.2014.07.014

Yadav U, Jonnalagadda A, Agrawal A. Derivation of extended-OBurnett and super-OBurnett equations and their analytical solution to plane Poiseuille flow at non-zero Knudsen number[J]. Journal of Fluid Mechanics, 2024, 983: A29. doi: 10.1017/jfm.2024.141

Yang J Y, Zeng S H, Zhao W W, et al. Numerical assessment of double lateral jets interaction in rarefied nonequilibrium crossflows via nonlinear coupled constitutive relations[J ] . Aerospace Science and Technology, 2025, 157: 109851. doi: 10.1016/j.ast. 2024.109851

Zeng S H, Yang J Y, Jiang Z Z, et al. Numerical assessment and application of simplified generalized hydrodynamic model in hypersonic near-continuum flows[J]. AIAA Journal, 2024, 62(6): 2370−2376. doi: 10.2514/1.j063031

Zeng S H, Yang J Y, Zhao W W, et al. Thermal non-equilibrium nonlinear coupled constitutive relations for hypersonic rarefied diatomic gas flows[J]. Aerospace Science and Technology, 2025, 161: 110147. doi: 10.1016/j.ast.2025.110147

Zeng S H, Yang J Y, Zhao W W, et al. Effect of bulk viscosity and relaxation model on nonequilibrium shock structure of diatomic gases[J]. Physical Review Fluids, 2025, 10(6): 063401. doi: 10.1103/8mtq-n1dn

Zeng S H, Yang J Y, Zhao W W, et al. Nonlinear coupled constitutive relations for hypersonic reacting flows with thermal nonequilibrium effect[J]. Physics of Fluids, 2025, 37: 016116. doi: 10.1063/5.0249391

Bird G A. Approach to translational equilibrium in a rigid sphere gas[J]. The Physics of Fluids, 1963, 6(10): 1518−1519. doi: 10.1063/1.1710976

Bird G A. Molecular gas dynamics and the direct simulation of gas flows[M]. Oxford: Oxford University Press, 1994. doi: 10.1093/oso/9780198561958.001.0001

Mieussens L. Discrete-velocity models and numerical schemes for the Boltzmann-BGK equation in plane and axisymmetric geometries[J]. Journal of Computational Physics, 2000, 162(2): 429−466. doi: 10.1006/jcph.2000.6548

Xu K. Direct modeling for computational fluid dynamics: construction and application of unified gas-kinetic schemes[M]. Hackensack, NJ: World Scientific, 2015.

Liu C, Xu K, Sun Q H, et al. A unified gas-kinetic scheme for continuum and rarefied flows IV: full Boltzmann and model equations[J]. Journal of Computational Physics, 2016, 314: 305−340. doi: 10.1016/j.jcp.2016.03.014

Wei Y F, Cao J Z, Xu K. Unified gas-kinetic scheme for reactive flow with multi-scale transport and chemical non-equilibrium[J]. Journal of Computational Physics, 2026, 546: 114514. doi: 10.1016/j.jcp.2025.114514

Tao S, Wang W H, Lin Y S, et al. Discrete unified gas kinetic scheme modeling for complex thermal flows using hierarchical Cartesian grids[J]. Engineering Analysis with Boundary Elements, 2025, 179: 106443. doi: 10.1016/j.enganabound.2025.106443

Liu C, Zhu Y J, Xu K. Unified gas-kinetic wave-particle methods I: continuum and rarefied gas flow[J]. Journal of Computational Physics, 2020, 401: 108977. doi: 10.1016/j.jcp.2019.108977

Liu Hongyu, Yang Xiaojian, Zhang Chuang, et al. UGKWP and IUGKP methods for multi-scale phonon transport with dispersion and polarization[PP]. arXiv, 2025. doi: 10.48550/arXiv.2506.16203

Cao J Z, Wei Y F, Long W P, et al. Adaptive criterion and modification of wave-particle decomposition in UGKWP method for high-speed flow simulation[J]. Computers Fluids, 2026, 305: 106896. doi: 10.1016/j.compfluid.2025.106896

Pu Z G, Xu K. Unified gas-kinetic wave-particle method for multiscale flow simulation of partially ionized plasma[J]. Journal of Computational Physics, 2025, 530: 113918. doi: 10.1016/j.jcp.2025.113918

Liu S, Zhong C W, Fang M. Simplified unified wave-particle method with quantified model-competition mechanism for numerical calculation of multiscale flows[J]. Physical Review E, 2020, 102: 013304. doi: 10.1103/physreve.102.013304

Yang S R, Liu S, Cao J Z, et al. Simplified unified wave-particle method for diatomic gases based on Rykov model[J]. Acta Mechanica Sinica, 2025, 42(2): 724678. doi: 10.1007/s10409-024-24678-x

Fan G C, Zhao W W, Yao S B, et al. The implementation of the three-dimensional unified gas-kinetic wave-particle method on multiple graphics processing units[J]. Physics of Fluids, 2023, 35(8): 086108. doi: 10.1063/5.0166092

Fan G C, Zhao W W, Jiang Z Z, et al. Study of shock structures using the unified gas-kinetic wave-particle method with various BGK models[J]. International Journal of Computational Fluid Dynamics, 2022, 36(1): 44−62. doi: 10.1080/10618562.2022.2082418

Liu W, Shu C, Teo C J, et al. A simple hydrodynamic-particle method for supersonic rarefied flows[J]. Physics of Fluids, 2022, 34(5): 057101. doi: 10.1063/5.0088946

Yang S R, Zhong C W, Ding N C, et al. A simplified unified wave-particle method for diatomic gases with rotational and vibrational non-equilibrium[J]. Aerospace Science and Technology, 2026, 168: 111261. doi: 10.1016/j.ast.2025.111261

Fei F, Zhang J, Li J, et al. A unified stochastic particle Bhatnagar-Gross-Krook method for multiscale gas flows[J]. Journal of Computational Physics, 2020, 400: 108972. doi: 10.1016/j.jcp.2019.108972

Feng K K, Tian P, Zhang J, et al. SPARTACUS: an open-source unified stochastic particle solver for the simulation of multiscale nonequilibrium gas flows[J]. Computer Physics Communications, 2023, 284: 108607. doi: 10.1016/j.cpc.2022.108607

An J Q, Fan G C, Zhao W W, et al. A joint hydrodynamics-particle method for efficient hypersonic simulation in all flow regimes[J]. Acta Mechanica Sinica, 2026, 42(5): 325150. doi: 10.1007/s10409-025-25150-x

Yang L M, Shu C, Yang W M, et al. An improved discrete velocity method (DVM) for efficient simulation of flows in all flow regimes[J]. Physics of Fluids, 2018, 30(6): 062005. doi: 10.1063/1.5039479

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

基于多尺度离散速度方法的高速飞行器跨流域流动数值模拟

基于修正N-S方程本构关系的气体动理论耦合方法

大型货机HWB布局气动性能与稳定性综合设计方法

航天飞行器热防护模拟感应耦合等离子加热器研究进展

考虑测量干扰抑制的风洞天平拓扑优化方法