旋转状态下气膜冷却模型的数值模拟

通过对带有90°倾角圆柱形交错孔排的涡轮叶片模型进行数值模拟,得到了不同主流雷诺数、旋转数和吹风比情况下前缘面与后缘面侧的气膜冷却流动与换热特性及各气膜孔流量系数的分配规律.结果表明,冷气受到离心力与哥氏力的共同作用向高半径处发生偏转,导致壁面冷却效率降低;雷诺数的增大会削弱气膜冷却效果,高吹风比则不利于气膜孔下游区域的冷却.各气膜孔的流量系数随吹风比的增大而增大,随旋转数的提高而减小.在后缘面侧,相同工况下各气膜孔的流量系数明显高于前缘面侧对应气膜孔的值.      

上游双扰动对圆柱气动特性影响的数值模拟

通过用商用计算流体力学软件FLUENT 6.0中的SIMPLE方法求解二维非定常层流模型,对雷诺数Re=200时上游双扰动体对下游圆柱气动特性影响进行了数值模拟研究.计算结果表明当双扰动体相隔较近(H/d=1.1,1.33)时,其后形成单一的Karman涡街,它们对下游圆柱的干扰作用类似于单扰动体时的情况,但是流动结构发生变化的临界间距Lcri较之有效宽度相同的单扰动体的要大,因此可以推断出在该种布置下双扰动体和圆柱绕流会出现2种流动模式.在空穴流动模式下面,下游圆柱的阻力和升力脉动最多减小94%和79%,系统阻力(包括圆柱和扰动体)也能减小33%.而当双扰动体尾流为双Karman涡街时(双扰动体相隔较远H/d=3.33,5.0),圆柱迎风面部分区域仍受到来流的冲击作用,气动力减小有限.     

微通道直角突缩弯道气体流阻特性

采用有限体积方法,研究了矩形截面微通道直角突缩弯道局部流场及压力损失机理,重点考虑了雷诺数和弯道进出口截面面积比对局部流场及压力损失的影响,雷诺数范围为1~200,特征尺寸为500 μm.研究表明,低雷诺数时压力损失主要由壁面摩擦产生,弯道压力损失系数与具有相同中心线长度的直通道的压力损失系数变化趋势相同,但在数值上较小,突缩弯道会引起弯道内速度型发生变化,影响弯道压力损失;而高雷诺数时,压力损失主要由弯道内角位置的速度分离及漩涡产生,突缩弯道会抑制速度分离及漩涡产生,进而降低弯道压力损失.根据计算结果,拟合了弯道压力损失系数与雷诺数及进出口截面面积比的关系式.     

宇宙线福布什下降的数值模拟

本文从银河宇宙线的太阳调制方程出发,认为激波对宇宙线的影响是由一扰动区产生的。在此扰动区中太阳风速度增加,扩散系数下降。由此进行了数值模拟,模拟结果表明:扩散系数的下降在产生福布什下降时要比太阳风速度的增加更为有效;福布什下降在近日球层内向外传播时,其幅度随径向距离的增大而衰减;两个无相互作用的激波同时存在时所产生的福布什下降为每个激波单独存在时的福布什下降的简单线性迭加。      

涡轮叶片尾缘复合通道隔板结构对流动的影响

涡轮叶片尾缘隔板结构的优化可以强化换热,其机理是破坏附面层的形成,并形成脉动效果.本试验通过流道中打入粒子,在激光片光源的照射下使用高速摄像头拍摄涡轮叶片尾缘带隔板的复合通道内的流动形式.主要试验了2种不同结构的隔板(矩形隔板与120°波形隔板)在射流孔全部打开和单个射流孔打开时的流动形式,以及波形隔板在不同条件下对流动形式的影响.结果表明:矩形隔板和波形隔板都对射流有所影响;波形隔板射流有其特有的震荡效果,在固定几何条件下,主要与进口雷诺数相关.     

矩形通道内气膜出流对内换热的影响规律

研究对象为矩形通道,包括2个直肋和9个气膜孔,主要研究气膜的出流比、开孔率对通道内换热的影响,温度场的测量采用了热色液晶测温技术.实验发现气膜出流对换热有强化作用,气膜孔对通道平均换热影响也与其在肋间的位置有关,其中肋后孔的影响最大.在不同的雷诺数下,开一个气膜孔时,出流比在6%~8%之间通道壁面内换热效果最好,低雷诺数下增强换热的效果显著;研究孔边的换热规律发现,在孔下游换热明显增强,从孔边到下游5倍孔直径处,强化换热比在1.1倍以上;此外还研究了开孔率对增强换热的影响规律.      

带气膜出流叶片前腔内表面冲击换热数值模拟

采用数值模拟方法对不同雷诺数下静止状态涡轮叶片前腔带气膜孔出流的冲击流动与换热特性进行了研究.分析了叶片前缘冲击流动产生的不同涡团对其内表面换热的作用机理.计算结果表明:相同雷诺数下,叶片前缘内表面气膜孔附近的换热强化比高于通道的平均值.随着雷诺数增加,换热强化比有所提高.冲击流动与通道流动耦合而形成的波浪形涡区,极大地扩展了冲击强化换热区域.气膜孔出流的抽吸作用对冲击流产生影响,进一步扩大了冷却空气在前缘内表面的覆盖范围.气膜腔叶根处纵向截面的涡团阻碍了冷气向叶根方向扩展,降低了冷却效率;而横向截面的涡团则促进冷气与壁面热气的掺混,提升了换热效果.      

中高层大气对低层大气脉冲扰动的响应

从基本的大气运动方程出发,采用二维全隐欧拉格式(FICE)建立极坐标下二维声重波传播的数值模式.对小振幅高斯形态声重波传播过程的模拟结果表明,该数值模式能够再现声重波波包向上的稳定传播过程.应用此数值模式模拟了在没有背景风场和存在背景风场条件下中高层大气对10km高空处近似脉冲点源扰动的响应过程,获得了中高层大气对低层大气脉冲扰动响应的图像.结果表明,电离层高度对重力波扰动的响应不但出现在激发源激发之后,还与激发源具有较大距离,扰动的幅度随高度的增大而增大;有背景风时,顺风情况下重力波波动传播路径较逆风情况下平缓,且传播的水平距离比逆风情况下远,波动的振幅较逆风情况下弱,逆向背景风场能够加速重力波向上传播.同时,模拟还表明脉冲扰动正上方可以激发产生一种以6 min为主要周期,具有声学分支特性的快速波动.本文模拟的扰动结果与2004年12月26日苏门答腊-安达曼地震引起的电离层TEC扰动具有类似的图像.      

近地微电离大气中阿尔芬波的截止

利用全电流双流体模型分析了地面磁扰动信号在近地层大气中的传播.结果表明磁扰动信号的传播模式受中性粒子elax-elax离子碰撞频率、磁扰动频率、磁粘滞系数以及大气离化系数等参数的影响.中性粒子碰撞效应导致近地层大气中不存在阿尔芬波传播通道.当磁扰动频率远小于0.001Hz时,磁扰动传播存在趋肤效应,且在近地大气中从高往下逐渐增强;当磁扰动频率远大于0.001Hz时,磁扰动以类似于真空电磁波的模式在近地大气中传播.      

来流局部干扰对三角翼气动特性的影响

通过在三角翼上游加入干扰圆柱的风洞实验方法,研究了来流干扰对微小型飞行器MAV(Micro Air Vehicle)气动特性的影响.结果表明,在刚性和弹性三角翼顶点上游加入圆柱干扰时,两者均出现缓失速,刚性翼产生缓失速与干扰圆柱尾流关系密切,弹性翼的缓失速不仅与此有关,还与弹性翼的振动有关.无干扰或在机翼顶点加入干扰时,在攻角为4°~18°内弹性翼的升力系数比刚性翼的要大,但升阻比相对要小.由于弹性翼的振动与机翼绕流结构、气动力之间的耦合,弹性翼顶点与翼尖振动的主频随着攻角增大呈规律性的变化,失速攻角附近翼尖的振动主频是其涡脱落频率.      

吸气式高超声速飞行器冷流试验设计及验证

对于吸气式飞行器而言,地面冷流试验是检验其进气道性能及气动特性的一项重要手段.以二元混压式进气道、机体/推进系统耦合为基本特征,设计了采用超燃冲压发动机为推进系统的内外流一体化巡航飞行器,针对其高超声速特性开展了冷流风洞试验,来流速度范围Ma=5.0~7.0,攻角范围α=-4°~8°.测压试验结果表明,随着来流马赫数的增大,进气道的总压恢复系数下降;而流量系数先上升,在设计点达到最大值;在一定攻角范围内,进气道的总压恢复系数和流量系数提高,但当攻角增大至巡航攻角时,随着攻角的增大,进气道的总压恢复系数和流量系数逐渐下降.测力试验验证了数值算法的有效性,除轴向力系数以外,其余气动特性系数的发展规律及数值基本吻合,可通过修正试验值的方式外推出飞行器的气动特性数据.     

降落伞开伞过程的多结点模型仿真

根据降落伞的结构和力学特征,在轴对称假设下创建了伞衣及回收物系统的多结点结构模型.通过考虑应力,重力和气动力的作用效果,建立了用于无气流攻角平面圆形降落伞充气模拟的多结点结构模型动力学方程组.对开伞过程中的流场变化引入准定常假设,利用simple算法数值模拟求解RNG(Renormalization Group)k-ε湍流模型下的雷诺平均N-S(Navier-Stokes)方程以获得选定时刻的伞衣表面压力分布.结合多结点模型动力学方程组的解算代码和计算流体力学程序,采用流固耦合的方法对选定的平面圆形降落伞模型的开伞过程进行了动态仿真,得到了开伞过程中降落伞外形和特性的变化.通过结果分析和比较,证明了多结点模型的可行性,发展出了一种用于降落伞流固耦合计算的新方法.      

几何和流动参数对上翘后体阻力系数的影响

利用N-S方程对9个上翘后体模型进行了气动力计算,主要研究后体几何参数和流动参数对上翘后体阻力系数的影响.研究结果表明, 后体的压差阻力系数分别随上翘角、收缩比的增加及迎角的减小而明显增加;后体摩擦阻力系数分别随后体的长细比的增加和雷诺数的减小而增加;后体越扁平,其压差阻力系数越大;在跨音速时,波阻系数也与上翘角有关,上翘角增加会导致波阻系数进一步加大.     

超临界层流翼型优化设计策略

针对超临界层流翼型设计问题,提出一种两轮优化策略。采用γ-Reθt转捩模型耦合剪切应力输运(SST)模式的湍流模型对翼型边界层转捩进行预测。翼型几何参数化建模采用形状分类函数转换(CST)方法,设计变量为描述翼型几何特征的参数。第1轮优化的目的是尽量提高层流区域的比例,气动分析模型为基于Kriging模型的代理模型,优化算法为遗传算法,通过优化获得满足约束要求的层流翼型。第2轮优化目的是对第1轮优化获得的翼型进行微调,进一步提高翼型的升阻比,气动分析直接采用CFD程序,优化算法采用基于梯度的优化算法。算例表明,应用本文提出的两轮优化策略,可将超临界翼型NASA SC(2)0412优化设计成超临界层流翼型,翼型的上下表面层流区比例分别达到了55.5%和47.0%,升阻比提高了38.1%。     

后掠角对后掠机翼边界层稳定性及转捩的影响

后掠机翼边界层流动稳定性及转捩对翼型的设计及优化有着重要的参考价值,而机翼后掠角是引起后掠机翼边界层横流失稳的关键参数之一.以NACA0012翼型为研究对象,通过求解三维可压缩Navier-Stokes方程计算了展向无限长后掠机翼的基本流场;通过求解Orr-Sommerfeld方程得到了扰动Tollmien-Schishting波演化的中性曲线及幅值曲线,研究了后掠角对后掠机翼边界层流动稳定性的影响;最后采用eN方法进行了转捩预测.研究发现,随后掠角的增大,横流强度和扰动幅值放大指数n均先增加后减小,且后掠角在40°~50°之间横流强度达到最大值.当后掠角在50°左右时,用转捩预测eN方法计算的幅值增长指数N值最大,导致转捩发生所需的初始扰动幅值最小,转捩最易发生.      

前体尾流对降落伞工作性能的影响

为研究前体尾流对降落伞工作性能的非定常影响,基于Realizablek-ε湍流模型采用PISO算法开展了物伞系统的非定常绕流数值计算,获得了精细的流场旋涡结构。在此基础上,研究了不同拖曳比下物伞系统的尾涡演变规律、流场分布规律以及伞衣气动特性变化。结果表明:前体尾涡导致伞衣入口处的涡量大小和方向时刻变化,随拖曳比增加,涡量黏性耗散增强,进入伞衣的旋涡强度逐渐减弱,伞衣入口形成稳定的负涡量区,伞衣尾涡脱离周期随之延长;拖曳比对尾涡区后端（伞衣入口处）流场压力的影响远大于前端,随拖曳比增加,流动形式逐渐由闭式转变为开式,流场的速度分布和压力分布更为对称,伞衣入口形成稳定的正压区,内外压差增加;当拖曳比大于9时,前体尾流对降落伞阻力系数和表面压强系数的影响减小。      

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

类X-51A飞行器非定常湍流精细模拟

针对类X-51A飞行器在超声速大迎角状态下存在的大范围非定常分离流动,开展了精细化湍流数值模拟研究。计算基于高阶格式下的延迟分离涡模拟方法（DDES）,来流马赫数为2.5,迎角为10°。分析了该复杂流场中存在的分离流动现象、分离流动诱导的气动特性变化规律以及压力脉动特点;其中重点研究了壁面压力脉动强度分布情况和监测点压力脉动频谱特性。分析结果表明:飞行器大迎角飞行时从侧缘诱导出明显的分离涡,并对尾部舵面产生干扰;受干扰尾舵表现出明显的非线性及非定常气动特性;分离涡的存在导致飞行器尾舵前缘等位置的壁面压力脉动显著增强,200~300 Hz的低频高幅值脉动可能会导致结构破坏。