壁面质量引射对高速飞行器减阻降热影响的研究

针对未来新型高速飞行器高升阻比和低热载荷的发展需求，本文选取贴近工程实际应用的三角翼为研究对象，开展了飞行状态下大面积壁面质量引射对减阻降热影响的数值模拟研究，对比分析了壁面质量引射在不同的飞行高度、来流马赫数、引射流量、引射温度和气体模型条件下对三角翼气动力热特性的影响。从高速飞行器减阻降热的工程需求出发，指出了质量引射流动控制在工程应用中存在的问题与后续应重点关注的方向。     

一种背负式无附面层隔道进气道的数值模拟研究与实验验证

采用CFD方法对背负式无隔道进气道/机身一体化流场进行了研究。主要分析了机身上表面附面层的发展情况、进气道进口鼓包排除附面层气流的特性以及进气道内部的流动特征,并将所得到的结果与实验数据进行了对比,比较了两种不同网格的计算准确度。研究发现进口鼓包能够有效地隔除机身上表面的附面层气流,进口段横向压力梯度是导致附面层气流"溢出"进气道的主要驱动力,另外进气道的流量系数对排除附面层气流的效果有着显著影响。     

“蚌式”进气道附面层扫除特性风洞试验

以某低、高速风洞为平台,设计搭建了“蚌式”进气道附面层扫除特性测量试验系统,进行了不同流量系数和来流马赫数下进气道鼓包表面附面层扫除特性的风洞试验,通过对试验数据的整理、计算和对比分析同型号的飞行试验结果,研究了“蚌式”进气道鼓包表面附面层扫除特性。研究结果表明:在相同的来流马赫数下,随着流量系数的增大,鼓包表面附面层的扫除能力逐渐减弱;在亚声速工况的绝大多数流量范围内,鼓包表面压力系数沿鼓包中心线对称分布、压力梯度变化明显,且在不同截面沿主流方向具有增大的特征,鼓包构型对附面层扫除效果较强;超声速工况下具有明显附面层扫除能力的流量范围明显小于亚声速工况,进气道唇口形成的弓形激波是影响鼓包表面不同位置压力梯度变化的主要因素,进而决定着附面层扫除特性。在接近来流马赫数18及以上飞行工况下,附面层的扫除能力减弱,附面层分离加强,进而会造成较大的进气压力损失和畸变。     

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

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

波瓣引射-混合器特性参数影响的数值研究和验证

运用不可压缩流动Navier—Stokes方程,对波瓣喷管引射-混合器的流场和引射特性进行了三维数值研究。在计算过程中,主流进口采用质量流量边界条件,二次流进口和混合流出口采用压力边界条件,均设置为环境大气压力。与相关的实验结果对比表明本文的计算方法可以有效地预测引射流量比和混合流场。针对混合管的结构参数开展了系列研究,获得了混合管截面比和长径比对于引射系数和热混合效率的影响趋势。     

来流条件对叶片表面附面层转捩影响的数值模拟研究

通过对二维压气机叶栅叶片表面附面层流态的数值模拟分析,研究了在不同来流马赫数和冲角下,叶片表面层流附面层向紊流附面层的转捩过程,转捩区域的大小等。归纳分析了部分来流条件对转捩过程的影响机制。结果表明:来流马赫数和冲角对转捩起始点和结束点的位置,即转捩所占区域的大小有重要影响。研究工作取得的进展不仅为进一步从机理上深刻认识来流条件对叶片表面附面层流态的影响机理提供了一定的科学依据,而且对今后指导叶片表面附面层流态模似实验工作的开展具有重要参考价值。      

二维引射—混合器流场的数值研究与验证

本文运用一维引射特性方程和不可压缩流动Navier-Stokes方程,对二维引射—混合流场进行了数值研究。首先用一维引射特性方程估算理想的引射气流流量;然后用N-S方程对混合管内的气流混合流动进行计算,用计算得到的混合管速度分布和壁面粘性剪切力来修正一维引射特性方程中的经验常数,逐次逼近真实的引射流量比。计算结果表明,本文的数值方法可以有效地预测引射流量比和混合流场。     

质量引射对边界层稳定性的影响

采用直接数值模拟（DNS）方法和线性稳定性理论(LST)研究了质量引射对马赫数为6的钝板边界层中第二模态扰动演化的影响。研究了质量引射及扰动的不同因素对扰动演化的影响,包括引射速度幅值、引射宽度、位置、组合以及体积流量。研究结果表明:质量引射对二维和三维扰动波都起不稳定作用,当引射位置接近扰动波中性曲线的下支界时,不稳定作用最明显;体积流量是决定质量引射对边界层作用的重要因素;当体积流量较小时,质量引射组合放置能起到单个质量引射叠加的效果;LST预测得到的N值修正与直接数值模拟的结果相近,可以定量地反映质量引射对边界层流动稳定性的影响。     

来流条件对叶片表面附面层转捩 影响的数值模拟研究

通过对二维压气机叶栅叶片表面附面层流态的数值模拟分析,研究了在不同来流马赫数和冲角下,叶片表面层流附面层向紊流附面层的转捩过程,转捩区域的大小等。归纳分析了部分来流条件对转捩过程的影响机制。结果表明来流马赫数和冲角对转捩起始点和结束点的位置,即转捩所占区域的大小有重要影响。研究工作取得的进展不仅为进一步从机理上深刻认识来流条件对叶片表面附面层流态的影响机理提供了一定的科学依据,而且对今后指导叶片表面附面层流态模似实验工作的开展具有重要参考价值。     

平衡气体对乘波体气动力热特性影响

针对乘波体在高空以高马赫数飞行时的特点,给出了一种考虑高温平衡气体效应的乘波体设计方法,并通过数值模拟的方法研究了高空高马赫数飞行条件下高温平衡气体对乘波体气动特性的影响。研究结果表明,相对于完全气体,平衡气体对乘波体的气动力特性影响不大,但对乘波体的前缘驻点区和下表面的热流及辐射平衡温度有较大的影响。另外在乘波体较为适宜的飞行攻角范围内,攻角的增大不会显著改变平衡气体效应对乘波体气动力特性的影响,也不会显著改变平衡气体效应对乘波体前缘驻点区的热流与辐射平衡温度的影响,但却会在一定程度上减小平衡气体效应对乘波体下表面的气动热及辐射平衡温度的影响区域。     

射流对压气机叶栅分离流控制的数值研究

现代先进轴流压气机级负荷不断提高的发展趋势导致流动分离日益严重.借助数值模拟分别对非定常射流和定常射流进行了参数优化研究.结果表明:基于射流的主动流动控制能有效弱化或消除流动分离,不同射流方式存在不同的最优射流参数(射流方向、位置、速度和频率等),这就为利用射流控制轴流压气机分离流动的工程应用奠定了一定的理论基础.     

前缘喷淋射流对导叶压力面冷却特性的影响

为了评估涡轮导叶的前缘喷淋射流对压力面多排气膜孔冷却特性的影响,在高亚声速风洞中进行了实验,获得了有无前缘喷淋射流时叶片表面的气膜冷却效率和传热系数。叶栅进口雷诺数(基于叶片弦长)范围为2.0×105~4.0×105,出口等熵马赫数为0.95,叶片前缘和压力面分别都包含6排圆形孔,质量流量比的范围分别为2.46%~4.57%和2.00%~3.71%。实验结果表明:在没有前缘喷淋射流时,压力面前半段的气膜冷却效率受质量流量比的影响较小,而后半段的气膜冷却效率随着质量流量比升高而增大。前缘喷淋射流使压力面多排气膜孔的冷却效率提高了20%~70%,并且使气膜冷却效率沿流向分布更均匀。不论是否有前缘喷淋射流,压力面的传热系数比都随质量流量比升高而增大,沿流向看,前缘喷淋射流提高了压力面前缘和尾缘区域的传热系数比而对压力面中间区域影响较小。     

加热表面水珠运动特性研究

为研究加热表面的水珠运动特性,提出了加热表面水珠的几何参数、受力及其运动过程的计算方法。试验获得了水珠的表面阻滞力、黏性阻力和气动力计算关系式中的相关系数,给出了不同风速条件下水珠运动的临界直径,进行了加热表面水珠运动试验并对其过程进行了数值计算。试验结果表明：水珠的无量纲表面阻滞力保持恒定,运动时受到的黏性阻力与其运动速度和宽度有关,在外流场作用下所受的气动力可采用修正圆球阻力公式进行计算。将试验获得的相关系数加入水珠运动模型中,对加热表面水珠运动过程进行数值计算,计算结果与试验结果吻合,说明本文所述方法能够准确地模拟加热表面水珠大小和位置随时间的变化。     

Coating deposition regularity depended on orientation difference in PS-PVD plasma jet

The YSZ coatings are prepared by the plasma spray-physical vapor deposition (PS-PVD) technology based on a specific experimental design. The structure, thickness and growth angle of YSZ coatings on the entire circumferential surface of the cylindrical sample are studied. The results indicated that the structure, thickness and deflection growth angle of YSZ coatings are related to the orientation of deposition location. The numerical simulation of the multiphase mixed fluid near the substrate is carried out and the deposition regularity and mechanism of YSZ coatings prepared by PS-PVD is deduced. The growth rate is related to the local characteristics of the plasma flow field, and is directly proportional to the field pressure and inversely proportional to the field velocity. The growth angle of the coating is generally affected by the flow direction of the plasma jet. Especially, the normal component of velocity vector, V_norm, mainly affects the speed at which the coating grows vertically upwards. The tangential component of velocity vector, V_tan, determines the degree that the coating growth direction deviates from the vertical direction. When V_tan ≠ 0, the coating forms a fine column with a certain deflection angle and finally develops into an oblique columnar structure.     

Semiempirical Models of the Slow and Fast Solar Wind

Coronal holes can produce several types of solar wind with a variety of compositional properties, depending on the location and strength of the heating along their open magnetic field lines. High-speed wind is associated with (relatively) slowly diverging flux tubes rooted in the interiors of large holes with weak, uniform footpoint fields; heating is spread over a large radial distance, so that most of the energy is conducted outward and goes into accelerating the wind rather than increasing the mass flux. In the rapidly diverging open fields present at coronal hole boundaries and around active regions, the heating is concentrated at low heights and the temperature maximum is located near the coronal base, resulting in high oxygen freezing-in temperatures and low asymptotic wind speeds. Polar plumes have a strong additional source of heating at their bases, which generates a large downward conductive flux, raising the densities and enhancing the radiative losses. The relative constancy of the solar wind mass flux at Earth reflects the tendency for the heating rate in coronal holes to increase monotonically with the footpoint field strength, with very high mass fluxes at the Sun offsetting the enormous flux-tube expansion in active region holes. Although coronal holes are its main source, slow wind is also released continually from helmet streamer loops by reconnection processes, giving rise to plasma blobs (small flux ropes) and the heliospheric plasma sheet.     

基于高温热管的超燃燃烧室热防护结构

提出了基于先进热管理思想的燃烧室热防护结构.面板采用腔体式平板高温热管，实现面板等温化，降低局部高温区的温度；在热管腔体内部设计燃油冷却通道，实现对超燃燃烧室面板的燃油主动冷却.对其各项性能进行了数值分析，给出了设计参数对系统性能的影响规律，并完成了结构样件研制及石英灯试验考核.典型设计状态下，其单位面积质量为无氧铜面板的35.4%，高温合金面板的38.2%.石英灯局部加热条件下，面板最高温度为1123K时最大温差为80K.相比于传统燃油冷却方式，该型防热结构能够有效提高超燃发动机燃烧室热防护的整体性能，是超燃发动机热防护的一种重要概念.      

新型无阻流板反推装置流场结构和影响参数研究

应用CFD数值模拟技术,分析了基于二次流喷射的无阻流板涡扇发动机反推控制机理,并详细分析了二次流喷射的压力、位置、角度及喷流缝宽度对流场结构和反推性能的影响。在计算模型中没有考虑反推导叶,为此,用与反推力成正比的参数,即反推质量流量比来衡量性能的优劣。计算结果表明:二次流喷射压力、位置和角度是影响反推性能的重要参数,并且在一定的风扇涵道流下,存在最佳的二次流喷射位置、角度和压力;而二次流喷射缝宽度在一定范围内对反推性能影响不大。     

同轴双剪切气-气喷嘴数值模拟

通过求解k-ε湍流模型的Navier-Stokes方程组,对以气氢/气氧为推进剂的同轴双剪切喷嘴燃烧室流场进行数值模拟研究.结果表明:气-气推进剂在燃烧室内形成两个燃烧面,能在大流量下实现较高的燃烧效率;同轴双剪切喷嘴的氧喷注速度和氢氧速度比是同轴双剪切喷嘴设计的关键参数,氧喷注速度越小使推进剂的燃烧位置越远离喷注面,而氢氧速度比越大使燃烧位置越靠前;与气-液同轴剪切喷嘴相比,同轴双剪切气-气喷嘴的设计可以取较大的氧喷注速度和适合的氢氧速度比.      

压气机转子出口流场的发展及三维紊流特性

用单斜热丝和高频压力探针仔细测量了单级压气机转子出口及下游的三维紊流流场, 揭示转子出口流场的变化, 分析不同流量状态下单级压气机转子出口的三维紊流特性。      

低雷诺数下吹吸气射流对NACA0012翼型气动性能的影响

基于Favre过滤的大涡模拟方法,对雷诺数Re=10^4,迎角α=6°下的NACA0012翼型上表面吹吸气射流进行了数值模拟,从翼型周围流场流线图、速度场云图、上下表面压力系数曲线以及上表面边界层位移厚度等多角度地分析了射流位置以及速度变化对翼型气动性能的影响。结果表明:射流位置对翼型气动性能影响较大,且吸气射流要明显优于吹气射流。对于吸气射流,前缘吸气要明显优于中后缘吸气,可有效增升减阻,并减小翼型尾部流动分离,抑制翼型气动参数扰动,其最佳吸气位置随着速度的增大逐渐向下游移动;而吹气射流对翼型气动系数的作用效果较差,但中后缘的吹气射流可减小飞行过程中的气动扰动量,且吹气越大,效果越明显。