移动边界浅水问题的数值研究ξ.pdfVIP

　A 辑第 17 卷第 4 期 　　　　　水 动 力 学 研 究 与 进 展 　　　Ser. A , Vol. 17 ,No. 4 　　2002 年 8 月 　　　　　　　JOURNAL OF HYDRODYNAMICS 　　　　　Aug. , 2002 ( ) 文章编号 2002 移动边界浅水问题的数值研究 毛献忠 , 　　潘存鸿 (浙江省水利水电河口海岸研究设计院 ,浙江杭州 310020) 　　摘 　要 : 　移动边界问题是水力计算中难点之一 ,本文用拉格朗日坐标系来描述浅水方程 ,并 用二阶 Godunov 算法求解有移动边界的浅水问题的数值解。在拉格朗日坐标系下 ,移动边界条件 十分容易处理 ,数值计算结果表明 ,该方法十分有效。 关 　键 　词 : 　浅水方程 ;移动边界 ;拉格朗日坐标 ; Godunov 格式 中图分类号: 　TV133 　　　文献标识码 :A 1 　引言 目前,不恒定流的数值计算是解决水力学工程实际应用问题中重要的研究手段之一 ,有限 差分法、有限元法和有限体积法是水力学计算中广泛采用的计算方法。对于有间断的流动 ,象 溃坝波、水跃等水力学现象 ,近二十多年来在空气动力学中提出的各种高性能计算格式已经广 ( ) [1 ] 泛应用到这些混合流 缓流和急流 问题计算中,并十分成功模拟这种混合流动现象 。 移动边界问题是水力计算中常见的边界条件 ,象干河床上的瞬时溃坝波问题、干河床上的 洪水演进问题、潮流涨落变化引起的陆地动边界问题、以及波浪爬高问题等等 ,都是移动边界 问题。目前处理这一问题最常用的方法是假设干地区域的初始条件被一层水深很小的薄水层 所覆盖 ,这种处理方法简单实用 ,但是它要求计算域必须覆盖到洪水可能到达的所有地方 ,这 在计算中造成很多额外的计算网格和计算时间 ,而且 ,对于某些问题 ,如溃坝波问题 ,这种处理 方法会导致完全不同的物理现象 ,文献 [2 ] 采用耗散有限元方法 ,并用网格变形技术及引入 计算边界条件计算洪水演进前锋面 ,计算过程中不断调整计算网格 ,这种网格的调整不但会浪 费计算时间 ,更重要的是新网格的信息是通过旧网格上的信息插值而得 ,这个过程势必会丧失 部分计算精度。 以上这二种移动边界的处理方法的缺点是欧拉坐标系下求解移动边界问题中无法避免 的 ,因为在欧拉坐标系下 ,计算网格在空间上是固定的 ,它不能象拉格朗日坐标系下的计算网 格随着流体质点的运动而随之相应变化。这就启发我们如果采用拉格朗日坐标系来求解移动 边界的浅水问题 ,移动边界问题就可以迎刃而解。 ( ) 把欧拉坐标系和拉格朗日坐标系统一起 文献 [3 ] 引入统一坐标系 ,用参量 h 0 ≤h ≤1 收稿日期 : ( ) 基金项目 : 　国家自然科学基金资助项目 ( ) 作者简介 : 　毛献忠 1968～ ,男 ,高级工程师 ,博士。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 508 水 　动 　力 　学 　研 　究 　与 　进 　展 　　　　　　　　　2002 年第 4 期 来, 当 h = 0 时为欧拉坐标系, 当 h = 1 时为拉格朗日坐标系, 并应用统一坐标系成功模拟二 ( ) 维浅水问题。本文采用文献 [