On the Validity and Sensitivity of CFD Simulations for a Deterministic Breaking Wave Impact on a Semi Submersible

Author(s):  
Henry Bandringa ◽  
Joop A. Helder

To assess the integrity and safety of structures offshore, prediction of run-up, green water, and impact loads needs to be made during the structure’s design. For predicting these highly non-linear phenomena, most of the offshore industry relies on detailed model testing. In the last couple of years however, CFD simulations have shown more and more promising results in predicting these events, see for instance [1]–[4]. To obtain confidence in the accuracy of CFD simulations in the challenging field of extreme wave impacts, a proper validation of such CFD tools is essential. In this paper two CFD tools are considered for the simulation of a deterministic breaking wave impact on a fixed semi submersible, resulting in flow phenomena like wave run-up, horizontal wave impact and deck impacts. Hereby, one of the CFD tools applies an unstructured gridding approach and implicit free-surface reconstruction, and uses an implicit time integration with a fixed time step. The other CFD tool explicitly reconstructs the free surface on a structured grid and integrates the free surface explicitly in time, using a variable time step. The presented simulations use a compact computational domain with wave absorbing boundary conditions and local grid refinement to reduce CPU time. Besides a typical verification and validation of the results, for one of the CFD tools a sensitivity study is performed in which the influence of small variations in the incoming breaking wave on the overall results is assessed. Such an analysis should provide the industry more insight in the to-be-expected sensitivity (and hence uncertainty) of CFD simulations for these type of applications. Experiments carried out by MARIN are used to validate all the presented simulation results.

2002 ◽  
Vol 456 ◽  
pp. 295-318 ◽  
Author(s):  
YING LI ◽  
FREDRIC RAICHLEN

The run-up of non-breaking and breaking solitary waves on a uniform plane beach connected to a constant-depth wave tank was investigated experimentally and numerically. If only the general characteristics of the run-up process and the maximum run-up are of interest, for the case of a breaking wave the post-breaking condition can be simplified and represented as a propagating bore. A numerical model using this bore structure to treat the process of wave breaking and subsequent shoreward propagation was developed. The nonlinear shallow water equations (NLSW) were solved using the weighted essentially non-oscillatory (WENO) shock capturing scheme employed in gas dynamics. Wave breaking and post-breaking propagation are handled automatically by this scheme and ad hoc terms are not required. A computational domain mapping technique was used to model the shoreline movement. This numerical scheme was found to provide a relatively simple and reasonably good prediction of various aspects of the run-up process. The energy dissipation associated with wave breaking of solitary wave run-up (excluding the effects of bottom friction) was also estimated using the results from the numerical model.


Author(s):  
Gerasimos A. Kolokythas ◽  
Athanassios A. Dimas

In the present study, numerical simulations of the free-surface flow, developing by the propagation of nonlinear water waves over a rippled bottom, are performed assuming that the corresponding flow is two-dimensional, incompressible and viscous. The simulations are based on the numerical solution of the Navier-Stokes equations subject to the fully-nonlinear free-surface boundary conditions and the suitable bottom, inflow and outflow boundary conditions. The equations are properly transformed so that the computational domain becomes time-independent. For the spatial discretization, a hybrid scheme with finite-differences and Chebyshev polynomials is applied, while a fractional time-step scheme is used for the temporal discretization. A wave absorption zone is placed at the outflow region in order to efficiently minimize reflection of waves by the outflow boundary. The numerical model is validated by comparison to the analytical solution for the laminar, oscillatory, current flow which develops a uniform boundary layer over a horizontal bottom. For the propagation of finite-amplitude waves over a rigid rippled bed, the case with wavelength to water depth ratio λ/d0 = 6 and wave height to wavelength ratio H0/λ = 0.05 is considered. The ripples have parabolic shape, while their dimensions — length and height — are chosen accordingly to fit laboratory and field data. Results indicate that the wall shear stress over the ripples and the form drag forces on the ripples increase with increasing ripple height, while the corresponding friction force is insensitive to this increase. Therefore, the percentage of friction in the total drag force decreases with increasing ripple height.


Author(s):  
Bülent Düz ◽  
René H. M. Huijsmans ◽  
Peter R. Wellens ◽  
Mart J. A. Borsboom ◽  
Arthur E. P. Veldman ◽  
...  

Numerical simulations of wave phenomena necessarily have to be carried out in a limited computational domain. This implies that incoming waves should be prescribed properly, and the outgoing waves should leave the domain without causing reflections. In this paper we will present an enhanced type of such generating and absorbing boundary conditions (GABC). The new approach is applied in studies of extreme hydrodynamic wave impact on rigid and floating structures in offshore and coastal engineering, for which the VOF-based CFD simulation tool ComFLOW has been developed.


Author(s):  
Dakui Feng ◽  
Xianzhou Wang ◽  
Zhiguo Zhang ◽  
Yanming Guan

The catamaran is composed of two monohulls, the flow fields between the inner and outer side of each monohull are different, the bodies must be considered as lifting bodies. So it is very important to know the lifting effect on hydrodynamic characteristics of catamaran hull at the preliminary design stage of its hull form. The pressure Kutta condition is imposed on the trailing-surface of the lifting body by determining the dipole distribution, which generates required circulation on the lifting part. The method is based on Green’s second theorem. Rankine Sources and dipoles are placed on boundary surfaces. Time-stepping scheme is adopted to simulate the wave generated by the catamaran with a uniform speed in deep water. The values of the potential and position of the free surface are updated by integrating the nonlinear Lagrangian free surface boundary conditions for every time. A moving computational window is used in the computations by truncating the fluid domain (the free surface) into a computational domain. The grid regeneration scheme is developed to determine the approximate position of the free surface for the next time step. An implicit implement of far field condition is enforced automatically at the truncation boundary of the computational window, Radiation condition is satisfied automatically. The influences on the wave making resistance of the distance between the twin hulls of the Wigley catamaran on the hydrodynamic characteristics are discussed. The numerical results are presented compared with the existing simulation result. The method can be used to simulate the flow fields around the foil near free surface.


2006 ◽  
Vol 128 (2) ◽  
pp. 89-99 ◽  
Author(s):  
J. Zang ◽  
R. Gibson ◽  
P. H. Taylor ◽  
R. Eatock Taylor ◽  
C. Swan

The objective of this research, part of the EU FP5 REBASDO Program, is to examine the effects of second order wave diffraction in wave run-up around the bow of a vessel (FPSO) in random seas. In this work, the nonlinear wave scattering problem is solved by employing a quadratic boundary element method. A computer program, DIFFRACT, has been developed and recently extended to deal with unidirectional and directional bichromatic input wave systems, calculating second order wave diffraction loads and free surface elevation under regular waves and focused wave groups. The second order wave interaction with a vessel in a unidirectional focused wave group is presented in this paper. Comparison of numerical results and experimental measurements conducted at Imperial College shows excellent agreement. The second order free surface components at the bow of the ship are very significant, and cannot be neglected if one requires accurate prediction of the wave-structure interaction; otherwise a major underestimation of the wave impact on the structure could occur.


Author(s):  
Csaba Pakozdi ◽  
Timothy E. Kendon ◽  
Carl-Trygve Stansberg

The slamming of breaking waves on the legs of large volume offshore platforms has received increased attention over recent years. To investigate this problem, MARINTEK’s Wave Impact Loads JIP has, in one of its sub-tasks, focused towards an idealised model test setup of a rectangular cylinder in breaking waves. The model consists of a vertical column with a fragment of a horizontal deck attached. The model is fixed at a distance L ahead of the wave maker. Physical scale model test experiments of the block in regular waves and in wave groups have been carried out in Phase 1 of the JIP (2008). The objective of this study is the CFD simulation of a long crested breaking wave and its impact on the aforementioned cylinder and deck structure in order to find out the feasibility of the numerical reconstruction of such events. The commercial CFD tool Star-CCM+ V5.03.0056 (www.cd-adapco.com) is used in this study. This paper considers results from the test setup, and compares the measured wave elevation against results from the CFD code. The position of the cylinder in relation to the breaking wave front is investigated in the numerical simulation in order to analyze its effect on the slamming force. Use of an unsteady wave boundary condition, matching the exact motion history of the wave-maker with the measured free surface elevation at the wave maker gives an almost exact matching between the computed wave profile and the measured wave profile. The improvement in the numerical tool of Star-CCM+ which makes it possible to use higher order time integration scheme for VOF significantly decreases the numerical diffusion of the wave propagation. This new scheme also enables the use of a time step 10 times larger than the first order scheme which reduces the computational time. Because a large time step can be chosen it is important that the time step is small enough to capture the correct time evolution of the physical phenomena of interest. Capturing the pressure evolution at a slamming event demands very high spatial resolution. Spatially averaged slamming pressures look fairly similar to the model test observations, while further work is needed for a more detailed comparison.


Author(s):  
Mayilvahanan Alagan Chella ◽  
Hans Bihs ◽  
Dag Myrhaug ◽  
Øivind Asgeir Arntsen

Wave loads from breaking waves on offshore wind turbine (OWT) substructures in shallow waters still remain uncertain. The interaction of breaking waves with structures is characterized by complex free surface deformations, instantaneous impact of the water mass against the structure, and consequently large wave forces on the structures. The main objective of the paper is to investigate wave impact pressures and kinematics during the interaction of breaking waves with a vertical cylinder using the open-source computational fluid dynamics (CFD) model REEF3D. The model is based on the Reynolds-averaged Navier–Stokes (RANS) equations coupled with the level set method and k–ω turbulence model. Three wave impact conditions are considered in this study. The numerically simulated free surface deformations around the cylinder during the breaking wave interaction are also presented for different wave impact conditions. For three wave impact conditions, the wave impact pressure and the horizontal and vertical components of the particle velocity are computed in front of the cylinder and analyzed. The pressure and velocity profile at their maximum values are also examined and discussed. In addition, the total force is calculated for three breaking conditions and they are correlated with the pressure and kinematics during the interaction.


Author(s):  
Ruosi Zha ◽  
Wei Qiu ◽  
Heather Peng

It is important to predict wave impact and structural responses of ship and offshore structures in extreme sea conditions. There are advantages to apply Lagrangian particle methods to simulate highly nonlinear breaking free surface flow and fluid-structure interactions (FSI). In this paper, an improved moving particle semi-implicit (MPS) method was developed to solve the FSI problems. At each time step, the fluid motions and the structural responses are solved. For flow computations, a modified mixed source term method and an improved free surface identification method were adopted to suppress pressure oscillations. Moreover, a particle collision model was used to enhance the numerical stability and avoid nonphysical solutions. The discretized Poisson equations for pressure were solved by a parallel version of the bi-conjugate gradient stabilized method based on the message passing interface (MPI) approach. For structural responses, solids were treated as isotropic elastic particles. Validation studies were carried out for cases of 2D dam breaking and its interaction with a rubber gate. The numerical solutions are in good agreement with experimental data and other published numerical results.


Geophysics ◽  
1991 ◽  
Vol 56 (2) ◽  
pp. 231-241 ◽  
Author(s):  
R. L. Higdon

Absorbing boundary conditions are needed for computing numerical models of wave motions in unbounded spatial domains. The boundary conditions developed here for elastic waves are generalizations of ones developed earlier for acoustic waves. These conditions are based on compositions of simple first‐order differential operators. The formulas can be applied without modification to problems in both two and three dimensions. The boundary conditions are stable for all values of the ratio of P‐wave velocity to S‐wave velocity, and they are effective near a free surface and in a horizontally stratified medium. The boundary conditions are approximated with simple finite‐difference equations that use values of the solution only along grid lines perpendicular to the boundary. This property facilitates implementation, especially near a free surface and at other corners of the computational domain.


Author(s):  
Erik Jan de Ridder ◽  
Tim Bunnik ◽  
Johan M. Peeringa ◽  
Bo Terp Paulsen ◽  
Christof Wehmeyer ◽  
...  

The objective of the Joint Industry Project Wave impact on Fixed foundations (WiFi JIP) was to increase the understanding of breaking and steep wave impact’s on fixed foundations of offshore wind turbines (OWT). The project was set-up as a Joint Industry Project (JIP) and in total 20 companies and research institutes participated in the project. In this paper a summary of the complete WiFi JIP project will be presented. At the start of the project the state of the art design methods and guidelines were reviewed (WP1). Thereafter a jacket and a monopile foundation were designed using these state-of-the-art tools that were available at the start of the project. This effort has been reported in WP2 , where design computations were carried out using the embedded stream function approach for several sea states. In this WP Siemens, ECN and Ramboll also calculated the impact response of the monopile to surging and spilling type wave breakers with their engineering tools. In the next phase the designed foundations were tested in MARIN’s shallow water model basin. The foundation for the monopile was modelled as a rigid and flexible foundation. The foundations were tested in regular waves, irregular sea states and so called focused waves. During the model tests the wave heights, wave run-up, accelerations, impact pressures and loads on the foundation and boat landing were measured. The model test results were reported in WP3 and 7 and used as validation for WP9 and 10. WP4 delivered more understanding of realistic design conditions for areas typical for OWT, like the North Sea. Particular attention was paid to the probability of occurrence of breaking and steep waves and the associated slamming load. For this an extensive 5 week experimental program was performed from September to October 2013 in the wide wave-current flume at Deltares (Atlantic Basin). During these tests both waves and current were simulated and two bathymetries. WP8 provided analyses of the performed full scale measurements on the response of a OWT. The full-scale measurements were done for a Vestas V90 3MW wind turbine in the Belwind windfarm which is located 46 km off the coast of Zeebrugge on the Bligh Bank. The CFD simulations performed in WP 9 showed that a good agreement is obtained between the CFD simulations and the model and full scale measurement. In work package 10, an improved methodology was developed based on the outcome of the previous WP’s to model the breaking wave impact of plunging type breakers. In WP11 and 12 this new approach is applied on different case study’s by ECN.


Sign in / Sign up

Export Citation Format

Share Document