Ma数对后掠机翼流动非线性失稳的影响

横流非线性失稳是导致后掠机翼流动湍流转捩的主要因素,而目前大部分的研究都集中在不可压缩流动.首先,本文采用非线性抛物化扰动理论分析了来流Ma=0.1,0.7,1.5,2.0条件下,后掠机翼流动的非线性失稳过程.计算结果表明来流Ma数激励了横流涡的发展,增加了饱和横流涡的流向涡量,并且来流Ma数增加,使得饱和横流涡趋近壁面.其次,通过对饱和横流涡二次失稳的分析,本文计算结果表明来流Ma数增强了饱和横流涡的二次失稳,并且在亚声速来流条件下,二次失稳"z"模态占优,而在超音速来流条件下,二次失稳"y","z"模态同时占优.     

可压缩钝楔边界层转捩到湍流的直接数值模拟

直接数值模拟了来流Mach数为6的钝楔边界层在特定扰动(音速点附近壁面吹吸)下的转捩到湍流的整个过程. 分析了平均速度剖面, 脉动速度均方根及剪切应力等统计量, 并与不可压理论及实验结果进行了比较. 展示了转捩过程中的涡结构并分析了压力梯度对转捩的影响.     

槽道流转捩中发卡涡演化与波增长的关系

通过对槽道流常规转捩的直接数值模拟,研究了转捩突变前后,发卡涡的演化与T-S波增长的对应关系,重点对占流场绝大多数、转捩前期增长缓慢的T-S波的急速增长在发卡涡演化中所处的阶段进行了分析.研究发现,在发卡涡头部形成阶段,流场中典型的未增长起来的波一直保持着缓慢的增长趋势;当发卡涡头部产生分离时,这些波的实际增长率开始爆发式增加,并在短时间内呈数量级增长,急速增长的过程一直延续到发卡涡头部混乱之后.     

边界层转捩过程中环状涡和尖峰结构的演化

采用高精度直接数值模拟方法, 研究在边界层转捩过程后阶段中的环状涡和尖峰结构的非线性演化及其对流动的影响. 以平板可压缩流边界层转捩为例, 精细的数值模拟清楚地展示了环状涡等典型的涡结构及其上喷和下扫的流动现象, 证实了边界层转捩过程中尖峰结构的形成与环状涡紧密相关, 特别是, 对于最近在实验中观测到的正尖峰结构, 文中的数值结果也发现同样的结构, 并进行了相应的机理分析.     

串列双方柱的风压特性及其流场机理

为了研究干扰条件下方柱的风压特性及其流场机理,以串列双方柱为研究对象,采用大涡模拟方法,在雷诺数Re=8×104、间距比P/B=1.1～5的条件下,研究了两个方柱的风压系数、气动力系数、风压非高斯特性、风压相关性随间距比的变化规律,重点探讨了双方柱流场特性及其与风压非高斯特性的内在联系.研究结果表明:随着间距比的增大,串列双方柱依次表现为3种流态,即单一钝体、剪切层再附和双涡脱流态,风压特性与流态密切相关.风压的非高斯特性和风压相关性随流态变化呈现为3种类型:在单一钝体流态下,柱间回流区附近的表面风压呈现明显的非高斯特性且风压相关性较强;在剪切层再附流态下,方柱尾流的涡脱强度低,风压相关性弱,但风压非高斯区域大;在双涡脱流态下,受上游方柱尾流旋涡的作用,方柱侧风面的风压相关性较强,下游方柱的侧风面和背风面出现大范围的风压非高斯区域.     

边界层流动中湍流斑形成机制的数值研究

边界层流动中层流向湍流的转捩是一个古老但仍没有解决的问题．虽然对转捩过程的探索已付出了艰苦的努力,但依然有一个重要的物理过程还没有弄清楚,即转捩过程中湍流斑形成的原因以及湍流斑的运动特征是什么,这些问题正有待于人们做进一步的研究．文中提出以壁面局部脉冲形式的小扰动来模拟边界层流动中湍流斑的初始扰动场,采用紧致有限差分的方法,数值模拟了边界层流动中湍流斑的生成和演化规律;结果显示,它们在好多方面反映出湍流斑的基本特征．     

边界层流动中湍流斑形成机制的数值研究

边界层流动中层流向湍流的转捩是一个古老但仍没有解决的问题.虽然对转捩过程的探索已付出了艰苦的努力,但依然有一个重要的物理过程还没有弄清楚,即转捩过程中湍流斑形成的原因以及湍流斑的运动特征是什么,这些问题正有待于人们做进一步的研究.文中提出以壁面局部脉冲形式的小扰动来模拟边界层流动中湍流斑的初始扰动场,采用紧致有限差分的方法,数值模拟了边界层流动中湍流斑的生成和演化规律;结果显示,它们在好多方面反映出湍流斑的基本特征.     

平面槽道流中层流-湍流转捩的“breakdown”过程的内在机理

层流-湍流的转捩是一个古老但依然没有解决的问题. 尽管对转捩的研究付出了巨大的努力, 但仍然有一个重要的问题还没有涉及到, 那就是最终导致转捩的breakdown的内在机理是什么? 传统的看法是转捩从扰动的放大开始. 由于非线性的作用, 当扰动增长时, 高次谐波将产生, 流动变得越来越复杂, 最终导致湍流. 这一过程似乎非常清楚, 但却漏掉了一个环节, 即转捩的breakdown过程究竟发生了什么. 我们对直接数值模拟的结果进行了分析, 分析的结果表明, 平均流剖面稳定性特性的改变在breakdown过程中起了关键性的作用.     

超音速平板边界层转捩中层流突变为湍流的机理——时间模式

采用时间模式对来流Mach数为4.5的平板边界层转捩过程做了直接数值模拟, 然后对结果进行分析. 发现在层流-湍流转捩的breakdown过程中, 层流剖面得以快速转变为湍流剖面的机理在于扰动增长后, 对平均剖面进行的修正导致了其稳定性特征的显著变化. 虽然在层流下第二模态T-S波更不稳定, 但在层流突变为湍流的过程中, 第一模态不稳定波却占主导作用.     

全风向角下小间距双方柱绕流的大涡模拟

采用大涡模拟方法研究了间距比（P/B）为1.5(P为柱心间距,B为方柱边长)、风向角α为0°～90°等条件下双方柱在均匀来流作用下的气动力、流态划分、表面风压和流场特性。研究发现：小间距双方柱流动干扰效应显著,下游柱平均气动力随风向角的变化规律与上游柱和单方柱差异较大,且可能受更大的升力绝对值。将小间距比双方柱绕流分为前角分离流态（α=0°～10°）、分离泡流态（α=20°～30°）、附着流流态（α=40°～60°）及间隙侧分离泡流态（α=70°～90°）四种模式。附着流流态和间隙侧分离泡流态的间隙区出现较强的负压。尾流负压区的强度随风向角先增强后减弱,在α=50°附近达到峰值。前角分离流态时双柱具有一个整体的尾流回流区,而在其他流态下上下游方柱均有独立的回流区。     

小攻角高超音速尖锥边界层的转捩预测及e^N法的改进

研究小攻角高超音速尖锥边界层的转捩预测问题．来流马赫数为6,锥的半锥角为5°,攻角为1°．采用传统的e^N方法得到的转捩位置很不令人满意．原因在于,在传统的做法中,只考虑了增长中的扰动波,而事实上对于不同的子午面,扰动在开始增长前的衰减过程是各不相同的．根据我们以前的一些工作,对e^N方法进行了新的诠释,并且提出实质上的改进．不仅考虑了扰动的增长过程还考虑了扰动的衰减过程．由线性理论计算得到的扰动从初始幅值增长到1％的位置被看作是转捩位置．这样得到的新的转捩预测结果相当令人满意．此外,对于基本流的计算,研究表明边界层方程可以用于小攻角的情况,其计算量远远小于直接数值模拟．     

小攻角高超音速尖锥边界层的转捩预测及e~N法的改进

研究小攻角高超音速尖锥边界层的转捩预测问题.来流马赫数为6,锥的半锥角为5°,攻角为1°.采用传统的eN方法得到的转捩位置很不令人满意.原因在于,在传统的做法中,只考虑了增长中的扰动波,而事实上对于不同的子午面,扰动在开始增长前的衰减过程是各不相同的.根据我们以前的一些工作,对eN方法进行了新的诠释,并且提出实质上的改进.不仅考虑了扰动的增长过程还考虑了扰动的衰减过程.由线性理论计算得到的扰动从初始幅值增长到1%的位置被看作是转捩位置.这样得到的新的转捩预测结果相当令人满意.此外,对于基本流的计算,研究表明边界层方程可以用于小攻角的情况,其计算量远远小于直接数值模拟.     

不同半径比Taylor-Couette湍流的直接数值模拟研究

采用谱方法求解柱坐标系下三维不可压缩流Navier-Stokes方程,直接数值模拟了不同半径比下的同心旋转圆筒间的Taylor-Couette(TC)湍流问题.由于TC湍流是由大尺度Taylor涡和湍流随机脉动的叠加而成,采用周向平均成功识别了TC湍流中的大尺度Taylor涡,并将其诱导的脉动流动与湍流随机脉动分离开来,对比分析了不同半径比下大尺度Taylor涡对脉动强度和湍动能的贡献;同时,通过计算雷诺应力输运方程,研究了半径比对湍动能生成、耗散和再分配等动力学特性的影响.计算结果表明,在小半径比的宽槽TC湍流中,Taylor涡诱导的脉动流相对较弱,湍流随机脉动更为强烈,其对湍流统计特性的影响占优;在大半径比的窄槽TC湍流中,流体脉动特性主要源于Taylor涡的贡献;另外,随着半径比的增大,近外壁面附近流体沿周向剪切作用增强,流场表现出类似于平板Couette流的流动特性.     

黏性流中的链环涡管演化与螺旋度分析

本文通过直接数值模拟研究了链环涡管在不可压缩黏性流中的演化过程.在初始时刻的链环涡管由两个变形的涡环组合而成,其螺旋度为依赖于涡轴参数方程的解析表达式,进而可利用该初始流场进行涡管演化研究与螺旋度分析.发现当初始链环涡管的涡量方向具有相同手性时,链环涡管和环面纽结涡管的演化具有类似的涡动力学过程;而当它们具有相反手性时,涡管间的强涡量梯度会使两涡环在短时间内产生剧烈的涡重联,从而导致涡环由快速的尺度级串过程达到类湍流状态.     

叉排玻璃球(20mm)多孔介质内流态特征的实验研究

利用叉排20 mm玻璃球床、折射率匹配技术、PIV技术截取流场中心图进行实验。用Tecplot软件得出Re=4.7至Re=857.8的流场、流线和等涡量线图。根据流场、流线及等涡量线的密集程度判断:流场在Re=4.7时,流动缓慢,从Re=23.3至Re=128.4,流动强度递增,形成四边形涡波区域,Re=128.4至Re=285.9流动逐渐增强;Re=285.9至Re=857.8,流动稳定;流线在Re=4.7至Re=81.7时,流线流动缓慢递增,在Re=128.4至Re=245.1,流线流动逐渐增强,形成海螺状图形,Re=326.8以后流动稳定;本实验结果对一般流动流态和多孔介质内流态演变研究都有促进作用,为流动区域划分研究起了重要参考作用。     

湍流边界层雷诺应力和壁脉动压强相位滞后分析

在风洞平板湍流边界层外层引入圆柱尾涡的周期性扰动,在尾流干扰下的湍流边界层内进行测量,对相位滞后关系进行了实验研究,给出了湍流边界层内的雷诺应力和壁面脉动压强相对于边界层内的大尺度结构变形率的相位滞后沿壁面法向的变化趋势,分析结果表明剪切湍流涡粘模式中的涡粘系数应该是随时均流梯度变化的复数形式.     

边界层人工转捩的大涡模拟及壁面脉动压力研究

为了验证边界层人工转捩技术实现船模边界层壁面脉动压力载荷相似的可行性,采用大涡模拟以及功率谱估计获得平板湍流边界层壁面脉动压力功率谱,并与试验值、半经验模型值进行了对比,验证了数值模拟的可靠性。采用大涡模拟获得了绊线上下游壁面脉动压力的均方值、自功率谱和波数频率谱;以脉动压力均方值、自功率谱特征为判据,对比了方形、锯齿形绊线的转捩效果。结果表明:绊线高度雷诺数大于吉宾斯雷诺数;两种绊线均可实现边界层转捩,且壁面脉动压力的功率谱特征、谱级相近;与同流速平板湍流边界层相比,即使在当地雷诺数较低时,绊线下游边界层壁面脉动压力自功率谱平台区仍然高出10 dB左右。设计的绊线可用于减弱船模边界层壁面脉动压力的尺度效应,并有助于实现壁面脉动压力载荷相似。     

平板边界层湍流的数值分析

对于不可压缩平板边界层的时间模式的直接数值模拟(DNS)而言,在转捩后继续计算,就可以得到湍流直接数值模拟结果.分析计算结果发现,当湍流充分发展后,其平均流剖面、各种量的脉动均方根值和雷诺应力等沿平板法向的分布都具有相似性.不过,从转捩完成至湍流充分发展之间有一过渡过程,其间上述相似性不成立.这一结果为简化但又较准确地计算工程技术问题中的湍流提供了一种可能性.此外,还分析了其他一些湍流特征,如形状因子和位移排移厚度随时间的演化规律.为了观察到相干结构,须将展向的计算域缩小一半,以增加分辨率,研究结果表明,湍流边界层近壁区存在相干结构,其主要表现是准流向涡或涡对.     

振荡流底层拟序结构运动理论模式

湍流拟序结构是一种可检测的有序运动,在壁面湍流结构中往往表现出复杂的三维结构.针对对称振荡流(纯振荡流),建立了描述对称振荡流湍流拟序结构的理论模式,并对其特性进行了研究.探讨了振荡流拟序结构能够保持的雷诺数范围、拟序扰动波数范围以及所能激发不稳定拟序波的波数频率范围,分析了振荡流激发的二维扰动波与三维扰动波的相互关系以及拟序波的对称性特征、拟序扰动流速、雷诺应力、附加涡量的时空分布等.结果表明,对称性振荡流拟序结构不仅与通常研究的雷诺数等因素有关,而且与振荡频率有很大的关系.     

硅油扭振阻尼器理论分析及实验研究

以涡－流函数法研究了硅油流体在扭振阻尼器内部的流动规律,在一定简化条件下,用有限元法对受叶片初始扰动下的流体流动性态进行了数值模拟,得到了流场规律．通过涡场强度的分析,揭示了有关阻尼机理,并利用频响函数法对典型阻尼减振系统进行了实验．结果表明,文中采用的方法在指导设计中是可行的．