fluent帮助文件第16章-部分预混燃烧的模拟

只适用于两个进口的情况，这 2进口可以是纯燃料、纯氧化剂， 或者燃料和氧化剂的混合物。two-mixture-fraction模型的情况下，进口数目限制延展到三个， 但是将带来较大的计算量。更多信息，请参考 14.1.1和 15.1.2。 16.2理论 部分混合模型，通过求解一个输运方程来求平均反应进度 （以决定焰峰的位置），和混合物组 份方程 和 。火焰前方（c＝0），燃料和氧化剂是混合的但未燃烧，火焰后边（c＝1）混合 物是燃尽了的。 16.2.1标量的计算 平均标量（如组份质量、温度和密度），用 表示示，是通过计算 f和 c的概率密度函数（PDF）： 在薄火焰的假设下，于是只有未燃反应物和已燃产物存在，平均标量取决于 这里下标 b和 u分别表示已燃和未燃。 已燃部分的标量，φb，是混合物的函数，通过组分燃料质量 f和氧化剂质量（1－f）并使混合 物平衡φb。当非绝热混合物和／或考虑层流 的时候，φb 依然是热焓和／或应力的函数，但 是这并不改变基本公式。未燃部分的标量，φu，类似地通过组合燃料质量 f 和氧化剂质量（1 －f）来计算，但是混合物没有反应。 就像非预混模型一样，化学计算和为已燃混合物 PDF，在 prePDF里面执行，FLUENT已经构建了 向上层工作台供使用。 在未燃混合物里忽视了湍流波动（PDF）和非绝热的影响，所以平均未燃物的标量 只是 的 函数。这些假设对大多数部分预混燃烧了流动是有效的，做这些假设是为了减少内存需求。未 燃的密度、温度、组份质量百分数、比热和热扩散（指示φ）在 prePDF里面用最小二乘法拟合 成 的三阶多项式。 由于未燃物标量是平滑变化、并是 的缓慢变化函数，这些多项式拟合通常很精确。在 prePDF 里面，用户有权改变这些多项式，如果你要改变它们。 16.2.2层流火焰速度 反应进度模型需要层流火焰速度（公式 15.2-4），层流火焰速度取决于混合物组份、温度以及 压力。对很好地预混合系统如 15章里的，反应流有一种混合物，层流火焰速度在整个火焰域里 近似为常数。然而，在部分预混系统里，层流火焰传播速度将随着反应混合物（平很比）的变 化而变化，这种变化必须考虑。 精确的层流火焰速度理论上很难确定，通常是通过试验或者一维摸拟计算出来。prePDF使用拟 合曲线 Goyygens［89］来获得层流火焰速度。这些曲线是为氢气、甲烷、乙炔、乙烯、乙烷和 丙烷燃烧火焰的。进口油气比从贫油极限到化学恰当比、未然物温度从 298k 到 800k，压 力从 1到 40bars的情况下，这些假设都是有效的。 prePDF把曲线拟合成为分段线性多项式。富油极限和贫油极限下的油气比也可以确定，并转化 为混合物百分比。混合物如果比贫油极限更贫油或者比富油极限更富油的话，将不会燃烧，将 出现 0 火焰速度。这要求输入数值 10 倍层流火焰速度。层流火焰速度的最小和最大的 极限 就是输入的第一个和最有一个 数值。 ！！这些火焰速度对于燃料是纯 H2, CH4, C2 H2, C2 H4, C2 H6, 和 C3 H8 是比较精确的。如果 氧化剂不是空气或者有其它的的燃料，那么这种曲线拟合是不正确的。虽然 prePDF缺损是甲烷 空气混合物，层流火焰速度多项式和富油或者贫油状态下最可能不正确。层流速度多项式应该 由其它项决定，比如别的一些相关文献介绍的测量方法或者详细的 1 维模拟，然后输入到 prePDF。 16.3使用部分预混燃烧模型 设置和求解部分预混燃烧问，结合部分非预混和部分预混的设置。16.3.1中给出了这个过程 的大致轮廓，同时还有相应的关于到预混和非预混章节里去查询细节的相关信息。部分预混燃 烧模型的特殊的输入，在 16.3.1和 16.3.2中介绍。 16.3.1设置和求解步骤 1、用 prePDF 生成一个 PDF 的查询表。你可以跟着非预混燃烧章节里，14.3.1 或者 14.4.6 讲 述的步骤进行设置，在 define case面板里，打开 partially premixed model选项。 Setup Case... ！！如果 prePDF 有警告你，在查询计算阶段，任何关于计算层流火焰速度的函数、计算未燃密 度、温度、比热和热扩散率的参数，超出范围的参数都是有效的，你必须在保存 PDF 文件前改 变多项式系数或者分段线性点的系数。 Display Partially Premixed Properties... 细节请参阅 16.3.2。 2、读入网格文件到 FLUENT，并设置好其它你采用的与部分预混燃烧模型相关的模型（湍 流、辐射等） 3、激活部分预混燃烧模型 （a） 打开 Species Model 面板里的 Partially Premixed Combustion 模型。 Define Models Species... （b） 如果必要，改变 Species Model 面板里的 Model Constants（模型常数）。这些对于预 混燃烧模型来说是永恒的常数，很多情况下，你不必改变他们的缺损值。详细情况参见 15.3.2 （c） 点击 Species Model 面板上的 OK 后，指定包含你在 prePDF 里生成的查询图表的 PDF 文件名，并把它读入 FLUENT。 4、定义计算域内未燃物质的物理特性 Define Materials... FLUENT 将自动地为层流火焰速度选择 prepdf-polynomial 函数，指示 PDF 的查询表的分段线 性多项式函数将用来计算层流火焰速度。有可以选择用户自定义函数而不是分段线性多项式函 数。参见 15.3.4，有关于设定未燃物的其它特性的详细资料。 5、设定进出口的火焰进度变量（c）和平均混合物百分数（ ）和其变异（ ）。（对于那些 包含多个进口的，你必须为平均次要进口部分分数定义边界条件以及其变异量） Define Boundary Conditions... 参见 14.3.3节，查看关于指导如何设定混合物分数和变异量条件，以及进口热流量和速度进口。 ！！如果你在 PDF里面定义燃料为进口预混组分，那么你应该设定 set 和 。如果 第定义的是纯燃料，那么你必须设定正确的化学恰当比( )，和进口处的 。 6、初始化进展变量 Solve Initialize Patch... 详见 15.3.6节。 7、求解物理问题并执行后处理 参见 14.3.3 节，详见关于设定求解参数的指导。（这些指导是为非预混燃烧计算的，但是也相 应适合于部分预混燃烧）。 16.3.2在 prePDF里改变未燃物组分的特性多项式 对于未燃物，prePDF计算温度、密度、热容量、热扩散作为平均组分百分数的函数， （公式 16.2-3）和层流火焰速度作为 的分段多项式函数。 未燃物温度，密度，热容，热扩散是通过 CHEMKIN数据库［125］里的热动力和输运特性计算的。 由于它们都是平滑的和 的缓慢变化的函数，所以多项式拟合的处理通常非常精确。 但是，就像在 16.2节里说的，层流火焰速度取决于化学动力学的细节和分子输运特性，不是直 接在 prePDF里计算得来。相反，详细模拟得到的火焰速度，做了曲线拟合。这些拟合只局限于 这些燃料 H2, CH4, C2 H2, C2 H4, C2 H6, 和 C3 H8，空气为氧化剂，贫油极限恰当比为 1，未燃 物温度从 298～ 800k，压力从 1～40bar。如果拟的参数调在这些范围外，prePDF 将会在你计 算查询表时警告你。在这种情况下，必须在部分预混模型特性面板（图 16.3.1）里改变分段线 性点的的取值，然后再保存 PDF文件。 Display Partially Premixed Properties... 图 16.3.1：prePDF里的部分预混模型特性面板 对每一个 的多项式函数（密度，温度，比热和热扩散），你可以指定 c0、c1、c2和 c3的值（公 式 16.2-3里的多项式系数）。 也可以对 10个不同点指定分段线性混合物百分数（f）和相应的层流火焰速度（sl）。第一个指 定的点是下限，最后一点是上限。在上下限范围以外，层流火焰速度是常数并等于极限值。 同时注意，你可以为层流火焰速度选择用一个用户自定义函数，这时分段线性拟合就是无关的 了。