边界层转捩的数值模拟

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

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

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

压气机静叶栅层流分离泡转捩与角区分离数值模拟与实验

采用γ-Reθ转捩模型对某可控扩散叶型(CDA)平面叶栅全攻角范围进行了三维数值计算，通过对比数值计算结果与叶栅实验吻合较好。在此基础上，分析了进口来流湍流度和雷诺数变化对叶栅表面层流分离、转捩以及角区分离的影响。结果表明:进口湍流度低于5%时，吸力面存在层流分离，当进口湍流度大于5%后，层流分离移除，但转捩会一直存在；随着进口湍流度或雷诺数增加，吸力面和压力面转捩位置均会前移；随着进口湍流度增加，吸力面角区分离会有所减小，雷诺数增加对角区分离的影响不大。      

低压涡轮叶型边界层相互作用的数值模拟

基于Lantry-Menter转捩模型,应用商用流体计算软件求解非定常雷诺平均N-S方程组,对进口雷诺数为2.5×104,来流湍流度为2.5%,尾迹折合频率为0.92状态下的尾迹/PAK-B叶型边界层的相互作用进行了数值模拟。数值计算揭示的尾迹在叶栅通道中的输运特性、尾迹诱导卷起涡形成、尾迹诱导边界层转捩等物理机制大致符合相关实验定性描述。     

基于浸入式边界法的振荡转子叶片数值模拟

对于流固耦合问题,传统的数值模拟方法由于需要不断的网格重构使得问题变得复杂,本文通过利用浸入式边界法建立了在单一坐标系和网格下求解振荡转子叶片的快速计算模型,从而避免了传统方法中由于需要不断地重构贴体网格造成的数值模拟的复杂性。为了验证该方法的正确性,分别对低KC(Keulegan-Carpenter)数下的振荡圆柱和两自由度振荡圆柱进行了数值模拟,计算结果与以往的试验结果和数值结果吻合得很好,证明了浸入式边界方法的可靠性。在此基础上,对振荡转子叶片进行了数值模拟。结果表明折合速度和叶栅稠度是影响振荡转子叶片的重要因素。并且,值得注意的是这种耦合过程并没有生成贴体网格,因此大大减小了计算时间,这样可以更加快速、准确地模拟真实的叶栅流动情况。     

用N-S方程计算三维叶栅流

发展了基于微元中心有限体积方法,采用雷诺平均N-S方程求解三维叶栅流动的程序,特别注意了格式中人工粘性的加法;程序采用Baldwin-Lomax模型。格式经过平板层流及湍流流动的校验,层流与Blasius解、湍流与实验经验公式及Wilcox的计算值符合很好,最后以Goldman涡轮导向器为例进行了叶栅流动计算校核,结果符合于实验以及Chima等人的计算。      

低雷诺数下涡轮叶栅流动分离实验与数值模拟

应用实验测量和数值模拟相结合的方法,研究了低雷诺数条件下高负荷涡轮叶栅吸力面的流动分离。通过对叶片表面压力系数、叶栅出口尾迹以及叶片表面气流分离位置和重新附着位置的比较发现,计算结果与实验结果吻合得相当好。应用本计算方法,对低雷诺数条件下雷诺数和来流湍流度对涡轮叶栅的流场的影响作了准确的模拟,对叶栅吸力面的气流分离、再附等做出了预测。实验研究和计算结果都表明,低雷诺数条件下叶栅损失的急剧增大是由于在低雷诺数条件下叶片吸力面发生了气流的分离,雷诺数越低或者进口湍流度越低,叶片吸力面的气流分离就越严重,由此导致的叶栅损失也就越大。     

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

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

绕翼型低雷诺数流动的数值分析研究

分离泡的产生是绕翼型低雷诺数流动的一个重要特征,分离泡通常都是非定常的,会对整个流场产生极大的影响,由于分离附面层的不稳定性,很快诱发转捩,并产生湍流再附。文中通过求解雷诺平均N-S方程,数值模拟了绕Eppler387翼型的低雷诺数非定常流动,并对两方程SSTk-ω湍流模型,代数B-L模型和层流的计算结果作了比较。N-S方程和两方程湍流模型的控制方程非耦合求解,时间推进均采用近似因子方法(AF);空间离散无粘项采用ROE格式,粘性项用中心差分方法计算;计算非定常流场时,采用伪时间子迭代(τt-s)方法保证二阶时间精度,湍流模型计算时都是在固定点转捩。最后分析了转捩对低雷诺数流动的影响、分离泡的存在引起的流动不稳定性和周期性的脱出涡,数值计算给出的时间平均升力系数、阻力系数和压力分布与实验结果比较吻合。     

二维涡轮叶栅流场的直接数值模拟

采用高精度有限差分格式求解非定常N-S方程组,对低雷诺数下二维涡轮叶栅流动进行了直接数值模拟,计算了雷诺数为10000,VKI涡轮叶栅在0°,8°以及-8°攻角下的流场,对涡轮叶栅非定常流动机理做了初步的探讨。计算结果表明:在叶栅尾缘处,逆时针方向和顺时针方向的主涡交替在壁面产生,并和主流相互作用产生二次涡,而当二次涡与主流连通发生掺混时,将引起主涡被分割并从叶片表面脱落;攻角在一定范围内的变化对VKI涡轮叶片表面边界层发展影响不明显。文中还对尾迹区的统计量特性和速度亏损特性等进行了研究。      

来流马赫数对附面层的影响

本文从附面层积分方程出发, 定性地分析了当表面压力系数 C_p分布相同时, 在保证Re数相等的条件下进口 Ma数对附面层发展的影响。并且用数值方法进行了模拟。与此同时用数值方法对由于表面曲率和自由流紊流度所引起的影响作出估计。结果表明; 在表面压力系数分布和 Re相同的情形下, 进口Ma数, 直到临界值附近时, 所引起的附面层的变化量是很小的, 一般在工程允许的误差范围之内。表面曲率所引起的变化同马赫数引起的变化相近, 而进口自由流紊流度引起的变化较之前两者较为明显。      

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。     

MHD控制激波诱导湍流边界层分离的机理分析

为了研究磁流体动力学(Magnetohydrodynamics:MHD)加速边界层对激波-湍流边界层相互作用的影响,用高阶有限差分法求解了小磁雷诺数近似的MHD湍流方程。其中,无粘通量采用WENN格式离散、粘性通量采用Roe平均中心差分离散,时间采用半隐式推进,并采取追赶法求解。计算给出了湍流、电场、磁场和电导率等参数对边界层分离的影响,数值结果显示:在同样的逆压梯度下,湍流边界层分离能更快地趋于稳态流场,且分离区比层流小;通过施加洛仑兹力加速,边界层速度型面变得更加饱满、位移厚度减小、分离点和再附点向激波与固壁的交点靠近,分离区尺寸减小甚至最终被消除。     

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

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

Steady and unsteady numerical investigation of transitional shock-boundary-layer-interactions on a fan blade

The present numerical study was performed to investigate the impact of both the Reynolds number variation and the used turbulence model to capture the boundary layer development on the characteristic of a BR710 fan blade. A one-equation model of Spalart Allmaras with an optional semi-empirical transition model of Abu-Ghanam Shaw has been applied. The transition model allows the boundary layer development from a laminar to a turbulent behaviour to be taken into account. This is of particular importance at low Reynolds numbers and thus high operating altitudes of airplanes when the transition location moves further downstream to the trailing edge. Therefore the interaction between the shock waves of the transonic fan blade row and the boundary layer leads to a significant change of the characteristic. The steady numerical studies of the 3-D blade passage have been carried out with the commercial 3-D Navier–Stokes solver NUMECA. Additionally an unsteady calculation has been applied to explain the penalty in efficiency on high operating altitudes.     

S-A模型在分离流动模拟中的改进

基于湍流不平衡性的分析,通过修正雷诺应力以及模型系数的计算,提高了S-A模型在二维边界层以及三维叶栅的分离流动上的预估能力。修正后的S-A模型给出了与实验吻合较好的结果,在对于分离流动的预估能力上优于原始模型。     

超临界层流机翼边界层及气动特性分析

高空长航时无人机设计巡航状态的雷诺数较小,黏性边界层对气动特性的影响较大。详细分析了雷诺数对机翼边界层和气动力的影响,用数值方法对超临界层流机翼三维层流-转捩-湍流混合边界层特性进行了研究,分析比较了高空小雷诺数和中空大雷诺数情况下机翼三维边界层的特性,尤其是边界层转捩点位置、表面摩阻和气动特性的雷诺数效应。研究表明雷诺数对于高空无人机机翼边界层厚度、摩擦阻力和升阻比影响较大;对层流机翼的转捩点位置和升力系数影响较小;自然层流机翼技术可以应用于高空无人机设计。     

湍流度和粗糙度对边界层转捩影响的研究

工程上计算边界层转捩位置,常用以线性稳定性理论为基础的e~N方法。实际上,层流边界层转捩受坏境因素(如湍流度T,粗糙度k和压力梯度等)影响较大。所以作为e~N方法的补充,须寻求一种能同时计算坏境因素影响的方法。现只有湍流模型方法能满足这种要求。     

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。