涡脱落撞击潜入式喷管引起低频压力振荡研究

通过不同障碍高和不同加质的冷流数值模拟,着重研究碰撞频率（碰撞处压力振荡频率）与强度(碰撞处p_rms/p_mean)对于SRM内压力振荡幅频特性的作用程度。研究结果表明,SRM内压力振荡幅度与碰撞处的压力振荡大小为相同量级。当碰撞强度较小时,碰撞只起到增强声场的作用;当碰撞强度较大时,碰撞产生的压力振荡随碰撞强度的增加逐渐覆盖淹没声场。SRM内压力振荡的幅频特性主要由涡脱落碰撞强度决定,一般只呈现低频且具有非声非线性特征。     

DES与DDES在湍流分离中的原理与性能研究

随着工程上流动结构的日益复杂，兼具雷诺平均Navier-Stokes(RANS)方法高效率与大涡模拟（LES）高精度的分离涡模拟（DES）类混合方法成为现阶段工程中最有效的湍流模拟方法之一。围绕DES类混合方法中的DES与延迟分离涡模拟（DDES）方法开展工作，分析二者开关函数构造上的不同，研究延迟因子作用机理，并考察DES与DDES方法的求解能力。研究表明：DES与DDES方法在模拟表现上存在一定差异，DDES方法通过引入延迟因子，保护RANS求解区域，改善模化应力不足，降低了DDES方法对交界面系数CDES敏感程度；DDES方法在计算过程中容易出现过度保护，导致求解瞬时涡结构能力不如DES方法，分析与延迟因子引入比重及开关函数构造形式有关。     

基于CFD/CSD耦合的高速射弹尾拍载荷特性研究

为提高射弹尾拍载荷的预测精度，建立了一套基于计算流体动力学（CFD）/计算结构动力学（CSD）双向耦合分析的计算方法和程序。射弹流体计算主控方程采用耦合SSTk-ω湍流模型和Schnerr-Sauer空化模型的Navier-Stokes方程，射弹结构计算采用基于模态叠加法简化的结构动力学方程，流固耦合界面插值采用径向基函数法，网格变形采用弹簧网格法。分别对泡型计算方法和流固耦合方法进行验证，在此基础上，计算对比1 000 m/s速度下射弹刚体和弹性体的尾拍泡型、结构变形和尾拍流体载荷特性差异。计算表明：弹性体尾拍过程，射弹泡型会产生弹身“二次拍击”、“局部沾湿”和沾湿面积增大等特殊现象，结构变形由弹性一弯模态主导，较大的变形引起流体载荷增大27%～105%，尾拍姿态角增大13%，尾拍频率增加20%，流固耦合效应对尾拍泡型、尾拍载荷和尾拍弹道均产生了较强的影响。     

气液针栓式喷注器液膜破碎过程和结构的数值研究

为了更好地认识针栓式喷注器雾化场的结构，基于网格自适应加密技术以及VOF (volume of fraction)方法追踪气液的分界面，采用realizable k-ε湍流模型模拟整个流动过程，还原了不同时刻气/液撞击的初次破碎过程，数值模拟结果与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。进一步对针栓式喷注器气/液撞击的初次破碎过程、内部流场涡结构、速度场进行分析，研究了初次破碎雾化的动力学过程和机理。研究结果表明：液桥的形成主要是由液洞的扩展和拉伸、合并而形成，而液滴主要是由中心液膜拉伸、液丝断裂以及液桥断裂而形成，液膜破碎阶段形成的涡结构是造成液膜断裂的主要原因。     

基于变权重-正态云模型的飞机轮胎滑水风险

针对飞机轮胎滑水行为的随机性和模糊性特征，提出基于正态云模型的滑水安全评价分析方法。引入变权理论动态调整权值，采用惩罚性变权函数降低常权权值对评价结果的主观影响。构建飞机轮胎滑水流固耦合仿真分析模型，选取飞机轮载、滑行速度、积水厚度、道面摩擦系数及刻槽深度作为风险影响因素，基于单因素云模型数字特征及变权向量求解综合隶属度，建立多元决策下的飞机滑水风险等级及划分标准。以某山区多雨机场轮胎滑水事件为实例进行验证，结果表明：以传统临界滑水速度指标进行条件判定仅得出允许起降的一般结论；对比常权和变权评价结果，工况1安全系数由1.09提高至1.17，工况2由2.09提高至2.94，可定量描述道面起降环境差异，滑水风险仍在可接受范围内，变权评价结果偏于保守；工况3安全系数由3.13提高至3.74，滑水风险上升至Ⅳ级，即使道面积水厚度符合上限要求，轮胎滑水发生几率仍有可能显著提高，与实际风险情况基本一致，对道面运行安全分级管理具备参照性。     

带自由表面水流的分离涡模拟

目前水动力学问题的数值求解仍然广泛采用RANS模型,特别对于大尺度的地表水流的数值模拟,该类模型计算效率较高,但其难以给出湍流场的小尺度涡结构.深入研究复杂地形条件下的水流运动,需要发展诸如LES等的高级数值模型.天然地表水流,地形、边界复杂,雷诺数通常较高,高级数值模拟受限于计算能力的限制,还难以工程应用.RANS和LES混合模型是目前具有工程应用的一种混合数值模型,其已经在CFD领域获得了长足发展.然而对于地表水流运动的数值模拟,相关工作还较少.本文将分离涡模型(DES),即一种RANS和LES的混合模型,应用于带自由表面的地表水流运动,建立了一套数值仿真模型.模型基于有限体积法,水平面内采用非结构计算网格,垂向为结构化网格,对流项离散格式采用二阶TVD格式,并行基于OpenMP语言库.算例表明DES模型有助于揭示复杂地形条件下带自由表面水流的大涡拟序结构.     

叶片反扭对跨声速大涵道比风扇性能的影响

使用基于流固耦合算法的叶片反扭程序,考虑了非定常气动力对叶片变形的非线性作用,研究了叶片反扭对跨声速大涵道比风扇性能的影响.以冷态叶型为起点,先计算离心力作用下的叶片变形,在此基础上使用流固耦合程序获得非定常气动力作用下的变形.考察了3个转速下叶片的动态变形对大涵道比风扇气动性能的影响.结果表明:在跨声速工况下,叶片表面激波位置的变化对叶片反扭角有很大作用,在考察的转速范围内,堵塞点使用设计叶型计算的流量大于动态叶型下的流量,数值达7%,将导致发动机起飞推力小于预测值.结果表明在大涵道比风扇设计阶段,预测气动性能使用准确叶型的重要性.      

细观周期性结构复合材料热固耦合双尺度渐进均匀化分析方法及有限单元法实现

针对周期性结构复合材料,根据稳态热固耦合平衡方程,将弹性矩阵和导热系数矩阵进行合并,将位移和温度作为空间场量同时进行小参数渐进展开,利用摄动理论和均匀化理论,推导并建立了周期性结构复合材料稳态热固耦合双尺度渐进均匀化分析方法。将复合材料力学分析从单一物理场推广到热固耦合场,为实现复合材料构件的热固耦合多尺度力学分析建立了基础。根据建立的偏微分方程式,利用变分原理推导得到其有限单元形式,用VC++6.0开发了热固耦合双尺度渐进均匀化分析程序。为验证该方法预测复合材料宏观等效性能的准确性,对纤维体积含量从20%到80%的单向纤维增强复合材料宏观等效参数进行了计算。当纤维体积含量低于65%时,得到的计算值与实验值吻合良好。更加真实地还原复合材料细观结构模型,计算精度将进一步提高。     

润滑系统管路流动与换热耦合计算研究

为研究沿程管路换热与滑油流动之间的耦合问题,通过分析主要元件的流动和换热特性,引入元件换热的能量方程,结合管路流动的连续性方程和动量方程,建立了润滑系统管路热流耦合计算模型。针对一段发动机供油管路,利用MATLAB/Simulink建模仿真工具包,建立其热流耦合网络计算模型,获得滑油流量、温度以及管壁温度等参数,采用商业软件Flowmaster对其进行验证。模型计算结果与Flowmaster计算结果非常接近,滑油流量和温度的最大误差均不超过3%,管道壁温的计算误差不超过2.4%,说明建立的热流耦合计算模型具有良好的仿真精度,适用于考虑沿程管路换热的发动机润滑系统计算分析。     

翅片位置对复合式减涡器减阻性能影响数值模拟

为了探索翅片-管复合式减涡器的翅片安装位置对共转盘腔径向内流压力损失的影响规律,对不同转速、翅片周向位置及安装角度下的去旋系统开展了数值研究,得到了不同工况下共转盘腔径向内流的流场结构及压力损失分布曲线。研究结果表明:减涡管能引导流体径向流入,并降低流体的旋流比;相比于管式减涡器,翅片-管复合式减涡器能明显降低盘腔内的总压损失;在不同旋转雷诺数下,翅片的周向安装位置α及安装角β均存在最佳值;在中、高旋转雷诺数下,最佳值分别为α=9°,β=30°,最佳结构下总压损失较基础模型低40%左右;改变翅片周向位置及安装角度可以明显改变气流进入减涡管的角度,在较优情况下,可以减小流体流入减涡管的阻力及在减涡管内的流动阻力,整体上减小了盘腔内总压损失。