首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 62 毫秒
1.
李锦  耿湘人  陈坚强  江定武  李红喆 《航空学报》2019,40(4):122176-122176
为评估直接模拟蒙特卡罗(DSMC)方法中新近发展的一种量子动理学(QK)化学反应模型,在自研高超声速稀薄气体流动数值模拟软件RariHV中实现该模型,并进行了单个绝热单元内的热泳平衡/非平衡反应测试。结果表明,数值解与理论解吻合良好,QK模型具有很好的预测化学反应流动的能力。为进一步验证QK模型在真实流动中的表现,还计算了高超声速绕圆柱的化学反应流动,QK模型计算的驻点线上温度和组分分布都与总碰撞能(TCE)模型结果吻合良好。相较于TCE模型,QK模型的优势是不再依赖反应速率系数的实验数据,在深空探测等化学反应数据缺乏的领域具有一定的应用前景。  相似文献   

2.
基于动理学模型的多尺度随机粒子方法   总被引:1,自引:0,他引:1  
传统随机粒子方法的时空离散步长受分子碰撞尺度(分子平均自由程和平均碰撞时间)的限制,当空间和时间离散尺度远大于碰撞特征尺度时,其输运系数的数值误差显著增加,因此如直接利用这类方法对跨流域流动精确求解,其计算效率往往是极低的(如DSMC方法,Fokker-Planck和BGK模型随机粒子方法)。通过对随机粒子方法输运系数离散误差的分析可知,这主要是因为传统算法将模拟分子的运动和碰撞解耦计算引起的。针对这一问题,本文介绍了适合于跨流域流动模拟的多尺度Fokker-Planck和BGK模型随机粒子方法。通过分子运动求解中耦合碰撞作用,在连续流区域,改进的随机粒子方法在较大的时间步长下仍能够满足宏观流体力学方程的输运性质。理论和计算结果显示,多尺度Fokker-Planck和BGK模型随机粒子方法可以高效准确地模拟从稀薄流到连续流的跨流域气体流动。  相似文献   

3.
基于Boltzmann-Shakhov模型方程,建立其沿特征线离散的一般形式,离散过程中对于碰撞项的处理采用显式和隐式加权平均的方法,其中权系数依赖于当地努森数,可根据当地流动尺度不同进行自适应调节。通过权系数的引入,对文献现有离散统一气体动理学格式进行改进,发展出具有尺度自适应特性的离散统一气体动理学格式(Scale Adaptive Discrete Unified Gas Kinetic Scheme,SADUGKS)。将SADUGKS格式应用于若干典型可压缩流动,对格式的有效性和尺度自适应特性进行了检验,所得数值结果与文献已有结果吻合较好,表明SADUGKS格式是一种求解宽范围努森数变化、跨流域多尺度流动问题的有效算法。  相似文献   

4.
首次提出了再入飞行过程中"小Knudsen数特征区(即[Kn1,Kn2]域)"的概念,并用直接模拟蒙特卡罗(DSMC)算法与求解Navier-Stokes方程两种模型详细研究了在这个区域中的流动.数值计算表明:采用王保国教授课题组自己编制的DSMC源程序和求解多组分、考虑热力学非平衡以及化学反应非平衡、弱守恒型Navi...  相似文献   

5.
真实流动环境下的真空羽流必然存在着多种不确定性,那么确定条件输入的数值模拟会存在偏差,因此需进一步研究不确定性对羽流流动特征的影响规律。本文采用直接模拟蒙特卡洛(DSMC)方法对不确定性输入的羽流流场进行模拟,采用稀疏的概率配置点方法对来流、壁面及模型参数等输入不确定量进行描述,对不确定性的传播和输出目标的平均值、方差及不确定度进行计算。研究表明,流场不确定性沿流线传播至流速最大处之后迅速增强,并在声速线前出现骤减的现象,传播至声速线之后,挡板壁面输入不确定性的影响凸显。最突出的是,压力不确定度在挡板驻点位置达到全场最大值,约为输入不确定量(3.54%)的2.1倍。此外,温度跳跃不确定度受到壁面温度不确定性输入的限制而近似保持为一个恒定值,约为输入不确定度的0.8倍。进而,壁面热流不确定度(5.54%)比壁面正应力不确定度(6.25%)略小,切应力不确定度最小(5.07%)。Sobol’全局敏度分析表明,喉道速度和喉道压力的输入不确定性对气动力/热不确定度的贡献是最大的,且远远超过了壁面温度和模拟分子直径不确定性输入的贡献。  相似文献   

6.
本文回顾了高阶气体动理学格式在湍流数值模拟中的应用.与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要.因此,基于两步四阶时间离散和高精度WE-NO重构,发展了具有四阶时间精度的气体动理学格式.该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力.目前...  相似文献   

7.
临近空间高超声速飞行器在高空长时间飞行,受气流加热影响,飞行器表面温度显著升高,依赖地面试验和传统DSMC仿真预测的热流值明显高于飞行观测值,导致飞行器防热系统的保守设计.本文发展了一种基于壁面辐射平衡的DSMC边界模型,通过热流值反算辐射平衡壁面温度,并以此温度作为下一个时间步DSMC计算的边界条件,迭代更新至给出壁面温度的收敛值.基于该温度边界条件,开发了适用于轴对称构型的DSMC求解器,并以钝锥构型对计算模型和求解器进行了验证.重点针对激波风洞试验条件下的双锥构型,开展数值模拟研究,结果表明:该构型恒温冷壁条件得到的壁面压力分布和热流与风洞试验结果吻合,两种温度条件下的压力峰值差异约为15.4%,但是整体气动力特性差异仅约为0.33%;相对于冷壁,辐射平衡计算得到的前缘处热流峰值降低约50%,再附点处的热流峰值降低约三分之二;两种条件相结合,可以给出壁面热流的预测范围.  相似文献   

8.
作为构造能够同时模拟稀薄和连续流动的混合算法的第一步,本文在前期工作的基础上,考察了气体动理学BGK-NS格式在复杂流动数值模拟中的能力,选取了高超声速再入飞行器大迎角侧向喷流及带栅格翼导弹超声速绕流两个算例.数值模拟结果与实验,以及文献中的计算结果都吻合很好.这充分说明了BGK-NS格式对复杂流动具有良好的模拟能力,也为下一步混合算法的建立打下了良好的基础.  相似文献   

9.
超/跨音风扇气动设计体系有关流动模型改进探讨   总被引:1,自引:0,他引:1  
当前超/跨音风扇准三维气动设计体系中,二维气动设计与三维流场分析在模型上不匹配。为加强二者之间的关联,提出了气动设计体系三角形概念。作为其实现方法的一部分,提出了内部气流脱轨角及其偏差作为定量指标,以对设计流场进行评价,并对下一轮设计进行指导。根据有限算例,探讨了内部以气流脱轨角及其偏差作为定量准则的可能性。   相似文献   

10.
神经网络在旋翼/机身气动干扰模型中的应用   总被引:1,自引:1,他引:1  
由于直升机自身的特点,旋翼/机身气动干扰呈现非线性,且受多种因素的影响,用神经网络来解决这一非线性问题是一个很好的办法。将旋翼/机身气动干扰试验数据构造的学习样本,对网络参数进行学习,可以得到旋翼/机身气动干扰神经网络模型,进而可以用该神经网络模型研究直升机机身受旋翼气流干扰时的空气动力特性。参数训练好的旋翼/机身气动干扰神经网络模型可直接用于直升机气动设计和实时仿真。笔者在对旋翼/机身气动干扰神  相似文献   

11.
高超声速稀薄流的气粒多相流动DSMC算法建模研究   总被引:1,自引:0,他引:1  
基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.  相似文献   

12.
高空羽流场的DSMC计算和实验研究   总被引:9,自引:2,他引:9  
以小推力发动机的高空羽流场为研究对象,实现了高空小喷管氮气流DSMC方法数值模拟,验证了DSMC方法在高空膨胀流中的可用性,并在现有设备基础上,进行一组高空羽流试验,给出试验与计算的对比结果。结果表明:喷管出口附近能得到满意结果,由于实现设备能力所限,远场试验没有达到预期效果。  相似文献   

13.
合并式电离—复合反应的DSMC模拟   总被引:1,自引:0,他引:1  
采用DSMC方法研究分子碰撞反应机理,在分析水平上模拟绝热封闭体系事并式电离-复合反应的驰豫过程。  相似文献   

14.
二维热化学非平衡流动的非结构网格DSMC方法及其应用   总被引:1,自引:0,他引:1  
研究了二维热化学非平衡流动非结构网格DSMC方法实现的过程。提出了一种新型的高效搜索算法,该算法不仅可以跟踪模拟分子在网格之间的迁移,而且在搜索过程中可以准确判别分子与物面是否相互作用,避免了原有算法中分子表面反射非确定论判据。为加快流场的时间发展历程,设计了适合DSMC方法的动态局部时间步长技术,使其可以应用到定常和非定常流场的计算。利用Fortran90的动态分配内存技术编制了计算程序。最后对过渡流域高超声速圆柱绕流进行了数值试验,计算结果初步验证了该算法的可行性。  相似文献   

15.
本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.  相似文献   

16.
稀薄流热化学非平衡效应的DSMC研究   总被引:2,自引:0,他引:2  
以钝体绕流驻点线流场为例,应用直接仿真蒙特卡罗法,模拟了五组元混合气体(N2,O2,N,O,NO)中的碰撞能松弛,离散反应,交换反应及复合反应,给出了碰撞分子各自由度间的能量交换模式及各类化学反应的选择模式。DSMC仿真结果与宏观粘性激波层方程计算结果,在飞行高度较低时基本相符,表明采用的热化学非平衡非模式是令人满意的,随着飞行高度增高稀薄气体效应增大,二者在激波附近有显著差异,表明高空飞行必须考  相似文献   

17.
本文发展了平衡粒子模拟方法(EPSM),建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型,并通过混合参数进行流区的自动识别,将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合,构造了可模拟化学反应流动的DSMC/EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟,将结果与DSMC方法的结果进行比较,验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较,发现混合算法能够大大提高计算效率。  相似文献   

18.
大瑞利数下二维Benard对流演化过程的DSMC仿真研究   总被引:2,自引:0,他引:2  
利用DSMC(Direct Simulation Monte Carlo)方法研究了大Rayleigh(瑞利)数条件下稀薄气体的二维Benard对流的特性。考察了流动的非定常演化过程及若干流动特征量随时间的变化特征,得到了与实验较为一致的结论。计算结果表明,当流动控制参数不断增大时,二维Berard对流从热传导状态过渡到热对流状态表现出明显的非定常演化特征。最后当瑞利数增大到一定程度后,注动显示出混沌的特性。  相似文献   

19.
采用DSMC(DirectSimulationMonte-Carlo)方法模拟球锥体过渡区流动。论文给出了不同Knudsen数球锥绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球锥阻力系数和底部阻力系数的分布,并将阻力系数模拟结果与桥函数一体化估算结果进行了比较分析。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号