Kinematics Under Extreme Waves

Author(s):  
Carl Trygve Stansberg ◽  
Ove T. Gudmestad ◽  
Sverre K. Haver

Nonlinear contributions in near-surface particle velocities under extreme crests in random seas can be important in the prediction of wave loads. Four different prediction methods are compared in this paper. The purpose is to observe and evaluate differences in predicted particle velocities under high and extreme crests, and how well they agree with measurements. The study includes linear prediction, a second-order random wave model, Wheeler’s method [1970, “Method for Calculating Forces Produced by Irregular Waves,” JPT, J. Pet. Technol., pp. 359–367] and a new method proposed by Grue et al. [2003, “Kinematics of Extreme Waves in Deep Water,” Appl. Ocean Res., 25, pp. 355–366]. Comparison to laboratory data is also made. The whole wave-zone range from below still water level up to the free surface is considered. Large nonlinear contributions are identified in the near-surface velocities. The results are interpreted to be correlated with the local steepness kA. Some scatter between the different methods is observed in the results. The comparison to experiments shows that among the methods included, the second-order random wave model works best in the whole range under a steep crest in deep or almost deep water, and is therefore recommended. The method of Grue et al. works reasonably well for z>0, i.e., above the calm water level, while it overpredicts the velocities for z<0. Wheeler’s method, when used with a measured or a second-order input elevation record, predicts velocities fairly well at the free surface z=ηmax, but it underpredicts around z=0 and further below. The relative magnitude of this latter error is slightly smaller than the local steepness kA0 and can be quite significant in extreme waves. If Wheeler’s method is used with a linear input, the same error occurs in the whole range, i.e., also at the free surface.

Author(s):  
Carl Trygve Stansberg ◽  
Ove T. Gudmestad ◽  
Sverre K. Haver

Four different methods for prediction of wave-zone particle velocities under steep crests in random seas are compared. The study includes linear prediction, a second-order random wave model, Wheeler’s method, and a new method proposed by Grue et al. (2003). Comparison to laboratory data is also made. The purpose is to observe and evaluate differences in predictions for high and extreme waves, and how well they agree with measurements. The whole range from below still water level up to the free surface is considered. It is found that the second-order random wave model works best at all levels of the water column under a steep crest in deep water, and is therefore recommended. Grue’s method works reasonably well in many cases for z &gt; 0, i.e. above the calm water level, but it overpredicts the velocities for z &lt; 0. Wheeler’s method, when used with a measured or a second-order input elevation record, predicts fairly well the velocities at the free surface z = ηmax, but it underpredicts around z = 0 as well as at lower levels. The relative magnitude of this underprediction is slightly lower than the local steepness kA0 and can be quite significant in extreme waves. If Wheeler’s method is used with a linear input, the same error occurs also at the free surface.


Author(s):  
Jørn Birknes ◽  
Øistein Hagen ◽  
Thomas B. Johannessen ◽  
Øystein Lande ◽  
Arne Nestegård

The present paper is concerned with the prediction of horizontal velocities underneath measured irregular wave surface elevations. The simple case of unidirectional waves in deep water is considered. The main challenge in calculating accurately the kinematics in the crest region is related to the treatment of the contribution from wave components with frequencies much higher than the frequencies near the spectral peak. When using linear or weakly nonlinear perturbation methods, the wave components are superimposed at the still water level and it is necessary to truncate the tail of the spectrum in order to calculate accurately the velocity in the crest region. In the present paper, results from three methods of calculating the crest kinematics are compared with the model test results of Skjelbreia et al. [1]: • The second-order model of Stansberg et al. [5] which truncates consistently the high frequency part of the spectrum. • The second-order model of Johannessen [13] which calculates the velocity directly at the instantaneous free surface. • The Wheeler [3] stretching method which stretches the linear velocity profile from the still water level to the instantaneous free surface. In addition to comparing the horizontal velocity profiles underneath the crest, time traces of horizontal velocity is compared at the free surface in the vicinity of a large crest. The latter comparison highlights the differences between the models and the challenge of accurate predictions close to top of crest. All three models show a reasonable agreement with model test results although it is clear that the first two methods are superior to the Wheeler method.


Author(s):  
Ole A. Hermundstad ◽  
Carl T. Stansberg ◽  
O̸yvind Hellan

A practical method for prediction of slamming loads and structural responses in the bow of an FPSO is presented. Incoming waves are simulated by a second-order random wave model, which describes the water elevation and kinematics. Vessel motions are calculated by linear analysis. The diffracted wave field is calculated taking into account linear 3D diffraction. Relative motions are then estimated by combining the linear vessel motions, second-order incoming waves and linear diffraction. The relative motions and velocities at the bow are used as input to numerical slamming calculations. The bow is divided into 2D sections and a boundary value problem is solved for each section applying the generalized Wagner-method of Zhao & Faltinsen (1993) and Zhao et al (1996). The 2D slamming calculations account for the local pile-up of water on each side of the section during impact. Structural responses are calculated from a finite-element model of the bow using the exact pressure distribution from the slamming calculations. This is achieved by automatic mapping of pressures onto the outer surface of the FE-model and performing a quasi-static structural analysis for each time-step. The methods are implemented into a package of computer tools, allowing the user to perform the various steps in the process with little manual editing of data. The system runs easily on a standard PC. Measurements on a 1:55 scaled model of an FPSO are used for validation of the bow slamming calculations. The model was equipped with five 3.85m × 1.65m (full-scale) panels in the upper part of the bow for slamming force measurements. The tests were run in storm conditions with steep waves. The calculated slamming force on a panel located at the foremost tip of the bulwark, 12.8 meters above the mean waterline, is compared with measured results for selected extreme slamming events. Considering the complexity of this problem and the relative simplicity of the approach, the agreement is very good.


Author(s):  
Xiang Yuan Zheng ◽  
Torgeir Moan ◽  
Ser Tong Quek

The one-dimensional Fast Fourier Transform (FFT) has been extensively applied to efficiently simulate Gaussian wave elevation and water particle kinematics. The actual sea elevation/kinematics exhibit non-Gaussianities that mathematically can be represented by the second-order random wave theory. The elevation/kinematics formulation contains double-summation frequency sum and difference terms which in computation make the dynamic analysis of offshore structural response prohibitive. This study aims at a direct and efficient two-dimensional FFT algorithm for simulating the frequency sum terms. For the frequency difference terms, inverse FFT and FFT are respectively implemented across the two dimensions of the wave interaction matrix. Given specified wave conditions, not only the wave elevation but kinematics and associated Morison force are simulated. Favorable agreements are achieved when the statistics of elevation/kinematics are compared with not only the empirical fits but the analytical solutions developed based on modified eigenvalue/eigenvector approach, while the computation effort is very limited. In addition, the stochastic analyses in both time-and frequency domains show that the near-surface Morison force and induced linear oscillator response exhibits stronger non-Gaussianities by involving the second-order wave effects.


2007 ◽  
Vol 129 (4) ◽  
pp. 327-334 ◽  
Author(s):  
Xiang Yuan Zheng ◽  
Torgeir Moan ◽  
Ser Tong Quek

The one-dimensional fast Fourier transform (FFT) has been applied extensively to simulate Gaussian random wave elevations and water particle kinematics. The actual sea elevations/kinematics exhibit non-Gaussian characteristics that can be represented mathematically by a second-order random wave theory. The elevations/kinematics formulations contain frequency sum and difference terms that usually lead to expensive time-domain dynamic analyses of offshore structural responses. This study aims at a direct and efficient two-dimensional FFT algorithm for simulating the frequency sum terms. For the frequency-difference terms, inverse FFT and forward FFT are implemented, respectively, across the two dimensions of the wave interaction matrix. Given specified wave conditions, the statistics of simulated elevations/kinematics compare well with not only the empirical fits but also the analytical solutions based on a modified eigenvalue/eigenvector approach, while the computational effort of simulation is very economical. In addition, the stochastic analyses in both time domain and frequency domain show that, attributable to the second-order nonlinear wave effects, the near-surface Morison force and induced linear oscillator response are more non-Gaussian than those subjected to Gaussian random waves.


2018 ◽  
Vol 6 (3) ◽  
pp. 105 ◽  
Author(s):  
Ankit Aggarwal ◽  
Csaba Pákozdi ◽  
Hans Bihs ◽  
Dag Myrhaug ◽  
Mayilvahanan Alagan Chella

The experimental wave paddle signal is unknown to the numerical modellers in many cases. This makes it quite challenging to numerically reproduce the time history of free surface elevation for irregular waves. In the present work, a numerical investigation is performed using a computational fluid dynamics (CFD) based model to validate and investigate a non-iterative free surface reconstruction technique for irregular waves. In the current approach, the free surface is reconstructed by spectrally composing the irregular wave train as a summation of the harmonic components coupled with the Dirichlet inlet boundary condition. The verification is performed by comparing the numerically reconstructed free surface elevation with theoretical input waves. The applicability of the present approach to generate irregular waves by reconstructing the free surface is investigated for different coastal and marine engineering problems. A numerical analysis is performed to validate the free surface reconstruction approach to generate breaking irregular waves over a submerged bar. The wave amplitudes, wave frequencies and wave phases are modelled with good accuracy in the time-domain during the higher-order energy transfers and complex processes like wave shoaling, wave breaking and wave decomposition. The present approach to generate irregular waves is also employed to model steep irregular waves in deep water. The free surface reconstruction method is able to simulate the irregular free surface profiles in deep water with low root mean square errors and high correlation coefficients. Furthermore, the irregular wave forces on a monopile are investigated in the time-domain. The amplitudes and phases of the force signal under irregular waves generated by using the current technique are modelled accurately in the time-domain. The proposed approach to numerically reproduce the free surface elevation in the time-domain provides promising and accurate results for all the benchmark cases.


Author(s):  
Dag Myrhaug ◽  
Carl Trygve Stansberg ◽  
Hanne Therese Wist

Statistics of the nonlinear free surface elevation as well as the nonlinear random wave kinematics in terms of the horizontal velocity component in arbitrary water depth are addressed. Two different methods are considered: a simplified analytical approach based on second-order Stokes wave theory including the sum-frequency effect only, and a second-order random wave model including both sum-frequency and difference-frequency effects. The paper compares results for the statistics of the nonlinear free surface, and the consequences of neglecting the difference-frequency effect in the first method are discussed.


Author(s):  
Jenny M. V. Trumars ◽  
Johan O. Jonsson ◽  
Lars Bergdahl

A phase averaging wave model (SWAN) is used to transform offshore sea states to the near to shore site of an offshore wind energy converter. The supporting structure of the wind turbine consists of a cylindrical monopile, and the wave forces and resulting base moments on it are calculated by Morison’s equation integrating from the bottom to the instantaneous free surface. For that purpose the wave-motion in the time domain at the monopile is realized by a second-order random wave model.


2006 ◽  
Vol 25 (5) ◽  
pp. 586-601 ◽  
Author(s):  
M. Onorato ◽  
A.R. Osborne ◽  
M. Serio ◽  
L. Cavaleri ◽  
C. Brandini ◽  
...  

Author(s):  
Carl Trygve Stansberg ◽  
Andreas Amundsen ◽  
Sebastien Fouques ◽  
Ole David Økland

The importance of including second-order nonlinear random wave kinematics in the numerical prediction of drag-induced shear forces and moments, at various levels on a bottom-fixed slender monopile in 40m water depth, is investigated. A vertical circular cylinder of diameter 0.5m is considered, representing typical dimensions of members in jacket type foundations of offshore wind turbines. The focus is here on the wave loads only, and wind and a propeller are therefore not included in this study. In particular, the main focus is on the effects from second-order random wave kinematics on the structural quasi-static time-varying loads due to drag forces in heavy storm wave conditions. Comparisons are made to the traditional use of Airy waves with various ways of stretching. An in-house numerical FEM code developed for structural analysis, NIRWANA, is used for this study. Thus one purpose of the present work is also to verify the implementation of the second-order random waves in the code. The results show significant effects, especially in the wave zone. Extreme crests are around 15%–20% increased, free-surface extreme particle velocities increase by around 30%–40%, while the velocities at levels below MWL are, on the other hand, somewhat reduced. The resulting peak shear forces, and in particular the moments, are thereby increased by typically 50%–100% in the upper parts of the column. At the base the peak shear forces are comparable to the traditional methods, while moments are still somewhat higher. Another effect is the generation of more high-frequency load contributions, which may be important to address further with respect to natural frequencies of such towers.


Sign in / Sign up

Export Citation Format

Share Document