Three-Dimensional Numerical Modeling of Pier Scour Under Current and Waves Using Level-Set Method

Author(s):  
Mohammad Saud Afzal ◽  
Hans Bihs ◽  
Arun Kamath ◽  
Øivind A. Arntsen

A three-dimensional (3D) computational fluid dynamics (CFD) model is used to calculate the scour and the deposition pattern around a pier for two different boundary conditions: constant discharge and regular waves. The CFD model solves Reynolds-Averaged Navier–Stokes (RANS) equations in all three dimensions. The location of the free-surface is represented using the level-set method (LSM), which calculates the complex motion of the free-surface in a very realistic manner. For the implementation of waves, the CFD code is used as a numerical wave tank. For the geometric representation of the moveable sediment bed, the LSM is used. The numerical results for the local scour prediction are compared with physical experiments. The decoupling of the hydrodynamic and the morphodynamic time step is tested and found to be a reasonable assumption. For the two situations of local pier scour under current and wave conditions, the numerical model predicts the general evolution (geometry, location, and maximum scour depth) and time development of the scour hole accurately.

Author(s):  
Mohammad Saud Afzal ◽  
Hans Bihs ◽  
Arun Kamath ◽  
Øivind A. Arntsen

Stability of offshore structures can be threatened by local scouring which could ultimately lead to their failure. As a consequence, knowledge of the scouring mechanism and the accurate prediction of the characteristic scour geometry is very important for the design of such structures. A three-dimensional computational fluid dynamics model is used to calculate the scour and the deposition pattern around a pier for two different boundary conditions: constant discharge and regular waves. The computational fluid dynamics (CFD) model solves Reynolds-Averaged Navier-Stokes (RANS) equations in all three dimensions. The location of the free surface is represented using the level set method, which calculates the complex motion of the free surface in a very realistic manner. For the implementation of waves, the CFD code is used as a numerical wave tank. For the geometric representation of the moveable sediment bed, the level set method is used. The numerical results for the local scour prediction are compared with physical experiments. The performance of the turbulence models, the formulations of the critical shear stress for the sloping bed and the effect of the variation of the sediment time stepping are investigated. The decoupling of the hydrodynamic and the morphodynamic time step is tested and found to be a reasonable assumption. For the two situations of local pier scour under current and wave conditions, the numerical model predicts the general evolution (geometry, location and maximum scour depth) and time development of the scour hole accurately.


2013 ◽  
Author(s):  
Gabriel Reliquet ◽  
Aurélien Drouet ◽  
Pierre-Emmanuel Guillerm ◽  
Erwan Jacquin ◽  
Lionel Gentaz ◽  
...  

The purpose of this paper is to present combination of the SWENSE (Spectral Wave Explicit Navier-Stokes Equations – [1]) method — an original method to treat fully nonlinear wave-body interactions — and a free surface RANSE (Reynolds Averaged Navier-Stokes Equations) solver using a single-phase Level Set method to capture the interface. The idea is to be able to simulate wave-body interactions under viscous flow theory with strong deformations of the interface (wave breaking in the vicinity of the body, green water on ship decks…), while keeping the advantages of the SWENSE scheme. The SWENSE approach is based on a physical decomposition by combining incident waves described by a nonlinear spectral scheme based on potential flow theory and an adapted Navier-Stokes solver where only the diffracted part of the flow is solved, incident flow parameters seen as forcing terms. In the single-phase Level Set method [2, 3], the air phase is neglected. Thus, only the liquid phase is solved considering a fluid with uniform properties. The location of the free surface is determined by a Level Set function initialised as the signed distance. The accuracy of simulation depends essentially on the pressure scheme used to impose free surface dynamic boundary condition. Comparisons of numerical results with experimental and numerical data for US navy combatant DTMB 5415 in calm water and in head waves are presented.


Author(s):  
Erlend Liavåg Grotle ◽  
Hans Bihs ◽  
Vilmar Æsøy ◽  
Eilif Pedersen

In this paper, numerical simulations of nonlinear sloshing in rectangular tanks are presented. Model implementations in the open source software reef3d are tested, and the results are compared with experimental data from three different conditions. The interface location is compared for both linear and nonlinear sloshing. The nonlinear sloshing is simulated in both two-dimensional (2D) and three-dimensional (3D). Video images from the SPHERIC project are compared with simulations for the interface. A condition with lateral wave impacts in sloshing, with a frequency close to the natural frequency of the first mode, can be found in this case. The numerical model is solving the Reynolds-averaged Navier–Stokes (RANS) equations with the k–ω turbulence model. The level set method is used to capture the interface. Higher order discretization schemes are implemented to handle time-evolution and convective fluxes. A ghost cell method is used to account for solid boundaries and parallel computations. It is found that the limiting factor for the eddy-viscosity has significant influence in the nonlinear sloshing cases. As the sloshing becomes more violent, the increased strain at the gas–liquid interface overproduces turbulence energy with unrealistically high damping of the motion. Three-dimensional simulations show slightly better comparison than 2D. Due to nonlinearities and small damping, the time to reach steady-state may take several cycles. The last case shows promising results for the global motion. As expected, the breakup of the liquid surface makes it difficult to resolve each phase. But overall, the numerical model predicts the sloshing motion reasonably well.


2009 ◽  
Vol 80 (12) ◽  
pp. 1520-1543 ◽  
Author(s):  
Qinglin Duan ◽  
Jeong-Hoon Song ◽  
Thomas Menouillard ◽  
Ted Belytschko

2008 ◽  
Vol 11 (4-6) ◽  
pp. 221-235 ◽  
Author(s):  
S. P. van der Pijl ◽  
A. Segal ◽  
C. Vuik ◽  
P. Wesseling

2009 ◽  
Vol 635 ◽  
pp. 1-26 ◽  
Author(s):  
NICOLAS PÉRINET ◽  
DAMIR JURIC ◽  
LAURETTE S. TUCKERMAN

We simulate numerically the full dynamics of Faraday waves in three dimensions for two incompressible and immiscible viscous fluids. The Navier–Stokes equations are solved using a finite-difference projection method coupled with a front-tracking method for the interface between the two fluids. The critical accelerations and wavenumbers, as well as the temporal behaviour at onset are compared with the results of the linear Floquet analysis of Kumar & Tuckerman (J. Fluid Mech., vol. 279, 1994, p. 49). The finite-amplitude results are compared with the experiments of Kityk et al (Phys. Rev. E, vol. 72, 2005, p. 036209). In particular, we reproduce the detailed spatio-temporal spectrum of both square and hexagonal patterns within experimental uncertainty. We present the first calculations of a three-dimensional velocity field arising from the Faraday instability for a hexagonal pattern as it varies over its oscillation period.


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.


Sign in / Sign up

Export Citation Format

Share Document