基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.      

先进飞行器气动隐身一体化设计方法研究现状与趋势

介绍了先进飞行器气动隐身一体化设计的研究现状与发展趋势。首先介绍了设计中的常用的数值模拟技术。针对气动力计算,重点分析对比了位势方程、Euler方程和N—S方程等气动力计算方法;针对隐身性能计算,重点分析对比了高频近似方法、时域法、频域法等雷达散射截面计算方法。在此基础上。讨论了先进飞行器气动隐身一体化设计的策略和步骤,最后对该设计方法的发展趋势进行了分析。     

激波诱导燃烧模拟中轴对称边界数值异常研究

对钝体激波诱导燃烧模拟中引起对称轴附近出现数值异常现象的几个影响因素进行了细致讨论.在飞行速度低于爆轰速度情形的模拟中,在一定条件下靠近对称轴的位置会出现明显的数值异常——燃烧阵面扭曲.分别从边界处理方法、计算格式、化学反应机制和计算网格等4个方面对该异常产生的原因进行了详细地讨论.研究表明4个因素都会导致数值异常产生,深入的分析表明计算格式是其中最根本的因素,而这种燃烧阵面扭曲的数值异常应该与粉刺问题(Carbuncle)类似.针对如何克服轴对称边界数值异常提出了一些合理的措施,并给出了验证算例的对比与分析.      

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

带舵旋成体侧向喷流流场特性分析

针对带舵旋成体外形,来流马赫数为2.0,迎角为0°,10°和20°,喷口位于×字舵中间,首先开展了喷主流干扰流场的数值模拟与风洞实验校验,得到相符的喷口附近子午线压力分布,验证了数值模拟方法的可行性;然后根据数值模拟得到的喷主流干扰流场结构和流动参数分布,分析了喷主流干扰流场结构随迎角的变化特点;最后从喷口附近旋成体和尾舵表面压力分布的变化分析了喷主流干扰对气动力/力矩放大因子的影响.研究结果表明:随迎角增大,喷流前激波向喷口靠近,分离区减小,喷主流压比增大,喷流穿透高度增大;气动力/力矩放大因子均大于1.0,随迎角增大,放大因子减小.喷口位于×字舵中间时,喷主流干扰对喷流直接力有增益作用;×字舵的存在减弱了喷流的包裹作用,对流场起到了有利干扰.      

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

不同进口马赫数下超声速燃烧流场特性研究

对有无楔板超燃冲压发动机模型内横向氢气喷流超声速燃烧流场进行了数值模拟,分析了进口马赫数对超声速燃烧流场特性及特征参数分布的影响特性。采用有限体积法求解多组元Navier-Stokes(N-S)方程,对不同进口马赫数下的燃烧流场进行了数值模拟,细致对比了流场激波结构、喷流穿透深度、燃烧阵面,燃烧效率及总压恢复系数等参数随进口条件的变化特征。结果表明:无论是否存在楔板结构,喷口后流场压强均随着进口马赫数的增加而减小,并且随进口马赫数的增加,氢气喷流穿透深度减小,楔板对喷流穿透深度基本无影响。较无楔板结构而言,设置楔板结构可以缓解燃烧室内流场对马赫数变化的敏感度,使燃烧更为稳定。在同一进口马赫数条件下,楔板布局有明显的促燃作用及激波点火效果,在一定程度上可增加此类发动机工作的马赫数范围,但以总压恢复系数略微降低为代价。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。