超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

0.3 m低温风洞液氮供给系统研制

基于系统级的一维热流体模拟分析,建立了适用于研究分析0.3 m低温风洞液氮供给系统的数学模型,并开展了系统漏热、两相流及缓冲罐中液氮容积等流体动力学分析;在系统现有控制策略及试验数据的基础上,基于该数学模型开展了系统压力动态响应分析,获得了在阀门动态调节过程中管网压力的瞬态响应,计算结果与试验值的总体误差控制在10%以内。喷射压力一致化改造避免了阀间干扰,添加的回流管道消除了供给末端的两相流现象,使喷射压力控制精度达到1.1%,调节时间减少到23 s,实现了风洞总温快速安全调节和精确控制。     

转子系统挤压油膜阻尼器设计与动力特性

鉴于挤压油膜阻尼器(SFD)设计必须要同时考虑转子系统的动力学特性,提出了一种转子系统与挤压油膜阻尼器耦合设计方法,给出了详细的设计流程。对所设计的阻尼器进行了CFD数值模拟、油膜压力测量以及减振效果实验,结果证明所提出的设计方法是有效的。相比未采用阻尼器,采用阻尼器后转子系统两个转盘的振幅分别下降了46%和39%。通过实验还研究了不平衡量、支承刚度、供油压力和滑油温度对挤压油膜阻尼器减振效果的影响。结果表明,相比于不平衡量和支承刚度对减振效果的影响,供油压力和滑油温度的影响并不显著。进行挤压油膜阻尼器设计时,重点应该关注转子上的不平衡量大小和支承刚度。     

退化量有缺失的无失效小样本轴承可靠性评估

针对有试验性能退化量参数记录、且(或)有部分数据缺失情况的小样本无失效轴承试验问题,通过由Taylor和Thompson提出的数据模拟法实现补全样本退化量,结合Bootstrap自助法扩大样本量,再根据基于性能退化轨迹的补充信息方法来进行其可靠性评估。选取7组受试轴承的振动退化量,对比在完整数据和带有缺失数据情况下的分析结果,发现可合理利用原舍弃不合规试验的部分有效信息,使之增加可靠性评估的样本数,从而得到较不用这些数据更为准确的结果,且所得结果较用完整数据结果绝对值相差在01以内。对比由极大似然估计法和加权E-Bayes法分析试验寿命数据的结果,发现该方法所得评估结果更优,与试验实际相比误差在10%以内,对于提高评估精度及降低试验成本有积极的实际意义。     

空心气冷低压涡轮动叶气热耦合数值模拟

以空心气冷低压涡轮动叶为研究对象,采用高质量的流体域和固体域网格控制技术,带转捩模型的双方程SST湍流模型,开展了基于CFD方法的叶片气热耦合问题研究。获得了不同冷气流量比(分别为1.0%,1.38%,1.8%和2.2%)、温比(分别为2.1,2.25,2.3,2.4和2.5)和压比(分别为1.4,1.6和1.8)对叶片换热特性的影响规律,设计状态中截面按弧长平均的叶片壁面金属温度计算值较试验值偏小0.3%,气热耦合计算的叶片壁面温度分布与试验结果吻合良好,验证了气热强耦合计算方法的精度,为涡轮叶片温度场分析提供了一种有效的方法。     

舰载直升机理论风限图飞行模拟获取方法

获取舰载直升机实用风限图,仅依靠实装飞行试验或建模仿真都是非常困难的。通过起降飞行模拟试验获取理论风限图,结合少量的海上实装飞行验证试验是一种综合较优的方法。起降飞行模拟试验需要真实飞行员和高等级飞行模拟器参与,飞行模拟器中的母舰气流场数据、母舰运动可由建模仿真或试验测试得到,且必须考虑母舰气流场的非定场特性;只有当飞行员视景或视觉反馈受限时(如能见度降低时),飞行模拟器基座的六自由度运动才是必需的;用DIPES或Bedford方法量化飞行员主观操纵负担,结合记录的客观操纵数据,评判起降难度;典型的飞行模拟试验程序包含了不同风况下的进近和甲板着舰,通常每隔15°和5kn一个风况。通过起降飞行模拟试验得到的理论风限图,可以更加有效和安全地指导海上实装飞行验证试验,最终获得实用风限图。     

基于并行CFD和优化技术的返回舱外形多目标优化设计北大核心CSCD

在给定的质心设计范围内,围绕球冠倒锥返回舱外形的高超声速气动单点静稳定性、配平升阻特性、质心横偏量的综合设计问题,提出了多点多目标优化设计数学模型。通过多目标优化设计方法结合并行数值模拟技术,对该多点多目标气动外形优化设计问题进行研究,为了加快多点数值计算进度,采用了嵌套并行方法,通过有效利用硬件资源来提高多个状态气动数值求解效率。根据以上方法给出的最优设计边界指出了返回舱单点静稳定性与配平升阻比和质心横偏量的矛盾关系,改善单点静稳定性会导致配平升阻比下降,使质心横偏量增加;反之,配平升阻比增加,质心横偏量减少都会使单点静稳定性变差。     

突加不平衡下GTF发动机转子系统瞬态响应分析

针对GTF发动机低压转子系统突加不平衡瞬态响应问题,基于有限元法建立了考虑啮合单元、行星架单元、膜盘联轴器单元的多体接触、多转子耦合动力学模型,介绍了考虑非线性因素下突加不平衡瞬态响应计算方法,求解了低压转子系统的突加不平衡瞬态响应,并分析了关键部件刚度对突加不平衡的影响规律。结果表明：当风扇转子发生突加不平衡时,各部件均表现为振动幅值突增,后迅速趋于稳定,低压涡轮转子发生明显拍振,突加不平衡载荷主要由行星架结构和第1.5支点刚性支承承担;行星架刚度主要影响风扇转子突加不平衡响应、第1.5支点和行星架处的外传力;膜盘联轴器刚度主要影响增压级突加不平衡响应、转子在突加不平衡之后稳定运转情况。