空间发展混合流的大涡模拟

采用大涡模拟方法数值模拟了空间发展混合流，通过在平均流中引入小扰动，研究了混合流大涡结构的产生和演化过程，捕捉了展向涡的卷起、配对、合并，以及二次流向涡的出现等大尺度的三维拟序结构。与实验及直接数值模拟的结果相比，吻合程度较好，表明本文所采用的LES方法对湍流的数值模拟是切实可行的。     

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景,采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型,采用高精度WENO格式进行空间离散,并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致,来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构,并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析,揭示了流场结构特征及其时空演化规律,并获得了流场密度脉动特性。      

旋流预混燃烧室流动混合的大涡模拟

采用大涡模拟的方法研究旋流预混燃烧室中的流动及燃料空气混合过程,选取亚格子湍动能模式为基准的涡黏模型,通过模拟涡核的运动过程定性的表征了拟序结构的发展;引入混合不均匀度的概率密度函数,研究入口边界条件尤其是入口速度对于涡团破碎及混合程度的影响规律;通过研究速度分量与燃料质量分数脉动的相关性,分析了不同方向的速度脉动对于混合过程的影响因子.这些均为进一步研究旋流预混燃烧室的燃烧不稳定性及污染物排放提供了理论上的指导.      

折流燃烧室燃烧流场大涡模拟

为深入了解航空发动机折流燃烧室内部复杂流场结构,对一种带有离心甩油盘的单头部环形折流燃烧室冷热态流场进行大涡模拟。数值计算模拟了从启动状态到稳定燃烧状态的完整非稳态过程,获得了该燃烧室流量分配、压力损失等参数以及冷热态流场结构。数值计算结果表明：（1）冷热态下燃烧区流场结构分为主回流区和次回流区两部分,主回流区冷态时呈现多涡结构,热态时回流区形状受燃油射流影响呈现对称的双涡结构;（2）燃烧室中各涡团结构由各进气孔射流相互作用形成,涡团结构促进燃烧室内部的能量和质量交换;（3）热态时燃烧室前后涡流板周围存在两个稳定的点火源。     

脉冲横向射流的大涡模拟

为了研究脉冲气膜冷却的流动特性,采用基于高阶对称加权本质无振荡(WENO)-Roe格式和隐式牛顿迭代时间推进、动态亚格子模型和预处理技术下的大涡模拟(LES)计算方法,对吹气比为1.0下的斜管横向射流进行了稳态和脉冲两种不同工况下的数值模拟.计算结果显示了射流流场中存在发卡涡结构,研究表明射流孔进口位置的脉冲改变了射流出口拟序结构,即改变了射流与主流的掺混过程.      

洁净室流场大涡模拟

应用湍流大涡模拟方法对矢流洁净室和全顶棚送风下部两侧回风式洁净室进行了数值模拟。计算中使用了隐含流线抑风耗散作用的具有三阶部分展开的Taylor-Galerkin离散格式，导出了经Gauss滤波后的基本方程组的有限元列式，应用涡粘性亚格子模式的湍流大涡模拟方法进行数值模拟。数值计算结果与实验结果符合较好，表明湍流大涡模拟能够较准确地捕捉到洁净室流场中旋涡的位置和尺度，在洁净室的数值仿真方面，与标准的K－ε湍流模式相比，湍流大涡模拟结果具有较高的可靠性。     

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

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

涡发生器控制平板边界层分离的大涡模拟

为了研究涡发生器(VGs)间距λ对控制边界层分离效果的影响,选取了4种涡发生器间距,λ/H(H为涡发生器高度)分别为5,7,9,11.采用大涡模拟(LES)方法对带逆压梯度的平板边界层分离流动及VGs控制分离流动进行了数值模拟.分析了有无VGs控制时,湍流场中大尺度相干结构及其演化规律,分别从旋涡间距、边界层内流体动能、压差损失等方面考察了VGs间距对控制流动分离效果的影响.研究结果表明当λ/H为5时,VGs间距过小抑制了旋涡的展向发展,λ/H为9,11时,VGs间距过大边界层内流体动能偏低,当间距λ/H为7时流动控制效果更优,此时计算域压差损失最小,相比较无VGs控制时,压差损失降低了30.95%.      

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.     

基于大涡模拟的圆转矩喷管尾喷流强化掺混机理研究

为揭示2元圆转矩喷管尾喷流强化掺混的内在机制,应用大涡模拟(LES)方程对2种相同进、出口直径的喷管模型(轴对称、2元圆转矩)在Ma=0.8、高雷诺数(2×10~5)条件下进行了数值模拟计算。结果表明:与轴对称喷管相比,圆转矩喷管射流掺混效应增强,速度衰减快,核心区长度和高温区域面积减小。同时尾喷流拟序结构变化说明:2种喷管主要拟序结构均包含涡环、涡辫、发卡涡、螺旋涡等相似结构;但圆转矩喷管在射流近场诱导出的涡旋更丰富,边角剪切涡发展更快,形成明显的CVP结构,导致其射流柱失稳时刻提前、距离缩短;同时,喷管形式的改变使得射流剪切层内雷诺应力增大,速度脉动增强。拟序结构发展及雷诺剪切应力变化说明在射流流场中涡旋发展耗散速度增大、速度边界层脉动增强、射流柱易失稳是导致射流掺混增强的本质因素,为异形喷管的强化掺混机理提供了依据。     

时间发展可压化学反应混合层直接数值模拟

基于一个简化的单步不可逆燃烧动力学模型,对时间发展二维可压化学反应平面混合流场,实现了高精度迎风紧致有限差分直接数值模拟。研究了化学反应放热与流体混合运动间的耦合作用,分析了前者对可压混合层厚度、拟序结构及产物生成率等方面的影响。     

NUMERICAL SIMULATION OF FREE MIXING LAYER WITH CROSS SHEAR

SupportedbytheNationalScienceFoundationofChinaSincethemiddleofthe7O'S,thestudyofevolutionofcoherentstructureinthenominallyplanemixinglayershasoccupiedtheeffortsofmanyresearchesinthefield-Alotofachieve-mentsaboutcounter-rotatingstreamwisevorticeshavebeenobtainedduringthelasttenyears.Mostnaturalshearflowsandmanyindustrialflowsinvolvecomplexities,suchas,withcrossshear.However,thescenarioofthecrossshearintroducedinitiallyisinherentlydifferentfromthatofthecrossshearintroducedinthemiddlestage,ass…     

交叉剪切混合层中拟序结构的数值研究

采用拟谱方法对时间模式的交叉剪切混合层进行了直接数值模拟。计算结果表明:与平面混合层一样,展向涡的拉伸作用是交叉剪切混合层中流向涡形成的主要机制。当展向剪切强度较大时(如两主流交叉角为40°),与初期展向KelvinHelmholtz相关的单向旋转流向涡在拉伸作用下很快增长起来,并“坍缩”成“肋状”涡。当交叉角为40°时,涡核区存在类似平面混合层中“方块状”涡的流向涡结构,展向涡辫区还存在一组符号相反的流向涡,不过与“肋状”涡对应的涡结构呈扁平状,始终没有“坍缩”。当交叉角为60°时,“肋状”涡非常强,以致完全抑制了平面混合层“对称模式”的发展。当交叉角小到20°时,流向涡结构更接近于对称分布,然而“肋状”涡却没有形成。另外,计算结果还证实:与二维混合层相比,大强度展向剪切的引入能够加强流场的混合,同时,适当增加展向扰动波初始强度和波数也是提高混合效率的有效手段。     

平直槽道内湍流运动的大涡模拟

为在中型机上实现对湍流运动的大涡模拟,采用Moin的基本模型,在数值计算中采用了诸如高斯滤波函数、速度压力交错网格、谱方法以及一种简单的压力计算的修正等方法和技巧,计算结果与Moin结果比较,表明在我国现有条件下也有可能进行大涡模拟方面的工作。     

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

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

DNS/LES方法在剪切湍流模拟中的应用

提出采用DNS/LES混合方法进行自由剪切类混合流动湍流研究,该方法兼顾了直接数值模拟(DNS)方法在转捩线性范围内的高精确度模拟和大涡模拟(LES)方法在转捩非线性及全湍流区的高效模拟两方面优点,通过对低对流Mach数流场计算表明该方法可以较好地模拟流场的转捩过程.同时,文中还对比实验分析了湍流场对光学传输效应的影响,给出了与实验一致的结论.     

渐扩通道内湍流流动与混合的数值模拟

应用标准形式的ｋ－ε湍流模型,对具有逆向压力梯度的渐扩通道内两股等温或不等温射流的湍流流动与混合进行了数值模拟。结果表明,时均速度场与温度场的计算结果与实验测量数据相符合。     

钝头双锥喷流致冷流场结构及密度脉动特性

采用大涡模拟方法对钝头双锥喷流致冷流场开展了数值模拟,研究了超声速喷流混合流场结构特征及密度脉动特性。大涡模拟方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runger-Kutta方法。数值模拟清晰地捕捉到了流场波系结构,精细地预测了流动发生失稳、转捩以及发展为充分发展湍流的物理过程,直接获得了流场密度脉动特性。通过有、无喷流状态对称面流场的对比,发现超声速喷流能够有效冷却光学窗口;喷流与主流形成的混合层不稳定,很快发生失稳和转捩,形成大尺度湍流结构,进而引起强烈的密度脉动。此外,获得了钝头双锥整体模型喷流致冷流场的空间发展形态特征。      

MT1涡轮导叶边界层拟序结构分析

采用自主开发的大涡模拟程序NUL-TURBO,对MT1高压涡轮导叶的跨声速流场进行了数值模拟。在对半叶高位置叶片表面等熵马赫数分布实验数据对比验证的基础上,研究了有无气膜冷却两种情况下叶片表面流场的拟序结构。结果表明:无气膜射流时,叶片吸力面近尾缘位置存在分离转捩过程,并在此过程中发现了发卡涡"森林"现象;有气膜射流时,射流出口位置叶片表面边界层直接转捩为湍流,并发现了叶片压力面逆转捩过程中存在的涡拉伸形态。