scholarly journals The Momentum Conserving Scheme for Two-Layer Shallow Flows

Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 346
Author(s):  
Putu Veri Swastika ◽  
Sri Redjeki Pudjaprasetya

This paper confronts the numerical simulation of steady flows of fluid layers through channels of varying bed and width. The fluid consists of two immiscible fluid layers with constant density, and it is assumed to be of a one-dimensional shallow flow. The governing equation is a coupled system of two-layer shallow water models. In this paper, we apply a direct extension of the momentum conserving scheme previously used for solving the one layer shallow water equations. Computations of various steady-state solutions are used to demonstrate the performance of the proposed numerical scheme. Under the influence of a given flow rate, the numerical steady interface is generated in a channel topography with a hump. The results obtained confirm the analytic steady interface of the two-layer rigid-lid model. Furthermore, the same scheme was used with an additional artificial damping to simulate the maximal exchange flow in channels of varying width. The numerical steady interface agreed well with the analytical steady solutions.

2012 ◽  
Vol 695 ◽  
pp. 341-365 ◽  
Author(s):  
Philip L.-F. Liu ◽  
Xiaoming Wang

AbstractIn this paper, a multi-layer model is developed for the purpose of studying nonlinear internal wave propagation in shallow water. The methodology employed in constructing the multi-layer model is similar to that used in deriving Boussinesq-type equations for surface gravity waves. It can also be viewed as an extension of the two-layer model developed by Choi & Camassa. The multi-layer model approximates the continuous density stratification by an $N$-layer fluid system in which a constant density is assumed in each layer. This allows the model to investigate higher-mode internal waves. Furthermore, the model is capable of simulating large-amplitude internal waves up to the breaking point. However, the model is limited by the assumption that the total water depth is shallow in comparison with the wavelength of interest. Furthermore, the vertical vorticity must vanish, while the horizontal vorticity components are weak. Numerical examples for strongly nonlinear waves are compared with laboratory data and other numerical studies in a two-layer fluid system. Good agreement is observed. The generation and propagation of mode-1 and mode-2 internal waves and their interactions with bottom topography are also investigated.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Thida Pongsanguansin ◽  
Montri Maleewong ◽  
Khamron Mekchay

A well-balanced scheme with total variation diminishing Runge-Kutta discontinuous Galerkin (TVD-RK DG) method for solving shallow water equations is presented. Generally, the flux function at cell interface in the TVD-RK DG scheme is approximated by using the Harten-Lax-van Leer (HLL) method. Here, we apply the weighted average flux (WAF) which is higher order approximation instead of using the HLL in the TVD-RK DG method. The consistency property is shown. The modified well-balanced technique for flux gradient and source terms under the WAF approximations is developed. The accuracy of numerical solutions is demonstrated by simulating dam-break flows with the flat bottom. The steady solutions with shock can be captured correctly without spurious oscillations near the shock front. This presents the other flux approximations in the TVD-RK DG method for shallow water simulations.


2011 ◽  
Vol 27 (3) ◽  
pp. 415-430 ◽  
Author(s):  
W.-D. Guo ◽  
J.-S. Lai ◽  
G.-F. Lin ◽  
F.-Z. Lee ◽  
Y.-C. Tan

ABSTRACTThis paper adopts the finite-volume multi-stage (FMUSTA) scheme to the two-dimensional coupled system combining the shallow water equations and the advection-diffusion equation. For the convection part, the numerical flux is estimated by adopting the FMUSTA scheme, where high order accuracy is achieved by the data reconstruction using the monotonic upstream schemes for conservation laws method. For the diffusion part, the evaluations of first-order derivatives are solved via the method of Jacobian transformation. The hydrostatic reconstruction method is employed for treatment of source terms. The overall accuracy of resulting scheme is second-order both in time and space. In addition, the scheme is non-oscillatory and conserves the pollutant mass during the transport process. For scheme validation, six advection and diffusion transport tests are simulated. The influences of the grid spacing and limiters on the numerical performance are also discussed. Furthermore, the scheme is employed in the simulation of suspended sediment transport in natural-irregular river topography. From the satisfactory agreements between the simulated results and the field measured data, it is demonstrated that the proposed FMUSTA scheme is practically suitable for hydraulic engineering applications.


2009 ◽  
Vol 66 (6) ◽  
pp. 1735-1748 ◽  
Author(s):  
W. T. M. Verkley

Abstract A global version of the equivalent barotropic vorticity equation is derived for the one-layer shallow-water equations on a sphere. The equation has the same form as the corresponding beta plane version, but with one important difference: the stretching (Cressman) term in the expression of the potential vorticity retains its full dependence on f 2, where f is the Coriolis parameter. As a check of the resulting system, the dynamics of linear Rossby waves are considered. It is shown that these waves are rather accurate approximations of the westward-propagating waves of the second class of the original shallow-water equations. It is also concluded that for Rossby waves with short meridional wavelengths the factor f 2 in the stretching term can be replaced by the constant value f02, where f0 is the Coriolis parameter at ±45° latitude.


2020 ◽  
pp. 2150138
Author(s):  
Hajar F. Ismael ◽  
Aly Seadawy ◽  
Hasan Bulut

In this paper, we consider the shallow water wave model in the (2+1)-dimensions. The Hirota simple method is applied to construct the new dynamics one-, two-, three-, [Formula: see text]-soliton solutions, complex multi-soliton, fusion, and breather solutions. By using the quadratic function, the one-lump, mixed kink-lump and periodic lump solutions to the model are obtained. The Hirota bilinear form variable of this model is derived at first via logarithmic variable transform. The physical phenomena to this model are explored. The obtained results verify the proposed model.


1999 ◽  
Vol 396 ◽  
pp. 223-256 ◽  
Author(s):  
B. S. BROOK ◽  
S. A. E. G. FALLE ◽  
T. J. PEDLEY

Unsteady flow in collapsible tubes has been widely studied for a number of different physiological applications; the principal motivation for the work of this paper is the study of blood flow in the jugular vein of an upright, long-necked subject (a giraffe). The one-dimensional equations governing gravity- or pressure-driven flow in collapsible tubes have been solved in the past using finite-difference (MacCormack) methods. Such schemes, however, produce numerical artifacts near discontinuities such as elastic jumps. This paper describes a numerical scheme developed to solve the one-dimensional equations using a more accurate upwind finite volume (Godunov) scheme that has been used successfully in gas dynamics and shallow water wave problems. The adapatation of the Godunov method to the present application is non-trivial due to the highly nonlinear nature of the pressure–area relation for collapsible tubes.The code is tested by comparing both unsteady and converged solutions with analytical solutions where available. Further tests include comparison with solutions obtained from MacCormack methods which illustrate the accuracy of the present method.Finally the possibility of roll waves occurring in collapsible tubes is also considered, both as a test case for the scheme and as an interesting phenomenon in its own right, arising out of the similarity of the collapsible tube equations to those governing shallow water flow.


2013 ◽  
Vol 15 (4) ◽  
pp. 1234-1257 ◽  
Author(s):  
Carmelo Juez ◽  
Javier Murillo ◽  
Pilar García-Navarro

Two-dimensional (2D) transient flow over an erodible bed can be modelled using shallow-water equations and the Exner equation to describe the morphological evolution of the bed. Considering the fact that well-proven capacity formulae are based on one-dimensional (1D) experimental steady flows, the assessment of these empirical relations under unsteady 1D and 2D situations is important. In order to ensure the reliability of the numerical experimentation, the formulation has to be general enough to allow the use of different empirical laws. Moreover, the numerical scheme must handle correctly the coupling between the 2D shallow-water equations and the Exner equation under any condition. In this work, a finite-volume numerical scheme that includes these two main features will be exploited here in 1D and 2D laboratory test cases. The relative performances of Meyer-Peter and Müller, Ashida and Michiue, Engelund and Fredsoe, Fernandez Luque and Van Beek, Parker, Smart, Nielsen, Wong and Camenen and Larson formulations are analysed in terms of the root mean square error. A new discretization of the Smart formula is provided, leading to promising predictions of the erosion/deposition rates. The results arising from this work are useful to justify the use of an empirical sediment bed-load discharge formula among the ones studied, regardless of the hydrodynamic situation.


Sign in / Sign up

Export Citation Format

Share Document