应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

基于大涡模拟的压缩拐角激波/边界层干扰研究

为了提高对激波/边界层相互干扰的基本理解,采用大涡模拟(LES)对来流马赫数Ma=2.9,转角为24°的压缩拐角激波与湍流边界层相互干扰进行了研究。采用回收/调节方法作为入口湍流生成技术,并在超声速平板湍流边界层上进行了验证。采用涡识别方法和数值纹影图等流动显示方法,研究了干扰区内激波与边界层相互干扰的结构变化特征。通过对比分析湍动能和雷诺正应力在不同流向位置的分布规律,研究表明:经过激波干扰后湍动能主要集中在边界层的外层,并在拐角附近形成一个低湍动能区;雷诺正应力流向分量和法向分量在边界层内的分布呈现为单峰结构,而展向分量呈现为双峰结构。运用间歇因子对分离激波的大尺度流向运动进行研究,发现激波围绕着平均分离点作前后运动,运动的尺度等于进口湍流边界层厚度的72%。证实了拐角下游G?rtler流向涡对的存在,并对其展向分布和空间演化特性进行了详细研究。     

沟槽表面边界层湍动能分布规律

 通过对不同风速下不同尺寸V型沟槽表面及光洁平板表面边界层内湍动能的测试,对比分析了沟槽表面边界层湍动能的分布规律。试验在一小型专用风洞中开展,流场测试中使用恒温式IFA300智能型流动分析仪,测试模型则采用有机玻璃材质的矩形平板结构;而在沟槽表面理论零点及壁面摩擦速度的计算中,采用基于湍流边界层Spalding壁面律公式同时计算壁面理论零点和壁面摩擦速度的改进方法。最终研究结果表明,沟槽结构主要影响边界层流场的近壁区,沟槽的存在提高了边界层中黏性底层内湍流脉动所具有的动能,有效降低了边界层中过渡区内的湍动能,最大相对降幅超过10%,较好地验证了基于“二次涡”的沟槽表面减阻理论。     

流线正交坐标下湍流分离再附流的实验研究

本文采用二维激光多普勒测速仪测量了二维非对称曲壁扩压器内的湍流分离再附流动。湍流边界层在强逆压梯度下于曲壁上分离,而后于后续的平直通道上再附。扩压器进口处的雷诺数为1.2×10~5,速度为25.2米/秒。在新定义的局部斜流线坐标下,给出了实验结果,并做了分析。实验结果表明:在瞬时分离点后的近壁反流中,雷诺剪应力为负值。采用最大雷诺剪应力为尺度,从分离至再附的流动过程,存在着雷诺剪应力的相似分布,同时Schofield-Perry速度分布在正向流动区成立。由雷诺剪应力相似分布及Schofield-Perry速度分布可假设一新的涡粘性湍流模型,其长度尺度来源于Schofield-Perry速度分布,其速度尺度由最大雷诺剪应力、自由流速度和Schofield-Perry速度尺度形成。     

气膜冷却恢复区的湍流模型

在 L am- Bremhorst低雷诺数 k- ε模型的基础上 ,发展了气膜冷却恢复区的湍流模型。原模型中的常数 cμ 和 c1代之以沿 x、y方向变化的参数 ,低雷诺数函数 fμ、f1和 f2 沿 x方向有所变化。修正后的模型可以在 0 .4～ 1.0的吹风比范围内 ,较好地描写恢复区内时均速度、湍动能和湍应力分布 ,并能正确模拟向不吹风情况过渡的恢复过程。     

马赫6柱-裙构型激波/湍流边界层干扰摩阻统计特性

为探究激波/边界层干扰表面摩阻统计特性,对马赫6柱-裙构型的激波/湍流边界层流场进行了直接数值模拟,推导了柱坐标形式的时均摩阻分解公式,将其与对流效应、流向不均匀性效应、分子黏性耗散效应、曲率效应和湍动能耗散关联起来,在此基础上对摩阻脉动和时均摩阻的统计特性进行了研究。概率密度分布结果表明干扰区内摩阻脉动概率密度明显偏离正态分布,摩阻脉动在该区域的间歇性较强,这与湍流脉动在激波干扰后存在显著内在可压缩性效应有关;谱功率密度结果表明摩阻脉动的能量在干扰前中频部分的能量占据主导,峰值对应频率约为0.14,干扰区后附近区域峰值位置转移到高频区域。采用摩阻分解公式对激波干扰前后的时均摩阻进行分解,结果表明时均摩阻的主要贡献来自分子黏性耗散效应和湍动能耗散效应;在激波干扰前分子黏性应力效应占据主导位置,干扰后湍动能耗散效应占据主导位置;激波干扰后对流效应明显增强,曲率效应的贡献在激波干扰前后基本保持不变;流向不均匀项贡献由正贡献变为负贡献,且所占比例有所上升,这主要是由压力梯度项的变化导致的。对分解后各项积分内函数沿法向分布的分析表明,分解后各项在边界层内分布呈现较大的差异,特别地,湍动能耗散项...     

基于Reynolds平均的高超声速二方程湍流模型

 在Reynolds平均的框架下推导了可压缩湍流Reynolds应力方程和湍动能方程。根据一定的假设和尺度分析简化并封闭了所推导的湍动能方程。在均匀湍流假设下,湍动能耗散率可分解成可压缩性耗散和旋度耗散,并对其中的可压缩性耗散进行了封闭;同时认为旋度耗散不受可压缩性影响,直接引用传统的Reynolds平均不可压缩湍动能耗散率模型方程。由此构造了适用于高马赫数的二方程可压缩湍流模型。应用所发展的模型计算了高超声速平板绕流,并与若干现有模型的计算、实验与半经验公式的计算结果进行了对比,验证了所发展的模型。在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对所发展的湍流模型,以及若干现有模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了湍流模型可压缩修正的必要性,得到了所发展模型的计算结果最接近实验结果的结论。     

湍流边界层大尺度相干运动的阵列TRPIV测量

使用高时间分辨粒子图像测速技术，研究湍流边界层中大尺度相干运动。由于大尺度运动的流向空间尺度与边界层厚度δ有关，因此沿流向排列4个高速相机进行拍摄，得到了约6.7δ×1.2δ的湍流边界层大视场，实验雷诺数Re_τ=422。针对流场中不同法向高度的流向脉动速度，采用沿流向方向进行空间小波变换的方法，得到不同空间尺度分量的脉动速度，并计算其占总流向脉动动能的比例，发现湍流边界层外区存在流向最大能量流向尺度，约为1δ。通过小波分解将湍流脉动速度场分为大尺度分量和小尺度分量。使用速度门限法，沿时间序列提取大尺度相干运动，利用泰勒冻结假设，将时间结构转化为空间结构，并与直接从空间得到的大尺度相干结构做对比。使用相位平均法测得大尺度相干结构的几何形态，发现从时间维度和直接从空间维度得到的喷射事件的流向尺度相近，而直接从空间提取的扫掠事件要比从时间提取的大。结果表明：流场中1δ尺度左右的大尺度运动是湍动能的主要贡献者；利用泰勒冻结假设可以从时间中提取出大尺度相干结构，与从流场空间直接提取的结果有着良好的一致性。     

应用一个新的具有转捩敏感性两方程低雷诺数湍流模型预测转捩流动

本文基于Launder-Sharma(LS)低雷诺数的 两方程湍流模型发展了一种新的具有转捩敏感性的低雷诺数湍流模型。针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,新模型进行了改进。新模型只使用当地物理量,不需要求解壁面距离、Y+和边界层积分参数。新模型能够适用于广泛的流动,且容易应用到通用的CFD程序中。对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,新模型能够预测转捩流动,并能对自由湍流变化给出合理的响应。 关键词：航空、航天推进系统;转捩模型;转捩流动;数值模拟     

沟槽壁湍流多尺度相干结构实验研究

应用热线测速技术,对沟槽壁面平板湍流边界层的减阻机理进行了实验研究.测量了风洞中并排放置的沟槽壁面平板及光滑壁面平板湍流边界层在不同雷诺数下不同法向位置的瞬时流向、法向速度分量的时间序列信号.运用流向速度分量信号的多尺度子波系数辨识壁湍流多尺度相干结构,用条件采样和相位平均技术提取了壁湍流多尺度相干结构喷射和扫掠时流向速度分量、法向速度分量和雷诺应力分量的相位平均波形.分析了沟槽壁面平板及光滑壁面平板湍流边界层中多尺度相干结构的持续时间、条件相位平均波形等特征.对沟槽壁面平板及光滑壁面平板湍流边界层多尺度相干结构的多种统计平均特征进行了比较,从壁湍流多尺度相干结构控制的角度研究了沟槽壁面平板湍流边界层的减阻机理.      

用子波系数概率密度函数研究湍流多尺度结构的间歇性

用热线风速仪在充分发展平板湍流边界层中测量流向速度信号的时间序列并进行湍流多尺度结构统计特性的实验研究。用热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量风洞中平板湍流边界层不同法向位置的流向速度分量的时间序列信号,用子波分析对湍流脉动速度信号进行多尺度分解,将推广的自相似性(ESS)应用于子波分析。用子波系数的概率密度函数验证了多尺度相干结构对间歇性的影响,并证明与距壁面距离有关。对比研究了提取湍流边界层多尺度相干结构前后湍流多尺度结构函数间歇性的变化规律及概率密度函数的变化规律,验证了多尺度相干结构是湍流边界层产生奇异标度律的原因。当这些相干结构被分离出后,由子波系数计算的ESS标度律就回到K olm ogorov线性标度律,这一现象在整个湍流边界层都是很明显的。      

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

周期性扰动对平板湍流边界层统计性质的影响

用IFA300热线风速仪和X形二分量热线探针精细测量了平板湍流边界层在施加不同频率的周期性吹吸扰动前后不同法向位置的瞬时流向、法向速度分量的时间序列信号.在平板湍流边界层的固壁表面沿展向切割一条5 mm宽的窄缝, 在平板底部用扬声器从壁面向平板湍流边界层内施加不同频率的周期性吹吸扰动, 研究不同频率的周期性扰动在平板湍流边界层内发展演化的规律及对平板湍流边界层统计性质的影响.      

湍流涡黏性模型方程中相位差的测量

基于理论上湍流相干结构复涡黏性模型对涡黏性系数的分析,应用热线测速技术测量了平板湍流边界层多尺度相干结构动力学方程中非相干结构成分对相干结构贡献的雷诺应力分量与相干结构流向速度流向变形率之间的相位差。分析了湍流边界层相干结构猝发过程中,非相干结构成分对相干结构贡献的雷诺应力分量与相干结构流向速度流向变形率之间的相位差沿湍流边界层法向的变化规律,为建立更加符合实际的湍流模型提供了实验依据。     

低湍流度下翼型边界层及近场尾流中湍流切应力和湍流动能间的相关

简述了在西北工业大学低湍流度风洞中，用热线风速仪对翼型边界层及近场尾流中雷诺正应力及切应力等的测量结果，着重讨论了雷诺切应力及湍流动能的分布规律及相关特性。结果表明，各区段的相关规律并不相同，与常规结果有明显差异。     

激波/平板湍流边界层干扰的数值模拟

采用理性GAO-YONG可压缩湍流模型,模拟了激波/平板湍流边界层干扰现象,结果表明:计算所得到的壁面压力分布、摩阻系数分布和速度型分布均与实验值吻合很好,并且比较准确地预报出了入射斜激波与平板湍流边界层干扰所引起的边界层分离点和再附点等流动特性.由此表明GAO-YONG可压缩湍流模型能够准确地用来模拟激波/平板湍流边界层的干扰流动.      

A study on turbulence transportation and modification of Spalart–Allmaras model for shock-wave/turbulent boundary layer interaction flow

It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent boundary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation(DNS). It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart–Allmaras(S–A) model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S–A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S–A model are obviously improved. Thus it can be concluded that the modification of S–A model with the pressure gradient can improve the predictive accuracy for simulating the shock-wave/turbulent boundary layer interaction.     

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

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

翼型尾缘噪声源空间分布与辐射特性关系

详细研究了翼型湍流边界层尾缘宽频噪声源空间分布与辐射特性的关系.采用基于雷诺平均流场的翼型尾缘宽频预测方法研究了NACA0012翼型湍流边界层尾缘宽频噪声在4种不同工况下的噪声源空间分布与辐射特性.首先计算了NACA0012翼型湍流边界层尾缘噪声源在不同频率下的空间分布.计算结果发现:边界层中湍流是翼型湍流边界层尾缘噪声声源.随着频率的增加,噪声源强度和噪声源空间尺寸都是先增加后减小,噪声源位置不断靠近翼型尾缘.同时也计算了边界层内不同位置处的噪声源对远场噪声的辐射特性,结果表明:边界层内层区域,其噪声频谱能量集中在高频;边界层外部区域,其噪声能量集中在中低频;攻角增大或者来流速度减小,噪声能量向低频转移.      

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。