Large-wave simulation of microscale breaking waves induced by a free-surface drift layer

Wave Motion ◽  
2007 ◽  
Vol 44 (5) ◽  
pp. 355-370 ◽  
Author(s):  
Athanassios A. Dimas
1999 ◽  
Vol 383 ◽  
pp. 285-305 ◽  
Author(s):  
MATTHEW MILLER ◽  
TOBIAS NENNSTIEL ◽  
JAMES H. DUNCAN ◽  
ATHANASSIOS A. DIMAS ◽  
STEPHAN PRÖSTLER

The effect of free-surface drift layers on the maximum height that a steady wave can attain without breaking is explored through experiments and numerical simulations. In the experiments, the waves are generated by towing a two-dimensional fully submerged hydrofoil at constant depth, speed and angle of attack. The drift layer is generated by towing a plastic sheet on the water surface ahead of the hydrofoil. It is found that the presence of this drift layer (free-surface wake) dramatically reduces the maximum non-breaking wave height and that this wave height correlates well with the surface drift velocity. In the simulations, the inviscid two-dimensional fully nonlinear Euler equations are solved numerically. Initially symmetric wave profiles are superimposed on a parallel drift layer whose mean flow characteristics match those in the experiments. It is found that for large enough initial wave amplitudes a bulge forms at the crest on the forward face of the wave and the vorticity fluctuations just under the surface in this region grow dramatically in time. This behaviour is taken as a criterion to indicate impending wave breaking. The maximum non-breaking wave elevations obtained in this way are in good agreement with the experimental findings.


2012 ◽  
Vol 1 (33) ◽  
pp. 65
Author(s):  
Gerasimos Kolokythas ◽  
Aggelos Dimakopoulos ◽  
Athanassios Dimas

In the present study, the three-dimensional, incompressible, turbulent, free-surface flow, developing by the propagation of nonlinear breaking waves over a rigid bed of constant slope, is numerically simulated. The main objective is to investigate the process of spilling wave breaking and the characteristics of the developing undertow current employing the large-wave simulation (LWS) method. According to LWS methodology, large velocity and free-surface scales are fully resolved, and subgrid scales are treated by an eddy viscosity model, similar to large-eddy simulation (LES) methodology. The simulations are based on the numerical solution of the unsteady, three-dimensional, Navier-Stokes equations subject to the fully-nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow boundary conditions. The case of incoming second-order Stokes waves, normal to the shore, with wavelength to inflow depth ratio λ/dΙ = 6.6, wave steepness H/λ = 0.025, bed slope tanβ = 1/35 and Reynolds number (based on inflow water depth) Red = 250,000 is investigated. The predictions of the LWS model for the incipient wave breaking parameters - breaking depth and height - are in very good agreement with published experimental measurements. Profiles of the time-averaged horizontal velocity in the surf zone are also in good agreement with the corresponding measured ones, verifying the ability of the model to capture adequately the undertow current.


Author(s):  
A. A. Dimas ◽  
A. S. Dimakopoulos ◽  
G. A. Kolokythas

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):  
Aggelos S. Dimakopoulos ◽  
Athanassios A. Dimas

A subgrid scale model is presented for the large-wave simulation (LWS) of spilling breaking waves over bottom of arbitrary shape and finite depth. According to LWS formulation, large velocity and free-surface scales are fully resolved, while subgrid scales are accounted for by an eddy viscosity model, similar to large-eddy simulation (LES). The LWS-based model is applied on the two-dimensional wave propagation over a constant-slope bed. Fluid motion is described by the Euler equations for inviscid but rotational flow, subject to the fully non-linear free-surface boundary conditions. The application of LWS is facilitated by a boundary-fitted transformation, which introduces free-surface elevation terms in the Euler equations and simplifies the numerical implementation. Subgrid velocity scales are modeled similarly to LES, while the effect of free-surface subgrid scales are modeled by wave SGS stresses model. The resulting equations are solved numerically by a two-stage fractional time-step scheme, while an absorption zone is placed in the outflow region to minimize reflection by the outgoing waves. The simulation is carried out for the propagation and breaking of waves over a flat bed with constant slope 1/35 and results are compared to available experimental data. The numerical predictions for the breaking height, the breaking depth and the free-surface elevation dissipation in the surf zone agree very well with the corresponding measurements. The model predicts the vorticity generation in the breaking face of the wave and the appearance of the undertow current in the surf zone. The predicted shear of the undertow current is higher than the measured one.


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 CFD model REEF3D. The model is based on the Reynolds-Averaged Navier-Stokes (RANS) equations coupled with the level set method (LSM) and k-ω turbulence model. Three wave impact conditions are considered in the present 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):  
Jannette B. Frandsen ◽  
Alistair G. L. Borthwick

Nonlinear effects of standing wave motions in fixed and vertically excited tanks are numerically investigated. The present fully nonlinear model analyses two-dimensional waves in stable and unstable regions of the free-surface flow. Numerical solutions of the governing nonlinear potential flow equations are obtained using a finite-difference time-stepping scheme on adaptively mapped grids. A σ-transformation in the vertical direction that stretches directly between the free-surface and bed boundary is applied to map the moving free surface physical domain onto a fixed computational domain. A horizontal linear mapping is also applied, so that the resulting computational domain is rectangular, and consists of unit square cells. The small-amplitude free-surface predictions in the fixed and vertically excited tanks compare well with 2nd order small perturbation theory. For stable steep waves in the vertically excited tank, the free-surface exhibits nonlinear behaviour. Parametric resonance is evident in the instability zones, as the amplitudes grow exponentially, even for small forcing amplitudes. For steep initial amplitudes the predictions differ considerably from the small perturbation theory solution, demonstrating the importance of nonlinear effects. The present numerical model provides a simple way of simulating steep non-breaking waves. It is computationally quick and accurate, and there is no need for free surface smoothing because of the σ-transformation.


2016 ◽  
Author(s):  
Janek Meyer ◽  
Hannes Renzsch ◽  
Kai Graf ◽  
Thomas Slawig

While plain vanilla OpenFOAM has strong capabilities with regards to quite a few typical CFD-tasks, some problems actually require additional bespoke solvers and numerics for efficient computation of high-quality results. One of the fields requiring these additions is the computation of large-scale free-surface flows as found e.g. in naval architecture. This holds especially for the flow around typical modern yacht hulls, often planing, sometimes with surface-piercing appendages. Particular challenges include, but are not limited to, breaking waves, sharpness of interface, numerical ventilation (aka streaking) and a wide range of flow phenomenon scales. A new OF-based application including newly implemented discretization schemes, gradient computation and rigid body motion computation is described. In the following the new code will be validated against published experimental data; the effect on accuracy, computational time and solver stability will be shown by comparison to standard OF-solvers (interFoam / interDyMFoam) and Star CCM+. The code’s capabilities to simulate complex “real-world” flows are shown on a well-known racing yacht design.


Author(s):  
Anne M. Fullerton ◽  
Thomas C. Fu ◽  
Edward S. Ammeen

Impact loads from waves on vessels and coastal structures are highly complex and may involve wave breaking, making these changes difficult to estimate numerically or empirically. Results from previous experiments have shown a wide range of forces and pressures measured from breaking and non-breaking waves, with no clear trend between wave characteristics and the localized forces and pressures that they generate. In 2008, a canonical breaking wave impact data set was obtained at the Naval Surface Warfare Center, Carderock Division, by measuring the distribution of impact pressures of incident non-breaking and breaking waves on one face of a cube. The effects of wave height, wavelength, face orientation, face angle, and submergence depth were investigated. A limited number of runs were made at low forward speeds, ranging from about 0.5 to 2 knots (0.26 to 1.03 m/s). The measurement cube was outfitted with a removable instrumented plate measuring 1 ft2 (0.09 m2), and the wave heights tested ranged from 8–14 inches (20.3 to 35.6 cm). The instrumented plate had 9 slam panels of varying sizes made from polyvinyl chloride (PVC) and 11 pressure gages; this data was collected at 5 kHz to capture the dynamic response of the gages and panels and fully resolve the shapes of the impacts. A Kistler gage was used to measure the total force averaged over the cube face. A bottom mounted acoustic Doppler current profiler (ADCP) was used to obtain measurements of velocity through the water column to provide incoming velocity boundary conditions. A Light Detecting and Ranging (LiDAR) system was also used above the basin to obtain a surface mapping of the free surface over a distance of approximately 15 feet (4.6 m). Additional point measurements of the free surface were made using acoustic distance sensors. Standard and high-speed video cameras were used to capture a qualitative assessment of the impacts. Impact loads on the plate tend to increase with wave height, as well as with plate inclination toward incoming waves. Further trends of the pressures and forces with wave characteristics, cube orientation, draft and face angle are investigated and presented in this paper, and are also compared with previous test results.


Sign in / Sign up

Export Citation Format

Share Document