(2)由于方程的非线性引起迭代值波动而导致发散。采取的措施为：对每个因变量选取 适当的TDMA横扫次数，使求得的变量值与相应的系数及源项相适应，迭代过程采松弛法， 以减小变量的变化幅度。旋流数大时，迭代开始各变量的波动很大，相应地应加强低松弛， 随着迭代过程的进行，低松弛可逐步调整。 (3)采用伪源项使计算过程稳定。伪源项的形式为： St=1mmt1(091d-op) 式中，m。。+为通过每个控制容积的净质量流量。 2.2收敛准则 节点差分方程的余数为： R,=∑Gabab+b-appp 工作中注意到各文献中采用SIMPLE方法时，所用的收敛准则并不一致，文献〔8)使 用引R,1/(Q×)1a,Q为质量流量，文献〔9〕使用lRs/N,N为节点总数，文献〔10) 中只使用压力校正方程的质量源项b为收敛准则，还有的文献没有提到所用的收敛准则。这 种现象不利于对湍流模型的比较和讨论，特别是对复杂的流场。本工作中将上述几种准则做 了比较。表1示出计算收敛时，用不同准则的打印结果。可见在相同条件下，若采用不同的 收敛谁则和标准，所得的结果不一致。本计算收敛标准为|R,mx<0.005。 表1不同准则当收敛时的打印结果 Tablel Results based on different convergence criteria Va Vr Va P! IR0|s/(Q×)1a 0,949 0,375 0.651 0,0896 0,145 」R#|max 0.00499 0.00152 0.00396 0,00017 0.00467 1R/N <0,00499 <0.00152 <0,00396 <0,00017 0.00467 IbImax 0,00017 3 计算及实验结果的分析与比较 将旋流叶片安装角为60°的工况，计算与实验结果进行比较，因篇幅关系，有的图未列 入，请查文献〔2)。图2所示计算与实验的轴向速度分布，两者定性基本相符，但计算值普遍 比实验值大，且计算结果中心线上轴向速度开始衰减较慢。分析其原因，可能是实验测量的 误差较大，对C.的调试也不充分，C值在全流场中应不一样。例如，燃烧室前部的C,值应 比后部的大，此外，对湍流模型也还待做其它改进。 另外就流场切向速度分布，计算与实验结果表明，切向速度分布表现为一个自由涡包围 着一个固体祸。值得注意的是，在未调试C。前切向速度在径向上的分布，在燃烧室后部仅 表现出固体涡的特征。 图3所示计算的子午面上无因次流函数流谱，从它及计算和实验的流场子午面上的速度 455由于方程的非线 性引起迭代值波 动而 导致发散 。 采 取 的措 施 为 对每个因 变量选取 适 当的 横扫次数 ， 使求得 的变量 值与相 应的 系数 及 源项相 适 应， 迭代过程采松 弛法 ， 以减小 变量的 变化幅度 。 旋流数大时 ， 迭 代开始各变量的波动很大 ， 相 应地 应加强低松弛 ， 随着迭代过程的进行 ， 低松 弛可逐 步调整 。 采用伪源项使计算过程 稳定 。 伪 源项 的形式为 ， “ 二 、 功 ’ 一 毋 式中 ， 二 、 为通过每个控制容积的净质量流量 。 收敛准 节点差分方程的余 数为 ， 。 、 功 。 、 一 诱 工作 中注意到各文 献中采 用 方法时 ， 所用的 收 敛准则并不 一致 ， 文献 〔 〕 使 用 】 ， 】 功 二 ， 为质量流 量 ， 文 献 〔 〕 使 用 ， ， 为 节 点总数 ， 文献 〔 〕 中只使 用压力 校正方程 的质量 源项 为收敛准则 ， 还有的文 献 没 有提 到所 用的收敛准 则 。 这 种现象不 利于 对湍流 模型的 比较和 讨论 ， 特别是对复杂的流 场 。 本工作中将上 述几种准则做 了比较 。 表 示 出计 算收 敛时 ， 用不 同准则的打 印结果 。 可 见在相 同条件下 ， 若采 用不同的 收敛准则和标准 ， 所得的 结果不 一致 。 本计算收敛标准 为 ， 二 二 。 豪 不 同准则 当收数 时的 打 印 结果 ‘ 户 ， 价 刃 户 价 币 万 ‘ 。 。 。 。 。 。 。 。 。 。 。 。 。 计算及实验结果 的分析与比较 将旋流叶片安装角为 。 的 工况 ， 计算与实验结果进行 比较 ， 因 篇 幅关 系 ， 有的 图未列 人 ， 请查文献 〔 〕 。 图 所示计算与实验的轴 向速度分布 ， 两 者定性基本相 符 ， 但计 算值普遍 比实验值大 ， 且计 算结 果 中心线上轴 向速度开始衰减较慢 。 分析 其原因 ， 可 能是实验测量的 误差较大 ， 对 ， 的调 试也不充分 ， ， 值在全流 场 中应不 一样 。 例如 ， 燃烧室前部的 ， 值 应 比后部的大 ， 此 外 ， 对湍流 模型也 还待做其它 改进 。 另外就流场切 向速度分 布 ， 计 算与实验结 果 表明 ， 切 向速度分布表现 为 一个 自由涡包 围 着一个 固体涡 。 值得 注意 的是 ， 在未调 试 ， 前 切 向速度在径 向上 的分布 ， 在燃 烧室后部仅 表现 出固体涡的特 征 。 图 所示计算的子午面上 无 因次流函数流谱 ， 从它及计算和 实验的流场 子午面上的速度