数值耗散对压气机流动分离涡模拟的影响研究

分离涡模拟DES是压气机流动模拟中常用的高保真湍流模式。为了使DES准确解析湍流,数值耗散必须限制在合理范围内。然而,当前的压气机流动DES类研究工作中仍然普遍采用高耗散的迎风格式。本文首先基于DES类方法计算的各向同性衰减湍流结果,定量比较了多种不同数值格式的耗散,证实了高耗散迎风格式严重低估中高波数湍流能量。高阶重构格式可以一定程度上改善该问题,但能量耗散仍然过高。本文在高阶重构的基础上,进一步引入自适应耗散函数修改Riemann求解器,构造了自适应耗散格式。该格式在全波数范围都能准确地预测湍流能谱。将该格式配合DES类方法模拟跨声速离心压气机流动,其预测的压比相比于三阶迎风格式,更加接近实验结果。此外,自适应耗散格式显著提高了中小尺度流动结构的分辨率。分析表明,在使用DES类方法模拟压气机流动时,有必要采用数值耗散较低的离散格式,以准确预测压气机总体性能和流动结构。本文构造的自适应耗散格式是一种良好选择。

基于高精度格式的翼型失速特性脱体涡数值模拟

该文采用脱体涡模拟方法(DES),引入低耗散近似黎曼解HLLC格式并结合五阶WENO高精度数值格式对NACA0015后缘失速翼型在大攻角及失速状态下的流场进行数值模拟。通过与风洞实验结果进行对比,讨论高阶精度数值格式对DES方法失速气动特性的预测效果。结果表明,HLLC+五阶WENO格式与二阶Jameson中心差分格式相比能够提高DES方法翼型失速特性的预测精度,HLLC+五阶WENO格式精度高、耗散性较低,能够捕捉更小尺度的湍流涡结构。

用DES数值模拟分离绕流中的旋涡运动

脱体涡模拟(Detached-Eddy Simulation,DES)是近年来出现的一种结合雷诺平均方法和大涡数值模拟两者优点的湍流模拟方法.采用基于Spalart-Allmaras方程模型的DES方法,数值求解Navier-Stokes方程,模拟绕流发生分离后的旋涡运动.其中空间区域离散采用有限体积法,方程空间项和时间项的数值离散分别采用Jameson中心格式和双时间步长推进方法.通过模拟圆柱绕流以及翼型失速绕流,观察到了与物理现象一致的旋涡结构,得到与实验数据相吻合的计算结果.

基于DES射流式涡发生器数值模拟研究

为了更好应用射流式涡发生器抑制风力机叶片流动分离,采用分离涡(DES)方法对平板上射流式涡发生器流进行数值模拟,分析射流式涡发生器旋涡运动规律以及吹风比对射流孔下游湍流场的影响。通过对射流孔下游大尺度湍流相干结构研究表明,发卡涡及流向涡是射流湍流场中的典型相干结构,成熟发卡涡涡腿外侧诱导出的次生流向涡对近壁区能量的交换有着重要作用。对比不同吹风比时射流孔下游流场信息表明,吹风比越大,射流孔下游旋涡尺度越大,射流扰动范围越大,射流壁面摩擦力损失降低,但流场中压差损失变大。

电磁力控制圆柱体三维绕流场的脱体涡模拟

利用Maxwell方程组直接数值计算表面包覆电极与磁极圆柱体产生的电磁力分布,将其加入到动量方程中,采用脱体涡模拟(DES)方法,在雷诺数Re=3 900时,对电磁力作用下圆柱体在弱电解质中的绕流场结构及其升阻力特性进行了数值模拟与分析.结果表明,电磁力作用可提高圆柱体边界层内的流体动能,抑制流动分离的产生,减弱圆柱绕流场的三维特性,在电磁力作用参数达到某个临界值后,在圆柱体后方产生射流现象;同时,随着电磁力作用参数的增大,圆柱体压差阻力及其总阻力减小,但摩擦阻力增大,而且电磁力的作用还可以显著减小升力脉动幅值.

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

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

脱体涡模拟方法研究及其应用

基于现有的非定常雷诺平均N-S方程的求解方法,研究与发展了脱体涡模拟方法(DES)。在紊流附面层内用RANS方法,湍流模型采用S-A模型,在其他区域结合S-A模型运用Smagorinski的LES模型,并将S-A模型中的物面耗散项进行适当的改进,使得改进的模型既充当了RANS中的S-A紊流模型,又充当了Smagorinski大涡模拟模型。为验证所发展方法的有效性,将研究的DES方法应用到一些典型流场计算中,并与已有结果或实验结果进行比较。

非定常等离子体激励对大攻角流动分离控制的数值模拟研究

应用交流电(alternating current,AC)介质阻挡放电(dielectric barrier discharge,DBD)等离子体流动控制由于其结构简单、响应频率快、可实现实时定量控制等优点,正在成为等离子体流动控制技术的重点研究方向。结合基于分离涡模拟(detached eddy simulation,DES)和等离子体唯象体积力模型的方法研究非定常等离子体激励对NACA0015翼型在攻角为20°情况下流动分离控制。结果表明:非定常等离子体激励在高雷诺数、大攻角下对翼型分离具有明显的控制效果,可以达到增升减阻目的,且流动控制效果比定常激励效率更高;非定常等离子体激励流动控制与定常等离子体激励流动控制机理不同,非定常等离子体激励通过促进分离区内速度脉动,对流场产生非定常的干扰,使得分离剪切层提前失稳,增强流场涡结构的掺混,从而抑制流动分离。

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA(Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.

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

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

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

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

压气机叶栅角区分离中速度PDF双峰现象的DES数值模拟

为了研究在实验中发现的角区分离区域内速度概率密度分布(PDF)双峰现象并探索其机理,利用脱体涡模拟(DES)对压气机叶栅中角区分离流动进行数值模拟,成功模拟出双峰现象,并讨论了双峰现象的可能机理。DES模拟结果中,双峰点位置以及内外峰随双峰点位置变化而此消彼长的发展趋势与实验结果一致。对来流攻角脉动特性的分析表明,角区分离中双峰现象与叶栅前缘附近非定常来流和角区内分离流动的相互作用有关。

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

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

来流攻角对叶片脱体涡脱落频率的影响

对于绕流问题的研究,涡的脱落频率是一个重要的参数,它是用来衡量流场的周期变化程度。为了研究来流攻角对叶片后脱体涡脱落频率的影响,分别对相同速度、相同的叶片长度、来流攻角分别为±30°、±45°、±60°的叶片进行数值模拟。计算表明:脱体涡的脱落频率随着叶片安装角度绝对值的增大而减小,安装角度的正负对脱落频率的影响很小,可以忽略不计。

基于两方程湍流模型的DES方法在超音速圆柱底部流动计算中的应用

采用基于M enter k-ωSST两方程湍流模型的脱体涡模拟(D etached Eddy S im u lation,DES)方法,求解N av ier-Stokes方程,数值模拟了超音速下圆柱底部的大分离流动。脱体涡模拟在近物面区采用雷诺平均方法,在其它区域采用Sm agorinsk i大涡模拟方法,兼具前者计算量小的优点和后者能模拟大分离湍流流动的优势。与雷诺平均方法的计算结果进行对比发现,DES方法可以更好地模拟分离涡的发展,得到的底部径向压力分布的时间平均值与实验值吻合。

涡结构对小颗粒在圆管背风面碰撞分布的影响

针对锅炉管背风面飞灰沉积过程,采用数值模拟的方法进行了研究.对绕柱气相湍流的模拟采用了较新的分离涡模拟方法（detached eddy simulation, DES）,通过对粒径为1～10 μm的颗粒在圆管背风面碰撞分布规律的计算发现背风面近壁区存在的主涡、二次涡和分离点附近小涡对颗粒在背风面的碰撞均有一定的影响.当速度较小时,随着粒径的增大,颗粒在背风面碰撞的质量逐渐增加,大颗粒在主涡和二次涡作用下均有较大的碰撞质量流量,而且二次涡的作用明显强于主涡;速度增加后,由背风面3种涡结构引起的颗粒碰撞量均有不同程度的减少,此时主涡的作用较强,而且二次涡引起的大颗粒碰撞量有明显减少.通过与实验测试结果进行对比检验,证明了数值模拟结果的准确性.

圆柱绕流噪声预报的流场与声场模拟方法对比研究

以三维圆柱为研究对象,使用Lighthill声类比法研究其绕流发声问题。第一步进行流场计算,分别用大涡模拟(LES)、脱体涡模拟(DES)和瞬态雷诺平均法(URANS)模拟声源区流场,通过对比流场压力和涡量等参数,据此选取合适的流场仿真方法;第二步用基于Lighthill方程FW-H积分法和边界元法预报远场直发声,通过和Revell试验结果比较,分析各种计算方法差别。研究表明:进行流场仿真时LES计算结果最好,IDDES法在保证计算精度条件下能有效减少流场网格数量,URANS法误差很大;进行辐射噪声预报时,FW-H积分法和边界元法基本相同。

DES模型在压气机亚音转子中的应用探讨

本文运用分离涡模拟(DES)方法研究了不同工况下压气机亚音转子的流动情况,分析了其时均与瞬时流场中顶部间隙泄漏流动和根部角区的流动分离。通过与S-A模型计算结果的对比表明,DES模型在模拟顶部泄漏流动及二次泄漏、泄漏流在转子下游与尾迹的干涉时能够捕捉到更强的旋涡结构,在模拟转子根部角区的分离时也能获得更为丰富的流动现象.对不同工况的DES计算表明负荷的上升会使泄漏涡形成的位置向上游移动,从而导致并加剧二次流动,并对叶栅下游泄漏涡与尾迹的干涉产生影响。对设计工况下瞬时流场的分析表明,泄漏涡在叶栅下游体现出周期性的强弱变化,近叶根分离区也体现出明显的非定常性。

一种风力机专用翼型气动性能的三维数值模拟

为了准确预测风力机专用翼型在大攻角状态下的气动性能,运用脱体涡模拟(detached eddy simulation)方法对瑞典的FFA-W3-241翼型较大攻角范围内的气动性能进行三维数值模拟,对该翼型前缘粗糙状态下的气动性能进行预测.计算结果表明:建立翼型的三维模型,运用DES模拟风力机专用翼型气动性能的方法在线性区有很高的预测精度,在失速发展区的计算精度达到工程实际与研究的要求,在深度失速区有一定的预测精度,可用于定性分析.前缘粗糙度对FFA-W3-241翼型的气动性能有重要影响,前缘粗糙度的增加使FFA-W3-241翼型的最大升力系数下降了27.8%,失速过程趋于缓和;翼型在线性区和深度失速区对前缘粗糙度不敏感,在失速发展区对粗糙度敏感.

偏航状态下空腔流动的数值模拟

采用基于S-A(Spalart-Allmaras)湍流模型的脱体涡模拟方法(Detached-eddy simulation,DES),数值模拟研究了偏航状态下的空腔流动。首先,运用DES方法对非偏航状态下的空腔流动进行数值模拟,并和文献结果进行比较,验证该方法模拟空腔流动特性的能力;然后,运用DES方法数值研究偏航状态下的空腔流动,给出声压级的变化。数值研究表明:不同偏航角对空腔底部中心线处声压级产生很大的影响;与非偏航状态相比,在空腔底面中心线的某些位置处,给出了声压级随偏航角的变化规律为:声压级先是随着偏航角增加而缓慢增加,然后大大减小。