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

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

跨音压气机确定应力场分析

采用全三维黏性定常、非定常数值模拟技术,对NASA 67第1级的固有非定常流动进行了研究.首先分析了定常流场和非定常时均流场间的差别以及压气机中典型的转、静干涉现象,然后对通道内确定应力的分布规律和转、静交界面确定应力各项的分布特点进行了较为深入的研究,揭示了现有确定应力模型的一些局限性,可为确定应力建模提供一定的指导作用.     

应用流固耦合数值方法研究机翼的射流减振技术

采用流固耦合的数值计算方法研究了NACA 0 0 15翼型在大迎角 ( 15°～ 6 0°)范围的颤振 ,以及在翼型背部部分引入射流的减振技术。在流体区域用高精度、高分辨率算法求解Farve平均的Navier Stokes方程 ,在固体区域用四阶Runge Kutta法求解振动方程 ,并且每一个时间步后都在两个区域之间传递边界条件。计算结果表明在翼型背部引入适当射流会降低翼型的振动 ,并提高升力。但如果引入射流的速度过高 ,会在叶背处形成新的分离流 ,升力反而会下降     

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

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

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

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

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

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

整体弹性结构法及在叶片流固耦合分析中的应用

针对三维叶片时域流固耦合振动响应计算普遍耗时的问题,采用一种假设整体结构的模态位移,求解固体时域响应和实现高效网格变形,发展了一种3维时域流固耦合分析的整体弹性结构方法,并应用于压气机叶片0°和180°相位角的气动稳定性分析中。结果表明,所发展方法的计算结果与传统双向时域算法和文献的结果较为吻合,而计算效率相比于传统算法显著提升;在所分析的两个相位角下,叶片振动的气动阻尼均随流量减小先增大后降低,相比于0°相位角,180°下叶片的气动稳定性大幅提高,表明该方法能有效应用于叶片的工程流固耦合研究。     

气流激励下叶片动力响应分析方法

利用在物理过程上的耦合代替数学方程上的耦合,完善了双向顺序耦合求解的理论,进一步给出了双向顺序耦合的数值求解方法,并描述了其在物理上对应的过渡过程.结合具体工程算例,分析得到某小型发动机二级静叶在气流激励下的流场分布形式和叶片动力响应的特点,指出了双向顺序耦合数值求解方法的优越性.      

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

采用流固耦合的实体元空心叶片鸟撞数值模拟

针对目前通用流固耦合算法在模拟实体元结构破坏上存在不足,以MSC.Dytran软件为平台,研究和验证了不考虑失效和考虑失效的实体元平板叶片流固耦合数值模拟方法;在此基础上,结合大涵道比航空发动机工作过程中较为常见且非常严重的鸟撞事故,建立了实体元空心叶片鸟撞瞬态动力学有限元模型并进行相应计算,结果表明:计算较好地模拟了叶片在遭受鸟体撞击后会产生巨大的瞬时冲击应力,以及叶片由此产生局部塑性变形.最后模拟了叶片遭受鸟撞发生失效的过程.