MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

液氢绕回转体非定常空化的流动特性

为研究航空发动机内部低温燃料的非定常空化流动特性,采用大涡模拟(LES)对液氢绕回转体的非定常流动特性进行了分析。通过与实验结果的对比可知,所采用的数值计算方法能够有效地模拟液氢绕回转体流动的非定常过程。讨论了空化演化过程和旋涡动力特性,研究结果表明:空化的非定常演变过程大致可以分为3个阶段:附着空穴生长阶段、大尺度空泡发展阶段和小尺度空泡发展阶段,且反向射流是造成空穴不稳定及脱落的主要原因;探讨了空化与旋涡之间的交互作用,明确了空穴前端和尾端处的涡量源于旋涡伸长项,空穴内部涡量源于旋涡扩张项,空穴交界面、反向射流头部区域涡量源于斜压扭矩项。      

大截面Z型材滚弯成形缺陷的数值模拟及试验研究

以控制大截面Z型材四轴滚弯成形中型材局部起皱、扭曲等质量缺陷为目的,利用ABAQUS有限元软件数值模拟大截面Z型材数控四轴滚弯过程,在分析滚弯成形缺陷产生原因的基础上,提出了侧滚轮非对称加载和二次滚弯方式控制成形缺陷的方法。结果表明:二次滚弯方式对起皱缺陷具有较好的改善作用;侧滚轮非对称式加载方式(左侧滚轮高于右侧滚轮),对起皱和扭曲等缺陷具有明显的抑制作用。滚弯试验结果证明,型材滚弯成形的扭曲量相比数值模拟结果误差在10%以内。     

钛合金宽弦风扇叶片的振动特性

通过模态仿真和模态试验技术分别对某型号TC4钛合金宽弦风扇叶片的振动特性进行研究。采用ANSYS-APDL软件对风扇叶片模型进行模态仿真分析。根据风扇叶片的坎贝尔图以及不同模态的振型分布,识别风扇叶片的弯扭振型。风扇叶片不同转速不同模态下的转速裕度和频率裕度分析显示:第1阶弯曲模态在红线转速下的转速裕度和频率裕度超过20%,表明风扇叶片满足共振裕度要求。在三轴向振动台上利用扫描式激光测振仪,采用敲击法和随机带宽激励的方法测量叶片的响应,分析了风扇叶片不同模态的冲击响应特性和多轴响应特性,发现风扇叶片在周向和轴向激励下的响应以低阶弯曲模态为主,而在径向激励下的响应以扭转和复合模态为主,不同激励位置不同模态的响应特性与风扇叶片的振型息息相关。      

热力学效应对液氢非定常空化流动的影响

为研究热力学效应对低温流体非定常空化流动的影响,以液氢为研究对象,采用大涡模拟计算了热力学和等温条件下液氢绕回转体的非定常空化流动。计算结果表明:热力学效应延长了液氢的空化周期,增强了空化的非定常特性,抑制了空化的发展。热力学条件下的空泡团内部包含更细、更小的气泡,呈现出多孔质特征,且使整个流域存在1.5K左右的温度波动。通过旋涡运动分析发现,热力学效应使旋涡结构由空穴交界面向空穴内部移动。基于涡量传输方程分析了热力学和等温条件下空穴和旋涡的交互作用,发现在热力学条件下,旋涡伸长项的作用位置主要位于空穴前端和尾端的交界面处,旋涡扩张项和斜压扭矩项主要位于空穴内部;在等温条件下,旋涡伸长项和旋涡扩张项主要位于空穴上方交界面处,斜压扭矩项主要位于空穴尾端。      

不同进口条件对OGV/ 前置扩压器流场大涡 模拟结果的影响

为了研究不同进口条件对压气机出口导向叶片(Outlet Guide Vane,OGV)/前置扩压器流场性能的影响,分别采用合成涡方法与白噪声方法生成进口条件,研究不同进口条件对流场数值模拟结果的影响,并与文献给出的试验结果进行对比。结果表明:采用合成涡方法生成的进口速度场满足试验测得的时均速度分布与脉动均方根分布,且其生成的进口湍流结构可传播到下游较远处;合成涡方法能够较好地预测OGV内的流场,与白噪声方法相比,所预测的OGV吸力面分离区较小;在OGV出口截面、扩压器内和扩压器出口截面,采用2种方法得到的速度分布相差不大,与试验相比分布趋势相同。     

鸟撞工况下HCA动态拓扑优化整级叶片的数值模拟

航空领域中飞机发动机在满足抗鸟撞性能的前提下,实现轻量化是关键问题之一。针对航空发动机冷端风扇叶片,基于HCA算法的动态拓扑优化方法对叶片实现进一步减重,为了验证该动态优化方法的有效性,建立了鸟撞航空发动机整级叶片冲击动力学有限元模型,模拟受气动与离心载荷作用下高速稳定旋转的发动机风扇叶片遭受鸟体撞击的瞬态响应过程。基于LS-DYNA软件平台,对考虑了鸟撞的多工况、多约束条件下叶片的动态优化结果与优化前叶片的多项动态响应指标进行了对比分析,证明了HCA动态拓扑优化方法比传统叶片轻量化方法更为优秀,在满足适航条例强度要求的同时可使叶片减重比达到37.9%,论证了在叶片轻量化设计上基于HCA算法的动态拓扑优化方法的可行性与优越性。     

推力矢量对飞机大迎角动态气动特性的影响

推力矢量是提高战斗机大迎角动态气动特性,提升其过失速机动能力和飞行品质的重要手段。新一代战斗机的高机动性要求也使气动和推力矢量的融合控制研究日益重要。针对中国空气动力研究与发展中心∅3.2 m开口低速风洞,研制了喷流模拟器和通气动态试验装置,建立了带推力矢量的大迎角动态试验技术。开展了不同减缩频率、不同落压比、不同喷管偏角时的大迎角俯仰振荡运动特性试验研究。结果表明：与无喷流试验相比,带喷流时模型的动态特性均随着落压比和喷管偏角的变化呈现规律性的变化;力和力矩系数形成的迟滞曲线面积随着落压比和偏角的增加而增加;减缩频率的变化对模型的动态特性影响小于无喷流时的影响。总的来说,推力矢量的影响未改变模型大迎角动态特性的基本规律,但是随着推力矢量角度和大小的变化,有规律地改变了模型动态气动力和力矩的变化幅度。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

燃气-液压混合驱动大惯量起竖系统内弹道建模与参数优化

针对传统起竖方式起竖缓慢、起竖动力源功率密度小等问题,基于燃气与液压混合驱动的起竖方案,设计了一种变截面的固体药柱。基于经典零维内弹道方程,考虑燃气气-粒两相流影响及点火药、固体装药两种燃气的共同作用,结合起竖机构的运动规律,以大惯量起竖装置为研究对象,建立了起竖系统的一体化动力学模型。参照目标起竖角度曲线,采用遗传算法对药柱结构参数进行匹配,并对燃气驱动方案进行优化,最终确定了前端直径0.1m,后端直径0.04m的变截面药柱及药柱燃烧1.5s,燃气驱动4.62s的起竖方案。正向内弹道计算结果表明,基于精细化模型的内弹道优化,可使系统在较短时间内(4.62s)以较小的过载(小于0.35g)起竖至目标角度(1.28rad)。