A mass-conservative centered finite volume model for solving two-dimensional two-layer shallow water equations for fluid mud propagation over varying topography and dry areas

2012 ◽  
Vol 40 ◽  
pp. 54-70 ◽  
Author(s):  
Alberto Canestrelli ◽  
Sergio Fagherazzi ◽  
Stefano Lanzoni
Keyword(s):  
Shallow Water ◽  
Finite Volume ◽  
Shallow Water Equations ◽  
Fluid Mud ◽  
Two Dimensional ◽  
Finite Volume Model ◽  
Volume Model ◽  
Dry Areas

scholarly journals A Finite-Volume Icosahedral Shallow-Water Model on a Local Coordinate

2009 ◽  
Vol 137 (4) ◽  
pp. 1422-1437 ◽  
Author(s):  
Jin-Luen Lee ◽  
Alexander E. MacDonald
Keyword(s):  
Shallow Water ◽  
Finite Volume ◽  
Spherical Surface ◽  
Second Order ◽  
Water Model ◽  
Numerical Dissipation ◽  
Shallow Water Model ◽  
Finite Volume Model ◽  
Volume Model ◽  
Cartesian Coordinate

Abstract An icosahedral-hexagonal shallow-water model (SWM) on the sphere is formulated on a local Cartesian coordinate based on the general stereographic projection plane. It is discretized with the third-order Adam–Bashforth time-differencing scheme and the second-order finite-volume operators for spatial derivative terms. The finite-volume operators are applied to the model variables defined on the nonstaggered grid with the edge variables interpolated using polynomial interpolation. The projected local coordinate reduces the solution space from the three-dimensional, curved, spherical surface to the two-dimensional plane and thus reduces the number of complete sets of basis functions in the Vandermonde matrix, which is the essential component of the interpolation. The use of a local Cartesian coordinate also greatly simplifies the mathematic formulation of the finite-volume operators and leads to the finite-volume integration along straight lines on the plane, rather than along curved lines on the spherical surface. The SWM is evaluated with the standard test cases of Williamson et al. Numerical results show that the icosahedral SWM is free from Pole problems. The SWM is a second-order finite-volume model as shown by the truncation error convergence test. The lee-wave numerical solutions are compared and found to be very similar to the solutions shown in other SWMs. The SWM is stably integrated for several weeks without numerical dissipation using the wavenumber 4 Rossby–Haurwitz solution as an initial condition. It is also shown that the icosahedral SWM achieves mass conservation within round-off errors as one would expect from a finite-volume model.

Sign in / Sign up

Export Citation Format

Share Document