Simulations of Dam-Break Flows Using Free Surface Capturing Method

2008 ◽  
Vol 24 (4) ◽  
pp. 391-403 ◽  
Author(s):  
W.-Y. Chang ◽  
L.-C. Lee ◽  
H.-C. Lien ◽  
J.-S. Lai

AbstractA model adopting the surface capturing method is developed for the simulation of dam-break flows by solving the Navier-Stokes equations of weakly compressible and variable density flows in open channels. Due to the characteristics of weakly compressible flow equations, a compressibility parameter describing the compressibility of fluid is determined to obtain the time-accurate flow fields in both liquid and gas regions simultaneously. Accordingly, the location of free surface can be captured as a discontinuity of the density field for dam-break flow simulations. The numerical algorithm in the proposed method is based on the framework of the finite volume method for discretization in space. To deal with the discontinuity property of fluid density near the free surface, the TVD-MUSCL scheme is adopted to overcome numerical oscillations and dissipation. For discretization in time, the explicit 4-stage Runge-Kutta scheme is employed in the model. Finally, several typical dam-break flow problems in open channel are simulated to demonstrate the validation and applicability of the proposed model.

Author(s):  
B. Elie ◽  
G. Reliquet ◽  
P.-E. Guillerm ◽  
O. Thilleul ◽  
P. Ferrant ◽  
...  

This paper compares numerical and experimental results in the study of the resonance phenomenon which appears between two side-by-side fixed barges for different sea-states. Simulations were performed using SWENSE (Spectral Wave Explicit Navier-Stokes Equations) approach and results are compared with experimental data on two fixed barges with different headings and bilges. Numerical results, obtained using the SWENSE approach, are able to predict both the frequency and the magnitude of the RAO functions.


2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.


1958 ◽  
Vol 8 ◽  
pp. 966-974
Author(s):  
H. E. Petschek

Analyses of aerodynamic dissipation in ordinary un-ionized gases are all based upon the Navier-Stokes equations. These equations relate the rate of dissipation to the local gradients in velocity and temperature through the viscosity and heat conduction coefficients. Although it is true that in many flow situations the magnitude of the total dissipation in the gas does not depend on the magnitude of the viscosity coefficient, this coefficient does determine the minimum scale of variations observed in the gas and the form of the Navier-Stokes equations determines the type of phenomena which are observed on a small scale. In order to discuss dissipation in an ionized gas in the presence of a magnetic field, it is therefore necessary to re-examine the derivation of the basic flow equations. This paper attempts to do this for a case of a completely ionized gas and demonstrates that the basic microscopic dissipation mechanism is appreciably different. For example, it is shown that the minimum length in which the properties of the flow field can change noticeably is appreciably less than one mean free path.


1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.


2008 ◽  
Vol 3 (1) ◽  
Author(s):  
R.K. Saket ◽  
Anil Kumar

This paper presents a convective dominated reliable diffusion process in an axi-symmetric tube with a local constriction simulating a stenos artery considering the porosity effects. The investigations demonstrate the effects of wall shear stress and recirculation flow on the concentration distribution in the vessels lumen and on wall mass transfer keeping the porosity in view. The flow is governed by the incompressible Navier-Stokes equations for Newtonian fluid in porous medium. The convection diffusion equation has been used for the mass transport. The effect of porosity is examined on the velocity field and wall stress. The numerical solutions of the flow equations and the coupled mass transport equations have been obtained using a finite difference method. This paper explains the reliable effects of flow porosity on the mass transport.


Author(s):  
Iraj Saeedpanah ◽  
M. Shayanfar ◽  
E. Jabbari ◽  
Mohammad Haji Mohammadi

Free surface flows are frequently encountered in hydraulic engineering problems including water jets, weirs and around gates. An iterative solution to the incompressible two-dimensional vertical steady Navier-Stokes equations, comprising momentum and continuity equations, is used to solve for the priori unknown free surface, the velocity and the pressure fields. The entire water body is covered by a unstructured finite element grid which is locally refined. The dynamic boundary condition is imposed for the free surface where the pressure vanishes. This procedure is done continuously until the normal velocities components vanish. To overcome numerical errors and oscillations encountering in convection terms, the SUPG (streamline upwinding Petrov-Galerkin) method is applied. The solution method is tested for different discharges onto a standard spillway geometries. The results shows good agreement with available experimental data.


1989 ◽  
Vol 16 (6) ◽  
pp. 829-844
Author(s):  
A. Soulaïmani ◽  
Y. Ouellet ◽  
G. Dhatt ◽  
R. Blanchet

This paper is devoted to the computational analysis of three-dimensional free surface flows. The model solves the Navier-Stokes equations without any a priori restriction on the pressure distribution. The variational formulation along with the solution algorithm are presented. Finally, the model is used to study the hydrodynamic regime in the vicinity of a projected harbor installation. Key words: free surface flows, three-dimensional flows, finite element method.


2012 ◽  
Vol 09 (03) ◽  
pp. 1250036 ◽  
Author(s):  
MOHAMED ABDELWAHED ◽  
MOHAMED AMARA

Due to ever increasing water demand, the preservation of water quality is becoming a very important issue. Eutrophication is amongst the particular problems threatening the quality of water. This paper begins with presenting a mathematical model for aeration process in lake used to combat water eutrophication. Two phases are numerically simulated to study the injected air effect on water by using a corrected one phase model described by Navier–Stokes equations with variable density and viscosity representing the mixture. This model is numerically studied by coupling characteristics scheme for time discretization and mixed finite element method for space approximation. An error estimates in space and time for the velocity are obtained. Numerical results are given firstly in support of the mathematical analysis and secondly to simulate a real application case of the studied problem.


Sign in / Sign up

Export Citation Format

Share Document