A new composite scheme for two-layer shallow water flows with shocks. (English) Zbl 1426.76380
Summary: This paper is devoted to solve the system of partial differential equations governing the flow of two superposed immiscible layers of shallow water flows. The system contains source terms due to bottom topography, wind stresses, and nonconservative products describing momentum exchange between the layers. The presence of these terms in the flow model forms a nonconservative system which is only conditionally hyperbolic. In addition, two-layer shallow water flows are often accompanied with moving discontinuities and shocks. Developing stable numerical methods for this class of problems presents a challenge in the field of computational hydraulics. To overcome these difficulties, a new composite scheme is proposed. The scheme consists of a time-splitting operator where in the first step the homogeneous system of the governing equations is solved using an approximate Riemann solver. In the second step a finite volume method is used to update the solution. To remove the non-physical oscillations in the vicinity of shocks a nonlinear filter is applied. The method is well-balanced, non-oscillatory and it is suitable for both low and high values of the density ratio between the two layers. Several standard test examples for two-layer shallow water flows are used to verify high accuracy and good resolution properties for smooth and discontinuous solutions.
MSC:
| 76M12 | Finite volume methods applied to problems in fluid mechanics |
| 76B70 | Stratification effects in inviscid fluids |
Keywords:
two-layer shallow water equations; nonconservative systems; finite volume method; Riemann solver; nonlinear filterReferences:
| [1] | Garratt, J.R.: Review of drag coefficients over oceans and continents. Mon. Weather Rev. 105, 915–929 (1977) · doi:10.1175/1520-0493(1977)105<0915:RODCOO>2.0.CO;2 |
| [2] | Bouchut, F., Morales, T.: An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment. Modél. Math. Anal. Numér. 42, 683–698 (2008) · Zbl 1203.76110 · doi:10.1051/m2an:2008019 |
| [3] | Castro, M.J., Chacón, T., Nieto, E.D., Parés, C.: On well-balanced finite volume methods for nonconservative nonhomogeneous hyperbolic systems. SIAM J. Sci. Comput. 29, 1093–1126 (2007) · Zbl 1148.65065 · doi:10.1137/040607642 |
| [4] | Castro, M.J., Macías, J., Parés, C.: A Q-scheme for a class of systems of coupled conservation laws with source term: application to a two-layer 1D shallow water system. Modél. Math. Anal. Numér. 35, 107–127 (2001) · Zbl 1094.76046 · doi:10.1051/m2an:2001108 |
| [5] | Dudzinski, M., Medvidova, M.L.: Well-balanced path-consistent finite volume EG schemes for the two-layer shallow water equations. Comput. Sci. High. Perform. Comput. IV, 121–136 (2009) |
| [6] | Abgrall, R., Karni, S.: Two-layer shallow water systems: a relaxation approach. SIAM J. Sci. Comput. 31, 1603–1627 (2009) · Zbl 1188.76229 · doi:10.1137/06067167X |
| [7] | Benkhaldoun, F., Elmahi, I., Seaid, M.: A new finite volume method for flux-gradient and source-term balancing in shallow water equations. Comput. Methods Appl. Mech. Eng. 199, 49–52 (2010) · Zbl 1225.76209 · doi:10.1016/j.cma.2010.07.003 |
| [8] | Castro, M.J., LeFloch Muñoz, P.G., Parés, C.: Why many theories of shock waves are necessary: convergence error in formally path-consistent schemes. J. Comput. Phys. 227, 8107–8129 (2008) · Zbl 1176.76084 · doi:10.1016/j.jcp.2008.05.012 |
| [9] | Lee, W.K., Borthwick, A.G.L., Taylor, P.H.: On mathematical balancing of a two-layer shallow flow model. In: 1st IAHR Eur. Div. Congress. IAHR, Edinburgh (2010) |
| [10] | Schijf, J.B., Schonfeld, J.C.: Theoretical considerations on the motion of salt and fresh water. In: Proc. of the Minn. Int. Hydraulics Conv., Joint Meeting IAHR and Hyd. Div. ASCE, pp. 321–333 (1953) |
| [11] | Macías, J., Parés, C., Castro, M.J.: Improvement and generalization of a finite element shallow water solver to multi-layer systems. Int. J. Numer. Methods Fluids 31, 1037–1059 (1999) · Zbl 0969.76045 · doi:10.1002/(SICI)1097-0363(19991215)31:7<1037::AID-FLD909>3.0.CO;2-V |
| [12] | Castro, M.J., Macías, J., Parés, C., Garciá-Rodríguez, J.A., Vázquez-Cendoón, E.: A two-layer finite volume model for flows through channels with irregular geometry: computation of maximal exchange solutions. Application to the Strait of Gibraltar. Commun. Nonlinear Sci. Numer. Simul. 9, 241–249 (2004) · Zbl 1136.76394 · doi:10.1016/S1007-5704(03)00115-1 |
| [13] | Liska, R., Wendroff, B.: Analysis and Computation with Stratified Fluid Models. J. Comput. Phys. 137, 212–244 (1997) · Zbl 0898.76076 · doi:10.1006/jcph.1997.5806 |
| [14] | Audusse, E., Benkhaldoun, F., Sainte-Marie, J., Seaid, M.: Multilayer Saint-Venant equations over movable beds. Discrete Contin. Dyn. Syst., Ser. B 15, 917–934 (2011) · Zbl 1230.37102 · doi:10.3934/dcdsb.2011.15.917 |
| [15] | Ambrosi, D.: Approximation of shallow water equations by Riemann solvers. Int. J. Numer. Methods Fluids 20, 157–168 (1995) · Zbl 0831.76063 · doi:10.1002/fld.1650200205 |
| [16] | Sahmim, S., Benkhaldoun, F., Alcrudo, F.: A sign matrix based scheme for quasi-hyperbolic non-homogeneous PDEs with an analysis of the convergence stagnation problem. J. Comput. Phys. 226, 1753–1783 (2007) · Zbl 1173.76367 · doi:10.1016/j.jcp.2007.06.017 |
| [17] | Engquist, B., Ltstedt, P., Sjgreen, B.: Nonlinear filters for efficient shock computation. Math. Comput. 186, 509–537 (1989) · doi:10.1090/S0025-5718-1989-0955750-9 |
| [18] | Spinewine, B., Guinot, V., Soares Frazão, S., Zech, Y.: Solution properties and approximate Riemann solvers for two-layer shallow water flow models. Math. Comput. 186, 509–537 (1989) · Zbl 1271.76194 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
