硬球-拟颗粒模拟高超声速稀薄气体流动

doi:10.12034/j.issn.1009-606X.219120

过程工程学报 ›› 2019, Vol. 19 ›› Issue (6): 1093-1100.DOI: 10.12034/j.issn.1009-606X.219120

硬球-拟颗粒模拟高超声速稀薄气体流动

赵祺^1,2，赵明璨^1,2，马琳博^1,3，葛蔚^1,2*

1. 中国科学院过程工程研究所多相复杂系统国家重点实验室，北京 100190 2. 中国科学院大学化学工程学院，北京 100049 3. 中国科学院大学化学科学学院，北京 100049

收稿日期:2019-01-23 修回日期:2019-04-15 出版日期:2019-12-22 发布日期:2019-12-22
通讯作者: 葛蔚 wge@ipe.ac.cn
基金资助:
国防基础科研科学挑战专题项目;高芳烃高含氮重油催化转化反应基础研究;中国科学院前沿科学重点研究项目;中国科学院信息化专项课题

Hard-sphere/pseudo-particle modeling (HS-PPM) for hypersonic rarefied gas flow

Qi ZHAO^1,2, Mingcan ZHAO^1,2, Linbo MA^1,3, Wei GE^1,2*

1. State Key Laboratory of Multi-phase Complex Systems, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China 2. School of Chemical Engineering, University of Chinese Academy of Sciences, Beijing 100049, China 3. School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

Received:2019-01-23 Revised:2019-04-15 Online:2019-12-22 Published:2019-12-22
Contact: Wei GE wge@ipe.ac.cn
Supported by:
;Key Research Program of Frontier Sciences of Chinese Academy of Sciences

摘要/Abstract

摘要： 发展了适应高超声速稀薄流动的硬球-拟颗粒模拟(HS-PPM)，并分别采用硬球分子动力学模拟(HS)、HS-PPM和直接模拟蒙特卡洛(DSMC)对马赫数为5、努森数为0.8的圆球绕流进行模拟，证明了HS-PPM可以得到更接近HS的模拟结果。对雷诺数100、马赫数5~19的三维圆球绕流，采用HS-PPM模拟得到了填充率0.01~0.08、完全热边界和完全滑移边界条件下的曳力系数，与HS模拟结果的一致性较好；模拟了马赫数24、努森数0.11~4.55的零攻角三维尖锥绕流，结果与文献中DSMC的模拟结果相符。研究验证了HS-PPM处理高超声速稀薄气体流动的可行性。

关键词: 硬球-拟颗粒模拟(HS-PPM), 高超声速流动, 稀薄气体

Abstract: Hard-sphere/pseudo-particle modeling (HS-PPM) has been demonstrated as an effective discrete simulation method for supersonic gas flow, but the simulated flow was limited to relatively low Mach number (Ma?3~5) yet. Recently, hypersonic rarefied gas flow has attracted great attention in aerospace and material engineering under extreme conditions. In this work, flow past some objects with simple geometry was simulated in HS-PPM. The geometric method was used to describe the wall and specular, and diffuse reflection was coupled to represent the slip thermal boundary condition. The tangential accommodation coefficient was introduced to adjust the proportion of slip and thermal accommodation. The simulated domain size was analyzed to determine the optimal values and the influence of tangential accommodation coefficient on the drag coefficient was considered. The hard sphere (HS) modeling, HS-PPM and direct simulation Monte Carlo (DSMC) method were used to simulate the flow past a sphere with a Mach number of 5 and a Knudsen number (Kn) of 0.8, which proved that the HS-PPM results were closer to those of the HS model. On this basis, 14 cases of the flow past a three-dimensional sphere with a Reynolds number (Re) of 100 were simulated, and Mach numbers change from 5 to 19. The upper and lower limits of the simulated drag coefficient were obtained by using the fully diffuse boundaries and the slip boundaries respectively, and were in good agreement with the corresponding results of the HS model. In addition, the simulated drag coefficient of cones with axial flow at a Mach number of 24 and Knudsen numbers from 0.11 to 4.55 was also obtained. The attack angle was zero and the results were consistent with the results of DSMC. This study demonstrated that HS-PPM was effective for hypersonic rarefied gas flow. Furthermore, the error caused by the invariable collisional cross of hard spheres in simulating hypersonic flow was found, indicating the direction of future improvement.

Key words: Hard-sphere/pseudo-particle modelling (HS-PPM), Hypersonic flow, Rarefied gas

赵祺赵明璨马琳博葛蔚. 硬球-拟颗粒模拟高超声速稀薄气体流动[J]. 过程工程学报, 2019, 19(6): 1093-1100.

Qi ZHAO Mingcan ZHAO Linbo MA Wei GE. Hard-sphere/pseudo-particle modeling (HS-PPM) for hypersonic rarefied gas flow[J]. Chin. J. Process Eng., 2019, 19(6): 1093-1100.

参考文献

[1] Bird G A. The DSMC method [M]. CreateSpace Independent Publishing Platform, 2013:13.
[2] Bird G A. Approach to translational equilibrium in a rigid sphere gas [J]. The Physics of Fluids, 1963, 6(10): 1518-1519.
[3] Antonov S, Ivanov M, Kashkovskii A, et al. Influence of atmospheric rarefaction on aerodynamic characteristics of flying vehicles [J]. Rarefied Gas Dynamics, 1990: 522-530.
[4] Et Defense A E, Les Mureaux F. Computation of Transitional Flows Around Three-Dimensional Re-Entry Bodies [J]. Rarefied Gas Dynamics: Space Science and Engineering, 1992: 161.
[5] Moss J N, Price J M. DSMC simulation of OREX entry [C]//Proc. Int. Symp. Rarefied Gas Dyn., 20th, Beijing, 1996.
[6] Ivanov M, Markelov G, Gimelshein S, et al. High-altitude capsule aerodynamics with real gas effects [J]. Journal of Spacecraft and Rockets, 1998, 35(1): 16-22.
[7] Ge W, Li J. Macro-scale phenomena reproduced in microscopic systems—pseudo-particle modeling of fluidization [J]. Chemical Engineering Science, 2003, 58(8): 1565-1585.
[8] 葛蔚, 麻景森, 张家元, 等. 复杂流动多尺度模拟中的粒子方法 [J]. 科学通报, 2005, 50(9): 841-853.
[9] 沈国飞. 纳微流动的分子动力学模拟 [D]. 中国科学院研究生院(过程工程研究所), 2010: 73-98.
[10] Zhang C, Shen G, Li C, et al. Hard-sphere/pseudo-particle modelling (HS-PPM) for efficient and scalable molecular simulation of dilute gaseous flow and transport [J]. Molecular Simulation, 2016, 42(14): 1171-1182.
[11] 张成龙. 硬球—拟颗粒—软球耦合模拟及其化工应用 [D]. 中国科学院大学, 2016: 23-65.
[12] Chapman S, Cowling T G. The mathematical theory of non-uniform gases [M]. Cambridge University Press, 1970.
[13] Chen F, Wei G. Statistical properties of pseudo-particle systems [J]. Particuology, 2010, 8(4): 332-342.
[14] Bird G A, Broadwell J E, Vogenitz F W. Leading edge flow by the Monte Carlo direct simulation technique [J]. AIAA Journal, 1971, 8(3): 504-510.
[15] Santos W. Aerothermodynamic characteristics of flat-nose power-law bodies in low-density hypersonic flow [C]//22nd Applied Aerodynamics Conference and Exhibit. 2004: 5381.
[16] Santos W, Lewis M J. Calculation of shock wave structure over power law bodies in hypersonic flow [J]. Journal of Spacecraft and Rockets, 2005, 42(2): 213-222.
[17] Santos W. Some Extentions to the Aerothermodynamic Performance Study of Flat-Nose Power-Law Bodies [J]. Nature Immunology, 2013, 2(2): 172-180.
[18] Bailey A, Hiatt J. Free-flight measurements of sphere drag at subsonic, transonic, supersonic, and hypersonic speeds for continuum, transition, and near-free-molecular flow conditions [R]. Arnold Engineering Development Center Arnold AFS TN, 1971.
[19] Takata G Y, Vogenitz F W. Rarefied hypersonic flow about cones and flat plates by Monte Carlo simulation [J]. AIAA Journal, 2015, 9(1): 94-100.
[20] Geiger R. Semi-Annual Report [R]. DAH C04-67-C, 1968.
[21] Horstman C C, Kussoy M I. Cone drag in rarefied hypersonic flow [J]. AIAA Journal, 2015, 8(2): 315-320.

硬球-拟颗粒模拟高超声速稀薄气体流动

Hard-sphere/pseudo-particle modeling (HS-PPM) for hypersonic rarefied gas flow

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 1

编辑推荐

Metrics