基于多项式混沌法的翼型不确定性分析及梯度优化设计

耦合伴随方法和非嵌入式多项式混沌法，发展了高效、可靠的不确定性梯度优化设计方法。利用伴随方程法求解目标函数对不确定性变量的导数，发展了一种梯度增强型多项式混沌法。通过亚声速和跨声速下等多种算例可以证明该方法可以提高不确定性分析的效率和精度。同时，利用基于方差分解的全局敏感性分析方法对不确定性变量的敏感性进行了量化。建立了多项式混沌耦合伴随方程的统计矩梯度求解方法，并结合梯度增强型多项式混沌法搭建不确定性梯度优化设计系统。基于该优化设计系统对二维低亚声速和跨声速翼型开展确定性及不确定性优化设计研究。优化结果显示，相比于确定性优化设计，不确定性优化设计通过合理权衡确定性性能和不确定性性能，可提高抵抗马赫数和迎角不确定性扰动的能力，同时优化性能均值和标准差。其中阻力系数均值最大可降低17%，阻力系数标准差最大可降低80%。而确定性优化设计可能导致性能鲁棒性的降低。     

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

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

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

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

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

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

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

火箭发动机超声速过膨胀射流气动噪声特性研究

采用大涡模拟LES方法计算了火箭发动机超声速过膨胀射流形态及近场声压分布,研究了入口温度与环境温度的比值(温度比)对声场的影响;将声源分解,基于Ffowcs Williams-Hawkings (FW-H)方程获取了不同位置噪声源的远场噪声,并根据声压级频谱和湍流形态分析了超声速射流噪声的产生机理。研究表明,超声速过膨胀射流气动噪声由湍流混合噪声和宽频激波噪声组成,近场噪声源以马赫波形式向大方位角辐射中高频噪声,下游大尺度湍流向低方位角范围辐射低频噪声,声压级峰值频率随观测角度增大而升高;随温度比升高,马赫波辐射角度增大,噪声指向性发生改变。该研究可为运载火箭发动机地面试车或火箭发射段声学环境设计提供参考。     

7°尖锥高超声速边界层转捩红外测量实验

为了推动高超声速边界层转捩研究的深入开展，给边界层转捩机理研究、物理模型验证、转捩数据库构建和转捩天地相关性的建立等提供基础风洞实验数据，在中国空气动力研究与发展中心的Φ1 m高超声速风洞开展了边界层转捩规律红外热图实验。针对半锥角7°尖锥模型，研究了不同单位雷诺数、迎角和马赫数对尖锥边界层转捩位置的影响规律。实验单位雷诺数（0.49~2.45）×107/m，迎角范围-10°~10°，马赫数5~7，模型头部半径0.05 mm。通过红外热图技术测量模型表面温度分布，获得了较为详细的转捩位置和转捩参数影响规律。实验结果表明:在马赫数5~7范围内，马赫数增大，尖锥转捩位置提前，分析认为是高马赫数时的雷诺数较大、自由流噪声水平较高引起；随着单位雷诺数的增大，边界层转捩位置前移，转捩雷诺数保持不变，约为3.0×106；小迎角时，随着迎角的增大，迎风面边界层转捩推迟，背风面边界层转捩前移，在10°大迎角时，迎风区中心线转捩前移，出现迎角"转捩逆转"现象，背风区出现了流动分离导致的低热流条带。     

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

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