Computational Fluid Dynamics Simulations of Regular and Irregular Waves Past a Horizontal Semi-Submerged Cylinder

Author(s):  
Shengnan Liu ◽  
Muk Chen Ong ◽  
Charlotte Obhrai ◽  
Sopheak Seng

Two-dimensional (2D) numerical simulations have been performed to investigate both regular and irregular waves past a fixed horizontally semisubmerged circular cylinder. The 2D simulations are carried out by solving Navier–Stokes equations discretized by finite volume method. Volume of fluid (VOF) method is employed to capture the free surface in the numerical wave tank (NWT). Validation studies have been performed by comparing the numerical results of free surface waves past the cylinder with the published experimental and numerical data. The present numerical results are in good agreement with both the experimental and the other numerical results in terms of hydrodynamic forces and free surface elevation. Subsequently, the effects of the wave height and the wavelength on wave–structure interaction are investigated by conducting numerical simulations on the regular and the irregular waves past a semisubmerged cylinder at different wave heights and the wavelengths. The averaged and maximum vertical wave forces on the cylinder increase with the increasing wave height. The numerical results for the irregular waves are compared with those induced by the regular waves in terms of the maximum and averaged vertical wave forces. When the significant wave height and the spectral peak period of the irregular waves are equal to the wave height and the wave period of the regular waves, the maximum vertical wave force induced by the irregular waves is larger than that induced by the regular waves, meanwhile, the average vertical wave forces have the contrary relationship.

Author(s):  
Shengnan Liu ◽  
Muk Chen Ong ◽  
Charlotte Obhrai ◽  
Sopheak Seng

Two-dimensional (2D) numerical simulations have been performed using OpenFOAM (an open source CFD software package [1]) and waves2Foam (an OpenFOAM based add-on library for wave generations and absorption [2]) to investigate free surface waves past one fixed horizontally semi-submerged cylinder. The 2-D simulations are carried out by solving Navier-Stokes equations which are discretized based on finite volume method (FVM). Volume of Fluid (VOF) method is employed to capture the free surface in the numerical wave tank. Validation studies have been performed by comparing the numerical results of Stokes first-order wave past a semi-submerged circular cylinder with the published experimental data at different incident wave properties. The numerical results are in good agreement with the experimental data. Subsequently, regular and irregular waves past semi-submerged cylinder at different wave heights and the wave lengths are computed numerically to investigate the effect of the wave height and wave length on wave-structure interaction. The numerical results for irregular waves are compared with those induced by regular waves.


1986 ◽  
Vol 1 (20) ◽  
pp. 131
Author(s):  
Chien-Kee Chang ◽  
Ching-Her Hwang

Wave pressure is the most important external force for the design of breakwater. During recent years, there has been considerable development in the technology of vertical face breakwater; however, there is no reliable method to compute wave forces induced by irregular waves. The purpose of this study is to obtain statistical characteristics of irregular wave pressure distribution from the data of model tests. The results of this study shown that vertical face breakwater under the action of irregular waves, some waves are reflected, so that the next wave breaks a critical distance resulting in a rapidly rising shock pressure on the breakwater. On the average, the wave pressure increase with incoming wave height, but the maximum wave force does not necessarily occur for the largest wave height. It can be occurred for serval larger wave group in an appropiate phase composition. The irregular wave pressure distribution on the breakwater is quite uniform; the ratio of tested and calculated wave pressures decreases with the reduction of relative crest height of breakwater. Coda formula can predict the total horizontal force of the upper part of breakwater quite well except exetreme shock pressure occurred by non-breaking waves. Wave forces calculated by Miche-Rundgren and Nagai wave force formula are about 10% cummulated exceeding percentage of wave force obtained from model test.


Author(s):  
Ankit Aggarwal ◽  
Mayilvahanan Alagan Chella ◽  
Arun Kamath ◽  
Hans Bihs ◽  
Øivind Asgeir Arnsten

In the present study, the irregular wave forces on a fully submerged circular cylinder are investigated using the open-source computational fluid dynamics (CFD) model REEF3D. A complete three dimensional representation of the ocean waves requires the consideration of the sea surface as an irregular wave train with the random characteristics. The numerical model uses the incompressible Reynolds-averaged Navier-Stokes (RANS) equations together with the continuity equation to solve the fluid flow problem. Turbulence modeling is carried out using the two equation k-ω model. Spatial discretization is done using an uniform Cartesian grid. The level set method is used for computing the free surface. For time discretization, third-order total variation diminishing (TVD) Runge Kutta scheme is used. Ghost cell boundary method is used for implementing the complex geometries in the numerical model. MPI is used for the exchange of the value of a ghost cell. Relaxation method is used for the wave generation. The numerical model is validated for the irregular waves for a wave tank without any structure. Further, the numerical model is validated by comparing the numerical results with the experimental data for a fully submerged circular cylinder under regular waves and irregular waves. The numerical results are in a good agreement with the experimental data for the regular and irregular wave forces. The JONSWAP spectrum is used for the wave generation. The free surface features and kinematics around the cylinder is also presented and discussed.


1992 ◽  
Vol 25 (9) ◽  
pp. 211-216
Author(s):  
A. Akyarli ◽  
Y. Arisoy

As the wave forces are the function of the wave height, period and the angle between the incoming wave direction and the axis of the discharge pipeline, the resultant wave force is directly related to the alignment of the pipeline. In this paper, a method is explained to determine an optimum pipeline route for which the resultant wave force becomes minimum and hence, the cost of the constructive measures may decrease. Also, the application of this method is submitted through a case study.


2021 ◽  
Vol 9 (5) ◽  
pp. 520
Author(s):  
Zhenyu Liu ◽  
Zhen Guo ◽  
Yuzhe Dou ◽  
Fanyu Zeng

Most offshore wind turbines are installed in shallow water and exposed to breaking waves. Previous numerical studies focusing on breaking wave forces generally ignored the seabed permeability. In this paper, a numerical model based on Volume-Averaged Reynolds Averaged Navier–Stokes equations (VARANS) is employed to reveal the process of a solitary wave interacting with a rigid pile over a permeable slope. Through applying the Forchheimer saturated drag equation, effects of seabed permeability on fluid motions are simulated. The reliability of the present model is verified by comparisons between experimentally obtained data and the numerical results. Further, 190 cases are simulated and the effects of different parameters on breaking wave forces on the pile are studied systematically. Results indicate that over a permeable seabed, the maximum breaking wave forces can occur not only when waves break just before the pile, but also when a “secondary wave wall” slams against the pile, after wave breaking. With the initial wave height increasing, breaking wave forces will increase, but the growth can decrease as the slope angle and permeability increase. For inclined piles around the wave breaking point, the maximum breaking wave force usually occurs with an inclination angle of α = −22.5° or 0°.


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 ◽  
Author(s):  
Chien Ming Wang ◽  
Huu Phu Nguyen ◽  
Jeong Cheol Park ◽  
Mengmeng Han ◽  
Nagi abdussamie ◽  
...  

<p>Floating breakwaters have been used to protect shorelines, marinas, very large floating structures, dockyards, fish farms, harbours and ports from harsh wave environments. A floating breakwater outperforms its bottom-founded counterpart with respect to its environmental friendliness, cost-effectiveness in relatively deep waters or soft seabed conditions, flexibility for expansion and downsizing and its mobility to be towed away. The effectiveness of a floating breakwater design is assessed by its wave attenuation performance that is measured by the wave transmission coefficient (i.e., the ratio of the transmitted wave height to the incident wave height or the ratio of the transmitted wave energy to the incident wave energy). In some current design guidelines for floating breakwaters, the transmission coefficient is estimated based on the assumption that the realistic ocean waves may be represented by regular waves that are characterized by the significant wave period and wave height of the wave spectrum. There is no doubt that the use of regular waves is simple for practicing engineers designing floating breakwaters. However, the validity and accuracy of using regular waves in the evaluation of wave attenuation performance of floating breakwaters have not been thoroughly discussed in the open literature. This study examines the wave transmission coefficients of floating breakwaters by performing hydrodynamic analysis of some large floating breakwaters in ocean waves modelled as regular waves as well as irregular waves described by a wave spectrum such as the Bretschneider spectrum. The formulation of the governing fluid motion and boundary conditions are based on classical linear hydrodynamic theory. The floating breakwater is assumed to take the shape of a long rectangular box modelled by the Mindlin thick plate theory. The finite element – boundary element method was employed to solve the fluid-structure interaction problem. By considering heave-only floating box-type breakwaters of 200m and 500m in length, it is found that the transmission coefficients obtained by using the regular wave model may be smaller (or larger) than that obtained by using the irregular wave model by up to 55% (or 40%). These significant differences in the transmission coefficient estimated by using regular and irregular waves indicate that simplifying assumption of realistic ocean waves as regular waves leads to significant over/underprediction of wave attenuation performance of floating breakwaters. Thus, when designing floating breakwaters, the ocean waves have to be treated as irregular waves modelled by a wave spectrum that best describes the wave condition at the site. This conclusion is expected to motivate a revision of design guidelines for floating breakwaters for better prediction of wave attenuation performance. Also, it is expected to affect how one carries out experiments on floating breakwaters in a wave basin to measure the wave transmission coefficients.</p>


Author(s):  
Malene H. Vested ◽  
Stefan Carstensen ◽  
Erik Damgaard Christensen

As the demand for offshore wind energy continues to grow, the strive to understand the wave forces acting on the substructure of the wind turbines continues. In regard to wind turbine design, it is vital to consider not only the total wave force, but also the local wave forces. Local forces are particularly important for the design of secondary structures as e.g. mooring platforms. Typically, however, experimental studies mainly concern total forces or idealized local forces. We present here a rather simple way to measure local forces along a model monopile. The study is conducted in a wave flume of 28 m in length, in which waves are generated by a piston-type wave maker at a water depth of 0.515 m and shoal onto a bed of slope 1:25. A model monopile is installed and subjected to forcing from a series of both regular and irregular waves. In the experimental set-up, the model monopile is fixed at the bottom and the top and consists of seven independent cylindrical sections. The cylindrical sections are connected by force transducers which measure local shear, and so the associated local forces may be determined. The measured local forces are compared to the force distribution given by Morisons equation combined with linear theory and Wheeler stretching, which is a force estimate commonly used in the industry. This study shows that the total force is rather well captured by Morison’s equation. The force distribution estimated from Morison’s equation, however, shows larger discrepancies from the measured forces. This encourages for further measurements. In this study, we show that it is possible to measure force distribution on a model monopile in a simple and cost-effective manner. The aim is here to demonstrate the method and we will later present a larger body of work associated with the outcome of the measurements.


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.


2016 ◽  
Author(s):  
Mayilvahanan Alagan Chella ◽  
Hans Bihs ◽  
Arun Kamath ◽  
Dag Myrhaug ◽  
Øivind Asgeir Arnsten

The main purpose of the study is to investigate the breaking wave interaction with a group of four circular cylinders. The physical process of wave breaking involves many parameters and an accurate numerical modelling of breaking waves and the interaction with a structure remain a challenge. In the present study, the open-source (Computational Fluid Dynamics) CFD model REEF3D is used to simulate the breaking wave interaction with the multiple cylinders. The numerical model is based on the incompressible Reynolds Averaged Navier-Stokes (RANS) equations, the level set method for the free surface and the k–ω model for turbulence. The model uses a 5th-order conservative finite difference WENO scheme for the convective discretization and a 3rd-order Runge-Kutta scheme for time discretization. The numerical model is validated with experimental data of large-scale experiments for the free surface elevation and the breaking wave force on a single cylinder. A good agreement is seen between the numerical results and experimental data. Two different configurations with four cylinders are examined: in-line square configuration and diamond square configuration. The breaking wave forces on each cylinder in the group are computed for the two cases and the results are compared with the breaking wave force on a single isolated cylinder. Further, the study investigates the water surface elevations and the free surface flow features around the cylinders. In general, the cylinders in both configurations experience the maximum forces lower than the maximum force on a single cylinder. The results of the present study show that the interference effects from the neighbouring cylinders in a group strongly influence the kinematics around and the breaking wave forces on them.


Sign in / Sign up

Export Citation Format

Share Document