多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

高超声速飞行体亚密湍流尾迹RCS特性的相似规律研究

再入或临近空间高超声速飞行器在大气层中飞行时,会形成等离子体绕流和尾迹,将严重影响飞行器的电磁散射特性,从而使得临近空间高超声速飞行器表现出与大气层内常规飞行器和大气层外轨道飞行器在电磁散射特性上的明显不同,对飞行器的探测、跟踪带来新的挑战。弹道靶是高超声速飞行器尾迹电磁散射特性研究的重要设备,为了由弹道靶试验结果外推出真实飞行条件下亚密湍流尾迹的RCS,本文从高超声速流场模拟的双缩尺率和亚密湍流尾迹RCS模拟的Born近似出发,推导出真实飞行条件下和地面弹道靶试验之间亚密湍流尾迹RCS模拟的一种相似规律,推导结果显示当D1/D2=η,V∞1=V∞2,ρ∞1/ρ∞2=1/η,ω1/ω2=λ2/λ1=1/η时,雷达散射截面服从σ1/σ2=η4的相似规律。推导出的相似规律可以为超高速飞行器亚密湍流尾迹RCS研究提供技术支持。     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

超声速条件下亚毫米液滴的变形破碎模态

以超声速气流中液滴变形破碎行为为研究内容,对水平激波管中承受激波冲击的亚毫米水液滴(0.44~1.09 mm)变形破碎过程进行了观测,实验激波马赫数范围为1.07~2.11。利用纹影法,结合高分辨率高速相机对不同破碎模态下液滴的变形破碎特征进行了记录,得到了袋状、多模态、剪切和爆炸式等破碎模式下的液滴纹影图像,分析了液滴运动参数的时空关系。得出了液滴变形阶段,液滴无量纲横向变形宽度以及液滴无量纲迎风面位移随无量纲时间的变化发展规律,并且得出在液滴初始直径相同时,不同液滴破碎模式的无量纲最大横向变形宽度的变化,其中袋状、多模态、剪切破碎模式的无量纲横向最大变形宽度均在1.15~1.61范围内变化,爆炸式破碎模式的无量纲横向最大变形宽度均在0.21~0.68范围内变化。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

边界层抽吸对超燃冲压发动机流场特性的影响研究

为了解边界层抽吸对超燃冲压发动机流场的影响,采用风洞试验和数值计算对隔离段激波串特性以及燃烧室燃烧特性进行了研究。结果表明,在发动机入口马赫数2.0,总温950K,总压0.82MPa的来流条件下,当量比为0.18先锋氢气与不同当量比煤油共同燃烧呈不稳定状态,激波串在隔离段内前后振荡传播。当煤油当量比为0.29时,激波串振荡前缘远离抽吸位置,边界层抽吸对发动机流场基本没有影响。随着煤油当量比逐渐增大,激波串前缘位置到达抽吸区附近,边界层抽吸开始产生影响,改变了隔离段内的激波串动态演化过程、形态结构以及位置分布,同时有效提高了隔离段抗反压特性,使得煤油最大当量比可以由0.38增大至0.42。此外,边界层抽吸对发动机内的亚燃/超燃区域分布也会产生影响。     

凹腔/支板结构亚燃冲压燃烧室性能

为了避免基于凹腔火焰稳定器的亚燃冲压燃烧室壁面喷注时燃料与主流空气掺混非均匀性问题和提高燃烧室的性能,提出在亚燃冲压燃烧室中使用支板喷注代替壁面喷注的方案,数值模拟了凹腔/支板结构亚燃冲压燃烧室中燃料分布及流场结构,并分析了支板结构对燃料空气混合及燃烧室性能的影响。研究表明:支板虽然使燃烧室出口的总压恢复系数相对于壁面喷注方式下的降低了63%,但能使燃料均匀分布于整个流道内,增强了燃料与空气掺混,使燃烧室出口的混合效率和燃烧效率分别提高了21.4%和20.5%。燃烧效率的提高弥补了采用支板导致的燃烧室内气流的额外总压损失所带来的机械能损失,使得支板喷注时燃烧室出口的比冲提高了39.6%。因此,在亚燃冲压燃烧室中设置凹腔/支板结构,有利于提高燃烧室整体性能。      

亚燃模态下释热分布对发动机性能的影响

为揭示释热分布对亚燃模态下超燃发动机性能的影响规律,对马赫数2.0,总温1100K,总压1.0MPa的来流,完成了燃烧室直连式实验和数值模拟研究。针对上游释热和下游释热条件,对比了发动机壁面压力分布,分析了内流道一维质量加权马赫数,获得了发动机部件和总体性能数据。结果表明,在总当量比相同情况下,上游释热能够获得更好的发动机性能,没有尾喷管时比下游释热获得的推力高出约18%,但在有尾喷管时只相差2.6%;对于本文构型,燃烧室和尾喷管是发动机推力的主要来源,两种释热分布下,二者产生的推力超过了发动机总推力的90%;但对于更高总当量比,上游释热可能会导致进气道不启动,需要增加下游释热获得更高的发动机性能。     

叶片丢失后发动机整机响应模拟试验与仿真

依据大涵道比涡扇发动机结构设计了叶片丢失试验台,开展了一系列模拟叶片丢失试验,并采用显式有限元方法进行了数值仿真,研究了发动机叶片丢失后整机结构响应与载荷传递规律。结果表明:丢失叶片与机匣存在叶尖与叶身两次撞击,对应的加速度曲线存在两个响应峰值;转子转速越高,加速度响应幅值越大。叶片飞断后转子不平衡载荷传递路径为前轴承支承-中轴承支承-中介支板-机匣结构;叶片撞击机匣导致的冲击载荷则由风扇机匣向后传递,最终传给吊装结构;发动机承受的载荷是由不平衡和冲击影响耦合得到,其中冲击载荷为主要部分。该研究为掌握真实发动机叶片丢失下整机响应规律提供了试验模拟方法与数值仿真分析工具。      

航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究

为深入了解真实航空发动机内燃烧流场,采用动态亚网格模型结合单步快速化学反应(FC)、火焰面(FLM)及反应进度变量(FPV)等三种燃烧模型对径向两级反向旋流燃烧室单头部构型进行了大涡模拟。为避免模型简化误差,对燃烧室包括全部气膜冷却孔在内的精细结构进行了完全仿真。在已达到统计定常状态的冷态流场基础上首先模拟了燃料喷注和掺混过程,约2.6ms后燃料到达真实的点火位置,随后采用FPV模型在半径3mm的球形区域数值模拟了点火,展示了在主燃孔横向射流作用下初始火焰沿化学恰当比混合分数等值线传播并充满整个火焰筒的发展过程,结果显示火焰到达燃烧室出口的耗时约为6～7ms。不同模型算法预测的平均温度场与CARS测量结果作了对比,LES-FPV,RANS-FPV,LES-FLM以及参考文献中RANS-FLM计算平均误差分别为3.47%,4.17%,7.76%和11.22%,表明LES改进了模拟精度,且FPV模型显著优于FLM模型。RANS-FPV预测的出口存在严重热斑,导致其给出的出口温度分布因子(OTDF)及最大径向温度分布因子(RTDF)值分别达到0.593和0.313;LES-FPV结果均匀性最好,其预测值分别为0.284和0.193。