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

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

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

超声速凹槽大涡数值模拟研究

在超声速燃烧室内常采用凹槽增强燃料与空气的混合,实现火焰稳定,本文采用大涡模拟方法对开式凹槽流场进行了数值模拟研究,计算采用Smagorinsky 亚网格模型模拟小尺度涡的作用.数值模拟结果给出了凹槽流动的自激振荡的发展过程,以及凹槽流动中拟序结构的演变过程,如涡的卷起、增长,涡与涡之间的合并、破碎等,计算得到的凹槽波动的峰值频率和理论结果一致.     

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

控制策略对无轴承开关磁阻电机

分析了无轴承开关磁阻电机(BSRM)的3种负载控制策略(方波控制策略、最小磁势控制策略和平均悬浮力控制策略),以及每种策略下由悬浮力方程和转矩方程求解超前角和绕组电流的方法。由于控制对象和求解过程所引入约束条件的不同,即便相同负载状况时,3种控制策略求取的超前角和绕组电流结果也不同。利用瞬态有限元法获得了电机的动态磁密分布,并研究了3种控制策略下电机分别表现出的磁场特性。通过双频法分离铁损,分别计算了3种控制策略在不同转速时的电机涡流损耗和磁滞损耗,给出了电机铁心各部分损耗随转速变化的关系,由此得到了3种控制策略对铁心损耗的影响。     

边界层厚度对腔体气动声学特性影响数值模拟

为了研究来流边界层厚度对开式腔体气动声学特性的影响,基于分离涡模拟方法,计算了来流马赫数为2.0条件下,不同来流边界层厚度与腔体深度比时,长深比为5.88的腔体流动特性,得到了该腔体声压级的频谱特性.计算结果表明:随着来流边界层厚度增加,形成的剪切层稳定性增强,失稳后上下摆动幅度减少,失稳生成的大尺度涡与超声速主流的相互作用减弱,使得大尺度涡发展到腔体后缘时所具有的平动动能和转动动能降低.大尺度涡撞击腔体后缘在腔体内形成的气动噪声的声压级降低,最大减小幅度达7.5dB.同时各阶模态的频率也发生偏移,偏移值在100Hz左右.基于新的假设重新推导了Rossiter公式,明确了经验常数的物理意义,并以此解释了频率偏移现象.      

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。     

引入压力梯度的针对角区分离流动的Spalart-Allmaras模型的改进

针对Spalart-Allmaras（S-A）模型在角区分离计算中的问题,将无量纲的压力梯度引入其涡黏性输运方程的生成项,得到了改进的S-A模型.通过对两套含角区分离的低速压气机叶栅进行验证计算发现:与实验结果相比,原始S-A模型所得的分离区偏大,分离区内壁面压力偏低;而改进模型得到了与实验一致的分离区尺寸以及吸力面、压力面压力系数分布等结果.针对S-A模型涡黏性生成项和耗散项的分析表明:引入的无量纲压力梯度有效的识别了角区分离,在分离区内改变了涡黏性的生成、耗散关系,增大了涡黏性,从而缩小了计算所得分离区,同时在主流区保留了原始S-A模型的计算结果,进而带来了良好的改进效果.      

乙醇-空气稀液雾值班火焰的大涡模拟研究

为加深对气液两相湍流燃烧现象的认知,检测火焰面模型在液雾燃烧中的适用性,本文在欧拉-拉格朗日架构下使用火焰面/反应进度变量模型(FPV)数值模拟了乙醇-空气稀液雾值班火焰。欧拉坐标系下的气相湍流场使用大涡模拟方法模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑了相间质量、动量、能量交换。燃烧模型中采用碳元素定义混合物分数,在混合物分数方程源项中体现液相对燃烧模型的影响。模拟得到的气相温度分布和液相统计值均和实验数据较好吻合,验证了该燃烧模型对稀液雾扩散类型火焰的适用性。分析瞬时图发现,该稀液雾火焰的最高燃烧温度往往位于当量混合物分数附近,在出口下游20倍直径处火焰完全点着,此处上游FPV模型能给出局部点火熄火现象。蒸发作用在剪切层和点火区域较强,而液滴和火焰作用较弱,单一液滴很少被火焰包围。     

S形进气道锤击波荷载特性分析

利用基于改进的延迟分离涡模拟（IDDES）方法,对亚声速和超声速来流条件下某S形模型进气道进行了非定常计算,研究了发动机喘振所产生的瞬时高压波形对锤击波传播规律的影响.结果表明:锤击波产生后沿进气道迅速向前传播,运动过程中锤击波的运动速度基本保持不变,但强度不断增强.同时受气流离心力的影响,S形进气道弯曲段半径较大一侧壁面受到的锤击波气动荷载值更大.发动机喘振所产生的瞬时高压的加载梯度增加使得锤击波传播速度及强度增强,而压力卸载方式对锤击波强度的影响不明显.在亚声速和超声速来流条件下,增加瞬时高压峰值均使得锤击波荷载强度显著增强,并近似符合二次函数分布规律,而且超声速来流条件下锤击波强度较亚声速来流更强.