带源项二维浅水方程正问题和反问题的求解 - 西安理工大学学报.pdf

　西安理工大学学报 ＪｏｕｒｎａｌｏｆＸｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ（２０１１）Ｖｏｌ．２７Ｎｏ．１　 １０７ 　　文章编号：１００６４７１０（２０１１）０１０１０７０５ 带源项二维浅水方程正问题和反问题的求解 １ ２ １ １ 刘文涛 ，闵涛 ，郝竹林 ，冯民权 （１．西安理工大学 水利水电学院，陕西 西安７１００４８；２．西安理工大学 理学院，陕西 西安７１００５４） 摘要：考虑到求解二维浅水方程的正问题及反问题的复杂性，建立高分辨率的有限元格式对正问 题进行模拟，运用最佳摄动量法对反问题进行了研究，通过对两类溃坝问题进行数值模拟及文献 对比，验证了有限元方法具有较高的离散精度同时又避免了数值解的伪震荡，并运用最佳摄动量 法对正问题算例进行糙率率定，模拟结果较为理想。 关键词：浅水方程；有限元法；最佳摄动量法；糙率率定 中图分类号：ＴＶ１３９．２　　　文献标志码：Ａ  ＴｈｅＳｏｌｕｔｉｏｎｔｏＦｏｒｗａｒｄａｎｄＩｎｖｅｒｓｅＰｒｏｂｌｅｍｓｏｆ２Ｄ ． ＳｈａｌｌｏｗＷａｔｅｒＥｑｕａｔｉｏｎｗｉｔｈＳｏｕｒｃｅＴｅｒｍｓ １ ２ １ １ ＬＩＵＷｅｎｔａｏ，ＭＩＮＴａｏ，ＨＡＯＺｈｕｌｉｎ，ＦＥＮＧＭｉｎｑｕａｎ （１．ＦａｃｕｌｔｙｏｆＷａｔｅｒＲｅｓｏｕｒｃｅｓａｎｄＨｙｄｒｏｅｌｅｃｔｒｉｃＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｘｉ’ａｎ７１００４８，Ｃｈｉｎａ； ２．ＦａｃｕｌｔｙｏｆＳｃｉｅｎｃｅ，Ｘｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｘｉ’ａｎ７１００５４，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｉｎｃｏｎｓｉｄｅｒｉｎｇｔｈｅｃｏｍｐｌｅｘｉｔｙｏｆｓｏｌｖｅｄｐｒｏｂｌｅｍｓａｎｄｉｎｖｅｒｓｅｐｒｏｂｌｅｍｓｏｆｔｈｅｔｗｏｄｉｍｅｎｓｉｏｎａｌ ｓｈａｌｌｏｗｗａｔｅｒｅｑｕａｔｉｏｎｓ，ａｈｉｇｈｒｅｓｏｌｕｔｉｏｎｆｉｎｉｔｅｅｌｅｍｅｎｔｓｃｈｅｍｅｉｓｅｓｔａｂｌｉｓｈｅｄｔｏｓｉｍｕｌａｔｅｔｈｅｆｏｒｗａｒｄ ｐｒｏｂｌｅｍ，ａｎｄｔｈｅｏｐｔｉｍａｌｐｅｒｔｕｒｂａｔｉｏｎｍｅｔｈｏｄｉｓｕｓｅｄｔｏｓｔｕｄｙｔｈｅｉｎｖｅｒｓｅｐｒｏｂｌｅｍ．Ｔｈｅｆｉｎｉｔｅｅｌｅｍｅｎｔ ｍｅｔｈｏｄｗｉｔｈｈｉｇｈｄｉｓｃｒｅｔｅａｃｃｕｒａｃｙａｎｄｐｓｅｕｄｏｓｈｏｃｋｔｏａｖｏｉｄｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｉｓｔｅｓｔｅｄｖｉａｔｈｅｎｕｍｅｒｉ ｃａｌｓｉｍｕｌａｔｉｏｎａｎｄｌｉｔｅｒａｔｕｒｅｃｏｎｔｒａｓｔｏｆｄａｍｂｒｅａｋｐｒｏｂｌｅｍｓｏｆｔｗｏｔｙｐｅｓ；ａｎｄｔｈｅｏｐｔｉｍａｌｐｅｒｔｕｒｂａｔｉｏｎ ｍｅｔｈｏｄｉｓｕｓｅｄｔｏｃａｌｉｂｒａｔｅｔｈｅｒｏｕｇｈｎｅｓｓｏｆｔｈｅａｌｇｏｒｉｔｈｍｃａｓｅｓｏｆｔｈｅｆｏｒｗａｒｄｐｒｏｂｌｅｍ．Ａｃｃｏｒｄｉｎｇｌｙ，ｔｈｅ ｓｉｍｕｌａｔｅｄｒｅｓｕｌｔｓａｒｅｒｅａｌｉｓｔｉｃ． Ｋｅｙｗｏｒｄｓ：ｓｈａｌｌｏｗｗａｔｅｒｅｑｕａｔｉｏｎ；ｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄ；ｏｐｔｉｍａｌｐｅｒｔｕｒｂａｔｉｏｎｍｅｔｈｏｄ；ｒｏｕｇｈｎｅｓｓｉ ｄｅｎｔｉｆｉｃａｔｉｏｎ 　　浅水方程属于非线性双曲型方程，不管初始值 ｏｄ）来计算流体时，其速度和压力需要以不同阶数的 ［１］ 如何光滑，其解都可能产生间断 。浅水方程正问 形函数（ｓｈａｐｅｆｕｎｃｔ