ITTC - Recommended Procedures and Guidelines: Preparation, Conduct and Analysis of Speed/Power Trials
ITTC - Recommended Procedures and Guidelines: Preparation, Conduct and Analysis of Speed/Power Trials
ITTC - Recommended Procedures and Guidelines: Preparation, Conduct and Analysis of Speed/Power Trials
5-04
-01-01.1
Procedures and Guidelines Page 1 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Procedure
Table of Contents
The present Recommended Procedure con- Normally, such stipulated conditions are not
cerns the preparation and execution of speed- experienced during the actual trials. In practice,
power trials, as well as the method of analysing certain corrections for the environmental condi-
the results. It has been defined by the 27th and tions, such as water depth, wind, waves, current
the 28th ITTC Specialist Committee on the Per- and deviating ship draught from specified
formance of Ships in Service. In this work the draught have to be considered. For this purpose,
Committee took into account: not only the shaft power and ship’s speed are
measured, but also relevant ship data and envi-
• ITTC 7.5-04-01-01.1 &.2, 2014, ronmental conditions shall be recorded during
• ISO 19019, 2002, the speed-power trials.
• ISO 15016, 2015, In case it is physically impossible to meet the
• Sea Trial Analysis JIP, 2006; (Boom, 2008). conditions in these Guidelines, a practical docu-
mented approach mutually agreed among
The descriptions for the calculation methods Owner, Verifier and Shipbuilder can be allowed.
of the resistance increase due to wind and waves,
as well as guidelines for analysis and speed cor- The applicability of these Guidelines is lim-
rections are based on relevant research results ited to commercial ships of the displacement
and modified from ITTC 7.5-04-01-01.2/2005 type.
to meet the IMO EEDI requirements.
All trial procedures and measurements shall
• The purpose of this document is to define be conducted in such a way that the speed at
and specify: Contract power and the speed at EEDI Power
• the responsibility of each party involved, are derived within 0.1 knots and the shaft power
• the trial preparations, within 2%.
• the vessel’s condition,
• the limiting weather and sea conditions,
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 6 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
1. Prior to dock-out the oil distribution The S/P trial agenda is a document prepared
mechanism showing the Propeller Pitch by the Shipbuilder, outlining amongst others the
shall be checked for zero pitch; scope of a particular Speed/Power trial. This
2. Check zero pitch reading in the measure- document contains the procedures on how to
ment system against the mechanical conduct the trial and table(s) portraying the runs
reading in the oil distribution box; to be conducted. It outlines the particular re-
3. Determine the maximum ahead pitch, sponsibilities of the Trial Leader, Trial Team,
design pitch, and maximum astern pitch ship’s crew/ Shipbuilder, and the Owner’s rep-
and then adjust the ship indicators to re- resentative. The scope of the S/P trials shall be
flect the measurements. Determine the in line with this document.
corrections to account for changes in
pitch due to shaft compression as thrust Preferably before the sea trials start, but at
increases and temperature effects on the the latest when the trial area is reached and the
Propeller Pitch control rod. environmental conditions can be studied, agree-
4. Verify the weight of the propulsor and ment between the Trial Team, Shipyard and
hub from the manufacturer’s specifica- Ship Owner and Verifier shall be obtained con-
tions for making thrust measurement cerning the limits of wind forces, wave heights
corrections. and water depths up to which the trials shall be
performed. Agreement shall be obtained con-
An important deliverable of this stage is a cerning the methods used to correct the trial data.
document describing the test set-up including The measured data, analysis process and the re-
evidence of the calibrations that have been car- sults shall be transparent and open to the Trial
ried out. Team.
placement determination shall be conducted ei- The objective during the S/P trial is to minimize
ther by reading the internal draught measure- the number of influencing factors.
ment system or by evaluating all tank soundings.
Although there are correction methods for
Displacement shall be derived from the certain deviations from the contract condition,
Bonjean data or using quadratic equations with these methods are only valid up to certain limits.
hydrostatic data, taking into consideration the
hog/sag using the draught data (forward, aft and In order to arrive at reliable S/P trial results
at half length) and the density of the water. the boundary conditions shall not exceed the
values given in this chapter.
The ship shall be brought into a loading con-
dition that is as close as possible to contract con- 6.1 Location
dition and/or the condition at which model tests
have been carried out. The loading condition High wind and sea state in combination with
shall be confirmed at zero ship’s speed. a heading deviating from head waves and fol-
lowing waves, can require the use of excessive
5.2 Trim rudder deflections to maintain heading, and thus
cause excessive fluctuations in propeller shaft
The trim shall be maintained within very torque, shaft speed and ship’s speed.
narrow limits. For the even keel condition, the
trim shall be less than 0.1% of the length be- The S/P trial shall be conducted in a location
tween perpendiculars. For trimmed trial condi- where the environmental conditions are ex-
tions, the forward draught shall be within ± 0.1 pected to be constant and to have only the small-
m of the ship condition for which model test re- est possible impact on the vessel in order to
sults are available. avoid unexpected environmental effects in the
S/P trial results.
5.3 Hull &propeller This means that the speed trial range shall be
located in a sheltered area (i.e. limited wind,
The ship shall have clean hull and propel-
waves and current). Furthermore, the area shall
ler(s) for the sea trial. Hull roughness and ma-
be free from hindrance by small boats and com-
rine growth can increase the resistance of the
mercial traffic.
ship significantly but are not corrected for in S/P
trials. Therefore, it is recommended that the hull
and propeller(s) be carefully inspected before 6.2 Wind
the sea trial, and cleaned as needed and as per
coating manufacturer's recommendation. The During the S/P trial the wind speeds shall not
dates of last docking and hull and propeller be higher than:
cleaning are to be recorded in the S/P trials re-
port.
where
h water depth [m],
B ship’s breadth [m],
Figure 1. Limits for allowable wave height TM draught at midship [m],
𝑉𝑉s ship’s speed [m/s],
6.3 Sea state G acceleration of gravity [m/s2].
The total significant wave height H1/3, de- The value of water depth to be used for cor-
rived from the significant wave heights of local rection shall not be less than the larger value ob-
wind driven seas (wind waves) HW1/3 and swells tained from the following two formulae:
HS1/3 , by 2
V
h = 2 B ⋅ TM and h = 2 S (3)
2
H1/3 =�HW1/3 2
+HS1/3 (1) g
shall satisfy the following criteria: Furthermore, areas with significant varia-
tions in the bottom contours shall be avoided.
H1/3 ≤ 1.5𝑥𝑥 when the wave height is derived The actual water depth during each speed run
from visual observations, shall be read from the ship’s instruments and
documented in the trial log.
H1/3 ≤ 2.25𝑥𝑥 when the wave spectrum encoun-
tered during the S/P trials is measured,
6.5 Current
with 𝑥𝑥 = �𝐿𝐿PP /100, Ideally S/P trials shall be conducted in a lo-
cation where current speed and direction are es-
See section 7.6.5 for definition of “observations” sentially uniform throughout the trial area.
respectively “measurements” in this context.
In cases of current time history deviating
The above limits are illustrated in Figure 1. from the assumed parabolic / sinusoidal trend
1
The Beaufort scale is given in Appendix B.
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 12 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
and the change of the current speed within the Table 1 Primary parameters
timespan of one double run is more than 0,5 Acceptable measurement Unit
knots/hour*timespan, neither of the correction devices
Ship Track DGPS [Latitude,
methods in Appendix H are applicable. Areas Longitude]
where this may occur shall be avoided for S/P or [m]
trials. Speed over DGPS [Knots]
Ground
Shaft Torque Torsion meter with [kNm],
or shaft power calibrated permanent
7. TRIAL PROCEDURES torque sensor or strain
gauges.
Power calculated from
7.1 Parameters that shall be recorded torque and RPM [kW]
Shaft RPM Pick-up, optical sensor,
In this chapter, an overview is given of the ship revs counter [RPM]
Propeller Pitch Bridge replicator
parameters that influence the trial speed. All Time GPS Time [hh:mm:ss]
these parameters shall be measured as accu- Water depth Ship echo sounder + [m]
rately as possible and recorded. nautical charts
Ship heading Gyro compass, or [deg]
compass- DGPS
For this purpose, a split has been made be- Relative wind, Ship anemometer, [m/s],
tween primary and secondary parameters. For speed and dedicated trial [knots],
each of the parameters the preferable measure- direction anemometer [deg]
ment methods are given. Height, period and Wave measuring device [m],
direction of wind such as wave buoy, radar, [sec],
waves and swell or lidar. Observation by [deg]
multiple Marinners.
7.2 Primary parameters Bow acceleration Calibrated acceleration [m/s2]
(for wave corr gauges
The primary parameters to be measured dur- method G.1))
ing each run and the accepted measurement de- Date [YYYY-
MM-DD]
vices are given in Table 1.
Table 2 Secondary parameters
It is recommended to record the wave height,
wave direction and period, absolute wind speed Acceptable measurement Unit
devices
and direction at station(s) in the vicinity of the Date [YYYY-
speed trial site. MM-DD]
Seawater density Salinity sensor, [kg/m3]
Conductivity Density
7.3 Secondary parameters Temperature (CDT) sensor
Seawater temperature Thermometer, CDT sensor [ºC]
Air temperature Thermometer [ºC]
The secondary parameters listed in Table 2
Air pressure Barometer [hPa],
shall be measured and recorded at the trial site [mBar]
at least once during the S/P trial. Sea trial area Geographical position by [Latitude,
DGPS Longitude]
Draughts, fore, Physical observation and / [m]
amidships and aft at or calibrated draught
zero speed gauges
Displacement According to draught [metric
readings and water density tons]
Torsion meter zero Torsion meter with kNm
setting calibrated torque sensor or
strain gauges
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 13 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
For the analysis of the S/P trials, i.e. to in- 7.6 Scope and conduct of the measure-
clude the effect of the propeller loading in non- ments
ideal conditions on the propulsion efficiency
and rpm, it is required that the model tests data
include the results of propeller load variation 7.6.1 Ship’s track and speed over ground
measurements as defined in [used in ITTC rec-
The ship’s position and speed shall be meas-
ommended procedures 7.5-04-01-xx (2017)]" )
ured by a global positioning system such as
or Appendix D.3.
DGPS. The DGPS system shall operate in the
Based on ITTC recommendations, the model differential mode to ensure sufficient accuracy.
tests shall be conducted according to the follow- Position and speed shall be monitored and
ing criteria: stored continuously.
7.6.4 Water depth height, period and direction into relevant equa-
tion (see section 10.3.2).
Measuring the water depth can be done by
using the ship echo sounder. It is important that 7.6.6 Density and temperature
the echo sounder is calibrated before the speed
run in combination with the check of the water The local seawater temperature and density
depth given on the charts and that the vessel’s at the trial site shall be recorded to enable the
draught is taken into account. Continuous re- calculation of the ship’s displacement and cor-
cording of water depth is recommended. rections with regard to viscosity. The water tem-
perature shall be taken at sea water inlet level.
7.6.5 Waves Air temperature and pressure shall be measured
at the trial location using a calibrated thermom-
The wave spectrum can be derived either by eter and barometer.
measurements or by observations.
7.6.7 Current
Wave measurements
Current speed in the direction of the ship’s
Preferably, the spectrum of waves induced
heading shall be derived as part of the evaluation
by local wind and swell originating from remote
of each run, either using the ‘Mean of means’
wind, shall be measured during the S/P trials.
method (Appendix H-2) or the ‘Iterative’
The spectrum is derived from a spectral analysis
method (Appendix H-1). See also section 10.2.4.
of the measured wave elevation as a function of
time. For this purpose, wave buoys in the speed
trial area or ship born equipment such as wave 8. CONDUCT OF TRIAL
radar or lidar can be used. The wave measure-
ment equipment shall be calibrated and the ac- On the day of and during the S/P trial, a num-
curacy shall be validated and documented. ber of prerequisites shall be met in order to ar-
rive at reliable trial results. In this chapter, an
The directions of the waves and swells may overview is given of the minimum requirements.
be derived from visual observations. Measure-
ment of directional wave spectrum is preferable.
8.1 Initiation
Wave observations
Prior to the S/P trials, the weather forecast
In case the wave spectrum encountered dur- shall be studied.
ing the S/P trials is not measured, the wave
height, direction and period shall be derived Whenever possible, the runs at EEDI power
from visual observations by multiple experi- shall be conducted in daylight to enable a clear
enced mariners, including the Owner’s repre- visual observation of the wave conditions. For
sentative and the Verifier. In addition to the trials in which the encountered wave spectrum
wave observations, wave now- or hind- cast data and the wave direction (both wind waves and
provided by an experienced and independent swells) are derived by measurements, these runs
weather office may be used. The wave spectrum may also be conducted without daylight.
is then obtained by entering the observed wave
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 15 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
It is important to check that the engine plant 8.3 Run duration and timing
configuration during the S/P trial is consistent
with normal ship operations. The S/P trial duration shall be long enough
to accommodate a speed/power measurement
Prior to the S/P trials, the following actions within the required accuracy. The run duration
shall be taken at the vessel's zero speed through shall be the same for all speed runs with a mini-
the water: mum of ten (10) minutes. The speed runs for the
same power setting shall be evenly distributed
1. draught reading as described in section in time.
5.1 and calculation of displacement,
2. measurement of wind speed and direc- 8.4 Trial direction
tion,
3. zero setting of shaft torque meter, The speed runs shall preferably be carried
4. measurement of water temperature and out by heading into and following the dominant
density. wave or wind direction, depending on which
effects the ship’s speed most.
8.2 Trial trajectory
Consequently, once the heading for the
The S/P trial runs need to be conducted over speed run and the reciprocal heading for the re-
the same ground area. For each base course, turn run are fixed, the selected heading shall be
each speed run will be commenced (COMEX) maintained very precisely throughout the S/P
at the same place (within reason). trial. However, if the ’Mean of means’ method
is used for current correction, the trial direction
can be changed between each power setting ac-
cording to change of weather condition.
8.5 Steering
No fixed approach distance can be given. In • Two (2) double runs (at the same power set-
order to verify that the vessel reached the steady ting) below the EEDI/Contract power,
ship’s condition the measured values of shaft ro- • Two (2) double runs (at the same power set-
tation rate, shaft torque (if available) and ship’s ting) around the EEDI/Contract power,
speed at the control position shall be monitored. • Two (2) double runs (at the same power set-
When all three values are stable the ship's con- ting) above the EEDI/Contract power.
dition shall be deemed "steady".
Two (2) double runs compensate for the ef-
8.7 Power settings fect of current and second order current varia-
tions. In order to obtain sufficient accuracy, the
A minimum of three (3) different power set- time intervals between each run at the same
tings are required. These shall be adequately dis- power setting shall be more or less the same
tributed within the power range of 65% MCR (time interval deviation of 25% between single
and 100% MCR. runs is acceptable).
All S/P trials shall be carried out using dou- If the results of the S/P trials of the first ship
ble runs, i.e. each run shall be followed by a re- of a series are acceptable, sister ships may be
turn run in the exact opposite direction per- subjected to a reduced speed trial program. The
formed with the same engine settings. runs shall comprise at least:
8.8.1 ‘Iterative’ method • One (1) double run below EEDI / Contract
power,
When the current correction is carried out • One (1) double run around EEDI / Contract
using the ‘Iterative’ method (Appendix H.1), the power,
runs shall comprise at least: • One (1) double run above EEDI / Contract
power.
• One (1) double run below EEDI / Contract
power, Additional runs – sister ships
• Two (2) double runs (at the same power set-
ting) around EEDI / Contract power, If ‘Mean of means’ method is used for sister
• One (1) double run above EEDI / Contract ships and a current variations of above 0.2 knots
power. within a double run are encountered, one (1) ad-
ditional double run at that power setting shall be
The EEDI / Contract power runs shall be conducted.
conducted not as the first or the last power set-
ting in the trial sequence. 8.8.4 Additional runs due to limiting wave
height
8.8.2 ‘Mean of means’ method
For the first of a series or a sister ship at any
When the current correction is carried out power setting, when the wave height is around
using the ‘Mean of means’ method (Appendix the limiting conditions and significant wave-in-
H.2), the runs shall comprise at least: duced ship motions are observed then one (1)
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 17 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
additional double run at that power setting shall shall always be to record as many parameters as
be conducted. possible by means of the measurement computer
in order to increase the level of accuracy of the
8.9 Test sequence S/P trials.
9.1.2 Location
Computed parameters:
ΒWT True wind angle in earth system [deg]
VWT True wind speed [m/s]
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 19 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Figure 5.Sign convention for wave directions Details of the methods are given in the fol-
lowing chapters. For wave and wind corrections
The wave direction is defined as the direc- the methods depend on the level of information
tion relative to the ship's heading from which the which is available to the conducting party of the
wave fronts are approaching. speed/power sea trials. The analysis and correc-
tion method to be followed is prescribed below
Input parameters:
and summarized in Table 3.
ψ Heading of the ship [deg]
HW1/3 Significant wave height (wind waves) Evaluation
[m], For the evaluation the Direct Power Method
HS1/3 Significant wave height (swell) [m], in combination with the propulsive efficiency
α Angle between ship heading and wave correction based on load variation tests (refer to
direction relative to the bow; 0 means ITTC 7.5-02-03-01 (2017)) shall be used.
head waves [deg]
VG Measured ship’s speed over ground Wind Correction
[knots] In calculating resistance increase due to
wind, four methods can be used, depending on
whether there are wind tunnel measurements
10. ANALYSIS PROCEDURE available or not:
This section describes the methods to ana- • Wind resistance coefficients from model test
lyse the results of speed/power trials as con- are used (Appendix F.1).
ducted according to the previous sections. The
method to be used depending on situation and If CFD simulations are available:
available data is given in Table 3.
• Wind resistance coefficients from simula-
tions are used (Appendix F.2).
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 20 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
If wind tunnel measurements or simulations • Under the condition that heave and pitch
are not available: motions are small, the direct correction wave
method based on wave reflection prescribed
• Wind resistance coefficients from standard in Appendix G.1 or G.2 shall be used.
data set (Appendix F.3) • In case significant heave and pitch is ob-
served during the trials, the empirical formu-
or lation of the response function prescribed in
Appendix G.2, shall be used for the analysis.
• Regression formula by Fujiwara et al. (Ap- This empirical transfer function covers both
pendix F.4). the mean resistance increase due to wave re-
flection and the motion induced added re-
are(is) used. sistance.
Wave correction
Table 3. Evaluation method to be followed. The numbers identify the method by the chapters or Appendix in which the
methods are described.
Evaluation / Correction Method
Condition Air Temp. & Water Displ. &
Evaluation Waves Wind Current
Resistance Density Depth Trim
Heading Yes H.1
changed be-
tween power
no H.1 or H.2
settings
Load Yes D
Variation
Test No D
available
G.1 or Included in
heave No
No G.2 method 10.3.3 10.3.4 10.3.5
Ship geome- and pitch
Yes G.2
try
available to G.1 or
Yes
Verifier G.2,G.3
Full Seakeeping
G.4
Model Tests available
Dataset of Wind Tunnel Tests
E.1/E.2
wind re- /CFD
sistance coef- Data set E.3
ficients
No E.4
Available
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 21 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
with
PSms shaft power measured for each run
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 22 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
ηS shaft efficiency (0.99 for conventional VS ship’s speed through water [m/s], see
shaft) 10.2.4
PBms brake power measured for each run 𝜉𝜉𝑛𝑛 , 𝜉𝜉𝑣𝑣 overload factors derived from load var-
ηM transmission efficiency iation model test (Appendix D)
Δ𝑉𝑉 speed correction due to shallow water
The corrected delivered power PDid is ob- [m/s], determined according to 10.3.4.
tained as follows (under condition 𝑃𝑃Dms −
∆𝑅𝑅𝑉𝑉S The extended analysis in Appendix J,
> 0):
𝜂𝜂
Did which is included informatively, is useful for
model test correlation purposes since it in-
1
PDid = volves the full-scale wake fraction.
2
∆RVS
2
∆RVS ∆RVS 10.2.4 Correction of the measured ship’s
PDms − + PDms − + 4 PDms ξP speed due to the effect of current
ηDid ηDid ηDid
(7) The ship's speed through water (VS) is the
measured speed over the ground (VG) corrected
VS ship’s speed through water [m/s], see
for the current speed (VC) at each run, VS=VG-VC.
10.2.4
ηDid propulsion efficiency coefficient in The current correction can be achieved by
ideal condition, from model test. two (2) different methods: either the ‘Iterative’
ξP overload factor derived from load vari- method or the ‘Mean of means’ method. The
ation model test. details are given in Annex H.
ΔR resistance increase due to wind, waves
and temperature deviations [N] (eq. 4 ). a) ‘Iterative’ method
PDid is the power in ideal conditions, i.e. no Based on the assumption that the current
wind, waves or other disturbances. For shallow speed varies with a semi-diurnal period, a cur-
water a speed correction is applied according rent curve as a function of time is created. In
to 10.3.4. Deviations in displacement are cor- the same process a regression curve represent-
rected for according to 10.3.5. ing the relationship between the ship’s speed
through the water and corrected power is de-
The correction of the propeller frequency
termined. The current curve and the regression
of revolution is also carried out considering
curve are created in one process. The regres-
load variation effect (Appendix D). The cor-
sion curve has no relation with the
rected shaft rate nid is
speed/power curve from the tank tests.
nms
nid = (8) The analysis of the direct power method as
PDms − PDid ΔV
ξn + ξv +1 described in Appendix D shall be repeated af-
PDid VS ter the value of VS has been derived by the cur-
rent correction analysis.
with
nms measured propeller frequency of revo- b) ‘Mean of means’ method
lution [1/s],
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 23 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
The conversion method to be followed to The wind speed and direction are usually
convert the trial results for trial condition to re- measured by the on-board anemometer, posi-
sults for the contractual or stipulated condition tioned mostly in the radar mast on top of the
is given in Appendix I. bridge. Both wind speed and direction at this
location may be affected by the geometry of
10.3 Calculation methods for resistance the vessel in particular the shape of the super-
increase and other corrections structure and the wheel house.
vector for the run-set at vertical position of an- α angle between ship heading and com-
emometer is found. The wind speed as meas- ponent waves; 0 means heading waves,
ured by the anemometer shall be corrected for VS ship’s speed through the water,
the wind speed profile taking into account the E directional spectrum.
vertical position of the anemometer and the
reference height for the wind resistance coeffi- If the directional spectrum is measured at
cients (normally 10 m) according to Appendix sea trials by a sensor and the accuracy is con-
E.2. This averaged true wind vector is then firmed, the directional spectrum is available. If
used to recalculate the relative wind vector for the directional spectrum is not measured it is
each speed-run of the set. This procedure is ex- calculated by the following relation:
plained in detail in Appendix E.1.
E = Sη(ω)G(α) (11)
The wind resistance coefficient shall be
based on the method according to Appendix F. with
G angular distribution function.
10.3.2 Resistance increase due to the effects Sη frequency spectrum.
of waves
The standard form of the frequency spec-
The most reliable way to determine the de- trum and the angular distribution function are
crease of ship’s increase of resistance in waves assumed for the calculation.
is to carry out sea keeping tests in regular
waves of constant wave height, at different For ocean waves the modified Pierson-
wave lengths and directions and at various Moskowitz frequency spectrum of ITTC 1978
speeds, and according to ITTC 7.5-02-07-02.2. is used:
22 s Γ 2 ( s + 1) 1
G (α ) cos 2 s (α − θ m ) (15) RT0 = ρ 0 SVS 2CT0 (19)
2π Γ(2 s + 1) 2
π π where
for − ≤ α − θm ≤
2 2 CF frictional resistance coefficient for ac-
tual water temperature and salinity,
where CF0 frictional resistance coefficient for ref-
s directional spreading parameter, erence water temperature and salinity,
Γ Gamma function, ΔCF roughness allowance associated with
θm primary wave direction; 0 means head- Reynolds number for actual water tem-
perature and salinity,
ing waves.
ΔCF0 roughness allowance associated with
For wind waves and swells, RAW is calcu- Reynolds number for reference water
lated for each run with the relevant wave temperature and salinity,
height, period and direction. CT0 total resistance coefficient for refer-
ence water temperature and salinity,
The resistance increase due to waves shall RAS resistance increase due to deviation of
be determined by tank tests or formulae shown water temperature and water density
in Appendix G. [N],
RF frictional resistance for actual water
temperature and salt content [N],
10.3.3 Resistance increase due to water tem- RF0 frictional resistance for reference water
perature and salt content temperature and salt content [N],
RT0 total resistance for reference water tem-
Both water temperature and salt content,
perature and salt content [N],
affect the density of the sea water and thus the
S wetted surface area [m2],
ship resistance. Usually, speed trials are cor-
VS ship’s speed through the water [m/s],
rected to a sea water temperature of 15°C and
a density of 1026 kg/m³. The effects of water 𝜌𝜌S water density for actual water tempera-
ture and salt content [kg/m3],
temperature and density that differs from these
values are calculated as follows: ρ0 water density for reference water tem-
perature and salt content.
ρ C + ∆CF0 CF, CF0, ΔCF and ΔCF0 are derived according
=
RAS RT0 S − 1 − RF F0 − 1 (16)
ρ0 CF + ∆CF to ITTC Recommended Procedures 7.5-02-03-
01.4, latest version, using the same roughness
with kS for ideal and actual condition.
1
=RF ρS SVS 2 ( CF + ∆CF ) (17) 10.3.4 Correction of the ship performance
2 due to the effects of shallow water.
1
=RF0 ρ 0 SVS 2 ( CF0 + ∆CF0 ) (18) Within the restrictions on water depth stip-
2 ulated in section 6.4, the results of
speed/power trials in restricted water depth
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 26 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Preparation
Conduct
Measured data
VG, Pms, nms
Correction of ship’s
speed for current
by ‘MoM’ or mean
Correction of power
for resistance increase
Wave data
‘Iterative’ method No
Wind data for current correc-
tion?
Yes
Water
temperature,
Water density Correction of ship’s
speed for current
Correction of power
for resistance increase us-
ing revised ship’s speed
Correction of power
Displacement for displacement
Final performance
Deng Rui, Li Chao, Dong Guoxiang, etc. "Op- Japan Ship Research Association, "SR208:
timization of the superstructure and fore- New Speed Trial Analysis Method", Re-
castle fairing of a multi-purpose ship" port of the SR208 Committee, 1993.
OMAE 2017-61256.
JTTC, "A proposal for a standard method of
Fujiwara, T., Ueno, M. and Ikeda, Y.: "A New speed trial analysis", Bull. SNAJ, No.262,
Estimation Method of Wind Forces and 1944.
Moments acting on Ships on the basis of
Physical Component Models", J. JTTC, "A tentative guide for the operation of
JASNAOE, Vol.2, 2005. speed trials with large vessels", Bull. SNAJ,
No. 442, 1966.
Hansa Int. Maritiem Journal 150th Year, No. 4,
April 2013, Hansa-online.de/STA-JIP.pdf. Kaiser, M., Ph.D. (2016), Results of aerody-
namic model tests for Handy Size Bulk
ISO 15016;”Guidelines for the assessment of Carrier, Technical Report No. RH-2016/T-
speed and power performance by analysis 104E, CTO.S.A.
of speed trial data”, 2002.
Maruo H., "On the increase of the resistance of
ISO 19019; “Guide for Planning, Carrying out a ship in rough seas (2nd report)", J. SNAJ,
and Reporting Sea Trials”, 2002. Vol. 108, 1960.
sistance”, 29th Symposium on Naval Hy- SRAJ, "A Study of ship speed trials", No. 2
drodynamics, Gothenburg, Sweden, Au- Standardization Panel, , Res. Rep. No. 12R,
gust 2012. 1972.
Raven, H.C.: “A New Correction Procedure Taniguchi, K. & Tamura, K., "On a new
for Shallow-Water Effects in Ship Speed method of correction for wind resistance
Trials”, Proceedings of PRADS2016, Co- relating to the analysis of speed trial re-
penhagen, 2016. sults", 11th ITTC, 1966.
Sea Trial Analysis JIP; “Recommended Prac- Tsujimoto, M, Kuroda, M. and Sogihara, N.,
tice for Speed Trials”, 2006, Public docu- “Development of a calculation method for
ment from www.marin.nl. fuel consumption of ships in actual seas
with performance evaluation”, Proceedings
The Ship Testing and Trial Trip Committee of of OMAE 2013, OMAE2013-11297, 2013.
the Association of Ship Technical Societies
in Norway, "Standardization code for trials World Meteorological Organization: Manual
and testing of new ships", 2nd Edition, 1971. on Codes, International Codes, Volume I.1,
Part A-Alphanumeric Codes, WMO-No.
SNAME, "Code for Sea Trials", 1989. 306 (1995 edition).
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 32 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Hull condition
Last date of cleaning hull
Hull appendages and Rudder
Geometry
Type
Rate of Movement during speed trials
Wind fetch
Height of anemometer above
waterline
Transverse Projected area above the
waterline including superstructures at
trial draught
Lateral projected area above the If F.3 method is used
waterline including superstructures at
trial draught
Propeller(s)
Type (FPP/CPP)
Pitch (FPP)
Direction of rotation
Number of blades
Shaft(s)
G modulus
Diameter (inside)
Diameter (outside)
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 33 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Probable Probable
Descriptive wave wave
term Mean veloc- height* in height* in
m s-1 km h-1 m.p.h. Land Sea Coast
ity in knots metres feet
0 Calm <1 0-0.2 <1 <1 Calm; smoke rises vertically Sea like a mirror Calm - -
1 Light air 1-3 0.3-1.5 1-5 1-3 Direction of wind shown by Ripples with the appearance of Fishing smack just 0.1 ¼
smoke drift but not by wind scales are formed, but without has steerage way (0.1) (¼ )
vanes foam crests
2 Light 4-6 1.6-3.3 6-11 4-7 Wind felt on face; leaves Small wavelets, still short but of smacks which then 0.2 ½
breeze rustle; ordinary vanes more pronounced; crests have a travel at about1–2 (0.3) (1)
moved by wind glassy appearance and do not knots
break
3 Gentle 7-10 3.4-5.4 12-19 8-12 Leaves and small twigs in Large wavelets; crests begin to Smacks begin to ca- 0.6 2
breeze constant motion; wind ex- break; foam of glassy appear- reen and travel about (1) (3)
tends light flag ance; perhaps scattered white 3–4 knots
horses
4 Moderate 11-16 5.5-7.9 20-28 13-18 Raises dust and loose pa- Small waves, becoming longer; Good working 1 3½
breeze per; small branches are fairly frequent white horses breeze, smacks carry (1.5) (5)
moved all canvas with good
list
5 Fresh 17-21 8.0-10.7 29-38 19-24 Small trees in leaf begin to
Moderate waves, taking a more Smacks shorten sail 2 6 (8½)
breeze sway; crested wavelets pronounced long form; many (2.5)
form on inland waters white horses are formed(chance
of some spray)
6 Strong 22-27 10.8-13.8 39-49 25-31 Large branches in motion; Large waves begin to form; the Smacks have double 3 9½
breeze whistling heard in telegraph white foam crests are more ex- reef in mainsail; care (4) (13)
wires; umbrellas used with tensive everywhere (probably required when fishing
difficulty some spray)
7 Near gale 28-33 13.9-17.1 50-61 32-38 Whole trees in motion; in- Sea heaps up and white foam Smacks remain in 4 13½
convenience felt when from breaking waves begins to harbour and those at (5.5) (19)
walking against wind be blown in streaks along the di- sea lie to
rection of the wind
8 Gale 34-40 17.2-20.7 62-74 39-46 Breaks twigs off trees; gen- Moderately high waves of All smacks make for 5.5 18
erally impedes progress greater length; edges of crests harbour, if near (7.5) (25)
begin to break into the spindrift;
the foam is blown in well-
marked streaks along the direc-
tion of the wind
9 Strong 41-47 20.8-24.4 75-88 47-54 Slight structural damage High waves; dense streaks of - 7 23
gale occurs (chimney pots and foam along the direction of the (10) (32)
slates removed) wind; crests of waves begin to
topple, tumble and rollover;
spray may affect visibility
10 Storm 48-55 24.5-28.4 89-102 55-63 Seldom experienced inland; Very high waves with long over- - 9 29
trees uprooted; considera- hanging crests; the resulting (12.5) (41)
ble structural damage oc- foam, in great patches, is blown
curs in dense white streaks along the
direction of the wind; on the
whole, the surface of the sea
takes on a white appearance;
the tumbling of the sea becomes
heavy and shock- like; visibility
affected
11 Violent 55-63 28.5-32.6 103-117 64-72 Very rarely experienced; Exceptionally high waves (small - 11.5 37
storm accompanied by wide- and medium-sized ships might (16) (52)
spread damage be for a time lost to view behind
the waves); the sea is com-
pletely covered with long white
patches of foam lying along the
direction of the wind; every-
where the edges of the wave
crests are blown into froth; visi-
bility affected
12 Hurricane 64 and 32.7 and 118 and 73 and - The air is filled with foam and - 14 45
over over over over spray; sea completely white with (-) (-)
driving spray; visibility very seri-
ously affected
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 34 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
* This table is only intended as a guide to show roughly what may be expected in the open sea, remote from land. It
shall never be used in the reverse way; i.e., for logging or reporting the state of the sea. In enclosed waters, or when
near land, with an off-shore wind, wave heights will be smaller and the waves steeper. Figures in brackets indicate the
probable maximum height of waves. (ref World Meteorological Organization, 1995)
Notes:
(1)*These values refer to well-developed wind waves of the open sea. While priority shall be given to the descriptive
terms, these height values may be used for guidance by the observer when reporting the total state of agitation of these
resulting from various factors such as wind, swell, currents, angle between swell and wind, etc.
(2) The exact bounding height shall be assigned for the lower code figure; e.g. a height to f4 miscoded as 5.
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 35 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Forwa rd Speed
Run No. Ti me Hea di ng UKC Speed Di recti on Hei ght Di recti on Peri od Hei ght Di recti on Peri od Torque Power Revs Torque Power Revs
/ Return (SOG)
[-] [-] [F / R] [deg] [kn] [m] [kn] [m/s ] [deg] [m] [deg] [s ] [m] [deg] [s ] [kNm] [kW] [RPM] [kNm] [kW] [RPM]
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 36 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
η ∆R
D.1. Propulsive efficiency correction =
Dms
ξP +1 (D-3)
ηDid Rid
The ship’s propulsive efficiency is affected
by the added resistance. This has to be taken where
into account when correcting the power. ξP overload factor derived from load vari-
ation model test, according to ITTC
The delivered power corrected to ideal con- Recommended Procedure 7.5-02-03-
dition is derived by 01.4(2017)
Rid resistance in ideal condition
P=
Did PDms − ∆ P (D-1)
This leads to the expression for the cor-
with rected delivered power:
ΔP correction of delivered power due to
the increased resistance and the ∆ RVS PDms
PDid =
PDms − 1 − ξ (D-4)
changed propulsive efficiency ηDid PDid P
∆ RVS η 1
∆P = + PDms 1 − Dms (D-2) PDid =
ηDid ηDid 2
∆ RVS
with PDms − η + (D-5)
Did
PDms delivered power derived from shaft 2
power or break power measured on + P − ∆ RVS + 4 P ∆ RVS ξ
Dms η
board for each single run [W], Did
Dms
ηDid P
VS ship’s speed through the water [m/s],
which can be obtained by the ‘Iterative’
method or the ‘Mean of means’ method, D.2. Correction of shaft rotation rate – ef-
Δ𝑅𝑅 resistance increase from relative wind, fect of added resistance and of shal-
waves and deviation of water tempera- low water
ture and water density for each run. The
value is computed according to section With the PDid found as described above the
10.3 in these Guidelines, [N], correction on shaft rate is:
ηDid propulsive efficiency coefficient in
ideal condition obtained from standard= ∆n P −P ∆V
ξ n Dms Did + ξ V (D-6)
towing tank test and interpolated to the nid PDid VS
speed VS,
ηDms propulsive efficiency coefficient dur- where
ing sea trial.
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 37 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
∆=
n nms − nid (D-7)
with
nm measured shaft rate [1/s],
nid corrected shaft rate [1/s],
ξn, ξv overload factors derived from load var-
iation model test,
ΔV speed correction due to shallow water
[m/s], determined by equation (K-10)
and Pd-Vs curve from model test.
nms
nid = (D-8)
P −P ∆V
ξ n Dms Did + ξ V +1
PDid VS
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 38 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
2
VWT ( i ) cosψ WT ( i ) + VWT ( i+1) cosψ WT ( i+1)
2
'
VWT(i/i+1) = 2
Figure E-1 True wind vectors and relative wind vec- VWT ( i )sinψ WT ( i ) + VWT ( i+1)sinψ WT ( i +1)
+
tors.
2
The true wind velocity and direction at the (E-4)
vertical position of the anemometer are calcu-
lated by:
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 39 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
VWT ( i )sinψ WT ( i ) + VWT ( i+1)sinψ WT ( i+1) The true wind velocity VWT( i ) , true wind di-
ψ WT(i/i+1)
'
= tan −1
VWT ( i ) cosψ WT ( i ) + VWT ( i+1) cosψ WT ( i+1) rection ψ WT(i) , relative wind velocity VWR ( i )
for VWT ( i ) cosψ WT ( i ) + VWT ( i+1) cosψ WT ( i+1) ≥ 0 and relative wind direction ψ WR(i) can then be
ψ =
( i ) , ψ WT(i) , VWR ( i ) and ψ WR(i) .
' ' ' ' '
WT(i/i+1) replaced by VWT
VWT ( i )sinψ WT ( i ) + VWT ( i+1)sinψ WT ( i+1)
tan −1 + 180
VWT ( i ) cosψ WT ( i ) + VWT ( i+1) cosψ WT ( i+1)
E.2. Correction for the vertical position of
for VWT ( i ) cosψ WT ( i ) + VWT ( i+1) cosψ WT ( i+1) < 0 the anemometer
(E-5)
The difference between the vertical posi-
tion of the anemometer and the reference
(i) =
'
VWR
height is to be corrected by means of the wind
(E-6)
= '2
VWT ( i ) + VG(i) + 2VWT ( i )VG(i) cos ψ WT(i) − ψ (i)
2 ' '
( ) speed profile given by formula (E-8).
1
ψ WR(i)
'
= tan −1
V '
WT ( i )
sin ψ( '
WT(i)
− ψ (i) ) VWTref
Z 9
= VWT ref (E-8)
VG(i) + V '
WT ( i ) (
cos ψ '
WT(i)
− ψ (i) ) Za
for VG(i) + VWT
' '
(
( i ) cos ψ WT(i) − ψ (i) ≥ 0 ) where
ψ WR(i)
'
tan −1
V '
WT ( i )
sin ψ( '
WT(i)
− ψ (i) ) + 180
VWTref true wind velocity at the reference
VG(i) + V '
WT ( i ) (
cos ψ '
WT(i)
− ψ (i) ) height [m/s];
for VG(i) + V '
WT ( i )
cos ψ ( '
WT(i)
− ψ (i) < 0 ) VWT true wind velocity at the vertical posi-
(E-7) tion of the anemometer in [m/s];
Z ref reference height for the wind resistance
where: coefficients in [m];
'
VWT averaged true wind velocity at the ver-
Za vertical position of the anemometer in
[m].
tical position of the anemometer [m/s];
'
VWR corrected relative wind velocity at the The reference height for the wind re-
vertical position of the anemometer sistance coefficients, 𝑍𝑍ref is selected as the cor-
[m/s]; responding height for the wind resistance coef-
ψ WT
'
averaged true wind direction at the ver- ficient from wind tunnel tests (normally 10m).
tical position of the anemometer [de-
The relative wind velocity at the reference
grees];
height is calculated by:
ψ WR
'
corrected relative wind direction at the
+ VG2 + 2VWTrefVG cos (ψ WT − ψ )
vertical position of the anemometer =
VWRref 2
VWTref
[degrees];
(i ) run number. (E-9)
VWTref sin (ψ WT − ψ )
ψ WRref = tan −1
VG + VWTref cos (ψ WT − ψ )
for VG + VWTref cos (ψ WT − ψ )≥0
ψ WRref = (E-10)
VWTref sin (ψ WT − ψ )
tan −1 +
180
VG + VWTref cos (ψ WT − ψ )
for VG + VWTref cos (ψ WT − ψ )<0
where:
VWRref relative wind velocity at the reference
height [m/s];
VWTref true wind velocity at the reference
height in [m/s];
ψ WRref relative wind direction at the reference
height [degrees];
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 41 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
F.3. Data sets of wind resistance coefficients LC = Loading Condition WT = Wind tunnel
L = Laden CFD = CFD computations
B = Ballast
Data sets of the wind resistance coefficient A = Average
CX are available for certain ship types shown in
Table F-1. These are based on projected frontal
area and wind speed at 10 m reference height.
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 42 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
0.0 60 0.31
70 0.47
80 0.48
-0.5 90 0.29
100 0.18
110 0.18
120 0.29
-1.0 130 0.37
140 0.47
150 0.56
-1.5 160 0.61
0 30 60 90 120 150 180 170 0.65
angle of attack [°] 180 0.63
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 43 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
LNG CARRIER
1.50
Angle Cx [-]
Prismatic integrated of
Prismatic
attack
integrated
1.00 [°]
0 -1.02
10 -0.99
0.50 20 -0.93
30 -0.83
40 -0.67
50 -0.48
Cx [-]
0.00 60 -0.30
70 -0.14
80 -0.03
-0.50 90 0.05
100 0.11
110 0.24
120 0.41
-1.00
130 0.56
140 0.68
150 0.79
-1.50 160 0.88
0 30 60 90 120 150 180 170 0.92
angle of attack [°] 180 0.91
LNG CARRIER
1.50
Prismatic extended deck Angle Cx [-]
of Prismatic
Spherical attack extended Spherical
1.00 [°] deck
0 -0.81 -1.11
15 -0.79 -0.79
0.50 30 -0.93 -0.70
45 -0.55 -0.56
60 -0.27 -0.37
75 -0.12 -0.12
Cx [-]
0.00
90 -0.20 0.30
105 0.00 0.58
120 0.19 0.83
-0.50 135 0.80 0.61
150 1.06 0.79
165 1.03 0.96
-1.00 180 0.60 0.74
-1.50
0 30 60 90 120 150 180
angle of attack [°]
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 44 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
90 -0.13 -0.25
0.0 105 0.28 0.12
120 0.77 0.49
135 1.09 0.91
-0.5 150 1.18 1.04
165 1.19 1.08
180 0.94 0.93
-1.0
-1.5
0 30 60 90 120 150 180
angle of attack [°]
90 -0.28 -0.26
0.0 105 0.04 -0.04
120 0.41 0.36
135 0.81 0.69
-0.5 150 0.99 0.90
165 1.00 0.85
180 0.89 0.67
-1.0
-1.5
0 30 60 90 120 150 180
angle of attack [°]
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 45 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
CAR CARRIER
1.50
Angle Cx [-]
Average of
1.00 attack Average
[°]
0 -0.52
15 -0.50
0.50
30 -0.40
45 -0.20
Cx [-]
60 -0.07
0.00
75 -0.04
90 0.12
105 -0.28
-0.50
120 0.17
135 0.68
150 0.90
-1.00
165 0.86
180 0.76
-1.50
0 30 60 90 120 150 180
angle of attack [°]
CRUISE FERRY
1.50
Angle Cx [-]
Average of
attack Average
1.00 [°]
0 -0.69
15 -0.73
30 -0.69
0.50
45 -0.34
50 -0.24
60 -0.26
Cx [-]
0.00 75 0.00
85 0.19
90 0.04
-0.50 95 -0.03
100 -0.04
110 0.10
120 0.24
-1.00
140 0.57
155 0.84
160 0.82
-1.50 170 0.73
0 30 60 90 120 150 180 180 0.67
angle of attack [°]
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 46 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 47 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
GENERAL CARGO
1.50
Angle Cx [-]
Average of
attack Average
1.00
[°]
0 -0.60
10 -0.87
0.50 20 -1.00
30 -1.00
40 -0.88
Cx [-]
0.00 50 -0.85
60 -0.65
70 -0.42
-0.50 80 -0.27
90 -0.09
100 0.09
-1.00 110 0.49
120 0.84
140 1.39
-1.50 150 1.47
0 30 60 90 120 150 180 160 1.34
170 0.92
angle of attack [°] 180 0.82
Angle Cx [-]
HANDY SIZE BULK CARRIER of Heavy
attack Ballast, Ballast,
1.50 [o] No
No
Cranes
Cranes
Heavy Ballast, No Cranes 0 -0.68 -0.59
1.00 10 -0.67 -0.63
Ballast, No Cranes
20 -0.66 -0.56
30 -0.62 -0.55
0.50
40 -0.58 -0.50
50 -0.45 -0.36
Cx [-]
60 -0.31 -0.21
0.00
70 -0.13 -0.06
80 -0.10 0.01
-0.50 90 -0.13 0.00
100 -0.05 0.02
110 0.08 0.04
-1.00 120 0.21 0.09
130 0.39 0.28
140 0.57 0.44
-1.50 150 0.65 0.56
0 30 60 90 120 150 180
160 0.71 0.61
Angle of attack [°] 170 0.71 0.64
180 0.80 0.64
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 48 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
Angle Cx [-]
HANDY SIZE BULK CARRIER of Heavy
Ballast,
attack Ballast,
1.50 Cranes
Cranes
[o]
0 -0.75 -0.71
Heavy Ballast, Cranes 10 -0.91 -0.84
1.00
Ballast, Cranes 20 -0.94 -0.89
30 -0.87 -0.82
0.50 40 -0.80 -0.70
50 -0.59 -0.49
Cx [-]
60 -0.34 -0.23
0.00 70 -0.18 -0.12
80 0.02 0.01
90 0.11 -0.02
-0.50
100 0.10 0.02
110 0.22 0.12
-1.00 120 0.44 0.36
130 0.65 0.55
140 0.78 0.71
-1.50 150 0.81 0.76
0 30 60 90 120 150 180
160 0.85 0.79
Angle of attack [°] 170 0.82 0.80
180 0.70 0.68
Angle Cx [-]
MULTI-PURPOSE CARRIER of
With partly With
attack containers containers
1.50 [°]
0 -0.81 -0.84
With partly containers 15 -0.76 -0.93
1.00
With containers 30 -0.71 -0.96
45 -0.59 -0.76
60 -0.47 -0.52
0.50
75 -0.22 -0.33
90 0.04 -0.14
Cx [-]
CDA ψ = 90(deg.)
=
WR
for ψ WR = 90 (deg.)
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 53 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
CMC
HC
LOA B
midship
Formula G-1 estimates the resistance in- G.2. Empirical correction method with fre-
crease in head waves provided that heave and quency response function for ships
pitch motions are small. The application is re- with heave and pitch during the
stricted to waves in the bow sector, within +/- speed runs (STAWAVE-2)
45 deg. off the bow. For wave directions out-
side this sector no wave correction is applied. The empirical method STAWAVE-2
(Boom, 2013) has been developed by STA-JIP
1 B to approximate the transfer function of the
RAWL = ρs gH W1/3
2
B (G-1)
16 LBWL mean resistance increase in heading regular
waves by using the main parameters such as
where ship dimensions and speed, see Figure G-2.
For this purpose, an extensive database of sea
B beam of the ship [m]
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 55 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
keeping model test results for a large popula- 11.0 for ω < 1
tion of ships has been used to derive parametric b1 = (G-7)
transformation functions. −8.50 elsewhere
1
f (Lpp, B) RAWRL = ρ S gζ A 2 Bα1 (ω ) (G-9)
2
2 1 1/2 1/4 1/8
Wave length / Ship Length
wave length / ship length
π 2 I12 (1.5kTM )
α1 (ω ) = f1 (G-10)
π 2 I12 (1.5kTM ) + K12 (1.5kTM )
Figure G-2 Parametric transfer function of mean re-
sistance increase in regular waves.
0.769
V
=
This empirical transfer function 𝑅𝑅wave co- f1 0.692 S + 1.81CB (G-11)
T g
vers both the mean resistance increase due to M
wave reflection 𝑅𝑅AWRL and the motion induced
resistance 𝑅𝑅AWML . where:
CB block coefficient,
=Rwave RAWRL + RAWML (G2) VS ship’s speed in m/s
k yy non-dimensional radius of gyration in
where,
lateral direction,
L PP ship length between perpendiculars,
RAWML = 4 ρS gζ A 2 B 2 / LPP raw (ω ) (G-3) TM draught at midship,
I1 modified Bessel function of the first
with kind of order 1,
K1 modified Bessel function of the second
raw (ω ) = kind of order 1,
b
(
b1
)
ω exp 1 1 − ω a1Fr1.50 exp ( −3.50 Fr )
b1
d1
with the following restrictions:
∞ Rwave (ω;VS ) 1
RAWL = 2 ∫ Sη (ω )dω (G-12) for τ ≥
0 ζA 2
4
−∞
∞
m4 ) H ( m)
1
2
(G-14)
tank tests in short waves or empirical ( m + k0τ )2 ( m + k cos α )
formula dm
( m + k0τ )4 − m2 k0 2
Applying the theoretical method, the mean
resistance increase in regular waves Rwave is 1
calculated from the components of the mean for τ <
4
resistance increase based on Maruo's theory
RAWM and its correction term which primarily
is valid for short waves RAWR. (
RAWM= 4πρS − ∫ + ∫ + ∫
m3
−∞
m2
m4
∞
m1 ) H ( m)
1
2
(G-15)
=
R RAWM + RAWR (G-13) ( m + k0τ ) ( m + k cos α )
2
wave dm
( m + k0τ )4 − m2 k0 2
where
RAWM: mean resistance increase in regular with
waves based on Maruo's theory (Maruo,
ωEVS
1960), which is mainly induced by ship τ= (G-16)
motion. g
RAWR: mean resistance increase due to wave
reflection for correcting RAWM. ω2
k= (G-17)
RAWR should be calculated with high g
accuracy because the mean resistance
increase in short waves is predominant g
one. k0 = 2
(G-18)
VS
This theoretical method is valid for all ship
types with the following restrictions: ωE= ω + kVS cos α (G-19)
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 57 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
m1 =
(
k 0 1 − 2τ + 1 − 4τ ) (G-20)
αT effect of draught and encounter fre-
2 quency,
ζA wave amplitude.
m2 =
(
k 0 1 − 2τ − 1 − 4τ ) (G-21) with
2
π 2 I12 ( keTdeep )
m3 = −
(
k0 1 + 2τ + 1 + 4τ ) (G-22)
αT =
π 2 I12 ( keTdeep ) + K12 ( keTdeep )
(G-26)
2
(
k 1 + 2τ − 1 + 4τ
m4 = − 0
) (G-23)
k=
e k (1 + Ω cos α )
2
(G-27)
2
ωVS
Ω= (G-28)
H 1 ( m) = ∫ σ ( x ) e imx
dx (G-24) g
L
sin 2 (α + β w ) sin β w dl +
where
1 I
∫
g gravitational acceleration, Bf = (G-29)
B + sin 2 (α − β w ) sin β w dl
H1(m) function to be determined by the distri-
II
∫
bution of singularities which represents
periodical disturbance by the ship,
where
VS ship’s speed through the water,
α encounter angle of incident waves (0 I1 modified Bessel function of the first
deg. means head waves), kind of order 1,
ρS density of fluid, K1 modified Bessel function of the second
ω circular wave frequency, kind of order 1,
ωE circular wave frequency of encounter. k wave number,
Tdeep draught; for a trim condition Tdeep is the
The expression of RAWR is given by Tsu- deepest draught,
jimoto et al. (2013). The calculation method βw slope of the line element dl along the
introduces an experimental coefficient in short water line and domains of the integra-
waves into the calculation in terms of accuracy tion (I&II ) are shown in Figure G-3.
and takes into account the effect of the bow
shape above the water.
1 fore
=RAWR ρS gζ A 2 BBf αT (1 + CU Fr ) (G-25) aft II
X
2 G α I
βw
waves
where Y
B ship breadth,
Bf bluntness coefficient, Figure G-3 Coordinate system for wave reflection.
CU coefficient of advance speed,
Fr Froude number, The coefficient of the advance speed in
oblique waves CU(α) is calculated on the basis
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 58 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
of the empirical relation line shown in Figure of advance speed CU is determined by the least
G-4 2, which has been obtained by tank tests of square method through the origin against Fr;
various ship types following to the procedures see Figure G-6.
in the next paragraph. When CU(α=0) is ob-
70
tained by tank tests the relation used in oblique CU
60
waves is shifted parallel to the empirical rela-
tion line. This is illustrated in Figure G-5 for 50
30
70
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
60 Bf
50 Exp. in heading waves empirical relation
relation used in oblique waves
40
30 70
CU
20 60
10
50
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 40
|Bf | 30
20
2
The empirical relation line in Figure G-4 was obtained as
follows. CU is derived from the result of tank tests and RAWM, In calculating RAWM the strength of the singularity σ is
as formula (G-30). calculated by the formulation of slender body theory as for-
mula (D-30) and the singularity is concentrated at depth of
CVPTM.
1 Rwave
EXP
( Fr ) − RAWM ( Fr )
CU − 1 (G-30) 1 ∂
Fr 1 ρ gζ 2 BB α − iωE − VS {Z r ( x ) B( x )} (G-31)
σ=
2
s A f T
4π ∂x
with with
EXP B(x): sectional breadth,
Rwave : mean resistance increase in regular waves meas- CVP: vertical prismatic coefficient,
ured in the tank tests. t: time,
TM: draught at midship,
x: longitudinal coordinate,
Zr: vertical displacement relative to waves in steady mo-
tion.
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 59 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
αU=CUFr
(ii) B f (α= 0) ≥ B fc and B f (α= 0) ≥ B fs
F=
S 68 − 310 B f (α ) (G-36)
= U (α
FC C= 0) (G-37)
and B fs = 68 − CU (α = 0) .
Fr 58
where B fc =
310 310
ship’s speed in
the sea trial con-
ditions in this
range
G.4. Seakeeping model tests
Figure G-6 Relation between effect of advance speed Transfer functions of the resistance in-
(αU=CUFr) and Froude number Fr. crease in waves (Rwave) may be derived from
the tank tests in regular waves. The tank tests
The tank tests should be conducted for at have to be conducted for the specific vessel ge-
least three different Froude Numbers Fr. The ometry at the trial draughts and trim, and at
Fr should be selected such that the speeds dur- contractual draughts if required. A minimum
ing the sea trials lie between the lowest and the of two different ship’s speeds VS covering the
highest selected Fr. speed range tested in the speed/power trials
have to be tank tested.
When tank tests are not carried out, the co-
efficient of advance speed in head waves CU(α As trials are not always conducted in head
=0) is calculated by the following empirical re- seas and following seas, the tank tests should
lations, formulae (G-31), shown in Figure G-4. not only comprise head and following waves
The formulae are suitable for all ships. but also the relevant oblique wave conditions.
( ) +
CU (α ) = sgn B f (α ) ⋅ CU (| B f (α ) |)
(G-32)
A maximum interval of incident wave angle
shall be 30° for head to beam seas (0°-90°) but
may be larger for beam to following seas (90°-
with 180°).
The ship’s speed through the water at the If two double runs, i.e. four runs, are con-
corrected power for each run VS is recalculated ducted, the following relationship is derived
as the updated one from the formula (H-3), and for each run from formula (F.7).
the processes of Stage 2 and Stage 3 are then
repeated until ∑ (P(VS′)i − Pidi ) is minimized.
2 { }
VG1 = VS + VC,2 (t + 3∆t ) 2 − VC,1 (t + 3∆t ) + VC,0 (H8)
{ }
VG2 = VS − VC,2 (t + ∆t ) 2 − VC,1 (t + ∆t ) + VC,0 (H-9)
H.2. ‘Mean of means’ method
{
VG4 = VS − VC,2 (t − 3∆t ) 2 − VC,1 (t − 3∆t ) + VC,0 } t is the start time of the first speed run of
a power setting,
(H-11) Δt is half of the elapsed time between two
successive runs.
where:
VS is the ship’s speed through the water in The current speed is accounted for by sub-
knots, stituting the above four formulae from (H-8) to
VG1: is the measured ship’s speed over the (H-11) for the formula (H-6). The ship’s speed
ground on the first of four runs in knots, through the water is the ‘Mean of means’ of the
VG2 is the measured ship’s speed over the two double runs.
ground on the second of four runs in
knots, The propeller shaft speed and power shall
VG3 is the measured ship’s speed over the be averaged over the two runs of each double
ground on the third of four runs in knots, run and then over the other double runs for the
VG4 is the measured ship’s speed over the same power setting.
ground on the fourth of four runs in
knots,
, ,
Converged?
DPM
Figure I-1 Example of conversion from trial condition to other stipulated load condition at 75% MCR
ITTC – Recommended 7.5-04
-01-01.1
Procedures and Guidelines Page 64 of 76
Preparation, Conduct and Analysis of Effective Date Revision
Speed/Power Trials 2017 05
The delivered power corrected to ideal con- The propulsive efficiency coefficient, ηD, is
dition, PDid, is derived by calculated as follows:
P= PDms − ∆ P (J-1) 1− t
Did
ηD = ηOηR (J-3)
1 − wS
with
PDms delivered power derived from shaft where:
power or brake power measured on ηO propeller open water efficiency, which is
board for each single run [W], derived from propeller open water char-
ΔP correction of delivered power due to the acteristics of the actual propeller, con-
increased resistance and the changed sidering the propeller load,
propulsive efficiency [W]. ηR relative rotative efficiency,
t thrust deduction factor,
The correction of delivered power, ΔP, can wS full-scale wake fraction.
be written as:
The self-propulsion factors relative rotative
∆ RVS η efficiency, ηR, thrust deduction factor, t, and
∆P = + PDms 1 − Dms (J-2)
ηDid ηDid model wake fraction, wM, are obtained from
model self-propulsion tests. Between full-scale
with wake fraction, wS, and model wake fraction, wM,
it is generally assumed that there is the follow-
ΔR Resistance increase [N], which is de- ing relationship:
rived from the data measured during sea
trial, 1 − wS =(1 − wM ) ei (J-4)
VS ship’s speed through the water [m/s],
which can be obtained by the ‘Iterative’
with:
method or the ‘Mean of means’ method,
ηDid propulsive efficiency coefficient, ηD , in ei: scale correlation factor of the wake frac-
ideal condition, tion.
ηDms propulsive efficiency coefficient, ηD,
during sea trial. Calculation of ηDms
The propulsive efficiency coefficients, ηDid The propulsive efficiency coefficient in the
and ηDms obtained as outlined in the following trial condition, ηDms, is obtained as follows, by
sections. rewriting the formula (J-3):
1 − tms
ηDms = ηOmsηRms (J-5)
1 − wSms
where:
ITTC – Recommended 7.5-04-01-01.1
Procedures and Guidelines Page 65 of 76
ηOms propeller open water efficiency in the Rid resistance in the ideal condition [N],
trial condition, which also can be derived from the
ηRms relative rotative efficiency in the trial measured data during sea trial.
condition,
tms thrust deduction factor in the trial condi- The self-propulsion factors in the ideal con-
tion, dition, ηRid, tid and wMid, are obtained from stand-
wSms full-scale wake fraction in the trial con- ard self-propulsion test and interpolated to the
dition. speed, VS.
Each self-propulsion factor in the trial con- The deviations of the self-propulsion factors
dition, ηRms, tms and wMms, is obtained by adding ΔηR, Δt and ΔwM are considered as the functions
the deviation of each factor between the trial and of ΔR/Rid. The details of the functions are de-
the ideal condition, ΔηR, Δt and ΔwM, to each scribed in Appendix J.2.
factor for the ideal condition, ηRid, tid and wMid,
respectively, as follows: It is acceptable that ΔηR, Δt and ΔwM are set
to zero, because these values are negligibly
= ηRid + ∆ηR ( ∆ R / Rid )
ηRms (J-6) small in comparison with the deviation of ηO due
to the load variation effect.
tms= tid + ∆t ( ∆ R / Rid ) (J-7) Propeller efficiency ηO and full-scale wake
fraction wS are determined using propeller open
w=
Mms wMid + ∆wM ( ∆ R / Rid ) (J-8) water characteristics for the ship's fitted propel-
ler, i.e. curves of thrust coefficient, torque coef-
where: ficient and load factor, according to the follow-
ing procedure.
ηRms relative rotative efficiency in the trial
condition, Thrust coefficient, torque coefficient and
tms thrust deduction factor in the trial condi- load factor can be written as follows:
tion,
wMms model wake fraction in the trial condi- KT = aT J 2 + bT J + cT (J-9)
tion,
ηRid relative rotative efficiency in the ideal
KQ = aQ J 2 + bQ J + cQ (J-10)
condition,
tid thrust deduction factor in the ideal con-
dition, τP =
aT + bT / J + cT / J 2 (J-11)
wMid model wake fraction in the ideal condi-
tion, where:
ΔηR(ΔR/Rid) deviation of relative rotative effi- KT thrust coefficient,
ciency corresponding to ΔR/Rid, KQ torque coefficient,
Δt(ΔR/Rid) deviation of thrust deduction fac- τP load factor equal to KT/J2,
tor corresponding to ΔR/Rid, J propeller advance coefficient,
ΔwM(ΔR/Rid) deviation of wake fraction corre- aT, bT, cT factors for the thrust coefficient
sponding to ΔR/Rid, curve,
ΔR resistance increase [N], which is derived aQ, bQ, cQ factors for the torque coefficient
from the data measured during sea trial, curve.
ITTC – Recommended 7.5-04-01-01.1
Procedures and Guidelines Page 66 of 76
KTms = aT J ms
2
+ bT J ms + cT (J-14) τ Pms =
KTms
(J-18)
2
J ms
where:
Jms propeller advance coefficient in the trial where:
condition, Jms propeller advance coefficient in the trial
condition,
Therefore, the propeller efficiency in the KTms thrust coefficient in the trial condition.
trial condition, ηOms, is:
ITTC – Recommended 7.5-04-01-01.1
Procedures and Guidelines Page 67 of 76
Then, the total resistance in the trial condi- ηRid relative rotative efficiency in the ideal
tion, Rms, is: condition,
tid thrust deduction factor in the ideal con-
Rms = τ Pms (1 − tms )(1 − wSms ) ρ V D
2
S S
2 2
(J-19) dition,
wSid full-scale wake fraction in the ideal con-
where: dition.
τPms load factor in the trial condition, The self-propulsion factor in the ideal condi-
tms thrust deduction factor in the trial condi- tion, ηRid, tid and wMid, is obtained from standard
tion, self-propulsion test and interpolated to the speed,
wSms full-scale wake fraction in the trial con- VS.
dition,
ρS water density in kilograms per cubic me- The full-scale wake fraction in the ideal con-
tre, dition, wSid, is calculated by the following for-
VS ship's speed through the water in metres mula obtained by rewriting the formula (J-4):
per second,
D propeller diameter in metres. 1 − wSid =(1 − wMid ) ei (J-20)
The total resistance in the ideal condition, Rid, The scale correlation factor of wake fraction,
is obtained by subtracting the resistance increase, ei, included in the above formula is obtained us-
ΔR, from the total resistance in the trial condi- ing the full-scale and model wake fractions in
tion, Rms: the trial conditions:
=
Rid Rms − ∆ R (J-22) 1 − wSms
ei = (J-21)
1 − wMms
where:
ΔR resistance increase [N], which is derived where:
from the data measured during sea trial.
wMid model wake fraction in the ideal condi-
The total resistance in the ideal condition, Rid, tion,
is also used when the self-propulsion factor in wSms full-scale wake fraction in the trial con-
the trial condition are calculated with the formu- dition derived from the formula (J-17),
lae (J-6) to (J-8). wMms model wake fraction in the trial condi-
tion derived from the formula (J-8).
Calculation of ηDid
The load factor in the ideal condition, τPid, is
The propulsive efficiency coefficient in the calculated by the following formula:
ideal condition, ηDid, is obtained as follows, by
Rid
rewriting the formula (J-3): τ Pid = (J-23)
(1 − tid )(1 − wSid ) ρSVS2 D 2
2
1 − tid
ηDid = ηOidηRid (J-5)
1 − wSid where:
Rid resistance in the ideal condition [N],
where tid thrust deduction factor in the ideal con-
ηOid propeller open water efficiency in the dition,
ideal condition,
ITTC – Recommended 7.5-04-01-01.1
Procedures and Guidelines Page 68 of 76
J id KTid
ηOid = (J-27)
2π KQid
where:
Jid propeller advance coefficient in the ideal
condition,
KTid thrust coefficient in the ideal condition,
KQid torque coefficient in the ideal condition.
VS (1 − wSid )
nid = (J-28)
J id D
ITTC – Recommended 7.5-04-01-01.1
Procedures and Guidelines Page 69 of 76
𝑃𝑃Dms VS, ΔηR, Δt, ΔwM are obtained from the result of the above process
𝐾𝐾𝑄𝑄ms = 3 × 𝜂𝜂Rms DPM
2𝜋𝜋𝜌𝜌S 𝑛𝑛ms 𝐷𝐷5
Current correction Final “VS” is determined
2 𝐽𝐽id 𝐾𝐾𝑇𝑇id
𝐾𝐾𝑇𝑇id = 𝑎𝑎 𝑇𝑇 𝐽𝐽id + 𝑏𝑏𝑇𝑇 𝐽𝐽id + 𝑐𝑐𝑇𝑇 𝜂𝜂Oid =
2𝜋𝜋 𝐾𝐾𝑄𝑄id
2
𝐾𝐾𝑄𝑄id = 𝑎𝑎𝑄𝑄 𝐽𝐽id + 𝑏𝑏𝑄𝑄 𝐽𝐽id + 𝑐𝑐𝑄𝑄
1 − 𝑡𝑡id
𝜂𝜂Did = 𝜂𝜂Oid 𝜂𝜂Rid
1 − 𝑤𝑤Sid
∆𝑅𝑅𝑉𝑉𝑆𝑆 𝜂𝜂Dms
∆𝑃𝑃 = + 𝑃𝑃Dms �1 − �
𝜂𝜂Did 𝜂𝜂Did
2
J.2. Application of load variation test results ∆R ∆R
∆ηR ξ R
= +ζR (J-25)
In order to determine each component of Rid Rid
propulsive efficiency coefficient ηD, propeller
2
open water tests, resistance and self-propulsion ∆R ∆R
tests are carried out at trial draught and evalu-= ∆t ξ t +ζt (J-26)
ated according to the tank's normal procedures. Rid Rid
In addition, a self-propulsion test with load var- 2
iation effect may be carried out at the trial ∆R ∆R
draught and, as a minimum, one speed close to= ∆w ξ w +ζw (J-27)
Rid Rid
the predicted EEDI speed (75% MCR). This
speed shall be one of the speeds tested in the where:
normal self-propulsion test.
ΔηR deviation of the relative rotative effi-
The self-propulsion test with load variation ciency,
effect includes at least 4 self-propulsion test Δt deviation of the thrust deduction factor,
runs, each one at a different propeller shaft ΔwM deviation of the wake fraction,
speed while keeping the model's speed constant. ΔR resistance increase in newtons,
The propeller shaft speed is to be selected such Rid resistance in the ideal condition in new-
that: tons,
ρS
∆=
R ( FD − FM ) λ 3 (J-30)
ρM
where:
ΔR resistance increase [N],
Rid full scale resistance at the actual speed
from resistance test [N],
FX external tow force measured during load
variation test [N],
FD skin friction correction force same as in
the normal self-propulsion tests [N],
λ scale factor,
ρS water density in full scale [kg/m3],
ρM water density in the model test [kg/m3].
The power correction is computed according d) The viscous resistance coefficient thus be-
to the following steps: comes:
CV 1.06CF (1 + k ) + ∆CF
= (K-4)
( h)
1.79
TM
ΔRV = RVdeep 0.57 (K-6) Assuming constant Admiralty coefficient,
the resistance increase caused by additional
sinkage is estimated as:
K.3. Estimate additional sinkage
Rsink= (1 + δ∇ )
2/3
(K-9)
The increase of the dynamic sinkage due to
shallow water is found from the formula given
in Raven 2016: K.6. Correct measured power
d ( sinkage ) = The power corrected for sinkage effect and
for shallow-water effect on viscous resistance is:
∇ Frh 2 Frhd 2 (K-7)
1.46 −
LPP 2 1 − Frh 2
1 − Fr hd
2
PDshallow ΔRVVS
=
PDdeep − (K-10)
rsink ηDid
with a minimum of 0.0 and with