低速高湍流度90°弯管流动数值模拟

分析了曲率对弯曲管道流动的影响 ,给出了曲率修正的双层 k-ε湍流模型 ,并数值模拟了非均匀来流低速高湍流度 90°弯管内流动。文中比较了曲率修正双层 k- ε湍流模型的计算结果与实验结果 ,以及双层 k- ε湍流模型和逆压力梯度修正双层 k- ε湍流模型的计算结果。通过比较 ,发现经过曲率修正后的双层 k- ε湍流模型更好地模拟了非均匀来流低速高湍流度 90°弯管内流动。同时比较了壁面参数的插值方法 ,发现在弯曲管道流动模拟中对壁面参数进行一阶插值比零阶插值能更好地模拟真实流场     

稳定分层湍流的PIV实验研究

在自行建立分层湍流实验装置的基础上,进行了一系列分层湍流的实验研究,运用激光诱导荧光(LIF)流动显示技术和PIV流场测量技术,定性研究了剪切和分层对湍流混合和湍流结构的影响.实验结果发现:剪切促进湍流混合,稳定分层抑制湍流混合,而且分层使得混合层内的湍流结构变得扁平和细长,同时还发现了混合层内涡旋的拉伸、合并、扭曲和变形等复杂的流动现象,特别是发现了其中的旋转涡对和反旋转涡对以及对湍流输运起重要作用的手指状结构.     

近壁湍流相干结构传播速度的研究

近壁湍流相干结构传播速度的问题,目前尚无理论结果。本文利用流动稳定性理论中的共振三波概念建立单个相干结构模型,采用紧致有限差分和Ｆｏｕｒｉｅｒ谱展开相结合的方法对近壁湍流相干结构进行了数值模拟,计算得到的相干结构流向传播速度与直接数值模拟的结果一致。     

基于大涡模拟的NACA0012翼型流动分离研究

本文针对翼型流动分离问题，利用求解RANS方程和大涡模拟数值计算方法，通过常用的三种湍流模型的对比：k-epsilon、SST k-omega和BSL Reynolds Stress，最终选取SSTk-omega湍流模型对NACA0012翼型进行定常流动和非定常流动数值计算，得到翼型的升力系数特性曲线和翼型上下翼面的速度和压力云图。研究结果表明，大涡模拟能够捕捉流动分离区域中主要涡结构，在一个周期的不同时间步能够观测到分离涡从附着、发展、脱落到再附着的完整演化过程，对翼型流动控制的研究具有一定的参考，为飞行器的增升减阻设计提供了理论支持。     

非正交曲线坐标湍流旋流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标系统，用坐标变换的方法处理锥形燃烧室内线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。紊流模型采用工程上常用的ｋ－ε模型。还运用贴坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提出了理论依据。     

透平叶片顶部间隙流动特性的实验和数值研究

利用粒子图像测速技术(PIV)捕捉透平叶片顶部泄漏流特征,并以此数据验证湍流模型和用商业软件CFX12.0进行的数值模拟方法.所研究的叶片为典型的GE-E3叶片,为了展示泄漏涡的生成和发展过程,用实验数据展示了3个不同截面的速度分布.数值计算中使用了混合网格生成技术及5种湍流模型.通过与实验数据的比较发现:RNG k-ε模型计算所得的泄漏涡与实验所拍摄的真实流动能较好地吻合.此模型和计算方法同样适用于研究叶顶射流对泄漏流的影响.计算结果显示:通过叶片顶部气膜孔射流产生的阻挡效应,最多能降低6.12%的主流泄漏.     

SST-DDES模型在大分离流动问题中的应用

采用基于SST k-ω湍流模型耦合延迟脱体涡模拟技术(Delayed detached-eddy simulation,DDES)的SSTDDES模型,数值求解Navier-Stokes方程。并针对低速流动问题,发展了基于隐式双时间步方法的低速预处理技术,用于数值模拟低速、高雷诺数湍流大分离流动问题。数值模拟了低速大迎角三角翼绕流及大迎角6∶1椭球绕流,观察到与物理现象相一致的旋涡特征,且得到的定量数值结果与实验数据相吻合。     

展向不同间距壁面粗糙元对近壁湍流拟序结构的影响

采用氢气泡流动显示技术,研究了壁面不同展向间距粗糙元对近壁湍流拟序结构的影响.实验中粗糙元布置在氢泡丝前,其间距按如下两种排列设置:粗糙元直径的0、1和3倍,等间距1、2和4cm;每种情况粗糙元直径均为2.3、4和6mm.基于平均流速和水力直径的流动雷诺数从14000到48000变化.实验观测了光滑壁面和粗糙元壁面的近壁湍流拟序结构条带特征.发现:对于两种不同排列,相同雷诺数和相同粗糙元间距条件下条带有量纲间距均随着粗糙元直径增大而减小,条带有量纲高度均随粗糙元直径的增大而增大;相同雷诺数和相同直径粗糙元条件下条带有量纲间距均随着展向间距增大而增大,条带有量纲高度均随展向间距的增大而减小.粗糙元壁面湍流条带有量纲间距均比光滑壁面小,条带有量纲高度均比光滑壁面大.研究结果对近壁湍流的流动控制具有重要指导意义.     

高强驻波声场中的湍流

反射端声压级高于140dB以后驻波管中的声流结构会发生明显变化。利用流动显示和LDV（Laser Doppler Velocimetry）技术对该问题进行了实验研究。流动显示结果表明，在160dB左右的高强驻波声场中随着旋涡结构的合并、破裂等现象的出现，声致流动有明显的湍流特性。LDV测量结果显示，随着驻波管中声压级的提高，一些测点的速度频谱明显加宽，出现高频速度脉动，而另一些点的速度频谱则没有显著变化。径向湍流强度测量结果显示，此时管中流动已是复杂的三维结构,而非轴对称的二维结构。     

后台阶流动的大涡模拟

利用大涡流模拟,对后台阶中的流动进行了数值计算,并和实验数据进行了比较,吻合较好,表明了所采取的亚网格湍流模型的合理性,但是需要改进。     

从壁面脉冲到湍流相干结构演化的数值模拟

发展了高效求解基于Navier-Stokes方程的扰动方程高阶紧致差分方法，研究了壁面强脉冲到湍流相干结构的复杂演化过程。该方法采用高精度、高分辨率的三维耦合中心差分格式，以及和高阶时间分裂法相结合因而更为有效，比一般谱方法更适用于较为广泛的流动区域和边界条件的湍流相干结构研究。中模拟了槽道流动中壁面强脉冲诱导的湍流相干结构的产生、发展和演化过程，显示了湍流相干结构的基本特征和发展规律．计算结果表明，所用壁面强脉冲模型比共振三波模型更有效．     

湍流边界层近壁区马蹄涡拟序结构的数值模拟与分析

采用大涡模拟对低雷诺数槽道湍流进行了数值模拟,进口边界条件采用了给定扰动波方法,经过短过渡段得到完全发展湍流.计算中采用了动态亚格子模型和预处理技术,进口采用的动量厚度雷诺数是670.大涡模拟对计算得到的平均速度剖面以及脉动分量与DNS解进行对比,验证了计算结果的可靠性.计算结果显示了湍流场中马蹄涡的形成及演化过程,其中包括形成单腿马蹄涡和形成对称涡腿的马蹄涡,同时发现流场中存在由亚谐波引起的拟序结构的交错现象,并在此基础上分析了湍流边界层近壁区马蹄涡结构的演化.     

NUMERICAL SIMULATION OF FLOW PATTERNS IN PLANAR JETS

INTRODUCTIONTurbulent jet flow is the most widely usedflow type in engineering applications especiallyfor combustion engines which provide most ener-gy sources for industrial production and trans-portation,and a large proportion of energy trans-portation and conversion processes is dominatedby transient evolution of large- scale structures,or so called coherent structures,in turbulentjets.As the flow control and optimization pro-cesses are closely related to the understanding ofdetailed in…     

低雷诺数下进气道异常起动现象及其影响因素探析

结合激波风洞实验和数值模拟分析,对一种二元混压式进气道在实验中低单位雷诺数下反而呈现出自起动特征的异常现象进行了研究。根据激波风洞的反复实验观察,表明随着来流单位雷诺数的降低,在继进气道进入不起动状态之后又会重新出现自起动特征的异常起动现象。该结果与层流模型计算得到的流场结构相符,而与湍流模拟结果差异显著;分析表明,层流情况下,由于分离区向前体压缩面大范围地延伸,缓解了进气道入口的逆压梯度,从而在喉道处可以形成主体为超声速的通畅流道,而湍流情况下,进气道入口处激波/边界层干扰形成过分集中的分离泡则呈现明显的壅塞状态;尽管层流情况下进气道流场结构呈现出较为通畅的类似起动的特征,但其流量系数仍明显低于湍流的情况。因此,实验上所观察到的这种异常起动现象严格地说并不属于真正意义上的起动状态。     

局部粗糙边界层流中二维T-S波的非线性演化

采用局部喷、吸、喷与吸组合的结构来模拟局部粗糙壁面,通过数值计算获得了稳定的三维边界层基本流的数值解.在此基础上,研究了最不稳定的和稳定的二维扰动T-S波在时、空间演化过程中非线性问题;探讨了局部粗糙的形式、强度大小及分布结构对二维T-S波的非线性演化的影响及其对增长速率贡献的大小.数值结果表明,三维局部粗糙壁面诱导产生的基本流是扰动波得以增长的关键因素;最不稳定的二维扰动T-S波在非线性演化过程中,由于非线性相互作用的不断增强,逐渐诱导激发出三维扰动波及高次谐波,其三维扰动波的流向波数和频率与二维扰动波的流向波数和频率大致相同.此外,展向速度的大小激励着二维扰动波的快速增长,随后发生初次失稳、流向涡形成、正负相间剪切层逐渐增强、诱导产生三维不稳定扰动波及高阶失稳等机制.上述结论与实验结果相吻合.     

强湍流下点火及火焰传播机理研究进展

可燃预混气点火过程研究是发动机燃烧领域最重要的课题之一。当前电火花强制点火广泛应用于各类发动机燃烧室中,其点火过程具有很强的瞬变性,从电火花产生到火焰完全形成的整个过程中,多种复杂因素联合影响点火火核的产生和发展。目前发动机高压、强湍流工况下的点火和火焰传播机理认识还不完善,亟需研究的科学问题是湍流和化学反应相互作用对点火和火焰传播的影响机制及其建模,包括湍流对点火的促进/抑制机制,湍流对火焰传播和火焰整体发展的影响规律,燃烧释热和火焰面不稳定性对湍流脉动速度(即火焰产生的湍流)的影响机制和对火焰传播速度的增强机制及由此导致的层流燃烧自加速转捩为湍流燃烧的理论,燃烧过程对标量通量输运(即反向或压力驱动输运)的影响机制。本文对强湍流下点火及火焰传播理论、实验和数值模拟方面的研究进展进行综述。     

非正交曲线坐标湍流族流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴你坐标系统，用坐标变换的方法处理锥形燃烧室曲线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。亲流模型采用工程上常用的ｋ－ε模型。还运用贴体坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提供了理论依据。     

采用CARS试验技术与UFPV数值方法研究航空发动机燃烧室

在自主开发的软件平台上,采用基于URANS的方法计算航空发动机燃烧室的三维两相燃烧流动,考虑了液态燃油从液膜-液滴-燃气-燃烧的完整物理化学过程。其中,颗粒相采用LISA一次破碎模型,KH-RT二次破碎模型和标准的蒸发模型,湍流燃烧模型采用可以考虑非稳态燃烧特性的非稳态火焰面/反应进度变量方法,得到了航空发动机燃烧室中温度、组分浓度和燃油液滴的颗粒直径分布规律。同时,采用CARS光学手段测量燃烧室主燃区的温度分布,并将数值计算结果与光学试验测量值进行比较,数值计算结果和试验值吻合较好,数值计算误差小于7.3%。说明了本文的数值计算方法和UFPV方法在计算航空发动机燃烧室的两相燃烧流动时具有较高的精度。     

封闭空间中不同压力下火焰过孔板加速机理研究

爆震或超级爆震发生时总会伴随着湍流火焰-冲击波相互作用,对其开展研究是揭示爆震或超级爆震机理的关键,研究火焰加速产生压力波的过程是火焰-压力波相互作用研究的基础性前提。基于自主设计的定容燃烧弹和Converge三维数值模拟方法,对封闭空间中火焰过孔板加速机理及影响因素开展了研究,讨论了初始压力对火焰过孔板加速的影响。依据火焰传播形态与速度,将火焰过孔板加速过程分为3个阶段:层流火焰阶段、射流火焰阶段和湍流火焰阶段。通过分析火焰过孔板过程中的流场情况,发现在火焰未到达孔板前,孔板附近存在强射流,火焰受强射流的驱动而急剧加速;但当火焰穿过孔板之后,火焰锋面前的流场速度沿着远离火焰的方向而逐渐下降,说明开始由火焰驱动未燃气体运动。比较不同压力下的火焰过孔板过程,发现湍流火焰传播速度和缸压振荡均随着初始压力的提高而升高。     

流动环境中三维圆形垂直浮力射流特性实验研究

为了揭示流动环境中三维圆形垂直浮力射流的规律,运用MicroADV流速测量系统和温度测量系统对不同射流——来流流速比情况下的圆形垂直热水浮力射流进行了实验研究。在对试验所得的速度场、温度场资料进行整理的基础上,对流场的流动结构,射流中出现的分叉现象,射流轨迹线,湍动能等时均特性和紊动特性进行了分析。认为基于最大速度和最高温度的射流轨迹线是不同的,射流横断面上两高温区的夹角的范围同前人在空气中所作试验是一致的。表明试验结果是可信的,可以用于揭示圆形浮力射流的规律和数值模拟的验证。