Mr. J. H. Michell on the Wave–Resistance of a Ship.
Philosophical Magazine, 1898, vol.45, Ser.5, pp.106–123.
John Henry Michell (1863–1940),
Australian mathematician and
Professor of Mathematics at
the University of
Melbourne.

XI. The Wave-Resistance of a Ship.
By
J. H. Michell.
[*]       (Russian)


THE object of this paper is to give a general solution of the problem of the waves produced by a ship of given form moving with uniform velocity in an inviscid liquid, and to determine the consequent wave–resistance to the motion of the ship. The only assumption made as to the form of the ship is that the inclination of the tangent plane at any point of its surface to the vertical median plane is small. This condition is not satisfied near the bottom of the middle body of a modern ship, but it seems probable that this will not much affect either the waves produced or the resistance, for the waves arise rather from the parts at the bow and stern at which the tangent to the surface is inclined to the direction of the ship's motion, than from the approximately cylindrical middle body. The neglect of friction is probably of little consequence. The eddying water close to the side will no doubt slightly alter the virtual shape of the ship, but the change in the inclination of the virtual tangent plane, on which the wave–making depends, will, almost certainly, be very small. Further, the effect of viscosity in destroying the waves produced by the bow will modify to some extent the interaction of bow and stern waves; but, seeing that it is the waves of length comparable with that of the ship which chiefly give rise to the resistance, the effect must be small. The conclusion is, therefore, that the course followed by W. Froude, of considering frictional resistance and wave–resistance separately and adding the two, will probably give a close approximation to the truth.

A summary of the experimental work on this question, as well as a sketch of (he theoretical work of Russell, Rankine, and the two Froudes, is given in White's interesting 'Manual of Naval Architecture' (1894), chap. xi. Recent mathe­matics on the subject has been devoted chiefly to explaining, in a general way, the interesting wave–patterns observed, but exception must be made of the papers of Sir W. Thomson (Lord Kelvin), Phil.Mag. (1886–7), in which the critical speed of a canal boat, examined experimentally long before by Russell, was mathematically discussed. Reference may be made to Lamb's 'Hydrodynamics' (1895), chap, ix., and to Sir W. Thomson's 'Popular Lectures on Navigation' (1891), p. 450, for the discussion of wave–patterns.

None of these writers has, however, considered the waves produced by an actual ship, and the present paper is intended to supply the necessary investigation. The formula obtained for the wave–disturbance is rather complex, but that for the resistance is much simpler, as the most complicated term in the former represents a local disturbance not leading to any resistance in an inviscid liquid. There will therefore be no difficulty to those who have the necessary apparatus in making comparisons with experiment. As to general results, for deep water the theory leads to a resistance which increases with the velocity, in an oscillating manner, to a maximum and then decreases to zero as the velocity is indefinitely increased. That the resistance is an oscillating function of the velocity has been experimentally found by Mr. W.Froude and his son [†] , to whom also we owe the general explanation in terms of the interference of bow and stern waves. But the ultimate vanishing of the resistance has not, so far as I know, been anticipated. From general considerations it is clear that, so far as the wave–form is concerned, the effect of increasing the velocity is the same as that of decreasing the acceleration of gravity, and, if gravity vanishes, there is no propagation of waves; but this is not quite the theorem to be obtained.

Fig. 1

Take the vertical median plane of the ship as y = 0, and the surface of the undisturbed water as z = 0, the axis Ox being in the direction of motion of the ship and Oz vertically downwards. We may suppose the ship at rest and the water moving backwards with uniform velocity v apart from the wave–disturbance. The motion is assumed steady and the velocity potential written –vx + f. Since the inclination of the ship's surface to the plane y = 0 is everywhere small, f will be small, and we shall neglect the squares of the velocities due to f in comparison with their first powers. At the surface of the water let z be the depression at (x, y) below the mean level. Then

  (1)

is the kinematic surface condition, and

the equation of pressure, which, since

,   ....    (q.p.)

gives

and, therefore, with (1)

   ....   (2)

On account of the symmetry of the ship with respect to the median plane y = 0, we have d f/dy = 0 when y = 0, except over the ship, where, if h is the semi–breadth at (x, z),

    (say),   ...   (3)

and this condition is taken to hold at the plane y=0, instead of at the surface of the .ship, the justification being the same as that for equation (I), Finally d f/dz = 0, at the bottom, z = h, of the water.

We now consider the solution for f, in the part of the water where y is positive, with the given boundary conditions at z=0, z=h, y=0.

The typical term in the solution is

a cos n ( z–h ) cos ( mx+ a ) cos ( py+ b ),

where m2+n2+p2=0. Here m must be taken real as the water extends from x = – ¥ to x = + ¥; n and p maybe either real or imaginary, but if p is imaginary [ =ip' ] the last factor must take the form e –p'y.

This term satisfies df/dz=0 at z=h, and it also satisfies equation (2), if

n tan nh = –v 2m 2/g   ........   (4)

This equation has an infinite number of real roots and one pure imaginary root given by

n' tanh n'h = v 2m 2/g ,     [n = in'].

We shall see that the imaginary root is alone responsible for the wave–making resistance. As for p it is always imaginary for the real roots of n, and is so for the imaginary root if m > n'.

The condition (3) will now require the expansion of the given function f(x, z) in the form

åå amn cos n( z h ) cos ( mx + a ),

where the summation with respect to m will take the form of an integral.

Suppose at first the function periodic in x so that

f ( x +l, z ) = f ( xl, z ),

and put

,

where r is a positive integer.

By Fourier’s method

where A0n is to be halved as usual.

Since the functions cos n( zh ) are all conjugate, as is easily proved, from these we get

and   

where A00 is to be halved; and the coefficients of the terms given by the imaginary roots, here, as always below, are got by pitting n=in'; so that

and so for Brn'.

Hence the theorem

Now let l become infinite, and putting

p r/l=m
p/l=dm,

we get

In particular, supposes the depth of water infinite, we have then

n h = r p + e ,
h dn = p ,
tan nh = tane,
cos n( z – h )= (–)rcos ( nz– e ),
2n + (sin 2n h ) / h = 2n,
n
tan e = – km 2................... [k=v2/g],

Also

tanh n'h = 1,
n = k m2,

The result receives some confirmation from what has been observed with torpedo–boats at high speeds. It has been found that the total resistance varies as a power of the velocity which at first is nearly the second, but which, increasing to a maximum, ultimately Incomes less than the second [‡] . A very simple investigation, given below, shows that in shallow water, if we neglect all but the long waves, the wave–resistance varies ultimately as the first power of the velocity.

I may mention that somewhat similar work to that of the present paper gives a theory of the damping of the oscillations of ships due to wave–making. This I hope to give in a subsequent paper.

Substituting, we get

which is the theorem on which the rest of the present paper is build. It is curiously easy to give an á posteriory proof of the theorem. Using the value of e given above we find

Integrate the last two terms with respect to n, viz.

and the quadruple integral becomes

The former integral is f( x,z ), and the latter disappears with the triple integral in the given formula.

Considering now, for simplicity, the water infinitely deep, it appears at once that the required solution for f is

   ....   (5)

for this gives

when y = 0.

In this expression attention must be called to the factor

in the second integral. This form is not required to satisfy the boundary conditions formulated above; and it is evident the solution is to a certain extent indeterminate with those conditions, for we may superpose any system of free waves symmetrical with respect to y = 0 on a particular solution satisfying them. The form of the factor in question is chosen in order to make the elementary diverging waves trail aft; in other words, to satisfy the condition that the ship advances into still water.

Leaving the reduction of the integrals on one side, for the present, we proceed to calculate the wave–resistance ( R ).

Let dp be the increase of pressure duo to the wave–disturbance. Then

,

the double integral extending over the median plane of the ship. Now measuring from the undisturbed surface

p = p + g r z – ½ r q 2 + ½ r v 2,

and therefore

     ...   (q.p.)

so that

 .

Substituting the value of d f/dx, we see that the first and third integrals in the expression for f add nothing to the resistance because

,

on account of the odd factor sin m ( x – x ) , and hence

(6)

where l = mv2/g , and

If the ship is similarly formed at bow and stern I = 0, the origin being at midship.

We can now prove that the resistance vanishes when the velocity is infinite.

Observe that

   ....  (7)

where F( x ) is less than the greatest value of f( x,z ) for a given value of x; and, therefore, if we substitute a large number t, instead ∞ as the upper limit of λ,, the part of R neglected is of order not greater than

or v 2/t 2,

and this vanishes when v = ∞ if we take

.

In the part of R retained l g/v 2 is small throughout, so that we may expand the circular functions and write

and   

Now

,

the ship being of finite length and h = 0 at both ends. Hence, using the formula (7)


and   

The successive terms are of orders

v–2 log v 2 , v –10/3 , …

and all vanish when v = ¥. The resistance therefore ultimately vanishes. Of course this result is only proved for a ship which is very short in comparison with [§] the depth of the water.

We now proceed to the reduction of the integral which gives the resistance duo to two elements of the surface.

Consider two elementary areas s, s' at ( x,z ), ( x',z' ) on the side of the ship, and lot θ, θ' be the inclinations of the horizontal lines in these areas to the axis of x. The resistance due to these two elements is


or, say,    

Now, writing for the moment z = x + i y, and taking the integral

around the circuit enclosing the region x and y positive, we get

,

and, therefore [§] , realizing,

..... .... (8)

where k = log 2 – g = .11593 … .

From which

,
since    J0"(s) + J0'(s)/s + J0(s) = 0
J0'(s) = J1'(s)

[**] and so for Yo (s).

Now taking

and putting

  ... ... (9)

from equation (8).

From which H can be readily calculated by mechanical quadrature in the case in which 4·r is not large, and. this is the case for ordinary ships. Elaborate tables of J0 and J1 are now available, and tables [**] of Y0 and Y1 have been calculated by Mr. B.A. Smith, who has kindly prepared tables of kJ0 – Y0 and kJ1 – Y1, appended to the present paper.

We now have


where     ,

and the expression for the mutual resistance is

,

where

r = g ( z + z' ) / v2,
s = g ( x – x' ) / v2.

For elements at opposite ends of the ship s will in general be large compared with unity and with ; and in this case we can put

   ...   q.p.

and so for ( sm ) and then approximately

and the resistance is

Now if l is the length of the free wave which travels with the velocity of the ship

v2 = g l / 2 p ,

and the formula can be written

.

This gives a maximum resistance, approximately, when

x – x' = ( n+7/8 )·l    ..........    [n an integer]

and a maximum assistance when

x – x' = ( n+3/8 )·l.

These formula correspond to the interference of the bow and stern waves, which has been so skillfully discussed by Mr. R.E. Froude. When the two elements are on the same vertical cross–section of the ship, another form of reduction may be given. Putting xx' = 0, the integral to be considered is

.

Put

           l2 = ½ ( 1+ m),
so that   ,

and        l2 = ½ (m – 1 );
whence    

or, if m = cosh f ,

,

where K0 is the Bessel function, so indicated by Gray and Mathews (pp. 67, 90).

Hence

,

since K0' = K'1 , and the corresponding term in the resistance is

The functions K do not appear to have been calculated, but their general course is known. It will be sufficient at present to refer to Gray and Mathews, chap. vii.

As an illustration of the process of calculating the resistance of a given ship at any speed, we may consider one of simple analytical form which is fairly ship–shape and easily realisable. Experimental confirmation of the result was not practicable, and the matter must be left in the hands of those who have the necessary apparatus at command.

Let the surface of the ship be

y = ±c ( 1 + cos ax )( 1 + cos bz ),

between

 x=± p /a,
z = 0   and   p /b;

Fig. 2.

so that, for y positive,

Here I = 0, and

where k = v2/g; and therefore

which is best calculated by mechanical quadrature.

Suppose, for example, in foot–second units

v = 20 (velocity of ship),
2 p /a = 200 (length of ship),
p /b = 20 (depth below water–line),
8c = 32 (greatest breadth),

then the integral is found to be 620, and the resistance is  R=940 lbs. wt. about.

This seems to be about what one would expect from the experimental results available; but I know of no formula with which to compare it, and experiment alone can decide whether the theory has numerical value. Of course the method of successive approximations can be applied if necessary.

To examine the case of the ship in shallow water in which all but the long waves are neglected, we may proceed as follows: – We make the motion steady as before and use the same set of axes. The pressure–equation gives at the surface

, as before ;

while the equation of continuity for long waves gives

where h is the depth of the water, and this is

... (q.p.)

Differentiating with respect to x and substituting for f  from the pressure–equation, we obtain

or   

where c is the velocity of free long waves.

The ship being wall–sided, and extending to the bottom of the water, the kinematic equation over the ship is

,

which on differentiation with respect to x and use of the pressure–equation gives

and, of course,    

over the rest of the plane y = 0.

Now if v > c, i.e. if the velocity of the ship is greater than that of the free wave, the equation


is solved in the form .... (10)

where the boundary condition gives

,
or      ... .... (11)

The form of solution (10) is employed in order to make the diverging waves trail aft.

The disturbance therefore consists of two bands at an angle with the line of the ship’s motion, the front of each band being a hump above the mean level and its back part a hollow, which is similar to the hump if the ship is similarly shaped fore and aft.

The resistance (R) is given by

,

so that it is infinite when the velocity of the ship is equal to that of the free wave, and ultimately varies as the velocity.

If v < c , the differential equation for z takes the potential form

 ,

putting .

The solution is now

,

where r' 2 = ( x' – x )2 + y' 2, and there is no wave–resistance.

If we add the solutions for an equally–spaced infinite number of ships moving abreast, we get the case of a ship in the centre of a canal.

In the Table appended there may be an error of 1 in the last place or possibly of 2 in the values for x > 3 or 4.


Melbourne University, August 9, 1897.

Tables of k Jo – Y o and k J1 – Y1. By Mr. B. A. Smith, M.C.E.

    x.      k Jo(x)–Yo(x) k J1(x)–Y1(x)    x.     k Jo(x)–Yo(x)k J1(x)–Y1(x)
.00.41.92432.7384
.014.7200100.0261.42.89722.6822
.024.027450.0453.43.87062.6286
.033.621533.3951.44.84462.5773
.043.333125.0767.45.81902.5282
.053.109120.0903.46.79402.4813
.062.925816.7695.47.76942.4362
.072.770514.4002.48.74532.3929
.082.635912 6255.49.72162.3514
.092.516811.2470.50.69832.3111
.102.409910.1457.51.67532 2729
.112.31339.2459.52.65282.2357
.122.22458.4971.53.63062.1999
.132.14287.8645.54.60882.1653
.142.06707.3230.55.58732.1319
.151.99616.8545.56.56612.0995
.161.92976.4450.57.54532 0681
.171.86716.0843.58.52482.0377
.181.80795.7642.59.50462 0083
.191.75175.4780.60.48161.9798
.201.69825.2209.61.46501.9521
.211.64724 9888.62.44561.9251
.221.59834.7779.63.42641.8988
.231.55154.5855.04.40761.8732
.241.50664.4091.65.38901.8483
.251 .46334.2476.66.37071.8241
.261.42164 0983.67.35251.8005
.271.38133.9603.68.33461.7775
.281.34243.8323.69.31691.7550
.291.30463.7131.70.29951.7329
.301.26803.6020.71.28231.7114
.311.23263.4982.72.26531.6904
.321.19813.4007.73.24851.6699
.331.16453.3091.74.23191.6496
.341.13193.2233.75.21551.6299
.351.10003.1423.76.19931.6105
.361.06903.0656.77.18331.5914
.371.03872.9932.78.16751.5728
.381.00922.9245.79.15181.5544
.39.98032.8593.80.13631.5365
.40.95192.7973.81.12111.5188

Table (continued).

    x.      k Jo(x)–Yo(x) k J1(x)–Y1(x)    x.     k Jo(x)–Yo(x)k J1(x)–Y1(x)
.82.10601.50154.8.4277–.3354
.83.09101.49134.9.4588–.2817
.84.07631.46745.0.4846–.2323
.85.06171.45095.1.5052–.1786
.86.04731.43465.2.5203–.1214
.87.03301.41855.3.5300–.0690
.88.01891.40265.4.5343.–.0159
.89.00501.38705.5.5332+.0373
.90–.00881.37155,6.5269.0892
.91–.02241.35635.7.5155.1393
.92–.03591.31135.8.4991.1873
.93–04931.32615.9.4782.2326
.94–.06251.31176.0.4527.2749
.95–.07561.29726.1.4232.3139
.96–.08851.28296.2.3900.3492
.97–.10121.26886.3.3535.3805
.98–.11381.25486.4.3140.4077
.99–.12631.21086.5.2721.4306
1.0–.13861.22716.6.2281.4488
1.1–.25481.09666.7.1826.4625
1.2–.3583.97576.8.1358.4715
1.3–.4501.86166.9.0881.4758
1.4–.5304.75267.0..0408.4755
1.5–.6008..61777.1–.0065.4704
1.6–.6604.51607.2–.0531.4609
1.7–.7101.41737.3–.0986.4471
1.8–.7500.35137.4–.1424..4290
1.9–.7805.25837.5–.1842.4071
2.0–.8017.16827.6–.2237.3814
2.1–.8142.08127.7–.2604.3524
2.2–.8180–.00237.8–.2910.3203
2.3–.8138–.08217.9–.3243.2854
2.4–.8017–.15808.0–.3510.2484
2.5–.7824–.22938.1–.3739.2093
2.6–.7560–.29598.2–.3928.1684
2.7–.7234–.35768.3–.4076.1206
2.8–.6848–.41398.4–.4182.0840
2.9–.6407–.46198.5–.4214.0408
3.0–.5920–.51008.6–.4264–.0020
3.1–.5390–.54928.7–.4242–.0438
3.2–.4823–.58238.8–.4177–.0851
3.3–.4227–.60928.9–.4071–.1251
3.4–.3606–.62999.0–.3926–.1637
3.5–.2970–.64449.1–.3744–.2001
3.6–.2320–.65259.2–.3527–.2343
3.7–.1666–.65469.3–.3277–.2657
3.8–.1014–.65059.4–.2996–.2940
3.9–.0368–.64069.5–.2688–.3194
4.0+.0266–.62519.6–.2359–.3412
4.1+.0881–.60429.7–.2012–.3594
4.2.1472–.57819.8–.1642–.3738
4.3.2035–.54729.9–.1261–.3845
4.4.2566–.512010.0–.0874–.3912
4.5.3058–.472810.1–.0181–.3911
4.63510–.430010.2–.0087–.3923
4.7.3917–.381210.3+.0302

[*] Communicated by the Author.
[†] 'On the Leading Phenomena of the Wave–Making Resistance of Ships'. Trans.Inst.Naval Architects, 1881
[‡] White, p. 470
[§] See Gray and Mathews, 'Bessel Functions' p.65, and Ex.18, p.230. The formula (8) was given by Weber.
[**] ‘Messenger of Mathematics,’ 1896.

Philosophical Magazine, 1898, vol. 45, Sec. 5, pp. 106–123.

Phil. Mag, S, .5. Vol . 45 No. 272. Jan. 1898. K