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

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

应用LDV研究二元管道内分离流动

应用二维激光多普勒仪测量了二元曲壁非对称扩张通道内的湍流附面层分离流动。得到了时均速度、雷诺应力、平均流动方向角及正反向间歇流动因子的分布。文中给出了所测结果的不确定度,并在激光测速的近壁测量上做了尝试。实验结果分析表明:间歇瞬时分离点对应于附面层内垂直压强梯度的最大值;Schofield的速度尺度和长度尺度可用于形成新的速度相似型;正反向间歇流动因子γpuo与平均流动方向角存在着简单的线性关系。     

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

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

三维弯管内分离流场的试验研究

通过热线测速系统对三维９０°弯管分离流场的气流特性进行了实验研究,通过对分离流场的测量及观察,给出了时均速度、脉动速度、雷诺应力的大小等数据。结果表明：在管道的中心线附近远离分离泡区域的紊流可以看作各向同性的;而分离泡与主流的过渡区域表现为强烈的各向异性;回流区的平均速度场分布与紊流度无关,回流的雷诺剪应力大小与当地平均速度无关,回流不是来自很远的下游。     

对称翼型近场尾流及后缘边界层流动特性的试验研究

本文简述了用热线风速仪测量“NACA 63-012”对称翼型近场尾流及后缘边界层流动特性,讨论了流场的时均特性及雷诺应力等湍流量的分布特性,提出了无因次速度分布的半经验公式、相关参数及相关曲线解析式,并得出了几点重要结论。     

吹吸扰动控制猝发过程中雷诺应力的实验

采用热线测速技术和双丝边界层热线探针,在风洞中以高于对应最小湍流时间尺度的分辨率,精细测量了平板湍流边界层施加不同频率周期性质量引射前后,不同流向和法向位置的瞬时流向、法向和展向速度分量的时间序列信号.采用流向脉动速度局部平均结构函数模极大值法检测了壁湍流猝发过程中拟序结构不同速度分量和雷诺应力分量的条件相位平均波形.发现在扰动源下游一定范围内,摩擦因数和雷诺应力幅值有所减小,扰动改变了法向或展向脉动速度与流向脉动速度的条件相位平均波形的相位差,使其保持在π/2的奇数倍左右,此时雷诺应力几乎为0,改变了湍流边界层原有雷诺应力的产生和扩散输运机制,抑制了湍流雷诺应力的产生和输运.      

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

三维边界层的湍流特性及其二阶矩模拟

根据现有实验数据,本文进一步分析了三维边界层内纵向涡对雷诺应力分布的影响,从而建立了一个足以模拟三维边界层湍流特性的二阶矩封闭模式。在新模式中,着重考虑了三维平均涡场引起的雷诺应力输运的非线性作用。以文[3]给出的模拟无限后掠翼实验为例,本文分别应用基准的二阶矩模式和新建立的非线性二阶矩模式对三维边界层内的雷诺应力分布作了数值预计。结果表明:应用基准二阶矩模式所得雷诺应力诸分量的预计值均比相应的测量值大得多,且完全不能预计切应力方向相对于速度梯度方向的滞后;而非线性二阶矩模式对雷诺应力分布的预计则基本上与测量数据符合,且能定性地模拟切应力方向相对于速度梯度方向的滞后。     

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

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

斜出口合成射流激励器横流输运特性与边界层控制

研制了一种斜出口合成射流(BSJ)激励器,应用粒子图像测速(PIV)技术获得了合成射流的瞬态和时均流动结构,分析比较了斜出口和常规平直出口激励器流场的异同。该斜出口激励器能将周围气体进行有方向的能量和质量输送控制,时均流场呈现沿壁面的横向流动输运特性。解释了横向流动输运的形成机理,并应用该激励器进行了边界层控制实验。研究结果表明:通过改变激励器的工作电压和频率,可方便有效地实现对平板边界层速度型的调控,使受控边界层速度型更"饱满",这对控制推迟边界层分离非常有利。     

一种关于静压气体轴承节流孔系数的计算方法

基于层流边界层方程的分离变量算法和雷诺方程的解析算法，提出了一种关于单节流孔静压气体止推轴承的节流孔系数的计算方法。该方法通过比较层流边界层方程计算获得的气体轴承的质量流量和雷诺方程计算获得的质量流量计算获得了节流孔系数。将计算获得的节流孔系数和节流孔系数为常数0.8代入单节流孔气体止推轴承的雷诺方程中，计算获得的承载力与分离变量算法求解层流边界层方程获得的承载力进行对比，可以发现，相对于采用节流孔系数为0.8来说， 采用该计算的节流孔系数求解雷诺方程的承载力与分离变量算法求解获得的承载力结果精度最大提高了8%。从而验证了该计算节流孔系数方法的正确性。      

高超声速有攻角锥裙直接数值模拟

采用直接数值模拟方法对有攻角的高超声速7°～34°锥裙开展了数值研究,通过对比0°、90°、180°周向子午面,评估了三维横流效应对激波/边界层干扰的影响规律和作用机制,包括壁面压力、摩阻、热流分布,分离泡非定常运动,再附边界层演化等。研究发现,不同周向站位均出现流动分离,横流区、迎风区内复杂激波结构与边界层相互作用导致壁面压力、摩阻、热流显著升高。热流与压力的比值在干扰区上升后由于再附降低,而热流与摩阻的雷诺比拟关系在分离区则完全失效。分离泡面积脉动的功率谱结果表明,分离泡非定常膨胀/收缩运动呈低频特征,且分离泡呼吸与激波低频振荡在横流区密切相关,在迎风区存在迟滞,而在背风区不相关。速度脉动场的本征正交分解结果表明,分离区的低频特征与低阶模态相应的剪切层附近大尺度结构相关。对下游再附边界层演化分析指出,攻角的存在导致雷诺应力在再附点附近大幅增强,其流向分量的恢复最为迅速,雷诺应力分量的峰值位置在背风区沿流向持续外移,而在迎风区、横流区已迅速向内层恢复。此外,雷诺应力各向异性不变量分布进一步表明干扰下游的近壁区湍流各向异性峰值在背风区弱于迎风区。     

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

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

圆转矩形喷管射流掺混特性试验研究

本文采用热线风速仪对圆转矩形收敛喷管的射流流场与掺混特性进行了试验研究,测量了喷口下游不同截面上,射流宽、窄对称面上射流的速度、雷诺剪应力、湍流强度,得到其沿径向和轴向的分布规律。发现在径向方向上,雷诺剪应力先增大后减小,最大值出现在射流与外流的交界面附近,湍流强度则逐渐减小;在轴向上,雷诺剪应力在近喷口区存在小幅波动,而后逐渐减小,湍流度则先略有增大而后变化逐渐减慢。射流特性在宽对称面与窄对称面上的分布规律相同。     

运用低雷诺数非线性k-ε模型研究湍流分离流动

给出了一个低雷诺数非线性ｋ－ε湍流模型，其中采用了一个非线性雷诺应力表达式并引入了两个考虑湍流近壁效应的修正函数。应用该模型和标准ｋ－ε模型计算了两种不同膨胀比的后台阶分离流。与实测结果比较表明，低雷诺数非线性模型结果较标准ｋ－ε模型结果有明显改进，它对流场的各特征量包括分离—再附长度、时均速度及各雷诺应力分量均给出了较好的预测。     

圆转矩形喷管射流掺混特性试验

采用热线风速仪对圆转矩形收敛喷管的射流流场与掺混特性进行了试验研究,测量了喷口下游不同截面上,射流宽、窄对称面上的速度、雷诺剪应力、湍流强度的分布规律.发现沿径向,雷诺剪应力先增大后减小,最大值出现在射流与外流的交界面附近,湍流强度则逐渐减小;在轴向上,雷诺剪应力在近喷口区存在小幅波动,而后逐渐减小,湍流度则先略有增大而后变化逐渐减慢.射流特性在宽对称面与窄对称面上的分布规律相同.      

SSG模型在压缩拐角激波/湍流边界层干扰模拟中的应用

采用非线性方案模化雷诺应力再分配项的二阶矩雷诺应力输运模型——SSG模型和两方程线性涡粘模型SST模型,对24°高雷诺数二维压缩拐角激波/湍流边界层相互干扰流动进行数值模拟。通过在壁面压力分布、壁面摩擦阻力系数分布、截面速度分布及分离区大小等方面与实验数据比较,综合评估了上述两种模型模拟该问题的能力。结果表明:两种湍流模型在上述各方面都取得了与实验较为一致的结果,但相比较而言,除SST模型在流动再附后对壁面摩擦阻力系数分布的模拟更接近实验值外,SSG模型结果在上述其他方面均更接近实验值。本文工作为压缩拐角激波/湍流边界层干扰这类流动数值模拟的湍流模型选择提供了参考。     

湍流边界层速度分布的显式表示

本文根据不可压缩湍流边界层内的Ｒｅｉｃｈａｒｄｔ立方律，得到了比现有的Ａ．Ｊ．Ｍｕｓｋｅｒ和Ａ．Ｌｉａｋｏｐｏｕｌｏｓ的公式更加符合Ｓｐａｌｄｉｎｇ公式及实验数据的湍流边界层内速度分布的显式表达式。