水平轴风力机偏航与俯仰尾流特性

风力机的尾流干扰效应是影响风场下游风力机性能和风场出力的主要因素，通过调整风轮姿态可有效改善尾流干扰效应，但同时也引入了更为复杂的非定常气动特性和尾流场分布。为了研究风轮姿态控制的尾流特性，基于升力线和自由涡尾迹模型建立了水平轴风力机气动和尾流仿真模型，以NREL 5MW风力机为研究对象，在额定风速条件下，研究了叶轮偏航和俯仰工况下的尾流结构特性。在此基础上，对存在风切变和塔架影响时的风力机偏航和俯仰工况产生的尾流干扰进行了分析。研究结果表明：风轮的偏航和俯仰会引起尾流涡结构失稳提前、远尾流区湍流度增大，同时会导致尾流中心发生偏转和尾流作用区域减小，并加快速度亏损恢复，降低尾流效应对下游的干扰。风切变会产生类似风力机俯仰的效果，使得尾流抬升，作用区域减小。塔架对近尾流区的干扰较为显著，而俯仰控制会降低塔架对尾流的干扰。     

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

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

上游尾迹入射角对压气机叶片负荷的影响

压气机叶片负荷的提高使得叶片表面边界层更容易分离,利用上游叶排产生的非定常尾迹能够抑制边界层分离.运动圆柱代替上游转子,在保证下游叶片进气速度大小及攻角不变的情况下,改变圆柱运动速度以获得上游尾迹与下游静叶吸力面不同的夹角,发现在低负荷小分离情况下静叶损失系数与上游尾迹入射角无关;高负荷大分离情况下静叶损失系数随上游尾迹入射角的增加而降低.分析尾迹作用下高负荷静叶通道内流场,当进入静叶通道的尾迹与叶片吸力面近似平行时,尾迹诱导边界层增厚.使得叶片表面分离泡随时间大幅值脉动,损失增加;当尾迹与叶片吸力面的角度逐渐垂直时,尾迹抑制了边界层分离,同时叶片表面分离泡位置近似不变,其原因是尾迹以负射流形式进入边界层内部,补充了边界层内部低能流体,使得损失减小.      

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

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

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

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

圆柱尾迹影响下近壁湍流低速条带的特征

利用氢气泡流动显示方法，对圆柱尾迹影响下的湍流边界层近壁区的低速流体条带特征进行了观察和分析。结果表明：与没有尾迹扰动的情况相比，在圆柱下游低速条带的平均展向间距减小，圆柱距离壁面较近时这种减小更显著，可达到22％；而当圆柱离壁面较远时，尾迹对条带平均间距的影响减弱，其最大减小量的出现向下游推迟。另一方面，尾迹的作用并未使条带间距的统计分布性质发生改变，它们仍然符合对数正态规律。     

计算风工程中κ-ε模型的一类边界条件

给定适当的湍流来流边界条件是计算风工程中一个重要的问题.然而通常情况下,由于问题的复杂性,很难给定一个既满足边界层自保持要求又与试验相符的边界条件.本文首先从模型方程本身出发,推导出一类近似满足κ-ε模型自保持边界条件要求的湍动能表达式.此表达式与有关文献建议的湍流边界条件比较,可以考虑湍动能随高度的变化,因而更符合实测规律;然后基于风洞试验数据,定义了一组κ-ε模型中的模型常数;最后通过模拟风洞试验的算例,验证并分析了此类湍动能边界条件配合试验常数,在标准κ-ε模型中模拟均匀大气边界层的适用性.     

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

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

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

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

动态微涡流发生器对激波/边界层干扰的影响研究

本文基于OPENFOAM数值仿真平台,采用动态网格技术和湍流离散涡(DES)模型,研究了微涡流发生器以一定速度向下游移动时,激波/边界层干扰(SWBLI)流场特性的变化,重点关注干扰区域内的流向和展向的流场特性。来流马赫数为4,微涡流发生器向下游移动速度为0m/s,20m/s和40m/s。研究表明：当MVG向下游移动时,SWBLI区域的“弓”形高压区会演化成“双弓”形；入射激波形成高压区的压力明显降低,同时,入射激波和反射激波形成高压区的峰值位置均会向下游移动；流场下游 “双圆弧”状高压区的高度逐渐降低；SWBLI区域边界层的高度逐渐降低,同时边界层底部的速度也有所降低；随着MVG移动速度的增加,对SWBLI流场的控制效果更加明显；动态MVG对流场的控制是通过尾迹涡和波系结构实现的。     

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

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

网格湍流CAARC模型风洞实验

对大气湍流边界层的真实模拟是风工程风洞模拟实验所要满足的基本条件之一，平均速度剖面，湍流度剖面，积分尺度和风谱是反映大气湍流边界层流动的四个最基本的因素。本文首先研究了不同孔隙率的均匀网格湍流场的流动特性，得到湍流度，积分尺度以及风谱的变化规律，然后在变化基中某一个因素而保持另外三个因素不变的条件下，研究了处于均匀网格湍流场中的CAARC模型风荷载响应的变化，从而确定每个因素的独立变化对模型所受风荷的影响。实验结果表明湍流度和积分尺度对结构物所受风载都有很重要的影响，在风洞大气湍流边界层模拟中应该予以充分的考虑，否则将会引起风荷载实验结果的偏差。     

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

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

风洞模拟近地面大气边界层

本文通过大气边界层风洞的模拟实验，介绍一种在风洞中模拟近地面大气边界层的方法。根据大风地区的现场实测资料，对大气的风剖面及湍流度等有关参数，提出了一些特殊的要求。常规的模拟方法难以保证与实际风环境的相似，为此专门研制了一种在保证平均风剖面相似的前提下，可使气流产生大湍流度的装置，简称湍流度调节器。通过对模拟边界层的平均风速，湍流度及风速谱的测量分析，证明与实际风区大风环境的风结构相当一致。为解决实     

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

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

温度边界层对壁湍流多尺度相干结构的影响

用热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中壁面加热平板湍流边界层不同法向位置流向速度分量的时间序列信号, 与常温情况相比, 热对流加速了壁湍流中流体的动量、质量和能量的交换, 加快了湍流边界层的发展;用子波分析方法研究了壁面加热对壁湍流多尺度相干结构喷射和扫掠条件相位平均波形的影响, 显示热对流明显增大了缓冲层相干结构扫掠过程和喷射过程的强度, 是增强壁湍流中流体的动量、质量和能量交换的根本原因.      

封闭式薄膜屋盖的参数共振分析

膜结构的非线性风振响应是流固耦合的共同作用问题,当脉动风频率与结构频率满足一定关系时,就有可能出现大幅参数共振。本文首先假定风为无粘性势流,利用室内外空气速度势函数、屋面涡旋以及尾流涡旋来模拟屋面与风场之间的大气湍流边界层。针对封闭式薄膜结构受纵向脉动风荷载作用下,推导了屋面所受的气动力表达式,建立了薄膜屋盖参数振动的非线性力学模型。然后通过数值计算分析,证明了薄膜屋盖在达到临界颤振风速之前有参数共振的可能性。因而,在实际工程中应该考虑出现参数颤振的可能性,通过调整结构参数或预应力,以降低参数共振的频率范围。     

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.