两种亚网格湍流模型的旋流扩散火焰大涡模拟

为研究旋流火焰结构,采用二阶矩亚网格(SOM-SGS)燃烧模型及Smagorinsky-Lilly和K方程亚网格湍流模型,对美国Sandia国家实验室测量的旋流火焰进行了大涡模拟,得到了与燃烧场速度、温度和温度脉动实测值相吻合的模拟结果。预报的瞬时温度分布云图与实际火焰的形状很相似。K方程亚网格湍流模型预报的瞬时温度分布比Smagorinsky-Lilly模型的预报结果更接近实际。燃烧火焰基本上位于回流区所在的位置,火焰在回流区被稳定。     

气固流动大涡模拟和两相湍流模型的评价

用基于气体Smagorinsky亚网格应力模型和颗粒动力论模型的双流体大涡模拟(LES)以及统一二阶矩两相湍流模型(USM-RANS),对轴对称突扩气固流动进行了数值模拟。结果表明:LES模拟的颗粒轴向平均速度和均方根脉动速度,以及USM-RANS计算得到的模拟结果都与实验结果吻合较好。LES得到的颗粒-气体轴向脉动速度关联值比USM-RANS模拟值更接近实验值,这表明USM-RANS模拟还有改进余地,并基本上得到了LES的验证。LES的瞬态模拟能显示气固湍流流动的各向异性两相湍流结构和颗粒弥散的发展过程,而USM-RANS无法实现。     

旋流稳焰瓦斯燃烧器湍流扩散火焰的数值模拟

为研究瓦斯燃烧器的火焰特性，采用κ-ε湍流模型、三气体扩散燃烧模型及综合辐射模型，对旋流稳焰瓦斯燃烧器湍流扩散火焰进行了数值模拟。研究结果表明，在旋流作用下，在火盆及其前端附近区域内形成了中心回流区，使燃气与助燃空气在此处发生剧烈掺混而被点燃，这有利于火焰稳定。在火焰前端仍存在一个高温尾流区，在选用和设计燃烧器时应充分考虑这一现象。燃气喷孔直径对火焰长度的影响非常显著，喷头顶部气孔位置和火盆锥口角对火焰长度影响较小。随着空气过剩系数和燃气中空气含量的增加，火焰长度明显减小。在实际应用中，应充分考虑燃烧器结构和操作条件的影响，以提高加热炉的效率和安全性。     

旋流稳焰瓦斯燃烧器湍流扩散火焰的数值模拟

为研究瓦斯燃烧器的火焰特性 ,采用k ε湍流模型、三气体扩散燃烧模型及综合辐射模型 ,对旋流稳焰瓦斯燃烧器湍流扩散火焰进行了数值模拟。研究结果表明 ,在旋流作用下 ,在火盆及其前端附近区域内形成了中心回流区 ,使燃气与助燃空气在此处发生剧烈掺混而被点燃 ,这有利于火焰稳定。在火焰前端仍存在一个高温尾流区 ,在选用和设计燃烧器时应充分考虑这一现象。燃气喷孔直径对火焰长度的影响非常显著 ,喷头顶部气孔位置和火盆锥口角对火焰长度影响较小。随着空气过剩系数和燃气中空气含量的增加 ,火焰长度明显减小。在实际应用中 ,应充分考虑燃烧器结构和操作条件的影响 ,以提高加热炉的效率和安全性。     

稳定分层湍流的大涡模拟

稳定分层湍流是大气和海洋流动中常见的自然现象, 是预测大气和海洋流动的重要物理机制. 我们用大涡模拟方法研究稳定分层流体湍流, 以期获得与实际大气和海洋中稳定分层湍流相符的性质. 我们在不同初始条件的流场和不同的流动特征参数（雷诺数Re和弗罗德数Fr）下用大涡模拟方法获得充分发展湍流场. 通过分析湍动能和湍流特征尺度随时间的演化考察稳定分层湍流演化过程; 并通过三维能谱、水平和垂直能谱来检验数值模拟是否和真实大气、海洋中实测结果一致. 研究结果表明, 大涡模拟可以用较少网格数获得和实测一致的稳定分层湍流性质. 不同的初始流场都能演化为湍流, 并且发现演化过程的后期总是内波能量占优势. 研究证实要获得和实测一致的稳定分层湍流的关键是特征参数必须满足一定条件： 特征雷诺数和弗罗德数不一定要和实际大气和海洋中的参数相等, 但是要求雷诺数足够大, 弗罗德数足够小（例如Re〉102, Fr〈0.1）; 特别重要的是组合参数ReFr2必须大于10.     

大涡模拟预报可靠性分析

针对明渠道湍流流动,研究了计算网格、计算域大小、亚格子（SGS）模型和雷诺数等因素对大涡模拟（LES）预报精度的影响。研究表明,即使简单外形的湍流流动,网格尺寸大小及匹配仍较大地影响大涡模拟计算结果,因此具体使用中应予以充分重视。     

离心泵叶轮内部清水湍流的动态大涡模拟

提出旋转坐标下的单相液体湍流的二阶双系数动态亚格子应力大涡模拟模型并加以验证。考虑亚格子应力的对称性和量纲一致性 ,在基于应变率张量和旋转率张量的不可约量及 Smagorinsky模型和亥姆霍兹定理的基础上 ,提出亚格子应力应表示为可解的应变率张量和旋转率张量的函数。在贴体坐标系下 ,利用交错网格系统和有限体积法离散湍流控制方程 ,采用 SIMPL EC方法求解方程。水泵叶轮中湍流的流速分布规律和压力分布等计算结果都与试验结果基本吻合 ,从而证明了此模型的可行性 ,以及计算方法和程序的可靠性     

180°方形弯管中湍流流动的大涡模拟

把大涡模拟应用于复杂几何形状的流域是一个重要研究课题。对１８０°方形弯管中的湍流流动作大涡模拟,以网格节点为中心的交错网格作为离散网格,用有限体积方法和中心差分进行空间离散,用Ａｄａｍｓ－Ｂａｓｈｆｏｒｔｈ格式进行时间离散,用简单的Ｓｍａｇｏｒｉｎｓｋｙ模型来作为次网格湍流模型。用多重网格方法求解压力方程以提高计算效率。数值模拟结果与实验结果取得了很好的一致。     

基于格子Boltzmann方程的大涡模拟对湍流时空关联性的研究

将格子Boltzmann方程和大涡模拟(LBE-LES)相结合,提出适应于格子Boltzmann方法(LBM)的涡黏性亚格子尺度模型,开展均匀各向同性湍流时空关联性的研究.采用D3Q19格式计算湍流的三维能谱、湍动能耗散率和其它高阶统计量,与实验和直接数值模拟结果的比较表明,该模型比传统涡黏模型有明显改进.考察了不同亚格子模型预测湍流频率波数能量谱的能力,结果表明,尺度涡产生的横扫作用是造成小尺度涡时间去关联的主要因素,不同波数的频率能量谱之间有一定的相似性,横扫速度是描述湍流频率波数能量谱的特征量.     

湍流横向射流大涡模拟及大尺度结构的涡环模型

采用大涡模拟方法数值模拟了流向椭圆喷嘴的湍流横向射流,重点研究其旋涡结构的产生、发展等动力学演化过程．结果表明：文献中所报道的横向射流基本涡结构,如反向旋转涡对、前缘涡、后缘涡、悬涡、肾涡、反肾涡等并非独立的物理实体,而是分别对应于新发现的横向射流基本涡结构——起始于喷嘴的三维拉伸涡环的局部结构．     

2.5MW煤粉氧燃烧旋流火焰数值计算

为了准确地数值预报煤粉旋流火焰特征和主要烟气成分,比较扩展的涡耗散模型在多湍流模型下的预测性和影响,针对IFRF(国际火焰研究基地)2.5MW煤粉氧燃烧和空气燃烧实验,利用组分输运方程,结合7步改进的总包反应,对炉内煤粉氧燃烧进行数值计算,并与实验结果和空气燃烧基况对比.结果表明:氧扩散率、发射率和比热容等物性参数修正后,湿循环、氧预混二次流的煤粉氧燃烧扩散火焰温度场与空气燃烧总体一致,火焰稳定且为Ⅱ型火焰结构;EDC模型(涡耗散概念模型)对组分预测更为准确,尤在可实现k-ε模型下的缓慢反应CO生成和IRZ(中心内回流区)高温预测;而FRED模型(有限速率涡耗散模型)在多湍流模型下对炉内温度、组分、火焰结构预测亦较准确.     

Large-eddy simulation of the compressible flow past a tabbed cylinder

The compressible flow past a tabbed cylinder has been studied numerically using large-eddy simulation for a free-stream Mach number M∞=0.75 and a Reynolds number based on the diameter Re=2×105. Because of the passive control of the flow past a tabbed cylinder, the mean drag coefficient of tabbed cylinder is less than that of a corresponding circular cylinder with a drag reduction up to 33%. The fluctuating lift coefficient is greatly suppressed to be nearly zero. Drag reduction due to the shearing process prevails over that due to the compressing process in this flow. Through investigating the mechanisms relevant to passive control of the tab, it is found that suppression of the shear layer instability can lead to a higher-base-pressure distribution, which can reasonably be associated with drag reduction and suppression of the lift fluctuations. The analysis of convective Mach number and self-sustained oscillations phenomenon inside the shear layers indicates that both the compressible effect and high-frequency forcing result in suppression of the shear layer instability of tabbed cylinder.     

多维流动湍流代数应力模型的比较

介绍了湍流传统代数应力模型和作者新建的代数应力模型的基本思想，用κ－ε模型和传统的及新的代数应力模型进行了多维数值计算，比较了正应力分量的各向异性和宏观湍流速度场。三种湍流模型之间存在明显差异，新的代数应力模型显示具有一定的优越性。     

H2/N2射流扩散火焰的一维全尺度湍流模拟

采用20步的反应机理并应用多分量扩散进行了模拟计算．讨论了Kolmogorov类型能量级串和多重涡效应导致的局部熄火现象．对于中心线发展量与条件统计量的计算符合实验数据．用一维全尺度湍流（ODT）模型对微量组分OH基和O基进行计算，给出了随着火焰发展自由基生成和消耗的信息．     

超声速旋流天然气分离器的旋流特性数值模拟

与传统的低温分离工艺相比，超声速旋流天然气分离器是天然气处理工艺技术的一大创新。在超声速旋流天然气分离器中，气流经过拉伐尔喷管绝热膨胀形成带液滴的超声速低温混合气流，在超声速翼的作用下混合气流由轴流转换成旋流，实现超声速旋流分离。超声速翼是实现气液分离的关键部件。设计了三角薄板型超声速翼，并利用CFD软件对超声速翼段内气流温度、压力、马赫数等特性参数的变化规律和翼段沿主流方向切向速度的变化情况进行了分析。结果表明，在所设计的超声速翼段内，气流能始终保持超声速，翼段出口马赫数为1．4，翼前无激波产生；分离器的旋流加速度最高在572000g，可实现良好的超声速气液旋流分离。     

用应力模型计算旋风分离器的流场

采用一种适合于强旋流的代数应力湍流模型对旋风分离器的流场进行了数值模拟，结果揭示了强旋流中湍流各向异性的特点，显示了强旋流所具有的涡结构及旋风分离器气动分离的机理     

基于RANS和LES湍流模型的二甲醚超临界喷雾数值模拟

分别采用RANS和LES湍流模型耦合KH-RT破碎模型对二甲醚超临界非燃烧喷雾实验进行数值模拟,研究两种湍流模型的适用性,并对比分析二甲醚亚临界喷雾和异辛烷超临界喷雾的发展过程。同时,基于LES湍流模型探讨了环境温度和压力以及喷射温度和压力的变化对二甲醚超临界工况下喷雾贯穿距的影响。结果表明,RANS和LES湍流模型均能较好地模拟二甲醚超临界喷雾,LES湍流模型的模拟结果与喷雾实验观察结果更吻合;与二甲醚亚临界喷雾相比,其超临界喷雾的雾化效果更好;在相同初始条件下,二甲醚超临界喷雾的雾化效果要优于异辛烷超临界喷雾;二甲醚超临界喷雾贯穿距随环境温度、燃油喷射温度和喷射压力的升高而增大,随环境压力的升高而减小。