高超声速非平衡流区域识别及网格自适应重构DSMC方法

李子玮; 王学德

doi:10.13224/j.cnki.jasp.20230379

高超声速非平衡流区域识别及网格自适应重构DSMC方法

doi: 10.13224/j.cnki.jasp.20230379

李子玮,
王学德^*, ,

南京理工大学能源与动力工程学院，南京 210094

详细信息

作者简介:
李子玮（1998-），男，硕士生，主要从事稀薄气体动力学方面研究

通讯作者:
王学德（1977-），男，副教授，博士，主要从事稀薄气体动力学方面研究。E-mail：wangxuede2000@njust.edu.cn

中图分类号: V448.2
计量
- 文章访问数: 576
- HTML浏览量: 131
- PDF量: 38
- 被引次数: 0
出版历程
- 收稿日期: 2023-06-08
- 网络出版日期: 2024-10-16

Hypersonic nonequilibrium flow region identification and adaptive grid reconfiguration DSMC method

LI Ziwei,
WANG Xuede^{*
, ,}

School of Energy and Power Engineering，Nanjing University of Science and Technology，Nanjing 210094，China

摘要

摘要:
采用分子动力学理论导出系统熵与熵产率的计算公式，识别高超声速稀薄流非平衡区域。根据熵产与非平衡现象的基本关系，确定了衡量气体非平衡态的判断方法。由于熵产率参数存在局限性，提出了一种熵产克努森数（Knudsen number）作为气体非平衡态衡量参数，与局部克努森数相结合作为非平衡态判据。基于两种非平衡判据相结合的策略，发展了一套高超声速流场非平衡区域识别与网格自适应重构的DSMC（direct simulation Monte Carlo）计算方法。对典型非平衡流场进行模拟，数值结果表明：所提出的识别参数能精准识别全流场中的非平衡效应，识别率达到了99%以上。在此基础上重构非平衡区域网格，计算得到的流场宏观量相对误差小于3%，壁面气动参数相对误差小于5%，证明了自适应重构方法的有效性。
- 高超声速 /
- 非平衡流动 /
- DSMC方法 /
- 熵产率 /
- 自适应网格
Abstract:
By using the molecular dynamics theory, the calculation formulas for system entropy and entropy production rate were derived to identify non-equilibrium regions in hypersonic rarefied flows. A method for measuring the non-equilibrium state of gases was determined based on the fundamental relationship between entropy production and non-equilibrium phenomena. Due to the limitations of entropy production rate parameters, a gas non-equilibrium measure parameter called the entropy production Knudsen number was proposed, which was combined with the local Knudsen number as a criterion for non-equilibrium states. By integrating two non-equilibrium criteria, a direct simulation Monte Carlo (DSMC) computational method was developed for identifying non-equilibrium regions and performing grid adaptive reconstruction in hypersonic flow fields. Typical non-equilibrium flow fields were simulated, the results showed that the proposed identification parameters can accurately identify non-equilibrium effects throughout the flow field, with an identification rate of over 99%. On the basis of this strategy, the grid in the non-equilibrium region was reconstructed, and the relative errors of the macroscopic flow field quantities were less than 3%, while the relative errors of wall aerodynamic parameters were less than 5%, which demonstrated the effectiveness of adaptive reconstruction method.
- hypersonic /
- nonequilibrium flow /
- DSMC method /
- entropy generation rate /
- adaptive grid

HTML

图 1 网格重构示意图

Figure 1. Schematic diagram of grid reconstruction

下载: 全尺寸图片幻灯片

图 2 双锥构型的外形（单位：mm）

Figure 2. Shape diagram of double-cone configuration （unit：mm）

下载: 全尺寸图片幻灯片

图 3 双锥构型非平衡区与网格重构

Figure 3. Non-equilibrium region and grid reconstruction of double-cone configuration

下载: 全尺寸图片幻灯片

图 4 主流方向速度的流线图

Figure 4. Streamline diagram of velocity in mainstream direction

下载: 全尺寸图片幻灯片

图 5 双锥的壁面气动特性分布曲线

Figure 5. Distribution curve of wall aerodynamic characteristics for double-cone

下载: 全尺寸图片幻灯片

图 6 圆柱构型非平衡区与网格重构

Figure 6. Non-equilibrium region and grid reconstruction of cylinder configuration

下载: 全尺寸图片幻灯片

图 7 圆柱绕流的密度与温度云图

Figure 7. Density and temperature contours of cylinder

下载: 全尺寸图片幻灯片

图 8 圆柱的壁面气动特性分布曲线

Figure 8. Distribution curve of wall aerodynamic characteristics for cylinder

下载: 全尺寸图片幻灯片

参考文献(19)

[1]	黄飞,张亮,程晓丽,等. 连续流失效对近空间飞行器气动特性的影响[J]. 空气动力学学报,2013,31(5): 623-628,640. HUANG Fei,ZHANG Liang,CHENG Xiaoli,et al. Effects of continuum breakdown on aerodynamics of near space vehicle[J]. Acta Aerodynamica Sinica,2013,31(5): 623-628,640. (in Chinese doi: 10.7638/kqdlxxb-2012.0002 HUANG Fei, ZHANG Liang, CHENG Xiaoli, et al. Effects of continuum breakdown on aerodynamics of near space vehicle[J]. Acta Aerodynamica Sinica, 2013, 31(5): 623-628, 640. (in Chinese) doi: 10.7638/kqdlxxb-2012.0002
[2]	王智慧,鲍麟,童秉纲. 尖化前缘的稀薄气体化学非平衡流动和气动加热相似律研究[J]. 气体物理,2016,1(1): 5-12. WANG Zhihui,BAO Lin,TONG Binggang. Similarity law of aero-heating to sharpened noses in rarefied chemical nonequilibrium flows[J]. Physics of Gases,2016,1(1): 5-12. (in Chinese WANG Zhihui, BAO Lin, TONG Binggang. Similarity law of aero-heating to sharpened noses in rarefied chemical nonequilibrium flows[J]. Physics of Gases, 2016, 1(1): 5-12. (in Chinese)
[3]	SCHROCK C,MCMULLAN R,CAMBEROS J. Calculation of entropy generation rates via DSMC with application to continuum/equilibrium onset[R]. AIAA2005-4830,2005.
[4]	CHIGULLAPALLI S,VENKATTRAMAN A,ALEXEENKO A,et al. Non-equilibrium flow modeling using high-order schemes for the boltzmann model equations[R]. AIAA2008-3929,2008.
[5]	BIRD G A. Breakdown of translational and rotational equilibrium in gaseous expansions[J]. AIAA Journal,1970,8(11): 1998-2003. doi: 10.2514/3.6037
[6]	WANG W L,BOYD I. Continuum breakdown in hypersonic viscous flows[R]. AIAA2002-651,2002.
[7]	BOYD I D,CHEN Gang,CANDLER G V. Predicting failure of the continuum fluid equations in transitional hypersonic flows[J]. Physics of Fluids,1995,7(1): 210-219. doi: 10.1063/1.868720
[8]	GARCIA A L,BELL J B,CRUTCHFIELD W Y,et al. Adaptive mesh and algorithm refinement using direct simulation Monte Carlo[J]. Journal of Computational Physics,1999,154(1): 134-155. doi: 10.1006/jcph.1999.6305
[9]	唐振宇,蔡国飙. 两种用于NS-DSMC耦合方法的连续失效参数对比[J]. 北京航空航天大学学报,2014,40(3): 389-393. TANG Zhenyu,CAI Guobiao. Comparison analysis between two continuum breakdown parameters applied in coupled NS-DSMC method[J]. Journal of Beijing University of Aeronautics and Astronautics,2014,40(3): 389-393. (in Chinese TANG Zhenyu, CAI Guobiao. Comparison analysis between two continuum breakdown parameters applied in coupled NS-DSMC method[J]. Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(3): 389-393. (in Chinese)
[10]	陈杰,赵磊. 高超声速流存在局部稀薄效应的一个判据及相应的流动特性[J]. 空气动力学学报,2018,36(1): 4-11. CHEN Jie,ZHAO Lei. A criteron for the existence of local rarefaction effect in a hypersonic flow field and the corresponding flow characteristics[J]. Acta Aerodynamica Sinica,2018,36(1): 4-11. (in Chinese doi: 10.7638/kqdlxxb-2017.0179 CHEN Jie, ZHAO Lei. A criteron for the existence of local rarefaction effect in a hypersonic flow field and the corresponding flow characteristics[J]. Acta Aerodynamica Sinica, 2018, 36(1): 4-11. (in Chinese) doi: 10.7638/kqdlxxb-2017.0179
[11]	CAMBEROS J,CHEN Poheng. Continuum breakdown parameter based on entropy generation rates[R]. AIAA2003-157,2003.
[12]	SCHROCK C,MCMULLAN R,CAMBEROS J. Continuum onset parameter based on entropy gradients using boltzmann’s H-theorem[R]. AIAA2005-967,2005.
[13]	EJTEHADI O,ESFAHANI J A,ROOHI E. Compressibility and rarefaction effects on entropy and entropy generationin micro/nano Couette flow using DSMC[J]. Journal of Physics: Conference Series,2012,362: 012008. doi: 10.1088/1742-6596/362/1/012008
[14]	CAMBEROS J. The production of entropy in relation to numerical error in compressible viscous fluid flow[R]. AIAA2000-2333,2000.
[15]	CAMBEROS J,SCHROCK C,MCMULLAN R,et al. Development of continuum onset criteria with direct simulation monte-carlo using boltzmann’s H-theorem: review and vision[R]. AIAA2006-2942,2006.
[16]	CARR R,BRANAM R,CAMBEROS J. Quantifying Non-Equilibrium using Entropy Generation[R]. AIAA2007-1127,2007.
[17]	BLAZEK J. Spatial discretisation[M]//Computational fluid dynamics: principles and applications. Amsterdam,Netherlands: Elsevier,2001: 129-179.
[18]	MOSS J,BIRD G. DSMC simulations of hypersonic flows with shock interactions and validation with experiments[J]. Portland,US: AIAA 37th Thermophysics Conference,2004.
[19]	谢东. 高超声速火星大气反应流场的DSMC方法研究及应用[D]. 南京: 南京理工大学,2023. XIE Dong. Research and application of DSMC method for hypersonic Mars atmospheric reaction flow field[D]. Nanjing: Nanjing University of Science and Technology,2023. (in Chinese XIE Dong. Research and application of DSMC method for hypersonic Mars atmospheric reaction flow field[D]. Nanjing: Nanjing University of Science and Technology, 2023. (in Chinese)