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

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

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

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

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

二维可压缩自由剪切层的理论分析

本文从"二维可压缩自由剪切层存在相对稳定的大尺度拟序涡结构,流场应该满足什么样的条件"这样一个反问题出发,应用常微分方程奇点理论对NS方程进行了分析,维持大尺度拟序涡结构的条件是流场内有交替出现的能量不平衡点,鞍点是温度极大值点,结点是温度极小值点;压缩性增加抑制了结点的生成和持续,是导致可压缩自由剪切层厚度增长速度比不可压情况变小的机理;得到的结论可以合理解释现有的物理现象.其次,对于目前研究可压缩自由剪切层广泛应用的对流马赫数进行讨论,认为在建立过程中采用鞍点是压力极大值点假设条件不合理,得出的公式不能反映介质特性,本文采用鞍点是温度极大值点作为假设条件,提出新的求解对流马赫数的公式.根据流动机理分析,提出通过在流场内建立交替出现的鞍点和结点,触发大尺度拟序涡结构,解决超燃冲压发动机混合问题;理论上给出大尺度拟序涡结构频率条件.     

可压混合层涡并放热效应的涡动力分析

 基于双元单步不可逆燃烧动力学简化模型 ,采用紧致差分算法直接求解二维可压 Navier-Stokes方程 ,研究了时间发展非预混扩散火焰平面反应混合层流动。拟序结构由数值线性稳定性分析方法设定的粘性可压最不稳定小扰动初场促生。随放热强度增大 ,得到了基频涡卷撕裂现象。涡动力分析表明 ,放热主要通过膨胀和斜压这两种强非线性过程抑制大尺度剪切涡的演化     

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

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

壁湍流调制雷诺应力和拟序结构应变率空间形态的层析PIV实验研究

采用高时间分辨率的层析PIV技术,测量了水洞壁湍流三维速度分量空间分布的时间序列,应用局部平均速度结构函数概念对流向速度信号进行多尺度分解,以流向速度分量的局部平均速度结构函数过零点作为特征量检测壁湍流中拟序结构猝发的喷射和扫掠过程,应用空间相位平均技术提取拟序结构猝发的喷射和扫掠过程各速度分量、涡量分量、拟序结构速度应变率分量以及调制雷诺应力分量的空间相位平均拓扑形态。为了引入平衡态涡粘模型假设模拟调制雷诺应力,研究了拟序结构猝发过程中调制雷诺应力分量和速度变形率分量的空间分布形态,发现两者之间的空间相位分布不一致。由于存在时空相位不同步性,说明需要考虑大尺度拟序涡结构引起动量传递的时空弛豫效应。应用经典的线性平衡态下的Boussinesq涡粘模型不能准确地描述壁湍流拟序结构动量传递非平衡现象的物理机理。对于壁湍流拟序结构动力学方程中调制雷诺应力的模拟,应采用包含时空相位信息的复涡粘张量模型。由于雷诺应力与速度变形率的时空不同步性,对非平衡非局部湍流场的数值模拟提出了一个挑战性的问题,建议采用包含相位信息的复涡粘张量模型来模拟雷诺应力张量,从而更加符合雷诺应力演化的物理机理,这一模型有可能成为一个很有发展前景的封闭模型,从而更加准确地预测工业领域中广泛存在的非平衡湍流。     

混合层中由次谐波导致的失稳过程的数值研究

用拟谱方法对随时间发展的二维混合层流场进行了直接数值模拟,给出了高雷诺数下在Ｋｅｌｖｉｎ－Ｈｅｌｍｈｏｌｔｚ波卷起后存在的限振幅次振波地大尺度拟序结构进一步的失稳演变过程,包括大涡的配对,撕裂以及二次配对,同时探讨了初始的次谐波强度,雷诺数相位差等对大涡配对和合并过程的影响。     

利用PIV技术对非光滑表面湍流边界层的实验研究

利用在线式 Ｐ Ｉ Ｖ 系统在低速风洞中对两种非光滑表面：阵列涡发生器表面和波纹壁面的湍流边界层进行了实验测量。观察到了壁面几何形状的改变对非光滑表面湍流边界层拟序结构的产生和发展的影响：阵列涡发生器表面（１０ｍ ／ｓ）湍流边界层内有明显的双剪切带状结构,外剪切带状结构接近边界层的外边界,小尺度的涡在内剪切带状结构的附近产生;波纹壁面（２０ｍ ／ｓ）湍流边界层内涡的尺度比较小。并在相同的壁面几何形状条件下,在不同的流动工况下,研究了非光滑表面对湍流边界层拟序结构的影响。实验结果表明,壁面几何形状的改变对外层的大尺度横向涡的产生和发展有明显的影响;而这种影响效果在不同的流动工况下相差很大。     

带凹坑圆孔气膜冷却流动传热特性

利用大涡模拟方法对平板带凹坑气膜孔冷气射流与主流的相互作用展开研究,并结合红外热成像壁温测试和粒子图像测速结果进行实验验证,揭示了凹坑对气膜冷却的强化机理.结果表明:凹坑在气膜孔下游诱发反肾型涡对,将冷气向壁面卷吸,改善了气膜贴壁性;与圆孔气膜冷却相似,凹坑气膜冷却流场也出现马蹄涡、发卡涡等大尺度拟序结构,但与圆孔下游...     

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

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

轴向旋涡强化矩形喷流掺混的机理研究

本文介绍了在低速射流风洞上开展的二维喷管来流带旋涡的模拟实验。通过流场观察，壁面静压和矩形喷流三维流场的测试。     

超声速混合层转捩过程实验研究的进展

超声速湍流机理的实验研究是一件十分困难的工作.在2000年以来,本研究小组在低噪声超声速混合层风洞研究、超声速流动精细结构测量技术研究方面取得了重要进展,这给超声速混合层湍流精细结构的研究奠定了基础.为了研究超声速混合层及其气动光学问题,在研制的超声速混合层风洞中,主要以基于纳米技术的平面激光散射技术(Nano-trace Planar Laser Scattering,简称NPLS)为基础,研究了几种对流马赫数的超声速混合层从层流到湍流转捩过程K-H不稳定涡的空间结构,以及K-H不稳定涡的空间结构随着时间的发展过程.实验结果清晰地反映了湍流混合的不稳定性与转捩的精细结构,以及转捩过程的展向精细结构.     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

大负荷低压涡轮叶型分离转捩流动的大涡模拟

应用动力模式大涡模拟数值方法,对来流无扰动、不可压、雷诺数为5×104(基于进口速度和轴向弦长),定常来流条件下大负荷低压涡轮叶型(Pak B)叶型吸力面非定常分离转捩流动进行了三维数值模拟。在与相关实验数据的对比基础之上,对非定常流动物理信息进行了详细的分析讨论,揭示了计算来流状态下的Pak B叶型吸力面非定常分离转捩流动机理。结果表明,由无粘Kelvin-Helmholtz机制产生的、空间线性增长的初始二维不稳定性在分离剪切内诱导展向旋涡形成并脱落,脱落过程中的展向涡在非线性增长的三维不稳定性作用下发生变形并最终破碎成湍流。计算得到的Kelvin-Helmholtz不稳定性特征频率处于相关实验测量范围内。     

V型尾缘喷管流场特性数值模拟

模拟计算了一系列V型尾缘(Chevron)收敛喷管的详细涡结构(包括流向涡和展向涡)、引射系数、推力损失等性能,并与相关实验结果进行了对比验证。采用有限体积QUC IK格式和标准k-ε湍流模型,结合壁面函数求解雷诺平均的N-S方程。通过模拟计算与实验结果的对比,验证了计算方法的可行性。通过计算得出:V型尾缘喷管齿弯角越大,核心区长度越短,喷管出口下游的流向涡和展向涡的强度和尺度也越大。V型尾缘的齿弯角每增加1°,对应的喷管的堵塞比约增加0.98%,喷管的推力损失约增加1.71%,引射系数约增加3.45%。     

受限空间二维超声速混合层混合增强数值研究

本文研究了具有双燃烧室冲压发动机超声速燃烧室工程背景的受限空间超声速混合层的混合增强措施——低速流入口加入沿流向振荡。研究结果表明：它对二维混合层的混合增强效果明显,加振荡后流场燃烧强度加大,但仍属于较弱燃烧强度;振荡频率决定着涡量聚集区域的位置和大小,振幅决定着涡量分配和聚集的均匀程度。     

双燃烧室冲压发动机超声速混合层混合增强数值研究

研究了具有双燃烧室冲压发动机超声速燃烧室工程背景的特定超声速混合层的混合增强措施——低速流入口加入沿流向振荡.研究结果表明:它对二维混合层的混合增强效果明显,加振荡后流场燃烧强度加大,但仍属于较弱燃烧强度;振荡频率决定着涡量聚集区域的位置和大小,振幅决定着涡量分配和聚集的均匀程度.      

波瓣混合器涡系结构及射流掺混机理的数值研究

借助流体力学软件ANSYS CFX,对波瓣混合器射流掺混流场进行了全三维定常数值模拟,研究了流场中各涡系结构的形成机理及发展过程,并详细探讨了其加速射流掺混过程的作用机制.结果表明:基于SST(shear stress transport)模型的封闭N-S方程能较好地模拟波瓣混合器射流掺混过程,波瓣特殊几何外形诱导产生的流向涡主要通过扭曲内外涵交界面的间接方式加速射流掺混过程,波瓣下游剪切层中K-H(Kelvin-Helmholtz)不稳定性发展而成的正交涡是直接加速射流掺混的关键因素,波谷附近二次流之间的相互作用所产生的通道涡对该区域内的射流掺混有明显的加速作用,受波瓣前缘切割的边界层在径向压力梯度作用下沿波瓣表面卷起而形成的马蹄涡对射流掺混的影响不是特别明显.