超声速光学头罩流场的PIV研究

在马赫数Ma=3. 8超声速风洞中.采用PIV(Particle Image Velocimetry,粒子图像测速)技术测量了超声速光学头罩流场的速度分布.PIV技术应用于超声速流场时,对系统的硬件配备、示踪粒子的跟随性以及PIV算法的精度有很高的要求.本文PIV系统选用高精度的同步控制器和高能量激光器;以纳米级粒径的粒子作为示踪粒子,通过斜激波响应实验分析了其在超声速流场中的跟随性;并采用多种高精度速度场算法对粒子图像进行处理.实验结果表明,示踪粒子在超声速流场中有很好的跟随性,采用的高精度速度场算法能够很好地反映超声速光学头罩流场的速度分布.     

DSMC算法气体壁面相互作用模型

研究了二维非结构网格高超声速稀薄流DSMC方法中气体与壁面相互作用的数值模型。发展了一种基于辐射平衡的壁面温度边界条件,与恒温壁面边界条件相比,该温度边界条件可以克服恒温壁面边界条件的自身缺陷并且能够在流场模拟中适时给出更贴近真实情况的壁面温度;采用了由完全漫反射和纯镜面反射模型组合而成的Maxwell壁面反射模型,该壁面反射模型能够更好地描述气体模拟分子在壁面散射的真实情况。采用本文方法对钝头体外形进行了数值模拟,结果表明本文采用的气体壁面相互作用模型能够提高钝头体背风面后部流场对温度变化的敏感度,并且随着飞行高度提升,敏感度有减弱的趋势。     

瑞利散射测速技术在高超声速流场中应用研究

采用基于法布里-珀罗干涉仪的干涉瑞利散射测速技术在Φ0.3m高超声速低密度风洞中进行了Ma5、Ma6、Ma12的流场速度和湍流度的测量,了解了瑞利散射速度和湍流度测量系统在高超声速流场中应用的情况,结果表明目前该风洞流场湍流度在1%以内,速度测量结果与流场校测偏差最大1.3%;对激波后返回舱模型绕流速度进行了测量,Ma6来流的测量结果与数值模拟结果吻合较好,而Ma12来流的测量结果与数值模拟结果相差69%,对原因进行了分析。在实验中发现目前Φ0.3m高超声速低密度风洞的流场存在一定程度的冷凝现象,并对后续研究工作提出了建议。     

犖犘犔犛技术及其在高速飞行器气动研究中的应用

近年来,与高速飞行器相关的超声速/高超声速流动受到了极大关注。这类流动所具有的非定常性、强梯度和可压缩性对试验研究提出了挑战。纳米示踪的平面激光散射技术(NPLS)是2005年由作者所在的研究团队研发的非接触光学测试技术。它能够获得超声速三维流场的某个剖面的瞬态流动结构,并且具有较高的时空分辨率。目前,许多研究结果表明NPLS是研究超声速湍流的一项非常有效的技术。近年来,作者应用 NPLS 技术在超声速湍流研究中取得了较大的进展,并且基于NPLS开发了其它几种技术,比如基于 NPLS 的密度场测量技术(NPLS-DT),能够获得超声速流动的密度场信息并还能进一步得到雷诺应力分布。本文介绍了NPLS技术并回顾了其在超声速边界层、激波/边界层相互作用等流动中的应用。由于能够获得雷诺压力和湍动能等统计量, NPLS技术有望在发展可压缩湍流模型的研究中发挥作用。     

气固两相自由射流的瞬态流场研究

为深入了解气固两相流的瞬态特性，对出口内径为10mm的气固两相自由射流流场进行了PIV实验研究，同时应用Fluent软件对其进行了数值计算加以对比。对实验和数值计算的气固两相速度场和固相浓度场的分析发现，气固两相射流中，固相的分散度小于气相，两相间存在明显的滑移速度；固相粒子集中在射流轴线附近，在射流中的扩散小于气相。     

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

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

电子束荧光技术用于测量高超声速飞行器流场密度

一种新发展起来的用电子束荧光技术直接测量高超声速模型周围流场密度的方法,首次在科学院力学所 JF-4B 高超声速炮风洞中对7°尖锥在 M_∞=7.8时进行了流场密度的测量,给出了模型表面边界层中分离与不分离两种情况下的平均密度和脉动密度剖面分布及电子束荧光的校准结果。验证了这项新的测量技术能成功地用于高超声速飞行器复杂流场中的密度测量。     

带气膜冷却结构的高超声速平板不同前缘形状下表面传热特性研究

当带红外成像制导系统的飞行器在稠密大气层内做高超声速飞行时,必须采取主动冷却方式防止严重气动加热造成的窗口材料热畸变以及复杂流场造成的气动光学畸变。本文根据成像窗口周围流动具有受高超声速钝头体绕流和气膜冷却结构(即背面为空腔的超声速后台阶)共同作用的特点,在 KD-01高超声速炮风洞中开展了带气膜冷却结构的高超声速平板在不同前缘形状下表面传热特性的试验研究,测量了 Ma8来流条件下喷缝下游表面传热系数,试验获得了2种前缘形状的带气膜冷却结构的高超声速平板喷缝周围瞬态流场 NPLS 图像。通过分析试验数据,得出以下结论:对于带气膜冷却结构(气膜不工作状态)的高超声速平板,模型前缘的形状对喷缝下游区域的表面热流整体分布有明显影响,在钝前缘情形下,表面热流分布接近相同前缘形状的平板边界层为层流状态时的表面热流分布;在尖前缘情形下,表面热流分布则表现出从层流边界层状态向充分发展湍流边界层状态变化的特性;喷缝下游分离和再附区表面传热特性和超声速后台阶流动类似,取决于喷缝上缘处边界层相对厚度。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

比热比和压比对高超飞行器尾喷流影响的实验研究

采用比热比为1.25的四氟化碳和空气的混合气体,模拟了超燃冲压发动机出口高温燃气的比热比。采用模型内喷管模拟发动机内喷流,风洞流场模拟飞行器外流。在0.5m常规高超声速风洞中,建立了模拟吸气式高超飞行器热态尾喷流干扰研究的实验手段,开展了喷流比热比对吸气式高超声速飞行器后体区域气动性能影响的实验研究。比较了相同外流和喷流落压比条件下,纯空气和混合气体喷流在喷流干扰区域的压力分布及流场结构。结果显示,混合气体喷流和空气喷流在喷流干扰区域的流场及表面压力分布差别明显。实验证实了喷流比热比是一个不可忽视的重要因素,在研究吸气式高超声速飞行器喷流干扰问题时应准确模拟。     

隔离段激波串精细结构与压力特性实验研究

在马赫数为2.5的等截面隔离段风洞中开展了无控制和安装T形涡流发生器两种情况的瞬态流场结构显示与压力测量的实验研究。运用常规纹影和基于纳米示踪的平面激光散射技术（NPLS）对两种不同状态的隔离段激波串三维流场精细结构进行了显示测量。结果表明:较传统纹影的测量结构而言,NPLS精细测量能够得到湍流边界层、激波串、分离区等细节结构。T形涡流发生器产生的展向涡与激波串相互作用,激波串前缘结构为分叉正激波,紧跟其后的第二道激波实际上结构与其类似。同时采用高频压力传感器对两种隔离段中激波串的壁面压力进行了测量,采用常规统计分析方法和差分平方累和方法对激波串压力分布、脉动及其上传特性进行了分析。分析表明,差分平方累和方法可以有效检测激波串的前缘位置。     

GHOST FLUID METHOD FOR RICHTMYER-MESHKOV INSTABILITY SIMULATION

INTRODUCTIONThe energy of the Richtmyer- Meshkov in-stability is excited by the interaction between ashock wave and a contact interface across thatfluids are different.Itis very important thatcap-ture the interface in multi- material flows exact-ly.Here we will apply the ghost fluid method tosimulate two dimensional Richtmyer- Meshkovinstability.The mathematical equation used to simulatethe Ritchmyer- Meshkov instability is the com-pressible Euler equation for two ideal gases.ρρuρvE…     

超声速弹头凹型光学头罩流动显示研究

在自行设计的Mα=3.8超声速风洞中,采用基于纳米技术的平面激光散射(NPLS)方法对超声速弹头凹型光学头罩流场进行了流动显示实验.高时空分辨率的NPLS图像再现了激波、膨胀波、边界层及尾迹等流场结构.观察到了边界层的产生、发展及转捩过程.通过对时间相关图像的分析,可以精确测定边界层内大尺度结构的几何特征和时间演化特征.     

二阶矩湍流模型在喷管两相流中的应用

本文应用一种考虑气体一颗粒湍流相互作用的二阶矩模型对喷管跨音速两相流动进行了数值模拟,研究了改变颗粒直径和载荷比对流动的影响。计算结果和实验数据及前人研究结果相吻合．表明二阶矩模型能够对跨音速两相湍流做出很好地预测。     

偏心底吹气体搅拌下冶金熔池水模型内三维湍流的实验研究

用一个圆柱形水模型模拟偏心底吹气体搅拌的冶金熔池,对模型内的三个典型剖面测量了三维速度场和湍流参数分布。方法是用二维激光多普勒测速仪在两个方向上测两次合成。面内的平行速度分布表明小流量偏心底吹熔池内的流谱为正循环流动。垂直速度分布表明,在气水两相流气泡柱附近有两个较强的圆周方向的环流,这种现象在中心底吹轴对称二维熔池流场中是不存在的。在气泡柱附近的区域内,速度模与湍动能值与较远的区域内的值相比要大得多,这就增强了这里的传热、传质和化学反应过程。在这里,气泡柱与壁面间的距离对湍流度的大小有影响,在速度值相近的情况下,较宽的区域内具有较大的湍流度值。     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

含化学反应流动的TVD格式及其在燃气射流中的应用

以Harten标准TVD格式为基础，结合固体火箭燃气射流的特点，以数学方法系统地推导了适合于高温、高压和高速流体流动的数值格式，给出了压力偏导数的合理计算公式；利用化学动力学知识，对火箭燃气射流流场中存在的多组分、含有限速率的化学反应问题进行了论述，阐明了化学反应质量源项的求解方法。以12组分9反应方程模型为例，利用编制的计算程序，对某火箭燃气自由射流流场进行了模拟。通过对结果的分析，肯定了数值格式的正确性。