绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

基于改进SST模型的分离流动数值模拟

对SST湍流模型进行了改进,并通过对典型分离流动进行数值模拟,来检验和提高模型预测分离流动的能力.基于亚声速分离流动,提出了提高雷诺应力的模拟精度和分离流修正的改进方法,并进行了对比研究得出结论;在亚声速分离流动分析结论基础上,采用了可压缩性修正方法,开展了跨声速、超声速和高超声速激波/边界层干扰分离流动的数值模拟研究.结果表明:提高雷诺应力的模拟精度和采用分离流修正明显地提高了SST湍流模型对分离流动的模拟能力;适当地可压缩性修正对超声速和高超声速分离流动的计算精度有改善作用.     

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。     

三角翼大迎角绕流的数值模拟

采用FINETM/HEXA非结构网格流场计算软件包对三角翼飞行器的低速大迎角绕流进行了数值模拟.结果表明,计算值与实验数据符合较好,采用Agglomeration 方法能大幅度提高收敛速度,而网格自适应方法用总体不大的网格单元数可较精细地模拟流场.本文还分析了前缘分离涡破裂前后的流动现象和旋涡横截面流线图谱的变化规律.     

k-ω SST两方程湍流模型中参数影响的初步分析

Menter的k-ωSST两方程湍流模型在流体力学计算中有良好表现。利用均匀试验设计方法对k-ωSST模型中八个参数在小攻角亚声速、跨声速和大攻角三种典型的翼型绕流流场计算中的影响规律进行了分析。分析结果表明在附着流中参数a1是对流场的影响最大的参数;而在分离流中,参数a1对流场的影响与湍流产生与耗散等项的作用相当,模型中的四个参数:σω2、β2、β*和a1的取值都会对计算结果有较大的影响。     

亚临界圆柱绕流的DES方法比较

针对地面结冰试验设施喷雾系统尾迹湍流模拟需求,考察了3种脱体涡模拟(DES)方法对三维亚临界圆柱绕流预测的准确性,比较和分析了瞬时流动特征和流场统计量。研究发现:1)从瞬时流动特征来看,k-ω雷诺平均NavierStokes(RANS)控制的区域对剪切应力输运k-ωDES(SSTk-ωDES)求解的准确性有显著影响;2)从流动参数统计量的误差范围来看,回流区长度和流向最小速度的预测质量能在一定程度上反应圆柱绕流数值模拟的准确性;3)综合比较,SSTk-ωDES预测值与试验值和大涡模拟(LES)数据最为吻合,具有应用于地面结冰试验设施喷雾系统尾迹湍流计算的潜力。     

两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

三角翼大攻角绕流数值模拟研究

采用Harten-Yee的二阶精度隐式TVD格式对N-S方程进行数值模拟,研究了后掠角为65°和70°的三角翼绕流的流场特性。计算结果表明：TVD格式中限制器函数的选取对于三角翼绕流结构的计算结果具有较重要的影响。结合张涵信的旋涡定性分析理论对旋涡的发展过程进行研究,此外,还对二次涡破裂现象和网格分布影响等问题进行了探讨。     

旋涡破裂引起的流动非定常特性

采用隐式近似因子分解法,通过求解三维、非定常、可压缩完全层流N-S方程,对绕大后掠三角翼旋涡破裂的流场进行了数值模拟。结果表明：当三角翼背风区流场发生旋涡破裂时,各流动参变量是非定常的,同时旋涡破裂点的位置沿轴向在一定范围内前后摆动。     

基于图像的测量方法在工业风洞跨声速流测量中的最新进展

风洞中非定常复杂流场的实验研究要求先进的测量技术.基于图像的测量技术中最重要的是测量平面流速度场、平面压强分布、模型位置和变形、模型温度以及定量的高速流可视化等技术.DLR(德国宇航研究院)对这些技术的应用包括从低速流到跨声速流、从增升装置到螺旋桨和旋翼、从弹射装置和水塔储水罐尾迹流旋涡到三角翼上涡破裂现象等的研究.由于跨声速风洞的特殊环境,将基于图像的测量技术用于跨声速流要求专门的技术开发和有经验的科学家.给出了DLR空气动力学和流动技术研究所将基于图像的测量技术应用于跨声速流研究的最新进展.     

三角翼后缘喷流对前缘涡破裂位置的影响

采用高分辨率的N-S方程数值模拟方法,研究了亚声速来流条件下对称喷流、差动喷流、矢量喷流对60°后掠角三角翼前缘涡破裂位置的影响。结果表明,对称喷流条件下,喷流与来流的速度比影响前缘涡破裂位置;差动喷流情况下,速度大的一侧喷流对涡破裂位置起主要影响;矢量喷流情况下,喷流方向对涡破裂位置影响显著,喷流方向与涡轴方向一致时,涡破裂位置延迟量最多。     

低速大迎角三角翼旋涡流动的有粘与无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Euler方程和雷诺平均Nayierstokes方程，应用Beam－Warming近似因子分解格式及其对角化形式，采用BaldwinLomax代数湍流模型。计算模型为70°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Euler方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

折流环形燃烧室流动过程数值模拟与分析

为深入研究折流式环形燃烧室燃烧过程,在任意曲线坐标系下发展了一套可对湍流燃烧进行数值分析的数值模拟软件。通过求解偏微分方程生成数值研究所需结构化网格,使用混合差分格式对气相控制方程进行离散并基于SIMPLE算法进行求解,分别采用RNG k-ε与Wilcox k-ω模型分析不同湍流模型对计算结果的影响,近壁区域流动采用三层壁面函数法加以处理。数值模拟结果表明,引入燃烧室的新鲜空气在扩压器出口流动方向发生折返,在扩压器下侧形成了较大的回流区域,导致较大的总压损失。不同紊流模型预测结果与试验数据符合较好,其中Wilcox k-ω最为准确,更适用于折流燃烧室燃烧过程的数值分析。