低压高负荷涡轮叶栅的大涡模拟

了解预测转捩与层流分离泡,对于低压高负荷叶片设计特别重要.为此,采用大涡模拟方法,对典型高负荷叶栅T106进行了模拟.计算描述了层流分离泡,给出了壁面压力、壁面摩擦力系数的分布,从而准确预测了T106压力面上的层流分离、转捩、湍流再附现象.计算同时刻画了层流分离泡的非定常性质,给出了瞬时压力,进行了频谱分析,确定主频为2500 Hz,并观察到了瞬时涡等有意义的现象.计算同时表明二维大涡模拟计算结果可以复现出湍流的许多重要现象并且精度可以满足工程需要.     

低雷诺数下周期性尾迹／层流分离泡相互作用的大涡模拟

本文采用大涡模拟(LES)方法,对低雷诺数条件下周期性尾迹/层流分离泡的相互作用进行了数值模拟.计算域入口尾迹采用独立计算的平面尾迹流给出.对周期性尾迹作用下的层流分离泡进行了统计特性及瞬态流场分析,并与相同条件下的无尾迹流动进行比较.计算结果表明:在尾迹作用下,时均层流分离泡长度明显缩短,分离泡变小;相位平均的分离点/再附点随时间不断变化;从瞬态结果容易看出,尾迹诱导的流向涡结构是导致分离泡提前转捩的主要原因.     

3维Mortar谱元法模拟Kelvin-Helmhtz不稳定性混合层

采用Mortar谱元法和多处理器并行计算技术模拟了Kelvin-Helmhtz界面不稳定性湍流的混合发展过程,通过对混合层动量厚度、能谱和总动能的计算,评估了Kelvin-Helmhtz混合层的演化机理。计算结果表明:3维Mortar谱元法具有高计算精度和光滑区域的指数收敛特性,可以有效模拟混合层流动的湍流混合和演化,能够捕捉到涡的合并现象和大涡到小涡的级联过程;初期的混合层层流运动发展成具有连续谱结构的湍流运动过程,实现了Kelvin-Helmhtz界面不稳定性混合层流动从2维发展到3维的转捩特征,总湍流统计动能的变化反映了粘性耗散过程的作用。通过对Kelvin-Helmhtz 3维界面不稳定性混合层流动和3维层流向湍流转捩过程的数值模拟,程序的有效性得到了验算,表明谱元法应用于湍流混合模拟是可行的。     

壁面换热和压力梯度条件下边界层演化的LES研究

本文采用大涡模拟(LES)方法对壁面换热边界条件下有压力梯度的分离边界层流动问题进行了数值模拟研究,对不同壁面热边界条件下分离边界层进行了统计特性及瞬态流场分析.结果表明:与绝热壁面条件相比,等温冷壁面条件(0.8倍来流总温)下附着边界层内的速度剖面比绝热壁面更加饱满,边界层形状因子降低、分离泡也被显著抑制,其时均分离泡尺寸无论是在流向还是法向均更小;等温冷壁面条件下分离剪切层展向涡脱落的位置提前,展向涡卷起的频率降低,分离剪切层内相干结构尺度整体减小,促使分离边界层转捩加速和再附位置前移。     

低速轴流涡轮内部流动的大涡模拟

在低进口雷诺数下,低速轴流涡轮内部可能存在复杂的边界层转捩和分离流动。准确模拟边界层转捩和流动分离对低速轴流涡轮的气动设计具有重要意义。本文以某单级低速轴流涡轮为研究对象,采用大涡模拟方法对其在进口雷诺数为20000情况下的内部流动进行了数值模拟研究,并与前期采用全层流模型、S-A模型、Abu-GhannamShaw(AGS)转捩模型的模拟结果进行了对比,对比分析发现,大涡模拟结果与实验结果吻合更好,可以准确模拟该涡轮叶片吸力面的流动分离和叶片通道内的二次流动。由大涡模拟结果可以得出,静叶尾迹和分离使尾迹区的流体流动速度降低,但尾迹对流动角的影响较小。动叶入口低速微团在做周向运动的同时沿径向运动;高速微团主要沿周向运动。静叶叶片表面的分离流存在较大的由叶顶向叶根的径向的运动;动叶吸力面叶顶处也存在较大的分离流动.     

基于大涡模拟下风力机翼型非定常转捩流动的动态模态分解分析

翼型绕翼流动对风力机整机性能产生重要影响.本文基于大涡模拟方法,得到风力机翼型非定常转捩流动.通过对压力流场的动态模态分解(DMD)分析,发现翼型层流分离泡的生成和发展为流动主要非定常特征,且该特征具有主要频率.预估得到高增长率T-S扰动波频率与DMD模态频率接近,发现高频变化的分离泡由边界层分离点下游不远处的T-S扰动波诱导K-H不稳定性而主导,以及中低频的DMD模态表现为T-S扰动波所引起的湍流结构.对比不同攻角下的DMD分解结果,发现上表面分离泡会逐渐向前移动,长度变短;而下表面分离泡略微后移,长度变化不大.     

B-C转捩模型的DDES方法模拟转捩大分离流动

根据零压力梯度平板的实验数据对B-C转捩模型中的经验关联式进行标定.标定后的B-C转捩模型在进口湍流度较低的平板算例上可以提供更合理的转捩预测结果.但受限于RANS方法自身特性，B-C转捩模型对转捩后出现的大分离流动计算可靠性较差.结合延迟分离涡模拟SA-DDES方法和B-C转捩模型，发展了转捩-延迟分离涡模拟BC-DDES方法.三维S-K平板的数值模拟表明，BC-DDES方法能得到与B-C转捩模型一致的转捩结果.三维圆柱绕流的数值模拟表明，BC-DDES方法得到与实验符合的平均阻力系数和表面压力分布，并且计算花费比tHRLES方法少.     

低Reynolds数翼型绕流主动控制技术

在低Reynolds数条件下,翼型绕流的上表面边界层由于抗逆压梯度能力变差容易发生流动分离,从而形成长层流分离泡.分离泡通常是非定常的,会诱发边界层的转捩、再附并形成湍流边界层.这个过程会使翼型的气动性能急剧下降,并伴随着强非线性效应.转捩后形成的湍流边界层也会产生高摩擦阻力.针对这种现象,文章以NACA0012翼型为例,通过隐式大涡模拟研究了有效的主动控制方案.为了统一分离控制技术和湍流边界层减阻技术,研究了在平板或槽道湍流中取得较好控制效果的壁面垂向反向控制方案.首先利用隐式大涡模拟研究了低Reynolds数条件下NACA0012翼型绕流的流场特征.其次分析并验证了反向控制方案在分离区控制流场的可行性,发现反向控制在分离区的作用相当于基于流场信息的壁面抽吸控制,且控制具有实时性和高效性,控制抽吸了前缘的低能流体,使得翼型前缘附面层变薄,并增强了其抗逆压梯度的能力,较大程度提高了翼型的气动性能.最后在湍流边界层验证了其减阻控制效果,发现反向控制阻断了流向涡的法向输运,抑制了涡结构的发展,并减弱了猝发过程,使得湍流的高摩阻力得到了有效降低.      

压缩拐角激波/湍流边界层干扰特性分析

用大涡模拟方法对Mach数3.0下的压缩拐角激波/湍流边界层干扰问题进行数值研究.对拐角上游平板区域边界层转捩及湍流进行模拟,设定平板区域长度,使得转捩过程于平板区域发生并充分完成,从而在拐角处产生激波/湍流边界层相互干扰,研究激波/湍流边界层的作用机理.研究表明:流场能够在非定常扰动激励下迅速转捩,并于平板区发展为完全湍流;湍流边界层与激波相互作用过程中,拐角附近分离区较层流情况明显减小;展向不同区域分离区大小差异较大,局部区域分离现象消失.     

飞行器低Reynolds数层流分离理论探讨

低Reynolds数流动由于自身特点导致气动特性严重恶化,非定常、非线性效应突出且预测困难,加之相关基础理论研究不足,给以临近空间低速飞行器和高性能微小型飞行器为代表的低Reynolds数飞行器的开发和研制带来了瓶颈和挑战.首先概述了飞行器低Reynolds数的范畴、低Reynolds数空气动力学的主要问题与挑战.随后从低Reynolds数层流分离基础理论出发,依次介绍了低Reynolds数层流分离经典理论、低Reynolds数层流分离非定常流动特性、低Reynolds数后缘层流分离泡.在此基础上,通过对经典长层流分离泡与后缘层流分离泡力学特性的差异以及随攻角和Reynolds数的演化规律的详细分析,逐步揭示了一些低Reynolds数复杂气动效应的本质,如小攻角升力系数的非线性效应,翼型随Reynolds数下降气动特性的二次恶化效应等.最后对低Reynolds数流动基础理论的发展过程进行了总结,并对层流分离诱导转捩及再附效应等复杂流动问题进行了展望.      

气泡-液体闭式多股射流的大涡模拟研究

采用气泡-液体两相流动的欧拉-拉氏大涡模拟,研究了矩形通道内多股射流形成的气泡-液体两相湍流流动,得到了气液两相湍流瞬态结构,产生和发展过程。研究结果发现气泡和液体都有瞬态大涡结构,气泡脉动比液体的强。大涡模拟统计结果给出了有无气泡两种情况下的液体湍流脉动速度均方根值分布。瞬态和统计结果都表明,气泡增大了液体湍流,液体湍流来源于其自身的剪切产生和气泡的作用。这与二阶矩模型模拟结果定性一致。     

非定常瞬态流动过程中的Lagrangian拟序结构与物质输运作用

从Lagrangian角度数值分析了圆柱瞬时起动过程中的非定常瞬态流动现象,如分离泡产生、破裂和涡脱落等及其产生的非定常效应,揭示了所列现象诱导的物质输运和迁移效应,首先采用双时间步长的特征线算子分裂算法数值模拟了圆柱起动过程中的瞬时流场,然后采用数值方法从流场中提取出Lagrangian拟序结构(LCSs),并根据非线性动力学理论研究了流动分离和旋涡演化过程中的物质输运作用,结果表明,圆柱瞬时起动后所产生的非定常阻力与相应瞬态现象中的物质输运有密切的关系:对称分离泡产生及其在流向方向的生长,能够使分离泡内压力升高且分布均匀,从而减小阻力;对称分离泡的失稳增强了分离泡与主流之间的物质输运作用,最终导致涡的脱落,并有利于推迟流动分离和减小分离区域,非定常流动中LCSs所描述的物质输运和迁移作用对流动控制的机理研究具有一定指导意义。     

Large eddy simulation of channel flow with and without periodic constrictions using the explicit algebraic subgrid-scale model

We analyse the performance of the explicit algebraic subgrid-scale (SGS) stress model (EASSM) in large eddy simulation (LES) of plane channel flow and the flow in a channel with streamwise periodic hill-shaped constrictions (periodic hill flow) which induce separation. The LESs are performed with the Code_Saturne which is an unstructured collocated finite volume solver with a second-order spatial discretisation suitable for LES of incompressible flow in complex geometries. At first, performance of the EASSM in LES of plane channel flow at two different resolutions using the Code_Saturne and a pseudo-spectral method is analysed. It is observed that the EASSM predictions of the mean velocity and Reynolds stresses are more accurate than the conventional dynamic Smagorinsky model (DSM). The results with the pseudo-spectral method were, in general, more accurate. In the second step, LES with the EASSM of flow separation in the periodic hill flow is compared to LES with the DSM, no SGS model and a highly resolved LES data using the DSM. Results show that the mean velocity profiles, the friction and pressure coefficients, the length and shape of the recirculation bubble, as well as the Reynolds stresses are considerably better predicted by the EASSM than the DSM and the no SGS model simulations. It was also observed that in some parts of the domain, the resolved strain-rate and SGS shear stress have the same sign. The DSM cannot produce a correct SGS stress in this case, in contrast to the EASSM.     

Large eddy simulation of compressible turbulent channel flow with spanwise wall oscillation

The influences of the modification of turbulent coherent structures on temperature field and heat transfer in turbulent channel flow are studied using large eddy simulation (LES) of compressible turbulent channel flows with spanwise wall oscillation (SWO). The reliability of the LES on such problems is proved by the comparisons of the drag reduction data with those of other researches. The high consistency of coherent velocity structures and temperature structures is found based on the analyses of the turbulent flow field. When the coherent velocity structures are suppressed, the transportations of momentum and heat are reduced simultaneously, demonstrating the same trend. This shows that the turbulent coherent structures have the same effects on the transportations of momentum and heat. The averaged wall heat flux can be reduced with appropriate oscillating parameters. Supported by the Key Subjects of National Natural Science Foundation of China (Grant No. 10732090), the National Natural Science Foundation of China (Grant No. 50476004), and the 111 Project (Grant No. B08009)     

Snapshot POD analysis of transient flow in the pilot stage of a jet pipe servo valve

The spatial evolution of a turbulent flow in the pilot stage of a jet pipe servo valve at the inlet pressure and deflection angle of the jet pipe is investigated using a large eddy simulation (LES). The pressure of the same flow field is measured by a high frequency dynamic pressure sensor in the experiments and is compared with the results of the LES, as well as their root-mean-square (RMS) and fast Fourier transform (FFT) results. The results of experiments and LES are in good agreement, indicating that LES is able to predict the flow dynamics. Velocity datasets based on LES are utilised to conduct the snapshot proper orthogonal decomposition (snapshot POD) technique. The snapshot POD analysis results of the first 4 modes show a full ability to directly visualise details of the coherent structures. The influences of the inlet pressure and deflection angle of the jet pipe are also discussed. Under different inlet pressures, the velocity eigenfunctions of the first mode are similar, while the locations and strengths of the vortices in high modes are different. The Lamb-Oseen vortices that affect the trajectory of jet streams are observed in the vicinity of the entrances of receiver channels only in the first mode, and several spindly vortices appear in the region of ?5?y/n?相似文献   

射流温度场与拟序结构间相互作用的大涡模拟

本文分别对常温空气和高温空气的平面自由射流进行了大涡模拟,采用分步投影法求解动量方程,亚格子项采用标准Ｓｍａｇｏｒｉｎｓｋｙ亚格子模式模拟,并同时求解了标志物输运方程以实现数值流场显示。模拟结果给出了湍流的瞬态发展过程以及流动中拟序结构的发展演变过程,通过对常温射流和高温射流瞬态流场的比较揭示了射流中温度场与拟序结构间相互作用的细节过程。     

Study of interactions between liquid vortices and bubbles by simultaneous particle image velocimetry and laser-induced fluorescence tracking

Extensive experimental research has been conducted using the particle image velocimetry (PIV), laser-induced fluorescence (LIF) imaging and backlit photographic recordings to study the complex interactions between coherent vortex structures created in the shear layer of jets and the bubbles. Triggering of the naturally-developing instabilities of the shear layer by a thin, pulsed annular flow surrounding the jets allowed the creation of large, orderly structures with controllable frequency and phase. Synchronization of the triggering with data acquisition permitted phase averaging of the data and revealed several interesting phenomena. In particular, the evolution of large vortices and bubble fields could be tracked and the interactions could be studied. The horizontal and vertical velocity components of the liquid and bubble fields and the vertical velocity of both the vortex and bubble rings that were created were measured by the PIV. LIF and image recordings have been combined to visualize bubble trapping inside large eddy structures.     

旋流和无旋突扩流动的LES和RANS模拟

本文用smagorinsky-Lilly亚网格尺度湍流模型对旋流突扩流动(s=0.53)和无旋突扩流动(s=0)进行了大涡模拟(LES模拟),同时分别用压力应变项为IPCM和IPCM+Wall模型的雷诺应力方程模型进行了RANS模拟,和LES的统计结果对比。LES的统计结果与雷诺应力模型的模拟结果及实验对照表明,LES结果与实验结果的吻合比雷诺应力模型的好,说明所用的亚网格尺度湍流模型对旋流流动是适用的,LES结果是可信的。LES的瞬态结果揭示出在旋流作用下,流场中存在复杂的旋涡脱落现象。大涡结构极易破碎成小涡,而在无旋突扩流动的情况下,由于剪切的作用更强,大涡结构的尺寸和范围比旋流流动的要大。     

Surface pressure characteristics of a highly loaded turbine blade at design and off-design conditions; a CFD methodology

The flow field passing through a highly loaded low pressure (LP) turbine cascade is numerically investigated at design and off-design conditions. The Field Operation And Manipulation (OpenFOAM) platform is used as the computational Fluid Dynamics (CFD) tool. In this regard, the influences of grid resolution on the results of k-ε, k-ω, and large-eddy simulation (LES) turbulence models are investigated and compared with those of experimental measurements. A numerical pressure undershoot is appeared near the end of blade pressure surface which is sensitive to grid resolution and flow turbulence modeling. The LES model is able to resolve separation on both coarse and fine grid resolutions. In addition, the off-design flow condition is modeled by negative and positive inflow incidence angles. The numerical experiments show that a separation bubble generated on blade pressure side is predicted by LES. The total pressure drop has also been calculated at incidence angles between -20° and +8°. The minimum total pressure drop is obtained by k-ω and LES at design point.