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

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

Gurney襟翼对某型客机流动控制数值模拟

为改善某型客机的起降性能,通过在机翼尾缘加装Gurney襟翼,对流场进行了数值模拟。对该客机机翼的控制翼型安装不同高度的Gurney襟翼进行数值模拟,结果表明安装Gurney襟翼可以提高多段翼型的升力系数和阻力系数,但会增强尾迹流动的不稳定性。将不同高度的Gurney襟翼应用于该客机的简化模型,机翼的大部分区域符合二维翼型研究得出的流动控制规律;在机翼外侧区域,Gurney襟翼使机翼附近流场中的翼尖涡发生了一定的变化。数值模拟的结果还表明,Gunney襟翼可以提高客机的升力系数,而且不会给飞机流场带来明显的改变。     

基于分离涡模拟的起落架气动噪声研究

用基于SA湍流模式的非定常RANS、分离涡模拟(DES)和延迟分离涡模拟(DDES)分别对四轮基本起落架模型进行了数值模拟,并根据所得的非定常流场计算了表面声压级分布和声压谱。三种方法所用的时间成本大致相同,URANS略低于其它两种,而DDES由于在附着流动区更好地保持RANS特性,故时间成本略低于DES。计算显示,非定常RANS在附着流动区能够得到合理的结果,但不能准确刻画分离流动的流动形态。DES和DDES都能较好地刻画起落架绕流的定常和非定常特性,DDES的结果略好于DES。因为起落架流动属于大分离流动,所以DES也能够得到相对正确的结果。研究结果验证了分离涡模拟在起落架大分离非定常流动预测与噪声预测中的可行性,对减小起落架噪声的方法研究具有一定的意义。     

后缘襟翼侧缘加装挡板降噪数值模拟研究

在飞机的起飞着陆过程中,机体噪声是主要的噪声源,后缘襟翼侧缘噪声作为机体噪声的重要组成部分之一,其降噪技术的发展一直受到广泛关注。采用分离涡(Detached Eddy Simulation)湍流模型模拟襟翼侧缘的非定常流场,分析并探究了襟翼侧缘流场中双涡结构的产生及其演变过程的主要特征。基于上述结果,采用了被动流动控制法,在襟翼侧缘加装最大高度为一倍襟翼厚度的不规则挡板结构,应用FW-H方法计算了远场噪声频谱特性以及指向性并分析了与基准构型的差别,探究了加装襟翼挡板结构后对流场以及声场特性的影响。模拟结果表明,加装挡板结构对襟翼的气动特性影响较小,挡板结构改变了侧缘的流场形态和侧缘涡的结构,延迟了侧缘涡系的发展进程,使得涡系的融合位置后移,涡系融合破裂产生脉动压力的区域远离襟翼吸力面,从而达到降低噪声的效果。此外,加装挡板结构后,噪声仍然具有一定的偶极子指向特性,偶极子轴垂直于襟翼弦线,在襟翼的大部分方位噪声幅值降低1 dB~4 dB,降噪效果主要体现在中高频段的宽带噪声。     

舰船空气流场数值模拟方法研究

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法,兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势。论文分别采用RANS和DES方法对舰船空气流场进行数值模拟,通过所得流场结果与实测数据的对比分析,说明DES方法对大分离流动具有更强的模拟能力,更加适用于舰船空气流场预报。     

基于SST-DDES方法的孤立气膜孔流动研究

采用基于SST湍流模型的延迟分离涡模拟(DDES)方法,对吹风比为0.5的平板孤立方孔横流射流进行数值模拟研究,并与实验测量结果进行对比;同时,基于捕捉到的流场大涡拟序结构,分析了射流与主流掺混及损失机理。研究结果表明:气膜冷却流场中存在复杂的大涡拟序结构,这种强三维流动对冷气覆盖效果和掺混损失起主导作用;与传统的RANS和URANS方法相比,SST-DDES方法不仅在射流下游平均流场的预测方面具有更高精度,而且还能捕捉更加丰富的射流/主流掺混过程流场细节及其时空演化特性,将SST-DDES类RANS/LES混合方法的应用领域扩展到横流射流类复杂流动的研究中具有较好的可行性。     

基于自适应重叠网格的三角翼跨声速流场计算

在自适应重叠网格系统下,数值求解非定常Navier-Stokes方程,开展了钝前缘三角翼跨声速流场的计算研究。目的在于考察交叠网格系统下,不同迎角的跨声速来流条件时流场细节的捕捉能力。其中,网格方面采用了贴体网格块精确描述机体外形,采用与之交叠的可自适应的直角网格捕捉脱体涡系的发展变化及涡与激波的干扰;求解涡黏性计算方面,采用了Spalart-Allmaras(SA)及其对应的Detached Eddy Simulation(DES)模型。基于以上描述的方法,针对钝前缘三角翼在来流马赫数0.85不同迎角状态进行了计算,比较了雷诺平均与DES计算的结果差异。计算结果表明,在重叠网格系统下,网格构建简便,适用性好,计算具备一定的数值精度;对于大迎角状态,DES方法能够有效地模拟脱体涡系的发展变化,获得更好的计算结果。     

冲击射流流动换热超大涡模拟研究

为了准确预测发动机热端部件中广泛采用的冲击射流冷却复杂的流动和换热特性,发展了基于BSL k-ω模型的超大涡模拟(VLES)高精度模拟方法,并对高雷诺数Re=4×104,两种不同射流距离2和6的单孔冲击射流及三孔冲击射流这一经典的流动传热问题进行三维非稳态高精度数值计算。同时,将分离涡方法 (DDES)和k-ωSST,RNG,Transition SST三种RANS方法的数值模拟和开发的超大涡模拟(VLES)方法进行对比。研究表明,VLES方法均能够准确捕捉冲击射流流场的复杂非稳态流动及传热特征,包括自由射流区、壁面射流区小尺度涡系和大尺度湍流结构的演化和破碎,同时冲击壁面的换热系数计算结果与实验值吻合较好。DDES方法未能准确捕捉流场复杂的小尺度湍流结构,壁面换热计算结果与实验值差异较大。RANS方法计算的换热结果与实验数据差异最大,基本未能预测到壁面换热特性。在相同的计算网格和计算方法下,VLES方法计算结果优于DDES方法,DDES方法一般好于RANS方法。这表明新开发的VLES方法能够准确地计算冲击射流相关的流动及换热问题。     

高升力翼型复杂流动模拟中湍流模型性能评估

采用软件Fluent中工程常用的7个涡黏湍流模型对某种高升力3段翼的降落阶段绕流进行了数值模拟.通过和试验及经过试验校验过的延迟脱体模拟（delayed detached eddy simulation,DDES）结果进行详细的对比分析,包括翼型压力系数、马赫数、涡量场和湍动能分布等,系统地研究了常用涡黏湍流模型对该高升力翼型的模拟性能.结果表明：对于平均流场,standard k-ω模型的模拟性能最好,能较好预测翼型压力系数、襟翼分离区位置和大小等;SST （shear stress transport）k-ω模型性能也较好,一方程SA （Spalart-Allmaras）模型和四方程v²-f模型具有相近的性能,而k-ε系列模型性能相对较差.对于湍流场的模拟性能,和上述平均流场对比所得结论相一致,但所有湍流模型均未预测出襟翼分离区附近湍动能“最大”的分布特征.     

应用等离子体进行翼型前缘分离控制的数值模拟研究

将等离子体对中性气体的作用模型化为彻体力矢量,求解带源项的Navier-Stokes方程,数值模拟了等离子体激励器在NACA0015翼型大迎角下的分离控制效果,彻体力为净电荷在外加电场作用下的电场力.解拉普拉斯方程得到外加电场分布,等离子体中的净电荷分布由泊松方程给出.为了较好地模拟分离涡的发展,采用了雷诺平均与大涡模拟相结合的脱体涡模拟(Detached Eddy Simulation)方法.通过与实验结果对比,发现该模型能较好地模拟等离子体激励器的控制作用.     

分区LES/DES混合方法及缝翼三维流场模拟

提出一种分区大涡/脱体涡模拟混合方法,采用高阶空间有限差分方法和二阶隐式LU-SGS时间推进方法,对有限展长缝翼的三维流场进行了数值计算。该混合方法吸收了脱体涡模拟方法和大涡模拟方法各自的优点,在非核心区域使用脱体涡方法进行计算,相对于大涡模拟方法具有较高的计算效率;相对于分区脱体涡模拟中采用的雷诺平均方法,在计算量没有明显增加的条件下提高了对流场各尺度流动的模拟能力。计算结果与风洞实验数据吻合良好,同时缝翼流场主要观测点的压力脉动数据为后续缝翼噪声分析和低噪声优化提供了基础。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

DDES方法在复杂旋翼流场计算中的应用

应用延迟脱体涡模拟（DDES）方法,进行了复杂分离流动旋翼流场的数值模拟研究。为控制计算规模,更好地模拟旋翼桨叶近壁面附近区域流动分离及远场尾迹发展,采用了贴体与背景自适应直角网格嵌套的网格系统;同时为提高计算速度,采用了基于分布式存储的数据并行模式,实现了贡献单元搜索、洞切割和流场计算的并行化。计算工作包含一个共轴双旋翼悬停和单旋翼下降状态,目的是分别讨论DDES方法在旋翼表面附近空间区域是否存在大范围分离时计算结果与雷诺平均Navier-Stokes结果之间的差异。计算表明,对于下降状态,DDES方法预测出了更强的尾迹区分离流动;对于桨叶表面基本为附着流动的双旋翼中等桨距角情况,两种方法预测的流场结构相似,仅在下旋翼桨叶气动力和旋翼桨根下方空间区域的诱导速度计算结果上存在微细差异。     

用DES数值模拟具有横向喷流的紊流流场

DES(DettachedEddySimulation)是把RANS(雷诺平均Navier Stokes方程)方法及LES(大涡模拟)方法结合起来的模拟有脱体涡的紊流流场的数值模拟方法。其主要思想是在物面附近解雷诺平均Navier Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文用DES及混合非结构网格数值模拟复杂紊流流场,算法采用Osher逆风格式,对超音速钝头体、三角翼、某型飞机及带喷流的导弹绕流用该程序进行了数值模拟。结果表明用DES模拟有分离的紊流流场是高效的。     

细长体大迎角湍流流场的数值模拟

通过数值方法对大迎角细长体超声速流场的模拟，探求正确模拟导弹大迎角绕流湍流流场的简单而有效的方法。数值模拟求解一般曲线坐标系下的三维可压缩Navier-Stokes方程，时间离散采用Euler向后差分，无粘项的空间离散采用二阶TVD格式，分别研究了Bald-win-Lomax代数模型及其修正形式(BLDS)对大迎角分离流动的模拟能力。数值试验表明修正的Baldwin-Lomax模型更精确地预测了流场的旋涡与分离情况，给出合理的表面压力分布，因而更准确地模拟了存在横向分离的导弹流场。     

基于IDDES方法和γ-Re_θ转捩模型的粗糙颗粒诱导高速边界层强制转捩模拟

采用γ-Re_θ转捩模型与IDDES (Improved Delayed Detached Eddy Simulation)相结合的方法对BAM6QT (Boeing/AFOSR Mach-6 Quiet Tunnel)中的马赫6来流条件下粗糙颗粒诱导转捩情况进行了数值模拟研究,通过与试验测量的压力脉动均方根值、脉动主频和边界层内的皮托压分布的定量对比及与文献中DNS (Direct Numerical Simulation)流场结构的定性对比,表明该方法可以捕捉粗糙单元诱导出的流向涡,能够模拟颗粒前缘分离激波和弓形激波之间的震荡现象,能够模拟流向涡向下游的发展失稳过程及其脉动发展过程。计算结果表明流向涡结构在粗糙颗粒下游40倍颗粒直径位置开始破碎,非定常扰动能也在该点附近增长达到最大值。粗糙单元诱导出了明显的条带涡结构,其中低速条带处于边界层上层,在向下游发展的过程中逐渐扩散至边界层外缘并耗散掉。高速条带位于边界层底部,在向下游发展的同时往展向两侧拓展,最终展向上多条高速条带接触并互相耦合,导致最后条带涡结构的破碎和尾迹区边界层的转捩。     

Gurney襟翼增升效应数值模拟

飞机的增升系统是影响其起降性的重要因素,为了系统地研究Gurney襟翼增升效应,本文采用数值模拟方法与结构网格求解不可压缩的Navier-Stokes方程对NACA0012翼型和机翼进行了研究。通过改变Gurney襟翼的高度、弦向位置来研究其对翼型气动特性的影响,并将结果应用于矩形机翼,计算结果表明:在一定范围内,Gurney襟翼越高、越靠近后缘则增升效果越好,在三维机翼上也具有相似的效果。     

支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

弹丸旋转空气动力效应非定常数值模拟

准确计算马格努斯力和力矩对旋转弹箭设计、弹道计算和稳定性研究都至关重要。采用非定常雷诺平均Navier-Stokes(RANS)方程对高速旋转的SOCBT弹丸进行了数值模拟,纵向气动特性在整个攻角范围内都较好地与试验值保持一致,而侧向气动特性马格努斯力和力矩在小攻角范围与试验数据吻合较好,在大攻角范围却存在一定的差异。采用延迟分离涡模拟(DDES)方法的计算结果有较为明显的改善,对比研究表明分离点位置对马格努斯效应有着显著影响。表明DDES方法对于提高旋转弹箭马格努斯效应的数值模拟精度有较大的潜力。     

Gurney襟翼对某水陆两栖飞机纵向气动特性的影响

为了掌握Gurney襟翼对飞机的影响机理及增升规律,改善某大型水陆两栖飞机的起降特性,利用数值计算方法分析研究了Gurney襟翼对单缝襟翼流场及气动特性的影响,研究了不同高度和偏度的Gurney襟翼对某大型水陆两栖飞机全机气动特性的影响。研究结果表明,采用Gurney襟翼明显提高了该型水陆两栖飞机的升力系数和抗浪性。