FLUENT算例 (9)模拟燃烧

1、计算流体力学作业FLUENT 模拟燃烧组分传输与气体燃烧 问题描述：长为2m、直径为0.45m的圆筒形燃烧器结构如图1所示，燃烧筒壁上嵌有三块厚为0.0005 m，高0.05 m的薄板，以利于甲烷与空气的混合。燃烧火焰为湍流扩散火焰。在燃烧器中心有一个直径为0.01 m、长为0.01 m、壁厚为0.002 m的小喷嘴，甲烷以60 m/s的速度从小喷嘴注入燃烧器。空气从喷嘴周围以0.5 m/s的速度进入燃烧器。总当量比大约是0.76(甲烷含量超过空气约28%)，甲烷气体在燃烧器中高速流动，并与低速流动的空气混合，基于甲烷喷嘴直径的雷诺数约为5.7×103。假定燃料完全燃烧并转换为:CH

2、4+2O2CO2+2H2O反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。1、 建立物理模型，选择材料属性，定义带化学组分混合与反应的湍流流动边界条件2、 使用非耦合求解器求解燃烧问题3、 对燃烧组分的比热分别为常量和变量的情况进行计算，并比较其结果4、 利用分布云图检查反应流的计算结果5、 预测热力型和快速型的NOX含量6、 使用场函数计算器进行NO含量计算一、 利用GAMBIT建立计算模型第1步 启动GAMBIT，建立基本结构分析：圆筒燃烧器是一个轴对

3、称的结构，可简化为二维流动，故只要建立轴对称面上的二维结构就可以了，几何结构如图2所示。(1) 建立新文件夹在F盘根目录下建立一个名为combustion的文件夹。(2) 启动GAMBIT (3) 创建对称轴 创建两端点。A(0,0,0)，B(2,0,0) 将两端点连成线(4) 创建小喷嘴及空气进口边界 创建C、D、E、F、G点CDEFGx00.010.0100y0.0050.0050.0070.0070.225 连接AC、CD、DE、DF、FG。(5) 创建燃烧筒壁面、隔板和出口 创建H、I、J、K、L、M、N点（y轴为0.225，z轴为0）。HIJKLMNx0.5000.5051.0001

4、.0051.5001.5052.000 将H、I、J、K、L、M、N向Y轴负方向复制，距离为板高度0.05。 连接GH、HO、OP、PI、IJ、JQ、QR、RK、KL、LS、ST、TM、MN、NB。(6) 创建流域 将以上闭合线段创建为面。 第2步 对空气进口边界进行网格划分（1） 划分甲烷进口边界为等距网格 点击Edges右侧黄色区域 按下Shift+鼠标左键，点击AC线段 Type选Successive Ratio,Radio 选1 在Spacing下面白色区域右侧下拉列表中选择Interval count 在Spacing下面白色区域内填入网格的个数5 保留其他默认设置，点击APPLY

5、（2） 划分空气入口边界为不等距网格 选择FG线时，若线段方向由F指向G，则按住Shift键，用鼠标中键点击FG线段，使线段方向由G指向F。 在Type项选择Exponet 在Ratio项输入0.38 Spacing选择Interval size 并输入0.005 点击APPLY（3） 划分小喷嘴壁面为等距网格 把CD、EF线段划分为网格数为4的等距网格 把DE线段划分为网格数为3的等距网格（4） 划分燃烧器出口边界为等距网格把燃烧器出口边界BN划分为35个等距离网格。（5） 划分燃烧器壁面为网格燃烧器壁面由GH、IJ、KL、MN组成 在Edges项选择GH、IJ、KL、MN 在Type项选择

6、Bi-exponent,在Ratio项输入0.55 在Spacing项选择Interval count，并输入62 Apply（6） 对壁筒上的三个隔板进行网格划分 把六个竖直边HO、IP、JQ、KQ、LS、MT分别划分为10个等距网格 把三个横边OPQRST分别化为2个等距网格（7） 对整个计算域进行面网格划分 点击Face右侧黄色区域 按下Shift+鼠标左键，点击面上的边线 在Elements选择Quad 在Type项选择Pave 在Spacing项选择Interval size,并输入网格间距0.008 Apply第3步 设置边界类型并输出文件（1） 设置甲烷速度入口边界 在Actio

7、n项为Add 在Name 项填入边界名inlet-fuel 在Type项选择WELOCITY_INLET 点击Edges右侧黄色区域 按住Shift键点击AC线段 Apply（2） 设置空气速度入口边界 在Name 项填入边界名inlet-air 在Type项选择WELOCITY_INLET 在Edges项选择FG线段 Apply（3） 设置压力出流边界 在Name 项填入边界名outlet 在Type项选择PRESSURE_OUT 在Edges项选择BN线段 Apply（4） 设置对称轴边界 在Name 项填入边界名axis 在Type项选择axis 在Edges项选择AB线段 Apply（

8、5） 设置小喷嘴的边界类型 在Name项填入边界名zozzle 在Type下选择WALL 在Edges项选择CD、DE、EF 点击apply（6） 输出网格文件 在File Name项确认文件名 选择Export 2-D(X-Y)Mesh Apply二、 利用FLENT-2d求解器进行模拟计算第1步 启动FLENT-2d求解器，读入网格文件。（1） 启动FLUENT-2d求解器（2） 读入网格文件combustion.msh（3） 检查网格（4） 网格信息 （5） 网格长度单位设置 （6） 显示网格 第2步 设置求解模型（1） 设置求解器 在Solver项选择Segregated 在Formu

9、lation项选择Implicit 在Space项选择Axisymmeric 在Time项选择Steady OK （2）选用k-湍流模型 在Model项选择k-epsilon OK （3）激活能量方程 选择Energy Equation OK （4）启动化学组分传输和反应 在Model先选择Species Transport在Reaction下选择Volumetric在Options下选择Diffusion Energy Source 在Mixture Material下拉列表中选择methane-air 在Turbulence-Chemistry Interaction 下选择Eddy-Di

10、ssipation OK 第3步 流体材料设置 在Denity下拉列表中选择incomprehensible-ideal-gas 在Cp项选择Constance,输入1000 点击Mixture Species 右边的Edit。 点击Cancel 在Material面板中，点击Reaction下拉列表右边的Edit OK 使用滚动条检查其余的物性 点击Chang/Create，接受材料物性的设置并关闭对话框 第4步 设置边界条件（1） 打开边界条件面板（2） 设定空气进口inlet-air的边界条件 在Zone项选择inlet-air 确定在Type项为velocity-inlet 在Velo

11、city Magnitude项输入空气入口速度0.5 在Turbulence Specification Method 项选Intensity and Hydraulic Diameter 在Turbulence Intensity项输入10 在Hydraulic Diameter项输入燃烧筒直径0.45 Species Mass Fractions 项均为常数，且在O2项输入0.22 点击OK（3） 设定燃料进口边界条件 在Zone 项选择inlet_fuel 确定Type项为velocity-inlet ,点击Set,打开燃料速度入口边界设置对话框 设置后点击OK（4） 设定压力出口边界条

12、件 在Zone项选择outlet 确定Type项为pressure-outlet点击set，打开压力出流白边界设置对话框 进行设置，点击OK（5） 设定燃烧筒外壁的边界条件 在Zone项选择wall 点击set，打开壁面边界条件设置对话框 在Thermal选项卡中的Thermal Conditions项选择Temperature 在Temperature项输入温度300 保留其他默认设置，点击OK（6） 设置燃料进口喷嘴壁面的边界条件 在Zone项选择nozzle 点击set，打开喷嘴壁面边界设置对话框 在Thermal选项卡中Thermal Conditions项 选择Heat Flux 在

13、Heat Flux项保留默认的零值 保留其他默认设置，点击OK 第5 步 初始化流场并求解（1）设置求解控制参数 打开求解控制参数设置对话框，在Under-Relaxation Factors项，设置每个组分的松弛因子为0.8 保留其他默认设置，点击OK（2）流场初始化 在Compute From下拉列表中选择all-zones 设置CH4为0.2 调整温度初始值到2000 点击Init（3）在计算期间打开残差图形监视器 打开残差监视器设置对话框在Options下，选择Print和Plot调整energy残差收敛标准为1e-05保留其他默认设置，点击OK（4）保存case文件 打开文件保存对话

14、框，键入文件名combustion1，点击OK（5）进行1000步迭代计算 打开迭代计算对话框。填入1000 （6）保存case和data文件 Case&Data 保存的文件名为combustion1.cas 和 combustion1.dat（7） 绘制温度分布云图 打开绘制分布云图设置对话框 在Option项选择Filled 在Contours of 下拉列表中选择Temperature.和Static Temperature 保留其他默认设置，点击Display Fluent是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%，凡是和流体、热传递和化学反应等有关的工

15、业均可使用。它具有丰富的物理模型、先进的数值方法和强大的前后处理功能，在航空航天、汽车设计、石油天然气和涡轮机设计等方面都有着广泛的应用。CFD商业软件FLUENT，是通用CFD软件包，用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术，因而FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型，使FLUENT在转换与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。目前与FLUENT配合最好的标准网格软件是ICEM，而不是早已过时的GAMBIT。

16、FLUENT软件包含三种算法：非耦合隐式算法、耦合显式算法、耦合隐式算法，是商用软件中最多的；FLUENT软件包含丰富而先进的物理模型，使得用户能够精确地模拟无粘流、层流、湍流。湍流模型包含Spalart-Allmaras模型、k-模型组、k-模型组、雷诺应力模型(RSM)组、大涡模拟模型(LES)组以及最新的分离涡模拟(DES)和V2F模型等。另外用户还可以定制或添加自己的湍流模型；适用于牛顿流体、非牛顿流体；含有强制/自然/混合对流的热传导，固体/流体的热传导、辐射；化学组份的混合/反应；自由表面流模型，欧拉多相流模型，混合多相流模型，颗粒相模型，空穴两相流模型，湿蒸汽模型；融化溶化/凝固

17、；蒸发/冷凝相变模型；离散相的拉格朗日跟踪计算；非均质渗透性、惯性阻抗、固体热传导，多孔介质模型（考虑多孔介质压力突变）；风扇，散热器，以热交换器为对象的集中参数模型；惯性或非惯性坐标系，复数基准坐标系及滑移网格；动静翼相互作用模型化后的接续界面；基于精细流场解算的预测流体噪声的声学模型；质量、动量、热、化学组份的体积源项；丰富的物性参数的数据库；磁流体模块主要模拟电磁场和导电流体之间的相互作用问题；连续纤维模块主要模拟纤维和气体流动之间的动量、质量以及热的交换问题；高效率的并行计算功能，提供多种自动/手动分区算法；内置MPI并行机制大幅度提高并行效率。另外，FLUENT特有动态负载平衡功能，确保全局高效并行计算；FLUENT软件提供了友好的用户界面，并为用户提供了二次开发接口（UDF）；FLUENT软件采用C/C+语言编写，从而大大提高了对计算机内存的利用率。在CFD软件中，Fluent软