Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-10T21:02:07.001Z Has data issue: false hasContentIssue false

A Q-scheme for a class of systems of coupled conservation laws with source term. Application to a two-layer 1-D shallow water system

Published online by Cambridge University Press:  15 April 2002

Manuel Castro
Affiliation:
Dpto. Análisis Matemático. Facultad de Ciencias. Universidad de Málaga, Campus de Teatinos s/n, 29080 Málaga, Spain. (grupo@anamat.cie.uma.es)
Jorge Macías
Affiliation:
Dpto. Análisis Matemático. Facultad de Ciencias. Universidad de Málaga, Campus de Teatinos s/n, 29080 Málaga, Spain. (grupo@anamat.cie.uma.es)
Carlos Parés
Affiliation:
Dpto. Análisis Matemático. Facultad de Ciencias. Universidad de Málaga, Campus de Teatinos s/n, 29080 Málaga, Spain. (grupo@anamat.cie.uma.es)
Get access

Abstract

The goal of this paper is to construct a first-order upwind schemefor solving the system of partial differential equations governing theone-dimensional flow of two superposed immiscible layers of shallow waterfluids. This is done by generalizing a numerical scheme presented byBermúdez and Vázquez-Cendón [3, 6, 27] for solving one-layer shallow water equations, consisting in a Q-scheme with a suitable treatment of the source terms.The difficulty in the two layer system comes from the coupling terms involving some derivatives of the unknowns.Due to these terms, a numerical scheme obtained by performing the upwinding of each layer, independently from the other one, can be unconditionally unstable.In order to define a suitable numerical scheme with global upwinding,we first consider an abstract system that generalizes the problem under study.This system is not a system of conservation laws but, nevertheless,Roe's method can be applied to obtain an upwind scheme based on Approximate Riemann State Solvers. Following this, we present some numerical tests to validate the resulting schemes and to highlight the fact that, in general, numerical schemes obtained by applying a Q-scheme to each separate conservation law of the system do not yield good results. First, a simple system of coupled Burgers' equations is considered. Then, the Q-scheme obtained is applied to the two-layer shallow water system.

Type
Research Article
Copyright
© EDP Sciences, SMAI, 2001

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Armi, L., The hydraulics of two flowing layers with different densities. J. Fluid Mech. 163 (1986) 27-58. CrossRef
Armi, L. and Farmer, D., Maximal two-layer exchange through a contraction with barotropic net flow. J. Fluid Mech. 164 (1986) 27-51. CrossRef
Bermúdez, A. and Vázquez, M.E., Upwind methods for hyperbolic conservation laws with source terms. Computers and Fluids 23 (1994) 1049-1071.
C. Berthon and F. Coquel, Travelling wave solutions of a convective diffusive system with first and second order terms in nonconservation form, in Hyperbolic Problems: Theory, Numerics, Applications, Vol. I of Internat. Ser. Numer. Math. 129, Birkhäuser (1999) 47-54.
M.J. Castro, J. Macías and C. Parés, Simulation of two-layer exchange flows through a contraction with a finite volume shallow water model, in Actas de las II Jornadas de Análisis de Variables y Simulación Numérica del Intercambio de Masas de Agua a través del Estrecho de Gibraltar, Cádiz (2000) 205-221.
M.J. Castro, J. Macías and C. Parés, Simulation of two-layer exchange flows through the combination of a sill and contraction with a finite volume shallow-water model. Internal Journal 1610, group on ``Differential Equations, Numerical Analysis and Applications'', University of Málaga (2000).
Coquel, F., El Amine, K., Godlewski, E., Perthame, B. and Rascle, P., Une méthode numérique decentrée pour la résolution d'ecoulements diphasiques. C. R. Acad. Sci. Paris, Sér. I 324 (1997) 717-723.
S.B. Dalziel, Two-layer Hydraulics Maximal Exchange Flows. Ph.D. thesis, University of Cambridge (1988).
Farmer, D. and Armi, L., Maximal two-layer exchange over a sill and through a combination of a sill and contraction with barotropic flow. J. Fluid Mech. 164 (1986) 53-76. CrossRef
P. García-Navarro and F. Alcrudo, Implicit and explicit TVD methods for discontinuous open channel flows, in Proc. of the 2nd Int. Conf. on Hydraulic and Environmental Modelling of Coastal, Estuarine and River Waters, R.A. Falconer, K. Shiono, and R.G.S. Matthew, Eds. 2 Ashgate (1992).
García-Navarro, P., Alcrudo, F. and Savirón, J.M., 1D open channel flow simulation using TVD McCormack scheme. J. Hydraul. Eng. 118 (1992) 1359-1373.
A.E. Gill, Atmosphere-Ocean Dynamics, Int. Geophys. Series 30, Springer-Verlag, San Diego (1982) 662 p.
E. Godlewski and P.A. Raviart, Numerical Approximation of Hyperbolic Systems of Conservation Laws, Appl. Math. Sc. 118, Springer-Verlag, New York (1996).
Harten, A., On a class of high resolution total-variation-stable finite-difference schemes. SIAM J. Numer. Anal. 21 (1984) 1-23. CrossRef
Harten, A., Lax, P. and van Leer, A., On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev. 25 (1983) 35-61. CrossRef
Helfrich, K.R., Time-dependent two-layer hydraulic exchange flows. J. Phys. Oceanogr. 25 (1995) 359-373. 2.0.CO;2>CrossRef
P.K. Kundu, Fluid Mechanics. Academic Press Inc., San Diego (1990) 638 p.
Lefloch, P.G. and Tzavaras, A.E., Existence theory for the Riemann problem for non-conservative hyperbolic systems. C. R. Acad. Sci. Paris, Sér. I 323 (1996) 347-352.
Lefloch, P.G. and Tzavaras, A.E., Representation of weak limits and definition of nonconservative products. SIAM J. Math. Anal. 30 (1999) 1309-1342. CrossRef
Roe, P.L., Approximate Riemann solvers, parameter vectors and difference schemes. J. Comput. Phys. 43 (1981) 357-371. CrossRef
P.L. Roe, Upwinding differenced schemes for hyperbolic conservation laws with source terms, in Proc. of the Conference on Hyperbolic Problems, C. Carasso, P.-A. Raviart, and D. Serre, Eds., Springer-Verlag, Berlin (1986) 41-51.
J.A. Rubal and M.E. Vázquez, Aplicación del método de volúmenes finitos y esquemas tipo Godunov a un modelo bicapa, in Actas de las II Jornadas de Análisis de Variables y Simulación Numérica del Intercambio de Masas de Agua a través del Estrecho de Gibraltar, Cádiz (2000) 223-239.
J.B. Schijf and J.C. Schonfeld, 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., Sept. 1953 (1953) 321-333.
J.J. Stoker, Water Waves. Interscience, New York (1957).
E.F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics. A Practical Introduction. Springer-Verlag, Berlin (1997).
M.E. Vázquez-Cendón, Estudio de Esquemas Descentrados para su Aplicación a las leyes de Conservación Hiperbólicas con Términos Fuente. PhD thesis, Universidad de Santiago de Compostela (1994).
Vázquez-Cendón, M.E., Improved treatment of source terms in upwind schemes for the shallow water equations in channels with irregular geometry. J. Comp. Physics 148 (1999) 497-526. CrossRef