Volume II
Volume II
Volume II
TOWINGTANKCONFERENCE
Copenhagen,DENMARK
p
g ,
Aug.31 Sep.5,2014
PROCEEDINGS VOLUMEII
Table of Contents
Proceedings of the 27th ITTC
Volume I
Preface
Organising committee
Supporting organisations
Report of the Executive Committee
10
14
60
128
195
263
332
414
Appendix 1
439
Appendix 2
444
Appendix 3
458
Appendix 4
475
Appendix 5
Member Organisations
496
Appendix 6
513
Volume II
Report of the Specialist Committee on CFD in Marine Hydrodynamics
522
568
585
639
680
726
INTRODUCTION
1.1
Membership
Chairman:
Prof. Takanori Hino
Yokohama National University, JAPAN
1.2
Secretary:
Prof. Pablo Carrica
The University of Iowa, USA
Meetings
Members:
Dr. Riccardo Broglia
Istituto Nazionale per Studi ed Esperienze di
Architettura Navale (INSEAN), ITALY
1.3
Tasks
1. Review,
from
an
interdisciplinary
perspective, the current status of CFD in areas
of importance to ITTC. Include resistance,
propulsion, propulsors, manoeuvring, steep and
breaking wave simulation, seakeeping, ocean
engineering and steady and unsteady flow field
prediction at model and full scale.
2. Review developments and identify need for
research in steady and unsteady CFD at full
522
2.
523
524
3.
525
workhorse
in
applications.
marine
hydrodynamics
526
4.
NEW DIRECTIONS
4.1
527
528
4.2
529
5.
WAKE SCALING
5.1
Introduction
5.2
530
10
531
11
5.3
draft, when using the RANS-predicted fullscale wake field in the cavitation test.
5.4
532
12
5.5
5.6
Starke & Bosschers (2012) discuss
RANS/BEM coupling errors in the prediction
of the effective wake especially at the propeller
root. They suggest calculating the induced
velocities only from the dipoles, not from the
sources.
Conclusions
533
13
6.
6.1
Introduction
6.2
Verification of calculations
534
14
6.3
535
15
6.4
Validation
6.5
7.
TRENDS IN NAVAL
ARCHITECTURE APPLICATIONS
Conclusions
536
16
Propulsion,
Propellers,
Seakeeping,
Manoeuvring, and Ocean Engineering.
7.1
Resistance
537
17
538
18
539
19
7.2
Propulsion
540
20
541
21
542
22
7.3
Propellers
543
23
24
545
25
546
26
7.4
Seakeeping
547
27
Carrica
et
al.
(2011)
presented
computations of the KCS model in head waves,
i.e. one of the cases proposed at the
Gothenburg 2010 CFD Workshop. Results are
compared with other submission to the
workshop and with experimental data; in
548
28
549
29
Manoeuvring
550
30
551
31
552
32
Manoeuvring in waves
RANS simulations of manoeuvring in
waves have also become possible. In Carrica et
553
33
554
34
555
35
7.6
Ocean Engineering
There is literature that deals with freesurface flows around offshore structures.
Problems of interest vary from green water
shipping to tidal stream energy conversion
systems.
FPSO's are generally operated in a specific
region and positioned to meet mostly head or
bow waves in order to reduce roll motions. In
Lim et al. (2012), experimental results for three
different FPSO bow shapes in regular head
waves were analysed and compared to each
other. Also CFD computations were carried out
as a sample validation case for the database
built for CFD code validation.
556
36
557
37
8.
CONCLUSIONS
9.
RECOMMENDATIONS TO THE
FULL CONFERENCE
558
38
10.
REFERENCES
559
39
560
40
561
41
Huang, B., Wang G., Yu Z., Shi S. Detachededdy Simulation for Time-dependent
Turbulent Cavitating Flows, J. Mech. Eng.,
Vol. 25, pp. 484-490, 2012.
562
42
563
43
564
44
Richardson,
L.F.,
The
approximate
arithmetical solution by finite differences of
physical problems including differential
equations, with an application to the
stresses in a masonry dam, Phil. Trans. R.
Soc. London, S. A 210, pp. 307357, 1911.
565
45
566
46
567
INTRODUCTION
1.1
Pusan National
February 2014
Membership
1.2
Korea,
SCOPE
3.
OVERVIEW
Meetings
The committee met four times:
University,
568
4.
4.1
Collaborate
with
the
Specialist
Committee on CFD to develop methods
for the validation of CFD codes using
detailed flow measurements.
STATE-OF-THE-ART REVIEW
569
570
571
0,2
0,1
0,1
-0,1
Vr(MARI
N)
0
0
200
400
-0,2
-0,1
Vr(CSSRC
)
-0,3
0,4
0,1
-0,2
200
400
-0,1
Vt(CSSRC)
-0,3
0,6
0,4
Vx(MARI
N)
0,2
Vx(CSSRC
)
0
200
200
400
Vr(MARI
N)
Vr(CSSRC
)
0,8
Vx(MARI
N)
0,6
0,4
Vx(CSSRC
)
0,2
0
0
400
Vr(MARI
N)
Vr(CSSRC
)
-0,2
-0,4
400
Vt(MARI
N)
200
-0,2
-0,3
0,2
200
400
ar/R=0.3
br/R=1.0
Figure 6. Comparison of SPIV data between MARIN and CSSRC (Wake of 320k VLCC).
4.2
572
573
4.3
574
575
5.1
5.
576
577
5.2
578
5.3
579
Component
error
estimation:
individual error sources are estimated
and propagated into an overall
uncertainty in the measurement.
6.
CONCLUSIONS
580
7.
RECOMMENDATIONS
7.1
7.2
analysis
581
8.
REFERENCES
582
583
Zhang, G., Lu, L., Xue, Q., and Wu, Y., 2011,
Comparison and Analysis on Hydrofoil
Wake Turbulence by PIV and TR-PIV,
PIV11, Tsukuba, Japan.
Zhang, G. and Zhang, N., 2013, The Test
Report of Fluctuation Pressure and Flow
Field of Cavity, CSSRC Technical
Report.
584
1.
INTRODUCTION
1.1
1.2
585
9. Review the technologies (hydrodynamic issues) for enhancement of the powering performance, such as speed reduction, energy
saving devices, hull form and propeller design, etc.
1.3
General Remarks
6. Take into account minimum power requirements for safe and effective manoeuvring with respect to the EEDI formula (sea
margin)
7. Describe the type of data (and the quality of
that data) that should be recorded during
full scale monitoring trials, including the issues of surface roughness.
8. Review the existing ITTC trial test procedures in this context. Review the existing
speed correction methods for Full Scale
Trial Measurements including ISO 15016,
and come up with recommendation if the
problems are identified, taking into account
the MARIN report as contained in document MEPC 62/5/5.
586
2.
2012 Guidelines on the method of calculation of the attained Energy Efficiency Design Index (EEDI) for new ships;
2012 Guidelines for the development of a
Ship Energy Efficiency Management Plan
(SEEMP);
2012 Guidelines on survey and certification
of the Energy Efficiency Design Index
(EEDI); and
Guidelines for calculation of reference lines
for use with the Energy Efficiency Design
Index (EEDI).
587
3.
The speed reduction coefficient fw was introduced in the 2012 Guidelines on the method
of calculation of EEDI for new ships, adopted
by MEPC.212(63). fw is a non-dimensional coefficient to compensate for the involuntary
speed loss in a representative sea condition of
wave height, wave frequency and wind speed.
In agreement with the AC as to the required action/deliverables on this task, it was
decided that the Committee should not develop
and deliver its version of fw prediction procedure, but rather just monitor and eventually cooperate with the Resistance, Propulsion and
Sea-keeping committees on this issue.
In this respect the Committee needs to admit that despite numerous attempts, liaison and
cooperation with the above committees was
poor. It seemed that the Sea-keeping committee
encountered difficulties in providing a sound
basis for the pertinent calculations.
4.
4.1
Regulatory Framework
The regulatory framework for GHG reduction from international shipping is based on
amendment of Annex VI of MARPOL Convention at IMO/MEPC62 and it started from 1st of
January, 2013 (Resolution MEPC.203(62)).
588
certification has been adopted, where it prescribes ITTC recommended procedures for
speed/ power trial part1 and part2 as preferable
standard (Resolution MEPC.234 (65)).
4.2
EEDI
At MEPC64, Oct. 2012, the Interim Guidelines for fw, non-dimensional coefficient for decrease in ship speed in a representative sea
condition,
has
been
approved
(MEPC.1/Circ.796). ITTC recommended procedure for speed/power trial analysis is approved as preferable standard (MEPC64/23).
For each ship type, the reduction rate is determined by ship size and phase of application.
The reduction rate is given in Table 2 where the
base of the reduction rate is determined by
Guidelines for reference lines.
589
Capacity is deadweight except that, for passenger ships and Ro-Ro passenger ships, Capacity is gross tonnage, and for container ships,
Capacity is 70 per cent of the deadweight,
Reduction of EEDI
after 1st,
after 1st,
Jan., 2013 Sep. 2015
X
X
**
X
X
X
X
X
X***
X
X
X
In the guidelines for EEDI survey and certification, EEDI calculation method based on
CFD may be accepted as equivalent to propeller open water test or used to complement the
tank tests conducted, such as evaluation of the
effect of energy saving device with approval of
verifier (Res. MEPC.214(63), 2012).
* not apply to ships of diesel-electric propulsion, turbine propulsion and hybrid propulsion system
** separate to Gas carrier and LNG carrier
*** only regulated for Passenger ship having
non-conventional propulsion
The CFD based methods for EEDI calculation will be reviewed when they are available.
Attained EEDI =
EEDI Numerator
EEDI Denominator
(1)
(2)
where
Vref is ship speed in a calm sea condition on
deep water,
590
The ITTC PSS Committee developed completely revised speed/power trials recommended procedures and these were proposed at
MEPC64. ITTC contribution to IMO/MEPC
has been acknowledged and the ITTC recommended procedures were denoted as preferable
standard for analysis of speed/power trials
(MEPC64/4/15, MEPC64/INF.6, MEPC65/
INF.7).
Size
20,000DWT
Bulk carrier
and above
10,00020,000DWT
10,000DWT
Gas carrier
Tanker
and above
2,00010,000DWT
20,000DWT
and above
4,00020,000DWT
15,000DWT
Container ship
and above
10,00015,000DWT
15,000DWT
and above
3,00015,000DWT
5,000DWT
Refrigerated cargo
carrier
Combination carrier
and above
3,0005,000DWT
20,000DWT
and above
4,00020,000DWT
LNG carrier***
Cruise passenger
ship*** having nonconventional propulsion
and above
1,0004,000GT
10,000DWT
and above
85,000GT
and above
25,00085,000GT
Phase 0
1, Jan.,
201331, Dec.,
2014
Phase 1
1, Jan.,
201531, Dec.,
2019
Phase 2
1, Jan.,
202031, Dec.,
2024
Phase 3
1, Jan.,
2025
and
onwards
10
20
30
n/a
0-10*
0-20*
0-30*
10
20
30
n/a
0-10*
0-20*
0-30*
10
20
30
n/a
0-10*
0-20*
0-30*
10
20
30
n/a
0-10*
0-20*
0-30*
10
15
30
n/a
0-10*
0-15*
0-30*
10
15
30
n/a
0-10*
0-15*
0-30*
10
20
30
n/a
0-10*
0-20*
0-30*
n/a
5**
15
30
n/a
5**
20
30
n/a
0-5* **
0-20*
0-30*
n/a
5**
20
30
n/a
0-5* **
0-20*
0-30*
n/a
10**
20
30
n/a
5**
20
30
n/a
0-5* **
0-20*
0-30*
IMO
prompted
ISO
to
update
ISO15016:2002 taking into account ITTC recommended procedures 2012. Currently ITTC is
cooperating with ISO to revise ISO15016
(MEPC66/4/4, MEPC66/INF.7).
4.3
EEOI/SEEMP
* Reduction factor to be linearly interpolated between the two values dependent upon
vessel size. The lower value of the reduction factor is to be applied to the smaller ship size.
** Phase 1 commences for those ships when the amendments to MARPOL Annex VI come
into effect
*** Reduction rate applies those ships constructed on or after [date of entry into force].
Note: n/a means that no required EEDI applies.
5.
Ship speed in calm water (Vref) can not be
directly measured in speed/power trial runs.
Therefore, corrections for wind effect, wave effect, current effect, displacement and shallow
water are performed. At MEPC62 consideration for these corrections in order to reduce the
vagueness and improve the robustness was requested from Norway (MEPC62/5/5).
PREDICTION OF POWER
INCREASE IN WAVES
The existing model testing procedures, related to the EEDI prediction and verification
process are reviewed in more detail in the next
section of the report and will not be discussed
here. The procedure for prediction of power increase in irregular waves from model tests is
developed by the Sea-keeping committee, and
591
6.2
6.
Above model test should be performed according to ITTC Recommended Procedure 7.502-02-01, 7.5-02-03-01.1 and 7.5-02-03-02.1.
Numerical calculations may be submitted to
justify derivation of speed power curves, where
only one parent hull form have been verified
with tank tests, in order to evaluate the effect of
additional hull features such as bulb variations,
fins and hydrodynamic energy saving devices.
ITTC, with AC chairman as its representative, and PSS committee contributed to Part III
- Verification of EEDI of the Industry Guidelines. The following sub-sections present the
outline of the Industry Guidelines, model-ship
correlation process and recommendations for
future work on this topic.
6.1
Verification Process
592
6.3
Quality System
6.5
2. Quality manual containing at least the information listed in the ITTC Sample quality
manual (2002 issue) Records of measuring
equipment calibration as described in Appendix 3 of Industry Guidelines.
It is recommended to use the graphical construction described in Figure 1 that can be described by the following general procedure, applied only to EEDI power reference point (75%
of MCR):
3. Standard model-ship extrapolation and correlation method (applied method and tests
description) .
6.4
Model-Ship Correlation
The principal steps of the Speed Trial prediction calculation, are given in ITTC Recommended Procedure 7.5-02-03-01.4 ITTC 1978
Performance Prediction Method (in its latest
reviewed version of 2011). The main issue of a
performance prediction is to get the loading of
the propeller correct and also to assume the
correct full scale wake. The right loading of the
propeller can be achieved by increasing the
593
Study for The Model-Ship Correlation Between Full and Ballast Condition. This study
was carried out based on the data from SRC,
Shipbuilding Research Centre of Japan. Number of data is 773 of all kind of ships.
312
461
It seems that further study is necessary, because the scatter of the model-ship correlation
is very large. All the model basins and ship
yards should gather their own data. Figure 2
may be the only data which shows the modelship correlation on different displacement. And
the number of the data is not small. Before the
594
6.6
Conclusion
7.
MINIMUM
POWER
REQUIREMENTS WITH RESPECT
TO EEDI
7.1
General
According to the above documents, the following definitions, applicability and assessment
procedures apply:
595
7.2
Discussion
a
b
0.0687 2924.4
0.0689 3253.0
See tankers above
596
the weather and sea conditions for the evaluation of maneuverability to milder ones, especially for the smaller vessels. IMO/MEPC also
modified slightly the minimum power lines for
both tankers and bulk carriers to improve fitting
with the available statistical data. Thus, a
member of the PSS Specialist Committee, updated the original study for four typical bulk
carriers of DWT 30000 (Handy), 57000 (SupraHandy-Max), 79000 (Kamsar-Max) and
176000 T (Cape-Size) and a VLCC 306000 T
to evaluate the most recent requirements. Lines
plans and sea trial data were used to derive the
performance of the aforementioned vessels in
adverse weather conditions. Actual operation
points were derived by matching the power requirements with the propeller characteristics
and the main engine operational diagram.
Long-crested head sea waves (worst case for
added resistance) and the worst peak wave period TP within the range specified by MEPC
232(65) were assumed.
On the basis of these results, useful conclusions with respect to the installed power margin
that is necessary to ensure safe operation of the
ships in the prescribed sea and weather conditions, are drawn, as follows:
All five ships studied, very easily satisfy
level 1 requirement, while some of them
satisfy only marginally the requirements of
level 2 simplified assessment. This constitutes a major failure of rationalism, dictating that level 1 should be the strictest one.
Since either level is sufficient to comply
with the requirement, it follows that the
simplest level 1 should not also be the easiest to fulfill.
The submerged lateral area of the hull, corrected for breadth effect, estimated using
the formulation proposed by IACS and incorporated in MSC-MEPC.2/Circ.11 and
MEPC 232 (65) is 55-70% higher than the
actual one, both at the full load and the
heavy ballast conditions and in the case of
the VLCC more than 80% higher.
The differences between the estimated
power requirements on the basis of MSCMEPC.2 / Circ.11 and the calculated results
are higher in the heavy ballast conditions
than in the full load ones. This can be attributed to the fact that the approximate relations in the supporting document of IACS
(IMO MEPC 64/INF.7, June, 2, 2012) to
the document IMO MEPC 64/4/13 have
been derived mainly on the basis of the full
597
8.
SPEED/POWER TRIAL
FUNDAMENTALS
8.1
Background
598
The IMO brief asked that a transparent, unambiguous and practical method had to be delivered which would be acceptable for all
stakeholders and that could be used for both
contractual agreements between yard and
owner as well as for the assessment of the IMO
EEDI for any new-built ship worldwide. At the
same time, the results of the speed/power trials
should be completely documented and traceable for the EEDI Verifier representing the flag
state of the vessel. This task was conducted by
the ITTC Committee for the Performance of
Ships in Service with the assistance of the
STA-Group which has been working in this
field since 2004.
8.2
In December 2011, the PSS Committee invited STA-Group to co-operate on the new
ITTC Guidelines for Speed/Power trials to be
submitted to IMO MEPC for EEDI verification.
STA-Group accepted the invitation and provided access for ITTC to STA-Group data and
method.
8.3
Approach
History
Speed/power trials are conducted to establish the performance of the vessel at design
draught and trim under stipulated weather conditions, usually deep water, no wind and no
waves. As the conditions encountered during
the trials often deviate from the contract conditions, corrections are applied during the analysis and reporting of the trial results. In the past,
institutes such as BSRA, NSMB, SNAME and
ITTC published methods for conducting and
analysing speed/power trials. Shipyards randomly selected and developed their own yard
standard from these methods. In 2002, the International Standard Organisation published
ISO 15016, which included a cumbersome
analysis method based on a wide choice of outdated correction methods and empirical data.
The analysis method is based on the old manual
599
The required number of double runs at various power settings was specified:
600
Corrected
Measured
VWT2
For sister ships, the programme can be reduced to one double run at contract power, at
EEDI power and at one other power setting between 65% and 100% of MCR. In adverse environmental conditions, additional double runs
are required. The measurements and recording
of all required signals during speed runs with a
minimum duration of 10 minutes have been
specified in detail in these Guidelines.
VWR1
VG1
V'WR1
VG2
V'WT
VWR2
VWT1
VG1
V'WT
VG2
V'WR2
Figure 5 Averaging of measured wind vectors over two counter runs to derive the true
wind vector
Wind drag coefficients for ships have been
published by many authors in the past; however
modern vessels are much larger and have a different geometry than ships used in well-known
wind resistance publications. Therefore, it is
important to use recent ship type and size specific data derived from proper wind tunnel
measurements or validated computational tools
such as LES-RANS CFD. For container ships,
it is crucial to distinguish the wind drag in ballast condition without containers on deck but
while taking into account the lashing bridges
(which are exposed to wind during trials) and
the design draught case where the vessel is
loaded with containers. Remarkably the wind
resistance coefficient of the loaded vessel is
normally smaller as the full container pack provides a better flow shape than the wheelhouse
and lashing bridges!
601
motion induced
added resistance
In all cases, it is possible to use the drag coefficients derived by means of qualified wind
tunnel tests or validated CFD analysis for the
specific ship geometry.
ra w
added resistance
due to reflection
More elaborate presentation of the recommended methods for wind drag correction follows in the next chapter.
1/2
1/4
1/8
/ Lpp
Figure 6 Added resistance in waves as function of wave length over ship length.
Wave correction. Even within the trial limits for wave height, the added resistance due to
waves can be a substantial part of the required
shaft power. The added resistance in waves increases quadratically with wave height and thus
even in low sea states the wave correction
method should provide an accurate prediction
of the added resistance for the specific ship and
the actual encountered wind driven sea and
swell conditions. At the same time, the method
should be practical requiring limited input; today, many yards refuse to deliver the body plan
to the shipowner and the encountered wave
spectrum is not normally measured.
Raw =
B
1
gH s2 B
Lb
16
Where:
B
= Beam of the vessel on the waterline[m]
Lb
= Distance of the bow to 95% of maximum beam on the waterline [m]
Hs
= Significant wave height [m]
The above expression is particularly practical for speed/power trials as only the ships
beam, the length of the bow section and the
significant wave height are required as input.
No other ship particulars such as parametric
coefficients or bluntness factors nor ship speed
602
or wave spectrum are required. It is simply assumed that the asymptotic short wave value of
the transfer function extends over the complete
range of wave frequencies and thus that the
vessel is not heaving and pitching, which can
be easily checked during trials.
250
STAWAVE1
STAWAVE2
200
150
100
50
0
Tanker loaded high speed
180
STAWAVE1
160
STAWAVE2
140
120
100
80
60
40
20
0
Container loaded
high speed
Container loaded
low speed
Container ballast
high speed
Container ballast
low speed
603
model tests in regular head waves for the specific ship geometry.
Prior to adopting the three wave correction
methods, ITTC PSS subjected these methods to
extensive correlation with model test results
made available for this purpose by HSVA,
MARIN, NMRI and SSPA. Some comparative
results are presented in the next Chapter 9.
Alternative to the use of the above prediction methods, the transfer function of added resistance can also be derived from seakeeping
model tests for the specific ship geometry and
loading conditions. This transfer function can
then be applied to the wave spectrum measured
during the trials.
604
The method has been validated by a comparison with model tests which were conducted
by SSPA specifically for ITTC PSS.
The documented procedure has been included in the final version of the ITTC recommended procedure 7.5-04-01-01.2. The
speed/power trial analysis thus requires the
load variation tests as a standard procedure to
be included in calm water ship power model
tests.
Corrections for temperature & density. The
usual corrections of power for temperature and
deviations are incorporated in the new procedures.
Corrections for Water Depth. In the new
ITTC speed/power trial procedure the speed
corrections for shallow water according to the
method published by Lackenby, 1963 has been
implemented. With the use of CFD analysis it
has been proven by Raven, 2013 that this
method strongly overpredicts the effect of shallow water. The reason for this is that the
method is based on systematic model tests in a
shallow water basin. The resistance of the
model is not only influenced by the water depth
but also by the horizontal restrictions of the
towing tank. Especially in shallow water the
horizontal restriction of the basin has a large
effect on the resistance. Work of the STAGroup is underway to develop and validate new
method by means of speed-power trials at a
range of water depths. This is considered to become one of the important improvements of the
new procedures in the near future.
Effect of surface roughness. The added resistance due to (hull/propeller) surface roughness is not addressed in the procedure. It is required that the ship should go on sea trial with
clean hull and propeller. In case some kind of
surface fouling is documented, the hull needs to
605
Wind effect. Using the wind tunnel database of NMRI for 54 ships, the following six
methods are compared in view of practical use.
The database contains contemporary ship
shapes.
Selected methods are (1) Fujiwara 2005
(Fujiwara et al., 2005), (2) Fujiwara 1998 (Fujiwara et al., 1998), (3) Isherwood (Isherwood,
1973), (4) Yamano (Yamano and Saito, 1997),
(5) Yoneta (Yoneta et al., 1992) and (6) STA
DataSet (Sea Trial Analysis Joint Industry Project, 2006). The methods (1) to (5) are regression formulae and method (6) is dataset for coefficient of wind resistance.
9.
9.1
Introduction
9.2
SE EST =
1 1 ns n
C AAij C AAij
nS n j =1i =1
(3)
Correction Methods
606
SEEST
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
SEEST
0.4
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
0.3
0.4
0.2
0.3
0.1
0.2
0.0
CAA
0.1
0.0
CAA
Figure 9 Averaged standard errors of longitudinal wind force coefficient (54 ships).
Car Carrier
a) Tanker: 16 ships
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
0.4
SEEST
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
0.3
0.2
0.4
0.1
0.3
0.0
0.2
CAA
Cruise Ferry
0.1
0.0
Tanker
CAA
SEEST
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
0.4
SEEST
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
0.3
0.2
0.4
0.1
0.3
0.0
0.2
CAA
0.1
0.0
CAA
Container Ship
607
LNG Carrier
0.5
Fujiwara2005
Yamano
Fujiwara1998
Yoneta
Isherwood
STA
0.4
0.3
0.2
0.1
0.0
CAA
Figure 10 Averaged standard errors of longitudinal wind force coefficient for each ship
type.
From the validation, it was found that Fujiwara 2005 method gives the best estimation.
608
STAWAVE2
Exp.
Faltinsen
STAWAVE1
NMRI
F-T
KAW
4.0
3.5
3.5
3.0
3.0
2.5
2.5
2.0
2.0
1.5
1.5
1.0
1.0
0.5
0.5
0.0
0.0
0.5
1.0
1.5
2.0
STAWAVE1
NMRI
F-T
KAW
4.0
0.0
2.5
0.0
0.5
STAWAVE1
NMRI
F-T
KAW
STAWAVE1
NMRI
F-T
4.0
STAWAVE2
Exp.
Faltinsen
2.5
STAWAVE2
Exp.
Faltinsen
3.5
3.0
3.5
2.5
3.0
2.0
2.5
1.5
2.0
1.0
1.5
0.5
1.0
0.0
0.5
0.0
2.0
1.5
/LPP
/LPP
4.0
1.0
STAWAVE2
Exp.
Faltinsen
0.0
0.5
1.0
1.5
2.0
0.0
0.5
1.0
1.5
2.0
/LPP
2.5
/LPP
STAWAVE1
NMRI
F-T
KAW
4.0
STAWAVE2
Exp.
Faltinsen
3.5
3.0
STAWAVE1
NMRI
F-T
KAW
4.0
2.5
STAWAVE2
Exp.
Faltinsen
2.0
3.5
1.5
3.0
1.0
2.5
0.5
2.0
0.0
1.5
0.0
0.5
1.0
1.5
2.0
/LPP
1.0
0.5
0.0
0.0
0.5
1.0
1.5
2.0
2.5
/LPP
609
STAWAVE1
800
700
600
500
400
300
200
100
0
100
200
300
400
Figure 17 Wave correction methods compared with NMRI model test data; note NMRI
model test results were used for the NMRI prediction method.
800
STAWAVE2
700
600
500
400
300
200
100
0
100
200
300
400
Figure 18 STAwave-1 and STAwave-2 correlated with model test results for various ship
types, loading conditions and speeds in irregular head waves
Comparison results for the NMRI method
combined with empirical relation for the parameter estimation applied to oblique/beam
waves has been provided by NMRI. Correlation
diagram is shown in Figure 19.
610
(1) STAWAVE1
Under the condition that the pitching and heaving are small/missing and head waves, the simplified estimation method can be used.
(2) STAWAVE2
In case only ship dimensions are available and
head waves, empirical correction method with
frequency response function can be used.
(3) NMRI method
611
Effect of Water Temperature and Salt Content. The effects of water temperature and salt
content are corrected considering difference of
water density and frictional resistance.
Tidal period is
about 12 hours
2 double runs
take 8 hours
Figure 20 Relation on tidal velocity and trials for large low speed ship as VLCC.
Load Correction. The load variation model
test has been selected to account for the influence of propeller loading on the propulsive efficiency. This method is also more transparent
and straight forward comparing to the previous
method (KT/J2 approach).
9.3
612
With the acceptance of these new procedures, the ITTC and IMO have established a
transparent, straightforward best practice and a
level playing field for the delivery of new
ships for all stakeholders.
10.
CONTEMPORARY TRENDS IN
OPTIMIZING SHIP PROPULSIVE
PERFORMANCE
Contemparary hull form design and optimisation is heavily supported by CFD (computational fluid dynamics) tools. Although the potential flow methods neglect all viscous effects
and cannot predict steep, breaking waves these
613
More and more RANS codes are used during the hull form design process giving much
more insight in local flow phenomena, boundary layer details or the complex wave formations emerging from partly submerged bulbous
bows in off-design conditions.
Van et al (2010) used a RANS code, coupled with a parametric hull form optimization
tool to develop hull forms that both minimize
resistance and improve the wake quality into
the propeller. In a sample case that examines
the afterbody of a tanker, the significant reduction in resistance was achieved with good
agreement with experimental data. However, it
is insufficient to judge the performance of afterbody on the basis of wake and resistance
only. Self-propulsion simulation and validation
are further needed.
Developments and applications of SBD
(Simulation-based Design to ship design was
also reported. Kim et al(2010) investigate on
the flexibility of use of some choices of design
variables including local and global ones in the
optimization of the KCS containership.
614
cepts can be developed prior to the modeltesting phase of a project. Yet, podded propeller has the disadvantage of shorter docking interval comparing with normal propeller which
should be taken into account.
615
616
617
Figure 32 PBCF
Combined devices. In 2012, an energy saving model test on a VLCC was carried out
in MARIC towing tank. The energy-saving devices consist of Simplified Compensative Nozzle(SCN), Rudder Bulb(RB) and Thrust
Fin(TF).
618
619
11.
Experimental methods to determine the effect of wind (wind force measurements in the
wind tunnel) and the effect of drift in side wind
conditions (force measurements using computerised planar motion techniques for a range of
specified drift angles, rudder angles and heeling
angles, if necessary) are state of the art and will
not be described in this context.
620
A. Mohseni, et al. (2012) presents the effects of waves, wind speed and direction, current speed and direction, and depth of water in
vessel voyage planning which is based on meteorology and satellite data and computer program based in the ISO/DIS 15016. The interpolation between satellite data, historical chart data and observed data can optimize voyage route
and cause reduction in sea passage time and
fuel oil consumption. Various analysis methods
for resistance increase due to ship motion,
wave diffraction, wind, steering, drifting, water
temperature, salt content, deviation of displacement, hull and propeller surface roughness
and shallow water effects are considered and
could be contained in computer program.
With the rapid progress of CFD technique,
numerical simulation has played more and
more important role in predicting wind resistance.
Zhu H., et al. (2009) performed the simulation of mean wind pressure distribution, form
factors and the wind loads of the platform
based on N-S equation.
621
Figure 42 Photo of free-running model experiment of a container ship (from behind the
ship).
Chuang (2013) performed a series of experiments on a model of 8000 DWT tanker in
a large towing tank and ocean basin. The model
was self-propelled and mainly running in moderate long wave conditions. Numerical simulation work was carried out in order to make
comparisons with experimental results. It was
concluded that time domain simulation is the
preferable method for the steering effect evelauation.
There have been large improvements in automated heading and steering control systems
technology from Resolution MEPC.213(63).
An integrated Navigation and Command System can achieve significant fuel savings by
simply reducing the distance sailed "off track".
The principle is simple; better course control
through less frequent and smaller corrections
will minimize losses due to rudder resistance.
Retrofitting of a more efficient autopilot to existing ships could be considered.
622
12.
623
624
white-box model for the purpose of ship routeing with respect to minimisation of emissions
was presented by Prpi-Ori & Faltinsen
(2012). Another example was published in
(Leifsson et al, 2008) where linear and nonlinear regression methods were implemented in
order to tune the general model to fit the characteristics of specific ship. It must be noted that
such models were developed for the purpose of
conceptual design. Therefore, use of white-box
models with their limitations and underlying
assumptions and uncertainties implicitly affects
their accuracies. These limitations were briefly
discussed by Petersen et al. (2012) revealing
that even large changes of ship performance
due to hull surface deterioration over one year
would not be detected. Similar conclusions
were found by Dinham-Peren & Dand (2010)
Some additional information with respect to
accuracy of white-box models were provided
by Leifsson et al. (2008).
white-box models,
black-box models and
grey-box models.
625
However, as the possible cause of the underperformance, the small range of speed change in
training data was suggested.
Use of so called grey-box models i.e. combination of semi-empirical (white) models and
black-box models were proposed in order to
overcome the limitations of pure black-box
models with reference to extrapolation. Both
parallel and serial combination of white and
black box models were tested by (Leifsson et
al, 2008) revealing similar performance. Similarly to Petersen et al. (2012) and Pedersen &
Larsen (2009) the ANN was used as the blackbox model. The white-box model was based on
Holtrop (calm water) and Isherwood (rough
water) methods. The developed model was
used for performance approximation of the 10k
DWT container ship sailing at design speed of
20 knots. The model was trained on the data set
registered during the quasi-static part of the
voyage (port approach and manoeuvres were
filtered) with use of MAREN energy management system installed onboard. Application of
grey-box model performed much better than
white-box model in terms of fuel consumption
prediction. However, it was noted that such result may be achieved due to largely simplified
white-box model used in comparison. Furthermore, it was stated that the same grey-box
model did not reveal satisfactory performance
in predicting the ship speed. Surprisingly, the
achieved results were only slightly better comparing to performance of pure white-box
model. Such results were difficult to explain.
626
its physical replicas built up with use of rapidprototyping based on two low-order models
(retaining 95% and 71% of the original surface
roughness). In order to trace the similarities and
differences of the replica surfaces, they were
scanned and the streamwise profiles of the
roughness amplitude were analysed with use of
probability density functions (PDFs). Both loworder representations preserved relatively well
flow characteristics outside roughness sublayer.
In case of the flow inside the boundary layer,
the differences are substantial. Only low-order
model, retaining 95% of the roughness details,
allowed to maintain the flow characteristics
close to the wall. It should be noted that in the
study of Mejia-Alvarez and Christensen (2010),
the roughness height was not the only scale
used to describe the surface. Beside the roughness height, the root-mean-square roughness,
skewness, flatness, streamwise and spanwise
surface gradients were presented. The need of
research within the field of correlation between
friction drag characteristics and surface texture
parameters was indicated by Flack and Schultz
(2010). They studied roughness parameters
presented in literature and the common surface
statistical parameters in order to identify hydraulically relevant roughness scales. Results
indicated that root-mean-square roughness
height (krms) and skewness of the surface elevation PDFs (sk) were the most effective parameters in terms of hydraulic performance of
the surface. A correlation between mentioned
parameters and commonly used sandgrain
roughness height (kS) was also provided:
(4)
627
13.
SURFACE ROUGHNESS
SIMULATION BY NUMERICAL
METHODS
13.1 Methods
The task to numerically simulate the effect
of surface roughness can be divided into two
parts: first the translation of a real roughness
condition as it appears in reality into simplified
parameters, like the equivalent sand roughness; and second to introduce the simplified
parameters into the numerical equations describing the near wall flow.
Sand-Grain Roughness in The Flow Equations. Thin boundary layer methods used together with a potential flow solver was until recently commonly used for ship flow simulations and is still relevant for many applications.
Surface roughness can be included in thin
boundary layer methods by introducing a velocity shift function. This applies a decrease in
the log-arithmetic layers mean velocity corresponding to the effect of roughness, see LeerAndersen & Larsson, 2003.
628
effects can be well simulated with these methods for a flat plate at Reynolds numbers corresponding to full scale ship applications. They
used the turbulence model SST k- which is
relevant and commonly applied in ship hydrodynamics.
13.2 Applications
Trial Speed-Power Prediction. The hull
surface of a newly built ship can be assumed to
be homogenous and well defined in terms of
size, texture and distribution (if we neglect the
fact that bio film growth can occur within a few
weeks in some locations). The step to translate
this kind of roughness to a single or dual parameter to be fed into the flow equations should
be possible, even though no examples thereof
have been found in the open literature. Castro
et al, (2011) demonstrate however how the surface roughness kS appearing in the roughness
allowance CF in the ITTC scaling procedure
can be translated to the surface roughness used
in wall functions and applied to the KRISO
container ship test case with good results.
13.3 Conclusions
Surface roughness is likely to affect not
only the skin friction on hull and propeller but
also the wake flow into the propeller.
Introducing
homogeneous
sand-grain
roughness into numerical methods for speed/
power prediction in trial condition seems to be
possible and several well documented methods
exist.
The possibility to study the effect of nonhomogenous roughness such as bio fouling in
operational condition is still limited. Progress
in this area would be helped by experiments
with consistent test techniques of a large number of realistic surface conditions (preferably
from the same laboratory). This could be used
to formulate models that bridge between real
roughness conditions and the simplified coefficients used in the numerical equations.
14.
CONCLUSIONS
The hull surface during operational condition is characterised by inhomogeneous roughness; barnacles, slime and corrosion that are
unevenly distributed over the hull and with
large variation in height, texture and density.
629
15.
REFERENCES
630
Andersen, P, Friesch, J., Kapple, J.J., Lundegaard, L. and Patience, G., Development of
a Marine Propeller with Nonplanar Lifting
Surfaces, Marine Technology, Vol. 42, No.
3, July 2005, pp. 144-158
Dang, Jie; Dong, Guoxiang; etal , An Exploratory Study On The Working Principles Of
Energy Saving Devices (ESDS) PIV,
CFD Investigations and ESD Design Guidelines. Proceedings of the 31st International
Conference on Ocean, Offshore and Arctic
Engineering OMAE2012 July 01-06, 2012,
Rio de Janeiro, Brazil
Bugalski T., Streckwall H., Szantyr J.A., Critical review of propeller performance saling
methods, based on model experiments and
numerical calculations, Polish Maritime Research, 4 (80), Vol. 20, 2013.
Candries M., Atlar M., Anderson C.D., Estimating the impact of new generation
antifoulings on ship performance: the presence of slime, Journal of Marine Engineering and Technology, 2003.
631
Eca, L., Hoekstra, M., 2011, "Numerical aspects of including wall roughness effects in
the SST k eddy-viscosity turbulence
model", Computers & Fluids, Vol 40 (2011)
pp299-314.
Garcia, V C, Thaler, G J. Comparison of Steering Control Algorithms for Optimized Autopilots. Proc. 7th Ship Control Systems
Symposium[C], 4:(1984):24-27.
Fujii, H. and Takahashi, T., 1975, "Experimental Study on the Resistance Increase of a
Ship in Regular Oblique Waves", Proc. of
14th ITTC, Vol.4, pp.351-360.
Fujiwara, T., Ueno, M. and Nimura, T., 1998,
"Estimation of Wind Forces and Moments
acting on Ships", Journal of SNAJ, Vol.183,
pp.77-90.
632
MA Y., Numerical Simulation and Experimental Study on the Aerodynamic Performance of the Sail. Wuhan University of
Science and Technology , 2009
Maruo, H., 1963, "Resistance in Waves", Research on Sea keeping Qualities of Ships in
Japan, SNAJ, Vol.8, 1963.
Mizokami, S., Kawakita, C., Kodan, Y., Takano, S., Higasa, S. and Shigenaga R., 2010,
"Development of Air Lubrication system
and verification by the full scale ship test",
Journal of the Japan Society of Naval Architects and Ocean Engineers, Vol.12, 6977.
633
Raven, H.C.; A computational study of shallow-water effects on ship viscous resistance, 29th Symposium on Naval Hydrodynamics, Gothenburg, August 2012.
Sasaki,
N.,
ZEUS and NOAH
Prord
jects of NMRI . Proceedings of the 3 International Symposium on Marine Propulsors, May 2013, Tasmania, Australia.
pp324-330
Raimo, Jaap van Heerd, Wave Damping Aftbody With Hybrid Podded Propulsors
SNAME2003
Sogihara, N., Tsujimoto, M., Ichinose, Y., Minami, Y., Sasaki, N. and Takagi, K., 2010,
"Performance Prediction of a Blunt Ship in
Oblique Waves", Journal of the Japan Society of Naval Architects and Ocean Engineers, Vol.12, pp.9-15.
634
Yu,
Long, etal
An Automatic Optimization Process for O
ptimal Ducted Propeller Design and its Application based on CFD Techniques. Proceedings of the 3rd International Symposium
on Marine Propulsors, May 2013, Tasmania, Australia. Pp314-319
van der Ploeg, A., Raven, H., 2010, CFDBased Optimization for Minimal Power and
Wake Field Quality Proc. of 11th International Symposium on Practical Design of
Ships and Other Floating Structures, Rio de
Janeiro, Brazil, pp.92-101.
van Manen, J.D., van Oossanen, P., 1988, Propulsion, Principles of Naval Architecture;
Volume II, SNAME.
Wei Jin-fang, et al. The Study On the Prediction Method Of the Ships Speed Loss In
the Wind Wave. The ninth National hydrodynamic Conference[C], 2009
635
Journal References
IACS PR No.38, 2013, "Procedure for calculation and verification of the Energy Efficiency Design Index (EEDI) ", IACS.
ITTC 7.5-04-01-01.1 : Recommended Procedures and Guidelines Speed and Power Trials, Part 1 Preparation and Conduct, 2012.
ITTC 7.5-04-01-01.2 : Recommended Procedures and Guidelines Speed and Power Trials, Part 2 Analysis of Speed/Power Trial
Data, 2012.
MEPC64/INF.22, 2012, "First version of industry guidelines on calculation and verification of the Energy Efficiency Design Index
(EEDI) ", Submitted by BIMCO, CESA,
IACS, ICS, INTERCARGO, INTERTANKO, ITTC, OCIMF and WSC.
MEPC65/WP.10, 2013, "Air pollution and energy efficiency -Report of the working
group on air pollution and energy efficiency", IMO.
MEPC66/4/4, 2014, "Work Progress on revision of ISO 15016:2002", Submitted by ISO
and ITTC.
MEPC66/INF.7, 2014, "Additional information
on revision of ISO 15016:2002", Submitted
by ISO and ITTC.
636
Resolution MEPC.233(65), 2013, "2013 guidelines for calculating reference lines for use
with the energy efficiency design index
(EEDI) for cruise passenger ships having
non-conventional propulsion", IMO.
Resolution MEPC.212(63), 2012, "2012 guidelines on the method of calculation of the attained energy efficiency design index
(EEDI) for new ships".
Nomenclature
ANN Artificial Neural Network
Resolution MEPC.213(63), 2012, "2012 guidelines for the development of a ship energy
efficiency management plan (SEEMP)".
Resolution MEPC.214(63), 2012, "2012 guidelines on survey and certification of the energy efficiency design index (EEDI)".
637
Schiffbau
HSVA
Hamburgische
Versuchsanstalt
638
OVERVIEW
1.1
Herbert Bretschneider
HSVA, Germany
Johan Bosschers (secretary)
MARIN, Netherlands
1.2
Theodore Farabee
NSWCCD, USA
Chiharu Kawakita
Mitsubishi Heavy Ind., Japan
Denghai Tang
CSSRC, China
639
INTRODUCTION
640
3.1
Noise source
Impact to Impact to
environment the ship
Propeller noise
non-cavitating
tonal components
BPFs
Low/ medium
Depend
on ship
Singing propeller
100 Hz 2 kHz
high
high
Propeller
non-cavitating
broadband
1 Hz 20 kHz
low
low
BPFs
high
high
high
high
low
high
Propeller cavitating
tonal
Propeller cavitating
broadband
Frequency range
10Hz - 20kHz
Propeller-hull
interaction
BPFs and
structural NF
Cavitation on
appendages
100 Hz 20 kHz
medium
medium
Wave breaking
100 Hz 10 kHz
low
low
Slamming
1 Hz 100 Hz
low
low
100 Hz 10 kHz
medium
medium
Main engines
1 Hz 500Hz
medium
high
Driving systems
10 Hz 1 kHz
low
medium
Auxiliary engines
and systems
10 Hz 2 kHz
low
medium
100 Hz 50 kHz
high
Medium
Active sonar
echo-sounder
10 Hz 30 kHz
low
low
Active sonar
navigation
10 Hz 100 kHz
low
low
1 Hz 100 Hz
high
low
Airguns
641
chant ships. The spectral levels of noise generated by the vessel at various speeds and propeller rotation rates are shown in Figure 3.1.
642
3.2
Figure 3.2 Typical noise levels for different types of ships (modified from McKenna et
al., 2012).
McKenna et al. (2012) present measured
source levels for several types of ship: (a) container ships and vehicle carriers, (b) bulk carriers and open hatch cargos, and (c) three types
of tankers. Figure 3.2 shows the 1/3 octave
band source levels with mean and standard
errors. Figure 3.3 shows the broadband (20 to
1,000 Hz) source level for these ships as a
function of ship speed. There is significant differences in both source level and spectral characteristics of underwater noise amongst the
ship types for which measurements were made.
643
The sound pressure level of a noncavitating propeller is less intense and of less
impact compared to a cavitating propeller. The
features of cavitating and non-cavitating propeller noise spectra are illustrated in Figure 3.4
(Frchou and Dugu et al., 2000).
f FS = f M
or :
644
nFS
nM
LFS
n
= LM FS
nM
GFS
n
= GM FS
nM
D
FS
DM
4
D
FS
DM
r
M
rS
7
r
M
rS
645
646
i ,ms
Re fs
=
Re ms
647
Figure
(HSVA).
3.6
Singing
model
The pressure fluctuations due to the turbulent boundary layer is a rather inefficient
(quadrupole) sound source when considered in
isolation, but it can become more efficient in
the presence of a rigid or especially a flexible
surface such as the hull plating of which the
vibrations generate sound (Blake 1986). The
radiation efficiency of the hull plates is strongly influenced by fluid loading and by the presence of ribs and stiffeners. Both spatial and
temporal characteristics of the turbulent boundary layer pressure fluctuations need to be taken
into account for the excitation of the hull vibrations. Unsteady surface pressure measurements
have been performed by Goody et al. (2007) on
the surface of a ship model hull in a towing
tank. The results compare well with an empirical model and, for low frequencies, with computational results using a Reynolds Averaged
Navier-Stokes Statistical Model. Similar measurements have been performed by Ciappi et al.
(2009) and De Jong et al. (2009). The scaling
parameters are strongly related to Reynolds
number and include the boundary layer displacement thickness and the wall shear stress.
In addition, hull conditions are critical.
propeller
648
model it is possible to establish a perfect correlation between model and full scale in term of
load bending moments and of the first bending
modes of the ship. A detailed overview of the
method and of the results so far achieved can
be found in Hirdaris et al. (2014). When local
response is considered, hydrodynamic loads
(pressure and acceleration) can be measured on
rigid models and the structural response calculated numerically or theoretically.
3.3
Generality
Machine vibrations can originate in the following ways (Urick, 1983): i) unbalanced rotating shafts, ii) repetitive discontinuities, e.g.
gear teeth, armature slots or turbine blades, iii)
reciprocating parts, e.g. combustion in engine
cylinders (piston slaps), iv) cavitation and turbulence in fluids flowing through pipes,
pumps, valves, condensers, and v) mechanical
friction as in bearings and journals
The phenomenon of slamming can be accurately tested at model scale if the model is
properly scaled to replicate global hydro-elastic
effects and tested at the correct Froude number.
It has been demonstrated that with this physical
649
Other diesel-powered ships employ medium speed, 4-stroke diesel engines, which connect to the propeller shaft via a reduction gear.
Typical speeds of these engines are 300 to
1,000 rpm. The engines can be rigidly or resiliently mounted.
The noise emission of medium speed engines can be separated in two bands. The lower
band covers the range between 6 Hz and approximately 150 Hz. The noise in this range is
generated by mass forces of the moving pistons, conrods and crankshafts and by gas forces
arising from the internal combustion process
and exhibits distinct frequencies which are integer or half integer multiples of the shaft rate
frequency. For the higher frequency band, engine noise is broader band, excited by the internal combustion process, and thump noise of
pistons, gear wheels, and valves.
650
tonal and broadband are lower for comparable power-to-weight ratios. Furthermore, the
tones produced by a gas turbine are much higher due to their higher rotation rate, which can
be as high as 3,600 rpm or 60 Hz. A comparison of representative diesel and gas turbine
vibration levels is provided in Figure 3.8.
Diesel Electric. The noise signature of diesel-electric ships typically contain energy contributions from the diesel generators and from
the electric propulsion motors in combination
with
the
frequency
converters
(synchroconverter or a cycloconverter). The
levels of electric propulsion motor noise, and
the frequencies at which they occur, vary by
ship and with propulsion shaft rpm. Noise
sources for electrical machines can be mechanical (angular and parallel shaft misalignment,
dynamically unbalance rotors, loose stator lamination, bearing), and electromagnetic (phase
unbalanced, slot opening, input current waveform distortion, magnetic saturation etc.). Even
large direct drive electric motors are quiet if
compared with reduction gears and piston engines.
Figure 3.8 Source vibration levels for Diesel and Gas Turbine (Fisher and Brown, 2005).
Auxiliary Machinery. Noise components
from rotating auxiliary machinery and other
shipboard equipment also contributes to a
ships overall noise signature, but usually at
lower levels than the propulsion systems. A
typical frequency range for vibration spectra
for auxiliary machinery is 1 Hz to 5 kHz.
For diesel-electric systems, the diesel generators operate at a constant rpm and therefore
their noise characteristics are not dependent on
ship speed. Moreover when used as a genset
they are usually elastically mounted. The same
consideration holds when a gas turbine is used
as a generator.
Sea water cooling pumps. Sea water cooling pumps are mainly of a centrifugal type and
the impellers cause tonals at impeller blade
passage rate and related harmonics. Source
levels in the pipes close to the pump reach up
to 180 dB (re: 1Pa) for non-cavitating condition and can be more that 200 dB in cases of
impeller cavitation.
Mitigation of blade tonals can be achieved
within the pump by increasing the tip clearance
651
of the impeller and accepting a reduced efficiency. Another measure is to introduce fluid
silencers up- and down-stream of the pump.
Bow and Stern Thrusters. Bow and stern
thrusters are mainly horizontal axis tunnel type
impeller systems and are strong noise sources.
Most of the noise from thrusters is caused by
cavitation on the impeller blades. The cavitation causes direct radiated noise and also structure-borne noise which propagates through the
hull structure and can radiate as underwater
noise. The spectrum of thruster noise is broadband with energy covering a very wide frequency range. Specialized thruster types such
as azimuths, pumpjets etc. have different noise
characteristics compared to conventional
thrusters. (Lloyds Register Consulting 2013).
3.3.2. Sonars
3.3.2.1.
652
3.3.2.2.
The operational purpose of the AMS dictates the sonars frequency range and source
strength. A majority of active military sonars
are used for anti-submarine warfare purposes
and thus operate in the low (~100 to 1kHz) to
medium (~1kHz to 8 kHz) frequency range so
that signal strength is not significantly impacted by acoustic absorption which increases with
frequency. Sonars operating in the high frequency range (~ 8 kHz and higher) are generally used as navigational aids or for mine hunting
where interest is in detecting the presence of
objects at shorter ranges and where higher (spatial) resolution is needed.
Fish finders operate similarly to echosounders except that the intent is for the acoustic pulses to reflect off fish instead of the ocean
bottom. They also operate at higher frequencies
to provide discrimination and at low source
levels as to not adversely disturb the fish that
are trying to be located. It is noted that the frequency of operation is potentially set at a frequency that provides maximum acoustic reflection for the fish of interest. Fish finders are
used both commercially and recreationally.
Searchlight sonars, which includes sidescan sonars, and acoustic cameras are examples
of high frequency sonar systems used for the
purpose of imaging underwater objects. These
sonars generally operate at lower source levels
to reduce issues with multiple reflections and at
quite high frequencies to provide high resolu-
653
tion capabilities. Acoustic Doppler current profilers have become common instruments for
high accuracy measurement of speed, either of
vehicles on which they are mounted or of currents passing over them. Water speed is measured based on the Doppler frequency shift of
pulses reflected back from particulates in the
water. To obtain highly accurate measurements
of speed, they typically operate at high frequencies.
NOISE REGULATION
4.1
Influence of
Environment
Noise
on
Marine
A final type of ship-board sonar system includes those used for underwater acoustic
communications. They typically operate in the
medium frequency range and have low to medium source strengths. This category includes
systems used for voice communication or as
underwater acoustic modems. Most often these
systems are used for communication between
surface vessels and submerged vehicles.
3.3.3. Airguns
Currently almost all marine seismic surveys
use arrays of airguns as a noise source for
seismic signals. An airgun is a twin piston steel
cylinder charged with high-pressure air (up to
200 bar). After triggering by an electric signal
the airgun suddenly releases the compressed air
to the lower outside pressure causing a transient high pressure peak like from explosives.
The peak pressure reaches values of about
230 dB (re: 1Pa at 1m), with a spectrum that
is of broadband type. Most airgun noise occurs
in the range below 1 kHz with increasing levels
at lower frequencies with a maximum typically
below 100 Hz.
4.2
Noise regulation
The problem of anthropogenic noise emissions in the sea has been assessed only in recent years. This problem has been analysed
mainly at a regional level, in particular for re-
654
ly or jointly as appropriate, all measures consistent with this Convention that are necessary
to prevent, reduce and control pollution of the
marine environment from any source, using for
this purpose the best practicable means at their
disposal and in accordance with their capabilities, and they shall endeavor to harmonize their
policies in this connection. Hence, by recognizing underwater sound as a pollutant by virtue of the substances or energy wording in
Article 1(1)(4), then Article 192 grants each
nation the authority to prevent, reduce and control such pollution, etc.
stricted areas where there is a higher concentration of species of marine mammals or fishes.
The national and international regulations reviewed to date often do not address underwater
noise quantitatively in the sense of specifying
acceptable underwater source levels but instead
restrict activities that are determined to harass
or harm marine animals.
4.2.1. International Framework
At an international level there are several
associations which deal with the protection of
marine mammals. In some of their regulations
or treaties they cover underwater sound. In the
following a short review of some of these regulations as they relate to underwater noise is
presented.
655
These non-mandatory Guidelines are intended to provide general advice about reduction of underwater noise and focus on the primary sources of underwater noise such as associated with propellers, hull form, onboard machinery, and operational aspects. Moreover, a
specific section addresses the use of numerical
tools for noise prediction indicating that CFD
can be used to predict the flow characteristics
around the hull and appendages, thus providing
the wake field in which the propeller operates
and propeller analysis methods, such as lifting
surface theory, or CFD, can be used for predicting cavitation. Finally, mention is made
that Statistical Energy Analysis (SEA) and Finite Element (FE) methods can be used to solve
the vibro-acoustic problem at high and low
frequency, respectively.
Other Organizations. Declarations regarding the impact of shipping noise have also been
made by many other international organizations, among them, the Convention on the Conservation of Migratory Species of Wild Animals (CMS), the International Union for the
Conservation of Nature (IUCN), the International Whaling Commission (IWC), the International Council for the Exploration of the Sea
(ICES), the International Fund for Animal
Welfare, the Whale and Dolphin Conservation
656
Starting in 2012, two multinational collaborative projects are partly funded by the 7th
Framework Programme of the European
Commission with the goal to develop tools to
investigate and mitigate the effects of underwater noise generated by shipping. These projects
are SONIC (www.sonic-project.eu) and AQUO
(www.aquo.eu).
The Agreement on the Conservation of Cetaceans of the Black Sea, Mediterranean Sea
and Contiguous Atlantic Area (ACCOBAMS).
The ACCOBAMS is a cooperative tool for the
conservation of marine biodiversity in the
Mediterranean and Black Seas. Its purpose is to
reduce threats to cetaceans in Mediterranean
and Black Sea waters and improve our knowledge of these animals. ACCOBAMS was concluded in the auspices of Convention on Migratory Species in 1996 and entered into force in
2001.
In 2010, under resolution 4.17, guidelines
to address the impact of anthropogenic noise
on cetaceans in the ACCOBAMS area were
adopted and a working group was established
that will focus on the mitigation of noise impact issues resulting from sonar, seismic surveys, coastal and construction works, and maritime traffic including commercial shipping.
657
In 2012 the ACCOBAMS and ASCOBAMS noise working groups joined. The
working group has produced among other
things a review of various international guidelines (Maglio, 2013).
The NEPA requires full disclosure of possible environmental impact, alternatives and
mitigation measures for any federal actions and
thus has direct impact on military activities in
the ocean.
Other US agencies that also have a regulatory or enforcement role related to anthropogenic sound include: the Marine Mammals
Commission (MMC); the Minerals and Management Service (MMS, under Department of
Commerce) and the US Navy.
Moreover, the US Coast Guard is one of the
five Armed Services of the US and enforces a
wide range of maritime safety, security and
environmental policies of the US. Issues related to habitability concerns due to ship-borne
noise levels would in part be handled by the
Coast Guard.
658
4.3
Hz-to-1 kHz and the higher covers from 1 kHzto-100 kHz and are given as:
ICES methodology.
The International
Council for the Exploration of the Sea (ICES)
derived noise limits for research vessels
(Mitson, 1995). The ICES proposal considers
the cod audiogram and selects the noise limit
by considering the lowest point of the curve at
30 dB over the threshold sensitivity. The curve
was interpreted as the limit over which behavioral effects (escaping) started to appear. The
selected point at 200 Hz represents the frequency of maximum sensibility for cod, the
corresponding level was set as the limit for
radiated underwater sound from a research vessel, free running at 11 knots, at a target distance
of about 20 m. The noise limit at 1 m was obtained applying the spherical dispersion law.
Hz
1 Hz
= 130 22 log10
kHz
1 kHz
The curves for the mentioned categories report maximum allowable noise levels (in dB re:
1Pa at 1m) versus frequency (1/3 octave resolution). In the case of the acoustic, fishery and
environmental categories two different curves
are given depending on the operational conditions of the ship. The curve relative to research
vessels corresponds to the ICES one except for
the format (third octave bands instead of narrowband (1 Hz) and the shape of the curve for
frequencies below 25 Hz, which contains less
659
restrictive limits.
5.1
of
Propeller
Predicting the low-frequency discrete spectrum of non-cavitating propeller noise has been
the subject of research for many years. Most of
the prediction methods are focused on obtaining the unsteady propeller forces or fluctuating
pressures of the flow field using CFD calculations. In turn, these pressures and forces are
used as the sources to acoustic methods of predicting the discrete tonal noise of a propeller.
Most methods for predicting marine propeller noise are for far field and as such various
higher-order terms are neglected in the numerical solutions. However, in the near field, the
higherorder terms play important roles. The
solutions of blade rate underwater noise induced by the unsteady force of a marine propeller in the time domain is described in the
paper by Kehr and Kao (2004). The method
can be used not only for computing far field
acoustic pressure, but also for near field pressures.
Seol, Suh and Lee (2002, 2004, 2005) present a numerical method to predict noncavitating tonal noise of an underwater propeller. The noise is predicted using a time-domain
acoustic analogy. The flow field is analyzed
with a potential-based panel method and the
Tonal noise of a propeller can also be assessed using the measured fluctuating forces of
660
661
5.2
662
of cavitation after the collapse of sheet cavitation. The size distribution of cavitation bubbles
is proposed to be represented by the beta distribution. The trailing thickness of sheet cavitation is treated as uniform along a radial direction of the propeller blade. The simulation result compares favourably to full scale measurement for the SEIUN MARU full scale propeller.
Salvatore et al. (2009) presented a numerical prediction method for non-cavitating and
cavitating noise of propellers operating in a
wake field. Propeller hydrodynamics is described by a BEM coupled with a nonlinear
unsteady sheet cavitation model. Hydroacoustic
models are based on a standard Bernoulli equation for incompressible flows and on the FW-H
equation with a Transpiration Velocity Model
to account for blade cavitation effects. While
fair agreement between results from the two
formulations is found for a non-cavitating propeller configuration, quantitative differences
are observed for cavitating flow noise predictions between the Bernoulli and FW-H models.
Numerical uncertainty in the evaluation of cavity pattern could have a strong impact on prediction of radiated noise levels.
Salvatore and Testa (2006) developed an integrated hydrodynamics/hydroacoustics approach for marine propeller sheet cavitation
noise. The hydrodynamics model is based on a
boundary element method (BEM) that is a potential flow formulation. A sheet cavitation
model using a surface tracking approach is
applied to estimate the transient cavity pattern
on the surface of propeller blades. Propeller
cavitation noise is studied through a general
hydroacoustics formulation based on the FW-H
equation. Hallander et al. (2012) presented
work of the EU project SILENV and compared results of underwater radiated noise of a
propeller obtained from different prediction
methods, including URANS/FW-H, SYSNOISE/BEM, lifting-surface/FW-H, and a
semi-empirical approach with sea trial data for
an LNG ship under cavitating conditions.
663
664
6.1
Full Scale
Method
Noise
Attention to hull and propeller surface conditions in support of acoustic testing varied
amongst organisations. Nominally half report
that hull conditions are checked, that propeller
conditions are checked and that the propeller
surface is polished prior to testing.
All organisations use commercially available omni-directional hydrophones that generally provide measurement capabilities to a frequency of 100 kHz or greater. Each organisation but one reports the type of hydrophone
calibration procedures used. Hydrophone arrays generally amounted to the use of 3 hydrophones with two organisations reporting using
arrays consisting of 10 or more hydrophones.
Nearly all organisations report making
noise measurements up to a frequency of approximately 50 kHz and as low as 10 Hz. All
but two report performing instrumentation
calibrations either and/or before testing. The
definition of acoustic source used during testing varied amongst the organisations and only
4 organisations reported using some type of
source localisation technique. A number of
different data reporting formats are used with
many organisations using multiple formats.
The formats ranged from one-third octave band
levels to narrow-band levels, each with or
without being range corrected. For the question of expected uncertainty for measurements, all 5 who provided an answer gave a
value of 3 dB or less. And, for the question on
confidence level on the quality of the test results, based on a scale of 1-to-10 (10 being
Measurement
665
Regarding whether measured levels are corrected for background noise and propagationrelated issues, 6 out of 10 reported measurements are corrected for background noise, 8 out
of 10 do not correct for environment-related
propagation losses (i.e., absorption), 8 report
that a 20log10(r/1m) adjustment for range correction is made. Only 3 report making an adjustment to correct for free-surface reflection.
For the question of expected uncertainty for
predicted noise source level the 5 who responded gave a value of 3 dB or less. And,
for the question of confidence level on the
quality of correction methods based on a scale
of 1-to-10, of the 8 who responded the highest
value was 9, the lowest 5, and the median being 7.
6.2
Model Scale
Methods
Noise
Measurement
The ship wake distributions used is presented in Figure 6.3. Almost half of the organisations are simulating both a model scale wake
and a full scale wake in their facility. The wake
field is measured by pitot tubes and/or LDV
while 5 respondents use PIV as well.
666
phones (frequency range, dimensions, operating pressure and temperature range) developed
for full scale measurements easily comply with
the requirements of a cavitation test facility.
0
6
11
Width [m]
3
2.5
2
1.5
1
0.5
0
0
10
12
Length [m]
Full scale
wake
calculated
using CFD
The propeller load condition is chosen according to the design point (11) or the model
test result (8) corrected to full scale (16). The
test condition is mostly adjusted by the thrust
coefficient KT (14) and rarely by the advance
ratio J (3) or torque coefficient KQ (1). The
667
90
(%) of respondents
80
70
60
50
40
30
20
10
0
Free surface
668
Reflections
due to wall
Background Reverberation
noise
Others
GUIDELINES
FOR
MEASUREMENTS
7.1
NOISE
Based on responses to the full scale measurements questionnaire and a review of existing full scale guidelines, both established and
in development, it is recommended that the
ITTC guideline for underwater noise measurement of full scale ships follow the standard
ISO/PAS 17208-1:2012(E). This standard was
drafted by international experts in the field of
underwater ship noise measurements and provides general requirements for measurements
in deep water. It is noted that guidelines for
measurement of ship underwater noise in shallow water are currently being pursued by numerous international organizations but such
standards will be slower to finalize due to the
much more complex nature of measurements in
shallow water where noise propagation issues
related to bottom affects must be accounted for.
7.2
The guidelines for model scale measurement of propeller cavitation noise are included
in the ITTC Quality Manual as guideline no.
7.5-02-01-05.
Many aspects of cavitation noise measurements are related to other procedures such as:
Procedure 7.5-02-03-03.1 on Model-scale
cavitation tests. Part of this procedure is
discussed by the Specialist Committee on
Wake Fields of the 25th ITTC, 2008 and
the Specialist Committee on Scaling of
Wake Field of the 26th ITTC, 2011.
Procedure 7.5-02-03-03.3 on Cavitation
induced pressure fluctuations, model scale
experiments. This procedure is discussed by
the Specialist Committee on Cavitation Induced Pressures of the 23rd ITTC, 2002.
669
Hydrophone arrays enable noise measurements with high directivity to scan the model to
identify local noise source regions and should
be used if permitted by facility capabilities and
testing budget. Examples of the setup of an
array are provided by Abott et al. (1993),
Chang and Dowling (2009), Park et al.
(2009b), and Lee et al. (2012).
For every test condition the background
noise of the facility has to be determined to
check the quality of the acquired noise data and
to correct the cavitation noise data if the difference between the two is less than 10 dB. Two
procedures to measure background noise can
be applied: Replacement of the propeller by a
dummy boss, or increase of tunnel pressure to
suppress propeller cavitation. Examples of
background noise sources in a cavitation tunnel
are pump cavitation, non-cavitating turning
670
fs
fs fs
=
ms ms ms
For the noise levels, different scaling formula have been derived which have been presented by the Cavitation Committee of the 18th
ITTC (1987) as:
fs ms
fs ms = 20 log10
ms
fs
fs fs fs
fs 2
ms
ms ms
ms
This formula is valid for noise in proportional band widths. For the case of constant
band width, frequency scaling should be included as well and the powers of cavitation
number and shaft rotation rate change. Assuming spherical spreading loss, we find for the
above, that x= 1 is generally applied. For the
other parameters in this equation, two different
approaches for determining their values have
been reported. One approach is based on linear
acoustics (e.g. Strasberg 1977) while the other
approach assumes constant acoustic efficiency
which is considered to be more valid for higher
frequencies (De Bruijn and Ten Wolde 1974,
Levkovskii 1980). The corresponding parameters are given in Table 7.1. As the cavitation
number is usually kept constant, it is only the
scaling of the tip speed that is relevant.
671
672
10
REFERENCES
RECOMMENDATIONS
The 27th Specialist Committee on Hydrodynamic Noise recommends adopting the following guidelines:
7.5-02-01-05 : Model Scale Noise Measurements
7.5-04-04-01: Underwater Noise from
Ships, Full Scale Measurements.
Ainslie, M.A., 2011, Potential causes of increasing low frequency ocean noise levels,
161st Meeting of the Acoustical Society of
America, Vol. 12.
Andrew, R.K., Howe, B.M., Mercer, J.A.,
Dzieciuch, M.A., 2002, Ocean ambient
sound: comparing the 1960s with the 1990s
for a receiver off the California coast.
Acoustics Research Letters Online, Vol.
3(2), pp 65-70.
673
Baiter, H.J., 1989, On Cavitation Noise Scaling with the Implication of Dissimilarity in
Cavitation Inception, ASME International
Symposium on Cavitation Noise and Erosion in Fluid Systems, San Francisco, USA.
Bark, G., 1985, Prediction of propeller cavitation noise from model tests and its comparison with full scale data, Journal of Fluids
Engineering, Vol. 107, pp 112-119.
Chang, N.A. and Dowling, D.R. (2009), Raybased Acoustic Localization of Cavitation
in a Highly Reverberant Environment,
Journal of the Acoustical Society of America, Vol. 125(5), pp 3088-3100 .
Bosschers, J., 2010, On the influence of viscous effects on 2-D cavitating vortices, 9th
International Conference on Hydrodynamics, Shanghai, China.
Bosschers, J., Lafeber, F.H., de Boer, J., Bosman, R., Bouvy, A., 2013, Underwater
Radiated Noise Measurements with a Silent
Towing Carriage in the Depressurized
Wave Basin, Symposium on Advanced
Measurement
Techniques,
AMT13,
Gdansk, Poland.
674
675
IMO/MEPC 59/19, 2009 Noise from commercial shipping and its adverse impacts on
marine life report of correspondence
group
IMO/MEPC.1/Circ.833, 2014 Guidelines for
the reduction of underwater noise from
commercial shipping to address adverse
impacts on marine life .
IMO/MEPC 66/17, 2013 Noise from commercial shipping and its adverse impacts on
marine life outcome of DE 57
676
Kehr, Y-Z., Kao, J-H., 2004, Numerical prediction of the blade rate noise induced by
marine propellers Journal of Ship Research, Vol. 48 (1), pp 1-14.
Kellett, P., Turan, O., Incecik, A., 2013, A
study of numerical ship underwater noise
prediction, Ocean Engineering, Vol. 66, pp
113120.
Mitson, R.B., 1995, Underwater Noise of Research Vessels. Review and Recommendations, ICES Cooperative Research Report
N. 209.
Lvik, A., 1981, Scaling of propeller cavitation noise, Noise sources in ships,
Nordforsk, Stockholm, Sweden.
677
Rolland, R.M., Parks, S.E., Hunt, K.E., Castellote, M., Corkeron, P.J., Nowacek, D.P.,
Wasser, S.K., Kraus, S.D., 2012, Evidence
that ship noise increases stress in right
whales. Proceedings of the Royal Society
B, Vol. 279, pp 2363-2368.
Salvatore, F., Testa, C., Greco, L., 2009, Coupled HydrodynamicsHydroacoustics BEM
Modelling of Marine Propellers Operating
in a Wakefield, 1st International Symposium on Marine Propulsor, Trondheim,
Norway.
678
Seol, H., Suh, J-C., Lee S., 2005, Development of hybrid method for the prediction of
underwater propeller noise, Journal of
Sound and Vibration, Vol. 288, pp 345
360.
Seol, H., Park, C., Park, Y., Kim, G., 2013,
Measurement of Propeller Cavitation
Noise in the MOERI Large Cavitation Tunnel, Symposium on Advanced Measurement Techniques, Gdansk, Poland.
679
depend on
location
No
No
Yes
HHI
No
DMO
Hull
mounted
1721
on a trestle 1 m
above bottom
50
water depth
10
30
100
50
800
Satellite
information
200
30
over 200
5, 20
30
0320
max 150
100 m or ship
length
60
80
From vessel
60
50
14, 20, 60
Yes
TSI
Spain
Max. 150
14, 400
18- 22
Yes
TNO
Netherlands
Holland,
Norway
No
Hydrophone
hull above the
propeller
KRISO
Sea State
3
2-3
2
A9: Typical Closest Distance between Noise Source and Reference Hydrophone (cpa).
Beaufort
Mitsui Lab.
A: SITE INFORMATION
Korea
Eckernfrder
Bucht
No
MHI
Japan
Haringvliet,
Heggernes
Yes
CETENA
Italy
Aschau
Yes
WTD 71
Germany
Others
Others
Location
Name
Others
DCNS
ORGANISATION
France
COUNTRY
4*
N/A
Permanent bottom
anchor, depth tide
compensated
30300
Yes
NSWC/CD
USA
DCNS
ORGANISATION
WTD 71
Germany
CETENA
Italy
C2:Hydrophone Operation
Omnidirectional (O),
flush mounted (F)
no
before meas.
by the ship
owner
yes
No
yes
Yes
0.022
0.093
0.0095,0.022
0.05,0.093
4dB
at 100kHz
270; +/2dB
at 15 kHz
0.024
0.12
0.04
0.35
2dB
at100kHz,
Omni
-214
-144
Omni; +/-1.5dB
at 50kHz
-170
-40
+80
40atm
B&K 8103,
8105
100atm
B&K 8103,
8100
200k
B&K 8105
C: HYDROPHONE INFORMATION
yes
diver's
observation
diver's
observation
10
90k
KRISO
-2
+ 50
4
100k
5
120k
no
No
Mitsui Lab.
no
No
HHI
0.021
0,12
-205
-30
+80
100atm
0.1
10k
B&K 8104
Korea
MHI
Japan
-2
+55
yes
C1:Hydrophone Manufacturer
and Type
Yes
No
B2: Do you check the hull surface condition?
A10: What kind of environmental factors are measured during noise measurements?
Water depth
Water density
Current velocity
Water temperature
Temperature distributions
Sea state
Air temperature
France
COUNTRY
-174, -205
-40
+80
150 m
5
160k
B&K
8106,8105,
8104
no
No
TNO
DMO
Netherlands
-170
20
200k
B&K
no
According to
vessel
maintenance
planning
According to
vessel
maintenance
planning
GPS location
TSI
Spain
0.05
0.35
+/- 1dB
at 10 kHz
+/- 0.5dB
at 56 kHz,
- 186
0
+25
up to 400 m
5
56 k
yes
NSWC/CD
USA
10 10 k
Others (specify)
Others (specify)
D8: Type of Sound Source Localization Methods
Beamforming
Crosspower spectrum
Others
No
No
108
10 50 k
Peak
Max. 3
Hydrophone
calibrator
yes
CETENA
Italy
BE
4 k - 80 k
1.6 - 10 kHz
20
yes
WTD 71
Germany
yes
DCNS
ORGANISATION
C4: Calibration
Do you have hydrophone calibration
procedures?
How often do you calibrate? Before /
After measurements, Regularly
France
COUNTRY
3-5
No
100 80 k
Yes
BE
40
- 20k
Depends on
application
Yes
1 - 100k
RMS
Parallel with
over side test
at sea
No
BE/EE
102
0-100k
RMS
No
B/E
144
10 - 80 k
Depends on
application
typically 3
(more possible)
B,A
yes
TNO
One or more if
necessary
B,A
yes
HHI
Hydrophone
calibrator
B,A
yes
KRISO
DMO
Netherlands
Hydrophone
Hydrophone Hydrophone
calibrator/
calibrator
calibrator biannual factory
calibration
yes
Mitsui Lab.
Korea
Water tank
yes
MHI
Japan
No
20-100 k
Postprocessing
Hydrophone
calibrator
yes
TSI
Spain
Yes
B/E
5 - 56 k
Depends on application
(for cavitation)
yes
NSWC/CD
USA
DCNS
ORGANISATION
No
Yes
Mitsui Lab.
Japan
No
Yes
KRISO
Yes
No
HHI
1-2
Time signal
Korea
hemispherical
diffusion
Yes
Numerical
computation
No
No
Yes
Yes
Yes
No
No
No
No
No
No
No
1-2
No
No
No
Yes
No
Yes
Yes
No
No
TNO
DMO
Netherlands
depends on
application.
Conversion to
'dipole source
level'.
No
1.5
No
Yes
No
MHI
No
CETENA
WTD 71
No
Italy
Germany
Others (specify)
Harmonics
1/3 octave converted to 1Hz
Normalization to 1 m distance
Others
D12: Expected Uncertainty (dB)
for Full Scale Sound Pressure
+/- 1
Level Measurements
D13: Confidence Level on the Quality of Test Results
Scale of 1~10: 1 very uncertain,
8
10 very confident
France
COUNTRY
Yes
No
Yes
No
Yes
TSI
Spain
+/- 3
No
No
Yes
No
No
NSWC/CD
USA
12
0
200
15
10
350
60
15
2.6
0.6
0.6
10.5
2.2
2
CSSRC
(2001)
SSSRI
(1985)
HSVA
(1990)
12
25
250
25
11
2.8
1.6
INSEAN
(1978)
5
3
101.
10
3.6
2.25
U. Genova
(1979)
8.5
10
101
20
1.2
0.57
0.57
U. Sharif
4
20
200
3.2
0.63
0.63
JMUC
(1980)
MHI
(1970)
12
20
100
20
2.8
0.71
0.71
14
9.81
196.
14
2.6
0.6
0.6
15
10
300
10
10
2
2
A: FACILITY INFORMATION
KRISO
(2009)
16.5
20
350
10
12.5
2.8
1.8
12.7
9.8
196
2.6
0.6
0.6
HHI
(1984)
Korea
Netherla
Norway
nds
6
2.5
101
240
18
8
MARIN
(1973)
Japan
18
10
600
1.2
2.2
MARINTEK
(1967)
Iran
MEGURO
(2006)
NB. A5b was introduced later following comment from HSVA and SSPA, answers correspond to answer a) in question A4.
ORGANISATION
Italy
Russia
10
30
300
10
6
1.3
1.3
KRYLOV
(1961)
Germany
Sweden
6.9
20
300
20
9.6
2.6
1.5
SSPA
(1970)
China
10
10
100
25
3.5
0.8
0.8
Rolls-Royce
COUNTRY
4
100
101
40
2
0.63
0.35
18
3.5
415
12.4
3.048
3.048
USA
Turkey
Istanbul TU
NSWC/CD
(1991)
B2: Appendages
Rudder
Stabilizer fins
Thruster openings
Energy saving devices
Shaft brackets
Bilge keels
Sonar dome
B3: Turbulence Stimulators
Ship hull (None, Studs, Trip Wire,
Other)
Propeller (None, Carborundum,
Paint)
B4: Propeller Manufacturing
Method (NC+ hands, NC only)
Accuracy (, mm)
B5: RPM measurement
Location (Dynamometer,Motor,
Shaft)
Number of pulses / rev
ORGANISATION
CSSRC
250
Al
250
brass
NC+
0.05
NC+
0.05
36
NC+
0.1
1024
360
3600
NC
0.02
Bro, Al
250
W/FRP
5-6
INSEAN
Brass
250
GRP
11
6 - 10
U. Genova
100
NC+
ITTC
various
200250
Italy
U. Sharif
NC
0.05
steel
300
MHI
JMUC
NC+
0.05
Al.
250
tacho
NC+
0.05
B: MODEL SETUP
MEGURO
360
NC+
1024
NC+
0.02
Al.
250
W / FRP
6-9
KRISO
Al.
250
FRP
HHI
100
NC+
0.05
Al.
250
Korea
360
NC+
0.05
Bro.
230250
7 - 12
MARIN
Netherla
Norway
nds
3600
NC+
0.02
Bro.
250
Various
1.2
MARINTEK
Japan
Russia
Sweden
Div
NC+
0.02
100
NC+
0.05
200300 200250
babbitt,
duralumi Bro.
n
Foam
KRYLOV
Iran
SSPA
Germany
HSVA
SSSRI
China
720
NC+
0.05
Bro.
250
Div
3.5
RollsRoyce
COUNTRY
NC+
0.05
Bro.
150
Turkey
Istanbul TU
15 or 60
NC+
0.075
Al.
300460
GRP
6.1 - 12.2
USA
NSWC/CD
C4: Calibration
Do you have hydrophone calibration
procedures?
How often do you calibrate? Before
/After measurements, Regularly
How do you calibrate?
Hydrophone calibrator / water tank
C5: Hydrophone Location
Acoustic chamber (AC)
Outside wall or windows
Flushed wall or windows
Rake / in flow
Inside basin
C6: Dimensions of Location
Height (m)
Width (m)
Depth (m)
Stand-off distance (mm)
C7: Number + Position
How many hydrophones are
installed?
Upstream of the propeller plane
At propeller plane
Downstream of the propeller plane
Others
C8: Hydrophone Array
Do you use an array?
Planar type, combo random, spiral
arm log spaced
Hydrophone spacing (mm)
What is the design frequency of the
hydrophone array (kHz)
What is the half power spatial
bandwidth?
What is the aperture size? (m)
HC
0.22
0.5
0.22
420
No
HC,WT
2
2.2
9,5
1000
Yes
0.512.5
10,20,40
PT
Yes
8104
B&K
SSSRI
Yes
B&K
8105
Hydrophone type
CSSRC
ORGANISATION
China
Reson
HSVA
No
2.4
11
2.8
930
HC
Yes
TC 4032
B&K
INSEAN
No
1
1
1
2
HC
Yes
8103
U. Genova
No
0.74
0.26
0.08
300
HC
Yes
TC 4013
Reson
Italy
MEGURO
JMUC
MHI
U. Sharif
No
4
1
HC
B, A
Yes
8103
B&K
No
WT
Yes
8103
B&K
CO
Yes
1
2
1
48
11.8
2
1.5
HC
Yes
B&K
8103,
8105
KRISO
to 5
1.6 x 1.6
to 100
random 120170
circle
Yes
2
1
32
32
8
1.4
1
300700
HC
Yes
B&K
8103,
8104
C: HYDROPHONE INFORMATION
HHI
No
0.08
0.08
0.068
HC
B, A
Yes
8104
B&K
Korea
Netherla
Norway
nds
No
HC
Yes
No
HC
Yes
Reson
B&K
TC 4014various
1
MARIN
Japan
MARINTEK
Iran
Russia
No
0.4
0.6
0.35
900
Yes
8103
B&K
KRYLOV
Germany
SSPA
No
1
1
1
HC
Yes
8103
B&K
Sweden
No
2
516
various
HC
B, A
Yes
B&K,PCB
8103,
112A22
RollsRoyce
COUNTRY
No
1
1
HC
B, A
Yes
8103
various
various
to 50
Other
S-A
Yes
various
56
13.71
2.14
1.15
550
No
TC 4056
Reson
Turkey
Istanbul TU
B&K
USA
NSWC/CD
CSSRC
ORGANISATION
No
Yes
18
1535
0.12
Yes
Yes
No
Yes
1722
1824
030
No
13.5
815
No
No
0.7R
12:00
7090
4060
58
2535
0.40.9
No
No
Yes
No
Yes
0.8R
12:00
No
0.8R
12:00
Kt
1-Wt
Kt
1-Wt
Kt
1-Wt
M
Yes
Yes
SSSRI
HSVA
Yes
No
SCL
No
Kt
1-Wt
Yes
U. Genova
No
1525
40
Yes
56
2030
0.11
Italy
MEGURO
JMUC
MHI
U. Sharif
No
0.63.6
333
No
Yes
0.7R
12:00
No
Yes
No
50
Yes
25
No
No
0.7R
12:00
No
Kt
1-Wt
No
3050
No
No
0.9..0.95
r
Yes
Kt
Yes
various
Yes
1723
3080
Yes
310
2030
0.32
No
Yes
0.7R
12:00
No
Kt
1-Wt
Yes
D, M
Yes
KRISO
Yes
No
0.7R
12:00
Kt
1-Wt
Yes
HHI
Netherla
Norway
nds
Yes
No
SCL
Yes
Yes
Yes
No
0.7R
12:00
No
Kt
1-Wt
Yes
Yes
1025
3070
Yes
Yes
1017
3640
No
No
1018
2540
Yes
No
1518
2530
Yes
712
2.55
23
25
25
3040
815
1825
0.42.5 0.30.6 0.050.1 0.20.3
No
No
0.7R
12:00
Yes
Kt
1-Wt
Korea
MARIN
Japan
MARINTEK
Iran
Russia
Yes
1725
3040
Yes
48
2540
0.33
No
Yes
SCL
No
Kt
1-Wt
No
KRYLOV
Germany
INSEAN
Yes
China
SSPA
No
1525
3040
Yes
49
2035
SCL
Yes
Yes
Sweden
No
1020
20
No
60
Yes
23
1733
1.06
No
Yes
0.7R
12:00
No
Kt
1-Wt
Yes
Yes
2135
6080
Yes
518
Yes
Yes
SCL
No
Kt
1-Wt
Yes
USA
Turkey
Istanbul TU
3040
Yes
45
2025
0.2
Yes
Yes
0.7R
12:00
No
Kq
1-Wq
Yes
RollsRoyce
COUNTRY
NSWC/CD
A
Yes
B
Yes
CSSRC
SSSRI
D, +
HSVA
A
Yes
INSEAN
A
Yes
D, +
Italy
A
Yes
*:
U. Sharif
MHI
B
Yes
MEGURO
B, A
Yes
JMUC
20 s
30 revs
Yes
Yes
1
100
Yes
10
20
Yes
No
No
HHI
B, A
Yes
R, I
Netherla
Norway
nds
B, A
Yes
R, I
D, +
B, A
Yes
R, I
Russia
B, A
Yes
SSPA
A
Yes
Sweden
A
Yes
D, +
B
Yes
Yes
3060 s
No
20
12.8
Yes
10 s
No
1001000
100
Yes
20 s
Yes
50
20
Yes
30 s
No
1 000
80
Yes
0
20
Yes
Yes
150300
7.530 s
s
No
1
22
Yes
No
20
12.5
Yes
> 10 s
Yes
1000
50
1m
1Pa/Hz
B
Yes
R, I
USA
Turkey
Istanbul TU
No
10
40
10
50
1Pa/Hz
RMS
A
Yes
R, I
KRISO
A
Yes
Korea
MARIN
Japan
MARINTEK
Iran
KRYLOV
Germany
U. Genova
Others (specify)
E8: Sound Source Localization Method
(a)Beamforming
(b)Matched-filter array processing
(c)Crosspower spectrum
(d)None
(e)Others (specify)
*
ORGANISATION
China
RollsRoyce
COUNTRY
NSWC/CD
Italy
SSSRI
CSSRC
HSVA
1-2
5
KRISO
3-5
TS
Netherla
Norway
nds
HHI
Russia
10
**: Evaluate vibration influence using coherence measurements (subtract coherent power)
***: Correction to far-field based on calculations or in-situ reciprocity measurement
identify
noise
source
location
MHI
JMUC
NB
U. Sharif
MEGURO
1/3 octave
Narrowband normalized to 1 Hz
Harmonics
1/3 octave converted to 1Hz
Normalization to 1 m distance
Others: time signal (TS),
narrowband (NB)
E12: The Reverberation of Facility
Have you investigated the
reverberation of facility
E13: What Corrections
(a) Free surface
(b) Reflections due to wall
(d) Reverberation
U. Genova
Korea
MARIN
Japan
MARINTEK
Iran
KRYLOV
Germany
INSEAN
vibration/acceleration
machiner
measurements(specify)
y noise
filtering
ORGANISATION
China
***
SSPA
Sweden
RollsRoyce
COUNTRY
**
Turkey
Istanbul TU
USA
NSWC/CD
CSSRC
empirical
Yes
SSSRI
Yes
HSVA
Yes
No
Italy
U. Genova
JMUC
MHI
U. Sharif
0
2
3
3
No
similar to
ITTC
1987
Yes
*: Depends on type of noise mechanism. In-house cavitation relationships & non-cavitating scaling relationship
No
Yes
MEGURO
No
F: SCALING METHODS
KRISO
No
Yes
Korea
HHI
No
Netherla
Norway
nds
Yes
Russia
Yes
10
10
No
Yes
1
1
2
1
No
Sweden
Yes
SSPA
0.5
1
1.5
1
No
de Bruijn own
own
and ten facility facility
Wolde scaling scaling
(1974) method method
Yes
MARIN
Japan
MARINTEK
Iran
KRYLOV
Germany
INSEAN
Others (specify)
ORGANISATION
China
No
Yes
RollsRoyce
COUNTRY
0.5
1
2
1
No
No
Yes
Turkey
Istanbul TU
Yes
USA
NSWC/CD
BACKGROUND
Membership
1.2
CNR-INSEAN, Italy
February 2012
Dr Aurlien Babarit,
cole Centrale de Nantes, France
Dr Arnold Fontaine,
Pennsylvania State University, USA
Prof Yanping He
Shanghai Jiao-Tong University, China
SCOPE
Dr Marek Kraskowski,
Centrum Techniki Okrtowej (CTO),
Poland
Dr Francesco Salvatore,
Istituto Nazionale per Studi ed Esperienze
This report addresses a number of key issues in the physical and numerical testing of
marine renewable energy systems, including
wave energy devices, current turbines, and offshore wind turbines. The report starts with an
overview of the types of devices considered,
and introduces some key studies in marine renewable energy research. The development of
new ITTC guidelines for testing these devices
is placed in the context of guidelines developed
680
3
3.1
3.2
OVERVIEW
Technology Readiness Level
In the case of the renewable energy industry, the following stages of Technology Readiness Levels (TRLs) are commonly considered
(e.g. Mankins (2009)). TRL 1-3 correspond to
research stages up to and including proof of
concept, TRL 4-5 correspond to component,
sub-system and system validation in laboratories and/or simulated operational environments
and TRL 6-9 correspond to prototype demonstration in operational environment through to
system proving via successful deployment.
681
overtopping
(a)
(b)
(c)
Figure 1. Overtopping Devices: (a) Schematic Diagram; (b) TAPCHAN built onshore in Norway
in the 1980s; (c) 1/4.5 scale model of floating device Wavedragon deployed in Denmark in 2003.
Turbine
Incident waves
Chamber
(a)
(b)
(c)
Figure 2. Oscillating Water Column Devices: (a) Schematic; (b) Pico shore-based OWC built in
the Azores in 1999; (c) Oceanlinx floating OWC deployed in 2010 in New South Wales Australia
(a)
(b)
(c)
(d)
Figure 3. Oscillating Body Devices: (a) Schematic; (b) Carnegies Ceto heaving buoy; (c)
Aquamarines Oyster device; (d) Pelamis Wave Powers P2 device. Devices c) and d) have
been tested at the European Marine Energy Centre (EMEC) full-scale test site in the UK.
682
Oscillating Bodies: In these devices, incident waves make one or several bodies oscillate. Relative motions between the bodies and
the sea bottom or between the bodies themselves are used to drive a Power Take Off system (PTO), often based on hydraulic components. The working principles and examples of
well-known prototypes are shown in figure 3.
WECs may also be classified using the location of installation. Some wave energy converters are designed to be installed at the shoreline, some in near-shore shallow-water regions
while other can be installed in deep water offshore.
A final approach to classification of WECs
uses considerations of size. Devices of small
dimensions with respect to wavelength are
called point absorbers. Examples are Carnegies Ceto device or the Aquamarines Oyster.
Large devices with the longest dimension parallel to the wave crests are called terminators, whereas devices with the longest dimension parallel to the wave propagation direction
are called attenuators. The Wavedragon is an
example of a terminator and the Pelamis is an
example of an attenuator.
683
work initiated a considerable amount of research on wave energy until the mid-1980s,
when funding stopped partly because of the
decline in the oil price. References to many
interesting studies which were conducted at
that time can be found in the review paper by
Falnes (2007) and in McCormick (1981), Berge
(1982) and Evans & Falcao (1985).
3.3
Current Energy
684
Energy conversion mechanisms are generally fully submerged and may be positioned at
various depths from the free surface. Existing
technologies include bottom-fixed devices as
well as floating and mid-water devices in
which one or more turbines or equivalent
mechanisms are fitted to surface platforms or
submerged moored structures. Figure 5 shows
some examples of turbine and non-turbine concepts. A review of marine current device technologies can be found in Khan et al. (2009).
(a)
(b)
(d)
(c)
(e)
(f)
(g)
(h)
(i)
(j)
685
Each of these two layouts presents advantages and disadvantages. In terms of power
efficiency, that is the amount of power that is
produced for given momentum of the water
mass processed by the device, HACTs have in
general higher performance than VACTs. Conversely, HACT power capturing capability is
very sensitive to the alignment of turbine axis
and current direction, whereas VACTs are insensitive to the direction of the onset flow in a
plane normal to the turbine axis.
For this reason, the latter are frequently preferred in case of bidirectional tidal currents,
whereas the operation of HACTs in tidal current sites implies that suitable blade/turbine
orientation solutions are implemented. A simple solution for HACTs is to adopt blades with
bidirectional profiles that, at a price of reduced
hydrodynamic efficiency, do not require to be
oriented according to tide direction.
In the case of deployment in shallow water
sites, bottom-fixed installations of HACTs are
preferred, but many examples of floating turbines exist. VACTs are typically appended to
floating platforms with the advantage that generation set and power control systems can be
placed on the platform above the water surface
yielding inherent advantages in terms of inspection and maintenance operations. A disadvantage of floating devices is vulnerability in
case of extreme weather conditions (e.g. hurri-
686
One of the projects most advanced in development is the MCT SeaGen shown in Figure
7, a 1.2MW horizontal axis tidal device deployed and operational since 2008 at Strangford Lough, Northern Ireland (Fraenkel
(2010)). The system consists of twin power
trains mounted on a crossbeam supported by a
bottom fixed tower. The cross beam can be
raised above the free surface for maintenance.
Each turbine has a diameter of 16 m and two
blades that can be pitched through 180 degrees
to operate in bi-directional flows with current
speed of 2.5 m/s. The company has plans to
deploy a 10MW turbine array in Wales and an
8MW tidal farm in Scotland. Other examples
of smaller horizontal-axis turbines are operating and grid-connected along the UK and Canadian coasts.
(a)
(b)
687
(a)
(b)
(c)
Figure 8 Examples of horizontal turbine prototypes
under assessment in field tests: (a) Voith Hydro
Test Turbine, Jindo, Korea; (b) Andritz Hydro
Hammerfest HS1000, EMEC, Scotland; (c) Open
Hydro Test Turbine, EMEC Scotland
688
689
piles, tripods, jackets and gravity bases. Floating Offshore Wind Turbines (FOWTs) utilise
floating support structures in many cases utilising technology originally developed for the
offshore oil and gas industries. In the windturbine literature these are often categorised as
ballast-stabilised (e.g. Spars), Mooringstabilised (e.g. TLPs) and buoyancy-stabilised
(e.g. barges and semi-submersibles). While
bottom-mounted systems are currently operating commercially in several European countries, FOWTs are still in a process of evolution
with several technology demonstration projects
under way around the world. Some hybrid concepts are under development, attempting to
exploit potential for synergy between wind and
current or wave energy systems. As of 2013 the
two largest offshore wind farms are located in
UK waters: the 630 MW London Array (Figure
11) is the largest offshore wind farm in the
world, with the 504 MW Greater Gabbard
wind farm the second largest. All turbines are
of the bottom-mounted monopile type.
3.4
Wind energy
690
The increasing size of full-scale wind turbines requires that experiments adopt relatively
small scale models. Devices currently being
deployed employ turbines up to 170m in diameter, and sizes may continue to increase as
offshore installations become more economically attractive both in land and offshore. Consequently it is difficult to keep the physical
similarities of not only the scale of the model
but also the external wind load conditions.
Most of the experiments of 2-7MW FOWTs in
a conventional ocean basin or test tank now use
1/50-1/200 scale models. Some particular challenges of these tests are addressed in Section 7.
(a)
(b)
691
(a)
(b)
(c)
(e)
(d)
(f)
692
GUIDELINE DEVELOPMENT
ers, sponsors/investors and/or regulators on the
processes required in order to assess device
performance at various stages of device design
and development, including hydrodynamic
model testing. However, much of the guidance
given is rather general and does not address the
specific nature of many critical challenges of
tank testing MRE devices.
693
Organisation
Type
Status(Reference)
Title
Infor
mal
Liaison
Yes
Organisation
EuropeanMarineEn
ergyCentre(EMEC)
Equimar
Marinet
SupergenMarine
OpenSeatest
centre(full
scale&nurs
ery)
EUfunded
ResearchPro
ject
EUfunded
ResearchPro
ject
UKresearch
project
Infor
mal
Liaison
Yes
TankTestingofWave Published2009
EnergyConversionSys http://www.emec.org.uk/tank
tems
testingofwaveenergy
conversionsystems/
Bestpracticefortank Published2010
Yes
testingofsmallmarine http://www.equimar.org/equi
energydevices
marprojectdeliverables.html
Ongoingdevelopment Underdevelopment
Yes
ofguidelinesforwave
andcurrentenergy
testing
Guidancefortheex
Published2008
No
perimentaltanktesting http://www.supergen
ma
ofwaveenergycon
verters
rine.org.uk/drupal/files/reports
/WEC_tank_testing.pdf
694
5
5.1
One of the key differences between the behaviour of a marine renewable energy device
and other fixed or floating structures subject to
fluid action is the presence of a Power Take
Off (PTO) system specifically designed to extract energy from the interaction between the
device and the fluid flow. The device performance, in terms of power capture as well as other
aspects of the device response such as motions
and/or hydrodynamic loads, thus depends upon
the behaviour of the PTO system. Hence appropriate simulation of the power take-off is
essential during small-scale model tests to determine the performance of the system.
5.2
695
The more complex strategies typically employ phase control with highly-tuned devices,
such as point absorbers, in order to broaden the
frequency range at which high response occurs.
In latching control (e.g. Bjarte-Larsson & Falnes (2007), Durand et al. (2007)) the relative
696
es are typically designed to operate over a specific range of conditions, typically RPM or
oscillation frequency, for optimal power generation. Maintenance of optimal conditions is
usually achieved by loading the device to hold
RPM or oscillation to a desired range. An unloaded device would typically operate under
freewheeling conditions often undesirable at
full scale and producing minimum shaft power.
The PTO system, comprised of drive train,
power generation and power electronics, is the
sub-system that provides the necessary device
loading. Proper small-scale device testing must
include some form of PTO modeling.
The important parameters of interest relative to PTO function in a small-scale model test
of a rotating device are shaft RPM and torque
at the rotor to shaft attachment. The smallscale PTO must be designed to absorb the delivered power by the shaft under controlled
conditions to minimize shaft torque and rpm
drift during testing. In tests of model current
devices, the PTO can be represented by direct
electrical power generation, by mechanical /
hydraulic / magnetic loading or by using a
speed or torque control drive to control rotor
RPM.
5.3
Generators, either permanent magnet or inductance, can be used in small scale model
testing using direct-drive or gear-box coupling.
While full-scale devices often have the power
generator installed in the nacelle, space limitations in small-scale model testing may require a
revised configuration where the generator is
installed in a downstream dynamometer or outside the facility. These adaptations often involve additional components like seals, bearings and gearbox configurations.
697
5.4
698
6.1
NUMERICAL MODELLING OF
MARINE RENEWABLE ENERGY
DEVICES
699
700
be adapted to deal with this feature, using approaches such as generalised modes. Dedicated
solvers may be developed in particular cases
(Renzi & Dias, 2012). In case of overtopping
devices, linear potential theory cannot be used
and one may rely on empirical laws (Borgarino
et al., 2007).
Although linear potential flow theory is often the only practical option, its limitations
must be acknowledged. WECs are often designed to have their natural frequencies within
the wave spectrum. Thus, the motion response
of the device may be large, which violates the
assumption of small amplitude motion. Significant discrepancies in comparison with experiments may be observed (Durand et al., 2007).
Some research groups are putting efforts in
developing medium fidelity models to address
this issue, while keeping the computational cost
moderate (Gilloteaux et al., 2007, Guerber et
al., 2012).
Some WECs are composed of many articulated bodies with many more degrees of freedom than the conventional 6 DOFs for conventional rigid body motions (Soulard & Babarit,
2012, Babarit et al., 2013). BEM solvers must
701
Linear Bf,M
Fundamental
Definition ol
hydrodynamics
Explicit coefficients
Nonlinear BEM
Boussinesq
Mild-slope
Not capable
Implicitly
lmplicitly
capable
capable
Not capable
Not capable
Implicit solver
Explicit inclusion
Explicit incusion
Nonlinear dynamics
Vortex shedding
WEC radiation
Implicitly capable
Explicitly capable
Diffraction
Implicitly capable
Explicitly capable
Not capable
lmplicity
Implicitly
capable
capable
Supra-grid
Sub-grid
CFD
Explicit
Explicit source
absorption laver
strength
Implicit fluid
flow
lmplicitly
pbase-
capable
averaged dynamics
Explicit
absorption layer
Explicit
source
strength
lmplicit inclusion
Explicit inclusion
Not capable
Explicitly
capable
Approximated by phase-decoupled
refracti on-diffracti
Implicit solver
lmplicitly capable
Impticitly
capable
Implicitly
capable
Implicitly
capable
ComDutation
Primary dependent
Secondary dependent
Complexity of
function
Number of panels
Determinate of array
Number ofpanels
and complexity
of equations
Number of panels
Number oltime-steps
"size"
Sover
Number of cells
Number of ceils
Number of cells
Number of time-steps
Number of time-
Complex and
stable
Medium
Low
steps
poss unstable
Usabil
Required skill
Software availability in
2012
Low
High
Medium
High
Commercial code
Research code
Commercial code
Research code
available
only
available
only
model required
t**
Limited to mid-slope
702
Low
Medium
High
Commercial and
open-source code
available
With few exceptions, turbine and nonturbine systems have in common that power is
generated by means of lifting surfaces subject
to rotary or oscillatory motions. According to
the device type, these surfaces consist in turbine blades, foils or sails. The prediction of
hydrodynamic forces generated on these lifting
surfaces can be obtained by computational
models with different levels of approximation
used to describe relevant fluid-dynamics
mechanisms. Due to similarities of energy capturing mechanisms and of device layouts to
some extent, computational methods for marine
current devices are closely related to models
used for wind turbine modelling.
In general, computational modelling provides analysis and design tools that are complementary to testing physical models in laboratories and field sites. In particular, aspects
that cannot be reproduced by physical tests can
be investigated by numerical simulations. This
includes, for instance, estimation of scale effects on results of laboratory tests carried out
on small models in order to provide reference
data for the design of full-scale prototypes to
be deployed in open water. The analysis of
multiple device operations in arrays and the
impact of energy capture on the environment
are other examples of applications of computational models.
703
The most popular approach among viscousflow solvers is the Reynolds-Averaged NavierStokes Equation (RANSE) method. This methodology allows to estimate the effects of flow
vorticity and turbulence on hydrodynamic
loads. Turbulence modelling is introduced to
avoid the explicit solution of eddies in the flow.
If more detailed descriptions of local flow perturbations and of blade/foil generated wakes
are necessary, Detached Eddy or Large Eddy
Simulation model (respectively, DES and LES)
are preferred to RANSE.
The problem of predicting the trailing vorticity pattern shed by VACT blades is addressed by Li and Calisal (2010a, 2010b, 2011)
using a 3D potential-flow vortex method with
viscosity-effects correction. The impact of the
wake shed by one turbine blades on the following blades is analysed for different configurations and the best performing layouts are determined.
Viscous-flow correction models are used to
improve turbine blade load predictions from
inviscid-flow models that cannot describe blade
stall and post-stall conditions when blades are
at high angle of attack, a common condition in
case of operation at small TSR. In particular,
704
water channel. To this purpose, the computational domain reproduces the test section of the
facility. The turbine region is enclosed into a
cylindrical grid block that rotates with respect
to the rest of the grid fixed with channel walls.
Suitable boundary conditions are imposed at
the interface between rotating and fixed grids.
Numerical results include pressure and loads
distributions on turbine blades and flowfield
quantities for given current speed and different
values of the Tip Speed Ratio (TSR).
705
Figure 15. Vorticity field and streamlines calculated across a VACT (Maitre et al., 2013)
Furthermore, Gaden and Bibeau (2010) present the results of a computational study using
the commercial RANSE code ANSYS CFX of
the flow around an axisymmetric shroud. Turbine induction is estimated by a simplified
momentum source model and is not explicitly
solved by RANSE. This turns into a computationally efficient methodology that can be used
to perform parametric studies to determine duct
geometries capable to maximise turbine power
output for a given current energy density. The
capability to achieve power output gains up to
a factor 3 with respect to unducted turbine op-
706
posite materials for blades. Dealing with composites, hydroelastic effects on blades are to be
taken into account because of the large deflections that blade loading determine. A hydroelastic analysis of a horizontal axis turbine is
described by Nicholls et al. (2013). The FluidStructure Interaction model combines a LSM to
predict blade loads and a structural model from
the Finite Element commercial code ANSYS.
Numerical results indicate that hydroelastic
tailoring consisting of blade bend-twist coupling has a potential to reduce blade loading
and to increase efficiency and reliability.
707
6.2
Offshore Wind Turbines (OWTs) are generally designed using simulation codes which
aim to model the coupled effects of wind inflow, aerodynamics, hydrodynamics, turbine
control systems and structural elasticity. These
are known as aero-hydro-servo-elastic codes;
some examples are shown in Table 1 (Robertson et al., 2014). Due to the high level of complexity of these codes it is desirable that they
are verified and validated to ensure their accuracy. A series of research tasks have been developed under the International Energy Agency
(IEA) Wind to address this need.
The first task was designated the Offshore
Code Comparison Collaboration (OC3) which
operates under Subtask 2 of the International
Energy Agency (IEA) Wind Task 23. The OC3
programme consisted of four phases each analysing a different problem. The basis of the
study was a code-to-code comparison. In Phase
I, the NREL 5-mW wind turbine was simulated
in 20m of water on a monopile with rigid foundation; in phase II the foundation was made
flexible using different models to represent
soil-pile interactions. In Phase III, the water
depth was changed to 45 m and the monopile
was replaced with a tripod substructure; whilst
in Phase IV, the wind turbine was installed on a
floating spar-buoy in deep water (320 m). Full
details may be found in Jonkman & Musial,
(2010).
708
7.1
Uncoupled Tests
Measurement of wave / current loads without simulation of wind loads offers a number of
advantages. Facilities without wind generation
may be employed, set-up and calibration time
will be reduced, whilst a larger model can be
used for the measurements of wave/current
loads. Furthermore, measurement uncertainty
may be reduced since the entire range of the
load measurement systems can be used for
measuring the wave/current loads. However,
where tests are aimed at investigating the coupled dynamic response of the structure in operational conditions, including realistic modelling of flexibility and aerodynamic damping,
then inclusion of the aerodynamic coupling due
to the rotor is necessary.
In studies of floating offshore wind turbines, model tests without the rotor can be carried out at preliminary stages of the tests or for
special purposes, such as comparison of differ-
709
ent support structures in terms of their responses to waves, or for the validation of numerical models. However, final tests aiming at
evaluating the global response of the system
from the concept validation stage to the prototype and demonstration stage should include at
least simplified modelling of the rotor due to
the strong coupling present between the platform dynamics and the rotor-generated forces
and moments.
7.2
Simplified Simulation with Wind Generation. Where a simulated wind field may be
created in a tank, it is possible that a solid or
porous disc may be used in place of the rotor in
conjunction with a battery of fans. The disc
should be sized to generate a drag load in the
simulated wind field corresponding to the predicted mean thrust on the turbine. If a rotating
disc is employed with Froude-scaled rotary
moment of inertia, and rotation speed, it is pos-
710
7.3
711
7.4
A second approach addressing low Reynolds number effects is the placement of studs
or other roughened materials as a turbulence
stimulator along the leading edge of a blade;
however this is unlikely to improve the turbine
performance adequately on its own to yield
comparable performance with the full-scale
device, and can yield unrealistic results if laminar separation occurs, as well as unrealistic
unsteady aerodynamic loads during flow reattachment.
712
wind speed. A further example of this technique is given by De Ridder et al. (2013).
Correct modelling of the gyroscopic moments introduced by the rotor can be achieved
by Froude scaling of the mass properties and
rotor speed. However, modelling the mass
properties of the rotor at small scale will require lightweight materials of high strength
potentially presenting some severe challenges
in manufacture. For example, the 1:80 scale
blades shown in Figure 18 have a target mass
of only 35g. It can be difficult to achieve the
necessary combination of accurate geometry
and very low mass using conventional model
making techniques, and novel approaches may
be required: for example, film-coated blades
may be manufactured from components created
by a 3D printer to achieve both mass and elasticity requirements.
7.5
Key issues for design load calculations include the choice of models for wind gradient
and the turbulence. Obrhai (2012) reviews the
current guidelines for wind modelling for offshore wind turbines.
The IEC standard 64100-3 (2009) recommend that the wind speed profile as a function
of height is given by the power law:
V z Vhub z zhub
Contribution of the wind load to global response of the FOWT is also strongly affected
by the pitch control strategy of the turbine. In
the near future, it is possible that 3D printers
will enable the manufacture of both very light
weight models of the RNA (rotor nacelle as-
713
UNCERTAINTIES IN PHYSICAL
MODEL EXPERIMENTS AND
EXTRAPOLATION TO FULLSCALE
However there are some very specific challenges in relation to the analysis of uncertainty
of power output of MRE devices which are not
addressed by existing procedures. These are
discussed further below.
Arguably the area in which greatest development is required is in the prediction of power
from wave energy devices. Very few examples
exist in the literature of uncertainty analysis of
wave energy device tests. Lamont-Kane et al.
714
REFERENCES
For all types of devices additional uncertainty results from extrapolation of model-scale
results to full-scale, particularly in cases in
which multiple scaling ratios are relevant. This
affects wave energy devices with PTOs which
do not behave according to Froude scaling,
such as oscillating water columns, floating or
mid-water current turbine systems and floating
offshore wind turbines, for which for which
both Froude and Reynolds scaling are relevant,
as discussed in sections 6 and 7.
Amelio, M., Barbarelli, S., Florio, G., Scornaienchi, N.M., Minniti, G., Cutrupi, G.,
Snchez-Blanco, M., 2012, Innovative tidal
turbine with central deflector for the exploitation of river and sea currents in on-shore installations. Applied Energy, Vol. 97, pp. 944
955.
Antheaumea, S., Matre, T., Achard, J.-L.,
2008, Hydraulic Darrieus turbines efficiency
for free fluid flow conditions versus power
farms conditions. Renewable Energy, Vol. 33,
pp. 21862198.
Babarit A. and Clment A.H., 2006, Optimal latching control of a wave energy device in
regular and irregular waves, Applied Ocean
Research, Vol 128 (1), pp 56-64
715
Babarit A., Hals J., Muliawan M.J., Kurniawan A., Moan T., Krokstad J., 2012, Numerical benchmarking study of a selection of
wave energy converters, Renewable Energy,
Vol 41, pp 44-63
Bellamy N.W., Peatfield A.M., 1986, Design and performance of the circular sea clam
wave energy converter, Proceedings, 3rd International Symposium on wave, tidal, OTEC and
small scale hydro energy
Babarit A., 2013, On the park effect in arrays of oscillating wave energy converter,
Renewable Energy, Vol 58, pp 68-78
Babarit A., Gendron B., Singh J., Melis C.,
Jean P., 2013, Hydro-elastic modelling of an
electro-active wave energy converter, Proceedings, 32nd International Conference on
Offshore Mechanics and Artic Engineering,
Nantes, France
716
Cruz J., 2008) Ocean wave energy: Current Status and Future Perspectives, Springer,
ISBN 978-3-540-74895-3
De Ridder, E.-J., Otto, W., Zondervan, G.J., Huijs, F., Savenije, F., 2013, State of the
Art Model Testing Techniques for Floating
Wind Turbines, Proceedings, EWEA Offshore 2013, Frankfurt, Germany
Durand M., Babarit A., Pettinotti B., Quillard O., Toularastel J.L., Clment A.H., 2007,
Experimental validation of the performances
of the SEAREV wave energy converter with
real time latching control, Proceedings, 7th
European Wave and Tidal Energy conference,
Porto, Portugal
Clment A.H. and Babarit A., 2012, Discrete control of resonant wave energy devices,
Proceedings, Philosophical Transactions of the
Royal Society A, Vol 370, pp 288-314
Cretel J., Lewis A.W., Lightbody G., Thomas G.P., 2010, An application of model predictive control to a wave energy point absorber, Proceedings, IFAC conference on control methodologies and technology for energy
efficiency, 267-272
Evans D.V., 1976, A theory for wavepower absorption by oscillating bodies, Journal of Fluid Mechanics, Vol. 77, pp 1-25
717
Evans D.V., 1979, Some theoretical aspects of three dimensional wave energy absorbers. Proceedings, 1st Symposium on Wave
Energy Utilization, Gothenburg, Sweden
Falcao A.F. de O., Henriques J.C.C., Candido J.J., 2012, Dynamics and optimization of
the OWC spar buoy wave energy converter,
Renewable Energy, Vol 48, 369-381
Falnes J., 1980, Radiation impedance matrix and optimum power absorption for interacting oscillators in surface waves. Applied
Ocean Research, Vol 2, 75-80
Falnes J., 2002, Ocean waves and Oscillating Systems, Cambridge University Press.
ISBN 0-521-78211-2
Falnes J., 2007, A review of wave energy
extraction, Marine Structures, Vol. 20,
pp 185-201
718
Fontaine, A., Straka, W., Meyer, R., Jonson, M., 2013, A 1:8.7 Scale Water Tunnel
Verification & Validation of an Axial Flow
Water Turbine Technical Report Pennsylvania
State University and Sandia National Laboratories (SNL). 122 pp
Forestier, J. M., Holmes, B., Barret, S. and
Lewis, A., 2007, "Value and validation of
small scale physical model tests of floating
wave energy converters," Proceedings, 7th
European Wave and Tidal Energy Conference,
Porto, Portugal, 11-14th Sep. 2007.
Goupee A.J., Koo B., Kimball R.W., Lambrakos K.F., and Dagher H.J., 2012, Experimental Comparison of Three Floating Wind
Turbine Concepts, Proceedings, 31th International Conference on Ocean Offshore & Arctic
Eng., OMAE2012-83645 Rio de Janeiro, Brazil
719
Khan, M.J., Bhuyan, G., Iqbal, M.T., Quaicoe, J.E., 2009, Hydrokinetic energy conversion systems and assessment of horizontal and
vertical axis turbines for river and tidal applications: A technology status review, Applied
Energy, Vol. 86, pp 18231835.
Jonkman, J., & Musial, W., 2010, Offshore Code Comparison Collaboration (OC3)
for IEA Task 23 Offshore Wind Technology
and Deployment, Technical Report NREL/TP5000-48191 71 pp
Li, Y., Calisal, S.M., 2010a, Threedimensional effects and arm effects on modeling a vertical axis tidal current turbine, Renewable Energy, Vol. 35, pp 2325-2334.
720
Mason-Jones, A., O' Doherty, D. M., Morris, C. E., O' Doherty, T., 2013, Influence of a
velocity profile and support structure on tidal
stream turbine performance, Renewable Energy, Vol. 52, pp 23-30.
Luquet, R., Bellevre, D., Frchou, D., Perdon, P., Guinard, P., 2013, Design and model
testing of an optimized ducted marine current
turbine, International Journal of Marine Energy, Vol. 2, pp 6180.
McCormick M.E., 1981, Ocean Wave Energy Conversion, New York: Wiley
Meyer N. I., McDonal Arnskov, M., Vad
Bennetzen, L. C. E., Burcharth, H. F., Bunger
J., Jacobsen, V., Maegaard, P., Vindelok, S.,
Nielsen, K., Sorensen, J.N., 2002, Bolgekraftprogram: Afsluttend rapport fra Energistyrelesens Radgivende Bolgekraftudval, Technical
Report, Bolgekraftudvalgets Sekretariat, Ramboll, Virum, Denmark
721
Mycek, A.B., Gaurier, B., Germain, G., Pinon, G., Rivoalen, E., 2013, Numerical and
experimental study of the interaction between
two marine current turbines, International
Journal of Marine Energy, Vol. 1, pp 7083.
Parmeggiani, S., Kofoed, J. P. & FriisMadsen, E., 2013, Experimental Update of the
Overtopping Model Used for the Wave Dragon
Wave Energy Converter Energies, Vol. 6,
pp 1961-1992;
Payne G. S., Taylor J. R. M., Bruce T.,
Parkin P., 2008, Assessment of boundaryelement method for modelling a free-floating
wave energy device. Part 2: Experimental validation, Ocean Engineering, Vol 35, 342-357
Myers, L.E., Bahaj, A.S., 2012, An experimental investigation simulating flow effects in first generation marine current energy
converter arrays, Renewable Energy, Vol. 37,
pp 28-36.
722
Popko, W., Vorpahl, F, Zuga, A., Kohlmeier, M., Jonkman, J., Robertson, A., Larsen,
T. J., Yde, A., Stertr, K., Okstad, K. M.,
Nichols, J., Nygaard, T. A., Gao, Z., Manolas,
D., Kim, K., Yu, Q., Shi, W., Park, H.,
Vsquez-Rojas, A., Dubois, J., Kaufer, D.,
Thomassen, P., de Ruiter, M. J., Peeringa, J.
M., Zhiwen, H., von Waaden, H., 2012, Offshore Code Comparison Collaboration Continuation (OC4), Phase I Results of Coupled
Simulations of an Offshore Wind Turbine with
Jacket Support Structure Proceedings,
Twenty-Second International Offshore and
Polar Engineering Conference Rhodes, Greece,
June 1722, 2012
Previsic, M., Bedard, R., Hagerman, G.,
2004, EI2 EPRI Assessment: Offshore wave
energy conversion devices, Technical Report,
E2I EPRI WP-004-US-Rev1, Electricity Innovation Institute, USA.
Renzi E., Dias F., 2012, Resonant behaviour of an oscillating wave energy converter in
a channel, Journal of Fluid Mechanics,
Vol. 701, pp 482-510
Salvatore, F., Greco, L., 2008, Development and Assessment of Performance Prediction Tools for Wind and Tidal Turbines,
RINA Conference on Marine Renewable Energy, London, UK, November, 2008.
Saulnier, J. B., Clment, A. H., Falcao, A.
F. de O., Pontes, T., Prevosto, M., Ricci, P.,
2011, Wave groupiness and spectral bandwidth as relevant parameters for the performance assessment of wave energy converters,
Ocean Engineering, Vol. 38 (1), pp 130-147
723
724
Yang, L., Hals, J., Moan, T., 2010, Analysis of dynamic effects relevant for the wear
damage in hydraulic machines for wave energy
conversion, Ocean Engineering, Vol 337 (13),
pp 1089-1102
Yu, Y-H, Li, Y., 2012, RANS Simulation
of the heave performance of a two-body floating-point absorbed wave energy system,
Computers and Fluids, Vol. 73, pp 104-114
725
1. INTRODUCTION
726
727
2.1. General
Following the recommendation of the
Advisory Committee the committee members
decided
to
review
the
following
Recommendations and Guidelines:
2.4. Resistance
Current procedures adopted in various ice
tanks for resistance tests have been reviewed.
Model towing methods practiced by ice model
tanks are described. A method to determine the
maximum force applied at model towing through
the ice ridges is given. Description is given of the
methods used for correcting ice resistance data
with minor deviations from the design ice
thickness and strength.
728
(3)
where f is the ice flexural strength and a and
b are the weight coefficients selected based on
experience (for example, one can use values
a=b=0.5).
(1)
where HI is the ice thickness.
(5)
(2)
where
Froude
where RIT,level ice and RIT,presawn ice are the total ice
resistance in level ice and
pre-sawn ice,
respectively.
number
729
Ice resistance
(6)
speed
(8)
(7)
730
&
(9)
where C is the correction coefficient.
The value of the correction coefficient C is
determined based on analysis of model and fullscale tests data. For them C may be determined
as given in the equations (10 12) (HSVA
Report):
- for level ice
(Alekseev
(10)
(13)
To correct model test data one can use ice
resistance computations made with one or
731
4,0
3,5
3,0
R/RT, %
2,5
2,0
1,5
1,0
0,5
0,0
4
10
12
14
16
HR,
R/R850
1,5
1,0
Aalto University
Aker Arctic Technology Inc. (AATI)
Hamburgische Schiffbau-Versuchsanstalt
GmbH (HSVA)
Japan Marine United Corp. (JMUC)
Kyrylov State Research Center new ice
tank (KSRC)
Krylov State Research Center old ice tank
(KSRC)
National Maritime Research Institute
(NMRI)
National Research Council / Ocean,
Coastal
and
River
Engineering
(NRC/OCRE)
0,5
750
800
850
, kg/m3
900
950
b)
4.1.
3. QUESTIONNAIRE
A questionnaire for investigating the
inventory of the ice basins was developed by
732
4.3.
733
In August 2011 the Ice Committee was reinstalled with 10 members. For various reasons
only 7 were consistently contributing to the
activities. In December 2011 the committee
started work with a Kick-Off meeting in
Hamburg, followed by 4 more physical meetings
and 6 video conferences.
After completion of a questionnaire and
compilation of a new inventory list of the
existing ice model basins the committee
concentrated on review and editing of some ice
basin work related guidelines and procedures.
Two of them were finished and accepted by the
Advisory Committee for publication.
734
The Recommendations
include:
and
Guidelines
6.2.
735
736
8.
8.1.
LIST OF PUBLICATIONS
737
738
8.2.
Haase A. and Jochmann P., 2013. DYPIC Dynamic positioning in ice second phase of
model testing. Proceedings of the ASME
2013 32nd International Conference on
Ocean, Offshore and Arctic Engineering
OMAE2013, Nantes, France.
739
8.4.
740
741
742
8.6.
743
8.7.
744
745
8.8.
Metrikin, I., Borzov, A., Lubbad, R., and Lset,
S., 2012. Numerical simulation of a floater in
a broken-ice field: Part II - Comparative
study of physics engines. Proceedings of the
31st International Conference on Offshore
Mechanics
and
Arctic
Engineering
OMAE2012, Rio de Janeiro, Brazil.
746
747
http://www.ittc.info