基于低矮建筑物实测数据的改进湍流物理模型

要准确模拟低矮建筑物壁面上分离涡位置的平均风压比较困难,这些区域在绕钝体流动中也是很关键的部位.本文在Spalart-Allamaras湍流物理模型的基础上,根据钝体绕流流场特性,结合低矮建筑物上的实测数据,调试并选取合适的参数,提高涡粘预测的准确性,提出了改进的S-A湍流物理模型.通过对6m立方体和Silsoe低矮建筑模型进行稳态数值模拟,由改进S-A湍流模型得到的屋面、背风和侧风面上平均风压结果,比目前公认的模拟钝体绕流效果较好的SST湍流模型有较大改善,比风洞实验更接近实测结果.由此可以看出,本文改进湍流模型的方法,将是一条提高预测建筑物表面平均风压精度很有效的途径.     

基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.      

运用SIMPLE方法计算大曲率弯道中的湍流流动

在三维任意贴体非正交曲线坐标系下应用同位网格技术,发展了适于求解复杂流动的压力修正算法,采用速度的协变物理分量为动量方程中速度的求解变量,湍流模型选用标准的κ-ε模型.计算表明,方法具有良好的收敛性、准确性和稳定性,计算值与试验数值吻合良好,而且能预示通道内复杂的流动现象.     

火烧炉内气相燃烧过程的三维数值模拟

本文对某工业火烧炉内的三维湍流气相燃烧场特性进行了数值模拟,湍流模型采用可实现k-ε双方程模型,湍流燃烧采用有限速率/涡耗散模型,辐射换热采用P1模型,数值方法采用SIMPLE算法,实现了炉内火烧试验的数值再现,特别考察了入口压力、燃料/空气比、多烧嘴不同组合工况等因素对炉内温度场和流场的影响,并对该火烧炉的运行提出了一些合理化建议,为其实际运行提供了有益的参考。     

基于小波分析的缩减湍流数据采集量的新方法

由于湍流实验的数据采集量的大小通常没有准则,实验中为了降低偶然事件的影响通常都使用较大的数据量,从而在数据的采集、存储和后处理过程中都造成资源和时间的浪费.笔者根据小波分析所具有的良好的时频双局域性的特点,在文献[4,5]的基础上,配合统计学方法,在理论上提出了以小波方法缩减湍流实验数据采集量和某些流场计算方法计算量的方法.以使用小波方法分析湍流边界层的湍动能方法为例,证明了这种方法的合理性和可行性.     

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

DENTON程序湍流模型改进研究

改进了DENTON程序所使用的混合长度模型,使之能够适应壁面发展以及有适度分离流动的情况,并减少对经验的依赖.通过对NASARotor67转子和Rotor37转子的对比计算分析,验证了改进的混合长度模型能够进一步提高计算的可靠性,并且计算时间增加很少,保持了DENTON程序计算速度快的优势.      

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

磁流体湍流中的结构和能量传输

磁流体湍流中存在多种相干结构,包括电流和涡旋结构。文章主要对近几年有关磁流体湍流相干结构及其相关的能量传输的一些工作进行了介绍。不同于中性流体湍流中的管状结构,磁流体中的相干结构多呈现为片状,强电流片附近也存在强涡片结构。磁流体湍流中的能量级串使能量跨尺度从大尺度传递至小尺度;与中性流体湍流能量传输主要集中在相近尺度上不同,磁流体湍流除了动能传输,还有磁能传输,以及动能与磁能之间的跨尺度互相转换。研究表明,磁流体湍流中的能量传输及耗散具有强间歇性,集中在较小区域范围,且与相干结构存在关联。在存在背景磁场的情况下,磁流体湍流中的相干结构在背景磁场方向上被拉长,而同样方向上的能量传输受到抑制。在高马赫数的情况下,激波的产生会使能量传输明显增强,同时带来很强的间歇性,结构函数标度律远远偏离线性标度关系,并且随着阶数的增加,标度指数会达到饱和值。     

厚覆冰导线升力突变机理及数值模拟研究

针对新月形厚覆冰导线的升力系数在风攻角15°附近存在突变的问题,分别采用基于k-ωSST湍流模型的雷诺时均法和大涡模拟(LES)的数值方法对新月形厚覆冰导线在风攻角10°～20°范围进行了模拟。通过对比两种数值方法计算得到的覆冰导线气动力系数、流场结构和表面风压,发现LES方法能够更好地捕捉新月形覆冰导线表面的小尺度涡结构,得到的覆冰导线气动力参数计算结果与风洞试验数据高度吻合;而k-ωSST湍流模型难以模拟壁面上小尺度涡,捕捉不到升力系数的突变。根据覆冰导线不同壁面区域的压力分布,发现上侧壁面处的涡结构影响整体流场,并在下侧壁面曲率、来流夹角和壁面切线方向共同作用下导致升力系数突变。LES的气动力参数模拟结果可为覆冰导线防舞提供参考。