High resolution finite-element analysis of shallow water equations in two dimensions

Abstract We present in this study the Taylor–Galerkin finite-element model to stimulate shallow water equations for bore wave propagation in a domain of two dimensions. To provide the necessary precision for the prediction of a sharply varying solution profile, the generalized Taylor–Galerkin finite-element model is constructed through introduction of four parameters. This paper also presents the fundamental theory behind the choice of free parameters. One set of parameters is theoretically determined to obtain the high-order accurate Taylor– Galerkin finite-element model. The other set of free parameters is determined using the underlying discrete maximum principle to obtain the low-order monotonic Taylor–Galerkin finite-element model. Theoretical study reveals that the higher-order scheme exhibits dispersive errors near the discontinuity while lower-order scheme dissipates the discontinuity. A scheme which has a high-resolution shock-capturing ability as a built-in feature is, thus, needed in the present study. Notice that lumping of the mass matrix equations is invoked in the low-order scheme to allow simulation of the hydraulic problem with discontinuities. We check the prediction accuracy against suitable test problems, preferably ones for which exact solutions are available. Based on numerical results, it is concluded that the Taylor–Galerkin-flux-corrected transport (TG-FCT) finite-element method can render the technique suitable for solving shallow water equations with sharply varying solution profiles.

[1]  Derek M. Causon,et al.  HIGH-RESOLUTION FINITE-VOLUME METHOD FOR SHALLOW WATER FLOWS , 1998 .

[2]  Tony W. H. Sheu,et al.  A NUMERICAL STUDY OF NONLINEAR PROPAGATION OF DISTURBANCES IN TWO-DIMENSIONS , 1996 .

[3]  E. F. Toro,et al.  Riemann problems and the WAF method for solving the two-dimensional shallow water equations , 1992, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[4]  H. Paillere,et al.  Multidimensional upwind schemes for the shallow water equations , 1998 .

[5]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[6]  Thomas Molls,et al.  Space-Time Conservation Method Applied to Saint Venant Equations , 1998 .

[7]  Reinaldo García,et al.  Numerical solution of the St. Venant equations with the MacCormack finite‐difference scheme , 1986 .

[8]  Xu-Dong Liu,et al.  Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes , 1998, SIAM J. Sci. Comput..

[9]  Pilar García-Navarro,et al.  Genuinely multidimensional upwinding for the 2D shallow water equations , 1995 .

[10]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[11]  M. Hsu,et al.  Numerical simulation of supercritical shock wave in channel contraction , 1995 .

[12]  Guillermo Hauke,et al.  A symmetric formulation for computing transient shallow water flows , 1998 .

[13]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[14]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[15]  Theodor Meis,et al.  Numerical solution of partial differential equations , 1981 .

[16]  N. Katopodes,et al.  Computing Two-Dimensional Dam-Break Flood Waves , 1978 .

[17]  J. Donea.  Generalized Galerkin Methods for Convection Dominated Transport Phenomena , 1991 .

[18]  D. Zhao,et al.  Approximate Riemann solvers in FVM for 2D hydraulic shock wave modeling , 1996 .

[19]  S. Sengupta,et al.  Boundary conformed co‐ordinate systems for selected two‐dimensional fluid flow problems. Part I: Generation of BFGs , 1986 .

[20]  R. F. Warming,et al.  The modified equation approach to the stability and accuracy analysis of finite-difference methods , 1974 .

[21]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[22]  J. Donea.  A Taylor–Galerkin method for convective transport problems , 1983 .

[23]  R. Maccormack.  The Effect of Viscosity in Hypervelocity Impact Cratering , 1969 .

[24]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[25]  M. Hanif Chaudhry,et al.  Explicit Methods for 2‐D Transient Free Surface Flows , 1990 .

[26]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[27]  Eleuterio F. Toro,et al.  Experimental and numerical assessment of the shallow water model for two-dimensional dam-break type problems , 1995 .