2D numerical simulation of shallow water and bedload transport in channel confluences by considering the non-hydrostatic pressure

2021 ◽  
Author(s):  
Behnam Balouchi ◽  
Nils Rüther ◽  
Mahmood Shafaei Bejestan ◽  
Kordula Valerie Anne Schwarzwälder ◽  
Hans Bihs
Keyword(s):  
Hydrostatic Pressure ◽  
Shallow Water ◽  
Flow Pattern ◽  
Stokes Equations ◽  
Maximum Velocity ◽  
Bedload Transport ◽  
Water Model ◽  
Effective Parameters ◽  
Complex Flow ◽  
Channel Confluence

<p>Channel confluence is one of the important sections of channel networks which is also common encountered in nature. Six different zones exist at a channel confluence: 1) stagnation zone, 2) flow deflection zone, 3) flow separation zone, 4) maximum velocity zone, 5) flow recovery zone and 6) shear layers between combining flows zone. Due to the complexity of flow pattern at channel confluence, this location is always interesting among researchers. Although there are a number of studies on the flow and sediment pattern at confluences, there are still some gaps to be studied. Hence, a calibrated numerical model should be a good tool for evaluating the various effective parameters on flow and sediment patterns. The numerical 2D shallow-water model used in this paper is SFLOW which was developed by NTNU. Besides, the model calibration part of the current study is done by using a set of data from laboratory experiments.</p><p>This study attempt to simulate bed changes at channel confluences with a 2D shallow-water modeling under non-hydrostatic pressure, and show the applicability of the SFLOW model for this complex flow pattern. SFLOW solving the depth-averaged Navier-Stokes equations which is equipped with cutting-edge solvers. Besides, SFLOW modeled turbulency with depth-averaged two-equation RANS. In comparison with other codes, one of the interesting features of SFLOW is solving the non-hydrostatic pressure besides of hydrostatic part. This leads to a more realistic representation of the complex flow and sediment patterns of channel confluences, and consider less computational power than full 3D models.</p>

scholarly journals Propagation of Solitary Waves over a Submerged Slotted Barrier

2020 ◽  
Vol 8 (6) ◽  
pp. 419 ◽  
Author(s):  
Yun-Ta Wu ◽  
Shih-Chun Hsiao
Keyword(s):  
Flow Pattern ◽  
Solitary Waves ◽  
Vortex Shedding ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Model Data ◽  
Free Surface Elevation ◽  
Complex Flow ◽  
Laboratory Observation ◽  
Navier Stokes Equations

In this article, the interaction of solitary waves and a submerged slotted barrier is investigated in which the slotted barrier consists of three impermeable elements and its porosity can be determined by the distance between the two neighboring elements. A new experiment is conducted to measure free surface elevation, velocity, and turbulent kinetic energy. Numerical simulation is performed using a two-dimensional model based on the Reynolds-Averaged Navier-Stokes equations and the non-linear k-ɛ turbulence model. A detailed flow pattern is illustrated by a flow visualization technique. A laboratory observation indicates that flow separations occur at each element of the slotted barrier and the vortex shedding process is then triggered due to the complicated interaction of those induced vortices that further create a complex flow pattern. During the vortex shedding process, seeding particles that are initially accumulated near the seafloor are suspended by an upward jet formed by vortices interacting. Model-data comparisons are carried out to examine the accuracy of the model. Overall model-data comparisons are in satisfactory agreement, but modeled results sometimes fail to predict the positions of the induced vortices. Since the measured data is unique in terms of velocity and turbulence, the dataset can be used for further improvement of numerical modeling.

scholarly journals Congested shallow water model: roof modeling in free surface flow

2018 ◽  
Vol 52 (5) ◽  
pp. 1679-1707 ◽  
Author(s):  
Edwige Godlewski ◽  
Martin Parisot ◽  
Jacques Sainte-Marie ◽  
Fabien Wahl
Keyword(s):  
Shallow Water ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Hyperbolic Equations ◽  
Mechanical Energy ◽  
Free Surface Flow ◽  
Surface Flow ◽  
Water Model ◽  
Shallow Water Model ◽  
Time Step

We are interested in the modeling and the numerical approximation of flows in the presence of a roof, for example flows in sewers or under an ice floe. A shallow water model with a supplementary congestion constraint describing the roof is derived from the Navier-Stokes equations. The congestion constraint is a challenging problem for the numerical resolution of hyperbolic equations. To overcome this difficulty, we follow a pseudo-compressibility relaxation approach. Eventually, a numerical scheme based on a finite volume method is proposed. The well-balanced property and the dissipation of the mechanical energy, acting as a mathematical entropy, are ensured under a non-restrictive condition on the time step in spite of the large celerity of the potential waves in the congested areas. Simulations in one dimension for transcritical steady flow are carried out and numerical solutions are compared to several analytical (stationary and non-stationary) solutions for validation.

On the necessity of non-hydrostatic pressure models for free surface flow over complex topography

2020 ◽  
Author(s):  
Isabel Echeverribar ◽  
Pilar Brufau ◽  
Pilar García-Navarro
Keyword(s):  
Free Surface ◽  
Hydrostatic Pressure ◽  
Shallow Water ◽  
Shallow Water Equations ◽  
Water Model ◽  
Finite Volume Scheme ◽  
Complex Topography ◽  
Source Terms ◽  
Pressure Correction ◽  
Wide Range

<p><span><strong>There is a wide range of geophysical flows, such as flow in open channels and rivers, tsunami and flood modeling, that can be mathematically represented by the non-linear shallow water 1D equations involving hydrostatic pressure assumptions as an approximation of the Navier Stokes equations. In this context, special attention must be paid to bottom source terms integration and numerical corrections when dealing with wet/dry fronts or strong slopes in order to obtain physically-based solutions (Murillo and García-Navarro, 2010) in complex and realistic cases with irregular topography. However, although these numerical corrections have been developed in recent years achieving not only more robust models but also more accurate results, they still might find a limit when dealing with specific scenarios where vertical information or disspersive effects become crucial. This work presents a 1D shallow water model that introduces vertical information by means of a non-hydrostatic pressure correction when necessary. In particular, the pressure correction method (Hirsch, 2007) is applied to a 1D finite volume scheme for a rectification of the velocity field in free surface scenarios. It is solved by means of an implicit scheme, whereas the depth-integrated shallow water equations are solved using an explicit scheme. It is worth highlighting that it preserves all the advantages and numerical fixes aforementioned for the pure shallow water system. Computations with and without non-hydrostatic corrections are compared for the same cases to test the validity of the conventional hydrostatic pressure assumption at some scenarios involving complex topography.</strong></span></p><p><span>[1] J. Murillo and P. Garcia-Navarro, Weak solutions for partial differential equations with source terms: application to the shallow water equations, Journal of Computational Physics, vol. 229, iss. 11, pp. 4327-4368, 2010.</span></p><p><span>[2] C. Hirsch, Numerical Computation of Internal and External flows: The fundamentals of Computational Fluid Dynamics, Butterworth-Heinemann, 2007.</span></p>

scholarly journals Towards a new friction model for shallow water equations through an interactive viscous layer

2019 ◽  
Vol 53 (1) ◽  
pp. 269-299 ◽  
Author(s):  
François James ◽  
Pierre-Yves Lagrée ◽  
Minh H. Le ◽  
Mathilde Legrand
Keyword(s):  
Shallow Water ◽  
Stokes Equations ◽  
Water Model ◽  
Finite Volume Scheme ◽  
Inviscid Fluid ◽  
Phase Lag ◽  
Viscous Layer ◽  
Viscous Effects ◽  
Shallow Water Models ◽  
Water Models

The derivation of shallow water models from Navier–Stokes equations is revisited yielding a class of two-layer shallow water models. An improved velocity profile is proposed, based on the superposition of an inviscid fluid and a viscous layer inspired by the Interactive Boundary Layer interaction used in aeronautics. This leads to a new friction law which depends not only on velocity and depth but also on the variations of velocity and thickness of the viscous layer. The resulting system is an extended shallow water model consisting of three depth-integrated equations: the first two are mass and momentum conservation in which a slight correction on hydrostatic pressure has been made; the third one, known as von Kármán equation, describes the evolution of the viscous layer. This coupled model is shown to be conditionally hyperbolic, and a Godunov-type finite volume scheme is also proposed. Several numerical examples are provided and compared to the Multi-Layer Saint-Venant model. They emphasize the ability of the model to deal with unsteady viscous effects. They illustrate also the phase-lag between friction and topography, and even recover possible reverse flows.

Modeling of Tarbela Reservoir and Water Flow Simulation Through Its Spillways

2010 ◽  
Author(s):  
Muhammad Abid ◽  
Muftooh Ur Rehman Siddiqi
Keyword(s):  
Stokes Equations ◽  
Numerical Study ◽  
Discharge Rate ◽  
Flow Simulation ◽  
Maximum Velocity ◽  
Warming Trend ◽  
Height Distribution ◽  
Complex Flow ◽  
Local Flow ◽  
Buoyant Force

This numerical study is performed to predict the flow patterns and characteristics in Tarbela dam which is a multipurpose dam during the summer season the flow of the dam reaches to its design capacity or near flooding due to weather changes resulting from the global warming trend. A 3D model was made in Pro-Engineer® and meshed in ICEM CFD®. Commercially known software, ANSYS CFX®, was applied to numerically solve the Navier-Stokes equations for solution domain. The calculated results such as pressure, velocities, flow rate, surface height, and water buoyant force were compared with the actual data where available. The numerical calculations show uneven discharge through each gate due to the complex flow pattern just upstream of the weir. Maximum velocity was observed along the spillways outlet. In conclusion, the results from numerical simulation are generally well agreed with the existing data, the flow information such as flow field patterns at increased flow, local flow disturbances, discharge rate and surface height distribution obtained used for the behavior of existing dam and can be used for engineering design purpose of future dams.

Axisymmetric gravity currents in a rotating system: experimental and numerical investigations

2001 ◽  
Vol 447 ◽  
pp. 1-29 ◽  
Author(s):  
MARK A. HALLWORTH ◽  
HERBERT E. HUPPERT ◽  
MARIUS UNGARISH
Keyword(s):  
Flow Field ◽  
Shallow Water ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Salt Water ◽  
Gravity Currents ◽  
The Body ◽  
Water Model ◽  
Maximum Radius ◽  
Theoretical Predictions

The propagation at high Reynolds number of a heavy, axisymmetric gravity current of given initial volume over a horizontal boundary is considered in both rotating and non-rotating situations. The investigation combines experiments with theoretical predictions by both shallow-water approximations and numerical solutions of the full axisymmetric equations. Attention is focused on cases when the initial ratio of Coriolis to inertia forces is small. The experiments were performed by quickly releasing a known cylindrical volume of dense salt water of 2 m diameter at the centre of a circular tank of diameter 13 m containing fresh ambient water of typical depth 80 cm. The propagation of the current was recorded for different initial values of the salt concentration, the volume of released fluid, the ratio of the initial height of the current to the ambient depth, and the rate of rotation. A major feature of the rotating currents was the attainment of a maximum radius of propagation. Thereafter a contraction–relaxation motion of the body of fluid and a regular series of outwardly propagating pulses was observed. The frequency of these pulses is slightly higher than inertial, and the amplitude is of the order of magnitude of half the maximum radius. Theoretical predictions of the corresponding gravity currents were also obtained by (i) previously developed shallow-water approximations (Ungarish & Huppert 1998) and (ii) a specially developed finite-difference code based on the full axisymmetric Navier–Stokes equations. The ‘numerical experiments’ provided by this code are needed to capture details of the flow field (such as the non-smooth shape of the interface, the vertical dependence of the velocity field) which are not reproduced by the shallow-water model and are very difficult for, or outside the range of, accurate experimental measurement. The comparisons and discussion provide insight into the flow field and indicate the advantages and limitations of the verified simulation tools.

scholarly journals ASYMPTOTIC DERIVATION OF THE SECTION-AVERAGED SHALLOW WATER EQUATIONS FOR NATURAL RIVER HYDRAULICS

2009 ◽  
Vol 19 (03) ◽  
pp. 387-417 ◽  
Author(s):  
ASTRID DECOENE ◽  
LUCA BONAVENTURA ◽  
EDIE MIGLIO ◽  
FAUSTO SALERI
Keyword(s):  
Shallow Water ◽  
Open Channel ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Classical Equation ◽  
Free Surface Flows ◽  
Water Model ◽  
Local Flow ◽  
Longitudinal Length ◽  
Friction Term

The section-averaged shallow water model usually applied in river and open channel hydraulics is derived by an asymptotic analysis that accounts for terms up to second order in the vertical/longitudinal length ratio, starting from the three-dimensional Reynolds-averaged Navier–Stokes equations for incompressible free surface flows. The derivation is carried out under quite general assumptions on the geometry of the channel, thus allowing for the application of the resulting equations to natural rivers with arbitrarily shaped cross sections. As a result of the derivation, a generalized friction term is obtained, that does not rely on local uniformity assumptions and that can be computed directly from three-dimensional turbulence models, without need for local uniformity assumptions. The modified equations including the novel friction term are compared to the classical Saint Venant equations in the case of steady state open channel flows, where analytic solutions are available, showing that the solutions resulting from the modified equation set are much closer to the three-dimensional solutions than those of the classical equation set. Furthermore, it is shown that the proposed formulation yields results that are very similar to those obtained with empirical friction closures widely applied in computational hydraulics. The generalized friction term derived therefore justifies a posteriori these empirical closures, while allowing to avoid the assumptions on local flow uniformity on which these closures rely.

scholarly journals Simulation of Ocean Circulation of Dongsha Water Using Non-Hydrostatic Shallow-Water Model

Water ◽  
2020 ◽  
Vol 12 (10) ◽  
pp. 2832 ◽  
Author(s):  
Shin-Jye Liang ◽  
Chih-Chieh Young ◽  
Chi Dai ◽  
Nan-Jing Wu ◽  
Tai-Wen Hsu
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Wave Propagation ◽  
Free Surface ◽  
Hydrostatic Pressure ◽  
Velocity Field ◽  
Shallow Water ◽  
Shallow Water Equations ◽  
Water Model ◽  
Shallow Water Model

A two-dimensional non-hydrostatic shallow-water model for weakly dispersive waves is developed using the least-squares finite-element method. The model is based on the depth-averaged, nonlinear and non-hydrostatic shallow-water equations. The non-hydrostatic shallow-water equations are solved with the semi-implicit (predictor-corrector) method and least-squares finite-element method. In the predictor step, hydrostatic pressure at the previous step is used as an initial guess and an intermediate velocity field is calculated. In the corrector step, a Poisson equation for the non-hydrostatic pressure is solved and the final velocity and free-surface elevation is corrected for the new time step. The non-hydrostatic shallow-water model is verified and applied to both wave and flow driven fluid flows, including solitary wave propagation in a channel, progressive sinusoidal waves propagation over a submerged bar, von Karmann vortex street, and ocean circulations of Dongsha Atolls. It is found hydrostatic shallow-water model is efficient and accurate for shallow water flows. Non-hydrostatic shallow-water model requires 1.5 to 3.0 more cpu time than hydrostatic shallow-water model for the same simulation. Model simulations reveal that non-hydrostatic pressure gradients could affect the velocity field and free-surface significantly in case where nonlinearity and dispersion are important during the course of wave propagation.

Sign in / Sign up

Export Citation Format

Share Document