Solution of shallow-water equations using least-squares finite-element method

• 来源期刊: 力学学报(英文版)
• 年，卷（期）: 2008,（5）
• 页码范围: 523-532
• 页数: 10页
• 语种: 英文
• DOI: 10.1007/s10409-008-0151-4

摘要

A least-squares finite-element method (LSFEM) for the non-conservative shallow-water equations is pre-sented. The model is capable of handling complex topogra-phy, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include: (1) sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; (2) upwind scheme is no needed; (3) a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) condition; and (4) the resulting system of equations is sym-metric and positive-definite (SPD) which can be solved effi-ciently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylin-der and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.