Numerical simulation of two-layer shallow water flows through channels with irregular geometry

This paper deals with the numerical simulation of flows of stratified fluids through channels with irregular geometry. Channel cross-sections are supposed to be symmetric but not necessarily rectangular. The fluid is supposed to be composed of two shallow layers of immiscible fluids of constant densities, and the flow is assumed to be one-dimensional. Therefore, the equations to be solved are a coupled system composed of two Shallow Water models with source terms involving depth and breadth functions. Extensions of the Q-schemes of van Leer and Roe are proposed where a suitable treatment of the coupling and source terms is performed by adapting the techniques developed in [J. Comput. Phys. 148 (1999) 497; Comput. Fluids 29(8) (2000) 17; Math. Model. Numer. Anal. 35(1) (2001) 107]. An enhanced consistency condition, the so-called C-property, introduced in [Comput. Fluids 23(8) (1994) 1049] is extended to this case and a general result providing sufficient conditions to ensure this property is shown. Then, some numerical tests to validate the resulting schemes are presented. First we verify that, in practice, the numerical schemes satisfy the C-property , even for extremely irregular channels. Then, in order to validate the schemes, we compare some approximate steady solutions obtained with the generalized Q-scheme of Van Leer with those obtained by using the asymptotic techniques developed by Armi and Farmer for channels with simplified geometries. Finally we apply the numerical scheme to the simulation of the flow through the Strait of Gibraltar. Real bathymetric and coast-line data are considered to include in the model the main features of the abrupt geometry of this natural strait connecting the Atlantic Ocean and the Mediterranean Sea. A steady-state solution is obtained from lock-exchange initial conditions. This solution is then used as initial condition to simulate the main semidiurnal and diurnal tidal waves in the Strait of Gibraltar through the imposition of suitable boundary conditions obtained from observed tidal data. Comparisons between numerical results and observed data are also presented.

[1]  L. Armi.  The hydraulics of two flowing layers with different densities , 1986, Journal of Fluid Mechanics.

[2]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[3]  H. Bryden,et al.  Exchange through the Strait of Gibraltar , 1994 .

[4]  D. Farmer,et al.  The flow of Atlantic water through the Strait of Gibraltar , 1988 .

[5]  María Elena Vázquez Cendón,et al.  Estudio de esquemas descentrados para su aplicación a las leyes de conservación hiperbólicas con términos fuente , 1994 .

[6]  M. Vázquez-Cendón.  Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .

[7]  R. LeVeque.  Approximate Riemann Solvers , 1992 .

[8]  J. B. Schijf,et al.  Theoretical considerations on the motion of salt and fresh water , 1953 .

[9]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[10]  E. F. Toro,et al.  Model Hyperbolic Systems with Source Terms: Exact and Numerical Solutions , 2001 .

[11]  G. Lawrence.  The hydraulics of steady two-layer flow over a fixed obstacle , 1993, Journal of Fluid Mechanics.

[12]  J. Backhaus,et al.  Control Variability and Internal Bore Evolution in the Strait of Gibraltar: A 2-D Two-Layer Model Study , 2001 .

[13]  Jean-Antoine Désidéri,et al.  Upwind schemes for the two-dimensional shallow water equations with variable depth using unstructured meshes , 1998 .

[14]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[15]  Carlos Parés,et al.  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 , 2001 .

[16]  P. García-Navarro,et al.  On numerical treatment of the source terms in the shallow water equations , 2000 .

[17]  J. G. Lafuente,et al.  Sea Level in the Strait of Gibraltar: Tides , 1990 .

[18]  Bram van Leer,et al.  Progress in Multi-Dimensional Upwind Differencing , 1993 .

[19]  J. Cunge,et al.  Practical aspects of computational river hydraulics , 1980 .

[20]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[21]  D. Farmer,et al.  Maximal two-layer exchange through a contraction with barotropic net flow , 1986, Journal of Fluid Mechanics.

[22]  A. E. Gill.  Atmosphere-Ocean Dynamics , 1982 .

[23]  Philip L. Roe,et al.  Upwind differencing schemes for hyperbolic conservation laws with source terms , 1987 .

[24]  E. Toro.  Godunov Methods: Theory and Applications , 2001 .

[25]  K. Helfrich.  Time-Dependent Two-Layer Hydraulic Exchange Flows , 1995 .

[26]  Alfredo Izquierdo,et al.  TIDES AND TIDAL ENERGETICS OF THE STRAIT OF GIBRALTAR : A MODELLING APPROACH , 1998 .

[27]  David M. Farmer,et al.  Maximal two-layer exchange over a sill and through the combination of a sill and contraction with barotropic flow , 1986, Journal of Fluid Mechanics.