三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

壁面喷流可控环量机翼的数值计算

本文提出一种适用于可控环量机翼研究的势流/粘性流的数值解法,包括钝体分离绕流的一种势流解,以及改进的Spalding统一边界层理论的解。通过对具有壁面切向喷流可控环量圆柱绕流和类椭圆翼绕流所作出的数值计算表明,在其计算模型确定后,就不必依赖实验数据来模拟实际问题,计算结果与实验数据符合较好。     

基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.     

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

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

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

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

折叠V形钝体火焰稳定器原型的稳燃机制研究

为了掌握折叠V形钝体这一新型火焰稳定器的稳燃机制,采用风洞实验与混合雷诺平均/大涡模拟结合的方法对常压条件下10组不同折叠角的折叠V形钝体火焰稳定器原型进行了冷态实验与热态数值模拟研究。其中冷态实验风速10~50m/s,热态数值模拟工况为进气温度700K、来流速度50~150m/s,当量比0.3~1。通过PIV实验测得了折叠V形钝体回流区冷态流动特征,并结合数值模拟方法研究并讨论了折叠V形钝体绕流尾迹区旋涡特征;通过数值模拟研究了燃烧条件下、折叠V形钝体下游槽向驻涡结构发展变化规律。研究初步探明折叠V形钝体总压恢复系数、阻力系数、抗吹熄能力与燃烧效率优于标准钝体的原因,并发现折叠角角度较小的折叠V形钝体,尾迹中心区域无典型回流区,但仍能通过槽向驻涡结构形成驻涡稳燃速度型,该速度型有利于冷态绕流流速的恢复与热态燃气的加速。     

圆柱尾流中的二次分离及环量衰减问题

在圆柱绕流的非定常流动中,相应于旋涡的周期性脱落柱后的二次分离是周期性间断发生。二次涡的影响不能忽略。但对长时间流动的模拟还必须考虑尾涡的环量衰减。模拟Bluff Body流动有效的离散涡模型应包括二次分离及尾涡衰减两种机因。本文数值计算结果证明了上述结论。     

钝体尾迹的稳定性研究及流动控制探讨

系统地回顾了研究钝尾尾迹稳定性的主要方法，并利用局部绝对／对流稳定性分析及整体稳定性分析两种方法，研究了圆柱绕流定常和时均尾迹流场的稳定性。采用低维Ｇａｌｅｒｋｉｎ方法研究了尾迹对圆柱流向振动的响应问题，揭示了频率锁定，涡街的激发和抑制等现象。     

突然起动圆柱的非定常绕流计算

本文用分区方法定量模拟高Re数情况下突然起动圆柱的初期对称流动。绕圆柱的二维不可压流场可分为无粘区、附着粘性区和分离区。无粘区有解析解，用非定常层流边界层方程解附着粘性区，用离散涡方法处理随时间变化的分离区。试用了混合的Rankine和Lamb涡模型及二次离散涡。给出了Re=300、9500两个算例，基本符合有关实验结果。     

圆柱突然起动的快速涡解法

 文中将Chorin的涡方法、格子涡(VIC)方法和快速Fourier变换(FFT)相结合,发展了一套快速有效的求解高Re数下非定常N-S方程的算法(快速涡方法)。通过对二维圆柱突然起动问题的数值计算,得到流动分离及其演化过程(包括二次涡及二次涡对的生成);此外还得到回流区对称轴线上速度分布。这些结果与实验数据及其它数值解相比,是令人满意的。该程序对其它高Re数非定常流动问题也是适用的。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

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

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

基于被动控制的D型钝体减阻大涡模拟研究

钝体减阻在航空航天等多个领域内具有重要的应用潜力,也是重要的基础研究热点问题。为了精确捕捉流动控制的流动细节,并发展效果优良的流动控制减阻方案,论文围绕抽象出的D型钝体,采用大涡模拟方法,开展了被动控制减阻高精度数值模拟研究。首先基于前期研究成果,对D型钝体尾迹区剪切层附近放置一个光滑小圆柱的被动减阻方法开展了数值模拟,发现总阻力减小17.7%,与试验结果吻合很好,同时数值预测的速度场分布也与试验结果吻合良好。在此基础上,进一步提出了采用齿槽型表面结构的小圆柱对D型钝体尾迹区进行扰动,并开展了数值验证,发现总阻力减小21.4%,优于前期的减阻方法。最后研究了在三种雷诺数工况下两种小圆柱扰动情况下的减阻效果,均表现出良好的减阻效果,两种小圆柱扰动下总阻力最大降幅分别为19.6%和23.1%,同时基于大涡模拟计算结果对减阻流动机理进行了探讨。上述研究结果表明,通过进一步优化流场结构,可以得到更优的流动减阻方案。     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。