解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究

推进技术 ›› 2014, Vol. 35 ›› Issue (7): 897-904.

解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究

唐振宇,贺碧蛟,蔡国飙

北京航空航天大学宇航学院，北京 100191;北京航空航天大学宇航学院，北京 100191;北京航空航天大学宇航学院，北京 100191

发布日期:2021-08-15
作者简介:唐振宇(1984—)，男，博士生，研究领域为稀薄气体动力学及发动机真空羽流仿真。

Investigation on Boundary Conditions in Decoupled N-S/DSMC Method for Vacuum Plume Simulation of Thrusters

School of Astronautics，Beijing University of Aeronautics and Astronautics，Beijing 100191，China;School of Astronautics，Beijing University of Aeronautics and Astronautics，Beijing 100191，China;School of Astronautics，Beijing University of Aeronautics and Astronautics，Beijing 100191，China

Published:2021-08-15

摘要/Abstract

摘要： 准确模拟羽流流场是合理评估空间推力器真空羽流效应的基础。因羽流流场同时包含连续和稀薄两种流动机理，通常采用解耦方式的Navier-Stokes(N-S)方程和直接模拟蒙特卡罗(DSMC)混合方法对其进行模拟。为保证解耦N-S/DSMC方法应用于羽流计算时的准确性并尽量提高其计算效率，对计算中的DSMC入口和喷管壁面等边界条件设置问题开展了研究。通过与文献中低总压工况的羽流试验数据比较，确定了合理设置DSMC入口边界位置、壁面反射类型、壁温等边界条件的方法。针对将此方法应用于实际推力器工况计算时效率过低的问题，通过多种数值试验，表明从喷管出口处开始沿KnGL为0.05的等值线作为DSMC入口界面可同时保证仿真精度和较高的计算效率；并证明实际工况下壁温设置对羽流场仿真结果影响不大。

关键词: 解耦N-S/DSMC方法；推力器真空羽流；边界条件；DSMC方法

Abstract: Accurate simulation of plume flowfield is the basis of plume effects assessment of space thrusters. Both continuum and rarefied flow regimes are contained in the plume flowfield of thrusters，therefore，decoupled Navier-Stokes(N-S)and Direct Simulation Monte Carlo(DSMC)method is usually employed to simulate this kind of flow problem. To ensure both high accuracy and high efficiency when using decoupled N-S/DSMC method to solve vacuum plume flowfield of thrusters，the DSMC inflow and nozzle wall boundary conditions are investigated. By detailed comparing with the experimental data for the plume flows of a low-total-pressure thruster in literature，reasonable methods for setting boundary conditions，including the location of DSMC inflow boundary，reflection type of nozzle wall and wall temperature，are set up. When using these methods to simulate practical high-total-pressure thrusters，computational efficiency reduces seriously. Series of numerical tests are conducted，which show that the DSMC inflow interface can be set by starting from the nozzle exit and laying along the KnGL contour line of 0.05，meanwhile in the practical cases，wall temperature setting makes little influence on the simulation results.

Key words: Decoupled N-S/DSMC method；Vacuum plume of thrusters；Boundary condition；DSMC method

唐振宇,贺碧蛟,蔡国飙. 解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究[J]. 推进技术, 2014, 35(7): 897-904.

TANG Zhen-yu,HE Bi-jiao,CAI Guo-biao. Investigation on Boundary Conditions in Decoupled N-S/DSMC Method for Vacuum Plume Simulation of Thrusters[J]. Journal of Propulsion Technology, 2014, 35(7): 897-904.

参考文献

[1] Bird G A.Molecular Gas Dynamics and The Direct Simulation of Gas Flow[M].Oxford:Clarendon, 1994.
[2] 蔡国飙, 王慧玉, 祖国君, 等. 低密度小喷管流场的DSMC直接模拟[J].推进技术, 1996, 17(5): 9-13.(CAI Guo-biao, WANG Hui-yu, ZU Guo-jun, et al.The DSMC Simulation of Low Density Flow Field in Small Nozzle[J].Journal of Propulsion Technology, 1996, 17(5): 9-13.)
[3] 蔡国飙, 王慧玉, 祖国君, 等. 真空小喷管羽流场的Monte Carlo直接模拟[J].推进技术, 1997, 18(2): 44-49 (CAI Guo-biao, WANG Hui-yu, ZU Guo-jun, et al.Monte Carlo Direct Simulation of Small Nozzle Plume in Vacuum[J].Journal of Propulsion Technology, 1997, 18(2): 44-49.)
[4] 王平阳, 程惠尔, 杨卫华. DSMC法模拟双组元姿探发动机喷管流动[J].推进技术, 2003, 24(2): 112-114.(WANG Ping-yang, CHENG Hui-er, YANG Wei-hua.DSMC Simulation of Nozzle Flow of Bipropellant Attitude Control Engine[J].Journal of Propulsion Technology, 2003, 24(2): 112-114.)
[5] 蔡国飙, 王慧玉, 庄逢甘.真空羽流场的N-S和DSMC耦合数值模拟[J]. 推进技术, 1998, 19(4): 57-61.(CAI Guo-biao, WANG Hui-yu, ZHUANG Feng-gan. Coupled Numerical Simulation with N-S and DSMC on Vacuum Plume [J].Journal of Propulsion Technology, 1998, 19(4): 57-61.)
[6] 王平阳, 鹿洪芳, 程惠尔.双组元姿控发动机反流特性的DSMC模拟[J].工程热物理学报, 2004, 25(4): 640-642.
[7] Bird G A.Breakdown of Translational and Rotational Equilibrium in Gaseous Expansions[J].AIAA Journal, 1970, 8(11): 1998-2003.
[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.
[9] Wang W-L, Boyd I D.Predicting Continuum Breakdown in Hypersonic Viscous Fows[J].Physics of Fluids, 2003, 15(1): 91-100.
[10] Glass C E, LeBeau G J.Numerical Study of a Continuum Sonic Jet Interacting with a Rarefied Flow [R].AIAA 97-2536.
[11] 严立, 王平阳, 欧阳华.月面环境发动机羽流冲击力效应模拟计算[J].上海交通大学学报, 2012, 46(8):1310-1314.
[12] 王平阳, 亓雪芹, 程惠尔.壁面条件对微喷管流动的影响[C].北京: 2003空气动力学前沿研究论文集:455-461.
[13] 沈青.稀薄气体动力学[M].北京: 国防工业出版社, 2003.
[14] Boyd I D, Penko P F, Meissner D L, et al.Experimental and Numerical Investigations of Low-Density Nozzle and Plume Flows of Nitrogen[J].AIAA Journal, 1992, 30(10):2452-2461.
[15] He B J, Zhang J H, Cai G B.Research on Vacuum Plume and Its Effects[J].Chinese Journal of Aeronautics, 2013, 26(1): 27-36.
[16] 侯凤龙, 蔡国飙.DSMC方法的边界采样与统计技术[J].北京航空航天大学学报, 2010, 36(2): 159-162.
[17] Schwartzentruber T E, Boyd I D.A Hybrid Particle-Continuum Method Applied to Shock Waves[J].Journal of Computational Physics, 2006, 215(2): 402-416.
[18] Wu J S, Lian Y Y, Cheng G, et al.Development and Verification of a Coupled DSMC-N-S Scheme Using Unstructured Mesh[J].Journal of Computational Physics, 2006, 219(2): 579-607.(编辑:梅瑛)收稿日期:2013-08-07；修订日期:2013-09-23。作者简介:唐振宇(1984—),男,博士生,研究领域为稀薄气体动力学及发动机真空羽流仿真。E-mail:rain_tang@sa.buaa.edu.cn