湍流边界层各向异性张量模式的实验研究

使用高时间分辨率粒子图像测速技术在水槽内精细测量平板湍流边界层内不同流向、法向位置的瞬时速度矢量空间场的时间序列信号。利用实验数据对湍流各向异性的涡黏张量模式进行了实验研究,给出了湍流边界层不可压缩槽道湍流中各向异性的涡黏张量各分量的空间分布规律。各向异性的涡黏张量分量的空间分布揭示了流向平均速度的法向梯度对雷诺应力三个分量的不同贡献。涡黏张量三个分量空间上的不同规律说明对涡黏张量的各向异性的考虑是合理的,而各向同性的涡黏模型不适宜描述剪切湍流雷诺应力的物理本质。     

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

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

吹吸扰动对壁湍流边界层摩擦阻力的影响

对风洞中平板湍流边界层施加局部周期性吹吸扰动,并用IFA-300热线风速仪测量了其近壁区域不同法向位置瞬时速度的时间序列信号.利用湍流边界层对数律平均速度剖面,计算湍流的壁面摩擦阻力.对平板湍流脉动速度信号用子波分析进行多尺度分解,通过条件采样以及相位平均的方法获得不同尺度下的相干结构波形和能量分布,来研究吹吸扰动对平板湍流边界层的摩擦阻力的影响.结果表明,扰动使得窄缝附近流动阻力增加,而下游区域流动阻力减小.      

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

GAO-YONG湍流模型对湍流传热的数值模拟

将GAO-YONG湍流模型应用于湍流传热的研究,分别计算了平板剪切湍流和二维平面冲击射流的湍流传热问题.边界层剪切湍流流动与换热的计算表明:与传统的湍流模型不同,GAO-YONG湍流模型不需要对近壁区域做任何特殊处理(比如壁面函数、低Reynolds数修正等)即可模拟出从壁面到主流区的全部流动与传热情况;另外,对于冲击射流Nusselt数的模拟也得到了与实验符合较好的计算结果,准确地捕捉到了2种冲击高度下流场换热的不同特征,表明了GAO-YONG湍流模型能够较高精度地计算湍流换热.      

边界层中湍动能和耗散能最大的尺度分量特征研究

阐述了利用离散正交小波理论分析 1 -D湍流测试信号的方法,使用热线风速仪,在 3个雷诺数下,分别测量了光滑壁面湍流边界层的速度和加速度信号。运用所述方法分析了边界层测试信号,给出了最大湍动能密度和最大耗散能密度及其对应尺度沿 y+ 的分布曲线,并且给出该尺度分量的发生频率沿 y+ 分布曲线。分析可见 :两种最大能量密度的无量纲尺度沿 y+分布具有雷诺相似性,其发生频率随着雷诺数增加而增加,但是最大湍动能密度尺度的频率沿 y+基本是一常数;无量纲的最大湍动能密度沿 y+在边界层小部分区域具有雷诺相似性,而无量纲的最大耗散能密度沿 y+ 在整个边界层内具有雷诺相似性。     

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

SST湍流模型的激波与可压缩效应改进

由于缺乏对某些重要流动特征的考虑,针对不可压流发展的标准SST湍流模型在描述超声速流场时存在明显的局限性。为改善SST模型在吸气式高超声速推进系统内流等复杂超声速流场中的预测精度,基于流动特征结构定向开展了激波和可压缩效应改进。通过激波/湍流边界层判别函数和可压缩效应判别函数定位标准SST模型参数或建模假设失效的区域,针对性地改进湍流模型。采用超声速平板边界层流动、超声速压缩拐角分离流动、超声速隔离段复杂激波串流动以及HIFiRE-2超声速内流等算例进行了测试,结果表明改进模型具有与标准SST模型一致的边界层预测能力,但显著提高了对激波干扰流动及逆压分离流的预测精度。     

温度边界层对壁湍流多尺度相干结构的影响

用热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中壁面加热平板湍流边界层不同法向位置流向速度分量的时间序列信号, 与常温情况相比, 热对流加速了壁湍流中流体的动量、质量和能量的交换, 加快了湍流边界层的发展;用子波分析方法研究了壁面加热对壁湍流多尺度相干结构喷射和扫掠条件相位平均波形的影响, 显示热对流明显增大了缓冲层相干结构扫掠过程和喷射过程的强度, 是增强壁湍流中流体的动量、质量和能量交换的根本原因.      

基于变分模态分解的壁湍流拟序结构形态研究

应用多分量准二维变分模态分解技术,对直接数值模拟的湍流边界层的流向-法向平面进行了分解。通过对从低雷诺数到中等雷诺数(Reτ=400～1750)各个工况的瞬时流场的观察,定性地得到了不同模态所对应的物理结构特征。通过两点相关云图形态和倾斜角度的比较,定量地对比了不同模态对应的拟序结构倾角随法向高度的变化。一阶模态主要对应大尺度扫掠/喷射事件,第二阶模态则对应近壁流向涡结构。两阶模态的拟序结构倾角随法向高度的变化趋势与雷诺数呈弱相关,表明大小尺度模态的几何特征在所研究的雷诺数范围内存在自相似性。     

模型的转捩流动计算分析简

 为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Re_θt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序（HUNS3D）。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Re_θt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing 高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明：γ-Re_θt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

基于大涡模拟分析气膜冷却的湍流场

用大涡模拟考察了单孔平板气膜冷却的湍流场,气膜孔沿流向倾斜30°,气膜出流的雷诺数为2600,吹风比为0.5.计算结果表明:①气膜孔内的流动分离增加了湍流动能,导致气膜冷却效率的降低,在设计中要尽量避免或减少流动分离;②与射流侧向扩展有关的涡黏性系数在气膜孔两侧存在峰值,而采用各向同性的湍流模型预测气膜冷却时,涡黏性系数的峰值出现在射流与主流的剪切区,因此会高估射流的垂直穿透而低估射流的侧向扩展;③可以用大涡模拟辅助建立各向异性的湍流模型,以便提高湍流模型的预测精度.      

一种近壁理论湍流模型初探

在传统的k-ε湍流模型基础上,在近壁区放弃求解ε方程而代之以根据共振三波理论模型在近壁区系统计算的ε值构成封闭,通过对NACA0012二维翼形流场的计算,初步验证并适当修正了新的理论模型,由于注意到了理论模型计算结果的无量纲性,从而改善了该流场的计算结果,得到了与实验更为相符的C_p分布。本文的分析还表明,逆压力梯度对ε的影响非常显著,验证了理论分析和计算得到的结论。      

Turbulent transition in solar surges

It has been suggested that a surge can be modelled as a jet travelling in a sheared magnetic field, and that the transition to turbulence of this MHD tearing jet can explain several key observed features. In this paper we present our preliminary results of the transition to turbulencevia secondary instabilities of the MHD tearing jet. Our results confirm that turbulent transition can decelerate the surge, with decay times which compare well with surge data. Furthermore, we find that the turbulent MHD tearing jet forms magnetic field-aligned velocity filaments similar to those often observed in the surge flow field.     

可压缩湍流的多尺度分析

作者的研究团队近几年在可压缩湍流的多尺度性质方面开展了系统的研究工作.通过多过程分解,研究了可压缩湍流中的速度和热力学量的剪切部分、胀压部分、伪声模态、声模态、熵模态的多尺度性质,并总结了各类可压缩条件下的速度和热力学量的谱的标度律.通过滤波方法,研究了动能和热力学量的多尺度传输现象,并重点分析了可压缩性对动能胀压部分...     

大气紊流频谱与飞机响应的几个问题

 本文推导了各向同性紊流速度及其任意阶梯度间的相关函数和频谱函数的一般表达式,并分析了频谱函数对空间频率的奇偶性。对现有文献和飞行品质规范中几个含糊不清的问题进行了阐述和澄情。对飞机纵向大气紊流响应进行了研究计算,得到了一些有意义的结论。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

基于湍流输运特性对S-A模型在压气机角区流动模拟中的改进研究

在非设计工况下,压气机叶栅角区的三维分离加剧,流场内的湍流处于较强的非平衡状态,因此必须对于非平衡状态的湍流的输运过程给出准确的预估.基于不同工况下压气机内角区分离流动的湍流输运过程的分析,提出了一种修正Spalart-Allmaras(S-A)模型中模拟湍流黏性生成和耗散关系的改进方法,经过数值研究发现改进的S-A模型对于静子角区分离的模拟更接近物理实际,从而提高了对非设计工况的预估能力.      

改进S-A湍流模型对横向射流的CFD模拟

通过在Spalart-Allmaras(S-A)湍流模型中加入应变率张量,来改善模型对横向射流(JICF)的计算精度。模拟了速度比分别为0.5和1.5工况下的横向射流。经实验值验证,结果表明:改进后的S-A模型预测精度有所提升;并对湍流模型改进部分对横向射流的作用机制做出了分析。将改进后S-A模型的计算结果与k-ε模型,雷诺应力模型(RSM)和大涡模拟(LES)的结果进行比较,发现S-A模型在速度比为0.5时精度最高,而在速度比为1.5工况下没有明显优势;并对这一现象做出了合理分析。      

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.