Hans Bihs

Numerical Modeling of Abutment Scour with the Focus on the Incipient Motion on Sloping Beds

A three-dimensional computational fluid dynamics model is applied to predict local scour around an abutment in a rectangular laboratory flume. When modeling local scour, steep bed slopes up to the angle of repose occur. To predict the depth and the shape of the local scour correctly, the reduction of the critical shear stress due to the sloping bed must be taken into account. The focus of this study is to investigate different formulas for the threshold of noncohesive sediment motion on sloping beds. Some formulas only take the transversal angle (perpendicular to the flow direction) into account, but others also consider the longitudinal angle (streamwise direction). The numerical model solves the transient Reynolds-averaged Navier-Stokes equations in all three dimensions to compute the water flow. Sediment continuity in combination with an empirical formula is used to capture the bed load transport and the resulting bed changes. When the sloping bed exceeds the angle of repose, the bed slope is corrected with a sand-slide algorithm. The results from the numerical simulations are compared with data from physical experiments. The reduction of the bed shear stress on the sloping bed improves the results of the numerical simulation distinctly. The best results are obtained with the formulas that use both the transversal and the longitudinal angle for the reduction of the critical bed shear stress.

Evaluating wave forces on groups of three and nine cylinders using a 3D numerical wave tank

The evaluation of the complex wave regime due to wave interaction with a large group of cylinders placed in proximity requires an efficient and accurate numerical model. This paper presents the application of a two-phase Computational Fluid Dynamics (CFD) model to carry out a detailed investigation of wave forces and flow around vertical circular cylinders placed in groups of different configurations at low Keulegan-Carpenter (KC) numbers. The 3D numerical wave tank is validated by comparing the numerical results with experimental data. Further, the hydrodynamic effects associated with three cylinders placed in tandem, side by side and in a 3 × 3 square array of nine cylinders are investigated. Wave forces are seen to reduce along the row in a tandem array. In a side-by-side arrangement, the central cylinder experiences the highest force. A combination of these effects is seen in the 3 × 3 square array. The variation of the wave forces on the cylinders in the array for different center-to-center distances and incident wavelengths is evaluated and the results show that the wave forces are the highest on the cylinders when the center-to-center distance is slightly less than half the incident wavelength.

Energy transfer due to shoaling and decomposition of breaking and non-breaking waves over a submerged bar

ABSTRACT Wave propagation over a submerged bar is simulated using the open source CFD model REEF3D with various incident wave heights to study shoaling, wave breaking features and the process of wave decomposition into higher harmonics for relatively long waves of kd=0.52. The computed free surface elevations are compared with experimental data and good agreement is obtained for both non-breaking and spilling breaking waves for both the wave phase and free surface elevation, which has been difficult to obtain in current literature. The differences in the mode of wave shoaling over the weather side slope and the wave decomposition over the lee side slope of the submerged bar are discussed. The evolution of spilling breakers and plunging breakers over the bar crest is also studied. It is found that the free surface elevation continuously increases due to shoaling in the case of non-breaking waves, whereas breaking waves propagate with much lower free surface elevations after breaking over the bar crest. The power spectra of the free surface elevations at various locations indicate that the wave energy in the fundamental frequency is reduced by 76 for the non-breaking wave with kA=0.015 and by about 90 in other cases with higher incident wave heights with kA=0.023−0.034 due to energy dissipation and energy transfer to higher harmonic components as the wave propagates over the submerged bar.

Numerical Simulation of Flows Past Partially-Submerged Horizontal Circular Cylinders in Free Surface Waves

Two-dimensional (2D) numerical simulations are performed to investigate the flows past partially-submerged circular cylinders in free surface waves. The 2D simulations are carried out by solving the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with the k-ω turbulence model. The level set method is employed to model the free-surface waves. Validation studies of a numerical wave tank have been performed by comparing the numerical results with the analytical results obtained from the linear-wave theory. Wave forces on the partially-submerged cylinders have been calculated numerically and compared with the published theoretical and experimental data under regular-wave conditions. The free-surface elevations around the cylinders have been investigated and discussed.Copyright

Calculation of Wave Forces on Cylindrical Piles Using a 3D Numerical Wave Tank

Offshore constructions generally include a large number of vertical cylinders in the support structure. The calculation of wave forces on a vertical cylinder and hydrodynamic effects on it in the presence of neighbouring cylinders is of practical interest. In this paper, a 3D numerical model is used to calculate wave forces on bottom fixed cylindrical piles. Two cases are considered in this study: a single cylinder and a pair of tandem cylinders. A scenario with multiple cylindrical structures in close proximity introduces complex wave-structure interactions and would be of great interest to observe this in detail in a three-dimensional simulation. The wave force exerted on a cylindrical pile is numerically calculated by integrating the pressure and the wall shear stress around the surface of the cylinder. In the case of the single cylinder, the force calculated by the model is compared to the force predicted by the Morison formula and MacCamy-Fuchs theory. In the second case, the pair of cylinders is aligned in the direction of the incoming waves. The numerically calculated inline wave force on each cylinder is compared to the analytical solution for this setup and a good agreement is seen.The Reynolds-Averaged Navier-Stokes equations are used as the governing equations for the fluid flow in the numerical model. The convective terms are discretized using a 5th-order conservative finite difference WENO scheme. A 3rd-order accurate TVD Range-Kutta scheme is used for time discretization. Chorin’s projection method is used to discretize the pressure. The Poisson equation for pressure is solved using a preconditioned BiCGStab algorithm. The level set method is used to obtain a sharp representation of the free water surface. Turbulence in the flow is simulated using the k-ω model. The numerical model is adapted to parallel processing using the MPI library to improve the computing performance of the code.Copyright

Numerical Modeling of Breaking Waves Over a Reef With a Level-Set Based Numerical Wave Tank

Wave breaking is a highly unsteady, non-linear and extremely turbulent phenomenon. During the wave breaking process, the energy of the wave system is focused close to the crest of the wave and a spatial spread of wave energy occurs. Thus, the description of such a physical phenomenon is highly complex and it requires a deep insight into the breaking wave process. The accurate assessment of breaking wave kinematics is essential for an accurate prediction of hydrodynamic loads on structures. Besides, the understanding of the transformation of waves propagating over an artificial or natural reef is important concerning the coastal processes.The numerical model used in this study is a two-phase model, which solves the flow problem for air and water simultaneously. The Navier-Stokes equations are solved on uniform Cartesian grids in two dimensions. The complex free surface is captured by the level set method. A staggered grid is used for the computation with the velocities defined at the cell edges and the pressure at the cell centres. This avoids unphysical pressure oscillations that can occur due to the coupling of pressure and velocity in the incompressible Navier-Stokes equations. The Ghost Cell Immersed Boundary Method is employed to handle the boundary conditions for complex boundaries. Turbulence modelling is carried out using the k-ω model. Discretization of the convective terms is performed using the 5th order Weighted Essentially Non-Oscillatory (WENO) scheme. In this study, a two-dimensional numerical wave tank is used to simulate waves propagating over steep slopes and wave dissipation. The main objective of the present study is to investigate the wave breaking process over a submerged reef. This is accomplished by examining the wave profile during wave breaking and the breaker indices. Also, the numerical results are compared to data from physical experiments and the numerical results exhibit reasonable agreement with experimental data.Copyright

Numerical Modeling of Berm Breakwater Optimization With Varying Berm Geometry Using REEF3D

Harbours are important infrastructures for an offshore production chain. These harbours are protected from the actions of sea by breakwaters to ensure safe loading, unloading of vessels and also to protect the infrastructure. In current literature, research regarding the design of these structures is majorly based on physical model tests. In this study a new tool, a threedimensional numerical model is introduced. The open-source CFD (Computational Fluid Dynamics) model REEF3D is used to study the design of berm breakwaters. The model uses the Volume averaged Reynolds Averaged Navier-Stokes (VRANS) equations to solve the porous flows. At first the VRANS approach in REEF3D is validated for flow through porous media. A dam break case is simulated and comparisons are made for the free surface both inside and outside the porous medium. The numerical model REEF3D is applied to show how to extend the database obtained with purely numerical results, simulating different structural alternatives for the berm in a berm breakwater. Different simulations are conducted with varying berm geometry. The influence of the berm geometry on the pore pressure and velocities are studied. The resulting optimal berm geometry is compared to the geometry according to empirical formulations.

Extreme Wave Generation, Breaking and Impact Simulations With REEF3D

An accurate description of extreme waves is necessary in order to estimate maximum wave forces on offshore structures. On several occasions freak waves have been observed in the past, some causing severe damage. In order to model such extreme wave conditions with a computational fluid dynamics (CFD) model, emphasize needs to be put on the wave generation. One possibility is to use focused waves of first or second order based on irregular sea state wave spectra. For focused waves, the wave phase is chosen, so that the waves focus in a predetermined location at a specified time. Numerical tests have shown, that generating extreme waves based on this method is somewhat limited. The individual wave components are steep enough, that they start to break before the focus location. In the current paper, transient wave packets are used for extreme wave generation. This way, extreme waves can be generated that are higher, but only break at the concentration point. The transient wave packets method is implemented in the open-source CFD software REEF3D. This model uses the level set method for interface capturing. For the hydrodynamics, the Navier-Stokes equations are solved in three dimensions. The code employs a staggered Cartesian mesh, ensuring tight pressure-velocity coupling. Complex geometries are handled with a ghost cell immersed boundary method. High-performance computing is enabled through domain decomposition based parallelization. Convection discretization of the different flow variables is performed with the fifth-order WENO (weighted essentially non-oscillatory) scheme. For the explicit time treatment a third-order Runge-Kutta scheme is selected. In order to validate the extreme wave generation, numerical tests in an empty wave tank are performed and compared with experimental data. Then, the extreme wave breaking on a vertical circular cylinder is investigated. ∗Corresponding author, hans.bihs@ntnu.no

Numerical study of irregular breaking wave forces on a monopile for offshore wind turbines

Abstract The substructures of offshore wind turbines are subjected to different types of hydrodynamic loads from sea states. The wave forces exerted by irregular breaking waves are one of the serious concerns because of the uncertainties involved in defining the breaking wave and the resulting force calculations. In the present study, irregular breaking wave forces on a vertical pile structure are investigated using an open-source Computational Fluid Dynamics (CFD) model REEF3D. The Level Set Method (LSM) is used for modelling the free surface. The Bretschneider spectrum is used for the irregular wave generation. This is validated in the numerical wave tank by comparing the numerical wave spectrum with the experimental wave spectrum. The wave free surface is calculated at three wave gauge locations and compared with experiments. It is observed that the peak of spectral wave density is higher for the wave gauge located besides the cylinder due to shoaling, wave run up and reflections from the cylinder and the peak of spectral wave density is lower for the wave gauge located behind the cylinder due to wave breaking. Further, simulations are performed to study the wave forces on a monopile due to the depth-limited breaking waves. A good match is observed with the experimental and numerical results. Numerical wave energy spectra at different locations along the tank are compared to study the changes in the wave surface elevations due to the interaction of irregular breaking waves with a monopile. The statistical parameters for free surface elevation and wave forces are further investigated. The free surface features around the monopile during its interaction with waves are also studied.

3D Numerical modelling of pile scour with free surface profile under waves and current using the level set method in model REEF3D

Wave action stirs up sediments and keeps it in suspension while currents wash it away from the coastal zone. The combined action of the waves and current close to the sediment bed may worsen the situation by creating excessive sediment transport leading to the failure of hydraulic structures. In this study, numerical modeling of local scour under waves and current is carried out using the open source CFD model REEF3D, which solves the Navier-Stokes equation using the finite difference method. The simulated flow field from the Navier-Stokes equations is coupled with sediment transport algorithms in a numerical waves tank. Further, the calculated bedload and suspended load are linked with the Exner formula to calculate bed elevation changes. The free surface and scoured bed surface are captured using the level set method. Two case scenarios, namely scour under waves and scour under current are run until the equilibrium scour condition is achieved. The simulated results are compared with experimental data of Link (2006) and Sumer & Fredsoe (2001). Good comparison between experimental data and simulated results is observed. It is observed that for equal flow velocity in the flume, sediment transport under current only condition is larger than under waves alone.