三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。     

带隔热屏加力燃烧室热态流场计算

在任意曲线坐标系下对带有横向波纹隔热屏、外冷却通道和尾喷口的涡喷加力燃烧室的三维热态流场进行数值模拟.利用椭圆型微分方程和区域法生成三维贴体网格, k-ε模型预估紊流特性, EBU-Arre-henius紊流燃烧模型估算化学反应速率, 六通量模型预估辐射通量以考虑其对壁面和气流温度的影响.在非交错网格系统下应用混合差分格式离散控制方程, SIMPLE算法求解离散方程, 数值分析两种进口气流参数分布对加力室热态流场和壁温分布的影响.所得计算值与试验数据基本相符, 表明计算方法合理.      

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

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

RNG k-ε模型数值模拟油雾燃烧流场

采用RNG(R enorm alization G roup)k-ε紊流模型对某环形燃烧室火焰筒内三维两相燃烧流场进行数值模拟。运用偏微分方程法和区域法生成贴体网格,EBU-A rrhen ius紊流燃烧模型估算化学反应速率,六通量热辐射模型估算辐射通量。在非交错网格下求解气相采用S IM PLE算法,液相采用颗粒轨道模型和PS IC算法。近壁区处理分别采用两层和三层壁面函数。计算数据与试验值比较表明,采用三层壁面函数的RNG k-ε模型更适用于模拟三维两相燃烧流场。      

模型加力室大涡模拟的并行计算

采用基于MPI(消息传递库)的并行算法,在贴体网格下对带V形槽稳定器模型加力燃烧室紊流化学反应流场进行数值模拟,湍流模型采用k方程亚网格尺度模型,燃烧模型采用亚网格EBU模型,采用热通量辐射模型估算辐射通量。在程序设计中,采用动态内存分配、分区算法和多点重合交错网格系统,并行计算的结果与单机运行结果的对比表明计算结果是正确的,可以明显的提高运算效率,是解决复杂燃烧流动大规模数值模拟的有效手段。      

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

回流燃烧室燃烧过程的三维数值模拟

在三维任意曲线坐标系下数值模拟回流燃烧室火焰筒内两相燃烧过程,采用RNG k-ε模型模拟紊流粘性,EBU-Arrhenius模型模拟燃烧反应速率、离散坐标模型以及六通量模型考虑辐射传热,液相采用颗粒轨道模型,气相采用SIMPLE算法求解,并用PSIC算法考虑气液两相之间的相互作用的影响,计算得到燃烧室内速度、温度等各气流参数分布.通过将计算与实验结果对比表明,计算方法可靠,离散坐标模型优于六通量辐射模型,更适用于模拟火焰筒内两相燃烧流场.      

二维气膜冷却火焰筒壁温数值分析

采用数值试验方法，研究了二维气膜冷却燃烧室火焰筒的壁温分布，预估了来流速度、温度和吹风比对壁温及辐射对换热规律的影响。分析中、采用二元热流通量模型模拟辐射换热，用“K-ε”两方程模型描述紊流特性。用壁面函数法修正处于流场内部的固壁。     

加力燃烧室单相燃烧模型及其应用

提出了描述航空发动机加力燃烧室内气相湍流燃烧的模型方程，并考虑了加力喷油效应。通过结合双时间尺度ｋ－ε湍流模型、热通量高温辐射模型及ＳＩＭＰＬＥ方法，对实际运行的加力燃烧室内的湍流燃烧进行了数值模拟，得到了数值流场图及温度、组分浓度的等值线分布，给出了燃烧效率和总压恢复系数，数值解证明了湍流燃烧模型的可解性与合理性。     

数值分析二级涡流器环形燃烧室的燃烧性能

采用多维经验分析法在任意曲线坐标系下对包括二级突扩扩压器、双级轴向涡流器及火焰筒在内的单头部环形燃烧室的燃烧性能进行计算.并对RNGk-ε紊流模型加以改进, 使其更适用于数值模拟非交错贴体坐标系下的三维两相化学反应流.采用EBU-Arrhenius紊流燃烧模型、六通量热辐射模型以及颗粒群轨道模型对油雾燃烧进行模拟.同时, 数值分析了不同涡流器几何尺寸对燃烧室整体流场和燃烧性能的影响, 数值计算结果与实验数据相当一致, 表明计算方法合理, 计算程序可靠, 可用来估算环形燃烧室的燃烧性能.      

双旋流器燃烧室NOx生成研究

通过试验和数值模拟方法研究不同进口油气比与温度以及主燃孔布局变化对双旋流模型燃烧室两相燃烧整体流场、出口温度丁以及污染物NO_x排放的影响.采用二维粒子图像测速仪(PIV)测量燃烧室内热态流场,用多点温度耙和取样管测量燃烧室出口温度T与NOx在非结构网格体系下,采用Fluent商业软件对燃烧室三维两相燃烧流场和NO_x分布进行计算.计算结果与试验数据符合较好,表明不同油气比、进气温度和主燃孔布局变化对模型燃烧室两相燃烧整体流场、出口温度T与NO_x排放都有一定的影响,并随着油气比与进气温度增加,燃烧室出口温度T与NO_x排放也相应地提高.      

燃烧室流动的数值计算

提出了多连通域流场统一计算的处理方法,采用三维贴体曲线网格,非交错网格法和k-ε湍流模型。对包括扩压器,内外环通道和火焰筒的整个燃烧室流场进行了数值计算。可以计算流过主燃孔、掺混孔和气膜孔的流量。计算结果定性合理。      

大涡模拟模型燃烧室燃烧性能计算

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 k 方程亚网格尺度模型模拟亚网格湍流黏性;亚网格EBU(eddy-break-up)燃烧模型预估化学反应速率;多维经验分析法计算燃烧性能;并在非交错网格体系下气相采用SIMPLE(semi-implicit method for pressure-linked equations)算法对控制方程进行求解,液相采用随机离散模型,两相之间的耦合采用PSIC(particle-source-in-cell)算法.通过大涡模拟瞬态及时均计算结果表明:与粒子图像测速仪(PIV)测量的瞬态速度场、出口温度分布试验数据吻合,表明在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法,数值模拟模型燃烧室两相喷雾燃烧流场,所采用的亚网格模型可以用于燃烧室气液两相喷雾燃烧流场的大涡模拟;燃烧性能计算结果与试验测量结果基本一致,说明所采用多维经验分析法可以用来数值模拟航空发动机燃烧室燃烧性能的计算,特别是污染物的预估,为设计低污染高性能航空发动机燃烧室提供有用的设计依据.      

带圆筒头部火焰筒两相燃烧流场及性能计算

在三维任意曲线坐标系下采用代数雷诺应力紊流模型模拟紊流粘性,Arrhen inus-EBU紊流模型模拟燃烧室内化学反应速率,采用随机轨道模型模拟气液两相之间的相互作用,数值分析不同进口气流参数对带双圆筒头部环形燃烧室两相紊流燃烧流场的影响,同时采用多维经验分析法预估燃烧室性能。通过两种工况计算表明:在最大工况下所得的燃烧流场和燃烧室性能均优于慢车工况。计算结果与试验数据符合较好,表明燃烧过程的数值模拟以及多维经验分析法可为燃烧室优化设计提供有用的设计依据,适合于工程应用。     

求解流体力学反问题的数学方法及应用

流体力学反问题是一类重要的非线性问题，这类问题在工程中有广泛的应用，如过坝水流问题，给定流量的射流和翼型设计问题。由于这类流动求解区域不确定，导致问题的非线性特征，使数值计算十分困难，本文将以亚、跨声速翼型设计问题为例，提出了一种求解流体力学反问题的新方法。