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

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

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

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

比熵增概念及其在湍流模型中的应用

通过对不同马赫数下(Ma=0.7、2.25、6)平板边界层进行直接数值模拟研究,提出比熵增的概念.比熵增表征单位机械能的耗散,并且在数值上表现出较好的单调性和马赫数无关性,因而能够可靠地表征边界层的范围.应用比熵增概念修正BL零方程湍流模型长度尺度,构造BL-entropy,同时采用原始的BL模型,SA一方程模型对后台阶(Ma=0.128)以及高超声速锥柱裙组合体(Ma=7.05)流动进行数值模拟,并与实验数据进行比较,结果表明BL-entropy能够准确地获取长度尺度进而得到均匀合理的涡粘性分布,极大地提升了原始BL模型模拟复杂流动的能力.     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

不可压缩湍流的色散模型

分析了广泛存在于湍流运动中的能量逆转现象,揭示了其产生的原因,提出了湍流的色散性质.在此基础上,修改了Boussinesq假设,建立了包含色散效应的新的雷诺应力封闭式和湍流色散系数,给出与不同模型相结合时,湍流色散系数所具有的不同形式,并阐述了湍流运动中能量的传递方向及条件.采用不可压缩平板边界层流动和平面后台阶流动验证了其可信性和优越性.平板摩擦阻力系数及边界层速度型与实验结果吻合良好,平面后台阶流动的流向再附长度、台阶边压力系数及湍流强度等参数均比标准k-ε模型更接近实验结果.结果表明:色散项的加入可以在不显著增加计算量的同时显著改善预测精度,模型具有一定的工程应用价值.      

带有前缘气膜冷却燃气轮机叶片的气热耦合研究(英文)

使用气热耦合数值计算方法对带有"淋浴喷头"式前缘气膜冷却的C3X燃气轮机叶片进行数值模拟。通过对使用不同湍流模型的计算结果与实验测量值进行对比,使用转捩模型能够比较准确的模拟带有气膜冷却射流的边界层内的流动和传热过程。气膜冷却射流在叶片吸力面上游层流边界层区域内表现出比较明显的三维湍流特性,但在吸力面下游边界层转捩为湍流以后,湍流使得气膜冷却射流与主流燃气掺混过程加速而使其三维流动特性迅速消失;在吸力面由于边界层流动始终为层流状态,冷却射流与主流掺混较缓慢而使其仍能保持较强的三维湍流特性。边界层内由高温燃气与被冷却叶片之间温度差形成的温度梯度有益于增强边界层流动的稳定性。     

求解可压缩流动的同位网格SIMPLE方法研究

在Rhie-Chou动量插值的基础上,推导了同位网格可压缩SIMPLE算法.经过无粘流超音速凸包算例和激波/湍流边界层干扰算例计算发现,如果对流项采用高阶有界HLPA格式,密度插值采用一阶迎风和中心差分的混合格式,这种算法能够很好地模拟凸包超音流的流动现象,在采用了新型GAO-YONG湍流模型后也能够较好地模拟激波/湍流边界层干扰.      

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

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

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

ALGEBRAIC TURBULENCE MODEL WITH MEMORY FOR COMPUTATION OF 3-D TURBULENT BOUNDARY LAYERS WITH VALIDATION

Additional equations were found based on experiments for an algebraic turbulence model to improve the prediction of the behavior of three dimensional turbulent boundary layers by taking account of the effects of pressure gradient and the historical variation of eddy viscosity, so the model is with memory. Numerical calculation by solving boundary layer equations was carried.out for the five pressure driven three dimensional turbulent boundary layers developed on flat plates, swept-wing, and prolate spheroid in symmetrical plane. Comparing the computational results with the experimental data, it is obvious that the prediction will be more accurate if the proposed closure equations are used, especially for the turbulent shear stresses.     

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.     

高超声速锥柱裙模型边界层转捩的弹道靶实验

为研究高超声速边界层转捩现象、给边界层计算提供可靠的对比数据,在中国空气动力研究与发展中心趟高速弹道靶上开展了锥柱裙模型高超声速边界层转捩的自由飞实验。所采用的锥柱裙模型全长105mm,飞行速度1．94km／s（Ma=5．65）,单位雷诺数4．32×10^7～1．20×10^8m-1。使用激光阴影成像技术,获得了锥柱裙模型边界层转捩和湍流边界层发展的图像,测得的湍流边界层厚度在0．6～2．2mm之间,湍流涡的流向尺寸与边界层厚度的比值介于0．3～0．8之间且沿流向呈下降趋势。实验结果表明：弹道靶实验能够获得给定飞行环境下的高超声速边界层转捩图像,从图像中可以清晰判断转捩位置或区域、测量边界层厚度和分析湍流涡的尺寸。     

离散气膜孔下游边界层的实验研究

“Ｘ”型热线测量了离散气膜孔下游边界层内的湍流量,将热线探头绕轴转三个角方位以感受三维湍流场信息。结果表明,边界层内湍动能分布呈反“Ｓ”形,吹风比和孔排数对参数分布有重要影响,这种影响在很长距离上不消失,所以,常用的ｋ－ε模型不经修正不能正确预计这类边界层的发展。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

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

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

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

旋涡运动与湍流模型

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

边界层转捩的数值模拟

在数值模拟中将Abu-Ghannam&Shaw(AGS)转捩模型引入到湍流模型中,以考虑流动中的层流和转捩现象,并在此基础上,分别对低雷诺数情况下的平板和涡轮叶栅流动进行了模拟,并与实验结果进行了比较,表明加入转捩模型可以较好地描述边界层的转捩过程。同时,结果还表明AGS转捩模型对压力梯度变化较大情况下的模拟稍有偏差。