钝后缘翼型的后缘减阻装置研究

采用数值模拟方法进行了FB-4000-1000与FB-4000-2000钝后缘翼型的后缘减阻装置研究.分析了后缘减阻隔板安装位置对翼型气动特性的影响.其中后缘台阶中部安装隔板,以升力的小幅损失为代价实现中小迎角范围内减阻约35％,而在较大迎角时减阻约20％.     

Detached-eddy simulation of supersonic flow past a spike-tipped blunt nose

Space vehicle in atmosphere travels mostly at supersonic speed and generates a very strong bow shockwave around its blunt nose. Oblique shock and conical separated flow zone generated by a forward disk-tip spike significantly reduce the drag by reducing the high pressure area on the blunt nose. This study employs improved delayed detached eddy simulation to investigate the characteristic flow structures around a spike-tipped blunt nose at Mach number of 3 and Reynolds number(based on the blunt-body diameter) of 2.72x10~8. The calculated time-averaged quantities agree well with experimental data. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. It is found that two distinct instability modes exist: oscillation mode and pulsation mode. The former is related to the foreshock/turbulence interaction with nondimensional frequency at around 0.004. The latter corresponds to the interaction between turbulence and shock structures around the blunt nose, with a typical coherent structure shedding frequency at 0.092.     

计算网格对气动力气动热数值模拟影响的研究

以二维圆柱钝头超声速流场为研究对象,应用计算流体力学方法研究了计算网格在气动力及气动热数值模拟中的影响.计算程序中对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果表明,壁面法向第一层网格间距、网格长宽比以及扩张度等因素都会对壁面摩阻、热流的计算产生影响,尤其是壁面法向第一层网格间距影响最为突出.边界层的计算结果与实验数据进行了仔细对比,包括壁面压力分布,摩阻系数Cf分布,斯坦顿数St分布等.高马赫数滞止气流导致壁面附近产生巨大速度梯度以及温度梯度,这些梯度的准确计算需要合适的壁面法向第一层网格间距以及合理的网格长宽比、扩张度.     

高超声速边界层转捩对旋转钝锥自由飞运动的影响

通过在钝锥模型表面上布置人工绊线促使边界层强迫转捩,采用运动自由度不受约束的风洞模型自由飞试验技术研究边界层转捩对高超声速旋转钝锥自由飞行运动特性和气动特性的影响规律,并与自然转捩的旋转钝锥风洞模型自由飞试验结果作对比分析,试验马赫数为5.0,以模型长为特征尺寸的自由流雷诺数为1.68×106。研究结果表明:有人工绊线的旋转钝锥在自由飞行过程中有"激励稳定"的绕流流场,产生动态稳定的自由飞运动(动稳定导数系数小于0),而无转捩绊线的旋转钝锥在自由飞行中则有"激励不稳定"的绕流流场,产生动态不稳定的自由飞运动(动稳定导数系数大于0)。     

相变热图技术在表面热流测量中的应用

本文详细地介绍了在中国空气动力研究与发展中心低密度风洞中，进行的相变热图技术的应用工作。实验模型是尖锥及球锥，实验条件是自由来流马赫数Ｍ∞＝１６，来流雷诺数Ｒｅ∞／ｍ＝０．５×１０５～２．０×１０５，来流滞止温度Ｔ０＝９４０Ｋ，滞止压力ｐ０＝０．３９ＭＰａ～１．１７ＭＰａ。实验结果与国外计算结果作了比较。最大相对偏差２６．０％。本文对实验结果作了误差分析，介绍了相变漆响应滞后分析的一些结果。同时本文还介绍了相变热固技术的原理及方法。     

大钝头短体跨音速大迎角不同流型的计算比较

 采用不同流型 (无粘,层流,湍流 )的数值模拟方法计算了绕大钝头短体外形的跨音速大迎角时的流态,并进行了比较分析。数值方法中选用 van Leer分裂格式离散无粘通量项,并构造了一种 Limiter函数以保证 TVD性质;采用中心差分格式来离散粘性通量项。湍流计算中采用了 B-L湍流模型。算例表明不同流型的流态共同特点是三维分离的非定常性,但不同流型计算所得流态的差异较大。     

V字形钝前缘激波反射迟滞现象

针对内转式进气道唇口在宽速域条件下面临的复杂激波干扰问题,将唇口模化为V字形钝前缘,采用数值模拟并辅以风洞实验,研究了典型V字形构型（根部倒圆半径R与前缘钝化半径r之比R/r=1,半扩张角β=18°）激波反射结构随来流马赫数Ma∞的演变过程。结果表明,随着Ma∞的增大或减小,V字形后掠前缘的脱体激波产生规则反射（Regular reflection,RR）和马赫反射（Mach reflection,MR）,并且两者的相互转变过程出现迟滞。初场为RR时,V字形根部产生大范围的流动分离和分离激波;随着Ma∞由5.7逐渐增大至6.5,脱体激波的交点向下游移动并与分离激波的交点重合,使RR转变为MR。初场为MR时,马赫杆下游存在大尺度的反转涡对;随着Ma∞由6.7逐渐减小至5.9,反转涡对不再影响脱体激波,使MR转变为RR。通过Ma∞=6的风洞实验证实,在相同来流条件下存在RR和MR双解。基于对脱体激波交点、分离激波交点和反转涡对尺度随Ma∞变化规律的认识,建立了RR?MR的转变边界。在双解区中,RR工况的壁面压力最大值约为MR工况的2～3倍,表明迟滞现象将导致唇口气动载荷突变。     

高超声速二维钝楔边界层参数分析

为了指导二元高超声速进气道型面设计,基于两类高超声速二维钝楔外形,采用高精度CFD方法计算了四种典型工况下二维钝楔边界层参数,分析了前缘半径、壁面温度、压力梯度对边界层参数的影响,得出了边界层位移厚度、边界层总压恢复及壁面压力梯度沿流向的变化规律。研究表明,前缘半径及壁面温度对边界层参数影响较大,在两者一定的情况下,逆压梯度使边界层位移厚度增速放缓,顺压梯度会使边界层位移厚度快速增长。     

带配平翼钝体高超声速粘性绕流的数值模拟

采用ＮＮＤ差分格式，通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带配平翼钝体的高超声速粘性绕流。文中首先以球头驻点压强和热流为参数，讨论了在差分方程左端采用ＬＵ－ＳＧＳ、ＬＵ－ＡＤＩ和对角化ＡＤＩ三种不同隐式算法收敛效率的异同。然后以带配平翼钝体的高超声速粘性绕流为模型，对所研制的程序进行了计算验证。在计算中采用了代数方法和求解椭圆型方程方法相结合的网格生成技术，针对配平翼外形给出了贴体性、     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。