模型燃烧室湍流亚网格尺度模型

分别采用Smagorinsky-Lilly和k方程亚网格尺度模型,在任意曲线坐标系下对模型燃烧室湍流流场进行大涡模拟。分别用QUICK和二阶隐式差分格式离散控制方程及其时间项,并用SIMPLE算法进行求解。计算结果与PIV测量数据比较表明:大涡模拟方法能较好模拟旋流器后面回流区的瞬态变化,两种亚网格尺度模型都能较好模拟湍流流动,但k方程亚网格尺度模型稍优于Smagorinsky-Lilly亚网格尺度模型,因此采用k方程亚网格模型更适用于大涡模拟实际燃烧室复杂湍流流场。

亚网格尺度稳定化有限元求解不可压黏性流动

从亚网格尺度稳定化方法的基本原理出发,提出了适合时间推进求解非定常Navier-Stokes方程获得定常解的SGS稳定化方法.基于一定程度的近似和简化,获得了与时间步长相关的稳定化参数,从而排除了传统SGS稳定化方法在求解高Re数、小时间步长问题时所引发的数值不稳定性.把SGS稳定化方法应用于求解不可压湍流,结合标准k-ε湍流模型和壁面函数法估计湍流黏性系数,详细讨论了壁面函数法的实施、湍流输运方程的求解和保证湍流变量非负性的限制策略,发展了时间推进求解不可压湍流的分离式算法.二维外掠后台阶层流和湍流计算结果表明,该方法求解不可压黏性流动是可行的,并且具有稳定性好、计算精度高的特点.

用二阶矩亚网格尺度燃烧模型对丙烷-空气旋流扩散燃烧的大涡模拟

用二阶矩(SOM)亚网格尺度燃烧模型对环缝进燃料的丙烷-空气旋流湍流扩散燃烧进行了大涡模拟(LES).模拟得到统计平均的热态3个方向的速度、湍流度、温度、丙烷、氧和CO2浓度分布,其值与实验数据符合很好.结果表明,二阶矩(SOM)亚网格尺度燃烧模型适用于大涡模拟.环缝进气使湍流脉动强度、各向异性程度和温度分布趋于均匀.

气粒流动两相亚网格尺度动能模型的大涡模拟

目前多数两相流动的大涡模拟采用单相流动的亚网格尺度应力模型,没有考虑或者没有完整地考虑两相亚网格尺度应力间的相互作用。本文提出一种两相亚网格尺度动能模型(LES-kkp),可以完整地考虑其相互作用。将该模型用于突扩气粒两相流动的大涡模拟,同时进行统一二阶矩(USM)两相湍流模型的雷诺平均模拟(RANS)。瞬态结果显示出两相流动中的气体形成了典型的拟序结构,但颗粒流动结构完全不同于气体的,没有形成涡结构。LES-kkp和RANS-USM模拟结果实验检验表明,对时间平均速度两种模拟给出结果与实验结果吻合良好,即RANS-USM得到了LES-kkp的验证。对两相均方根脉动速度和两相脉动速度关联量,LES-k-kp的模拟结果优于RANS-USM的模拟结果。

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性；热通量辐射模型估算辐射换热,分别采用亚网格EBU燃烧模型(E-A model)、亚网格二阶矩输运方程模型(SOM)和亚网格二阶矩代数模型(SOM-A)估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解,液相采用Lagrange处理,并用PSIC算法对其进行求解．通过实验结果和计算结果的比较,表明在三维任意曲线坐标系下对燃烧室火焰筒两相湍流油雾燃烧流场进行大涡模拟,3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程,而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．

环形燃烧室两相喷雾燃烧的大涡模拟

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对环形燃烧室火焰筒气液两相紊流瞬态反应流进行数值模拟。采用椭圆偏微分方程生成三维贴体网格,计算中所采用的数学度模型有:k方程亚网格尺度模型估算亚网格紊流粘性;亚网格EBU燃烧模型估算化学反应速率;热通量辐射模型估算辐射换热。并在非交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相采用随机离散模型(Stachastic Separated Flow,简称SSF),在拉格朗日坐标系下追踪各油珠群沿各自轨道运动、质量损失及能量变化。通过计算结果与实验数据相比较,表明在三维贴体坐标系下对燃烧室火焰筒两相紊流油雾燃烧流场进行大涡模拟,采用欧拉-拉格朗日两相大涡模拟方法能反映两相紊流化学反应流流动及实际燃烧过程。

加力燃烧室热态流场的大涡模拟

本文采用k方程亚网格尺度模型对带V形槽稳定器模型加力燃烧室紊流化学反应流进行大涡模拟的研究,采用亚网格EBU燃烧模型估算化学反应速率,用热通量辐射模型估算辐射通量。数值计算表明,在稳定器后面开始在短时间内出现涡的交替脱落和逐渐消失过程,然后形成稳定的回流区,数值计算结果与实验比较吻合,表明可用大涡模拟研究实际燃烧室瞬态燃烧流场。

算子分裂法求解旋转坐标系下不可压黏性流动

提出了求解旋转坐标系下的不可压黏性流动问题的θ格式算子分裂算法.通过算子分裂,把不可压缩性、非线性和哥氏力占优三大耦合困难分割开来.采用亚网格尺度稳定化方法消除了Galerkin方法求解时由于不可压缩性和哥氏力占优所引发的数值振荡.结合最小二乘和共轭梯度法间接求解非线性子问题,排除了强对流作用所引发的数值振荡,避免了引入迎风格式离散对流项的必要性.同时该算法保证了迭代过程中有限元总刚度矩阵正定不变的特性,为求解线性方程组采用高效的求解器提供了可能.数值试验表明,该算法具有稳定性好、收敛速度快、计算精度高的特点.

加力燃烧室污染特性的大涡模拟

采用k方程亚网格尺度模型对模型加力燃烧室的污染特性进行大涡模拟研究,其中污染物CO亚网格尺度模型是以二步化学反应为基础,NO采用亚网格动力学模型,并用亚网格EBU燃烧模型估算化学反应速率,用热通量辐射模型估算辐射通量,通过数值计算结果与实验数据的比较,表明采用的亚网格污染物生成模型可以用于实际加力燃烧室污染物排放的预估.

模型燃烧室紊流燃烧的大涡模拟

采用了两种不同的亚网格尺度燃烧模型对带 V型稳定器的模型燃烧室紊流化学反应流动进行了大涡模拟 ,用 k-ε方程亚网格尺度模型确定亚网格紊流粘性 ,为了考虑热辐射对燃烧室壁温和气流温度的影响 ,运用热流法辐射模型估算热辐射通量 ,用 SIMPLE算法和混合差分求解大涡模拟各守恒方程 ,通过对两种不同亚网格尺度燃烧模型数值模拟结果与实验值的比较表明 ,两种燃烧模型都与实验值较吻合 ,但 G方程小火焰模型要比亚网格

贴体坐标系下模型加力室的大涡模拟

本文利用贴体网格对带V形槽稳定器模型加力燃烧室素流化学反应流流动进行大涡模拟的研究。采用区域法生成模型加力燃烧室的二维贴体网格,并采用多区域耦合法进行区域之间的数据传递,求解加力室整体流场。采用k方程亚网格尺度模型和亚网格EBU燃烧模型分别估算其亚网格紊流粘性和化学反应速率,用热通量辐射模型估算辐射通量,并用交错网格下SIMPLE算法和混合差分格式求解离散方程,壁面函数处理固壁边界条件。计算结果显示了稳定器后面的回流区气流结构,所得的热态流场模拟结果与实验比较吻合,表明采用贴体网格对模型加力燃烧室进行大涡模拟能真实反映流体流动及燃烧过程。

大涡模拟在加力燃烧室冷态流场中的应用

对带 V型稳定器模型加力燃烧室二维冷态流动进行了大涡模拟。采用了代数和 k方程两种亚网格尺度模型 ,并用 SIMPLE算法和混合差分格式求解离散方程 ,用壁面函数处理固壁边界条件 ,得到了模型加力燃烧室二维瞬态流场的模拟结果。通过对不同亚网格尺度模型和不同入口速度下的瞬态流场的比较 ,揭示了不同亚网格尺度模型和不同入口速度分布对带 V型稳定器燃烧室二维冷态瞬态流场的影响。

燃烧室污染物生成大涡模拟模型的研究

采用亚网格动力学模型和联合概率密度函数(PDF)亚网格模型对污染物NO的生成进行大涡模拟.计算中采用的数学模型有:k方程亚网格尺度模型来估算亚网格湍流粘性;亚网格EBU燃烧模型估算化学反应速率;在交错网格体系下,采用SIMPLE算法和混合差分格式来求解差分方程.数值分析两种不同亚网格污染物生成模型和两种进口气流温度对NO排放的影响,计算结果与实验数据比较表明:两种亚网格污染物生成模型都能较好的预估燃烧室污染物NO的生成,但亚网格动力学模型的计算工作量要比联合概率密度函数亚网格模型少得多,更适合工程应用.

模型燃烧室内瞬态紊流的大涡模拟

本文采用三种不同亚网格尺度模型对带有V型稳定器的模型燃烧室二维瞬态紊流流动进行了大涡模拟。并在交错网格系下用SIMPLE算法和混合差分格式求解离散方程。数值研究拟不同型式入口速度分布和不同亚网格尺度模型下模型燃烧室二维瞬态紊流流场。计算结果表明不同入口速度分布和不同亚网格尺度模型对瞬态流场和出口速度分布有一定的影响。本文通过数值模拟,揭示了V型稳定器后旋涡的产生和脱落过程。通过计算结果及实验数据的比较可知,本文采用的亚网格尺度模型可以用来模拟模型燃烧室紊流流场及稳定器后面回流区的流动情况。

多区域湍流燃烧大涡模拟

为了提高湍流大涡模拟方法在计算复杂工业燃烧设备中反应流体的效率和精度,提出了可以对结构复杂燃烧设备喷嘴内部无反应区域只求解基本流体控制方程,对喷嘴下游燃料反应区域加入燃烧模型进行分区域计算的多区域求解方法,并采用该方法和传统单区域求解方法分别对甲烷-空气预混合湍流燃烧进行大涡模拟.研究结果表明：多区域和单区域求解方法计算所得的轴向和径向时均速度、脉动速度以及OH质量分数一致;不过,多区域求解方法在计算速度上有明显优势,当CPU数量不超过64个时,计算速度可提升14.7%~20.6%;多区域和单区域方法的数值模拟结果与实验值进行了对比,三者吻合较好.

基于喷流拟序结构预测的SGS模型比较研究

拟序结构是湍流机理研究的重要对象,而亚网格尺度(sub-grid scale,SGS)模型的准确性是高精度大涡模拟的关键,对正确预测湍流拟序结构有重要影响.本文对马赫数0.9的可压缩喷流开展大涡模拟,使用4阶空间截差和3阶时间截差的高精度有限差分算法,亚网格尺度模型分别采用了Smagorinsky模型(Smagorinsky model,SM)、拟序结构动能模型(CKM)、选择多尺度模型(selective mixed-scale model,SMSM)、局部动态Smagorinsky模型(localized dynamic Smagorinsky model,LDSM)和拟序结构模型(coherent-structure Smagorinsky model,CSM)以及对平均耗散、瞬时涡结构、脉动速度POD主导模态的分析表明,u′的主导模态呈带状,在下游分岔或破碎,呈现出多尺度特征,v′的主导模态则呈流向排列的肋状,矢量(u′,v′)两端呈环流模式,表征流动卷吸,w′的主导模态呈流向排布的脊状,表征流场受展向拉伸的模式.脉动速度的POD模态对亚网格尺度耗散敏感,其中CKM模型预测的峰值耗散区对应u′模态的低谷区,因而预测的环流模式不明显;SM模型则未能预测出w′的脊状模态.SMSM,CSM和LDSM模型均较好地反映了湍流中涡结构的多尺度特性,揭示了流动卷吸、脊状的展向拉伸等流动模式,其中CSM模型的计算效率较高.