绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

非定常来流压力下基于DMD方法的预冷器换热特性

采用大涡模拟方法对非定常来流压力条件下叉排管束预冷器换热特性进行了研究,同时运用动力学模态分解方法对流场主控流动结构进行了识别,探讨了来流压力周期性变化频率对预冷器内部流动、换热性能和熵产的影响。结果表明:来流压力变化频率对预冷器时均和瞬态换热性能影响均不显著,但当来流压力变化频率增大至流场固有频率950 Hz时,流场发生共振,换热性能发生剧烈振荡;管束壁面剪切层运动和绕流脱落涡结构为主控流动结构,其时空演化过程对瞬时换热性能起决定作用;当流场发生共振时,剪切层的生长和演化与来流速度的脉动密切相关,前排管束的绕流涡脱落周期与来流压力/速度变化周期一致,而壁面剪切层的生长周期则为来流压力/速度变化周期的两倍。此外,叉排管束流场的换热熵产决定于主控流动结构,其时空演化特征与主控流动结构演化规律完全一致。     

阀体后90°圆形弯管内部流场PIV测量及POD分析

利用激光粒子图像测速技术对阀体后90°圆截面弯管的内部流场进行测量,获取了弯管流场在不同来流流速及观测截面上的高分辨率瞬态速度场数据。分析了图像相关计算所获得的瞬态及时均速度场,对流场进行本征正交分解,并进一步分析重构脉动场中的涡结构特性。结果表明：流体流经蝶阀在弯管中形成卡门涡街,涡结构信息包含在本征正交分解低阶模态中,脉动场的模态重构可以还原涡的流动特性。     

三维钝体结构风压分布与绕流特性的大涡模拟研究

针对风场作用下三维钝体结构的绕流特性与空间风压分布特性,采用大涡模型对其进行了数值模拟研究,以平均风速和脉动风速作为入口边界条件,在稳态分析的基础上对9种不同来流方向的工况进行了瞬态动力分析,再现了钝体绕流的冲撞、分离、重附着等现象,揭示了钝体绕流主、次涡循环区的形成、发展机理,提出了空间风压场分布的一般规律。研究表明:钝体绕流对流场风速剖面影响的强弱程度依次为尾流区、屋顶区、来流区,且对中心区域流场风速剖面的影响要大于两侧流场,并由中心区域向两侧逐渐减小;来流区钝体绕流特性与空间风压分布特性较屋顶区和尾流区稳定,尾流区涡旋的非对称生成、脱落对钝体左右两侧流场的冲撞、分离、重附着等具有重要的影响;来流方向与钝体呈45°或135°夹角时,对尾流风场的影响范围最小,但随着来流方向逐渐垂直于钝体,尾流影响区域也逐渐增大。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

低亚声速和跨声速矩形柱绕流的大涡模拟研究

采用基于非结构网格的有限体/有限元混合格式和大涡模拟的方法求解可压缩Navier-Stokes方程,研究了不同长宽比矩形柱低亚声速和跨声速绕流的流动特性。在雷诺数为22000时,对来流马赫数等于0.1和0.75,截面长宽比分别为1∶1、2∶1、3∶1和4∶1的矩形柱绕流进行了大涡模拟,以研究长宽比和压缩性对矩形柱绕流流场的影响。马赫数为0.1时,Strouhal数随着长宽比的增大先降低再增大然后再降低;长宽比为3∶1和4∶1时会有流动的再附产生;柱体上表面的三维特性在长宽比大时更明显。马赫数为0.75时,Strouhal数随着长宽比的增大逐渐减小;湍流脉动和涡脱落受到抑制;方柱的近尾迹区域,有两种形成机制不同的局部超声速区。     

覆盖多孔介质的圆柱尾迹实验研究

利用风洞试验研究了多孔介质对圆柱绕流的影响,通过热线、烟线和粒子图像测速对有/无覆盖多孔介质柱体尾迹的频率特征、平均流场和瞬时流场进行了测量和分析,进而探究多孔介质与流场之间的作用机理。结果表明:引入多孔介质后,尾迹的大尺度涡脱落结构得到有效抑制,圆柱涡脱落频率降低;多孔介质减弱了近尾迹场的速度功率谱密度峰值,增大了远场的速度功率谱密度峰值,涡脱落位置延后;多孔介质减小了尾迹区的速度,拓宽了尾迹的宽度,使得两侧剪切层变得细长,削弱了它们之间的相互作用,减缓了两侧涡相互作用的速率;同时,多孔介质减小了尾迹区的速度脉动以及雷诺应力,雷诺应力和平均涡量的峰值更加远离圆柱,减弱了剪切层的不稳定性,增大了涡形成长度。     

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。     

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

直升机雷达天线非定常气动特性数值模拟

基于Realizable κ-ε湍流模型对直升机飞行过程中雷达天线的非定常气动特性进行了CFD数值模拟分析,并对全尺寸和缩比雷达天线进行了不同雷诺数下的绕流计算,得到了雷达天线的非定常流场分布和气动载荷,通过处理得到了雷达天线脱落涡的频率和对应的斯脱拉哈数St.计算结果与风洞试验结果进行了比较分析,给出了基本分析结论.     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

圆自由射流多尺度涡结构迁移特性的实验研究

在流场背景湍流度比较大或流场中涡结构的尺度范围比较宽的情况下，为了研究湍流场中不同尺度涡结构的迁移速度，用IFA300双通道恒温热线风速计和两个热线探针测量圆自由射流中心线上距离一定的空间两点同一时刻的流向速度分量，用子波分析提取了每一个探针输出的不同尺度的脉动速度信号，计算了两个探针输出的同一尺度脉动速度信号的互相关函数，利用互相关函数达到最大值对应的延迟时间研究了湍流多尺度涡结构的迁移特性。实验结果表明：不同尺度涡结构具有不同的迁移速度，小尺度涡的迁移速度较大，大尺度涡的迁移速度较小。其中，一部分尺度的涡结构是湍流中的主要结构，该尺度涡结构对湍流的贡献也是主要的，其迁移速度与实验测得的湍流平均速度较为接近。     

180°矩形弯管流场的LDV测量

采用激光多普勒测速仪（简称LDV）对180°矩形弯管内流场进行了测量,得到时均速度、湍流强度等数据。除靠近内壁r^＋=0．1位置,弯管纵截面上的切向速度沿轴向基本不变,但靠近弯管上下壁面的切向速度逐渐减小直至为零。在弯管的主流区域,0°～60°之间的纵截面上,内侧切向速度增大,外侧切向速度减小;60°～180°之间的纵截面上,内侧切向速度减小,外侧切向速度增大。在整个弯管段内,内侧切向速度总是大于外侧的切向速度。由于受到边界层分离和二次流的影响,90°～180°纵截面上r^＊ =0．1位置的切向速度产生明显的变化。轴向速度值远小于切向速度值,并且沿轴向变化不大。轴向速度的正、负之分,说明了二次流的存在,并且二次流的旋转中心从外壁向内壁移动。切向和轴向湍流强度的数量级一样,基本在0．1V。左右。切向湍流度在150°～180°纵截面r^＊ =0．1位置的变化很大;但是轴向湍流强度分布比较平稳,其值沿轴向和径向变化不大。     

Statistical Analysis of Velocity Measuring Systems Employing the Photon Correlation Technique

Velocities of fluid flows and solid objects can be measured by performing photon count correlations on laser light scattered from the measuring region, the most common optical configuration being the laser Doppler system. A statistical analysis for the counting process is presented for both the constant velocity case and for turbulent fluid flows with formulations for specific laser Doppler configurations. For turbulent flow, relationships are derived for the shape of the count correlation in terms of the mean velocity and rms turbulent intensity, taking into account the shape of the turbulence correlation curve. Experimental results recorded for both a constant velocity solid object and for an unseeded turbulent air flow in a wind tunnel are presented, and, in the latter case, comparisons are made with hot-wire measurements.     

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

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

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

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

密相、湍流、可压、气固两相流理论及在管道中的应用

就密相、湍流、可压、气固两相流的基本方程用双流体模型进行了较说尽地归纳。与单相流体相比，在两相流的基本方程中，多了一个变量（体积浓度）及两相间的耦合作用，并用CDIPSA方法对可压密相上升管道湍流流动进行了计算，结果表明：随颗粒直升的增加，使得单个颗粒脉动速度减小；但是气体及颗粒的速度差却随之增大。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

旋流器后火焰流场的试验研究

 本文给出了用LDA测量旋流器后火焰流场的试验研究结果。试验表明,中心开孔的平底碟形火焰稳定器后方形成双回流区,湍流脉动的均方根速度在回流区边界区域上发生突变,且湍流脉动速度的概率密度分布图上存在双峰;用热电偶测出的温度分布表明内外两支火焰在双回流区的内外边界附近开始,并向下游伸展。试验结果为应用旋流器燃烧室性能改进提供了依据,也为发展带回流的湍流反应流数值计算方法提供了验证依据。     

带有横向扩展出口的气膜冷却孔的紊流流场研究

报道了用热线风速仪和流动可视化技术对带有横向扩展型孔的燃气轮机气膜冷却叶片紊流流场进行的研究。试验了三种射流 /主流质量流量比 ,即 M分别等于 0 .5 ,0 .89,1 .5的情况。测量了主流速度、紊流强度以及雷诺剪应力。通过分析沿主流方向和展向方向各紊流量的空间分布 ,可以得知横向扩展型孔对燃气轮机气膜冷却叶片紊流流场的改善状况。研制了发生干冰气雾的实验装置 ,并将其应用于观察射流与主流复杂的相互作用。流动可视化结果与热线风速仪测量结果相一致