在直机翼与平板连接处测定分离区内外的低、亚音速湍流边界层特性

本文介绍了在低速和亚音速来流情况下,对一竖立着机翼的平板,在其连接处的分离区内外所作的三维湍流边界层测量。低速实验测出沿二条不同站位线共十个测量站位上的三个平均速度分量分布和六个湍流应力分量分布。亚音速情况下,对分离区内外的二个站位作了同样测量。测量结果表明,气流的分离导致产生一个马蹄形旋涡,它是影响分离区内外三维湍流边界层特性的重要因素。文中还讨论了实验数据的处理方法,分析了有压力梯度和流线曲率影响的实验结果。     

亚音速扩压器内湍流边界层特性的试验研究

本文介绍了小展弦比二元直壁亚音速扩压器内的湍流边界层特性.精确地测量了扩压器沿程几个截面上的平均速度分布、壁面摩擦系数C_f分布以及由此而得出的边界层积分参数.并指出C_f要比Ludwelg-Ti-Ilmann经验公式所计算的结果稍低,所以要对C_f稍作修正.文章还分析了逆压力梯度对湍流边界层内外区的壁面律与尾迹律的影响,并得出无因次压力梯度参数与尾流强度参数之间的关系曲线.     

绕扁长椭球体三维分离的积分法

本文利用三维积分边界层法并结合流线法决定涡层型流动分离。由三维不可压湍流边界层方程求得边界层参数后,按流线法计算粘性流线。作为算例,研究了一带迎角的扁长椭球体在不可压缩湍流里的绕流,经与水洞里所观察到的α=15°流型的对比,二者基本一致。     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

流向弯曲壁超声速湍流边界层研究进展

流向弯曲壁面上的超声速湍流边界层特性是流体力学的重要科学问题,掌握其规律对提升高超声速飞行器设计水平有重要意义。国内外在该问题上的研究已经持续了50余年,近年来随着实验测量和数值模拟能力的提高,在流向曲率的影响机理上取得了突破性进展。本文系统梳理了流向凸曲壁及凹曲壁的超声速湍流边界层的演化机制,介绍了流向曲率、压力梯度和体积膨胀/压缩等关键因素的影响规律,总结了流向凸曲壁边界层湍流衰减的特性与成因,以及凹曲壁边界层中G?rtler不稳定和逆压梯度增强湍流的特性与成因。最后归纳其发展趋势,为下一步开展弯曲壁超声速湍流边界层研究提出了建议。     

超音速绕凹角湍流分离流动的计算

本文用统一的Levy-Lees变换以及正算法与逆算法相结合,求解了超音速绕凹角湍流分离流动。 对附着流区用边界层正算法,压强分布用流过尖劈统一的高超音速与超音速公式,湍流模型取代数涡粘性模型;对凹角分离区用边界层逆算法,给定位移厚度δ~*分布,湍流模型取代数松弛模型;边界层计算采用Cebeci-Keller Box方法;计算成功地算得分离流场,较好地预估了分离点与重附点位置以及壁面压强分布与表面摩擦应力分布。     

低雷诺数κ-ε封闭模式在三维湍流边界层计算中的应用

本文采用低雷诺数形式的k-ε湍流模型,计算具有大横向流的不可压缩三维湍流边界层。在湍流模型中,应用了涡粘性的各向异性系数。本文所得结果与用混合长模型所得结果以及量测结果的比较表明:对于边界层的不同特征量,湍流模型和涡粘性各向异性对计算结果的影响程度是不同的。总的来说,应用考虑涡粘性各向异性的k-ε模型,较之应用混合长模型,所得计算结果有所改善,但还不是根本性的改善。     

无限翼展后掠机翼的三元湍流边界层计算

本文应用零方程的湍流模型,考虑到三元湍流边界层中涡团粘性的各向异性分布,分别对Michel-Schneider模型和Cebeci-Smith分层模型,采用有限差分和牛顿迭代方法,对有逆压梯度的无限后掠平板机翼湍流边界层进行了计算。文章着重讨论了混合长度(在Cebeci-Smith模型中包括内层参数和外层参数)和涡团粘性各向异性系数的影响,并与风洞实验测量值进行了比较。     

翼型近尾迹流动的PIV研究—运动学特性

利用在线式PIV系统(ParticleImageVelocimetry),在低速风洞中对NACA0012翼型在雷诺数2.39×105,0°和4°攻角下的近尾迹流动进行了实验研究。实验结果表明,在较高的雷诺数下翼型近尾迹流动是一种以旋涡的运动学和动力学特性为主导的湍流剪切流。在测量范围内,翼型的尾缘处是近尾迹涡街的形成区;尾缘后0.5倍弦长的区域存在类似于卡门涡街的有序结构,是旋涡发展区域,旋涡具有较好的稳定性;距翼型尾缘0.5倍弦长至1倍弦长的区域,是翼型近尾迹流动由有序走向无序区域,旋涡开始破裂。翼型表面边界层对翼型近尾迹湍流剪切流的演化有重要影响。实验结果还给出了近尾迹流动的平均速度、湍流强度和剪切应变变化率,以及速度脉动量的二阶关联量u'u',u'v'和v'v' 的分布。      

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

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

主流压力梯度对气膜冷却效率影响的数值研究

分别对收缩通道、扩张通道和直通道中亚声速主流条件下的气膜冷却进行数值模拟，对比分析了不同主流压力梯度、次流吹风比条件下的主流和次流流场、温度场特征。研究结果表明:引起气膜冷却效率变化和不同发展趋势的因素可归结为主流边界层厚度、主次流自由剪切混合程度、肾形涡的强度和位置等因素。相对于零压力梯度的主流条件，在吹风比较小(M=0.25)的情况下，主流的逆压力梯度一方面增厚边界层、增强了气膜射流对主流的穿透，另一方面减小了肾形涡的强度，综合作用的结果是气膜平均冷却效率提高了4.91%。在吹风比较大(M=2)的情况下，主流的顺压力梯度扼制主流边界层的发展、抑制气膜射流的穿透能力，降低肾形涡涡核的位置，从而提高气膜冷却效率达17.40%。      

扩压器分离流场控制的研究

 本文研究中采用加置局部柱形扰源的方法,改变流场中的时均速度分布和紊流度分布,从而达到抑制分离,改善扩压器出口流场的目的。研究表明,局部柱形扰源可明显地延迟流场分离的形成,其主要机理在于近壁部分时均气流动量的变化。文中也介绍了紊流场的测定结果,并对附着流和分离流情况以及有无扰源条件下的紊流分布进行了分析讨论。     

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

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

The laminar-turbulent transition zone in the boundary layer

The flow during transition from the laminar to a turbulent state in a boundary layer is best described through the distribution of the intermittency. In constant-pressure, two-dimensional flow, turbulent spots appear to propagate linearly; the hypothesis of concentrated breakdown, together with Emmons's theory, leads to an adequate model for the intermittency distribution over flow regimes ranging all the way from low subsonic to hypersonic speeds. However, when the pressure gradient is not zero, or when the flow is not two-dimensional, spot propagation characteristics are more complicated. The resulting intermittency distributions often show peculiarities that may be best viewed as ‘subtransitions’. Previous experimental results in such situations are reviewed and recent results and models are discussed. The problem of predicting the onset of transition remains difficult, but is outside the scope of the present article.Although this paper is intended to be chiefly a survey, several new results in various stages of publication are also included.     

跨音速弯曲壁面正激波与湍流边界层干涉的计算结果

本文介绍弯曲壁面上湍流边界层与正激波干涉的某些计算结果。壁面曲率对激波前后的压力比有重要影响。在干涉区中,边界层位移、动量和能量厚度都有显著的增加。这种增加主要决定于激波前后的压力比。(TBTF)     

EXPERIMENTAL STUDY OF TURBULENT WINGBODY JUNCTION AND WAKE FLOW

ＥＸＰＥＲＩＭＥＮＴＡＬ　ＳＴＵＤＹ　ＯＦ　ＴＵＲＢＵＬＥＮＴ　ＷＩＮＧＢＯＤＹ　ＪＵＮＣＴＩＯＮ　ＡＮＤ　ＷＡＫＥＦＬＯＷＫａｎｇＨｏｎｇｗｅｎ；ＦａｎＭｅｎｇ；ＷｅｉＺｈｏｎｇｌｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＭｅｃｈａｎｉｃｓ，ＰｅｋｉｎｇＵｎｉ...     

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.     

后向台阶湍流边界层分离、再附及发展

应用激光测速仪对后向台阶湍流边界层分离、再附及发展进行了测量,得到时均流速、流向湍流度、平坦系数和偏斜系数等沿程分布。依据时均流速分布特性,得到分离区边界及再附点位置ｘ＿ｒ＝５．２ｈ＿ｓ（ｈ＿ｓ为台阶高度）,再附后湍流边界层恢复速度型相似的位置ｘ＿Ｔ＝１３．５ｈ＿ｓ。在详细分析湍流边界层区域时均及紊动特性的基础上,与光滑平板湍流边界层进行了比较。结果表明,除流速分布对数公式中积分常数较大外,无其它差别。此外,湍流边界层区存在尾流律,尾流参数π＝０．２９。     

静叶时序对压气机叶片附面层流动影响的数值研究

采用数值方法对两级低速压气机中径处的非定常流场进行模拟,针对压气机第2排静叶两个典型周向位置对动静叶干扰下的叶片附面层流动进行研究.建立尾迹与附面层干扰分析模型,结合叶片壁面摩擦力和近壁面附面层湍动能,详细分析了尾迹和势流干扰下静叶时序改变对叶片附面层流动产生的影响.对第2排静叶附面层的研究结果表明:静叶时序改变了尾迹在其叶排中的输运特征,能够降低壁面摩擦力和近壁面湍动能及其非定常最大波动幅值,影响吸力面附面层内动叶尾迹后沉寂区的宽度.在非定常条件下,尾迹能够诱导静叶层流附面层在尾迹干扰的局部范围内转捩发展为湍流状态,同时高湍流度尾迹的干扰具有抑制逆压梯度下附面层分离的作用,并能够延长层流区的范围.     

改善二元亚音扩压器性能的试验研究

本文按文献[1]提出的“失速裕度”设计法计算了二元亚音扩压器的几何尺寸和性能,其结果与作者所做的试验数据相差甚大.该扩压器是具有很大初始扩压角的突扩扩压器,在初始扩压段即出现严重的气流分离.此方法不适用的原因在于所用的附面层动量积分方程没有考虑壁面曲率项的影响.作者还对双扩压角直壁扩压器的性能进行试验研究,结果表明具有较大扩压角且短的初始扩压段的直壁扩压器性能优于相同面积比和长度的单扩压角直壁扩压器.此外,试验表明在分离点前安装适当型式和结构参数的叶片式涡流发生器能够控制扩压器内气流分离,从而使压力恢复系数C_p值提高25～34%.