高阶精度格式WCNS-E-5在亚跨声速流动中的应用研究

采用高阶精度非线性紧致加权格式WCNS-E-5和Baldwin-Lomax模型,求解雷诺平均Navier-Stokes方程,开展了典型翼型与机翼的湍流流动数值模拟研究.对方程中粘性项采用的四阶精度差分近似以及网格导数求解与边界格式的四阶精度,保证了高精度算法的实现.计算结果表明:本文算法能够准确地模拟这些翼型与机翼的亚跨声速流场,得到与实验测量十分吻合的壁面压力分布,计算结果对网格的依赖性小.     

Pω增强型k-ω湍流模型在三角翼旋涡流动的应用

在三角翼旋涡绕流数值模拟中,标准 Wilcox k-ω湍流模型生成项未考虑旋度的影响而导致预测的旋涡强度较弱。通过引入探测因子区分剪切层和涡核,在旋涡流动的高旋度区域增加ω方程生成项的方法,基于结构化网格上的 RANS 求解器,加入了 Pω增强型 k-ω湍流模型,对绕尖前缘三角翼亚声速和跨声速旋涡流场进行了数值模拟。计算结果与 NASA 的 NTF 风洞和 DLR 的 DNW-TWG 风洞试验数据进行了对比分析,结果表明：不论在亚声速还是跨声速自由来流条件下,Pω增强型 k-ω湍流模型计算的压力分布、涡破裂位置均与试验数据吻合良好,准确地预测出了三角翼上翼面的主涡、二次涡结构,特别是跨声速条件下激波干扰导致的涡破裂的临界迎角及涡破裂位置,表明 Pω增强型 k-ω湍流模型在绕三角翼旋涡流动数值模拟中具有良好的适用性。     

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

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

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

两种湍流模型及其可压缩修正在高超压缩拐角中的评估(英文)

采用在高阶精度差分格式对SA湍流模型和SST湍流模型及相应的可压缩修正模型在高超拐角流中进行了评估。可压缩修正方式考虑了密度梯度、压力膨胀和湍流马赫数等方法。为了减小数值误差与模型误差之间的混淆,控制方程的对流项采用了5阶精度的加权紧致非线性格式(WCNS-E-5),粘性项采用了一种半结点/结点交错的4阶中心格式。通过对马赫数为9.22的15度拐角和34度拐角湍流的模拟,考察了原始湍流模型及其修正模型的效果。计算表明:原始SST模型对高超拐角湍流的预测比原始SA模型准确,这种准确主要体现在对分离区预测、再附点附近压力和热流峰值预测上。通过混合采用Catris和Shur等的方法对 SA模型进行可压缩修正可以大大改进模拟效果。在SST模型的可压缩修正方法中,Catris的修正方法最好;考虑压力膨胀修正后得到的分离区远远偏大,比本原始SST模型更差;考虑湍流马赫数的修正方法得到的分离区偏小。本文还给了部分低阶格式的计算结果,高阶格式与低阶格式相比,对分离区大小、再附点附近的压力和热流峰值等的预测准度有所改进。     

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

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

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

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

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

CFD空间精度分析方法及4种典型畸形网格中WCNS格式精度测试

计算流体力学中的精度阶数体现的是数值解随着网格加密向精确解收敛的速度,计算格式、网格质量和边界处理等都会影响计算精度。5阶非线性加权紧致格式(WCNS-E-5)是一种具有较高的精度和分辨率的有限差分格式,对强激波和复杂网格的适应性强。但是,格式的模板自适应选择器和网格的复杂性等都会对最终计算精度造成影响。为了分析精度,首先给出了两种通过数值解反推计算格式精度的方法:一种方法适合已知精确解的情况,另一种方法适合未知精确解的情况。然后在斜交网格、拉伸网格、拐折网格和扭曲网格等4种典型畸形网格中考察了WCNS-E-5格式的精度,考察算例涉及线性波动问题、等熵涡传播问题和超声速边界层问题等。计算结果表明WCNS-E-5在这四种网格中都体现出高阶精度。     

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

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

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

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

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.     

机身前体横截面形状对三角翼涡破裂位置的影响

通过求解三维定常雷诺时均方程,采用剪切应力输运(SST)湍流模型,在亚声速范围内,分别对融合体前体-三角翼组合体和旋成体前体-三角翼组合体流场中翼涡破裂现象进行了数值模拟.模拟结果显示:与旋成体前体相比,在中、大迎角时,融合体前体的分离涡,涡量集中、强度高,进入机翼上方流场后,能够与翼涡密切耦合,彻底地改变翼涡强度的分...     

后缘喷流对三角翼绕流影响的N-S方程数值分析

本文用拟压缩性方法求解不可压流雷诺平均拟压缩Ｎ－Ｓ方程组,对带有后缘喷流的三角翼粘性绕流进行了数值模拟,求解中采用了Ｂｅａｍ－Ｗａｒｍｉｎｇ隐式近似因子分解格式以及ＭＭＬ代数湍流模型。计算结果说明,后缘喷流使涡核压强降低,使涡核速度增大,从而对三角翼前缘分离涡有稳定作用,并能增大上翼面的负压值和下翼面的正压值,从而可以增加部分升力。计算结果还说明,喷口面积或喷流下偏会使上述作用增强。     

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

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

旋涡运动研究的某些进展

 本文简要介绍研究旋涡运动在以下问题上的某些结果:低速不同后掠角三角翼在各个迎角下的九种分离流类型及其边界;应用微分方程定性论与拓扑学对三维分离流与旋涡流的分析;旋涡破裂形态,对三角翼前缘涡破裂的实验研究与理论分析;受控分离与旋涡的干扰,二旋涡的位移、绕转与合并等。     

大后掠翼前缘涡对其颤振特性的影响

 大迎角三角翼的前缘涡不仅可以改善其气动力特性，也会显著影响机翼的气动弹性特性。运用基于Euler方程的非定常气动力降阶模型（ROM）方法，耦合结构运动方程，在状态空间内建立了气动弹性分析模型，研究了70°削尖三角翼的大迎角颤振特性。研究结果显示前缘涡对该机翼颤振特性的影响不可忽略。颤振速度随迎角的增加而大幅降低，迎角α=20°时的颤振速度比α=0°时降低了22％。发现了颤振特性随迎角变化时出现的不连续现象，并揭示了该现象是由于系统颤振分支随着静态迎角的增加发生转移所致。     

NS-DBD激励控制非细长三角翼前缘涡仿真研究

通过在三角翼前缘施加纳秒脉冲介质阻挡放电（NS-DBD）激励唯象学模型，进行了47°后掠角钝前缘三角翼流动控制的仿真。分析了不同迎角下升力和阻力系数的变化、流场结构的变化、以及激励诱导旋涡的演化过程。研究表明:施加无量纲激励频率F+=1.44的NS-DBD激励后，可明显提高三角翼失速前后的升力系数；同时阻力系数也有所增加，变化趋势与实验结果一致。激励在前缘分离剪切层处诱导产生流向涡，改变了前缘剪切层结构，使其向内卷吸；激励后时均流场形成了明显的负压峰值，前缘涡附着线外移，吸力面回流区减小。      

TUNNEL INTERFERENCE IN UNSTEADY POST-STALL EXPERIMENTS

ＴＵＮＮＥＬＩＮＴＥＲＦＥＲＥＮＣＥＩＮＵＮＳＴＥＡＤＹＰＯＳＴ┐ＳＴＡＬＬＥＸＰＥＲＩＭＥＮＴＳＺｈａｎｇＷｅｎｈｕａ，ＤｉｎｇＫｅｗｅｎ，ＨｕａｎｇＤａ，ＬｉＺｈｉｑｉａｎｇ，ＺｈａｎｇＱｉｎｇｌｉ（６ｔｈＤｅｐｔ．ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉ...