数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.     

超高速流动模拟及热化学反应模型对比研究

超高速流动是飞行器再入大气层时所面临的高速高温流动环境,膨胀管是少数几种能模拟超高速流动的地面设备之一。采用数值模拟方法对超高速试验进行辅助分析诊断,流动模拟时热化学反应模型的选择对流场特性影响较大,分别选择5组分、11组分热平衡及5组分热非平衡模型,对比研究3种不同热化学反应模型对双楔试验模型数值模拟结果的影响,以进一步评估超高速流动模拟时热化学反应模型的适用范围。结果表明,试验气流条件下5组分化学模型即可满足要求,加速气流条件则必须采取11组分化学模型,而对于流动中热非平衡效应显著时,热化学非平衡模型更为适用。     

喷流干扰气动热数值模拟的若干影响因素

工程上对喷流气动热进行近似计算时,常常对真实的全化学反应流近似处理,这时需要了解喷流气体和喷管外形的近似方法对计算结果的影响规律。采用数值试验对比分析的方法,探讨了近似计算的2个影响因素:其一是使用量热完全的空气喷流近似燃气喷流进行数值模拟时,喷管几何形状对热流分布的影响;其二是使用空气喷流近似燃气喷流进行数值模拟的方法,分析不同的喷流热力学参数匹配方法对喷流形态及热流分布的影响。通过对不同方法进行对比计算,揭示了各个方法对喷流干扰气动热计算的影响规律,可以用于指导工程分析。      

高空高速无人飞行器热控制系统设计

针对飞行时间短、速度和高度变化快、表面温度波动大的无人飞行器UAV(Unmanned Aerial Vehicles)热控制系统设计难题,提出了一种可解决实际工程问题的热分析计算方法.即把热天工况、冷天工况和标准天工况作为设计/试验工况;采用参考温度法、高超音速工程预测法或计算流体动力学CFD(Computational Fluid Dynamics)数值模拟法,确定了飞行器表面温度分布,并把其作为后续热分析数学模型的外边界条件;分析结构热容量对瞬态热载荷的影响,建立与之相应的边值问题方程,并采用有限差分法求解;根据高空高速飞行特点及瞬态热载荷值,确定仪器设备舱调温系统方案.     

高超音速钝头体绕流反应速率与壁面催化问题

反应速率与壁面催化是反应流中两个重要问题。本文以高超音速化学反应钝头体绕流为例对这两个问题做了初步的分析和讨论。反应速率问题实质上是反应速率常数实验数据的误差及其影响,本文根据实验数据的分散度大小以样本数据组和等效速率常数不确定性因子这两个概念建立反应速率常数分析数据库,通过对数据误差和计算结果的分析与比较总结出高超音速钝体绕流高温特性与化学反应特性计算结果依赖于反应速率数据库的某些规律。壁面催化问题是固壁对与之接触的气体化学反应的催化影响,通过引入一个新的壁面催化率参数来考虑壁面催化性质对高超音速钝体绕流计算结果的影响。由于引入和强化了上述概念,给化学反应速率与壁面催化问题的分析与研究提供了有力的手段,因而提高了数值模拟结果的可信度。     

球头锥-柱再入飞行器的动力环境预示

考察了超声速和高超声速流场中一类球头锥-柱再入飞行器表面稳定压力和脉动压力的分布特性,利用CFX(Computation Fluid Explorer)软件建立球头锥-柱的外流场域网格,对其进行定常、非定常流场分析,给出一套预测表面稳定压力和脉动压力环境的工程新方法.将该方法与理论和经验估算方法进行了对比分析,验证了该方法的可行性和适用性.同时利用该方法研究了马赫数、来流攻角等因素对再入飞行器表面脉动压力环境的影响.利用ANSYS大型软件建立了再入飞行器导弹头部的动力学分析模型,根据确定的再入飞行器表面脉动压力环境和随机振动分析理论,采用有限元谱分析方法计算得到该飞行器弹头部件的随机振动响应加速度的均方根值和功率谱密度曲线,计算结果满足工程应用的要求.      

火星再入飞行器风洞试验与真实飞行之间相关性的探讨

由于风洞试验条件限制,难以完全模拟火星再入飞行器真实飞行环境,因此需要建立火星再入飞行器风洞条件与真实飞行之间的关联关系。基于国外文献公开数据,采用数值方法和对比分析方法探讨了类"探路者号"外形的火星再入飞行器的风洞试验与真实飞行之间的外推方法。结果表明,在高焓空气风洞和常规空气风洞试验条件下,可以将模型驻点附近的无量纲压力和压力系数作为相关性参数,将风洞条件与飞行条件相关联起来,但是不能直接利用风洞试验的热流、无量纲热流和Stanton数作为关联参数;在高焓CO₂风洞试验条件下,可以利用模型驻点附近的无量纲压力、压力系数和Stanton数作为外推参数,但是不能直接将风洞试验的热流、无量纲热流作为相关性参数,将风洞条件下的风洞数据通过外推获取飞行条件下飞行器的性能参数。      

高超声速三维化学非平衡流动的数值模拟

高温非平衡化学反应对高超声速再入飞行器绕流流场的气动热特性有较大影响.通过采用有限体积NND(Non-oscillatory, containing No free parameters and Dissipative scheme)格式求解三维N-S方程,分别计算了小攻角条件下球锥的完全气体粘性流场和化学反应非平衡气体流场.比较并分析了化学非平衡作用对三维高超声速钝头体绕流流场的影响.所得结果与国外文献结果符合的较好,证明了计算方法与程序的有效性和实用性.      

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

栅格翼和机身组合体的气动特性计算与分析

给出计算栅格翼和机身组合体亚超音速气动力干扰理论方法.亚音速采用涡格法计算栅格翼气动力,用一个二维偶极子和圆柱映像涡系模拟机身的上洗流场和栅格翼迁移到机身上气动载荷.超音速假定栅格翼处在机身横向流动的上洗流场中,用翼片理论计算机身对栅格翼气动力干扰的修正因子.用本文方法计算了几个栅格翼和机身组合体的气动特性的到了比较满意的结果.      

轴对称DCR进气分流流场数值模拟

用矢通量分裂法，采用ＭａｃＣｏｒｍａｃｋ二步格式，对ＤＣＲ（ＤｕａｌＣｏｍｂｕｓ－ｔｏｒＲａｍｊｅｔ，简称ＤＣＲ）进气分流流场进行了数值模拟，提供了全场清晰的波系结构和物理信息，计算结果表明，通过改变反应，可以有效地控制结尾激波的位置，从而改变亚燃室内回流区的大小，对组织亚燃室燃烧和火焰稳定创造良好的条件。     

NUMERICAL MODELLING OF VELOCITY AND TEMPERATURE DISTRIBUTIONS OF THE BUOYANCY CONVECTION EFFECT IN KNbO3 MELT

Numerical modelling of velocity and temperature fields in high-temperature KNbO3 melt of a loop-shaped Pt wire heater is carried out by using the commercial com putational code ANSYS for the mathematical solution of the governing equations.Based on the experimental boundary conditions and the Boussinesq approximation,the numerical modelling of a steady and two-dimensional model is applied to study the process under consideration of the buoyancy-driven convection condition. The result is compared with the previous experimental and theoretical data obtained in our laboratory, and the former is in agreement with the latter. Thus a theoretical guide for reasonable growth conditions is provided by studying in depth the real fluid flow effects in the crystal growth from the melt.     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

二维非结构网格的可压缩Euler方程流场解

讨论了非结构网格的生成方法,用Delaunay三角化方法和推进阵面法相结合的方法生成了非结构网格,提高了生成非结构网格的自动化程度和网格的质量.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼形和两段翼型的算例,并与实验结果进行了比较和分析.     

栅格翼空气动力特性的计算与分析

将涡格法推广应用于计算外形复杂的蜂窝式栅格翼的气动特性，给出了亚、跨、超音速流中，栅格翼升力，力矩和阻力特性以及载荷分布，并分析了空气动力特点及几何参数的影响。计算结果表明，栅格翼作为一种新型的承力稳定面和控制面，在某些方面优于传统单翼，修正提出了空气动力等效条件，并通过计算进行了验证。     

CARET进气道亚音速气动特性研究

通过对一种Caret型进气道进行试验与流场数值计算研究,给出该型进气道亚音速的基本气动特性.在低速情况下,进/发匹配点处由于流量系数较大,管道内往往存在分离,这是总压恢复降低的主要原因.试验表明在进气道出口的左上角存在低能量区,而数值计算显示在进口外下角处首先出现流动分离,分离在管道内发展时并未沿周向旋转.