ITTC Modelling of Wave Recommanded
ITTC Modelling of Wave Recommanded
5-02
-07-01.2
Procedures and Guidelines Page 1 of 27
Guideline
Disclaimer
All the information in ITTC Recommended Procedures and Guidelines is published in good faith. Neither ITTC
nor committee members provide any warranties about the completeness, reliability, accuracy or otherwise of this
information. Given the technical evolution, the ITTC Recommended Procedures and Guidelines are checked reg-
ularly by the relevant committee and updated when necessary. It is therefore important to always use the latest
version.
Any action you take upon the information you find in the ITTC Recommended Procedures and Guidelines is
strictly at your own responsibility. Neither ITTC nor committee members shall be liable for any losses and/or
damages whatsoever in connection with the use of information available in the ITTC Recommended Procedures
and Guidelines.
Table of Contents
1. PURPOSE OF GUIDELINE
2.1.2 Non-linear effects – analysis, control
The purpose of this recommended guideline
is to help laboratories in conducting tests on Real water waves are not exactly linear.
waves. The guideline addresses, among other They are influenced by nonlinear effects which
things, nonlinear effects, extreme events and are of higher order with respect to the wave am-
wave breaking, wave generation and control, plitude. Higher-order solutions may be derived
and interaction with wind and currents. by using the Stokes’s expansion as shown in
Newman (1977) or Dean and Dalrymple (2000),
or they may be derived by using fully nonlinear
2. WAVE MODELLING
methods, as e.g. presented in Rienecker and
Fenton, (1981). In a non-dimensional form, the
2.1 Regular waves wave amplitude is better represented by the
wave steepness ε = kA. With increasing wave
Ideally, regular waves are periodic unidirec-
steepness ε, the wave profile progressively de-
tional progressive wave trains, with a single
viates from a pure sinusoidal wave form and are
(monochromatic) basic harmonic.
characterised by higher crests and shallower
troughs, Toffoli et al. (2005). The asymmetry
2.1.1 Linear waves with respect to the horizontal axis (in this case a
vertical asymmetry) is generally referred to as
For most regular wave applications, the av- skewness, which is defined as (e.g. Babanin
erage wave height H and the average period T (2011))
are of main interest. Amplitudes A, defined by
H/2 or by crests AC and troughs AT, and the av- 𝐴𝐴𝐶𝐶
𝑆𝑆𝑘𝑘 = −1
erage steepness kA, are also used, where k is the 𝐴𝐴 𝑇𝑇
angular wave number.
and is positive, unless some other phenomena,
Theoretically, properties should be constant such as wave breaking, occurs. In the above
throughout time and in space, but in physical equation AC and AT denote the crest height and
generation in a wave basin, there is always a cer- the trough depth, respectively.
tain level of variation. Time windows for analy-
sis of the generated wave field are selected Nonlinear effects on wave height distribu-
based on criteria such as minimum variations, tions are discussed in Tayfun and Fedele (2007).
minimum transient effects in the model test set- The nonlinear effects in the water waves and, in
up, or minimum reflections from the beach or particular, the vertical asymmetry are also re-
from the side walls. Normally a minimum of 10 sponsible for wave drift.
wave cycles are selected. Parameters are defined
Non-linear regular wave characteristics are
by a time-domain (zero crossing) approach or by
defined by the significant components at higher
a frequency domain (Fourier or harmonic) ap-
harmonics. The asymmetric wave geometry,
proach, mainly with focus on the basic har-
with increased crests (and decreased troughs)
monic. Simple root-mean-squared (RMS) anal-
ysis of elevation records is also applied.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 4 of 27
and associated local steepness, may have im- tank model with wavemaking and absorption ca-
portant consequences in practice, for instance in pabilities, and HOS-Ocean, an open boundary
the case of stability tests. Recent studies also solver to approximate wave field transformation
show that the bound harmonics are also respon- from an initial state. The HOS method owes its
sible for the generation of rogue waves, Fedele accuracy and efficiency to the pseudo-spectral
et al. (2016). approach in solution procedure and has been
proven to be more efficient than finite-base dis-
To distinguish these ‘real’ non-linear effects cretization models, Ducrozet et al. (2012 and
in open-sea wave fields from ‘parasitic’ labora- 2016). The solution procedure details and for-
tory-induced ones, comparisons of laboratory mulation of the problem can be found in the
generated wave profiles with theoretical/numer- original works of Dommermuth et al. (1987),
ical nonlinear reference models are helpful. In and West et al. (1987). As in the Schäffer (1996)
this regard, although referring to a two-dimen- case, nonlinear regular waves of arbitrary order
sional wave system, Henderson et al. (2006) can be generated using the HOS method by se-
showed the relevance of accounting for nonlin- lecting a single wave component. High degree
ear effects in wave generation in order to of accuracy can be achieved using the HOS
achieve a wave pattern that propagates with a model, Ducrozet et al. (2007).
time independent form.
Over the past three decades, significant pro-
For second-order nonlinear regular (peri- gress has been made on the developed of a non-
odic) wave generation for intermediate water linear Fourier analysis methodology (NLFA) to
depth and deep water conditions, the method characterize nonlinear ocean waves in arbitrary
presented by Schäffer (1993) can be used by se- water depth. The work by researchers on this
lecting a single frequency component only. The topic was documented comprehensively in a
numerical algorithm will automatically generate full-length text by Osborne (2010). For waves in
a “clean” second-order Stokes wave containing coastal regions, the Korteweg-de Vries (KdV)
the fundamental component plus its phase- equation governing shallow water wave behav-
locked second-order frequency component only. iour was used. For waves further offshore, the
The unintended and undesirable second-order nonlinear Schrodinger (NLS) equations govern-
free wave induced by incompatible geometry of ing intermediate to deep water was applied (Os-
the wave board will be removed from the wave borne, 2010). Both the KdV and NLS equations
field by a potential wave equal in amplitude but admit both regular waves as well as irregular
with a 180-degree phase shift, thus having a can- wave solutions, and thus are applicable to anal-
celling effect on the second-order free wave (see ysis and simulation of nonlinear waves in both
Schäffer, 1996). categories (Bruhl and Oumeraci, 2016; Mohtat
et al. 2018b)). For regular waves, the nonlinear
Another approach, which is more general periodic waves are cnoidal waves (in the shal-
and can accurately generate nonlinear Stokes low water case) and Stokes waves (in the inter-
waves of arbitrary order, is the higher-order mediate to deep water case) (Mohtat et al.,
spectral (HOS) method developed by the 2020).
LHEEA-ECN (Hydrodynamics, Energetics and
Atmospheric Environment Laboratory of Ecole A note of caution is in order here. Regular
Centrale de Nantes) group. Their High-Order plane waves are stable only up to a limited steep-
Spectral method (HOSM) includes two pack- ness ratio. For long-duration regular wave gen-
ages, namely HOS-NWT, a numerical wave
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 5 of 27
eration, as the wave steepness increases, stabil- absorption mechanism is related to wave break-
ity of the one spatial-dimension plane waves ing and viscous dissipation. These issues were
will deteriorate with time and the wave field generally addressed by Chakrabarti (1994) and
may become two dimensional and converge to a USACE (2002).
stable wave surface consisting of rectangular
cells. This nonlinear phenomenon was ex- An important aspect to improve the effi-
plained in Henderson et al. (2006) via the stabil- ciency of the wave absorption device is porosity
ity of the one spatial dimension (third-order) on the surface and inner layers of the devices.
nonlinear Schrodinger equation. Usually, sloped and parabolic beach with porous
layers can have good absorption properties for
long waves and a wide range of wave steepness
2.1.3 Confinement effects
by an optimization of the shape and porosity.
Wave generation in water of finite depths in- Studies on this topic are provided by Straub et
troduces additional effects relative to that in al. (2011) and Lean (1967).
deep water. Dispersion is depth-dependent, with
Most existing wave tanks are of rectangular
shorter wavelengths and reduced speed in de-
shaped with straight corners. In case of wave
creased depths, Newman, (1977), Toffoli et al.
tanks with wave generators on multiple sides,
(2005). This may lead to spatial variations due
two adjacent sides of wave paddles are used to
to refraction effects unless the bottom is per-
generate oblique waves for a wide range of inci-
fectly horizontal and flat. Fully nonlinear solu-
dence angles. The combined movement of the
tions for wave propagation over topography are
paddles changes the wave crests directions for a
provided in Kennedy and Fenton (1997).
continuum wave field generation. However, in
Non-linear wave-wave interactions increase the corner there is at least one paddle with re-
with reduced depth, with sharper peaks but also stricted movements or blocked to prevent colli-
larger set down effects and corresponding return sions. This discontinuous surface (or singular-
currents. ity) at the corner creates a border effect on the
wave field that propagates to the test section of
In a wave basin, reflections on side walls and the wave tank. The simplest solution to prevent
wave generator are present, resulting in unde- a discontinuity of the wave field is adopting a
sired influence on the required wave field. gradual movement fader in the paddles corner
Waves reflected by the wave absorber at the op- but this results in an incorrect wave amplitude
posite end of the tank limit the total duration of propagating to the test section. A more complex
the test. Once the residual energy reflected by solution is provided by Matsumoto & Hanzawa
the (imperfect) absorbers and accumulated in (2001) who modified the paddle movement to
tank becomes unacceptably high compared to maintain the wave field stability. An optimiza-
the wave energy of the spectrum, the experiment tion of paddle movements is obtained numeri-
shall be terminated. In computing the accumu- cally before the test execution in the tank.
lated energy, the effects of the side wall also
have to be accounted for.
Another option to overcome a discontinuity
Passive wave absorbers suppress the reflec- at the corner is to use a circular layout of wave
tions and presents a cost effective feasibility for paddles. In this case, the paddle angle changes
a wide range of wave tanks. The physical wave gradually connecting the two adjacent sides of
the tank. In order to prevent collision between
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 6 of 27
paddles during movements, a gap between pad- culation may develop. In such cases, the only so-
dles is introduced. Normally, a mechanism lution is to wait enough time before starting a
made by metal strips and fabric connection is new test. The time interval depends on the spec-
used to cover the gap between the paddles. ifications by the client as well as on the specific
test and dimensions of the facility.
It is also important to make the transition
from the wave generator to the wave absorption 2.1.5 Radiation and reflection from model,
parts as smooth as possible, regardless if the beaches, etc.
wave absorption is active or passive.
As anticipated, wave reflected from the side
2.1.4 Wave frequency and low frequency re- walls and/or waves radiated or reflected by the
flections model and subsequently reflected by either the
side walls and/or the damping beaches may en-
In order to limit the effects of confinement ter the measuring field.
on wave generation, wave absorbers are usually
employed to suppress or substantially reduce the For instance, in a two-dimensional wave
reflections for the side walls of the tank. As al- tank, three-dimensional waves, albeit of small
ready discussed, sloped and parabolic beaches amplitude, can be generated due to the interac-
with porous layers are generally used but their tion with the structure and/or tank walls. The
effectiveness in absorbing the incoming wave phenomenon is noticeable for single flap wave-
changes with frequency and steepness. makers: in such cases the effect can be sup-
pressed, or at least significantly reduced in a
Theoretical approaches for the design of given range of wavelength, by using straighten-
wave absorption devices in shallow water are ers in the upstream portion of the tank. Multi-
provided in Lean (1967) for different bottom segmented wavemakers, with paddles of small
shapes. Laboratory tests to measure the absorber size compared to the tank width, can be properly
performances for different wave absorption so- controlled to suppress the phenomenon. The use
lutions when varying wave frequencies and of control system also allows suppression of
steepness are shown in Straub et al. (2011). other disturbances (e.g. reflection from the
beach, transverse standing wave and reflections
Different techniques to measure the reflec- from the model test).
tion coefficients based on two or three wave
probe signals are provided in Isaacson (1991). A Active absorption by the wave generator
technique to separate incoming and reflected contributes to reduce these effects. When the
waves in a deep water towing tank based on the model reflects or irradiates waves, it creates
Doppler shift in frequency domain from a standing waves between the wave generator and
known velocity of a single measuring probe is the model. The interaction of incidence and re-
presented in Drzewiecki and Sulisz (2019). flected waves changes the wave field by non-
linear effects, even induces wave breaking.
Particularly for long wave tanks, long and Thus, active absorption in the wave generator is
very low frequency oscillations may occur. often used to decrease the resonance of stand-
Such long waves cannot be effectively damped ings waves, as presented in a study by de Mello
out neither by passive nor active absorbers. Be- et al. (2013). Similarly, cross–tank waves could
sides the seiche waves, in long wave tanks, cir-
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 7 of 27
occur as standing waves. To prevent such trans- ity, e.g. Tulin and Waseda (1999). The phenom-
verse waves, one of the tank sides could be used enon is partly reduced by the increased role
for passive or active wave absorbers. played by the viscous dissipation, Ma et al.
(2012) but a careful check of accuracy and re-
The time interval of 20 minutes between peatability of the wave quality is needed when
subsequent repeats is suggested in ITTC (2014) using relatively small scales.
as acceptable for a typical facility. However, the
time interval can be tank specific and it has to be
2.2 Irregular waves
long enough to pitch down the residuary tank
disturbance beneath 1% of the next target wave
height, as pointed out in the document. 2.2.1 Wave spectra
In Van Essen et al. (2014 and 2016) methods Sea states are generally specified by the
to avoid/suppress long standing waves in test fa- short-term variance spectrum S(f) or S(ω),
cilities are described. where f and ω are the frequency and the angular
frequency respectively. Primary spectral param-
eters are the significant wave height Hs, defined
2.1.6 Deviations from ideal conditions
as 𝐻𝐻𝑠𝑠 = 𝐻𝐻𝑚𝑚0 = 4�𝑀𝑀0 and a characteristic
Ideally, a regular wave model would require wave period, e.g. the peak period Tp or the zero-
a unidirectional periodic wave field with ampli- crossing period Tz defined from the spectrum as
tude, period and direction constant throughout 𝑇𝑇𝑧𝑧 = 𝑇𝑇𝑚𝑚02 2𝜋𝜋�𝑀𝑀0 ⁄𝑀𝑀2 , the mean wave period
time and space. In practice, deviations from the 𝑇𝑇1 = 𝑇𝑇𝑚𝑚01 2𝜋𝜋 𝑀𝑀0 ⁄𝑀𝑀1 where Mi is the i-th spec-
ideal situation are observed, for various reasons, tral moment (DNV, 2011) defined as
which are associated with wavemaker, basin and
∞
wave absorbing devices. 𝑀𝑀𝑛𝑛 = ∫0 𝜛𝜛 𝑛𝑛 𝑆𝑆(𝜛𝜛)𝑑𝑑𝑑𝑑 (1)
Model testing procedures must take these ef- 𝑆𝑆(𝜛𝜛) being the power spectral density.
fects into account, in one or several of the fol-
lowing ways: a) avoiding them, b) reducing Many widely-used models for the spectrum
them, c) documenting them and interpreting of waves measured at a point (without regard to
their effect on vessel and offshore structure re- wave direction) are of the Bretschneider (1959)
sponses. form
If the disturbances cannot be suppressed, 𝐴𝐴
𝑆𝑆(𝑓𝑓) = 𝑓𝑓 5
exp(− B⁄𝑓𝑓 4 ) (2)
their effects on the measurements can be re-
duced by choosing a proper combination of lo-
cation and time windows. In all cases, reflec- where f is the frequency and A and B are con-
tions or other disturbances have to be measured stants. Among this class are those referred to as
and documented. Pierson & Moskowitz (one- and two-parameter
forms), see ITTC (1999, 2002, 2017a, and
Another important aspect concerns the 2017b) and Liu et al. (2018). These presenta-
model scale. When reducing the scale of the tions differ only with respect to the parameters
problem, and thus the wavelength, the waves are that are used in determining A and B.
keener to develop natural modulational instabil-
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 8 of 27
Other spectra include those related to the In case the wave theory in use is not able to
basic Bretschneider form, for cases where there capture all non-linear phenomena and basin ef-
is a limited fetch (JONSWAP), a finite water fects and the wave spectrum is not correctly re-
depth (TMA), or a combination of a known produced with the required accuracy, a Black-
wind speed and limited fetch, Mitsuyasu (1972). Box approach can be used to account for the ef-
A specified duration of random simulations is fects of nonlinearities. The Black-Box approach
important to achieve stationary irregular wave is based on calibration of the control signal in
conditions: typically, a duration of 1 hour (full the frequency domain using the ratio between
scale) for seakeeping and 3 hours (full scale) in the required and realized wave energy spectrum
offshore engineering is applied. as presented in Drzewiecki (2018).
this case, the irregular waves for the shallow wa- wave crests become less important when the di-
ter case are a superposition of cnoidal waves rectional spreading increases. For very large
plus their nonlinear interaction terms, and the ir- spreading the distribution of extremes ap-
regular waves for the intermediate to deep water proaches the second-order distribution, DNV-
case are a superposition of Stokes waves plus GL (2015) and Toffoli and Bitner-Gregersen
their nonlinear interaction terms (Osborne, (2011).
2010).
2.2.5 Short-crested wave modelling
2.2.4 Distribution of extremes
Most likely, the multi-modal sea is also mul-
From linear wave theory under an ideal ran- tidirectional. Generally, the swell component is
dom sea state, for a narrow-banded wave spec- unidirectional with very little directional spread-
trum, the wave crests follow the Rayleigh distri- ing whereas the wind sea may be characterized
bution (see Longuett-Higgins, 1952). However, by a more substantial directional spreading with
because of the nonlinearity of real waves, the central (peak) direction oriented with a cer-
higher-order effects should be accounted for and tain angle to the swell. A proper account for all
thus the Tayfun distribution is more representa- the peaks of the frequency spectrum as well as
tive for wave heights, Tayfun and Fedele of the directions of the different wave systems is
(2007). A probability distribution based on sec- fundamental for an accurate prediction of the
ond-order wave simulations was proposed by ship response.
Forristall (2000), which represents the wave
crest distributions of moderately steep short- A more general and detailed discussion of
crested waves quite well. the laboratory modelling of multi-modal and
multi-directional wave spectra is provided in the
However, there may be cases in which ef- Guideline 7.5-02-07-01.1, whereas a discussion
fects of third and higher order are relevant. It is on the methods generally employed for the
shown in Buchner et al. (2011) that the tail of measurement of directional seas are reported in
the distribution rises substantially when increas- the following.
ing the wave steepness. This result confirms
what was found by Onorato et al. (2006 and For modelling purposes, the directional
2009) who observed that for long-crested, steep characteristics of waves are sometimes assumed
and narrow banded waves, the probability of to be uncoupled from their spectral properties,
finding an extreme wave can be underestimated and then the spectrum of waves travelling within
by more than one order of magnitude if only sec- a given range of headings is taken to be some
ond-order theory is considered. proportion of that measured at a point. On this
basis, the directional spectrum is usually pre-
Considering the effect on the wave height sented in the form
and particularly the crest height, a proper ac-
count for the higher order effects in large waves 𝑆𝑆(𝑓𝑓, 𝜃𝜃) = 𝑆𝑆(𝑓𝑓) 𝐺𝐺(𝜃𝜃), (3)
is essential for a correct estimate of air gaps, im-
pact loads as well as green water conditions. where the spreading function G depends only on
the direction θ. Its most common form is
Results based on experiments or numerical
simulations indicate that higher order effects on 𝐺𝐺(𝜃𝜃) = 𝐹𝐹(𝑠𝑠) cos2𝑠𝑠 [(𝜃𝜃 − 𝜃𝜃1 )⁄2], (4)
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 10 of 27
where θ1 is the predominant wave direction, and affect both the breaking occurrence as well as
s is an index that determines the width of the di- bubble entrainment.
rectional spread. In other forms of G(θ), the
power 2s is replaced by s, or the argument of the 2.2.7 Geographical consistency of wave spec-
cosine may omit the factor 1/2. In another ap- trum selection
proach it can be expressed just in terms of its
angular harmonics, see e.g. Frigaard et al., In order to achieve an accurate prediction of
(1997) for further details. The function the responses of ships or offshore structures to
wave loading, it is important that wave spectrum
22𝑠𝑠−1 Γ2 (𝑠𝑠+1)
𝐹𝐹(𝑠𝑠) = , (5) characteristics are correctly reproduced in the
𝜋𝜋 Γ(2𝑠𝑠+1)
laboratory. When performing tests for ships or
ensures that the total variance of the directional offshore platforms and devices required to oper-
spectrum S(f,θ) is the same as that of the point ate in certain regions, it is strongly recom-
spectrum S(f). mended that the specific wave conditions are
employed as long as they are producible.
Beside the changes in the distribution of ex-
tremes discussed in the previous section, the di- There is evidence that some places in the
rectional spreading also changes some statistical world are characterized by specific wave condi-
properties of the waves like skewness and kur- tions which may be seasonal. Examples of the
tosis, Toffoli et al. (2010). inter-annual variability of the wind-sea and
swell significant wave heights are presented in
Semedo et al. (2011) for different Oceans and
2.2.6 Wave breaking – influence on statistics noticeable differences in Pacific and Atlantic
and kinematics Oceans have been observed. During winter in
the Northern Hemisphere, a strong north–south
In irregular seas, the occurrence of wave
swell propagation pattern is observed in the At-
breaking limits the maximum steepness and re-
lantic Ocean.
duces the probability of occurrence of extreme
events, Buchner et al. (2011). Another very important example is repre-
sented by the West Africa region investigated in
In some cases, the occurrence of breaking
the WASP-JIP (West Africa Swell Project JIP).
prevents the achievement of the specific spec-
The wave spectra at sites off West Africa are
trum. In order to overcome this problem, differ-
dominated by the constant presence of one or
ent seeds are used to generate the random phases
more swell wave systems (Olagnon et al. 2013),
until a reduction of the breaking occurrence is
with swell spectra exhibiting a strong narrow-
obtained. However, it can also indicate that the
ness, both in frequency and direction. Again
selected theoretical spectrum is not an appropri-
concerning the seasonal variations Prevosto et
ate representation of the actual wave spectrum
al. (2013) show that during the austral winter
at that steepness, i.e. the theoretical input spec-
swells approach the West Africa coast from
trum might not be suitable to describe a very
south to south-westerly direction whereas dur-
steep sea state with a lot of occurrence wave
ing the austral summer north-westerly swells are
breaking, as also described in Huang et al.
observed.
(2018). When performing tests in model scale
the increased role played by surface tension Several examples are briefly summarized
changes some of the breaking feature and may below:
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 11 of 27
In the above equations, a0, ω0 and k0 are the Ocean waves often appear in sequences of
amplitude, frequency and wave number, respec- wave groups. The response of marine structures
tively, and cg is the group velocity which can be on successive wave groups is quite different
estimated as half the phase velocity. from regular waves and an isolated wave group.
The above solution is valid for deep water. An interesting study is under development at
More recently, the solution has been generalized the University of Michigan in which a large da-
to finite water depth. The solution, which is pro- tabase of field data has been analyzed with the
vided in Onorato et al. (2013), is basically the aim of identifying the wave groups in Seyffert
same as above but the equation of the envelope et al. (2016). These real time series could then
is a little more complicated: be the basis for generating an ensemble of wave
time series, all of which contain wave groups of
𝑞𝑞𝑝𝑝 (𝑥𝑥, 𝑡𝑡) = 𝑎𝑎0 exp(−𝑖𝑖𝑖𝑖𝑎𝑎02 𝑡𝑡) known runs and probability of occurrence. The
4𝛼𝛼(1−𝑖𝑖2𝛽𝛽𝑎𝑎02 𝑡𝑡) (8) study is not yet mature and consolidated to be
×� 2 − 1� used as a practical approach for testing, but it
𝛼𝛼+𝛼𝛼�2𝛽𝛽𝑎𝑎02 𝑡𝑡� +2𝛽𝛽𝑎𝑎02 𝑥𝑥 2
seems to be a very interesting and promising di-
where: rection to follow.
1
𝛼𝛼 = �𝑐𝑐 2 − 𝑔𝑔ℎ sech2 [𝑘𝑘0 ℎ](1 An envelope approach to simulate wave
2𝜔𝜔0 𝑔𝑔 groups was proposed by Xu et al. (1993). Two
− 𝑘𝑘0 ℎ tanh2 [𝑘𝑘0 ℎ])� parameters, i.e. the groupness factor height
(GFH) and groupness factor length (GFL) were
and introduced to represent the group. An empirical
wave envelope with given GFH and GFL is pro-
𝜔𝜔0 𝑘𝑘02 8 + cosh[4𝑘𝑘0 ℎ] − 2tanh[𝑘𝑘0 ℎ]2
𝛽𝛽 = � − posed in Liu et al. (2013) together with a gener-
2 8sinh[𝑘𝑘0 ℎ]4 ation procedure for random waves with deter-
2
�2𝜔𝜔0 cosh[𝑘𝑘0 ℎ]2 + 𝑘𝑘0 𝑐𝑐𝑔𝑔 � mined wave groupness.
�
(𝑘𝑘0 sinh[2𝑘𝑘0 ℎ])2 �𝑔𝑔ℎ − 𝑐𝑐𝑔𝑔2 �
2.4 Measurement and analysis of long- and
short-crested waves
2.3.2 Wave groupness
For unidirectional waves with assigned di-
Although extreme events are very important
rection, one wave probe is sufficient. If direc-
for estimation of the extreme loading, there is
tional measurements are required, an array of
difficulty in using them for design as there is not
wave probes is needed.
much indication about the probability of occur-
rence in DNV-GL (2015). Moreover, not neces- Several approaches can be used to estimate
sarily one single event, although extreme, repre- short-crested waves and directional spectrum,
sent the most critical situation. Considering the differing from each other by the minimum num-
coupling of the wave system with the response ber of wave probes needed and the method of
of the ship for instance, the passage of a wave analysis employed. The Maximum Entropy
group with a sequence of a sufficient number of Method (MEM) exploits the similarity between
waves of critical height and length may be even the directional distribution function and a prob-
more dangerous than a single wave, even if of ability density function. There are different ver-
higher amplitude.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 13 of 27
sions developed which seem to reproduce mul- The capacitive probe consists of two metal
tiple spectral peaks reasonably well, but as for rods immersed in water in parallel, where at
parametric models (and all other methods as least one rod is covered with an electrical insu-
well), the estimates will be affected by the actual lator and the capacitance measured between two
method and the way it is applied. Another mutually insulated rods is linearly dependent on
widely used group of methods is based on the the immersion depth, Clayson (1989).
principle of Maximum Likelihood. It is nor-
mally easier to use and more computationally The ultrasonic probe emits the ultrasonic
efficient than MEM. wave vertically downward, where the time
elapsed from emission to the absorption of the
Generally, MEM is more accurate but less ultrasonic wave reflected from the water eleva-
robust: it does not always converge to an esti- tion, is linearly dependent on the distance to wa-
mate and in some cases the estimate is far off. ter elevation, Payne (2008).
The MLM method is more robust (its estimates
are never far off) and takes far less time to com- The optical probe consists of the camera and
pute, Van den Berg (2011). The accuracy im- the laser, where the laser spot at the water eleva-
proves if wave slopes and elevation instead of tion is captured by the camera and images are
the point elevation measurements are employed, processed into a value of wave height based on
particularly for a small footprint array. the triangulation method, Payne (2009).
Another method, based on a complex Morlet The finger probe, consist of the servo-drive
wavelet transform, which can be used with just mechanism tracking the water elevation, where
three probes located arbitrarily, was proposed the motions of the mechanism, that correspond
by Donelan et al. (1996). As the data are first to the water elevation, are measured using a
resolved in time and frequencies by the wavelet servo-drive encoder, Payne (2008).
transform, the stratification of the spectral con-
tent of the data in the time domain can simplify Each of the indicated solution, has certain
the directional estimation procedure. Compared advantages and disadvantages. The rods of the
with the estimation methods of maximum like- resistive and capacitive probes are immersed
lihood (MLM) and maximum entropy (MEM), into water. Therefore, the rods can start vibrat-
the wavelet based method could give the best es- ing under the wave loads and disturb the meas-
timation of the directional spectra, Donelan et ured wave elevation, especially at forward
al. (2015). speed. Direct contact with water causes the sed-
iments to settle on the surface of the immersed
Several types of probes can be used to meas- rods which need to be cleaned regularly. The re-
ure the waves: resistive probe, capacitive probe, sistance and capacity between the rods are de-
optical probe, ultrasonic probe, finger probe, as pendent on the temperature and chemical com-
described in Payne (2008). Each of these probes position of water and thus it has to be included
has a different design and operation principle: in the calibration of the wave probes. Conse-
quently, it is recommended to calibrate the
The resistive probe consists of two thin probes with cleaned rods in their operating en-
metal rods immersed in water in parallel, where vironment and check the calibration regularly.
the resistance measured between the rods is lin- The calibration determines the relation of the
early dependent on the immersion, Clayson probe output signal to the height of the wave
(1989). profile permeating the rods. The height of the
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 14 of 27
wave profile for calibration is simulated by the order 6 HOSM model with formulation of West
vertical immersion of the rods into the water. et al. (1987) corresponds to a nearly fully non-
The calibration has to be performed by measure- linear model, Clamond et al. (2006).
ment of the output signal of probe at enough
points of reference immersion depth. The num- In the pseudo-spectral method of solution,
ber of points should be large enough to deter- computation of some of the equations, usually
mine a relationship based on a mathematical including products of variables, are performed
function with a satisfactory accuracy. The pre- in the physical domain and the rest, usually in-
ferred relationship is given by a linear function. cluding the derivatives of the variables, in the
The scope of the calibration has to be deter- Fourier space. This approach exhibits some
mined in accordance with the measurement amazing convergence properties. The addition
scope of the probes and generation scope of the of a second-order wave maker theory in the
wave profiles. HOS-NWT makes this package even more suit-
able for validation and prediction of the wave
The tips of the ultrasonic and optical probes conditions. Details of the additional wave maker
are not in contact with water. Thus, they are not boundary condition and the solution procedure
exposed to the water and do not intrude the is well defined in Ducrozet et al. (2012).
measured waves. However, the ultrasonic
probes have a specified operating range of the Nonlinear ocean wave governed by the KdV
wave steepness to be reliable, Payne (2008) and for shallow water and the NLS for intermediate
the speed of sound in the air at given tempera- to deep water will appear as non-stationary in
ture and humidity has to be included into cali- the linear spectral sense as they contain Stokes
bration. waves and breathers (phase-locked Stokes wave
modes), both of which contain phase-locked
harmonic components. However, the wave
2.5 Non-stationary power spectrum (time
modes are stationary when viewed from a non-
and space)
linear Fourier analysis perspective (Osborne,
In recent years, there have been an improved 2010).
implementation of HOS packages that provide a
large database of validation available and the 2.6 Calibration
code can be considered matured for engineering
practices, Ducrozet et al. (2016). Some of the The wave environment needs to be cali-
applications that could be more appealing to the brated prior to the test to ensure the correct en-
energy industry and deep-water wave modelling vironment modelling has been achieved. The ef-
are, but not limited to, modulational instabili- fective duration of the wave field with transient
ties, Fernandez et al. (2014) and Toffoli et al. part removed needs to be sufficiently long as
(2010) and freak waves, Ducrozet et al. (2007) specified. The required tolerances for wave cal-
and Xiao et al. (2013). It should be noted that ibration are usually ±5% for both significant
HOSM is a strong nonlinear model can include wave height and peak wave period.
the wave-wave interactions up to a defined or-
der, which is called the order of HOSM. For in- The shape of the wave spectrum can be cali-
stance with an order of 3 HOSM model, an ac- brated comparing the measured wave spectrum
curacy equal to that of Zakharov equation in with the theoretical wave spectrum at the target
Zakharov (1968) can be achieved, Sergeeva et location in the basin. Correction factors that can
al. (2013). Further investigation revealed that an be applied to the wavemaker steering signal can
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 15 of 27
be calculated based on the frequency-by-fre- wave heights and reduces the wavelength. The
quency ratio between the measured and required problem becomes more complex in case of a
spectrum. depth-varying current (Ellingsen and Li (2017))
and in the open ocean the wave-current interac-
DNV-RP-C205 (2010) recommends to use tion can affect the modulation of the spectrum,
the Forristall distribution for stationary short- the wave breaking rates, and wave statistics
term storms noting that the distribution is based (Romero et al. (2017)).
on second-order wave simulations and extremes
predicted by this distribution may be on the low In the ocean, negative horizontal velocity
side. The criteria of the maximum wave crest for gradients (i.e. an accelerating opposing current
wave calibration test has not been established to or a decelerating following current) make waves
date. Huang (2017) and Liu and Zhang (2018) shorten and heighten which enhances wave
adopted a CFD and a Higher-Order Spectral steepness. As a result, a nonlinear mechanism
(HOS) method, respectively, to investigate this known as modulational instability develops,
subject. They found that third-order nonlinearity leading to the formation of large amplitude
should be taken into account when wave calibra- waves (the so-called rogue waves), even if they
tion is performed for long-crest waves in a wave would otherwise be unexpected. In Toffoli et al.
basin with maximum wave crest being a key cri- (2019), laboratory experiments and numerical
terion. Liu and Zhang (2018) also proposed an simulations with a current-modified version of
empirical model to estimate the maximum wave the Euler equations are presented to assess the
crest height for wave calibration. It must be role of an opposing current in changing the sta-
noted that the resulting crest height in a 3 hour tistical properties of unidirectional random
realization is dependent on the randomly distrib- wave fields. Results demonstrate in a consistent
uted phases. Especially, for very steep sea states and robust manner that an opposing current in-
the crest distribution of different 3 hour realiza- duces a sharp and rapid transition from weakly
tions can vary from one another and from theo- to strongly non-Gaussian properties with a con-
retical distributions. However, this does not sequent increase of the probability of occurrence
mean that a 3 hour realizations with large devi- of rogue waves. The tests were conducted with
ations from high order distributions can be ne- irregular unidirectional waves in a wave flume
glected, but is the result of the random 3 hour and a directional wave basin at Plymouth Uni-
realization of a certain sea state. A better repre- versity. The initial conditions at the wave maker
sentation of the crest distributions of a certain were given in the form of an input JONSWAP-
sea state can be achieved by running multiple 3 like wave spectrum to model waves in the fre-
hour realizations, to fill in the tail of the distri- quency domain. As the wave field entered into a
bution, see also Huang et al. (2018). region of opposing current, the wave height was
observed to increase. Evident breaking dissipa-
2.7 Interactions with and current and wind tion was observed for very strong current fields
for U/Cg > 0.3 (breaking appeared with even
less strong currents in the wave basin). The pres-
2.7.1 Waves on currents ence of the current also accelerated the down-
shift of the spectral peak, with energy migrating
Theoretically, a perfectly steady current that
from high to low frequencies bands within
is collinear with the waves slightly reduces the
scales of tens of wavelengths, in agreement with
wave heights and increases the wavelength.
modulational instability effects, Onorato et al.
Similarly, an opposing current increases the
(2009). The analysis of the statistical properties
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 16 of 27
of extreme (rogue) waves with an aid of kurto- Aside from specific studies on wind gener-
sis, the fourth-order moment of the probability ated waves, e.g. Buckley and Veron (2018), in
density function of the surface elevation. The the most practical applications the wind fetch in
kurtosis expresses the probability of extreme laboratory is too short to affect the wave dynam-
events in a record (this assumes the value of 3 ics substantially. In such a case, it is recom-
for Gaussian sea states). The sea state rapidly mended that the wave spectrum has to be chosen
transitioned from a weakly to a strongly non- to match the open ocean conditions, Babanin et
Gaussian condition as current speed increased; al. (2019).
maximum values of kurtosis were detected to
reach 3.5 or higher, which are remarkably high
for water waves. These features were evident in 3. GENERATION TECHNIQUES
both facilities. However, the wave basin exhib-
ited much higher kurtosis (> 4) than the wave 3.1 General aspects of wave generation tech-
flume. Agreement with numerical simulations niques
confirms that this transformation can be at-
tributed to quasi-resonant nonlinear interactions Many different type of wave makers are in
triggered by the background current. use, typically piston or flap type wave genera-
tors with hydraulic or electrical power units suit-
Generally, the specified wave conditions are able for servo-controleld actuation of the wave
calibrated with the current, if current is required boards. Through the motion of the wave board
for the tested location. If waves are calibrated on at prescribed amplitude and frequencies result-
current, changes in the wave height and wave ing in a displacement of the water in front of the
length are accounted for and embedded in the wave board, progressive waves are generated.
resulting wave field. Nonlinear wave-current in-
teraction effects influence the resulting hydro- The hydrodynamic transfer function re-
dynamic forces, such as wave drift damping quired to calculate the required wave board mo-
(and corresponding modification in slow-drift tions has been deduced by Biesel (1951) and ex-
excitation), wave-induced currents and wave perimentally confirmed, for a piston type wave
kinematics. maker, by Ursell (1959). An overview of the ap-
plication of the Biesel hydrodynamic transfer
2.7.2 Waves on wind function is also given in Schmittner et al. (2013).
In most of the cases, waves are the ultimate Piston type wave makers are typically used
effect of the action of the wind over the ocean in water of finite depth, as the entire water col-
surface. The wind modifies the kinematics and umn is moved by the wave maker. And flap type
dynamics of the wave group. In the specific case wave maker are typically used in intermediate to
of extreme waves evolving under the collinear deep water. The hinge depth or underwater
wind, the focus point shifts downstream, Kharif height of the wave board in combination with
et al. (2008). The airflow separation at the wave the maximum amplitude or angle determine the
crest, beside contributing to the momentum working range of the wave maker. Wave maker
transfer, sustains the steep waves, that subse- with large hinge depth can generate larger
quently breaks reaching larger steepnesses com- waves. On the other hand, when generating
pared to the no-wind condition, Iafrati et al. short waves, the mismatch between orbital mo-
(2019). tions and wave board geometry can become too
large, resulting in spurious waves.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 17 of 27
The water depth has to be considered accu- analysis, stability in time should also be docu-
rately to determine the adequate transfer func- mented, as well as stability in space whenever
tion. Depending on the particular frequency, a relevant.
slight reduction of water depth of about 5 %
might result in a wave height reduction up to 50 When regular waves are generated in wave
%. basins, some irregularities may appear due to
discontinuities in the paddle distribution. This
Second order or higher order mitigation can be either due to the end of the wave maker
methods can be applied to reduce the effect of a or to the transition from the wave making to the
mismatch between required orbital motions and wave absorbing paddles. A solution to suppress
flap motions. References describing these meth- or mitigate the problem is provided by Matsu-
ods are given in the following sections. moto & Hanzawa (2001) who introduce a spe-
cific correction to each paddle movement in or-
Directional waves can be generated with der to achieve a more uniform wave field in
multi-flap wave makers. For more details on di- front of the wave maker. The method is well
rectional wave generation refer to Section 3.3.3 documented and can be used to different wave
of this guideline. A more general and detailed maker configurations.
discussion of the laboratory modelling of multi-
modal and multi-directional wave spectra is pro- Regular waves of arbitrary order can be gen-
vided in the Guideline 7.5-02-07-01.1 erated using the HOS method as described in a
series of papers by Ducroset and collaborators
3.2 Regular wave generation (see Ducroset et al. 2007, 2012 and 2016). For
the NLS-based regular wave generation in inter-
The quality of generated regular waves mediate to deep water conditions proposed by
should be carefully monitored as quality degra- Mohtat et al. (2020) by restricting the total num-
dation may occur quickly during propagation, ber of nonlinear modes to just one.
Benjamin and Feir (1967) and Stansberg (1993).
Also, attention has to be paid to the passive or 3.3 Irregular wave generation
active wave absorption of the beach and other
boundaries. For the generation of regular waves, 3.3.1 Unidirectional wave generation
second-order generation techniques for irregular
waves, Sulisz and Hudspeth (1993) and Schäffer Unidirectional (or long-crested, 2D) irregu-
(1996) can be used. Although second-order cor- lar waves are frequently used in most model ba-
rections are very important in the case of steep sins not only because they represent real sea-
waves, the second-order correction of the flap states in a very simple form, but also because it
motions compensate for waves generated due to is assumed that they result in conservative load-
a mismatch of the flap motion with the wave or- ings and responses compared to short-crested
bital velocity. This is particularly relevant for (directional) seas. It is also easier to define a sea
very long and shallow water waves. state in a unique manner.
Depending upon the application, documen- In the generation of 2D irregular waves, it is
tation of possible deviations from ideal condi- important to maintain “true” randomness that
tions such as reflections should be made availa- will prevent unrealistic repetition of the waves.
ble from the tests, see the 22nd ITTC Report on
Environmental Modelling (1999). In the wave
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 18 of 27
Also, careful attention should be given to the ef- waves. The longer the duration of wave genera-
fects of the frequency range covered by the tion (determined by the nature of the model
servo system. The test duration and the number tests), the more frequency components are
of frequency components adopted are also im- needed. The specified duration of random simu-
portant if proper natural statistics of the wave lations is important, normally 3 hours (full
field are to be reproduced. Wave reflection from scale) for modelling a full storm, which is most
the beach and diffraction by the basin walls often used in offshore engineering tests. It can
should be monitored carefully. be changed depending on the phenomena the
test is focusing on, however it must be long
Linear irregular random waves can be mod- enough to realize statistical properties if nonlin-
elled as a summation of sinusoidal wave com- earities and extremes are to be studied. For sea-
ponents or as filtering of white noise. keeping tests, at least 100~200 waves have tra-
ditionally been used (typically 0.5~1 hour, full
The linear long-crested wave model with the scale), which is often defined as satisfactory if
summation of sinusoidal wave components is linear effects only are considered.
given by:
Given the characteristic function of the wave which represents the Fourier transform of the
generator, H(f), the problem is to find a function wave maker control time history.
y(t) to be used as input to the wave maker in or-
der to obtain the desired spectral density func- The white noise approach has the advantage
tion to be realized, Sz(f). The idea behind the of generating non-repeating records.
white noise generation approach is that the func-
tion y(t) can be obtained from w(t) through a 3.3.3 Directional wave generation
specific filter Q(f). The filter Q(f) may be
viewed as the inverse of that needed for whiten- Many basins now use multi-directional wave
ing the function y(t). generators to achieve more realistic wave envi-
ronments. Wave generators usually consist of
Hence, if Z(f) is the desired frequency con- many small wave boards, which can be con-
tent of the wave system to be generated, it is ob- trolled independently by electric or electric-hy-
tained as: draulic actuators.
𝑍𝑍(𝑓𝑓) = 𝐻𝐻(𝑓𝑓) ∙ 𝑄𝑄(𝑓𝑓) ∙ 𝑊𝑊(𝑓𝑓)
Due to the effects of the Biesel limit on the
and thus the frequency content of y(t) is size of the wave board, Biesel (1954) and reflec-
tion from the walls, wave characteristics in the
𝑌𝑌(𝑓𝑓) = 𝑄𝑄(𝑓𝑓) ∙ 𝑊𝑊(𝑓𝑓) test region need to be carefully determined. As
a reference, the relation of the length of each el-
Correspondingly, the spectral density func- ement of wave makers and the angle of emission
tions are related by: of waves is practically shown as follows:
𝑆𝑆𝑧𝑧 (𝑓𝑓) = |𝐻𝐻(𝑓𝑓)|2 |𝑄𝑄(𝑓𝑓)|2 𝑆𝑆𝑤𝑤 (𝑓𝑓) 𝜆𝜆
𝑙𝑙 =
As already stated, Sw(f) can be assumed to be √2 + sin𝜃𝜃
unity and then:
Here, 𝑙𝑙 is the length of each element of wave
𝑆𝑆𝑧𝑧 (𝑓𝑓) = |𝐻𝐻(𝑓𝑓)|2 |𝑄𝑄(𝑓𝑓)|2 makers, 𝜆𝜆 is the generated wave length and 𝜃𝜃 is
the angle of emission of waves.
which leads to:
𝑆𝑆𝑧𝑧 (𝑓𝑓) Modelling directional spectra in the labora-
|𝑄𝑄(𝑓𝑓)|2 = tory is generally associated with a significant
|𝐻𝐻(𝑓𝑓)|2
random scatter, especially in the finer features
By introducing the additional constraint that of the measurements. This reflects features of
Q(f) has to be a real function, the above equation real sea data, reflecting natural statistical scatter,
finally provides Stansberg (1998). Therefore, a robust descrip-
tion of the directional sea conditions is often re-
�𝑆𝑆𝑧𝑧 (𝑓𝑓) stricted to a few parameters only, such as mean
𝑄𝑄(𝑓𝑓) =
|𝐻𝐻(𝑓𝑓)| direction, directional spread, and a simple shape
parameter that expresses the bimodality such as
and then: skewness and kurtosis, see Kuik et al. (1988)
and Stansberg (2001), or simplified smooth par-
�𝑆𝑆𝑧𝑧 (𝑓𝑓)
𝑌𝑌(𝑓𝑓) = 𝑊𝑊(𝑓𝑓) ametric models such as the cos2s(θ) model.
|𝐻𝐻(𝑓𝑓)|
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 20 of 27
For second-order directional wave genera- Biesel, F., 1954, “Wave Machines”, Proceed-
tion the method by Schäffer and Steenberg ings of the First Conference on Ships and
(2003) can be used. Ducrozet et al. (2007, 2016) Waves, Hoboken, New Jersey, pp.288-304.
also have developed open-source HOS numeri-
cal codes for directional wave generation. Cur- Bretschneider, C. L., 1959, “Wave Variability
rently there are no known wavemaker theory for and Wave Spectra for Wind-generated Grav-
extensions of the NLS equation to directional ity Waves”, Beach Erosion Board, U.S.
wave generation, which is governed by the Army Corps of Engineers, Washington
Davey-Stewartson equation (see Davey and (DC), USA. Technical Memorandum No.
Stewartson (1974) and Osborne (2010)). 118.
Babanin, A.V., Rogers W.E., de Camargo, R., Bruhl, M. and Oumeraci, H., 2016, “Analysis of
Doble, M., Durrant, T., Filchuk, K., Ewans, Long-Period Cosine-Wave Dispersion in
K., Hemer, M., Janssen, T., Kelly-Gerreyn, Very Shallow Water Using Nonlinear Fou-
B., Machutchon, K., McComb, P., Qiao, F., rier Transform Based on KdV Equation,”
Schulz, E., Skvortsov, A., Thomson, J., Vi- Applied Ocean Research, Vol.61, 81-91.
chi, M., Violante-Carvalho, N., Wang, D.,
Waseda, T., Williams, G. and Young, I.R., Buchner, B., Forristall, G., Ewans, K., Christou,
2019, “Waves and Swells in High Wind and M. and Hennig, J., 2011, “New Insights in
Extreme Fetches, Measurements in the Extreme Crest Height Distributions (A Sum-
Southern, Ocean Front”, Maritime Science, mary of the ‘CREST’ JIP)”, Proceedings of
Vol. 6, art. 361. DOI: the 30th International Conference on Ocean,
10.3389/fmars.2019.00361. Offshore and Arctic Engineering
OMAE2011-49846, Rotterdam, The Nether-
Belenky, V., Bassler, C. C. and Spyrou, K. J. , lands.
2011, “Development of Second Generation
Intact Stability Criteria”, Hydrodynamics Buckley, M.P. and Veron, F., 2018, “The Tur-
Department Report, Naval Surface Warfare bulent Airflow Over Wind Generated Sur-
Center. face Waves”, European Journal of Mechan-
ics/ B Fluids, Vol. 73, pp. 132-143. DOI:
Benjamin, T.B., and Feir, J.E., 1967, “The Dis- 10.1016/j.euromechflu.2018.04.003.
integration of Wave Trains on Deep Water.
I Theory”, Journal of Fluid Mechanics, Vol. Chabchoub, A., Hoffmann, N. and Akhmediev,
27, pp. 417-430. N., 2012, “Observation of Rogue Wave
Holes in a Water Wave Tank”, Journal of
Biesel, F., Suquet, F., 1951, ‘Les appareils gen- Geophysical Research, Vol. 117, C00J02.
erateurs de houle en laboratoire.’, La Houille
Blanche, (2):147-165, Mars-Avril 1951.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 21 of 27
Chakrabarti, S. K., 1994, “Offshore Structure C., Nishimoto, K., 2013, “A Control and Au-
Modeling”, World scientific. tomation System for Wave Basins”, Mecha-
tronics, Oxford, Vol. 23, pp. 94-107.
Chen, J., Li M., Wang W., 1990, “A Study on
Typhoon Wave Spectrum in the Sea Area DNV, 2010, Environmental Conditions and En-
Adjacent Hong Kong”, Tropic Oceanology, vironmental Loads, Recommended Practice,
9(4): 1-8. (In Chinese) DNV-RP-C205.
Clauss, G.F., Hennig, J., 2003, “deterministic DNV, 2011, Modelling and Analysis of Marine
Analysis of Extreme Roll Motions and Sub- Operations, Recommended Practice, DNV-
sequent Evaluation of Capsizing Risk”, In RP-H103.
STAB 2003 – 8th International Conference
on the Stability of Ships and Ocean Vehi- DNV-GL, 2015, “Rogue Waves: Impact on
cles, Madrid, Spain. Ships and Offshore Structures”, DNV-GL
Strategic Research and Innovation, Position
Clamond, D., Francius, M., Grue, J. and Kharif, Paper 05-2015.
C., 2006, “Long time Interaction of Enve-
lope Solitons and Freak Wave Formations”, Dommermuth, D.G. and Yue, D.K., 1987, “A
European Journal of Mechanics-B/Fluids, High-order Spectral Method for the Study of
25(5), pp.536-553. Monlinear Gravity Waves”, Journal of Fluid
Mechanics, 184, pp.267-288.
Clayson, C., 1989, “Survey of Instrumental
Methods for the Determination of High-fre- Donelan, M.A., Babanin, A., Sanina, E. and
quency Wave Spectrum”, Technical Report Chalikov, D., 2015, “A Comparison of
267, Institute of Oceanographic Sciences Methods for Estimating Directional Spectra
Deacon Laboratory. of Surface Waves”, Journal of Geophysical
Research Oceans, 120, 5040–5053. Doi:
Cuong, H.T., Troesch, A.W. and Birdsall, T.G., 10.1002/2015JC010808.
1982, “The Generation of Digital Random
Time Series”, Ocean Engineering, Vol. 9, Donelan, M.A., Drennan, W.M., and Magnus-
pp. 581-588. son, A.K., 1996, “Non-stationary Analysis
of the Directional Properties of Propagating
Davey, A. and Stewartson, K., 1974, “On Three- Waves”, J. Physical Oceanography, 26,
Dimensional Packets of Surface Waves,” 1901–1914.
Proc. Royal Society of London A, Vol.338,
101-110. Drzewiecki, M. and Sulisz, W., 2019, “Genera-
tion and Propagation of Nonlinear Waves in
Dean, R.G. and Dalrymple, R.A., 2000, “Water a Towing Tank”, Polish Maritime Research,
Wave Mechanics for Engineers and Scien- 1(101), Vol. 26, pp. 125-133, doi:
tists”, World Scientific. 10.2478/pomr-2019-0014.
de Mello, P. C., Carneiro, M. L., Tannuri, E. A., Drzewiecki, M., 2018, “Control of the Waves in
Kassab, F., Marques, R. P., Adamowski, J. a Towing Tank with the Use of a Black-Box
Model”, Zeszyty Naukowe Wydziału El-
ektrotechniki i Automatyki Politechniki
Gdańskiej, No. 59, 37-42.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 22 of 27
Ducrozet, G., Bonnefoy, F., Le Touzé, D. and Depth”, Natural Hazards and Earth System
Ferrant, P., 2007, “3-D HOS Simulations of Sciences, 14(3), pp.705-711.
Extreme Waves in Open Seas”, Natural Haz-
ards and Earth System Science, 7(1), pp.109- Forristall, G. Z., 2000, “Wave Drest Distribu-
122. tions: Observations and Second-order The-
ory”, J. Physical Oceanography, Vol. 30, pp.
Ducrozet, G., Bingham, H.B., Engsig‐Karup, 1931-1943.
A.P., Bonnefoy, F. and Ferrant, P., 2012, “A
Comparative Study of Two Fast Nonlinear Frigaard, P., Helm-Petersen, J., Klopman, G.,
Free‐surface Water Wave Models”, Interna- Stansberg, C.T., Benoit, M., Briggs, M.J.,
tional Journal for Numerical Methods in Flu- Miles, M., Santas, J., Schäffer. H.A., and
ids, 69(11), pp.1818-1834. Hawkes, P., 1997, “IAHR List of Sea State
Parameters – an Update for Multidirectional
Ducrozet, G., Bonnefoy, F., Le Touzé, D. and Waves”, Proceeding of the IAHR Seminar
Ferrant, P., 2016, “HOS-ocean: Open-source on Multidirectional Waves and their Interac-
Solver for Nonlinear Waves in Open Ocean tion with Structures, 27th IAHR Congress.
Based on High-Order Spectral Method”,
Computer Physics Communications, 203, Henderson, D.M., Patterson, M.S. and Segur,
pp.245-254. H., 2006, “On the Laboratory Generation of
Two-dimensional, Progressive, Surface
Ellingsen, S. Å., and Li, Y., 2017, “Approxi- Waves of Nearly Permanent Form on Deep
mate Dispersion Relations for Waves on Ar- Water”, Journal of Fluid Mechanics, Vol.
bitrary Shear Flows”, Journal of Geophysi- 559, pp. 413-427.
cal Research: Oceans, 122, pp. 9889–9905,
https://doi.org/10.1002/2017JC012994. Hemer, M., Katzfey, J. and Hotan, C., 2011,
“The wind-wave climate of the Pacific
Fan S. and Wang T., 1992, “On the Characteris- Ocean”, Report for the Pacific Adaptation
tics of Ocean waves—an Investigation of Strategy Assistance Program, Department of
Waves of the Western Pacific Ocean”, Climate Change and Energy Efficiency, The
Transactions of Oceanology and Limnology, Centre for Australian Weather and Climate
(1):1-11. (In Chinese) Research.
Fedele, F., Brennan, J., de León S.P., Dudley, J. Huang, Z.J., Guo Q., 2017, “Semi-empirical
and Dias, F., 2016, “Real world ocean rogue Crest Distributions of Long-crest Nonlinear
waves explained without the modulational Waves of Three-hour Duration”, Proceed-
instability”, Scientific Report, Vol. 6, 27715. ings of the 36th International Conference on
Ocean, Offshore and Arctic Engineering,
Fedele, F. and Tayfun, M.A., 2009, “On Nonlin- OMAE 2017- 61226, Trondheim, Norway.
ear Wave Groups and Crest Statistics”, Jour-
nal of Fluid Mechanics. 620, 221–239. Huang, Z., Zhang, Y., 2018, “Semi-Empirical
Single Realization and Ensemble Crest Dis-
Fernández, L., Onorato, M., Monbaliu, J. and tributions of Long-Crested Nonlinear
Toffoli, A., 2014, “Modulational Instability Waves”, Proceedings of the 37th Interna-
and Wave Amplification in Finite Water tional Conference on Ocean, Offshore and
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 23 of 27
ITTC, 1999, “The Special Committee on Envi- Lan, C., 1984, “A Preliminary Study on Ty-
ronmental Modelling Final Report and Rec- phoon Generated Waves in the Sea Area
ommendations to the 22nd ITTC”, Proceed- Near Hong Kong”, Tropic Oceanology,
ings of the 22nd ITTC, pp. 1-38. 3(3):10-17. (In Chinese)
ITTC, 2002, “The Special Committee on Waves Lean, G. H., 1967, “A Simplified Theory of Per-
Final Report and Recommendations to the meable Wave Absorbers”, Journal of Hy-
23rd ITTC”, Proceedings of the 23rd ITTC, draulic Research, 5:1, 15-30, DOI:
pp. 505-736. 10.1080/00221686709500187.
ITTC, 2014, “Seakeeping Experiments”, Proce- Li, H., 1992, “Application of Wave Spectrum in
dure 7.5-02-07-02.1. Ocean Engineering”, Research and Develop-
ment of South China Sea, (2): 32-38. (In
ITTC, 2005, “The Special Committee on Stabil-
ity in Waves Final Report and Recommen- Chinese)
dations to the 24th ITTC”, Proceedings of
Liu, S., Zhang, and X., 2018, “Extreme Wave
the 24th ITTC, pp. 369-408
Crest Distribution by Direct Numerical Sim-
ITTC, 2017a, “The Special Committee on Mod- ulations of Long-crest Nonlinear Wave
elling Environmental Conditions Final Re- Fields”, App. Ocean Res., 86, 141-153.
port and Recommendations to the 28th
Longuet-Higgins , M. S., 1952, “On the Statisti-
ITTC”, Proceedings of the 28th ITTC, Vol-
cal Distribution of the Heights of Sea
ume II.
Waves”, Journal of Marine Research (11),
ITTC, 2017b, “Laboratory Modelling of Multi- 245–266.
directional, Irregular Wave Spectra”, Guide-
line 7.5-02-07-01.1.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 24 of 27
Ma, Y., Dong, G., Perlin, M., Ma, X. and Wang, Newman, J.N., 1977, “Marine Hydrodynam-
G., 2012, “Experimental Investigation on the ics”, MIT Press.
Evolution of the Modulation in Stability
with Dissipation”, Journal of Fluid Mechan- Ochi, M.K. and Hubble, E.N., 1976, “Six Pa-
ics, Vol. 711, pp. 101-121. rameter Wave Spectra”, Proceedings of the
15th Coastal Engineering Conference, pp.
Matsumoto, A., and Hanzawa, M., 2001, “New 301-328.
Optimization Method for Paddle Motion of Olagnon, M., Ewans, K., Forristall, G. and
Multi-directional Wavemaker”, Coastal En- Prevosto, M., 2013, “West Africa swell
gineering Proceedings, 1(25). doi: spectral shapes”, Proceedings of the 32nd
https://doi.org/10.9753/icce.v25.%p. International Conference on Ocean, Off-
shore and Arctic Engineering,
Mitsuyasu, H., 1972, “The One-dimensional OMAE2013-11228, Nantes, France.
Wave Spectra at Limited Fetch”, Proceed-
ings of the 13th Coastal Engineering Confer- Onorato, M., Proment, D., Clauss, G. and Klein,
ence, Vol. 1, pp. 289-306. M., 2013, “Rogue Waves: from Nonlinear
Schrödinger Breather Solutions to Sea-keep-
Mohtat, A., Yim, S.C., and Osborne, A.R. ing Test”, PLoS ONE, Vol. 8, e54929.
2018a, “Energy Content Analysis of Closed
Basin Wave Simulation Using Linear and Onorato, M., Osborne, A.R., Serio, M., Cava-
Nonlinear Fourier Analysis,” Proceedings of leri, L., Brandini, C. and Stansberg, C.T.,
the International Conference on Ocean, Off- 2006, “Extreme Waves, Modulational inSta-
shore and Arctic Engineering, OMAE2018- bility and Second-order Theory: Wave
78754, Madrid, Spain. Flume Experiments on Irregular Waves”,
European Journal of Mechanics-B/Fluids,
Mohtat, A., Yim, S.C., Osborne, A.R. and Chen, Vol. 25, pp. 586-601.
M., 2018b, “Nonlinear Deepwater Extreme
Wave Height and Modulation Wave Length Onorato, M., Cavaleri, L., Fouques, S., Gram-
Relation,” Proceedings of the ASME 397th stad, O., Janssen, P. A. E. M., Monbaliu, J.,
International Conference on Ocean, Off- Osborne, A.R., Pakozdi, C., Serio, M.,
shore and Arctic Engineering, OMAE2018- Stansberg, C.T., Toffoli, A. and Trulsen, K.,
18317, Fort Lauderdale, USA. 2009, “Statistical Properties of Mechani-
cally Generated Surface Gravity Waves: a
Mohtat, A., Yim, S.C., Adami, N. and Lomon- Laboratory Experiment in a Three-dimen-
aco, P., 2020, “A General Nonlinear Wave- sional Wave Basin”, Journal of Fluid Me-
maker Theory for Intermediate to Deepwater chanics, Vol. 627, pp. 235-257.
waves Using Inverse Scattering Transform,”
Proceedings of the ASME 37th International Osborne, A.R., 2010, Nonlinear Ocean Waves
Conference on Ocean, Offshore and Arctic and the Inverse Scattering Transform, aca-
Engineering, OMAE2018-78755, Madrid, demic Press, New York.
Spain.
Parente, C. E., Nogueira, I. C. M., Martins, R.
Mori, N., Liu, Paul C., Yasuda T., 2002, “Anal- P., Ribeiro, E. O., Grossman-Matheson, G.
ysis of Freak Wave Measurements in the Sea S., 2017, “Regional Environmental Charac-
of Japan”, Ocean Engineering, 29(11). terization of the Campos Basin, Southwest
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 25 of 27
Shu, F., Cai F., He J., Chen Z., Liao K., 2012, Tayfun, M. A. and Fedele, F., 2007, “Wave-
“Wave Spectrum Analysis of Guangdong height Distributions and Nonlinear Effects”,
Yangxi Offshore Area”, Journal of Ocean- Ocean Engineering, Vol. 34, pp. 1631-1649.
ography in Taiwan Strait, 31(4): 557-564. (
In Chinese) Toffoli, A., Babanin, A., Bitner-Gregersen, E.,
Magnusson, A.K., Monbaliu, J., Portilla, J.,
Spyrou, K. J., 2001, “The Nonlinear Dynamics Dumon, G. and Onorato, M., 2005, “On the
of Ships in Broaching”, The Annals of Ma- Shape of Nonlinear Wave Trains”, Proceed-
rie Curie Fellowship Association, Volume 1. ing of the Fifth International Symposium on
Ocean Waves Measurements and Analysis.
Stansberg, C.T., 1993, “On Spectral Instabilities
and Development of Nonlinearities in Prop- Toffoli, A., Gramstad, O., Trulsen, K., Mon-
agating Deep-water Wave Trains”, Proceed- baliu, J., Bitner-Gregersen, E. and Onorato,
ings of Coastal Engineering, pp. 658-671. M., 2010, “Evolution of Weakly Nonlinear
Random Directional Waves: Laboratory Ex-
Stansberg, C.T., 1998, “Modeling and Analysis periments and Numerical Simulations”,
of Multidirectional Waves: Statistical Varia- Journal of Fluid Mechanics, 664, pp.313-
bility and Directional Resolution”, Proceed- 336.
ings of the Real Sea Conference, KRISO
(KR), pp. 22-39. Toffoli, A., Bitner-Gregersen, E., 2011, “Ex-
treme and Rogue Waves in Directional
Stansberg, C.T., 2001, “Random Waves in the Wave Fields”, The Open Ocean Engineering
Laboratory – What is Expected for the Ex- Journal, 4, 24-33.
tremes”, Proceeding of the Rogue Waves
conference 2000, IFREMER (FR), pp. 289- Toffoli, A., Ducrozet, G., Waseda, T., Onorato,
301. M., Abdolahpour, M., Nelli, F., 2019,
“Ocean Currents Trigger Rogue Waves”, In
Straub, L., Bowers, C., and Herbich, J., 2011,” Proceedings of the Twenty-ninth Interna-
Laboratory Tests of Permeable Wave Ab- tional Ocean and Polar Engineering Confer-
sorbers”, Coastal Engineering Proceedings, ence, Honolulu, Hawaii, USA.
1(6), 44. doi:
https://doi.org/10.9753/icce.v6.44. Torsethaugen, K., 1993, “A Two Peak Wave
Spectrum Model”, Proceedings of the 12th
Suh, K.D., Park, W.S. and Park, B.S., 2001, International Conference on Ocean, Off-
“Separation of Incident and Reflected shore and Arctic Engineering, Vol. II, 175-
Waves in Wave-current Flumes”, Coastal 180.
Engineering, 43, pp. 149-159.Tayfun, M. A.,
1980, “Narrow-band Nonlinear Sea Waves”, Tulin, M. P. & Waseda, T., 1999, “Laboratory
Journal of Geophysical Research-Oceans, Observation of Wave Group Evolution, in-
Vol. 85, pp. 1548-1552. cluding Breaking Effects”, J. Fluid Mechan-
ics, Vol. 378, pp. 197-232.
Sulisz, W., Hudspeth, R., 1993, “Complete sec-
ond-order solution for water waves gener-
ated in wave flumes. Journal of Fluids and
Structures, 7, 253–268.
ITTC – Recommended 7.5-02
-07-01.2
Procedures and Guidelines Page 27 of 27
Ursell, F., Dean, R. G., Yu, Y. S., ‘Forced small- Xiao, W., Liu, Y., Wu, G. and Yue, D.K., 2013,
amplitude water waves: a comparison of the- “Rogue Wave Occurrence and Dynamics by
ory and experiment.’, Journal of Fluid Me- Direct Simulations of Nonlinear Wave-field
chanics, 1959. Evolution”, Journal of Fluid Mechanics,
720, pp.357-392.
USACE, 2002, “Coastal Engineering Manual”,
U. S. Army Corps of Engineers. Xu, D., Hou, W., Zhao, M. and Wu, J., 1993
,”The statistical simulation of wave groups”,
Van den Berg, J., 2011, “Comparison of Meth- Applied Ocean Research, Vol. 15, pp. 217-
ods for Short Crested Wave Analysis”, Pro- 226.
ceedings of the ASME 30th International
Conference on Ocean, Offshore and Arctic Yang Z., and He H., 1984, “The Wind Wave
Engineering, OMAE2011-49443, Rotter- Spectrum in Lianyun Harbor”, Chinese
dam, The Netherlands. Ocean Engineering, (3):46-58. (In Chinese
)Yang Z., He H., and Gao Z., 1986, “Sta-
Van Essen, S., Pauw, W. Schmittner, C., 2014 tistical Analysis and Spectral Analysis of
“Improvemnt of Shallow-water Basin Wave Wave in Chiwan Port”, Port & Waterway
Generation Using ‘Anti-waves’”, Proceed-
Engineering, (12): 5-10. (In Chinese)
ings of the ASME 33rd International Confer-
ence on Ocean, Offshore and Arctic Engi- Yim, S.C., Osborne, A.R. and Mohtat, A., 2017,
neering, OMAE2014-23449, San Francisco, “Nonlinear Ocean Wave Models and Labor-
USA. atory Simulation of High Seastate and
Rogues Waves,” Proceedings of the ASME
Van Essen, S., Pauw, W., van den Berg, J., 2016,
36th International Conference on Ocean,
“How to Deal with Basin Modes When Gen-
Offshore and Arctic Engineering,
erating Irregular Waves on Shallow water”,
OMAE2017-62706, Trondheim, Norway.
Proceedings of the ASME 35th International
Conference on Ocean, Offshore and Arctic Zakharov, V.E., 1968, “Stability of Periodic
Engineering, OMAE2016-54134, 2016, Waves of Finite Amplitude on the Surface of
Busan, South Korea. a Deep Fluid”, Journal of Applied Mechan-
ics and Technical Physics, 9(2), pp.190-194.
West, B.J., Brueckner, K.A., Janda, R.S.,
Milder, D.M. and Milton, R.L., 1987, “A Zhang, H., Schäffer, H.A., 2007, “Approximate
New numerical Method for Surface Hydro- Stream Function wavemaker theory for
dynamics”, Journal of Geophysical Re- highly non-linear waves in wave flumes,”
search: Oceans, 92(C11), pp.11803-11824. Ocean Engineering, 34, 1290-1302.
Wu X., and Teng X., 1995, “Application of Four
Typical Directional Spectrum Models to Bo-
hai Sea Area”, Acta Oceanologica Sinica,
17(3): 18-23.Wu X., Liao K., Lin Y., and Lu
F., 1998, “Wave Spectrum of Quanzhou
Bay”, Journal of Oceanography in Taiwan
Strai, 17(3): 282-288. (In Chinese)