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

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

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

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

壁湍流调制雷诺应力和拟序结构应变率空间形态的层析PIV实验研究

采用高时间分辨率的层析PIV技术,测量了水洞壁湍流三维速度分量空间分布的时间序列,应用局部平均速度结构函数概念对流向速度信号进行多尺度分解,以流向速度分量的局部平均速度结构函数过零点作为特征量检测壁湍流中拟序结构猝发的喷射和扫掠过程,应用空间相位平均技术提取拟序结构猝发的喷射和扫掠过程各速度分量、涡量分量、拟序结构速度应变率分量以及调制雷诺应力分量的空间相位平均拓扑形态。为了引入平衡态涡粘模型假设模拟调制雷诺应力,研究了拟序结构猝发过程中调制雷诺应力分量和速度变形率分量的空间分布形态,发现两者之间的空间相位分布不一致。由于存在时空相位不同步性,说明需要考虑大尺度拟序涡结构引起动量传递的时空弛豫效应。应用经典的线性平衡态下的Boussinesq涡粘模型不能准确地描述壁湍流拟序结构动量传递非平衡现象的物理机理。对于壁湍流拟序结构动力学方程中调制雷诺应力的模拟,应采用包含时空相位信息的复涡粘张量模型。由于雷诺应力与速度变形率的时空不同步性,对非平衡非局部湍流场的数值模拟提出了一个挑战性的问题,建议采用包含相位信息的复涡粘张量模型来模拟雷诺应力张量,从而更加符合雷诺应力演化的物理机理,这一模型有可能成为一个很有发展前景的封闭模型,从而更加准确地预测工业领域中广泛存在的非平衡湍流。     

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

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

变湍流度时翼型边界层及近场尾流的法向湍流特性初步试验 …

阐述了湍流度对边界层及近场尾流时均特性和边界层转捩位置的影响，讨论了流向和法向湍流特性分布规律，并着重分析了湍流度法法向湍流特性的影响，结果表明，法向湍流强度及雷诺正应力沿流向分段性更为明显，且法向与流向湍流参量之比的特征值有一定规律性。     

超疏水壁面湍流边界层减阻机理的TRPIV实验

利用高时间分辨率粒子图像测速（TRPIV）技术,开展超疏水壁面材料湍流边界层减阻机理的实验研究.在循环水槽中,对超疏水壁面和亲水壁面湍流边界层瞬时速度矢量场的时间序列进行了实验测量.得到了同一来流速度（0.17m/s）下超疏水壁面和亲水壁面湍流边界层的平均速度、湍流度及雷诺切应力沿法向的分布规律.提出了空间多尺度局部平均涡量的概念,并以此为特征量检测壁湍流发卡涡展向涡头的中心位置.用条件采样及空间相位平均技术提取了不同法向位置发卡涡展向涡头周围流向脉动速度和流线的空间拓扑,对发卡涡展向涡头的俯仰角进行了对比,并从鞍点-焦点动力系统的角度分析了发卡涡展向涡头附近的流线拓扑特征.研究表明：雷诺数约为13500时,相比亲水壁面,超疏水壁面实现了10.1%的减阻.超疏水壁面平均速度明显增大,雷诺切应力减小,流向湍流度减弱,发卡涡展向涡头俯仰角较小,近壁区相干结构的发展受到抑制.     

槽道湍流在确定分布的展向电磁力调制下的壁面减阻机理

为了研究槽道湍流中确定分布的展向电磁力的流动控制与减阻问题,采用直接数值模拟的方法,对槽道湍流的确定分布展向电磁力控制后的流场中流向脉动速度、法向脉动速度以及雷诺应力分布的变化规律进行了研究。研究结果表明,槽道湍流流场存在大量分布紊乱的准流向涡和发夹涡结构,经过确定分布的展向电磁力调制之后,该流场的发夹涡结构基本消失,仅存在规律分布的准流向涡结构。控制后的流场相比于湍流流场,其流向脉动速度和法向脉动速度均受到了束缚作用,尤其是近壁区域更为明显。整个控制过程的实质是雷诺应力离散点向原点集中的过程,这个过程导致了流场雷诺应力和壁面阻力的下降。     

沟槽壁湍流多尺度相干结构实验研究

应用热线测速技术,对沟槽壁面平板湍流边界层的减阻机理进行了实验研究.测量了风洞中并排放置的沟槽壁面平板及光滑壁面平板湍流边界层在不同雷诺数下不同法向位置的瞬时流向、法向速度分量的时间序列信号.运用流向速度分量信号的多尺度子波系数辨识壁湍流多尺度相干结构,用条件采样和相位平均技术提取了壁湍流多尺度相干结构喷射和扫掠时流向速度分量、法向速度分量和雷诺应力分量的相位平均波形.分析了沟槽壁面平板及光滑壁面平板湍流边界层中多尺度相干结构的持续时间、条件相位平均波形等特征.对沟槽壁面平板及光滑壁面平板湍流边界层多尺度相干结构的多种统计平均特征进行了比较,从壁湍流多尺度相干结构控制的角度研究了沟槽壁面平板湍流边界层的减阻机理.      

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

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

周期性扰动对平板湍流边界层统计性质的影响

用IFA300热线风速仪和X形二分量热线探针精细测量了平板湍流边界层在施加不同频率的周期性吹吸扰动前后不同法向位置的瞬时流向、法向速度分量的时间序列信号.在平板湍流边界层的固壁表面沿展向切割一条5 mm宽的窄缝, 在平板底部用扬声器从壁面向平板湍流边界层内施加不同频率的周期性吹吸扰动, 研究不同频率的周期性扰动在平板湍流边界层内发展演化的规律及对平板湍流边界层统计性质的影响.      

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

Numerical study of transonic separation flow using an anisotropic turbulence model

A transonic turbulent separation flow in a converging-diverging transonic diffuser was studied, when there existed a separation bubble on the top wall of the diffuser triggered by strong shock-wave-boundary-layer-interaction (SWBLI). To capture the essential behavior of this complex flow, the current study utilized an anisotropic turbulence model developed on the basis of a statistical partial average scheme. The first order moment of turbulent fluctuations, retained by a novel average scheme, and the turbulent length scale, can be determined from the momentum equations and mechanical energy equation of the fluctuation flow, respectively. The two physical quantities were readily used to construct the nonlinear anisotropic eddy viscosity tensor and to significantly improve the computational results. Comparisons between the computational results and experimental data were carried out for velocity profiles, pressure distribution, skin friction coefficient, Reynolds stress as well as streamline vectors distribution. Without using any empirical coefficients and wall functions, the numerical results were in good agreement with the available experimental data, further confirming that the nonlinear anisotropic eddy viscosity tensor is the decisive factor for the success of the computational results.     

基于大涡模拟分析气膜冷却的湍流场

用大涡模拟考察了单孔平板气膜冷却的湍流场,气膜孔沿流向倾斜30°,气膜出流的雷诺数为2600,吹风比为0.5.计算结果表明:①气膜孔内的流动分离增加了湍流动能,导致气膜冷却效率的降低,在设计中要尽量避免或减少流动分离;②与射流侧向扩展有关的涡黏性系数在气膜孔两侧存在峰值,而采用各向同性的湍流模型预测气膜冷却时,涡黏性系数的峰值出现在射流与主流的剪切区,因此会高估射流的垂直穿透而低估射流的侧向扩展;③可以用大涡模拟辅助建立各向异性的湍流模型,以便提高湍流模型的预测精度.      

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

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

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

旋转光滑直通道湍流流动二维热线实验

采用二维热线测量了旋转光滑直通道内不同流向位置上的平均速度和雷诺应力。实验结果表明:较高的当地旋转参数使得旋转对平均速度的影响区域扩大并且导致了前缘面附近湍流流动的再层流化。后缘面附近无量纲主流平均速度型严格按照旋转数顺序依次排列，且在对数律区呈现对数律分布。与此同时，所有无量纲雷诺应力分量在后缘面附近基本不受旋转影响。再层流化导致了前缘面附近无量纲主流平均速度型无法在对数律区维持对数律分布，且所有雷诺应力分量都随着旋转数和流向位置半径增大不断衰减；经过u-v象限分析，再层流化现象的直接原因被归于湍流脉动生成减少。      

有限体积时间推进法在跨音速扩压器湍流流场中的应用

本文提出了一种有限体积时间推进法,将其应用于均化N-S方程的求解,模拟跨音扩压器的湍流流场,并采用空间变时间步长和多重网格技术加快收敛速度。计算结果与有关文献的实验数据吻合良好。     

圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。