ITTC – Recommended 7.

5 – 02
03 – 01.4
Procedures Page 1 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

CONTENTS

1. PURPOSE OF PROCEDURE

2. DESCRIPTION OF PROCEDURE
2.1.1 Introduction for the Original 1978 ITTC Performance Prediction Method
for Single Screw Ships
2.1.2 Introduction for the 1978 ITTC Performance Prediction Method as Modified
in 1984 and 1987
2.2 Model Tests
2.3 Analysis of the Model Test Results
2.4 Full Scale Predictions
2.4.1 Total Resistance of Ship
2.4.2 Scale Effect Corrections for Propeller Characteristics.
2.4.3 Full Scale Wake and Operating Condition of Propeller
2.4.4 Model-Ship Correlation Factors
2.5 Analysis of Speed Trial Results
2.6 Input Data
2.7 Output Data
2.8 Test Example
3. PARAMETERS
3.1 Parameters to be Taken into Account
3.2 Recommendations of ITTC for Parameters
3.3 Input Data
4. VALIDATION
4.1 Uncertainty Analysis
4.2 Comparison With Full Scale Results
5. ITTC- 1978 PERFORMANCE PREDICTION METHOD (COMPUTER CODE)

COMMENTS OF PROPULSION COMMITTE OF 22nd ITTC

In its original form the ITTC 1978 Performance Prediction Method offers a valuable and rea-
sonably accurate prediction tool for reference purposes and conventional ships.
Edited by 22nd ITTC QS Group 1999 Approved
15th ITTC 1978 pp388 – 402
15th ITTC 1978, 17th ITTC 1984
17th ITTC 1984 pp326 - 333 th
and 18 ITTC 1987
18th ITTC 1987 pp266 - 273
Date Date
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 2 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

1978 ITTC Performance Prediction Method

1. PURPOSE OF PROCEDURE more convenient use of the program. These

extensions are summarized as follows.
The method predicts rate of revolution and
delivered power of a ship from model results. (1) Inclusion of prediction of propeller revo-
lutions on the basis of power identity.

2. DESCRIPTION OF PROCEDURE (2) Temporary measure for wTS > wTM

2.1.1 Introduction for the Original 1978 (3) Extension to twin screw ships
ITTC Performance Prediction Method
for Single Screw Ships (4) Addition of speed trial data

The method predicts rate of revolution and (5) Extension for the case of a stock propel-
delivered power of a ship from model results. ler in the self-propulsion test
The procedure used can be described as fol-
lows: (6) Adaptation to the input of the non-
dimensional resistance coefficient and
The viscous and the residuary resistance of the self-propulsion factors.
ship are calculated from the model resistance
tests assuming the form factor to be independ- In recent years, many member organizations
ent of scale and speed. have been asked by their customers for a gen-
eral description of the method, viz., model test
The ITTC standard predictions of rate of revo- and analysis of their results, calculation of full-
lutions and delivered power are obtained from- scale power and rate of propeller revolutions,
the full scale propeller characteristics. These and the model-ship correlation factors used.
characteristics have been determined by cor- Considering the above, it was decided to pre-
recting the model values for drag scale effects pare a user's manual of the 1978 ITTC method
according to a simple formula. Individual which includes all of the extensions and modi-
corrections then give the final predictions. fications made.

2.1.2 Introduction for the 1978 ITTC Per- 2.2 Model Tests
formance Prediction Method as
Modified in 1984 and 1987 Model tests required for a full scale com-
prise the resistance test, the self-propulsion test
The 1978 ITTC Method developed to pre- and the propeller open-water test.
dict the rate of propeller revolutions and deliv- In the resistance test the model is towed at
ered power of a single screw ship from the speeds giving the same Froude numbers as for
model test results has been extended during the the full scale ship, and the total resistance of
last two terms of the ITTC for a better and the model RTM is measured. The computer pro-
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 3 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

gram accepts either RTM in Newton, or in a non- ler and the model similar to the full-scale pro-
dimensional form of residuary resistance coef- peller should be tested in open water.
ficient CR assuming the form factor 1 + k. In
the latter case, the friction formula used can
then be either of the ITTC 1957, Hughes, 2.3 Analysis of the Model Test Results
Prandtl-Schlichting or Schönherr's formulae.
Resistance RTM measured in the resistance
The form factor 1 + k is usually determined tests is expressed in the non-dimensional form
from the resistance tests at low speed range or
by Prohaska’s plot of CFM against Fn4 R TM
C TM =
1
The ship model is not in general fitted with ρSV 2
2
bilge keels. In this case the total wetted surface
This is reduced to residual resistance coef-
area of them is recorded and their frictional
ficient CR by use of form factor k,
resistance is added in calculating the full-scale
resistance of the ship. viz.,

In the self-propulsion test the model is CR = CTM - CFM (1 + k)

towed at speeds giving the same Froude num-
bers as for the full-scale ship. Generally a tow-
ing force FD is applied to compensate for the Thrust, T, and torque Q, measured in the
difference between the model and the full-scale self-propulsion tests are expressed in the non-
resistance coefficient. dimensional forms

During the test, propeller thrust (TM), torque T Q

K TM = and K QM =
(OM) and rate of propeller rotation (nM) are ρD 4 n 2 ρD 5 n 2
measured.
With KTM as input data, JTM and KQTM are read
In many cases, stock propellers are used off from the model propeller characteristics,
which are selected in view of the similarity in and the wake fraction
diameter pitch and blade area to the full-scale
propeller. Then the diameter and the open- J TM D M
water characteristics of the stock propeller wTM = 1 −
V
have to be given as input data in the program.
In the open-water test, thrust, torque and rate of and the relative rotative efficiency
revolutions are measured, keeping the rate of K QTM
revolutions constant whilst the speed of ad- ηR =
K QM
vance is varied so that a loading range of the
propeller is examined. are calculated. V is model speed.
The thrust deduction is obtained from
In the case when a stock propeller is used in
the self-propulsion test, both the stock propel- T + FD − RC
t=
T
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 4 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

with  1

 k 3
∆C F = 105 S  − 0.6410 − 3
1   LWL  
FD = ρ M S M V M2 [C FM − (C FS + ∆C F )]  
2

where RC is the resistance corrected for differ- where the roughness kS=150.10-6 m and
ences in temperature between resistance and
self-propulsion tests: - CAA, is the air resistance
A
C AA = 0.001. T
(1 + k ).C FMC + C R S
RC = RTM
(1 + k ).C FM + C R
If the ship is fitted with bilge keels the total
resistance is as follows:
where CFMC is the frictional resistance coeffi-
cient at the temperature of the self-propulsion
S + S BK
test. C TS = [(1 + k )C FS + ∆C F ] + C R + C AA
S

2.4 Full Scale Predictions

2.4.2 Scale Effect Corrections for Propeller
2.4.1 Total Resistance of Ship Characteristics.

The total resistance coefficient of a ship The characteristics of the full scale propel-
without bilge keels is ler are calculated from the model characteris-
tics as follows
CTS =(1+k)CFS +CR+ ∆ CF +CAA
K TS = K TM − ∆K T
Where
K QS = K QM − ∆K Q
-k is the form factor determined from the where
resistance test
P c.Z
- CFS is the frictional coefficient of the ship ∆K T = − ∆C D .0.3.
according to the ITTC-1957 ship-model D D
correlation line
c.Z
∆K Q = − ∆C D .0.25.
- CR is the residual resistance calculated from D
the total and frictional coefficients of the
model in the resistance tests: The difference in drag coefficient ∆C D is
C R = C TM − (1 + k )C FM
∆C D = C DM − C DS
-. ∆C F is the roughness allowance
where
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 5 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

(1 − wTS )V S
nS = (r/s)
  J TS D
 t   0.04 5 
C DM = 2 1 + 2  −
 c   (R ) 16 (R ) 23 
 nco nco 
and - the delivered power:
−2.5 K QTS
 t  c  PDS = 2πρD 5 n S3 10 − 3 (kW)
C DS = 21 + 2 1.89 + 1.62. log  ηR
 c  kp 

- the thrust of the propeller:
In the formulae listed above c is the chord K
TS = T2 . J TS ρ .D 4 .n S2
2
length, t is the maximum thickness, P/D is the (N)
J
pitch ratio and Rnco is the local Reynolds num-
ber at x=0.75. The blade roughness kp is put
- the torque of the propeller:
kp=30.10-6 m. Rnco must not be lower than 2.105
K QTS
at the open-water test. QS = ρD 5 n S2 : (Nm)
ηR

2.4.3 Full Scale Wake and Operating Con- - the effective power:
dition of Propeller PE = C TS 1 / 2 ρ .V S3 .S .10 −3 (kW)

The full scale wake is calculated from the - the total efficiency:
model wake, wTM, and the thrust deduction, t: P
η D = DS
(1 + k )C FS + ∆C F PE
wTS = (t + 0.04 ) + (wTM − t − 0.04 )
(1 + k )C FM
- the hull efficiency:
1− t
where 0.04 is to take account of rudder effect. ηH =
The load of the full scale propeller is obtained 1 − wTS
from

KT S C TS 2.4.4 Model-Ship Correlation Factors

= .
J 2 2 D 2 (1 − t )(1 − wTS )2
Trial prediction of rate of revolutions and de-
livered power with CP - CN corrections
With this K T / J 2 as input value the full
scale advance coefficient JTS and the torque if CHOICE=0 the final trial predictions will be
coefficient KQTS are read off from the full scale calculated from
propeller characteristics and the following
quantities are calculated nT = CN.nS (r/s)
- the rate of revolutions:
for the rate of revolutions and
PDT = CP.PDS (kW)
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 6 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

2.5 Analysis of Speed Trial Results

for the delivered power.
The analysis of trials data is performed in a
Trial prediction with ∆CFC - ∆wC corrections way consistent with performance prediction but
starting PD and n backwards, i.e. from
If CHOICE=1 the final trial predictions are
calculated as follows: PD
KQ = .η RM .10³
2π . ρ .D 5 .n 3
KT S C TS + ∆C FC
= .
J 2 2 D 2 (1 − t )(1 − wTS + ∆wC )2 JS is obtained from the full-scale open-water
characteristics KQ ≈ JS then
With this KT/J² as input value, JTS and KQTS
are read off from the full scale propeller char- wT = 1 − J S .n.D / V
acteristics and
Further from KT ≈ JS characteristics
(1 − wTS + ∆wC )V S
nT = (r/s)
J TS .D T = K T . ρ .n ² D 4

K QTS T .(1 − t )
PDT = 2π . ρ .D 5 .n T3 . .10 − 3 (kW) CT =
η RM 1
. ρ .V ².S
2

Trial prediction with CNP correction Then we obtain

∆C FC = C T − C TS
If CHOICE = 2 the shaft rate of rotation is pre- ∆wC = wTS − wT
dicted on the basis of power identity as fol-
lows.
2.6 Input Data
 KQ  1000.C P .PDS
  =
 T 2π . ρ .D ²V S (1 − wTS )³
3 Input data sheets are given in ENCL.1
 J³

K Q0  KQ  2.7 Output Data

=   .η RM
J³  J T
- Output data I gives ITTC Standard Pre-
diction with CP = CN = 1.0, together with
n S = V S (1 − wTS ) / J TS .D
model and full scale propulsive coeffi-
cients (ENCL. 4).
nT = C NP n S - Output data II gives the final ship predic-
tion (ENCL. 5).
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 7 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

- Output data III gives the analysis of the propeller diameter D = 8.2m
speed trial results (ENCL. 6).
Calculations were carried out with the
ITTC Trial Prediction Test Program with:
2.8 Test Example
CP = 1.01
To illustrate the program a prediction was CN = 1.02
made for a hypothetical ship with the following
particulars: The input data were taken as shown in
length between ENCL. 1 and the printout of the input data and
perpendiculars Lpp = 251.5m results are given in ENCL. 4 - 6.
breadth B = 41.5m
draft T = 16.5m
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 8 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 9 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 10 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 11 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 12 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

3. PARAMETERS "F" format requires the data to be input with a

decimal point; the number can appear
3.1 Parameters to be Taken into Account anywhere in the field indicated.
"A" format indicates that alphanumeric char-
acters must be entered in the appropriate
Froude scaling law
card columns.
ship-model correlation line ,friction line
kinematic viscosity
mass density The card order of the data deck must fol-
low the order in which they are described
blockage
below.
form factor
propeller loading
hull roughness
Card No. 1 Identifications
see also 3.3 Input Data
Card Form CC Definition
column at Symbol
1- 8 A - Project No.
3.2 Recommendations of ITTC for Pa- 9-16 A - Ship model No
rameters Propeller model No.
see 4.9-03-03-01.1 Propulsion Test 17-24 A -

1987 p.263 In using the 1978 ITTC Method 25-32 F SCALE Scale ratio
it is recommended that the rudder(s) be fitted
in hull resistance experiments for barge type Card No. 2 Ship particulars
forms where inflow velocity is relatively
large. Card For- CC Definition
column mat Symbol
9-16 F LWL Length of waterline
17-24 F TF Draft, forward
3.3 Input Data
25-32 F TA Draft, aft
33-40 F B Breadth
All data are either non-dimensional or 41-48 F S Wetted surface, with-
given in SI-units. out bilge keels
49-56 F DISW Displacement
Every data card defines several parameters 157-64 F SBK Wetted surface of
which are required by the program; each of bilge keels
these parameters must be input according to a 65-72 F AT Transverse projected
specific format. area of ship above
waterline
"I" format means that the value is to be input 72-80 F C3 Form factor deter-
without a decimal point and packed to the mined at resistance
right of the specified field. tests
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 13 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

Card No. 3 Particulars of full scale Mean values of the trial correction figures,
Card For- CC Definition Cp and CN can be obtained from the trial test
column mat Symbol material of the individual institutions by run-
8- 8 I NOPROP Number of propellers
ning the ITTC Trial Prediction Test Program.
should be 1 since method
is valid only for single If an institution wishes to give predictions
screw ships with a certain margin the input CP-CN-values
15-16 I NPB Number of propeller must be somewhat higher than these mean
blades values.
17-24 F DP Diameter of propeller
25-32 F PD075 Pitch ratio at x=0.75
33-40 F CH075 Chord length of Propeller Cards Nos. 5-14 Result of resistance and self-
blade at x=0.75 propulsion tests and model propeller charac-
41-48 F TMO75 Maximum blade thick- teristics.
ness of propeller at
x=0.75 Card Format CC Definition
49-56 F RNCHM Reynolds number at column Symbol
open-water test based on 1- 8 F VS Ship speed in knots
chord length and local 9-16 F RTM Resistance of ship
velocity π .0.75

V = VA 1 + 


2
model
 J 
17-24 F THM Thrust of propeller
at x-0.75. 25-32 F QM Torque of propel-
Card No. 4 General ler:QM:100
Card For- CC Sym- Definition 33-40 F NM Rate of revolution
column mat bol 41-48 F FD Skin friction correc-
tion force
2.- 4 I NOJ Number of J-values in the 49-56 F ADVC Advance coefficient,.
open-water characteristics open water
(J ≤ NOJ ≤ 10) 57-64 F KT Thrust coefficient,
7- 8 I NOSP Number of speeds in the open water
self- propulsion tests 65-72 F KQ Torque coefficient,
(NOSPmax=10) open water
9-16 F RHOM Density of tank water
17-24 F RHOS Density of sea water
The J-margin in the open-water character-
25-30 F TEMM Temperature of resistance
test
istics must be large enough to cover the
31-36 F TEMP Temperature at self- model and full scale J-values with some mar-
propulsion test - gin.
36-41 F TEMS Temperature of sea water Input data sheets are given in ENCL. 1.
48-48 I CHOICE CHOICE=0 C − C
P N
trial corr.
CHOICE==1:
∆C FC − ∆wC trial corr.
49-56 F CP Trial correction for shaft
power.
57-64 F CN Trial correction for rpm
65-72 F DELT Trial correction for ∆C F
CFC
72-80 F DELTWC Trial correction for ∆w
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 14 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

(13th 1972 pp.155-180) Empirical Power

4. VALIDATION Prediction Factor ( 1+X )
4.1 Uncertainty Analysis 2) Propeller Dynamics Comparative Tests
not yet available (13th 1972 pp.445-446 )

3) Comparative Calculations with the ITTC

4.2 Comparison With Full Scale Results Trial Prediction Test Programme
(14th 1975 Vol.3 pp.548-553)
The data that led to t ITTC-78 method can
be found in the following ITTC proceedings:
4) Factors Affecting Model Ship Correlation
(17th 1984 Vol. 1, pp274-291)
1) Proposed Performance Prediction Factors
for Single Screw Ocean Going Ships
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 15 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

5. ITTC- 1978 PERFORMANCE PREDICTION METHOD (COMPUTER CODE)

C
C ****************************************************************************************************
C * *
C * 1978 ITTC PERFORMANCE PREDICTION METHOD FOR SINGLE SCREW *
C * SHIPS *
C * (REVISED 1983 TO INCLUDE TRIAL ANALYSIS AND TWIN SCREW SHIPS* *
C * *
C ****************************************************************************************************
C
C DECLARATIONS
C
COMMON /A/ FILE(2),MODELS(2), MODELP(2), LPP,LWL,TF,TA,B,S,
* SCALE,RNCHM,DISW,NOPROP,NPB,DP,PD075,CH075.
* TM075,C3,SBK,AT,CP,CN,DELCF,DELWC,KSI,KPI,
* RHOM,RHOS,TEMM,TEMP,TEMS,VS(10),RTM(10),THM(10),
* QM(10),NM(10),ADVC(10),KT(10),KQ(10),THD(10),
* FD(10),IC,NOJ,NOSP,PI
C
COMMON /B/ ETARM(10),ETAO(10),ETAH(10),ETAD(10),AWTM(10),
* AWTS(10),ACFM(10),ACTM(10),AVS(10),AVM(10),
* ATS(10),AQS(10),APDS(10),APE(10),APDT(10),
* ANS(10),ANT(10),BPDT(10),BNT(10),KTSJ2(10),
* KQS(10),KTS(10),ACTS(10)
DIMENSION FILE1(2),MODLS1(2),MODLP1(2)
C
REAL LPP, LWL, KS1, KS, KP1, KP, NM1, NM, KT, KQ, KTM, KQ0, JTM,
* KTSJ2, JTS, NS, KQTS, KTS, KQS, KQM
DATA TRIAL /‘TRIA‘/
500 FORMAT(6A4,F8.0)
501 FORMAT(10F8.0)
502 FORMAT(2I4,9F8.0)
503 FORMAT(2I4,2F8.0,3F6.0,I6,4F8.0)
504 FORMAT(9F8.0)
600 FORMAT(/5X,’NUMBER OF ADV,KT AND KQ POINTS =’,15/
* 5X,’NUMBER OF SPEEDS =’,15/
* 5X,’NUMBER OF SPEEDS OR ADVC POINTS >10’/)
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 16 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C CONSTANTS
C
G=9.81
PI=3.14159
KP1=30.0
KS1=150.0
KS=1.5E-4
KP=0.3E-4
C
C READ INPUT DATA
C
1000 CONTINUE
READ(5,500,END=999) FILE,MODELS,MODELP,SCALE
READ(5,501) LPP,LWL,TF,TA,B,S,DISW,SBK,AT,C3
READ(5,502) NOPROP,NPB,DP,PD075,CH075,TM075,RNCHM
READ(5,503) NOJ,NOSP,RHOM,RHOS,TEMM,TEMP,TEMS
* IC,CP,CN,DELCF,DELWC
NMAX=MAX0(NOJ,NOSP)

IF(FILE(1).EQ.TRIAL) GOTO 100

READ(5,504)(VS(I),RTM(I),THM(I),QM(I),NM(I),FD(I),
* ADVC(I),KT(I),KQ(I);I=1,NMAX)

C
C WRITE INPUT DATA
C
CALL OUTPUT(1)
C
C CHECK
C
IF(NOJ.LE.10.AND.NOSP.LE.10) GOTO 2
WRITE(6,600) NOJ.NOSP
GOTO 1000
2 CONTINUE
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 17 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C RECALCULATION OF INPUT DATA
C
DO 3 I=1,NOJ
KT(I)=KT(I)*0.1
KQ(I)=KQ(I)*0.01
....3 CONTINUE
DELCF=DELCF*0.001
RNCHM=RNCHM*100000.
VISCP=((0.585E-3*(TEMP-12.0)-0.03361)*(TEMP-12.0)+
* 1.2350)*1.0E-6
VISCM=((0.585E-3*(TEMM-12.0)-0.0361)*(TEMM-12.0)+
* 1.2350)*1.0E-6
VISCS=((0.659E-3*(TEMS-1.0)-0.05076)*(TEMS-1.0)+
* 1.7688)*1.0E-6
C
C CORRECTION OF PROPELLER CHARACTERISTICS
C
CDM=2.0*(1.0+2.0*TM075/CH075)*(0.044/RNCHM**0.16667-
* 5.0/RNCHM**0.66667)
CDS=2.0*(1.0+2.0*TM075/CH075)/(1.89+1.62*ALOG10(CH075
* /KP))**2.5
DCD=CDM-CDS
DKT=-0.3*DCD*PD075*CH075*NPB/DP
DKQ=0.25*DCD*CH075*NPB/DP
DO 4 I=1,NOJ
KTS(I)=KT(I)-DKT
KQS(I)=KQ(I)-DKQ
KTSJ2(I)=KTS(I)/ADVC(I)**2
4 CONTINUE
DO 5 I=1,NOSP
VS1=VS(I)*0.15444
VM1=VS1/SQRT(SCALE)
NM1=NM(I)
C
C
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 18 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C CALCULATE ROUGHNESS ALLOWANCE AND SHIP TOTAL RESISTANCE

C
RNLP=LWL*VM1/(VISCP*SCALE)
RNLM=LWL*VM1/(VISCM*SCALE)
RNLS=LWL*VS1/VISCS
CFMC=0.075/(ALOG10(RNLP)-2)**2
CFM=0.075/(ALOG10(RNLM)-2)**2
CFS=0.075/(ALOG10(RNLS)-2)**2
CTM=RTM(I)*SCALE**3/(0.5*RHOM*VS1**2*S)
CR=CTM-(1.0+C3)*CFM
RTMC=RTM(I)*(1.0+C3)*CFMC+CR)/((1.0+C3)*CFM+CR)
THD(I)=(THM(I)+FD(I)-RTMC)/THM(I)
DELCF=(105.0*(KS/LWL)**0.33333-0.64)*0.001
CAA=0.001*AT/S
CTS=((1.0+C3)*CFS*DELCF)*(S+SBK)/S+CR+CAA
C
C MODEL PROPULSIVE COEFFICIENTS
C
FNOP=NPROP
KTM=(THM(I)/FNOP)/(RHOM*(DP/SCALE)**4*NM1*NM1)
KQM=(QM(I)*0.01/FNOP)/(RHOM*(DP/SCALE)**5*NM1*NM1)
JTM=APOL(0,KT,ADVC,NOJ,KTM,IX)
KQ0=APOL(0,ADVC,KQ,NOJ,JTM,IX)
WTM=1.0-JTM*DP*NM1/(VM1*SCALE)
C
C FULL SCALE WAKE
C
IF(JRUDER) 6,5,6
5 WTS=(THD(I)+0.04)+(WTM-THD(I)-0.04)*((1.0+C3)*CFS+DELCF)/
* ((1.0+C3)*CFM)
GOTO 7
6 WTS=(THD(I) )+(WTM-THD(I) )*((1.0+C3)*CFS+DELCF)/
* ((1.0+C3)*CFM)
GOTO 7
7 IF(WTS.GT.WTM) WTS=WTM
ETARM(I)=KQ0/KQM
C
C SAVE AREAS
C
ACTM(I)=CTM
ACFM(I)=CFM
AWTM(I)=WTM
AWTS(I)=WTS
ACTS(I)=CTS
AVS(I)=VS1
AVM(I)=VM1
8 CONTINUE
C
C ITTC STANDARD PREDICTION
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 19 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
CALL IP
C
C RETURN FOR NEW INPUT
C
DO 20 I=1,2
FILE1(I)=FILE(I)
MODLS1(I)=MODELS(I)

20 MODELP1(I)=MODELP(I)
SCALE1=SCALE
GOTO 1000
C
100 CONTINUE
DO 110 I=1,2
FILE(I)=FILE1(I)
MODELS(I)=MODLS1(I)
110 MODELP(I)=MODLP1(I)
SCALE=SCALE1
C
CALL ANLSYS
C
C RETURN FOR NEW INPUT
C
C
GOTO 1000
999 STOP
END
C
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 20 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C
*****************************************************************************************************
***
C
C OUTPUT IS USED FOR PRINTING INPUT DATA AND RESULTS
C
C IOUT= 1 INPUT DATA IS PRINTED
C 2 RESULT PAGE 1
C 3 RESULT PAGE 2
C
C
*****************************************************************************************************
***
C
SUBROUTINE OUTPUT(IOUT)

COMMON /A/ FILE(2),MODELS(2),MODELP(2),LPP,LWL,TF,TA,B,S

* SCALE,RNCHM,DISW,NOPROP,NPB,DP,PD075,CH075,
* TM075,C3,SBK,AT,CP,CN,DELCFC,DELWC,KSI,KPI,
* RHOM,RHOS,TEMM,TEMP,TEMS,VS(10),RTM(10);THM(10),
* QM(10),NM(10),ADVC(10),KT(10),KQ(10),THD(10),
* FD(10),IC,NOJ,NOSP,PI
C
COMMON /B/ ETARM(10),ETA0(10),ETAH(10),ETAD(10),AWTM(10),
* AWTS(10),ACFM(10),ACTM(10),AVS(10),AVM(10),
* ATS(10),AQS(10),APDS(10),APE(10),APDT(10),
* ANS(10),ANT(10),BPDT(10),BNT(10),KTSJ2(10),
* KQS(10),KTS(10),ACTS(10)
C
REAL LPP,LWL,KS1,KS,KP1,KP,NM1,NM,KT,KQ,KTM,KQ0,JTM,
KTSJ2,JTS,NS,KQTS,KTS,KQS
DIMENSION TEXT (16)
DATA TEXT /’INPU’,’T DA’,’TA ‘,’ ‘,
* ‘OUTP’,’UT D’,’ATA ‘,’1 ‘,
* ‘OUTP’,’UT D’,’ATA..’,’2 ‘;
* `TRIA`,`L AN`,ÀLYS`,ÌS `/
600 FORMAT(‘1’,19X,’1978 ITTC PERFORMANCE PREDICTION’,10X,
* ‘ENCL:’/
C?? * 20X,’METHOD ‘,8X,
* ‘REPORT:’/20X,4A4/)
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 21 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

601 FORMAT(5X,’IDENTIFICATION :’,18X,’SHIP:’//

* 5X,‘PROJECT :’,2A4,
* 10X,’LENGTH PP :’,F8.2,’ (M)’/
* 5X,’SHIP MODEL’ :’,2A4,
* 10X,’LENGTH WL :’,F8.2,’ (M)’/
* 5X,’PROPELLER MODEL :’,2A4,
* 10X,’DRAFT FWD :’,F8.2,’ (M)’/
* 5X,’SCALE FACTOR :’,F8.2,
* 10X,’DRAFT AFT :’,F8.2,’ (M)’/
* 43X,’BREADTH :’,F8.2,’ (M)’/
* 5X,’PROPELLER:’,
* 28X,’WETTED SURFACE :’,F8.0,’ (M**2)’/
* 43X,’DISPLACEMENT :’,F8.0,’ (M**3)’)
602 FORMAT(5X,’NUMBER OF PROPELLERS:’,I8/
* 5X,’NUMBER OF BLADES :’,I8,
* 6X,’FRICTION COEFFICIENT CF’/
* 5X,’DIAMETER :’,F8.3,’ (M)’,
* 2X,’CALCULATED ACCORDING TO ITTC-57’/
* 5X,’PITCH RATIO 0.75R :’,F8.4,
* 6X,’FORM FACTOR :’,F6.3,’ (BASED ON ITTC-57)’/)
603 FORMAT(5X,’HULL ROUGHN.*10**6 :’,F6.1,’ (M)’,
* 2x,’BILGE KEEL AREA :’,F6.1,’ (M**2)’,
* 5X,’PROPELLER BLADE ROUGHN.*10**6:’,F6.1,’ (M)’,
* 2X,’PROJ.AREA ABOVE WL. :’,F6.1,’ (M**2)’/)
604 FORMAT(5X,’CHORD LENGTH OF PROP.BLADE AT X=0.75:’,
* F7.4,’ (M)’/
* 5X,’THICKNESS OF PROP.BLADE AT X=0.75:’,
* F7.4’ (M)’/)
605 FORMAT(5X,’DENSITY OF WATER (TANK ) :’F7.1,
* ‘ (KG/M**3)’/
* ’DENSITY OF WATER (SEA ) :’F7.1,
* ‘ (KG/M**3)’/
* 5X,’TEMP. OF WATER (RESISTANCE TEST) :’F7.2,
* ‘ (CENTIGRADES)’/
* 5X,’TEMP. OF WATER (SELF PROP. TEST) :’F7.2,
* ‘ (CENTIGRADES)’/
* 5X,’TEMP. OF WATER (SEA ) :’F7.2,
* ‘ (CENTIGRADES)’//
* 5X,’MODEL TEST RESULTS:’,
* 30X,’OPEN WATER CHARACT.;’/
* 54X,’RNC :’’F5.2,’*10**5’/)
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 22 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

606 FORMAT(5X,’SHIP RESIS- FRICT. THRUST TORQUE RATE OF ‘,

* 2X,’ADVANCE THRUST TORQUE’/
* 20X,’REVS. RATIO COEFF. COEFF.’/
* 5X,’KNOTS N N N NM RPS ’,
* 7X,’J 10*KT 100*KQ’/)
607 FORMAT(1X)
608 FORMAT(‘+’,3X,F5.1,1X,F7.1,1X,F7.2,2X,2F7.1,F9.2)
609 FORMAT(‘+’,49X,F10.3,F7.3,F8.3)
610 FORMAT(5X,’SHIP MODEL:’//
* 8X,’SPEED RES. COEFF. FRICT. COEFF. THRUST DED.’,
* 2X,’MEAN REL.ROT.’/
* 6X,’VS VM TOTAL’,32X, ‘WAKE EFFIC.’/
* 5X,’KNOTS M/S CTM*1000 CFM*1000’,8X,’TM’,
* 7X,’WTM ETARM’/)
611 FORMAT(4X,F5.1,F7.3,F8.3,6X,F7.3,7X,F7.3,3X,F7.3,F8.3)
612 FORMAT(/5x,’ITTC STANDARD PREDICTION CP=CN=1.0 :’//
* 5X,’SPEED EFF. POWER DELIV. POWER RSATE OF REVS’,
* 2X,’ THRUST TORQUE’/
* 6X,’VS’,7X,’PE’,10X,’PD’,12X,’N’,10X,’T’,8X,’Q’/
* 5X,’KNOTS’,5X,’KW’,10X,’KW’,11X,’RPS’,9X,’KN’,
* 6X,’KNM’/)
613 FORMAT(4X,F5.1,F10.0,3X,F9.0,4X,F9.3,3X,F9.0,F8.0)
614 (FORMAT(/5X,’SPEED TOT. EFF. PROP.EFF. HULL EFF. SHIP WAKE’,
* 3X,’OPEN WATER CHAR. FULL SCALE:’/
* 5X,’KNOTS ETAD ETA0 ETAH’,/X,’WTS’,
* 9X,’J 10*KT 100*KQ’/)
615 FORMAT(‘+’,3X,F5.1,F8.3,3(3X,F7.3))
616 FORMAT(‘+’,50X,3F7.3)
617 FORMAT(/5X,’SHIP DELIVERED POWER RATE OF REVS.’/
* 5X, ‘SPEED --------------------------- ---------------------‘/
* 5X,’KNOTS KW HP RPS RPM’/)
618 FORMAT(4X,F5.1,2X,2F8.0,3X,F7.3,F8.2)
619 FORMAT(/5X,’SHIP TRIALS PREDICTION CP=’,F7.3,’ CN=,F7.3)
620 FORMAT(/5X,’SHIP TRIALS PREDICTION DELCFC*1000=’,
* F6.3,’ DELCW=’,F6.3)
ITEX=ICUT*4-4
WRITE(6,600) (TEXT(ITEX+1),I=1,4)
WRITE(6,601) FILE,LPP,MODELS,LWL,MODELP,TF,SCALE,TA,B,S,DISW
WRITE(6,602) NOPROP,NPB,DP,PD075,C3
C
GOTO(10,20,30,40) , IOUT
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 23 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C INPUT DATA IS LISTED
C
10 CONTINUE
WRITE(6,603) KS1,SBK,KP1,AT
WRITE(6,604) CH075,TM075
WRITE(6,605) RHOM,RHOS,TEMM,TEMP,TEMS,RNCHM
WRITE(6,606)
NMAX=MAX0(NOJ,NOSP)
DO 1 I=1,NMAX
WRITE(6,607)
IF(I. LE. NOSP) WRITE(6,608) VS(I);RTM(I);FD(I),THM(I),
QM(I),NM(I)
IF(I. LE.NOJ) WRITE(6,609) ADVC(I),KT(I),KQ(I)
1 CONTINUE
RETURN

C
C RESULTS PAGE 1
C
20 CONTINUE
WRITE(6,610)
DO 21 I=1,NOSP
CFM=ACFM(I)*1000.0
CTM=ACTM(I)*1000.0
WRITE(6,611) VS(I),AVM(I),CTM,CFM,THD(I),AWTM(I),ETARM(I)
21 CONTINUE
WRITE(6,612)
DO 22 i=1,NOSP
WRITE(6,613) VS(I),APE(I),APDS(I),ANS(I),ATS(I),AQS(I)
22 CONTINUE
WRITE(6,614)
DO 23 i=1,NMAX
WRITE(6,607)
IF(I.LE.NOSP) WRITE(6,615) VS(I),ETAD(I),ETA0(I),ETAH(I);
AWTS(I)
XKTS=KTS(I)*10.0
XKQS=KQS(I)*100.0
IF(I.LE.NOSP) WRITE(6,616) ADVC(I),XKTS,XKQS
23 CONTINUE
RETURN
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 24 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C RESULTS PAGE 3
C
30 CONTINUE
DCFC=DELCFC*1000.0
IF(IC.EQ.1) WRITE(6,620)DCFC,DELWC
IF(IC.NE.1) WRITE (6,619) CP,CN
WRITE(6,617)
DO 31 I=1,NOSP
WRITE(6,618)VS(I),APDT(I),BPDT(I),ANT(I),BNT(I)
31 CONTINUE
....40 RETURN
END

C
C
*****************************************************************************************************
***
C
C IRAT= 0 INTERPOLATION WITH A 2:ND DEGREE POLYNOMIAL
C =1 INTERPOLATION WITH A RATIONAL FUNCTION OF 2:ND DEGREE
C X = ARGUMENT ARRAY
C Y = VALUE ARRAY
C N = NUMBER OF ARGUMENTS
C EX = ARGUMENT
C IFEL = ERROR RETURN CODE
C
C
*****************************************************************************************************
***
C
REAL FUNCTION APOL(IRAT,X,Y,N,EX,IFEL)
DIMENSION X(1),Y(1)
C
C CHECK NUMBER OF POINTS > 2
C
IFEL=0
IF(X(1).GT.X(N)) GOTO 2
IF(X(1).GT.EX.OR.X(N).LT.EX) GOTO 7
DO 1 I=1,N
L=1
IF(EX-X(I)) 4,4,1
1 CONTINUE
GOTO 4
2 CONTINUE
IF(X(1).LT.EX.OR.X(N).GT.EX) GOTO 7
DO 3 I=1,N
L=I
IF(EX-X(I)) 3,4,4
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 25 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

3 CONTINUE
4 CONTINUE
M=2
IF(L.EQ.1) M=1
IF(L.EQ.3) M=3
LM=L-M
X1=X(LM+1)
X2=X(LM+2)
X3=X(LM+3)
Y1=Y(LM+1)
Y2=Y(LM+2)
Y3=Y(LM+3)
C
C INTERPOL. 2:ND DEGREE POLYNOMIAL
C
X21=X2-X1
X31=X3-X1
X32=X3-X2
IF(IRAT.EQ.1) GOTO 6
C1=Y1
C2=(Y2-C1)/X21
C3=(Y3-C1-C2*X31)/(X31*X32)
APOL=C1+(EX-X1)*(C2+C3*(EX-X2))
RETURN
6 CONTINUE
C
C INTERPOL. RAT. FUNCTION
C
Y21=Y2*X2*X2-Y1*X1*X1
Y32=Y3*X3*X3-Y2*X2*X2
A0=(Y32-X32*Y21/X21)/(X32*X31)
B0=(Y21/X21-A0*(X1+X2)
C0=((Y1-A0)*X1-B0)*X1
APOL=(C0/EX+B0)/EX+A0
RETURN
7 CONTINUE
WRITE(6,8)
8 FORMAT(/5X,’INCREASE THE J-RANGE’)
STOP
END
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 26 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C ********************************************************************
C
C ITTC PREDICTIONS
C
C ********************************************************************
C
SUBROUTINE IP

COMMON /A/ FILE(2),MODELS(2),MODELP(2),LPP,LWL,TF,TA,B,S,

* SCALE,RNCHM,DISW,NOPROP,NPB,DP,PD075,CH075,
* TM075,C3,SBK,AT,CP,CN,DELCFC,DELWC,KSI,KPI,
* RHOM,RHOS,TEMM,TEMP,TEMS,VS(10),RTM(10),THM(10),
* QM(10),NM(10),ADVC(10),KT(10),KQ(10),THD(10),
* FD(10),IC,NOJ,NOSP,PI
C
COMMON /B/ ETARM(10),ETA0(10),ETAR(10),ETAD(10),AWTM(10),
* AWTS(10),ACFM(10),ACTM(10),AVS(10),AVM(10),
* ATS(10),AQS(10),APDS(10),APE(10),APDT(10),
* ANS(10),ANT(10),BPDT(10),BNT(10),KTSJ2(10),
* KQS(10),KTS(10),ACTS(10)
C
REAL LPP,LWL,KS1,KS,KPI,KP,NM1,NM,KT,KQ,KTM,KQD,JTM,
* KTSJ2,JTS,NS,KQTS,KTJT2,KQOS,KQS,KTS
DO 3 I=1,NOSP
VS1=AVS(I)
CTS=ACTS(I)
WTS=AWTS(I)
C
C CALCULATE THE FULL SCALE LOAD ADVANCE COEFF: AND
C TORQUE COEFF.
C
FNOP=NOPROP
KTJT2=S*CTS*0.5/((DP*(1.0-WTS))**2*(1.0-THD(I))) /FNOP
JTS=APOL(1,KTSJ2,ADVC,NOJ,KT,KTJT2,IX)
KQOS=APOL(0,ADVC,KQS,NOJ,JTS,IX)
C
C THE RATE OF REV. AND THE DELIVERED POWER
C
NS=(1.0-WTS)*VS1/(JTS*DP)
APDS(I)=2.0*PI*RHOS*DP**5*NS**3*KQOS/ETARM(I)*0.001
ANS(I)=NS
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 27 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C THE THRUST AND TORQUE OF THE PROPELLER
C
ATS(I)=KTJT2*JTS**2*RHOS*DP**4*NS*NS*0.001
AQS(I)=KQOS*RHOS*DP**5*NS*NS/ETARM(I)*0.001
C
C THE EFFECTIVE POWER, TOTAL AND HULL EFFICIENCY
C
APE(I)=CTS*0.5*RHOS*VS1**3*S*0.001
ETAD(I)=APE(I)/APDS(I)
ETAH(I)=(1.0-THD(I))/(1.0-WTS)
IF(IC.EQ.1) GOTO 1
C
IC1=IC-1
IF(IC1)10,11,12
C
C TRIAL PREDICTION WITH CP-CN CORRECTIONS (ITTC1978 ORIGINAL)
C
10 ANT(I)=CN*NS
BNT(I)=ANT(I)*60.0
APDT(I)=CP*APDS(I)
BPDT(I)=1.36*APDT(I)
GOTO 100
C
C TRIAL PREDICTION WITH CP-CN CORRECTIONS
C CN BASED ON POWER IDENTITY
C
12 APDT(I)=CP*APDS(I)
BPDT(I)=1.36*APDT(I)
KQJ3T=1000.0*APDT(I)/(2.0*PI*RHOS*DP**2) /FNOP
KQJ3T=KQJ3T/(VS1**3*(1.0-WTS)**3)
KQ0J3=KQJ3T*ETARM(I)
JTS=APOL(1,KQSJ3,ADVC,NOJ,KQ0J3,IX)
NS=(1.0-WTS)*VS1/(JTS*DP)
ANT(I)=CN*NS
BNT(I)=ANT(I)*60.0
GOTO 100
11 CONTINUE
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 28 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
C TRIAL PREDICTION WITH DELCF-DELWC CORRECTIONS
C
KTJT2=S*(CTS+DELCFC)/(2.0*(1.0-THD(I))*(DP*
* (1.0-(WTS-DELWC)))**2)
JTS=APOL(1,KTSJ2,ADVC,NOJ,KTJT2,IX)
KQOS=APOL(0,ADVC,KQS,NOJ,JTS,IX)
ANT(I)=(1.0-WTS+DELWC)*VS1/(JTS*DP)
BNT(I)=ANT(I)*60.0
APDT(I)=2.0*PI*RHOS*DP**5*ANT(I)**3*KQOS/ETARM(I)*0.001
BPDT(I)=1.36*APDT(I)
2 CONTINUE
ETAD(I)=KTJT2*JTS**3/(2.0*PI*KQOS)
3 CONTINUE
C
C WRITE OUTPUT
C
CALL OUTPUT(2)
CALL OUTPUT(3)
RETURN

SUBROUTINE ANLSYS
C
C***********************************************************************************************************
****
C * *
C * ANALYSIS ACCORD1NG TO 1978 ITTC PREDICTION METHOD *
C * *
C***********************************************************************************************************
****
C
C
DIMENSION VST(10),XNT(10),XPD(10),
* THDT(10),WTMT(10),WTST(10),ETART(10),CRWT(10),
* YNT(10),YPD(10),CPT(10),CNT(10),CNPT(10),ZNT(10)
* DCFT(10),WTSS(10),DWT(10),DCFM(10),DWM(I0),
* KQJ3(10)
C
COMMON /A/ FILE(2),MODELS(2),MODELP(2),LPP,LWL,TF,TA,B,S,
* SCALE,RNCHM,DISW,NOPROP,NPB,DP,PD075,CH075,
* TM075,C3,SBK,AT,CP,CN,DELCFC,DELWC,KS1,KP1,
* RHOM, RHOS,TEMM,TEMP,TEMS,VS(10),RTM(10),THM(10),
* QM(10),NM(10),ADVC(10),KT(10),KQ(10),THD(10),
* RA(10),IC,NOJ,NOSP,PI
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 29 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

C
COMMON /B/ ETARM(10), ETA0(10),ETAH(10),ETAD(10),AWTM(10),
* AWTS(10),ACFM(10),ACTM(10),AVS(10),AVM(10),
* ATS(10),AQS(10),APDS(10),APE(10),APDT(10),
* ANS(10),ANT(10),BPDT(10),BNT(10),KTSJ2(10),
* KQS(10),KTS(10),ACTS(10)

C
REAL LPP,LWL,KS1,KS,KP1,KP,NM1,NM,KT,KQ,KTM,KQ0,JTM,
* KTSJ2,JTS,NS,KQTS,KTJT2,KQOS,KTS,KQS,KQM,
* KQJ3,KQJ3T
C
C
DO 5 I = 1,NOJ
5 KQJ3(I) = KQS(I) /ADVC(I)**3
C
NOST=10

READ(5,510) (VST(I), I=1,NOST)

READ(5,510) (XNT(I), I=1,NOST)
READ(5,510) (XPD(I), .I=1,NOST)

510 FORMAT (10F8.0)

C
C COUNT NO. OF TRIAL RUNS
NOST = 0
DO 8 I = 1, 10
IF (VST(I).GT.0. ) NOST=NOST+1
8 CONTINUE
IF(XNT(1).GT.20.) GOTO 20
DO 10 I=1, NOST
XNT(I) = XNT(I)*60.0
10 XPD(I) = XPD(I)*1.36
20 CONTINUE
DO 50 I=1, NOST
VST1=VST(I)*1852.0/3600.0
CTST = APOL(0,AVS, ACTS, NOSP,VST1, IX)
THDT(I)= APOL(0,AVS, THD, NOSP,VST1, IX)
WTMT(I)= APOL(0,AVS, AWTM, NOSP,VST1, IX)
WTST(I)= APOL(0,AVS, AWTS, NOSP,VST1, IX)
ETART(I)= APOL(0,AVS, ETARM,NOSP,VST1, IX)
CF =APOL(0,AVS, ACFM, NOSP,VST1, IX)
CT =APOL(0,AVS, ACTM, NOSP,VST1, X)
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 30 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

CRWT(I)= CT - (1.0+C3)*CF
FNOP =NOPROP
KTJT2 =S*(CTST/FNOP )*0.5 / ((DP*(1.0-WTST(I)))**2*(1.0-THDT(I)))
JTS =APOL(1, KTSJ2, ADVC, NOJ, KTJT2, IX)
KQOS=APOL (0, ADVC, KQS, NOJ, JTS, IX)
NS=(1.0-WTST(I))*VST1/(JTS*DP)
PDS = 2.0*PI*RHOS*DP**5*NS**3*KQ0S/ETART(I)*0.001*FNOP
YNT(I)= NS*60.0
YPD(I) = PDS*1.36
CPT(I)= XPD(I)/YPD(I)
CNT(l)=XNT(I)/YNT(I)
PDT1 = XPD(I) /1.36
XNT1 = XNT(I) / 60.0
FKQ = PDT1*START(I)*1000.0 / (2.0*PI*RHOS*DP**5*XNT1**3) / FNOP
FJT = APOL(0,KQS,ADVC,NOJ,FKQ,IX)
FKT = APOL(0,ADVC, KTS,NOJ,FJT,IX)
KQJ3T=FKQ * (DP*XNT1)**3 / ((1-WTST(I))*VST1)**3
FJQ= APOL( 1,KQJ3,ADVC,NOJ,KQJ3T,IX)
ZNT(I)=(1.0 -WTST(I)) * VST1 / (FJQ*DP) * 60.0
CNPT(I)=XNT(I) / ZNT(I)
THS= FKT * RHOS * DP**4*XNT1**2
CTS=THS*(1.0 - THDT(I)) / (0.5*RHOS*VST1**2*S) * FNOP
DCFT(I)=(CTS - CTST)*1000.0
WTSS(I)= 1.0 - FJT*DP*XNT1/VST1
DWT(I) = WTST(I) - WTSS(I)
DWM(I) = WTMT(I) - WTSS(I)

C
C CALCULATION OF FRICTIONAL RESISTANCE ~COEFF. OF SHIP
C
T = TEMS
FNU = ((0.659E-3*(T-l.0)-0.05076)*(T-1)+1.7688)*1.0E-6
RNLS= ALOG10(LWL*VST1/FNU)
CFS = 0.075 / (RNLS-2.0)**2
C
DCFM(I) = CTS - (l.0+C3)*CFS - ( CRWT(I)+0.001*AT / S )*S / (S+SBK)
DCFM(I) = DCFM(I) * 1000.0
CRWT(I) = CRWT(I) * 1000.0
50 CONTINUE
C
CALL OUTPUT(4)
WRITE(6,600)
 ITTC – Recommended 7.5 – 02
03 – 01.4
Procedures Page 31 of 31
Performance, Propulsion
1978 ITTC Performance Prediction Effective Date Revision
1999 00
Method

600 FORMAT(' ',19X,'TRIAL ANALYSIS ACCORDING TO ITTC 1978 METHOD',///)

WRITE(6,610) ( VST(I), I=1, NOST)
610 FORMAT(5X.. ' SHIP SPEED - TRTAL',7(F10.2, 2X) /)
WRITE(6,620) ( XNT(I), I=1, NOST)
620 FORMAT(5X, ‘ PROP, RPM –TRTAL ',7(F10.2, 2X) /)
WRITE(6,630) ( XPD(I), I=1, NOST)
630 FORMAT(4X, 'DELIV.POWER-TRIAL ',7(F11.0,1X) //)
WRITE(6,640) ( YNT(I), I=1, NOST)
640 FORMAT(/5X, ‘ PROP. RPM -CN=1 ',7(F10.2,2X) /)
WRITE~(6,650) ( ~YPD(I), I=1,NOST)
650 FORMAT(4X, ' DELIV. POWER -CP =1',7(F11.0,1X) /)
WRITE(6,660) ( ZNT(I), I=1, NOST)
660 FORMAT(5X, ‘ PROP. RPM -CNP=1 ',7(F10.2,2X), //)
WRITE(6,670) ( CPT(I), I=1, NOST)
670 FORMAT(/5X, ‘ CP ‘,7(F10.3,2X) /)
WRITE(6,680) (CNT(I), I=1, NOST)
680 FORMAT(5X, ‘CN ‘,7(F10.3,2X) /)
WRITE(6,690) (CNPT(I), I=1,NOST)
690 FORMAT(5X, ‘CNP ',7(F10.3,2X) //)
WRITE(6,700) (DCFT(I), I=1,NOST)
700 FORMAT(/5X, ‘DCFC*1000 -CP=CN=1’,7(F10.3,2x) /)
WRITE(6,710) ( DWT(I), I=1, NOST)
710 FORMAT(5X, ' DWC CP=CN=1’,7(F10.3,2X) //)
WRITE(6,715) ( DCFM(I), I=1, NOST)
715 FORMAT(/5X, 'DCF *1000 ITTC-57’,7(F10.3,2x) /)
WRITE(6,717) ( DWM(I), I=1,NOST)
717 FORMAT(5X, ‘DW = WM-WTRIAL ',7(F10.3,2X) //)
WRITE(6,720) ( CRWT(I) ,I=1, NOST)
720 FORMAT(/5X, ‘ CR*1000 ‘,7(F10.3,2X) /)
WRITE (6,730) ( THDT(I), I=1, NOST)
730 FORMAT(5X, ‘ THDM ',7(F10.3,2X) /)
WRITE(6,740) ( WTMT(I), I=1, NOST)
740 FORMAT(5X, ’ WTM ',7(F10.3,2X) /)
WRITE(6,750) ( WTST(I), I=1, NOST)
750 FORMAT(5X, ‘ WTS CP=CN=1 ’,7(F10.3,2x) /)
WRITE(6,760) ( WTSS(I), I=1, NOST)
760 FORMAT(5X, ’ WTS TRIAL ’,7(F10.3,2X) /)
WRITE(6,770) ( ETART(I), I=1, NOST)
770 FORMAT(5X, ‘ ETARM ‘ ,7(F10.3,2X) /)
RETURN
END