参数和系数的数值见表1 k,e. P.Y e 图3交错网格示意图 图2计算用座标系示意图 表1 计算中所用的参数和系数 名称 CD C1 C: ke×104 ecx104 (m2/S2)(m2/s3) dopt △T 数值 0.09 1.0 1.31.42 2.23 9.95/19.51.09/2.97 0.997 0.06 名称N JN KN d(m) X(m) (m) Z(m) Ve(m/s) L(m) 数值·10 10 201.96X10-3 95X10-3952X10-3665X10-32.23/3.1295X1013 在使用这些系数下，尤其使用了ao:最佳人工阻尼因子（在程序中自动实现），本 程序的收敛速率是很快的。若以残差下降三阶为收敛标准，本问题仅迭代146次即可收敛。 在M150计算机上计算，CPU时间仅需8分多钟。 4、计算流程 本文采用的计算流程框图见图4 12参数和系数的数值见表 片 圣 一一门… 一 — 一 一 — ， 产 了 ’川 叮 闷 自 ， 二二刁尸二 】 护二 ，二 一户丁 七己 】 ， 盆 盈 ，， 了 盈 卜 ， 厂 孟 厂 尸 勺 ’ 厂 图 交错 网格示 意图 图 名 计算用 座 标系示意 图 表 计 算 中 所 用 的 参 数 和 系 数 ， 。 “ 任 “ 泛 又 ‘ ， △ ‘ 功 功 ， ‘ 火 火 数名一 一称值 在使用这些 系数下 ， 尤 其使用 了 。 最佳人工 阻尼 因子 在程序中 自动实 现 ， 本 程序的收敛速率是很快的 。 若以残差下降三 阶为收敛标准 ，本问题仅迭代 次即可 收 敛 。 在 计算机上计算 ， 时间仅需 分 多钟 。 、 计算流程 本文采用 的计算流程框 图见图