不同迎角下脊形前体绕流数值模拟研究

脊形前体有较强的背风涡流场,不同的前体形状对前体涡流场和气动力有很大的影响。本文针对脊形前体飞行器大迎角湍流大分离流动计算的困难,采用IDDES混合湍流模型,以及与之匹配的非定常算法,研究了不同来流迎角下脊形前体的气动特性,以及背风涡非定常演化、破裂的细致流动结构。选取了不同脊形角,以及不同上、下高宽比的脊形前体进行计算。计算结果表明,在迎角较小时,随着迎角的增大,前体主涡会逐渐增强,在迎角较大时,前体主涡破裂;在相同迎角下,脊形角较小时,前体涡较强,涡升力也更大;对于相同脊形角的前体,当上半截面高宽比较小时,前体主涡强度较大,前体涡破裂临界迎角较小,即会提前破裂。     

非恒定增压式微尺度拉瓦尔喷管流动特性

针对微尺度喷流在航天器运动状态切换时出现的非恒定增压变化,采用直接模拟蒙特卡洛(DSMC)方法对阶跃式增压和线性式增压两种模式下的微尺度拉瓦尔喷管流场进行了模拟,并对其变化过程中的流动特性进行了对比分析。结果显示：阶跃式增压会导致流动特性出现较大幅值的峰谷式波动,而线性式增压下的流动特性则呈现出线性变化的特点;黏性力对微尺度喷流的非恒定增压变化产生了重要的黏滞作用,在喉部扩张段至出口的流场中尤为明显;在设定的条件下,阶跃式增压过程中喷流产生的总冲量较线性式增压高59.5%,质量流量高74.7%,单位工质提供的冲量低约8.6%,波动性也高于线性式模型,阶跃式增压适用于系统需要较大推力改变运动状态且推进剂充足的情况,而线性式增压在系统精确微调或需要推进剂产生更高效能时具有明显的优势。     

基于拉格朗日方法的通气云状空泡三维非定常脱落特性研究

采用均相流模型并结合FBM湍流模型,对绕轴对称回转体通气云状空泡流动特性进行了三维数值模拟,基于实验结果对数值方法进行验证,同时利用基于拉格朗日体系的有限时间李雅普诺夫指数(FTLE)、拉格朗日拟序结构(LCS)和粒子追踪方法分析了其三维非定常脱落特性。研究结果表明:纵截面上空泡覆盖区域的拉格朗日拟序结构整体呈椭球状分布,内部为不规律的复杂拟序结构;不同横截面上拟序结构分布存在很大差异。空泡内部的非对称流动结构和周向流动导致空泡呈现很强的三维运动特性。反向射流在周向上推进的不同步性,是造成空泡呈现不规则断裂和大尺度U型空泡团脱落的主要原因。     

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

叶片冷却的静叶时序对涡轮的性能影响

燃气轮机工作过程中上游叶片的尾迹会周期性地扫掠下游叶片,这使得下游叶栅通道内的流场呈现显著的非定常性。而时序效应作为一种典型的非定常流动现象,已受到国内外学者的广泛关注。本文研究双级气冷涡轮,利用全三维非定常数值仿真的方法分析了第1级与第2级静叶之间的时序位置对涡轮流量和效率的影响。计算结果表明处于不同时序位置下的涡轮第2级静叶叶身冷气流量的差异可达15.6%,涡轮41截面效率和热效率的波动幅值则分别达到了0.33%和0.26%。涡轮冷气流量的变化一定程度上影响了41截面效率,但时序位置对下游叶栅通道流动损失的影响才是造成涡轮效率差异的主要原因。     

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

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

钝体稳定射流火焰的火焰面模型数值研究

基于稳态火焰面(SLFM)和交互式非稳态欧拉颗粒火焰面(EPFM)模型对Sydney大学CH4/H2钝体稳定湍流扩散火焰行了数值研究,采用修正的雷诺应力模型(RSM),同时对两种不同规模的甲烷详细化学反应动力学机理进行研究,比较了燃烧模型和反应机理对湍流火焰结构、活性自由基以及氮氧化物预测精度的影响,与实验数据对比结果表明:两种反应机理得到的温度场和主要组分分布基本相同;SLFM模型能对速度场和标量场的分布进行较为准确预测,采用EPFM模型修正后,部分区域OH预测结果更加靠近实验结果;采用EPFM模型对SLFM模型耦合GRI-Mech 211的计算结果修正后,NO量级降低近2倍,预测精度明显改善,与实验结果实现较好的符合,验证了化学反应动力学机理以及非稳态效应对氮氧化物预测精度的影响。     

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

高速颤振模型设计中颤振主要模态的判断

颤振模型设计时难以实现完全的动力学相似,需要对颤振主要模态进行合理选择。采用数值分析方法,对颤振模型设计时主要模态的选取问题进行研究。通过各阶模态振型下气动刚度系数的比较、指定运动形式下广义非定常气动力的计算和不同模态截断下颤振结果的收敛特性分析,研究了颤振分析时不同模态运动之间的相互影响,对模态运动引起的气动力和颤振特性变化进行评估。以高超翼面模型为研究对象的数值算例结果表明,几种分析方法所判断的颤振主要模态基本一致。其中基于振型的广义气动刚度系数参数,避免了非定常气动力的计算,可作为颤振模型设计或颤振分析时主要模态选取的快速判断方法。     

轴向间距对涡轮时序效应影响的数值研究

时序效应可以改进多级叶轮机械的效率,为了探讨轴向间距对时序效应的影响,并应用到叶轮机械气动设计中,对1.5级静/动/静布局的涡轮叶片进行了3维非定常数值模拟。结果表明:保持轴向间距不改变,转子居中时涡轮效率整体最高,时均效率最大值为0.866,转子前移时次之,时均效率最大值为0.861,转子后移时效率最低,时均效率最大值仅为0.859,且转子相对位置发生改变并不影响效率-时序位置曲线;总轴向间距变大,时序效应变弱,效率提升幅度仅为0.7个百分点,效率-时序位置曲线随之改变。