可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

转盘轴承集聚效应分析及其对回转性能的影响

针对航天发射台转盘轴承工作时出现滚动体集聚导致摩擦力矩增大的问题，研究转盘轴承中滚动体之间间隙变化对集聚效应的影响。通过对典型工位处进行有限元分析，对比不同接触类型转盘轴承的滚道变形和滚动体接触作用力的分布规律。在确定滚动体与隔离块的几何关系后，进一步分析了滚道变形和滚动体接触作用力的时变性对滚动体间隙变化的影响，并将间隙变化和集聚效应建立了联系。有限元计算得到的发射台框架变形与实验测量数据一致，验证了有限元计算的可靠性。仿真结果显示，由于下滚道发生轴向变形，滚动体在运动过程中不断交替进行“爬坡”和“下坡”；滚道径向变形导致滚动体之间增大了较多的间隙；各滚动体的驱动摩擦力不断变化，运动状态具有时变性，当转盘轴承中间隙增加时，运动状态不断变化的滚动体导致增加的间隙出现累积，引起滚动体集聚。为消除滚动体集聚效应，探索其对摩擦力矩的影响，进行了对比实验。结果表明:集聚效应导致发射台的摩擦力矩增加，且变化剧烈。      

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

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

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

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

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

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

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

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

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

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

微流控芯片技术在血细胞变形和流动性分析研究中的应用进展

近年来，随着微流控芯片技术的快速发展，微流控芯片在生物医学研究领域得到了广泛关注。由于其具有高通量、高灵敏度、集成化、低消耗及可控化等诸多特点，为在多细胞水平研究细胞迁移和分选动力学提供了新的技术平台。利用微流控芯片微通道结构设计灵活的特点，可在实验条件下模拟正常的生理和病理条件下的复杂血管；其微米尺寸的微通道也适于单细胞引入、操纵及检测。因此，用微流控芯片技术在单细胞层面对细胞生物力学性能表征也引起了广泛关注。以健康和疾病中的血细胞为例，从单细胞变形、流动、黏附、机械疲劳等力学性能表征到多细胞迁移及分离动力学等方面归纳目前微流控芯片技术在细胞力学分析和表征方面的研究进展。     

基于滑动子空间的杂波样本筛选方法

针对受目标污染的非均匀杂波辅助样本筛选问题，本文利用超低旁瓣子空间平滑加窗处理，克服已有的非均匀杂波样本筛选方法面临的初始协方差矩阵包含目标信号而引起的目标相消问题。一方面，采用超低旁瓣子孔径滑窗方法大幅提升旁瓣区信杂噪比，结合拟合方差偏离量提取目标大信号能量并予以剔除。另一方面，结合广义内积（Generalised Inner Product,GIP）等算法在子空间平滑的基础上进一步剔除非均匀杂波样本，从而得到独立同分布的杂波协方差矩阵。基于子孔径滑窗的杂波样本剔除方法能够大幅降低系统运算复杂度，提高系统应用效能。通过对机载实测数据的验证分析，该算法能够改善目标的输出信杂噪比，保留微动效应的离散弱目标信号，具备很高的工程应用价值。     

短螺旋型燃烧室旋流结构及气动边界发展规律数值模拟研究

短螺旋型燃烧室的头部轴线沿发动机周向与发动机轴线偏转一定角度,可以在保证燃烧效率不变的条件下有效缩短燃烧室轴向长度。短螺旋型燃烧室流场的最大特点是旋流流动单侧受限。为了研究头部安装角α变化对燃烧室旋流流动特性的影响,基于数值方法对短螺旋型燃烧室进行了计算分析。结果表明:随着α变化,旋流器下游旋涡依次出现对称环状、马蹄状、环状结构;随着α增加,气动边界逐渐出现并抑制旋流的周向扩张,导致流场出现不同的旋涡结构;不同α下的切向角动量随轴向距离增加而衰减,但α为35°和45°时,气动边界在非受限侧出现并对旋流产生约束,角动量衰减变慢;当α为0°,15°,35°,45°时,燃油液滴依次集中分布在旋流器下游两旋涡边缘、侧壁面和头部端壁、非受限旋涡边缘。本文研究了不同α下的旋涡结构及气动边界沿轴向的演化过程,为短螺旋型燃烧室进一步的设计与优化提供基础。