来流参数对防热瓦横缝旋涡结构及热环境的影响

针对高超声速飞行器表面缝隙内部流动,通过求解可压缩Navier-Stokes方程,自主研发了一套能够较好模拟缝隙流动特性的计算流体力学(CFD)软件。利用该软件研究了来流参数对防热瓦横缝旋涡结构及热环境的影响。计算结果表明:随着来流雷诺数的增加,缝内旋涡结构呈现主涡个数增多形态趋于饱满的变化趋势,缝隙壁面绝对热流和无量纲热流增加;随着来流马赫数的增加,缝内主涡个数、形态基本不变,但主涡旋转速度增加,缝隙壁面绝对热流增加,无量纲热流基本不变;随着来流迎角的增加(迎角较小时),缝内旋涡结构和热流变化规律基本与增加来流雷诺数相同。由此分析可知,涡量向下传递并形成旋涡的距离,即形成所谓"死水区"的深度,主要由来流雷诺数和来流迎角决定。     

高超声速钝双楔绕流流动转捩与分离流动的壁温影响

为研究壁温对吸气式高超声速飞行器进气道转捩流动的影响,选取钝双楔这一典型外形,基于德国Aachen工业大学Thomas与Herbert所开展的双楔高超声速风洞试验,分析了一些已有的计算流体力学(CFD)研究内容,并结合本文不同方法的CFD数值模拟结果,讨论了不同壁面温度对该双楔模型高超声速绕流流动转捩与流动分离的影响。对于双楔模型,流动分离一般发生在拐角附近,由于流动分离旋涡的剪切作用会诱发流动转捩,转捩又会改变流动分离强度、分离涡尺寸,若分离流动存在非定常特征则将导致非定常旋涡运动与流动转捩的复杂相互作用。通过比较已有文献的CFD数值模拟结果与本文计算结果,表明只有按照转捩思路开展的数值模拟才能够反映该风洞试验情况。计算结果与试验数据的比较显示,文献中按照第一压缩面层流与第二压缩面湍流状态计算得到的结果能够在一定程度上与风洞试验数据相符,本文使用MUSCL格式、剪切应力输运(SST)湍流模型与γ-Reθ关联转捩模型这种计算方法,得到的结果与试验数据符合较好,正确地反映了风洞试验情况。分析还表明,在分离流动之前的区域,随着壁面温度的升高,壁面热流会下降,近壁区域黏性系数变大,边界层内速度剖面不饱满,速度边界层较厚,厚的速度边界层容易发生流动分离现象。     

稀薄流区高超声速飞行器表面缝隙流动结构及气动热环境的分子模拟

针对高空稀薄流区的高超声速飞行器表面缝隙或缺陷结构导致的局部气动加热问题,采用直接模拟Monte Carlo(DSMC) 方法研究了70、75、80km和90km等4个飞行高度下稀薄流区高超声速缝隙流动问题,考虑稀薄气体效应和三维效应对缝隙内部流场结构和热流的影响。结果表明:上述飞行高度下,外部流动的分离和再附在缝隙内部形成一个充满腔体的单涡结构;稀薄气体效应对缝隙内部流动结构和壁面热流影响明显,随着高度的增加,主涡涡心上移,其形状逐渐变得“扁长”,右上角逐渐变尖,热流越来越集中分布于缝隙下游侧面的顶部区域;三维缝隙效应阻碍来流气体分子进入缝隙,导致主涡涡心上移,二维缝隙假设会高估缝隙表面的热流。      

临近空间高超声速飞行器表面非规则缝隙流动及传热特性研究

高超声速飞行器表面的缝隙会对其整体的气动传热和烧蚀特性产生重要影响,研究缝隙内部流动与传热机理对高超声速飞行器热防护系统的设计具有重要意义。针对临近空间高超声速飞行器表面缝隙的非规则变形问题,在热流梯度集中处构建倾斜式形变边界,并采用非结构网格DSMC方法对其内部的流动和传热特性进行模拟,分析非规则形变对缝隙内部涡旋流动及传热特性所造成的影响。结果表明:倾斜式非规则变形会同时产生高温引流和加速导流作用,而高温引流会大幅提高缝隙底部的气动热负荷,是缝隙底部热防护的难点问题。     

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

网格对高超声速钝头体表面热流数值模拟结果的影响

高超声速钝头体热流的精确预估，对高超声速飞行器热防护的设计具有重要意义，而网格划分尺度对预估飞行器表面热流精度具有重要的影响。应用理论分析及数值模拟方法，提出了一种基于壁面网格雷诺数及基于钝体特征长度的来流雷诺数网格划分方法，并给出了壁面网格雷诺数的取值范围。应用所提出的网格划分方法，对不同来流条件的高超声速半圆柱及球头表面热流进行了数值模拟，结果表明，应用所提出的网格划分方法，及给出的可供参考的壁面网格雷诺数的取值范围划分网格，在满足高超声速钝体表面热流精度要求的同时，有利于网格合理分布并提高数值计算效率。     

高超声速飞行器瞬态热试验

为了进行高超声速飞行器热防护系统的初步设计和数值计算的验证,设计开发了高超声速飞行器瞬态气动加热地面试验系统及其控制软件.试验系统能够根据飞行器的飞行轨迹和外形参数加载瞬态热流,实时测出结构表面的热流值和温度,得到飞行器的表面试验热流曲线和温度曲线.试验系统采用真空舱模拟飞行环境,并为此设计了冷却床,在真空环境下能比较真实的模拟热防护系统的下表面热环境,使瞬态热试验的原理更加合理,精度进一步提高.      

高超飞行器内流道激波振荡问题的数值研究及试验验证

采用数值模拟和风洞试验的方法,对吸气式高超声速飞行器盲腔状态的流动特性进行了研究。采用＂双时间步＂方法进行了内外流一体化的非定常数值模拟,利用彩色纹影系统对高超声速飞行器前体流场进行显示,并采用动态压力传感器测得了飞行器内流道的壁面压力随时间的变化。结果表明：盲腔构型在高超声速飞行中会出现周期性的激波振荡现象。数值模拟所得流场变化特征、内流道壁面压力振荡周期和壁面压力变化趋势与试验结果吻合良好。     

LENS风洞试验返回器模型气动热特性模拟

国家数值风洞高超声速流动模拟软件HyFLOW的研制对打破国外同类软件的技术壁垒具有重要意义。与国外DPLR软件进行了对比研究，同时系统介绍了HyFLOW软件求解器的数值方法、物理化学模型以及壁面催化特性计算模型等主要方法，采用典型算例对有限催化模型进行了数值验证，最后基于LENS风洞试验146 mm返回器标模外形开展了高超声速气动热特性数值模拟。研究结果表明，HyFLOW软件在高超声速热化学非平衡流动模拟与评估方面的气动力计算精度高，与国外同类软件DPLR相当，同时其壁面催化条件下的气动热计算精度可靠，可信度高。     

壁面催化对高超声速飞行器气动特性影响

针对高超声速流动中的高温真实气体效应,采用数值模拟求解三维N-S方程的方法研究了壁面催化对高超声速飞行器气动特性的影响规律。研究发现:对于文中所选两类高超声速飞行器———大钝头CEV再入飞行器和仿HTV2高升阻比升力体飞行器,壁面催化对表面压力影响均较小,对剪切应力的影响在飞行器不同部位表现不同:在头部和前缘等强压缩区域,壁面催化对表面剪切力影响明显,在大面积和背风区位置,壁面催化对表面剪切力影响微弱。这是因为壁面催化使得具有更大惯性的大分子气体在近壁处聚集,从而导致更高的速度梯度。由于大钝头外形波阻在整体气动特性中占优,而高升阻比外形头部强压缩区域面积较小,头部强压缩区域剪切力对整体气动特性贡献较小,最终体现为壁面催化对整体气动力影响微弱。     

A swept fin-induced flow field with different height mounting gaps

In order to apply the air fin successfully and ensure the maneuverability of hypersonic vehicle, a key problem to be studied urgently is the heat flux brought by the fin mounting gap. The appearance of mounting gap and fin shaft can induce many complex flow structures which need more attentions to be investigated. Under Ma 6, Nano-tracer-based Planar Laser Scattering (NPLS) and Temperature Sensitive Paints (TSP) were applied to visualize and measure transient flow structures and heat flux distribution of a swept fin-induced flow field with different height mounting gaps. Complementarily, Reynolds-averaged N-S equations were solved with k-ω SST turbulent model. The heat flux distribution results of numerical simulation and TSP observed the change of high heat flux region with different mounting gap, both in position and magnitude. The streamlines based on Computational Fluid Dynamics (CFD) and flow visualization results obtained by NPLS revealed the cause of high heat flux region. The high heat flux region in this flow field is mainly related to the reattachment of vortex and flow stagnation. The increase of gap height can lead to stronger gap overflow and shaft-induced horseshoe vortex, which are source of the high heat flux around the fin. The case with the highest mounting gap (4 mm) en-counters the most severe aerodynamic heating, both on the surface of fin and plate. Thus, under the premise of ensuring the flexibility of the fin, the gap should be set as small as possible.     

尖楔结构低速高/中温双路气流组合热试验方法

在典型高超声速飞行工况下,数值模拟分析了高温合金尖楔前缘结构沿气流方向大温度分布梯度将带来严重的热强度问题,产生大温度分布梯度的根本原因是尖楔结构头部区域平均热流密度与后段平板区域平均热流密度之差,而受头部区域热流密度具体分布的影响不大;进而提出了一种低速高/中温双路气流组合热试验方案,并通过数值模拟方法证明了该方案具有两股气流参数可以独立调节分别满足尖楔结构头部驻点区域及后段平板区域大、小两种热流密度的优点,进而解决单喷口低速高温燃气流热模拟试验难题,满足尖楔结构高超声速飞行工况下大温度梯度模拟要求.同时,该方案通过高/中温气源的合理组合搭配可以大大降低尖楔结构热试验所需高温气源发生功率,推广应用于电弧风洞可拓展其热试验范围.      

基于CFD的机翼结冰过程分析

采用一种适用于分析机翼结冰过程的数值模拟方法,该方法在飞机结冰气象条件下采用计算流体力学（Computational Fluid Dynamics,CFD）对机翼流场进行计算,获得流场中二维机翼翼型的水滴撞击特性,并遵循质量和能量守恒方程,计算出机翼表面的结冰量,预测出机翼的积冰形状。所得冰形与文献中的试验数据吻合良好,表明数值模拟方法有效。另外还分析了结冰对机翼升力、阻力和压力系数的影响,其结果可以指导飞机机翼防／除冰系统的设计研究     

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

高超声速前缘空腔数值模拟研究

介绍了一种前缘空腔高超声速被动流动控制方法并将其应用于轴对称钝锥外形,采用高分辨率的CFD数值模拟方法对不同参数的前缘空腔外形进行了对比研究,阐述了高超声速前缘空腔降低前缘壁面热流的机理,总结了不同空腔参数对空腔外形气动力、热的影响规律。结果表明,合理设计空腔参数可以使空腔外形的前缘壁面热流有较大程度的降低,并且不影响外形的升阻比等气动特性。     

悬停状态共轴刚性双旋翼非定常流动干扰机理

基于运动嵌套网格方法,建立了一套适合于悬停状态下共轴刚性双旋翼非定常干扰流场分析的计算流体力学(CFD)方法。首先,基于高效的运动嵌套网格技术,采用积分形式的可压雷诺平均Navier-Stokes(RANS)方程作为双旋翼非定常流场求解控制方程,湍流模型选用Baldwin-Lomax模型,时间推进采用双时间方法。在CFD方法的验证基础之上,对干扰过程中的桨尖涡涡核位置及强度演变规律进行了细致分析,揭示了共轴双旋翼非定常干扰流场中上、下旋翼桨尖涡与双旋翼桨叶之间的贴近干扰、碰撞现象,以及上、下旋翼桨尖涡之间的相互干扰机理。然后,进一步研究了不同总距角下的共轴旋翼系统中上、下旋翼的非定常气动特性以及影响规律。计算结果表明:上旋翼桨叶的桨尖涡会直接与下旋翼桨叶发生碰撞,导致下旋翼桨叶拉力损失;上旋翼桨叶的桨尖涡和下旋翼桨叶的桨尖涡相互干扰,改变了桨尖涡的强度和轨迹;上、下旋翼桨叶相互靠近时,上、下旋翼桨叶的拉力均会上升,之后相互远离时上、下旋翼桨叶拉力均会先下降再上升。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

平板空气舵干扰区烧蚀凹陷对热环境影响

针对高速飞行条件下空气舵干扰区烧蚀产生的局部凹陷对气动加热的影响问题,建立了平板-空气舵流动模型,针对典型高速飞行状态,采用高温热化学非平衡数值模拟,研究了空气舵缝隙区的流动结构和气动加热规律,并对舵缝干扰区的烧蚀外形进行了模化,分析了干扰区烧蚀凹陷对流动结构和气动加热的影响,结果表明:烧蚀凹陷改变了干扰区压力分布规律...     

轴承腔内壁与油膜换热的数值模拟与试验

航空发动机后轴承腔内壁与滑油的换热分析是轴承腔热防护结构设计的基础。对航空发动机轴承腔内壁换热模拟试验件开展试验与数值模拟研究,得到了滑油油膜对轴承腔内壁的换热影响。通过测量试验件外壁面、内壁面以及滑油油膜的温度得到了试验件内壁面换热热流密度与换热系数的分布;结合CLSVOF(Coupled Level Set and Volume Of Fluid)油/气两相流以及热-流-固耦合计算方法对试验件进行了换热分析,并将内壁对滑油的局部热流量的计算结果与试验结果进行了对比,结果显示两者在各个工况下均吻合较好。通过将局部换热系数计算值与当地的滑油流动雷诺数Rel进行对比分析,结果显示内壁局部努赛尔数Nuw与Rel的0.7次方成正比关系。另外,对转速对换热的影响进行分析得到Nuw与旋转雷诺数Re_(rot)的0.345次方呈正比关系。     

基于对流换热的低温空化流动数值模拟研究

为了考虑低温介质的热力学效应对空化发展的影响,基于气泡表面对流换热平衡建立了温降与气泡生长的关系,引入夹带理论估计了对流换热系数,并对一种输运型空化模型进行修正,将修正后的空化模型以二次开发的形式嵌入至商业软件中,同时引入能量方程源项以及物性参数随温度变化关系,对二维翼型表面空化流动进行数值仿真,通过与实验结果的对比,发现计算结果与实验结果符合较好,修正后的空化模型能够更好地预测空化区内温度的分布,最大温降偏差由62.18%降低至7.14%,平均温度偏差由0.59%降低至0.28%。考虑热效应之后,空化区主要由气液混合组成,来自主流的液体一部分经对流传递至空化区,一部分在翼型头部和气液界面处发生空化形成蒸汽,导致空化区气相体积分数显著减小,空化区与主流之间的界面变得模糊,最大温降和压降均发生在翼型头部位置。