多孔底面轴对称二相重力流的数值模拟.pdfVIP

( ) 第 43 卷 　第 3 期 复旦学报 自然科学版 Vol. 43 No. 3 2004 年 6 月 Journal of Fudan University (Natural Science) J un. 2004 　　文章编号 (2004) 多孔底面轴对称二相重 力流的数值模拟 胡　越 , 周　亮 , 张慧生 (复旦大学 力学与工程科学系 , 上海 　200433) 摘 　要 : 对于具有多孔介质底面的轴对称二相重力流 , 引进基于浅水近似的控制方程和相应的边界条件 ,采用 贴体坐标变换使运动边界问题化为固定边界问题 ,提出了基于特征插值并结合使用梯形积分公式和 Newton Raphson 迭代法在时间和空间都具有二阶精度的数值边界条件. 为检验格式的性能和避免编写程序时可能出现 的错误 ,对类似的方程构造了一类精确解. 在空间上采用了二步 Lax 格式、二阶 TVD 格式、三阶 ENO 格式及五 阶 WENO 格式 ,在时间上采用了二阶及三阶的 TVDRungeKutta 方法对该问题进行数值模拟. 数值结果表明 , 在解的光滑区域 ,这几种格式的精度都很高 ,但是在大梯度区 ,二步 Lax 格式将会产生强烈的数值振荡 ,且振荡 不会随网格宽度的减小而减小 ,而其他 3 种格式将不会或仅会产生幅度要小得多的数值振荡 ,且振荡会随网格 宽度的减小而趋向于零. 对实际应用目的来说 ,结合使用二阶 TVDRungeKutta 方法的二阶 TVD 格式是一个 经济而又适当的选择. 关键词 : 数值模拟 ; 轴对称二相重力流 ; 多孔介质底面 中图分类号: O 357. 3 　　　　文献标识码 : A ( ) 当一种密度较大的流体 主要沿水平方向 流入另一种密度较小的流体时就会产生重力流. 这种现象 广泛存在于自然界及工业生产中1 ,2 . 通常情况下 ,成分及温度的差异 ,流体中的颗粒悬浮物 ,或几者的综 合都可能产生重力流3 . 过去对于这一问题的研究主要集中在不可渗透底面上的重力流. 然而 ,有许多 实际问题涉及到重力流体在可渗透底面上流动并伴随流体质量损失的现象. 例如 , 当可燃或有毒流体和 压缩气体从容器中泄漏以后在沙石底面上流动 ,便形成了这种重力流. 真实的物理过程还应包括流体的 蒸发及与周围物质的化学反应4 ,5 . 而对于这类问题最感兴趣的是被沙石底面吸收的质量及流体的传播 距离. 对于这些知识的了解将有助于降低灾害所造成的损失. 典型的问题是在一个放置着较轻流体的容器中突然释放一定体积的大密度流体 ,容器的水平底面是 可渗透的. 现阶段所见到的实验工作似乎只有 Thomas ,Marino 和 Linden 对含盐流体在充满淡水并带有 多孔底面的矩形容器中流动6 . Ungraish 和 Huppert 7 使用两步 Lax 格式8 ～10 对二维问题和轴对称问题 进行了数值模拟. 文献 7 的轴对称情形的结果带有明显的数值振荡. 由于轴对称问题更具有实际应用 价值 ,所以我们应寻找合适的数值格式以减少或避免轴对称情形计算结果中的数值振荡. 1 　控制方程、边界条件和精确解 考虑如图 1 所示的无量纲模型 :在多孔底面 z = 0 上有一层密度为 ρ 的静止的无限流体. 重力沿 z