高压涡轮变叶顶蜂窝角度二次流动特性研究

为减少叶顶泄漏流带来的气动损失，本文对高压涡轮叶顶复合蜂窝的排布角度进行寻优，并分析其气动性能。研究过程保持叶顶蜂窝几何形状不变，改变复合蜂窝在叶顶的排布角度，降低叶顶二次流的总压损失系数和叶顶相对泄漏比。以叶栅出口下游30%轴向弦长位置的面平均总压损失系数为目标参数，利用Isight软件嵌套图形-网格自动生成流程，对0～57°旋转角度内的蜂窝排布方式进行遍历寻优，得到低总压损失的蜂窝排布方式。研究表明，最优排布结果与平叶顶相比，叶栅总压损失降低5.21%，与基准角度蜂窝相比降低1.34%。最优排布方案对叶顶泄漏流的阻碍效果更明显，增大了蜂窝对气流的耗散能力，降低了跨叶顶的横向驱动力，减少了泄漏涡的损失。     

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

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

复杂内筋铝筒段旋压变形规律和再结晶组织演变数值仿真

流动旋压是一种先进的内筋薄壁筒段整体近净成形方法。本文运用数值仿真方法研究多级筋铝合金筒段热旋压过程中材料流变规律和组织演变特性。将内变量型2219铝合金本构模型嵌入ABAQUS仿真平台建立了准确反映宏观流变和微观晶粒演变的热旋压仿真模型，完成了具有多级特征的内筋筒段旋压成形规律分析。研究表明在减薄率为50%、温度为300℃的旋压参数下，随内筋宽度增加材料由“挤压”式转变为“塌陷”式填充内筋，筋槽内材料变形量下降，筋槽内材料晶粒大小不均匀性增加；窄筋在250℃下的晶粒尺寸明显小于其在350℃下的，而宽筋在250～350℃下的晶粒尺寸较为接近，成形温度对窄筋晶粒大小有明显影响，而对较宽内筋的影响有限。     

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

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

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架，设计开发了LES(Large Eddy Simulation)湍流模型，主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法，给出设计的软件框架和计算流程，尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例，验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型，具备高精度数值格式、亚格子模型、湍流统计等通用模块，可为国内学者提供一个LES湍流模拟研究的开放平台。     

同轴旋转倒置圆台环隙间流动特性分析

针对同轴旋转倒置圆台环隙间流体复杂流动问题, 对其环隙间流动特性进行了实验研究。重点进行染色液流动显示实验和PIV流场测速实验, 对实验结果做定性及定量分析, 研究内筒转速和环隙宽度对环隙间流动特性的影响。染色液流动显示实验和PIV流场测速实验分别定性和定量地展现了环隙间螺旋涡的产生及变化过程。对不同内筒转速和环隙宽度下的螺旋涡涡心运动周期进行分析, 结果表明, 内筒转速升高, 周期减小；环隙宽度增大, 周期增大。运用瞬时流动和时均流场解析了环隙间螺旋涡运动产生机制, 探究内筒转速和环隙宽度对3种雷诺应力大小的影响与分布情况。内筒转速变化, 雷诺径向正应力始终最大；环隙宽度变化, 雷诺切应力始终最小。      

基于新型丝状电极等离子体激励器的前体涡控制

本文首次将新型丝状暴露电极DBD等离子激励器应用于大迎角下细长体非对称涡控制。丝状暴露电极的材料的选择对DBD推力以及推力效率至关重要,通过地面精细推力测量对丝状暴露电极等离子体激励器进行了优化,结果表明,本文研究材料中采用钨丝作为暴露电极,其推力效率最优;且随着电极直径从d=0.3 mm减小到d=0.08 mm,DBD推力效率显著提升。基于优化后的DBD激励器,将其应用于前体非对称涡控制:未施加等离子体控制时,压力测量以及PIV结果均表明细长体背风区流场为明显的非对称涡结构;在等离子体激励下,该非对称涡结构可变为对称甚至反向非对称,且非稳态激励控制能力明显优于稳态激励。研究发现,大迎角下细长体非对称涡控制与背风区原始涡系结构有关,其中包含对称涡系和非对称涡系。本文研究为大迎角下细长体非对称涡控制提供了一种新思路,同时也为丝状暴露电极DBD等离子体激励器的应用提供参考。     

电磁力作用下槽道流速度响应的放大机理

在电解质溶液中,电磁场产生的电磁力可以控制流体的运动,从而达到很好的减阻、增升、减振等效果。但由于所施加的电磁力较大,导致控制效率很低,因此以较小的电磁力诱导出大的流动响应成为提高流动控制效率的关键。以层流槽道流动作为研究对象,在槽道的下壁面施加沿展向余弦分布的展向电磁力,推导了线性条件下流向响应速度的解析解,并通过直接数值模拟对非线性条件下的响应进行了计算。结合解析解和数值解,揭示了流场中速度响应的放大机制,讨论了电磁力和流场参数对响应放大效果的影响。结果表明:当振幅较小时,速度响应处于线性范围内,其放大倍数与Re~2成正比,随着渗透深度的增大,先迅速增大后缓慢减小;随着展向波数Kz的增大单调减小。随着振幅的增大,放大倍数进入非线性范围,其值逐渐减小,但速度响应值先增大后减小。在振幅处于10~(-3)～10~(-2)量级时,速度响应可达到的最大值超过0.2,此时的放大倍数在102量级。因此,利用流场的放大效应,是实现高效流动控制的重要环节。     

高压捕获翼构型跨声速流动特性初步研究

高压捕获翼是一种可以在高超声速条件下同时获得高升阻比、高容积率和高升力系数的新型布局概念。为初步分析该类新型布局的宽速域气动特性,以一种圆锥-圆台体组合高压捕获翼原理性构型为计算模型,以典型跨声速条件(马赫数0.92、0°攻角)为计算工况,进行了计算和分析。流场分析结果表明,在跨声速流动条件下,机体与高压捕获翼之间存在比较强烈的气动干扰,且干扰的强烈程度与高压捕获翼-机体间的垂向距离大小直接相关。与不带捕获翼的参考构型相比,增加捕获翼会导致机体尾部分离区范围增大,并在机体与捕获翼之间的开放通道内形成类似于变截面收缩管的流动,致使沿流向方向出现了明显的激波串,进而导致捕获翼下表面壁面压力出现较为明显的波动。同时,由于机体和捕获翼间的垂直距离沿展向方向逐渐增加,导致该波动在对称面附近最为剧烈,然后随展向位置逐渐增加,压力波动逐渐减弱。     

运载火箭推进剂复杂流动传热问题数值模拟中的模型简化方法

在运载火箭改进优化或新型运载火箭设计过程中,会遇到各种复杂的流动与传热问题。按照真实产品和物理过程直接建立的数值计算模型过于复杂,不同尺度的结构、流动与传热、长时间的飞行过程相互耦合,出现计算时间太长(数年)、收敛困难、程序调试困难等问题。针对初始条件复杂、飞行过程复杂、边界条件复杂3类具体问题,总结了3种相应的模型简化方法,分别是工程经验法、极限参数法、低维近似法,为推进剂复杂流动与传热的数值模拟提供参考。