任意曲线坐标系下二维浅水方程的数值模拟.pdf

２０１２年９月 第３０卷　第５期 　 Ｓｅｐｔ．２０１２ Ｖｏｌ．３０　 Ｎｏ．５ ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７４－８５３０．２０１２．０５．０２１ 任意曲线坐标系下二维浅水方程的数值模拟 吴卫国，薛世峰 （中国石油大学（华东）工程力学系，山东 青岛２６６５５５） 摘要：采用拟合边界曲线的方法来构造天然河道、湖泊、港口河口以及海湾的复杂边界问题，建 立任意正交曲线网格以克服由于复杂边界而引起的计算困难．在此基础上，推导出任意曲线坐 标系下的二维浅水方程、湍流动能方程和湍流动能耗散方程；应用有限差分方法对该方程组进 行数值离散，并用交错方向隐式格式实现在计算区域内对任意曲线坐标系下的二维浅水方程进 行数值求解．为了验证在任意曲线坐标系下二维浅水方程数值求解方法的可靠性、正确性，以Ｄｅ Ｖｒｉｅｎｄ的１８０°平面弯道水槽试验物理模型为例进行数值模拟，结果表明，数值计算的结果和Ｄｅ －２ Ｖｒｉｅｎｄ的试验结果相当吻合，最大绝对误差值约为１０ ，因此，数值计算方法合理可行，可为任 意复杂边界的天然河道、湖泊等水域的水动力研究提供有效的计算方法． 关键词：浅水方程；曲线坐标；交错方向隐式格式；复杂边界；数值模拟 中图分类号：Ｓ２７７９；ＴＶ１４５２１　文献标志码：Ａ　文章编号：１６７４－８５３０（２０１２）０５－０５９８－０５ Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆ２Ｄｓｈａｌｌｏｗｗａｔｅｒｅｑｕａｔｉｏｎｉｎ ａｒｂｉｔｒａｒｙｃｕｒｖｉｌｉｎｅａｒｃｏｏｒｄｉｎａｔｅｓ ＷｕＷｅｉｇｕｏ，ＸｕｅＳｈｉｆｅｎｇ （ＤｅｐａｒｔｍｅｎｔｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＣｈｉｎａＵｎｉｖｅｒｓｉｔｙｏｆＰｅｔｒｏｌｅｕｍ，Ｑｉｎｇｄａｏ，Ｓｈａｎｄｏｎｇ２６６５５５，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅｐｒｏｂｌｅｍｓｗｉｔｈｃｏｍｐｌｅｘｂｏｕｎｄａｒｙｓｈａｐｅｓ，ｓｕｃｈａｓｎａｔｕｒａｌｒｉｖｅｒｃｈａｎｎｅｌｓ，ｌａｋｅｓ，ｅｓｔｕａ ｒｉｅｓａｎｄｂａｙｓ，ｗｅｒｅｓｏｌｖｅｄｂｙｕｓｉｎｇｔｈｅｂｏｕｎｄａｒｙｃｕｒｖｅｆｉｔｔｉｎｇｍｅｔｈｏｄｓ．Ａｒｂｉｔｒａｒｙｏｒｔｈｏｇｏｎａｌｃｕｒｖｉｌｉ ｎｅａｒｇｒｉｄｗａｓｅｓｔａｂｌｉｓｈｅｄｔｏｏｖｅｒｃｏｍｅｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｄｉｆｆｉｃｕｌｔｉｅｓｃａｕｓｅｄｂｙｔｈｏｓｅｃｏｍｐｌｅｘｂｏｕｎｄａｒｉｅｓ． Ｔｈｅｎ，ａｓｅｔｏｆｄｅｒｉｖｅｄｅｑｕａｔｉｏｎｓｉｎｔｈｅａｒｂｉｔｒａｒｙｃｕｒｖｉｌｉｎｅａｒｃｏｏｒｄｉｎａｔｅｓ，ｉｎｃｌｕｄｉｎｇ２Ｄｓｈａｌｌｏｗｗａｔｅｒｅ ｑｕａｔｉｏｎ，ｔｕｒｂｕｌｅｎｃｅｋｉｎｅｔｉｃｅｎｅｒｇｙｅｑｕａｔｉｏｎａｎｄｄｉｓｓｉｐａｔｉｏｎｒａｔｅｅｑｕａｔｉｏｎｅｔｃ．ｗｅｒｅｎｕｍｅｒｉｃａｌｌｙｄｉｓｃｒｅ ｔｅｄｂｙｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄ．Ｉｎａｄｄｉｔｉｏｎ，ｔｈｅ２Ｄｓｈａｌｌｏｗｗａｔｅｒｅｑｕａｔｉｏｎｗａｓｎｕｍｅｒｉｃａｌｌｙｓｏｌｖｅｄ ｗｉｔｈｉｎｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｄｏｍａｉｎｂｙｕｓｉｎｇｔｈｅａｌｔｅｒｎａｔｉｎｇｄｉｒｅｃｔｉｏｎｉｍｐｌｉｃｉｔ（ＡＤＩ）ｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅ． Ｉｎｏｒｄｅｒｔｏｖｅｒｉｆｙｔｈｅｒｅｌｉａｂｉｌｉｔｙａｎｄｃｏｒｒｅｃｔｎｅｓｓｏｆｔｈｅｍｅｔｈｏｄ，ｔｈｅＤｅＶｒｉｅｎｄ＇ｓ１８０°ｐｌａｎｅｃｕｒｖｅｆｌｕｍｅ