ALE方法求解圆柱的涡致振动

数值求解基于拉格朗日-欧拉(ALE)描述的不可压缩流体的N-S方程,计算了较小雷诺数下圆柱的涡致振动现象.N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散,计算网格采用H-O非交错网格系统,并结合分块耦合方法.柱体运动简化为弹簧-阻尼-质量系统,柱体运动方程采用经典龙格-库塔方法求解.通过模拟柱体和流体之间的非线性耦合作用,成功地捕捉到了"锁定"、"拍"和"相位开关"等现象,并与试验数据相吻合.另外,本文详细分析了Re＝200时圆柱尾涡形态、升力、阻力以及圆柱位移等随圆柱自振频率变化的过程,捕捉到尾涡结构变化的频率转折点.     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

高对流Mach数三维混合层转捩特性分析及小激波结构模拟

采用空间大涡模拟方法对超/超混合、超/亚混合两类三维可压缩平面混合层转捩及其全湍流流场进行了研究,认为混合层动量厚度饱和点可作为流场转捩完成的标志.计算所得到的线性扰动波激励下流场转捩拟序结构与随机扰动下自然失稳结构以及文献结果进行了比对,其结果是一致的,表明了引入线性扰动激励来研究流场转捩结构是合理的.同时,本文还在较高对流Mach数流动下得到了三维流场动态小激波结构,其分布具有非对称特性,且形状与实验及直接数值模拟结果相似.不同条件混合层转捩计算表明:高对流Mach数下混合层转捩以Λ涡结构的形成和发展为主导机制,受扰动及对流条件的影响Λ涡结构不尽相同,某些情况下流场出现二维与三维涡结构共存现象.充分发展湍流区域,流场脉动速度分量量级相同,湍流压缩效应随着对流Mach数提高而明显增强.     

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。      

可压缩混合层的涡结构演化与流质混合

为研究可压缩混合层的流动结构,采用七阶精度广义紧致格式离散对流项和用显式八阶精度的中心格式离散粘性项,数值求解了非定常三维可压缩Navier-Stokes方程。用约4亿规模的网格,直接数值模拟了对流马赫数为0.7的超声速可压缩混合层的空间发展流动,获得了自初始流动失稳直至充分发展湍流流动结构的精细演化历程,所得结果表明:大尺度涡结构的生成使得混合层的动量厚度快速增长,并主宰了由被动标量质量分数展示的可视混合厚度的量级,充分发展湍流的小尺度的结构主要使得该量级厚度内的流质混合趋于均匀,对可视厚度增长的贡献非常有限。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

立楔诱导的高超声速流动分离数值分析

高超声速飞行器进气道等关键部件引起的激波与边界层相互作用将导致流动分离,从而改变当地压力分布与局部受热情况,影响飞行稳定性与飞行安全,因此需要对高超声速流动的分离现象进行细致研究。采用高精度5阶特征型WENO格式与3阶TVD型Runge-Kutta方法,求解三维Navier-Stokes方程,对立楔诱导的高超声速激波与边界层相互作用引起的分离流动流场结构进行了细致的数值模拟与分析。结果表明,5阶特征型WENO格式分辨率远高于类TVD格式;Ma=6时得到清晰的激波结构、分离涡结构及其演化过程和壁面极限流线的拓扑结构,证明了WENO格式应用于高超声速分离流动的可行性与高分辨率;对不同来流Mach数的对比证明Mach数的增大抑制流动分离,导致分离涡减小。     

非均匀网格湍流大涡模拟高精度有限体积解法

为准确预测不可压复杂湍流,提出了一种可用于大涡模拟均匀或非均匀网格上的高精度有限体积法。该方法空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动证实了该方法具有近四阶的空间精度;并在此基础上,采用动态Smagor-insky亚格子应力模式,成功地实现了充分发展槽道湍流和后台阶湍流流动的大涡模拟计算,所得结果与直接数值模拟结果吻合良好,且采用非均匀网格可在比均匀网格数少的离散系统上得到同样满意的结果。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.      

三维超声速混合层内标量混合的大涡模拟

对三维对流马赫数0.62的超声速混合层流动的标量混合进行大涡模拟, 控制方程对流项采用五阶精度的WENO格式求解, 小尺度涡的作用采用一方程LDKM亚格子模型处理, 过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了混合层流场大尺度拟序结构以及标量场的演化过程, 研究表明标量混合过程受混合层内涡系演化所控制, 标量场具有明显的三维特征.模拟得到的速度、组分及其脉动的统计时均结果和实验结果相符较好.      

可压缩混合层密度梯度场特性的实验研究

利用纹影技术研究了不同对流马赫数（Mc=0．28,0．38,0．72）可压缩混合层密度梯度场的特性。由纹影图像可以看出：具有可压缩效应（高M。数）的混合层存在不同程度的摆动。将纹影图像对比以后还发现：当Mc〈0．3时,在可压缩混合层中发现了类似于Brown—Roshko的结构,混合层看不出明显的摆动;随着可压缩效应的增强（对流马赫数Mc〉0．3）,混合层开始摆动,并且随着对流马赫数的增大,混合层开始摆动位置向前移动。此外,还观察到不同对流马赫数下大尺度结构出现的区域有不同的特点。     

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程，对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法，亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析，给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析，计算结果与理论及实验吻合。     

可压缩混合层中的大尺度结构及其作用初探

利用粒子示踪图像流动显示技术和粒子图像测速技术(PIV技术)研究了对流马赫数Mc=0.38时可压缩混合层发展早期出现的大尺度结构的形状和特性,并将这些混合现象与湍流脉动场结合起来研究了这些大尺度结构对流动混合效率的作用.将粒子示踪图像和湍流脉动强度场对比以后发现:混合层流动中早期出现的大尺度结构沿着顺时针方向旋转着向前移动;该结构上的流体微元具有高涡量.由于该结构旋转和高涡量产生了较大的诱导速度,加上该结构的非定常性导致在这些大尺度结构频繁出现的区域具有很高的脉动强度,从而可能带来较高的混合效率.     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

双涡/贫油驻涡燃烧室的贫油熄火特性试验

对在燃烧室进口高速气流条件下工作的贫油双涡结构驻涡燃烧室,通过改变主流马赫数与温度和改变其凹腔几何尺寸进行贫油熄火性能试验.试验结果表明:①随着主流马赫数的增加,贫油熄火余气系数减小,主流速度是影响贫油熄火性能的最主要因素.②主流温度的升高对拓宽燃烧边界有利,但是在达到600K后,对贫油熄火余气系数的影响逐渐减弱.③双涡试验件后体进气量的增加可以使贫油熄火性能变差.④不同凹腔宽深比对贫油熄火余气系数的影响很大,综合各个试验件结果,深宽比为0.88的双涡试验件的贫油熄火性能最好.      

非线性扰动方程研究平面混合层发声机制

在流动过程中由脉动引起的噪声问题在自然界和工业界中广泛存在,研究流动过程中不稳定波的发声机制对于理解、预测并最终控制气动噪声有着重要意义。本文以二维 Blasius 相似性解作为基本流场,以超声速混合层作为研究对象,采用非线性扰动方程(NLDE)研究超声速混合层不稳定波的近场动力学特性与远场声辐射之间的内在联系。针对不同类型的对流马赫数,分别研究对流马赫数 Mc =0.5(<1)和 Mc =1.2(≥1)两种情况下的扰动发声机制,结果表明：当对流马赫数 Mc <1时,流场中的发声机制主要由大尺度结构的涡条发声,且基本波与亚谐波之间的非线性作用能增强辐射强度;当对流马赫数 Mc ≥1时,根据扰动相速度是相对于上层还是下层自由流速度为超声速,可以进一步分为快慢两种模态,分别对快慢两种模态以及其相互作用模态进行了数值模拟研究,计算得到流场中的发声机制是以马赫波形式辐射向远场,即马赫波辐射,其辐角的计算结果与理论值相符。     

超声速粘性剪切流空间稳定性对称紧致差分数值分析

用一类高精度对称紧致差分格式数值离散一阶改型三维可压粘性扰动方程,对导出的非线性离散特征值问题采用二阶修正Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ边值迭代局部解法,实现了超声速剪切流的线性空间稳定性分析。基于对流Ｍ数、Ｒｅ数、速度比、密度比等参数研究,讨论了压缩性效应、粘性效应、超声速扰动快／慢模态等,结果显示超声速粘性剪切流的弱不稳定性同多种影响因素密切相关。     

低排放驻涡燃烧室冷态流场特性试验

利用粒子图像测速仪(PIV)对低排放驻涡燃烧室模型进行冷态流场测量,获得该燃烧室流场的变化规律和压力损失的变化情况.试验结果表明:低排放驻涡燃烧室涡系结构稳定,在主流后方的驻涡区存在主副双涡结构,联焰板后方的驻涡区存在单涡结构,两者之间特征截面呈现出双涡涡系逐渐向单涡涡系过渡,涡心位置也随截面的变换而变化.随着进口马赫数(0.15~0.30)的增大,主涡面积随之增加,副涡面积在进口马赫数为0.2时最大,而各特征截面上主副双涡及单涡的涡心位置基本不变.驻涡区涡系强度及其边缘的气流速度,以及主燃区的气流速度均随进口马赫数增大而提高.总压损失随进口马赫数(0.15~0.30)的增大而增加.      

驻涡加力燃烧室贫油熄火性能的影响

为了研究驻涡加力燃烧室的贫油熄火性能,设计了三种不同方案的凹腔结构,进行了贫油熄火（LBO）性能试验。采用航空煤油作为燃料,试验中的主要研究参数如下:外内涵进气压比在0.97~1.07之间变化;主流马赫数在0.13~0.20之间变化。研究结果表明:驻涡加力燃烧室贫油熄火油气比随着外内涵进气压比的增加而增加,随着主流马赫数的增大而增加,并且在马赫数较大时,变化幅度增大;在驻涡区上部分区域,轴向平均绝对速度越小,以及主涡涡心位置在一定范围内更贴近凹腔后壁面,驻涡加力燃烧室的贫油熄火性能更优。      

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。