基于无结构化网格浅水方程的隐式解法-人民长江.PDF

第４６卷 第５期 人　民　长　江 Ｖｏｌ．４６，Ｎｏ．５ ２０１５年 ３月　 　 Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　 　 Ｍａｒ．， ２０１５ 　　文章编号：１００１－４１７９（２０１５）０５－００８１－０４ 基于无结构化网格浅水方程的隐式解法 唐岳 灏 （长江工程职业技术学院，湖北 武汉４３０２１２） 摘要：为提高有限体积法计算浅水方程的数值稳定性，采用Ｒｏｅ方法近似 Ｒｉｅｍａｎｎ解计算界面通量，利用ＴＶＤ －ＭＵＳＣＬ格式对守恒变量进行重构，推导并建立了高效的隐式计算格式。该格式基于无结构化网格，将计算 精度提高到二阶，使用面积加权计算流速梯度，对坡底梯度的处理满足静水平衡要求。为在时间积分上采用 隐格式，又导出了控制方程的雅可比矩阵全解析形式，并采用牛顿 －拉夫逊方法进行迭代求解。通过对溃坝 算例的对比分析，验证该方法的稳定性、和谐性和高效性，同时也证明该算法具有捕捉溃坝激波问题的能力。 关　键　词：浅水方程；Ｒｏｅ格式；非结构网格；有限体积法；隐格式 中图法分类号：ＴＶ１３５　　　文献标志码：Ａ ＤＯＩ：１０．１６２３２／ｊ．ｃｎｋｉ．１００１－４１７９．２０１５．０５．０１９ 　　用有限体积法计算浅水方程时，目前绝大部分基  Ｉ Ｖ Ｕｄ ＋ （Ｆ － Ｆ ）ｄＳ＝ Ｈｄ （１） Ω ν Ω ∫ ∮ ｔ ∫ ｔ 于无结构网格的数值方法都使用显格式。显格式比隐 Ω Ｓ Ω ｕ ｕ ｖ 格式在程序编制上要简单，但由于受到 ＣＦＬ条件限 Ｔ Ｖ    式中，Ｕ＝［ｈ，ｈｕ，ｈｖ］ ；Ｆ ＝［０，ｈ ｎ＋ｈ ｎ，ｈ ｘ ｙ ｘ ｙ ｘ    制，计算的时间步长一般较小，计算效率受到影响。另 ｖ  Ｔ Ｔ 外，由于源项的非线性性，显格式的数值稳定性也受到 ｎ＋ｈ ｎ］；Ｈ＝［０，ｇｈ（Ｓ －Ｓ），ｇｈ（Ｓ －Ｓ）］； ｘ ｙ ０ｘ ｆｘ ０ｙ ｆｙ ｙ  影响［１