一种godunov算法及其在非平底流污染物传输中的-长江科学院院报.pdfVIP

　 第３０卷第８期 长　江　科　学　院　院　报 Ｖｏｌ．３０　Ｎｏ．８　　 　　２０１３年 ８月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ａｕｇ．２０１３　 　 　　ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００１－５４８５．２０１３．０８．０２１ ２０１３，３０（０８）：９７－１０１ 一种 Ｇｏｄｕｎｏｖ算法及其在非平底流污染物 传输中的应用 １ １ １ ２ 杨金波 ，段文刚 ，黄国兵 ，李订芳 （１．长江科学院水力学研究所，武汉４３００１０；２．武汉大学数学与统计学院，武汉４３００７２） 摘要：采用和谐的加权平均通量（ＷｅｉｇｈｔｅｄＡｖｅｒａｇｅＦｌｕｘ，简称ＷＡＦ）算法，研究了浅水波方程的间断解及污染物传输 问题。该算法采用ＷＡＦ格式和ＨＬＬＣＲｉｅｍａｎｎ求解器近似单元边界数值通量，中心差分格式离散地形源项，然后理 论上证明了该算法是和谐的。最后利用ＷＡＦ算法对非平底地形上浅水波间断解及污染物传输问题进行数值计算， 精确地捕捉到了间断解和污染物运动过程，结果表明该算法满足守恒性，具有高分辨率、无振荡及捕捉污染物运动边 界的能力。 关　键　词：浅水波方程；ＷＡＦ格式；溃坝；污染物传输 中图分类号：Ｘ５２　　　文献标志码：Ａ　　　文章编号：１００１－５４８５（２０１３）０８－００９７－０５ 非平底河床上浅水流的间断解及污染物传输问题， 将计算结果同精确解或其他文献比对，证明该方法 １　研究背景 的计算结果精度高、分辨率高，能计算真实的溃坝和 河流、湖泊和近海岸的污染问题越来越严重，一 污染物传输问题。 方面破坏了生态环境；另一方面对人类生命安全和 经济发展带来致命性的危害，因此有必要研究污染 ２　数值格式 物传输方式的数学模型。 污染物传输的数值算法包括浅水波和污染物传 ２．１　浅水波及污染物传输方程 输扩散的计算。在浅水波方程计算方面，近年来空 控制方程由浅水波方程和污染物运输方程组 气动力学中的 ［１］ 成，其中，污染物部分主要研究污染物传输过程，忽 Ｇｏｄｕｎｏｖ 型有限体积算法越来越普 ［２］ 略污染物扩散项，如 遍地应用到浅水波方程中，ＷＡＦ方法 是 Ｔｏｒｏ基 Ｕ Ｆ   于Ｇｏｄｕｎｏｖ思想提出的显示、二次精度的守恒有限 ＋ ＝Ｓ　。 （１） ｔ ｘ   体积格式。非平底河床是浅水方程计算的一个关键 其中： ［３］ 问题，为了构造和谐的数值格式，Ｒｏｅ 提出迎风离 Ｈ  ｑ    散双曲型方程中源项的思想；Ｂｅｒｍｕｄｅｚ