L一液体 G一气体 图3边界条件定义图 (2) 在侧壁r=R0<Z<H中=k=3=0 (22) DN高6t52ap-t0g/B p-中P (23) (3) 在对称轴r=0、中=k=司3 Oror 0 (24) （）-c2 (25) ④在自由面2=H:<R◆-8登-8配=0 (26) =(e+p8)/ (H-0.5Na) (27) (5)在气液两相区Z=H0<r<rc Nat 中=N2N:（N,-N2)CN2中2-N2,24s+N1(N:中。-N中) (28) 5=-2、2Ψ3-NΨ2 pr Na2Na3 (Na3-No2) (29) K=N(N:NK3N, K3 (30) 32 (31) -=N/八N4） (32) 8。数值求解方法 对以上拟线性椭圆型偏微分方程及有关代数方程组采用Spalding等〔2〕计算湍流回 流的方法，用差分法数值求解，采用13×13的网格，对某些变址使用了超亚松驰法，以 FORT RANN算法语言编制程序，在M150计算机上，编译约需26秒，计算约需7分钟， 计算收敛标准ε-0.005其表达式为： lp-p-|/Σlm|≤e (33) 式中p-,中、K、3,k一达代次数，三为内点总和。 79图 边 界条件定义 图 在侧壁 二 中 省 二 协 一 中， ， 。 ， ， 在对称 旅击 一 阳击 轴 、 中二 气、 乡 。 。 艺 在 自由面 小 口 一 口口卜八 一 · 井资晶 ‘ 一 ” · 一 ， 在气液两相 区 ’ 中 话一恶 一飞 、一 〔 。 。 劝 ， 一 ’ 。 中 。 、 ， 、 。 、 八 一 工、 。 中 一 。 、卜 〕 占二 甲 一 。 甲 。 。 。 一 。 ， 。 一 ， 。 一 。 旦鱼 一 二 ，。 一 。 。 · 落警 · 一 、 、 坠 一 ，· 一 。 一 一 。 橄值求解方法 对 以上拟线性椭圆型偏微分方程及有关代数方程组 采用 等 〔 〕 计算湍 流 回 流的方法 ， 用差分 法数值求解 ， 采用 火 的网格 ， 对某些变量使 用了 超 亚 松 驰 法 ， 以 算法语言编制程序 ， 在 计算机上 ， 编译约需 秒 ， 计算约需 分 钟 ， 计算收敛标准￡ 其表达式为 万 甲 一 甲 一 ’ 】 艺 甲 ￡ 式 中甲一 右 厂 。 、 、 ，二 ，、、 、 。 ， 、 、 ， 二 。 、 ， 。 。 。 ， ， 甲 、 、 口 长 一 达 划人 戈又， 一 刀 尸 切才 ‘ 、 、 沙 马