基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

左行运动激波主导的管内流动特性研究

为了分析左行运动激波主导的管内流动特征，本文采用非定常数值仿真方法，对亚声速进口条件下等直管道内左行运动激波传播与演化特性、左行运动激波/边界层干扰特征开展研究。研究结果表明：在出口周期性强压力脉动干扰下管内存在连续的左行运动激波，该左行运动激波传播特征具有相似性，激波强度、传播速度按幂函数规律衰减。气流经过左行运动激波后总压、总温、静压阶跃式升高，随后受膨胀波影响气流总压、总温、静压下降；左行运动激波/边界层干扰诱发形成翼型回流区，该回流区随运动激波强度衰减逐渐减小。理论与数值分析表明存在左行运动激波后速度为零和运动激波两侧总压相等的两个临界状态。波前马赫数低于临界值或左行运动激波强度高于临界值时，左行运动激波后为倒流、波后总压高于波前。     

翼身融合运输机分布式电推进系统设计及油耗评估

针对翼身融合运输机开展了分布式电推进系统的总体设计与油耗评估。通过数值计算完成了70t载质量翼身融合飞机的气动设计与优化。在巡航马赫数为0.80和10km高度的设计点,最大升阻比达到了24。通过求解积分边界层方程组,完成了电推进系统的总体设计。电推进系统包含10个推进风扇,风扇直径为1.45m，压比为1.35，巡航功率为2.94MW。建立了考虑燃烧过程的发动机一维性能模型,对发动机油耗进行了评估,获得了不同发动机循环参数下燃油消耗。建模结果表明,基于翼身融合布局和分布式电推进技术,可使运输机的油耗较C-17节省近50%。      

基于等离子体激励的湍流边界层减阻控制

为了探究介质阻挡放电（dielectric barrier discharge, DBD）等离子体气动激励对平板湍流边界层的减阻情况，在控制来流速度为10.7 m/s的低速风洞中进行等离子体平板湍流边界层减阻控制实验。本实验重点研究了不同激励频率对湍流边界层的减阻控制效果，使用热线风速仪系统采集流向速度信号，获得边界层平均速度分布和脉动速度分布。对实验结果进行对比分析发现，在施加不同频率的等离子体激励之后，边界层内对数区速度明显减小；随着激励频率的增加，局部减阻率呈现出先增大后减小的趋势，在激励频率为200 Hz时，减阻率达到最大为7.4%。     

壁面温度分布的高超声速平板边界层稳定性分析及转捩预测

针对来流马赫数为4.5、6.0和7.0的高超声速平板边界层，取30km高空处的气体参数，壁面为等温、绝热和温度分布等3种不同条件，采用eN方法进行转捩预测。其中，壁面温度分布条件下，在等温壁（冷壁）和绝热壁之间，给出4种流向分布进行分析。取扰动的初始幅值为0.3‰，以幅值达到1.5%作为转捩判断依据。结果表明:当温度为来流温度时，等温壁面条件的转捩位置最靠前，并随马赫数增大更加靠前；绝热壁面条件的转捩位置随马赫数增大而推后；壁面温度分布条件下，在相同时刻，马赫数越大，转捩位置越靠前。相同马赫数下，壁面温度较高者转捩位置较靠后（马赫数为7.0时，不完全满足此规律）。在马赫数为4.5和6.0时，绝热壁面条件转捩由第一模态主导，其余情况主导转捩的都是第二模态。     

温度对高速平板边界层转捩雷诺数的影响

为了深入研究温度对高速边界层稳定性和转捩的影响,并对不同温度条件下的转捩结果进行关联,采用线性稳定性理论和eN方法,对不同壁温、总温条件下的高速平板边界层进行了稳定性分析和转捩预测。温度条件对转捩的影响体现在两方面,一方面会改变快、慢模态的分枝特性,由此改变扰动的线性增长路径;另一方面可以改变边界层本身的不稳定性,引起扰动增长幅值和转捩位置的变化。结果表明,壁温或总温的变化,会产生不同的模态分枝类型,且总温的增加对边界层起稳定作用。另外,发现对于相同马赫数,在不高于1000K的壁温条件下,扰动的增长与壁温存在一致性的变化规律,结合eN方法,获得了转捩雷诺数与壁温比和N值的函数关系式。利用该公式,可以对不同初始扰动幅值和壁温条件的边界层转捩进行预测和关联。     

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

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

带轴对称台阶的圆锥高超声速边界层转捩试验

在Ma 6风洞内，通过高频脉动压力测试技术和基于纳米粒子示踪的平面激光散射(NPLS)技术，分别对带前向、后向轴对称台阶的圆锥高超声速边界层转捩进行了试验研究。采用功率谱密度分析和互相关计算等方法对脉动压力数据进行分析，得到了边界层中扰动波的发展规律，定量分析了第二模态波的相关参数。结果显示：两种模型中第二模态波在沿流向向下游发展的过程中，其幅值均先增大再衰减、特征频率均逐渐减小；特征频率和传播速度整体上均随雷诺数的增大而增加(后台阶模型中特征频率由100 kHz增至196 kHz，前台阶中则由 97 kHz 增至174 kHz)、波长变化规律则与之相反(后台阶中由6.35 mm降至4.54 mm，前台阶由7.35 mm降至 4.66 mm )；后台阶模型中第二模态波初次出现位置比前台阶中更靠近上游，边界层转捩位置较前台阶前移。将NPLS结果与高频脉动压力测试结果进行对比，两者吻合较好。     

基于MSFLE与Liutex的涡合并实验研究

应用具有拉格朗日性质的运动单帧长曝光(MSFLE)图像测量方法,对矩形管内湍流边界层涡结构合并现象开展了实验研究。研究了雷诺数Reθ在97到194之间管道湍流边界层流向-法向平面内涡合并现象。在实验中,测量系统以与涡运动速度相近的速度匀速移动,采用连续长曝光记录示踪粒子运动轨迹,以捕捉涡随时空的演变过程,骨架提取图像处理算法获得示踪粒子速度,并应用Liutex物理量表征涡的旋转强度。研究表明:MSFLE方法测量装置简单,对实验条件要求低,可从拉格朗日视角直观测量湍流边界层涡结构及周围流场的时空演变过程,MSFLE图像测量方法与Liutex涡识别算法结合可以很好地应用于湍流边界层涡结构的可视化与量化。矩形管内湍流边界层涡合并的条件是两个涡相邻、强度和尺寸基本相同且为同向旋转涡,合并中两个涡的强度呈反向变化,合并生成的新涡的强度和尺寸基本为初始合并时两个涡之和,且旋转方向与两个涡同向。