流场计算中自由表面的处理.pdf

维普资讯 华 北 水 利 水 电 学 院 学 报 1g04年 12月 Jou~ lofNorthCbi∞ [n~ituteofW aterc。 yandHydroelectricPower 流场计算中自由表面 的处理 ／ 马福喜 王克修 王金瑞 ，f／’1 — — — — — — — — — — — — — 一 4／一 f (木刺系) (院办) (^事娃) j一 [提 要] 本文针对水流计算中不规则边界和 自由表面问题 ，从流固两相漉理论出发 ，推 导得出具有不规则边界时的紊流控制方程，引入通度概念处理不规则计算建 莽；采用MAc方法的思路，直接从物质输运方程 出发 ，推得流体体积输运方 程，用VOF法追踪随时问起状的 自由表面。利用所建立的数学模型计算了实 用重力坝厦其下游戽式消力池中水流流场，并对一个收缩式消力池中三维术 跃进行数值模拟 。计算结果较好地与实验结果符舍。 [关键词] 通度、不规则边界、自由表面 前 言 水l蟊嚣 一，f氪响 不规则计算边界和 自由表面问题是计算流体力学中经常遇到的复杂 问题 。对于不规则 边界的处理 ，通常采用阶梯化近似 。但由于它难 以对边界精确而光滑地拟台，对计算精度 有一定影响。为克服这点，美国的ThompsonJ．F．等人提出了边界拟台坐标生成方法 。该 法采用一组偏微分方程将物理边界转换成为简单的计算边界 ，使不规则的物理网格转换成 为规则的计算 网格 。但 由于引^坐标变换 ，使得流体控制方程相应地复杂化，使该法实际 应用受到 了艰 制。本文从流固两相流理论出发，推导得出具有不规则边界时的紊流控制方 程，引人通疫概念处理不规则边界。计算表 明它概念简单，计算程序易于编制，计算结果 明显优于阶梯化近似方法。通常对 自由表面的处理，对于大区域平缓的自由表面水流 问题 采用静压假定或刚盖假定 对于水面起伏较大 的水流采用标高函数法 ，它采用一个高度函 数反映 自由表面的位置，适用于 自由表面变化不太剧烈并且没有重叠的单值 自由表面 本 文引用 MAc方法的思路 ，直接从物质输运方程出发，推得流体体积输运方程，采用VOF 法追踪三维复杂的自由表面，能追踪 自由表面起伏较大或具有重叠的多重 自由表面。 1 控制方程 两相流理论认为，两种连续介质的流动，分别 由各 自的质量 动量、能量守恒方程控 制 n ，其一般形式为 质量守恒方程 — et*,t,i+ ·(R )一 0 (1) 本文于 1994年 5月收到 维普资讯 第4捌 马福喜等 ， 藏墒计算中 自由表酉的处理 动量 守恒方程 1 + (P )一一 O~7P+ [ ( )]+ _-+ (2) , 式中 H——介质的粘性系数； P——两相混合压力 s．——两相之间的相互渗透作用项} ——是 f种介质所 占体积浓度 。 =1，2分别代表两种介质。 具有不规则边界的水流可认为是两相流的特倒。假定两相漉中t一相的各漉动变量在 时间和空间上保持恒定不变，该相即为固体}假定另一相为漉动相．由 (1)和 (2)得不 可压缩流体连续方程和动量方程