Bible Encyclopedias
Hydromechanics

13. Theorem

The resultant vertical thrust on any portion of a curved surface exposed to the pressure of a fluid at rest under gravity is the weight of fluid cut out by vertical lines drawn round the boundary of the curved surface.

Theorem

The resultant horizontal thrust in any direction is obtained by drawing parallel horizontal lines round the boundary, and intersecting a plane perpendicular to their direction in a plane curve; and then investigating the thrust on this plane area, which will be the same as on the curved surface.

The proof of these theorems proceeds as before, employing the normality principle; they are required, for instance, in the determination of the liquid thrust on any portion of the bottom of a ship.

In casting a thin hollow object like a bell, it will be seen that the resultant upward thrust on the mould may be many times greater than the weight of metal; many a curious experiment has been devised to illustrate this property and classed as a hydrostatic paradox (Boyle, Hydrostatical Paradoxes, 1666).

Consider, for instance, the operation of casting a hemispherical bell, in fig. 2. As the molten metal is run in, the upward thrust on the outside mould, when the level has reached PP', is the weight of metal in the volume generated by the revolution of APQ; and this, by a theorem of Archimedes, has the same volume as the cone ORR', or rya, where y is the depth of metal, the horizontal sections being equal so long as y is less than the radius of the outside FIG. 2. hemisphere. Afterwards, when the metal has risen above B, to the level KK', the additional thrust is the weight of the cylinder of diameter KK' and height BH. The upward thrust is the same, however thin the metal may be in the interspace between the outer mould and the core inside; and this was formerly considered paradoxical.

Analytical Equations of Equilibrium of a Fluid at rest under any System of Force. 14. Referred to three fixed coordinate axes, a fluid, in which the pressure is p, the density p, and X, Y, Z the components of impressed force per unit mass, requires for the equilibrium of the part filling a fixed surface S, on resolving parallel to Ox, f flpdS = f f fpXdxdydz, (I) where 1, m, n denote the direction cosines of the normal drawn outward of the surface S.

But by Green's transformation f flpdS = f f PPdxdydz, (2) thus leading to the differential relation at every point = dy dp The three equations of equilibrium obtained by taking moments round the axes are then found to be satisfied identically. Hence the space variation of the pressure in any direction, or the pressure-gradient, is the resolved force per unit volume in that direction. The resultant force is therefore in the direction of the steepest pressure-gradient, and this is normal to the surface of equal pressure; for equilibrium to exist in a fluid the lines of force must therefore be capable of being cut orthogonally by a system of surfaces, which will be surfaces of equal pressure.

Ignoring temperature effect, and taking the density as a function of the pressure, surfaces of equal pressure are also of equal density, and the fluid is stratified by surfaces orthogonal to the lines of force; n ap, dy, P d z , or X, Y, Z (4) are the partial differential coefficients of some function P, =fdplp, of x, y, z; so that X, Y, Z must be the partial differential coefficients of a potential -V, such that the force in any direction is the downward gradient of V; and then dP dV (5) ax + Tr=0, or P+V =constant, in which P may be called the hydrostatic head and V the head of potential.

With variation of temperature, the surfaces of equal pressure and density need not coincide; but, taking the pressure, density and temperature as connected by some relation,such as the gas-equation, the surfaces of equal density and temperature must intersect in lines lying on a surface of equal pressure.

15. As an example of the general equations, take the simplest case of a uniform field of gravity, with Oz directed vertically downward; employing the gravitation unit of force, 1 dp i dp t dp dp/dz = p = pzn (4) z n+I pz 1 /n p-p n-H ?t), (5) supposing p and p to vanish together.

These equations can be made to represent the state of convective equilibrium of the atmosphere, depending on the gas-equation p = pk =RA (6) where 0 denotes the absolute temperature; and then d9 d p R dz - dz (p) n+ 1' so that the temperature-gradient deldz is constant, as in convective equilibrium in (I I).

From the gas-equation in general, in the atmosphere n d dp _ I dp 1 de _ d0 de i de (8) z p dz-edz-p-edz-k-edz' which is positive, and the density p diminishes with the ascent, provided the temperature-gradient de/dz does not exceed elk. With uniform temperature, taking h constant in the gas-equation, dp / dz= =p / k, p=poet/ k , (9) so that in ascending in the atmosphere of thermal equilibrium the pressure and density diminish at compound discount, and for pressures p 1 and 1, 2 at heights z 1 and z2 (z1-z2)11? = loge(P2891) =2.3 logio(p2/p1) (io) In the convective equilibrium of the atmosphere, the air is supposed to change in density and pressure without exchange of heat by conduction; and then PIN = (e/e0) n+1, d5 -(n-{--I) P -(n+I)R ' y - where is the ratio of the specific heat at constant pressure and constant volume.

In the more general case of the convective equilibrium of a spherical atmosphere surrounding the earth, of radius a, (1-1?-=(n+ I) Po --a 2 dr, (12) gravity varying inversely as the square of the distance r from the centre; so that, k = po/po, denoting the height of the homogeneous atmosphere at the surface, 0 is given by (n+I)k(I -9/6 0) =a(I -a/r), (13) or if c denotes the distance where 0=o, 0 _a (14) 0 r c -a' When the compressibility of water is taken into account in a deep ocean, an experimental law must be employed, such as p - po=k(P - Po), or P/po=I+(p-p0)/A, A=kpo, (15) so that A is the pressure due to a head k of the liquid at density under atmospheric pressure po; and it is the gauge pressure required on this law to double the density. Then dp/dz=kdp/dz = P, = Poe ik, p - po= kpo(ez Ik -1); (16) and if the liquid was incompressible, the depth at pressure p would be (p - po) 1po, so that the lowering of the surface due to compression is ke h I k -k -z= 1z 2 /k, when k is large. (17) For sea water, A is about 25,000 atmospheres, and k is then 25,000 times the height of the water barometer, about 250,000 metres, so that in an ocean 10 kilometres deep the level is lowered about 200 metres by the compressibility of the water; and the density at the bottom is increased 4%.

On another physical assumption of constant cubical elasticity A, dp = Ad p /P, (p - po)IA= lo g (P/Po), (18) dp _ A dp (I 1 zd p dz - P ' A Po-p -z, I - p -k, A kPo ' (19) (3) P dx Pdy Pdz -., (I) When the density p is f un dp/ iform, this becomes, as before in (2) § 9 P pp ==Pzz++paoconstant. (2) (3) Suppose the density p varies as some nth power of the depth below 0, then (7) and the lowering of the surface is 2 ° - z=klog po - z= - k log(1 - k) - zt12 k (20) Po as before in 17).

16. Centre of Pressure. - A plane area exposed to fluid pressure on one Side experiences a single resultant thrust, the integrated pressure over the area, acting through a definite point called the centre of pressure (C.P.) of the area.

Thus if the plane is normal to Or, the resultant thrust R =f fpdxdy, (r) and the co-ordinates x, y of the C.P. are given by xR = f f xpdxdy, yR = f f ypdxdy. (2) The C.P. is thus the C.G. of a plane lamina bounded by the area, in which the surface density is p. If p is uniform, the C.P. and C.G. of the area coincide.

For a homogeneous liquid at rest under gravity, p is proportional to the depth below the surface, i.e. to the perpendicular distance from the line of intersection of the plane of the area with the free surface of the liquid.

If the equation of this line, referred to new coordinate axes in the plane area, is written xcos a+y sin a - h=o, (3) R = f f p(h - x cos a - y sin a)dxdy, (4) zR= f fpx(h - xcos a - y sin a)dxdy, (5) yR = f f py(h - x cos a - y sin a)dxdy. Placing the new origin at the C.G. of the area A, f f xdxdy = o, ffydxdy = 0, (6) R = p hA, (7) xhA = - cos a f f x 2 dA - sin affxydA, fxydA, (8) yhA = - cos a ff xydA - sin ail y 2 dA. (9) Turning the axes to make them coincide with the principal axes of the area A, thus making f f xydA = o, xh = - a 2 cos a, y h = - b 2 sin a, (io) where ffx2dA=Aa2, ffy 2 dA= Ab 2 , (II) a and b denoting the semi-axes of the momental ellipse of the area. This shows that the C.P. is the antipole of the line of intersection of its plane with the free surface with respect to the momental ellipse at the C.G. of the area.

Thus the C.P. of a rectangle or parallelogram with a side in the surface is at a of the depth of the lower side; of a triangle with a vertex in the surface and base horizontal is 4 of the depth of the base; but if the base is in the surface, the C.P. is at half the depth of the vertex; as on the faces of a tetrahedron, with one edge in the surface.

The core of an area is the name given to the limited area round its C.G. within which the C.P. must lie when the area is immersed completely; the boundary of the core is therefore the locus of the antipodes with respect to the momental ellipse of water lines which touch the boundary of the area. Thus the core of a circle or an ellipse is a concentric circle or ellipse of one quarter the size.

The C.P. of water lines passing through a fixed point lies on a straight line, the antipolar of the point; and thus the core of a triangle is a similar triangle of one quarter the size, and the core of a parallelogram is another parallelogram, the diagonals of which are the middle third of the median lines.

In the design of a structure such as a tall reservoir dam it is important that the line of thrust in the material should pass inside the core of a section, so that the material should not be in a state of tension anywhere and so liable to open and admit the water.

17. Equilibrium and Stability of a Ship or Floating Body. The Metacentre. - The principle of Archimedes in § 12 leads immediately to the conditions of equilibrium of a body supported freely in fluid, like a fish in water or a balloon in the air, or like a ship (fig. 3) floating partly immersed in water and the rest in air. The body is in equilibrium under two forces: - (i.) its weight W acting vertically downward body, and (ii.) the buoyancy of the of the displaced fluid, and acting B, the C.G. of the displaced fluid; for equilibrium these two forces must be equal and opposite in the same line.

The conditions of equilibrium of a body, floating like a ship on the surface of a liquid, are therefore: (i.) the weight of the body must be less than the weight of the total volume of liquid it can displace; or else the body will sink to the bottom of the liquid; the difference of the weights is called the " reserve of buoyancy." (ii.) the weight of liquid which the body displaces in the position of equilibrium is equal to the weight W of the body; and (iii.) the C.G., B, of the liquid displaced and G of the body, must lie in the same vertical line GB.

18. In addition to satisfying these conditions of equilibrium, a ship must fulfil the further condition of stability, so as to keep upright; if displaced slightly from this position, the forces called into play must be such as to restore the ship to the upright again. The stability of a ship is investigated practically by inclining it; a weight is moved across the deck and the angle is observed of the heel produced.

Suppose P tons is moved c ft. across the deck of a ship of W tons displacement; the C.G. will move from G to G 1 the reduced distance G1G2 = c (P/W); and if B, called the centre of buoyancy, moves to B 1, along the curve of buoyancy BB 1 , the normal of this curve at B 1 will be the new vertical B1G1, meeting the old vertical in a point M, the centre of curvature of BB I, called the metacentre. If the ship heels through an angle 0 or a slope of I in m, GM =GG 1 cot 8=mc(P/W), (r) and GM is called the metacentric height; and the ship must be ballasted, so that G lies below M. If G was above M, the tangent drawn from G to the evolute of B, and normal to the curve of buoyancy, would give the vertical in a new position of equilibrium. Thus in H.M.S. " Achilles " of 9000 tons displacement it was found that moving 20 tons across the deck, a distance of 42 ft., caused the bob of a pendulum 20 ft. long to move through ro in., so that GM (2) also cot 0 = 24, 8 =2° 2 4 '. (3) In a diagram it is conducive to clearness to draw the ship in one position, and to incline the water-line; and the page can be turned if it is desired to bring the new water-line horizontal.

Suppose the ship turns about an axis through F in the water-line area, perpendicular to the plane of the paper; denoting by y the distance of an element dA if the water-line area from the axis of rotation, the change of displacement is EydA tan 8, so that there is no change of displacement if EydA = o, that is, if the axis passes through the C.G. of the water-line area, which we denote by F and call the centre of flotation.

The righting couple of the wedges of immersion and emersion will be ZwydA tan 8.y =w tan 0Zy 2 dA =w tan 0.Ak 2 ft. tons, (4) w denoting the density of water in tons/ft.', and W =wV, for a displacement of V ft.3 This couple, combined with the original buoyancy W through B, is equivalent to the new buoyancy through B, so that W.BB 1 =wAk 2 tan 8, (5) BM =BB 1 cot B=Ak e /V, (6) giving the radius of curvature BM of the curve of buoyancy B, in terms of the displacement V, and Ak e the moment of inertia of the water-line area about an axis through F, perpendicular to the plane of displacement.

An inclining couple due to moving a weight about in a ship will heel the ship about an axis perpendicular to the plane of the couple, only when this axis is a principal axis at F of the momental ellipse of the water-line area A. For if the ship turns through a small angle 0 about the line FF', then b1, b 2, the C.G. of the wedge of immersion and emersion, will be the C.P. with respect to FF' of the two parts of the water-line area, so that b 1 b 2 will be conjugate to FF' with respect to the momental ellipse at F.

The naval architect distinguishes between the stability of form, represented by the righting couple W.BM, and the stability of ballasting, represented by W.BG. Ballasted with G at B, the righting couple when the ship is heeled through 0 is given by W.BM. tan 0; but if weights inside the ship are raised to bring G above B, the righting couple is diminished by W.BG. tan 0, so that the resultant righting couple is W.GM. tan 8. Provided the ship is designed to float upright at the smallest draft with no load on board, the stability at any other draft of water can be arranged by the stowage of the weight, high or low.

19. Proceeding as in § 16 for the determination of the C.P. of an area, the same argument will show that an inclining couple due to K FIG. 3.

through G, the C.G. of the fluid, equal to the weight vertically upward through the movement of a weight P through a distance c will cause the ship to heel through an angle 0 about an axis FF' through F, which is conjugate to the direction of the movement of P with respect to an ellipse, not the momental ellipse of the water-line area A, but a confocal to it, of squared semi-axes a 2 -hV/A, b 2 - hV/A, (I) h denoting the vertical height BG between C.G. and centre of buoyancy. The varying direction of the inclining couple Pc may be realized by swinging the weight P from a crane on the ship, in a circle of radius c. But if the weight P was lowered on the ship from a crane on shore, the vessel would sink bodily a distance P/wA if P was deposited over F; but deposited anywhere else, say over Q on the water-line area, the ship would turn about a line the antipolar of Q with respect to the confocal ellipse, parallel to FF', at a distance FK from F FK= (k2-hV/A)/FQ sin QFF' (2) through an angle 0 or a slope of one in m, given by P sin B= m wA FK - W'Ak 2V hV FQ sin QFF', (3) where k denotes the radius of gyration about FF' of the water-line area. Burning the coal on a voyage has the reverse effect on a steamer.

Hydrodynamics 20. In considering the motion of a fluid we shall suppose it non-viscous, so that whatever the state of motion the stress across any section is normal, and the principle of the normality and thence of the equality of fluid pressure can be employed, as in hydrostatics. The practical problems of fluid motion, which are amenable to mathematical analysis when viscosity is taken into account, are excluded from treatment here, as constituting a separate branch called "hydraulics" (q.v.). Two methods are employed in hydrodynamics, called the Eulerian and Lagrangian, although both are due originally to Leonhard Euler. In the Eulerian method the attention is fixed on a particular point of space, and the change is observed there of pressure, density and velocity, which takes place during the motion; but in the Lagrangian method we follow up a particle of fluid and observe how it changes. The first may be called the statistical method, and the second the historical, according to J. C. Maxwell. The Lagrangian method being employed rarely, we shall confine ourselves to the Eulerian treatment.

The Eulerian Form of the Equations of Motion. 21. The first equation to be established is the equation of continuity, which expresses the fact that the increase of matter within a fixed surface is due to the flow of fluid across the surface into its interior.

I n a straight uniform current of fluid of density p, flowing with velocity q, the flow in units of mass per second across a plane area A, placed in the current with the normal of the plane making an angle 0 with the velocity, is oAq cos 0, the product of the density p, the area A, and q cos 0 the component velocity normal to the plane.

Generally if S denotes any closed surface, fixed in the fluid, M the mass of the fluid inside it at any time t, and 0 the angle which the outward-drawn normal makes with the velocity q at that point, dM/dt = rate of increase of fluid inside the surface, (I) =flux across the surface into the interior _ - f f pq cos OdS, the integral equation of continuity.

In the Eulerian notation u, v, w denote the components of the velocity q parallel to the coordinate axes at any point (x, y, z) at the time t; u, v, w are functions of x, y, z, t, the independent variables; and d is used here to denote partial differentiation with respect to any one of these four independent variables, all capable of varying one at a time.

To transfer the integral equation into the differential equation of continuity, Green's transformation is required again, namely, (++) dxdydz= ff (l +mr t } ndS, (2) or individually dxdydz = f flldS, ..., (3) where the integrations extend throughout the volume and over the surface of a closed space S; 1, m, n denoting the direction cosines of the outward-drawn normal at the surface element dS, and, 77, any continuous functions of x, y, z.

The integral equation of continuity (I) may now be written l f fdxdydz+ff (lpu+mpv+npdso, (4) which becomes by Green's transformation (dt +d dz dy dx (p u) + d (p v) + d (p w) l I dxdydz - o, dp leading to the differential equation of continuity when the integration is removed.

22. The equations of motion can be established in a similar way by considering the rate of increase of momentum in a fixed direction of the fluid inside the surface, and equating it to the momentum generated by the force acting throughout the space 5, and by the pressure acting over the surface S.

Taking the fixed direction parallel to the axis of x, the time-rate of increase of momentum, due to the fluid which crosses the surface, is - f'fpuq cos OdS = - f f (lpu 2 -+mpuv+npuw)dS, (1) which by Green's transformation is (d(uiu 2) dy dz dxdydz. (2) y The rate of generation of momentum in the interior of S by the component of force, X per unit mass, is fffpXdxdydz, f pXdxdydz, (3) and by the pressure at the surface S is -f. flpdS= _ dx dxdydz, (4) by Green's transformation.

The time rate of increase of momentum of the fluid inside S is )dxdydz; (5) and (5) is the sum of (I), (2), (3), (4), so that /if (dpu+dpu2+dpuv +dpuw_ +d p j d xdyd z = o, (b)` dt dx dy dz dx / leading to the differential equation of motion dpu dpu 2 dpuv dpuv _ X_ (7) dt + dx + dy + dz with two similar equations.

The absolute unit of force is employed here, and not the gravitation unit of hydrostatics; in a numerical application it is assumed that C.G.S. units are intended.

These equations may be simplified slightly, using the equation of continuity (5) § for dpu dpu 2 dpuv dpuw dt dx + dy + dz =p Cat +uax+vay+waz? dp dpu dpv dpw -z)' reducing to the first line, the second line vanishing in consequence of the equation of continuity; and so the equation of motion may be written in the more usual form du du du du d dt +udx+vdy +wdz =X -n dx' with the two others dv dv dv dv i dp dt +u dx +v dy +w dz - Y -P d y' dw dw dw Z w dw i d p dt +u dx +v dy +wd - -P dz.

23. As a rule these equations are established immediately by determining the component acceleration of the fluid particle which is passing through (x, y, z) at the instant t of time considered, and saying that the reversed acceleration or kinetic reaction, combined with the impressed force per unit of mass and pressure-gradient, will according to d'Alembert's principle form a system in equilibrium.

To determine the component acceleration of a particle, suppose F to denote any function of x, y, z, t, and investigate the time rate of F for a moving particle; denoting the change by DF/dt, DF = 1t F(x+uSt, y+vIt, z+wSt, t+St) - F(x, y, z, t) dt at = d + u dx +v dy+ w dz and D/dt is called particle differentiation, because it follows the rate of change of a particle as it leaves the point x, y, z; but dF/dt, dF/dx, dF/dy, dF/dz (2) represent the rate of change of F at the time t, at the point, x, y, z, fixed in space.

(5) (8) (I) The components of acceleration of a particle of fluid are consequently Du dudu du du dt = dt +u dx +v dy + wdz' Dr dv dv dv dv dt -dt+udx+vdy+wdz' dt v = dtJ+udx+vdy +w dx' leading to the equations of motion above.

If F (x, y, z, t) =o represents the equation of a always the same particles of fluid, DF =o, or dF {-u dx-rzd { w d d = o, Trt y _ which is called the differential equation of the bounding surface. A bounding surface is such that there is no flow of fluid across it, as expressed by equation (6). The surface always contains the same fluid inside it, and condition (6) is satisfied over the complete surface, as well as any part of it.

But turbulence in the motion will vitiate the principle that a bounding surface will always consist of the same fluid particles, as we see on the surface of turbulent water.

24. To integrate the equations of motion, suppose the impressed force is due to a potential V, such that the force in any direction is the rate of diminution of V, or its downward gradient; and then X= -dV/dx, Y= -dV/dy, Z= -dV/dz; (I) and putting dw dv du dw dv du Ty - dz -2 ' dz - dx -2n ' dx - dy2?, d -{- d ' v ? d? = dx dy dz the equations of motion may be Written du - 2v? 2wr { a 0, dt2WE+2UC+ dz = o, dw dt - 2un+2v+ dH = 0, where H = fdp/p +V +1q 2, (7) 2 2 +v 2 2 (8) and the three terms in H may be called the pressure head, potential head, and head of velocity, when the gravitation unit is employed and Zq 2 is replaced by 1q 2 1 g. Eliminating H between (5) and (6) DS du dv dw (du dv d1zv dt u dx n dx udx' 5 -, dzi =°' and combining this with the equation of continuity Dp du dv dw p iit dx+dy+ dz = °' (10) D i du n dv dw_ dt (p p dx p dx p dx - o, with two similar equations. Putting (12) a vortex line is defined to be such that the tangent is in the direction of w, the resultant of, n, called the components of molecular rotation. A small sphere of the fluid, if frozen suddenly, would retain this angular velocity.

If w vanishes throughout the fluid at any instant, equation (I I) shows that it will always be zero, and the fluid motion is then called irrotational; and a function 4) exists, called the velocity function, such that udx+vdy-{-wdz =-d, (13) and then the velocity in any direction is the space-decrease or downward gradient of cp.

25. But in the most general case it is possible to have three functions 0, ?, m of x, y, z, such that udx+vdy +wdz = -dcp-mdl,G, (I) as A. Clebsch has shown, from purely analytical considerations (Crelle, lvi.); and then = Z d(?G, m), ? = 1d(1G, m), ?= 1 d(?', m) d(y, z) d(z, x)' y), and , a ?+n d +s m - O,, d +n dinn +? d o, dx dy dz dx dy dz so that, at any instant, the surfaces over which tk and m are constant intersect in the vortex lines.

Putting do dipdp _ Hdt - dt -K, the equations of motion (4), (5), (6) § 24 can be written -2u +2wn - d(x',t))o,...,. ... and therefore dK dK dK cWcmn d y - -? dz =o. Equation (5) becomes, by a rearrangement, dK dmdm dm din dx dt +u dx + dy +Zee dz + dx (dt +u dx +v dy +w d) = o,. .

dK d1. Dm dm Dl (8) ,dx - dx dt + dx dt = °' ...' ...' and as we prove subsequently (§ 37) that the vortex lines are composed of the same fluid particles throughout the motion, the surface m and satisfies the condition of (6) § 23; so that K is uniform throughout the fluid at any instant, and changes with the time only, and so may be replaced by F(t).

26. When the motion is steady, that is, when the velocity at any point of space does not change with the time, dK dx-2v{ +2wn = o ,.. .,. .. d - K dK dK _ dK dK dK ?dx n dyd °, udx dz - ° and K=fdp/o+V+2q 2 =H (3) is constant along a vortex line, and a stream line, the path of a fluid particle, so that the fluid is traversed by a series of H surfaces, each covered by a network of stream lines and vortex lines; and if the motion is irrotational H is a constant throughout the fluid.

Taking the axis of x for an instant in the normal through a point on the surface H = constant, this makes u = o, = o; and in steady motion the equations reduce to dH/dv=2q-2wn = 2gco sin e, (4) where B is the angle between the stream line and vortex line; and this holds for their projection on any plane to which dv is drawn perpendicular.

In plane motion (4) reduces to dH = 2q"= q /av q? if r denotes the radius of curvature of the stream line, so that I dp + dV - dH _ dq 2 q2 (6) p dv dv dv dv - r ' the normal acceleration.

The osculating plane of a stream line in steady motion contains the resultant acceleration, the direction ratios of which are du du, du d i g d g 2 _ dH dx +v dy + dz - 2v? }-2w ? = dx dx' ... , (7) and when q is stationary, the acceleration is normal to the surface H = constant, and the stream line is a geodesic.

Calling the sum of the pressure and potential head the statical head, surfaces of constant statical and dynamical head intersect in lines on H, and the three surfaces touch where the velocity is stationary.

Equation (3) is called Bernoulli's equation, and may be interpreted as the balance-sheet of the energy which enters and leaves a given tube of flow.

If homogeneous liquid is drawn off from a vessel so large that the motion at the free surface at a distance may be neglected, then Bernoulli's equation may be written H = PIP--z - F4 2 / 2g = P/ p +h, (8) where P denotes the atmospheric pressure and h the height of the free surface, a fundamental equation in hydraulics; a return has been made here to the gravitation unit of hydrostatics, and Oz is taken vertically upward.

In particular, for a jet issuing into the atmosphere, where p=P, q 2 /2g = h - z, (9) or the velocity of the jet is due to the head k-z of the still free surface above the orifice; this is Torricelli's theorem (1643), the foundation of the science of hydrodynamics.

27. Uniplanar Motion. - In the uniplanar motion of a homogeneous liquid the equation of continuity reduces to du dv dx' dy-O' u= -d,y/dy, v = d i t/dx, (2) surface containing so that we can put _ (6) (9) we have (I) (2) (5) (I) where 4 is a function of x, y, called the streamor current-function; interpreted physically, 4-4c, the difference of the value of 4, at a fixed point A and a variable point P is the flow, in ft. 3 / second, across any curved line AP from A to P, this being the same for all lines in accordance with the continuity.

Thus if d,/ is the increase of 4, due to a displacement from P to P', and k is the component of velocity normal to PP', the flow across PP' is d4 = k.PP'; and taking PP' parallel to Ox, d,, = vdx; and similarly d/ ' = -udy with PP' parallel to Oy; and generally d4,/ds is the velocity across ds, in a direction turned through a right angle forward, against the clock.

In the equations of uniplanar motion = dx - du = dx + dy = -v 2 ?, suppose, so that in steady motion dx I +v24 ' x = ?' dy I +v2" dy = 0' d4' Y' =o, and 2 must be a function of 4'. Y If the motion is irrotational, u=-x-- dy' 2' d y = dx' y y so that :(, and 4' are conjugate functions of x and y, 0+4,i = f(x + y i), v 2 4 =o, v 2 0 =o; or putting 0+0=w, +yi=z, w=f(z).

The curves 0 = constant and 4, = constant form an orthogonal system; and the interchange of 0 and 4, will give a new state of uniplanar motion, in which the velocity at every point is turned through a right angle without alteration of magnitude.

For instance, in a uniplanar flow, radially inward towards 0, the flow across any circle of radius r being the same and denoted by 27rm, the velocity must be mfr, and 0=m log r, ,y=m0, +4,i =m log re ie , w=m log z. (7) Interchanging these values =m log r, 4, = mO, 4,+4,i =m log rei e (8) gives a state of vortex motion, circulating round Oz, called a straight or columnar vortex.

A single vortex will remain at rest, and cause a velocity at any point inversely as the distance from the axis and perpendicular to its direction; analogous to the magnetic field of a straight electric current.

If other vortices are present, any one may be supposed to move with the velocity due to the others, the resultant stream function being = gy m log r =log IIrm; (9) the path of a vortex is obtained by equating the value of 1P at the vortex to a constant, omitting the rm of the vortex itself.

When the liquid is bounded by a cylindrical surface, the motion of a vortex inside may be determined as due to a series of vorteximages, so arranged as to make the flow zero across the boundary.

For a plane boundary the image is the optical reflection of the vortex. For example, a pair of equal opposite vortices, moving on a line parallel to a plane boundary, will have a corresponding pair of images, forming a rectangle of vortices, and the path of a vortex will be the Cotes' spiral r sin 20 = 2a, or x-2+y-2=a-2; (io) this is therefore the path of a single vortex in a right-angled corner; and generally. if the angle of the corner is jr/n, the path is the Cotes' spiral r sin n0=na. (II) A single vortex in a circular cylinder of radius a at a distance c from the centre will move with the velocity due to an equal opposite image at a distance a 2 /c, and so describe a circle with velocity mc/(a 2 -c 2)in the periodic time 21r(a 2 -c 2)/m. (22) Conjugate functions can be employed also for the motion of liquid in a thin sheet between two concentric spherical surfaces; the components of velocity along the meridian and parallel in colatitude 0 and longitude A can be written d¢_ i _ d4, I dip _ dy (13) d8 sin - 0 dX' sin 0 dX de' and then = F (tan O. eAi). (14) 28. Uniplanar Motion of a Liquid due to the Passage of a Cylinder through it.-A stream-function 4, must be determined to satisfy the conditions v24 =o, throughout the liquid; (I) I =constant, over any fixed boundary; (2) d,t/ds = normal velocity reversed over a solid boundary, (3) so that, if the solid is moving with velocity U in the direction Ox, d4y1ds=-Udy/ds, or 0 +Uy =constant over the moving cylinder; and 4,+Uy=41' is the stream function of the relative motion of the liquid past the cylinder, and similarly 4,-Vx for the component velocity V along Oy; and generally 1,1'= +Uy -Vx (4) is the relative stream-function, constant over a solid boundary moving with components U and V of velocity.

If the liquid is stirred up by the rotation R of a cylindrical body, d4lds = normal velocity reversed dy = - Rx- Ry ds (5) ds 4' + 2 R (x2 + y2) = Y, (6) a constant over the boundary; and 4,' is the current-function of the relative motion past the cylinder, but now V 2 4,'+2R =o, (7) throughout the liquid.

Inside an equilateral triangle, for instance, of height h, - 2Ra/3y/h, (8) where a, 13, y are the perpendiculars on the sides of the triangle.

In the general case 4'=11.+Uy-Vx+2R(x 2 +y 2) is the relative stream function for velocity components, U, V, R.

29. Example z.-Liquid motion past a circular cylinder. Consider the motion given by w=U(z+a2/z), (I) 4,=U(r+- r ) cos 0= U + a1x, so that (2) = U (r-)sin 0= U(i -¢) y. Then 4, =o over the cylinder r = a, which may be considered a fixed post; and a stream line past it along which 4, = Uc, a constant, is the curve (r - ¢2) sin 0=c, (x2 + y2) (y - c) - a 2 y = o, (3) a cubic curve (C3).

Over a concentric cylinder, external or internal, of radius r=b, 4,'=4,+ Uly =[U I - + Ui]y, (4) and 4" is zero if U 1 /U = (a 2 - b2)/b 2; (5) so that the cylinder may swim for an instant in the liquid without distortion, with this velocity Ui; and w in (I) will give the liquid motion in the interspace between the fixed cylinder r =a and the concentric cylinder r=b, moving with velocity U1.

When b = o, U 1 =00; and when b = oo, U 1 = -U, so that at infinity the liquid is streaming in the direction xO with velocity U. If the liquid is reduced to rest at infinity by the superposition of an opposite stream given by w = - Uz, we are left with w = Ua2/z, (6) =U(a 2 /r) cos 0= Ua2x/(x2+y2), (7) 4, = -U(a 2 /r) sin 0= -Ua2y/( x2+y2), (8) giving the motion due to the passage of the cylinder r=a with velocity U through the origin 0 in the direction Ox.

If the direction of motion makes an angle 0' with Ox, tan B' = d0 !dam _ ?xy 2 = tan 20, 0 =-10', (9) dy/ y and the velocity is Ua2/r2.

Along the path of a particle, defined by the of (3), _ c) sine 2e, - x 2 + y2 = y a 2 ' (Io) sin B' de' _ 2y-c dy 2 ds ds' on the radius of curvature is 4a 2 /(ylc), which shows that the curve is an Elastica or Lintearia. (J. C. Maxwell, Collected Works, ii. 208.) If 01 denotes the velocity function of the liquid filling the cylinder r = b, and moving bodily with it with velocity Ul, 41 = -U1x, (12) and over the separating surface r =b 4, = I U (+- a2) a2 +b2 (13) 1 Ul b2 - a 2 ... b2' and this, by § 36, is also the ratio of the kinetic energy in the annular 4,1 interspace between the two cylinders to the kinetic energy of the liquid moving bodily inside r = b. Consequently the inertia to overcome in moving the cylinder r=b, solid or liquid, is its own inertia, increased by the inertia of liquid (a2+b2)/(a2,..b2) times the volume of the cylinder r=b; this total inertia is called the effective inertia of the cylinder r =b, at the instant the two cylinders are concentric.

With liquid of density p, this gives rise to a kinetic reaction to acceleration dU/dt, given by 7rp b 2 a 2 b l b d J = a 2 +b2 M' dU, if M' denotes the mass of liquid displaced by unit length of the cylinder r =b. In particular, when a = oo, the extra inertia is M'. When the cylinder r =a is moved with velocity U and r =b with velocity U 1 along Ox, = U b e - a,1 r +0 cos 0 - U ib2 - 2 a, (r +Q 2 ') cos 0, = - U be a2 a2 (b 2 - r) sin 0 - Uib2 b1)a, (r - ¢2 sin 0; b and similarly, with velocity components V and V 1 along Oy a 2 b2 ?= Vb,_a,(r+r) sin g -Vi b , b2 a ,(r+ 2 ) sin 0, (17) = V b, a2 a, (b2 r ) cos 0+Vi b, b, a, (r- ¢ 2) cos h; (18) and then for the resultant motion z 2zz w= (U 2 + V2)b2a a2U+Vi +b a b a2 U z Vi -(U12+V12) b2 z a2b2 Ui +VIi b 2 - a 2 U1 +Vii b 2 - a 2 z The resultant impulse of the liquid on the cylinder is given by the component, over r=a (§ 36), X =f p4 cos 0.ad0 =7rpa 2 (U b z 2 + a 2 Uib.2bz a2); (20) and over r =b Xi= fp? cos 0. bdo 7rpb 2 (u, b 2 a2 Uibb +¢z), and the difference X-X 1 is the component momentum of the liquid in the interspace; with similar expressions for Y and Y1.

Then, if the outside cylinder is free to move - 2a2 2 X 1 = 0, T T = b2 a2, 7rpa 2 Ub 2+ a 2. (22) But if the outside cylinder is moved with velocity U1, and the inside cylinder is solid or filled with liquid of density v, 2 U i 2pb2 and the inside cylinder starts forward or backward with respect to the outside cylinder, according as p> or < v.

30. The expression for w in (i) § 29 may be increased by the addition of the term im log z =-m0 + im log r, (1) representing vortex motion circulating round the annulus of liquid.

Considered by itself, with the cylinders held fixed, the vortex sets up a circumferential velocity m/r on a radius r, so that the angular momentum of a circular filament of annular cross section dA is pmdA, and of the whole vortex is pm7r(b2-a2).

Any circular filament can be started from rest by the application of a circumferential impulse 7rpmdr at each end of a diameter; so that a mechanism attached to the cylinders, which can set up a uniform distributed impulse rpm across the two parts of a diameter in the liquid, will generate the vortex motion, and react on the cylinder with an impulse couple-pmira 2 and pm7rb 2 , having resultant pm7r(b 2 -a 2), and this couple is infinite when b = oo, as the angular momentum of the vortex is infinite. Round the cylinder r=a held fixed in the U current the liquid streams past with velocity q' =2U sin 0+m/a; (2) and the loss of head due to this increase of velocity from U to q' is q' 2 -U 2 - (2U sin e to space filled with liquid, and at rest at infinity, the cylinder will experience components of force per unit length (i.) -27rpmV, 27rpmU, due to the vortex motion; 2 dU 2dV (ii.) -71-pa 2 w,, -7rpa dt , due to the kinetic reaction of the liquid; (iii.) o, -7r(a-p)a 2 g, due to gravity, taking Oy vertically upward, and denoting the density of the cylinder by a; so that the equations of motion are 71-0-a 2 - di r = - 7pa2- -- 22rpmV, (4) aa 2 - = -7rpa 2 dV +27rpmV - 7r(cr - p) a2g , (5) 7r or, putting m = a 2 w, so that the vortex velocity is due to an angular velocity w at a radius a, (o+p)dU/dt+2pwV =o, (6) (a+ p) dV /dt - 2 pwU + (v - p)g = o. (7) Thus with g=o, the cylinder will describe a circle with angular velocity 2pw/(a+p), so that the radius is (a+p)v/2pw, if the velocity is v. With v=o, the angular velocity of the cylinder is 2w; in this way the velocity may be calculated of the propagation of ripples and waves on the surface of a vertical whirlpool in a sink.

Restoring v will make the path of the cylinder a trochoid; and so the swerve can be explained of the ball in tennis, cricket, baseball, or golf.

Another explanation may be given of the sidelong force, arising from the velocity of liquid past a cylinder, which is encircled by a vortex. Taking two planes x = =b, and considering the increase of momentum in the liquid between them, due to the entry and exit of liquid momentum, the increase across dy in the direction Oy, due to elements at P and P' at opposite ends of the diameter PP', is pdy (U - Ua 2 r2 cos 20 +mr i sin 0) (Ua 2 r 2 sin 2 0+mr 1 cos 0) + pdy (- U+Ua 2 r 2 cos 2 0 +mr1 sin 0) (Ua 2 r 2 sin 2 0 -mr 1 cos 0) =2pdymUr '(cos 0 -a 2 r 2 cos 30), (8) and with b tan r =b sec this is 2pmUdo(i -a 2 b2 cos 30 cos 0), (9) and integrating between the limits 0 = 27r, the resultant, as before, is 27rpmU. 31. Example 2. - Confocal Elliptic Cylinders. - Employ the elliptic coordinates n,, and -=n+Vi, such that z=cch?, cchncos,y=cshnsin-; (1) then the curves for which n and are constant are confocal ellipses and hyperbolas, and -d(n,) =c 2 (ch 2 n - cost) = 2c 2 (ch2n-cos2) = r i r 2 = OD 2, (2) if OD is the semi-diameter conjugate to OP, and ri, r 2 the focal distances, rl,r2 = c (ch n cos 0; r 2 = x2 +y2 = c 2 (ch 2 n - sin20 = 1c 2 (ch 2 7 7 +cos 2?).

Consider the streaming motion given by w =m =a+si, (5) 4=m ch (n -a)cos(-0), p=m sh(n-a)sin(-13). (6) Then =o over the ellipse n = a, and the hyperbola t = (, so that these may be taken as fixed boundaries; and %,1. is a constant on a C4. Over any ellipse n, moving with components U and V of velocity, =i+Uy-Vx=[msh(n-a) cos (3+Ucshn] sin k -[msh(n-a) sin (3+Vcchn] cos h; (7) so that ' =o, if U c sh n cos R, V = c ch n sin a, (8) m sh(n - a) m sh(n - a). having a resultant in the direction PO, where P is the intersection of an ellipse n with the hyperbola 13; and with this velocity the ellipse n can be swimming in the liquid, without distortion for an instant.

At infinity U = -me a cos (i = a m b oos (3, V= -me a sin 1 3 - C7,1 sin 0, (9) a and b denoting the semi-axes of the ellipse a; so that the liquid is streaming at infinity with velocity Q = m/(a+b) in the direction of the asymptote of the hyperbola (3. An ellipse interior to n = a will move in a direction opposite to the exterior current; and when n = o, U = oo, but V = (m/c) sh a sin 13. Negative values of n must be interpreted by a streaming motion on a parallel plane at a level slightly different, as on a double Riemann sheet, the stream passing from one sheet to the other across a cut SS' joining the foci S, S'. A diagram has been drawn by Col. R. L. Hippisley.

U -U? (p -a)(b2 -a2) U i p(b2 +a2) +0-(b2-a2)' (23) (19) (21) (14) (15) (16) X =-7rUa U, U -p(bz +a2) +0_(b2-a2), 2g 2g , (3) so that cavitation will take place, unless the head at a great distance exceeds this loss.

The resultant hydrostatic thrust across any diametral plane of the cylinder will be modified, but the only term in the loss of head which exerts a resultant thrust on the whole cylinder is 2mU sin Olga, and its thrust is 27rpmU absolute units in the direction Cy, to be counteracted by a support at the centre C; the liquid is streaming past r=a with velocity U reversed, and the cylinder is surrounded by a vortex. Similarly, the streaming velocity V reversed will give rise to a thrust 27rpmV in the direction xC. Now if the cylinder is released, and the components U and V are reversed so as to become the velocity of the cylinder with respect +m /a) 2 - U2 The components of the liquid velocity q, in the direction of the normal of the ellipse n and hyperbola t, are -mJi sh(n--a)cos(r-a),mJ2 ch(n-a) sin (E-a). (io) The velocity q is zero in a corner where the hyperbola a cuts the ellipse a; and round the ellipse a the velocity q reaches a maximum when the tangent has turned through a right angle, and then q _ (Ch 2a-C0s 2(3). (II) a e sh 2a and the condition can be inferred when cavitation begins. With #=o, the stream is parallel to xo, and 4)=m ch (n-a)cos = - Uc ch (n-a) sh n cos /sh (n-a) (22) over the cylinder n, and as in (12) § 29, =-Ux =-Uc ch n cos t, (23) for liquid filling the cylinder; and _ th n (14) 01 th (7 7 - a) ' over the surface of n; so that parallel to Ox, the effective inertia of the cylinder n, displacing M' liquid, is increased by M'thn/th(n-a), reducing when a= oo to /If' th n = M' (b/a). Similarly, parallel to Oy, the increase of effective inertia is NT'/th n th(n-a), reducing to M'/th n=M' (a/b), when a= oo, and the liquid extends to infinity.

32. Next consider the motion given by = m ch 2(77a)sin 2E, tii= -m sh 2(na)cos 2E; (I) in which > ' =o over the ellipse a, and =1'+IR(x2+y2) =[ -m sh 2(7 7 -a)+4Rc 2 ]cos 4Rc2 ch 2n, (2) which is constant over the ellipse n if 4Rc 2. msh2(n-a); (3) so that this ellipse can be rotating with this angular velocity R for an instant without distortion, the ellipse a being fixed.

For the liquid filling the interior of a rotating elliptic cylinder of cross section x2/a2+y2/b2 = 1, /(4) = m l (x 2 / a2 - - y2 /b 2) (5) with V21G1' =-2R =-2 m 1 (I / a2 + 21b2), 214 = m l (x2 / a2 + y2 / b2) - IR(x2+y2) = I R (x2 - y2) (a 2 - b2)/(a2+b2), cp 1 = Rxy (a 2 - b2)/(a2 +b2), w1 = cb1 + % Pli = - IiR(x +yi)2(a2b2)/(a2+b2).

The velocity of a liquid particle is thus (a 2 - b 2)/(a 2 +b 2) of what it would be if the liquid was frozen and rotating bodily with the ellipse; and so the effective angular inertia of the liquid is (a 2 -b 2) 2 /(a 2 +b 2) 2 of the solid; and the effective radius of gyration, solid and liquid, is given by k 2 = 4 (a 2 2) , and 4 (a 2 For the liquid in the interspace between a and n, m ch 2(0-a) sin 2E 4) 1 4Rc 2 sh 2n sin 2E (a2_ b2)I(a2+ b2) = I/th 2 (na)th 2n; (8) and the effective k 2 of the liquid is reduced to 4c 2 /th 2 (n-a)sh 2n, (9) which becomes 4c 2 /sh 2n = s (a 2 - b 2)/ab, when a =00, and the liquid surrounds the ellipse n to infinity.

An angular velocity R, which gives components - Ry, Ix of velocity to a body, can be resolved into two shearing velocities, -R parallel to Ox, and R parallel to Oy; and then ik is resolved into 4'1+1'2, such that 4/ 1 -R-Rx 2 and 1//2+IRy2 is constant over the boundary.

Inside a cylinder ?1 = - I iR (' x + yi) 2a2 /(a2 +b 2 '), (Io) 92+1 i 2 i = I ZR ( x + yi) 2b2 I (a2 +b 2), (21) and for the interspace, the ellipse being fixed, and a l revolving with angular velocity R (1 +11/11= - giRc 2 sh 2 (n-a+ i)(ch 2a+I)/sh 2 (a i -a), (12) 42+1,'21 = *iRc 2 sh 2(na+Ei)(ch 2a-1)/sh 2(a i - a), (13) satisfying the condition that 4/ 1 and ' // 2 are zero over n = a, and over n =a 1 constant values.

In a similar way the more general state of motion may be analysed, given by w =r ch2('-y), y =a+, i, (26) as giving a homogeneous strain velocity to the confocal system; to which may be added a circulation, represented by an additional term in w. x+yi =c1,1 [sin(+ 7 ni)] (17) i ' =Qc sh((n-a)sin((E-,6) (18) 'will give motion streaming past the fixed cylinder n = a, and dividing along t =43; and then x 2 -3/ 2 = c 2 sin ch n, 2xy = c 2 cos sh n.

In particular, with sh a =1, the cross-section of n = a is x 4 +6x 2 y 2 4= 2c 4, or x 4 -{-y =c 4 (20) when the axes are turned through 45°.

33. Example 3.-Analysing in this way the rotation of a rectangle filled with liquid into the two components of shear, the stream function 1//1 is to be made to satisfy the conditions v 2 /1 =0, 111+IRx 2 = IRa 2, or /11 =o when x= = a, +b1+IRx 2 = I Ra2, y ' 1 = IR(a 2 -x 2), when y = b Expanded in a Fourier series, 2 232 2 cos(2n+ I)Z?rx/a a -x 7r3 a Lim (2n+I) 3 ' (1) so that '?" 16 ? cos (2n + I) Iirx /a. ch (2n +1)I 7 ry /a yl-R3ct (2n +I)3.ch(2n+I)17b /a ' 16 cos(2n+I) 2 7 z /a w1=4,i+ 4, ii = iR ?3a2 an elliptic-function Fourier series; with a similar expression for 1,'2 with x and y, a and b interchanged; and thence 4, = '1 +h.

Example 4.- Parabolic cylinder, axial advance, and liquid streaming past.

The polar equation of the cross-section being rI cos 19 =al, or r + x = 2a, (3) the conditions are satisfied by = Ur sin g -2Uairi sin IB = 2Uri sin 10(14 cos 18a'), (4) 1J/ =2Uairi sin IO = -U1/ [2a(r-x)], (5) w =-2Uaiz1, (6) and the resistance of the liquid is 2lrpaV2/2g.

A relative stream line, along which 1/,' = Uc, is the quartic curve y-c=?![2a(r-x)], x = 4a2y2-(y g)4, r- 4a2y2 +(y c) 4 , 7) 4 a (y-c) 4a(y and in the absolute space curve given by 1', dy= (y- c)2, x= 2ac_ 2a log ( y -c) (8) 2ay y - c 34. Motion symmetrical about an Axis.-When the motion of a liquid is the same for any plane passing through Ox, and lies in the plane, a function ' can be found analogous to that employed in plane motion, such that the flux across the surface generated by the revolution of any curve AP from A to P is the same, and represented by 2s-4 -11'o); and, as before, if d is the increase in due to a displacement of P to P', then k the component of velocity normal to the surface swept out by PP' is such that 274=2.7ryk.PP'; and taking PP' parallel to Oy and Ox, u= -d/ydy, v=dl,t'/ydx, (I) and 1P is called after the inventor, " Stokes's stream or current function," as it is constant along a stream line (Trans. Camb. Phil. Soc., 1842; " Stokes's Current Function," R. A. Sampson, Phil. Trans., 1892); and dly/yds is the component velocity across ds in a direction turned through a right angle forward.

In this symmetrical motion =o, n=O' 2s,=27c dx ) d (?' z) = I (ax - I d ? 1 = yv24, (2) y 2 y y y suppose; and in steady motion, + y 2 dx v-t ' = o, dH +y 2dy0 2P = o, so that 2 "/ y = - y2, 7 2 1,G = dH/d is a function of 1,G, say f'(> '), and constant along a stream line; dH/dv = 2qi', H -f (1/.) = constant, throughout the liquid.

When the motion is irrotational, dq_ _I d deId> G =o, a=-dxy dy, v dy ydx' v 21 ,' = o, or dx + dy -y chi, '1/4724, 4 1 1+1 Rx2 = $Rc 2 (ch 2 a1 +I), +h+I Ry2 = 8Rc 2 (ch 2a 1 - I), (6) (7) b2)2/(a2 + b2).

Similarly, with the function (19) (2n+ I) 3 ch (2n+ I) ITrb/a' (2) Changing to polar coordinates, x =r cos 0, y = r sin 0, the equation (2) becomes, with cos 0 =µ, r'd + (I -µ 2)-d µ = 2 ?-r3 sin 0, (8) of which a solution, when = o, is = (Ar'+) _(Ari_1+) y2, ,?

0= {(n +I)Ar" - where P. denotes the zonal harmonic of the nth order; also, in the exceptional case of =Ao cos 0, 4) = Ao/r; 4'= Bor, 49 = - Bo log tan 2B sh - lx/y. (II) Thus cos 0 is the Stokes' function of a point source at 0, and Papb of a line source AB.

The stream function ,y of the liquid motion set up by the passage of a solid of revolution, moving with axial velocity U, is such that y Glib = - 15 42, iI ' + Uy 2 =cons t ant, (12) over the surface of the solid; and 4, must be replaced by41' =1l.-1-1-Uy2 in the general equations of steady motion above to obtain the steady relative motion of the liquid past the solid.

For instance, with n = I in equation (9), the relative stream function is obtained for a sphere of radius a, by making it ,y' =1y+2Uy 2 = 2U(r 2 -a 3 /r) sin? 0, 1, = -ZUa 3 sin e B/r; (13) and then =Ux(I+1a3/r2), 4=ZUa 3 cos 0/r2, -d r = Ua3 cos B, -d9=ZUa3 sin 0, so that, if the direction of motion makes an angle >G with Ox, tan (4y-0) =Z tan 0, tan =3 tan 0/(2-tan 2 e). (16) Along the path of a liquid particle 4)' is constant, and putting it equal to 2Uc2, (r 2 - a 3 /r) sin 2 0 = c 2 , sin 2 0 = c2r/(r3 - a3), (17) the polar equation; or y 2 = c2r3/(r3 - a 3), r3 = a3y2 /(y2._ c2), (18) a curve of the 10th degree (C10). In the absolute path in space cos Ili = (2 - 3 sin 2 6)/1/ (4-sin 2 6), and sin 3 B = (y 3 -c 2 y)/a 3 , (19) which leads to no simple relation.

The velocity past the surface of the sphere is dC r sin 0 dy 2U (2r+ a 2) r sin g z U sin e, when r =a; (20) so that the loss of head is (! sin e 0 - i) U 2 /2g, having a maximum a U 2 /2g, (21) which must be less than the head at infinite distance to avoid cavitation at the surface of the sphere.

With n =2, a state of motion is given by ty = - ZUy 2 a 4 µ/r 4, ' _ 0' = Ux+.4), 4'--3U(a 4 /r 3) P 2, P2 = Aµ 2 2, (23) representing a stream past the surface r4=a4µ.

35A circular vortex, such as a smoke ring, will set up motion symmetrical about an axis, and provide an illustration; a half vortex ring can be generated in water by drawing a semicircular blade a short distance forward, the tip of a spoon for instance. The vortex advances with a certain velocity; and if an equal circular vortex is generated coaxially with the first, the mutual influence can be observed. The first vortex dilates and moves slower, while the second contracts and shoots through the first; after which the motion is reversed periodically, as if in a game of leap-frog. Projected perpendicularly against a plane boundary, the motion is determined by an equal opposite vortex ring, the optical image; the vortex ring spreads out and moves more slowly as it approaches the wall; at the same time the molecular rotation, inversely as the cross-section of the vortex, is seen to increase. The analytical treatment of such vortex rings is the same as for the electro-magnetic effect of a current circulating in each ring.

36. Irrotational Motion in General.-Liquid originally at rest in a singly-connected space cannot be set in motion by a field of force due to a single-valued potential function; any motion set up in the liquid must be due to a movement of the boundary, and the motion will be irrotational; for any small spherical element of the liquid may be considered a smooth solid sphere for a moment, and the normal pressure of the surrounding liquid cannot impart to it any rotation.

The kinetic energy of the liquid inside a surface S due to the velocity function 4' f i (s given by T=2p + (d) 2+ (t) dxdydz, pff f 75 4 dS (I) by Green's transformation, dv denoting an elementary step along the normal to the exterior of the surface; so that d4ldv = o over the surface makes T = o, and then ( d4 2 d4) 2 'x) + (dy) + (= O, dd? = o, dyd? = 0, dz = O. (2) If the actual motion at any instant is supposed to be generated instantaneously from rest by the application of pressure impulse over the surface, or suddenly reduced to rest again, then, since no natural forces can act impulsively throughout the liquid, the pressure impulse W satisfies the equations I do = I d i dos - ax -u, - - y = -v, Pdz = -t, a =p4)-}-a constant, (4) and the constant may be ignored; and Green's transformation of the energy T amounts to the theorem that the work done by an impulse is the product of the impulse and average velocity, or half the velocity from rest.

In a multiply connected space, like a ring, with a multiply valued velocity function ¢, the liquid can circulate in the circuits independently of any motion of the surface; thus, for example, 4) =mB=m tan - l y/x (5) will give motion to the liquid, circulating in any ring-shaped figure of revolution round Oz.

To find the kinetic energy of such motion in a multiply connected space, the channels must be supposed barred, and the space made acyclic by a membrane, moving with the velocity of the liquid; and then if k denotes the cyclic constant of 0 in any circuit, or the value by which 4) has increased in completing the circuit, the values of 0 on the two sides of the membrane are taken as differing by k, so that the integral over the membrane dS= k f ? d ddS, and this term is to be added to the terms in (I) to obtain the additional part in the kinetic energy; the continuity shows that the integral is independent of the shape of the barrier membrane, and its position. Thus, in (5), the cyclic constant k = 27rm.

In plane motion the kinetic energy per unit length parallel to Oz T 2p J J [ (d4)) 2+ ( d dy (P)1dxdy=lpfl[ a) 2+ (=zp 4d ds=zp f ,ydvds. (7) For example, in the equilateral triangle of (8) § 28, referred to coordinate axes made by the base and height, = -2 Ra,(3y/h= - 1. Ry[(h - y) 2 ]/h (8) 4' =4G'- zR[(ahy)2+x2] - ZRRh 3 +3h 2 y+h)x 2 -y 2) -3 x2y + y3 ]/ h (9) and over the base y =o, dx/dv= -dx/dy = +Z R(111 2 - 3x 2)/h4 = - 1 R (s h2 -1-x 2). (to) Integrating over the base, to obtain one-third of the kinetic energy T, 3T = 2 pf '3 4R2(3x4-h4)dx/h 3 = pR2h4 / 1 35 V 3 (II) so that the effective k 2 of the liquid filling the trianglc is given by k 2 = T/Z p R 2 A = 2h2/45 = (radius of the inscribed circle) 2, (12) or two-fifths of the k 2 for the solid triangle.

Again, since d4)/dv =d /ds, d4)/ds= - d4y/dv, (13) T = 1 p f(1 9 d = - 2 p f4' d (14) With the Stokes' function ,y for motion symmetrical about an axis.

T =140 y d 27r y ds =7rp f c, 5 4 . (15) 37. Flow, Circulation, and Vortex Motion.-The line integral of the tangential velocity along a curve from one point to another, defined by s v as + u'a s) ds =f (udx+vdy-}-zdz), (I) is called the " flux " along the curve from the first to the second point; and if the curve closes in on itself the line integral round the curve is called the " circulation " in the curve.

With a velocity function 49, the flow -f d 4 = 4)142, (2) (9) (to) (6) (22) Z Uy (I -a4,ic /r4), so that the flow is independent of the curve for all curves mutually reconcilable; and the circulation round a closed curve is zero, if the curve can be reduced to a point without leaving a region for which 4 is single valued.

If through every point of a small closed curve the vortex lines are drawn, a tube is obtained, and the fluid contained is called a vortex filament. By analogy with the spin of a rigid body, the component spin of the fluid in any plane at a point is defined as the circulation round a small area in the plane enclosing the point, divided by twice the area. For in a rigid body, rotating about Oz with angular velocity the circulation round a curve in the plane xy is x ds yds) ds = times twice the area.

,In a fluid, the circulation round an elementary area dxdy is equal to dv du udx + (v+dx) dy- (u+dy) dx-vdy= () dxdy, so that the component spin is dv du (5) 2 dx - dy) in the previous notation of § 24; so also for the other two components and n.

Since the circulation round any triangular area of given aspect is the sum of the circulation round the projections of the area on the coordinate planes, the composition of the components of spin,, 7,, is according to the vector law. Hence in any infinitesimal part of the fluid the circulation is zero round every small plane curve passing through the vortex line; and consequently the circulation round any curve drawn on the surface of a vortex filament is zero.

If at any points of a vortex line the cross-section ABC, A'B'C' is drawn of the vortex filament, joined by the vortex line AA', then, since the flow in AA' is taken in opposite directions in the complete circuit ABC AA'B'C' A'A, the resultant flow in AA' cancels, and the circulation in ABC, A'B'C' is the same; this is expressed by saying that at all points of a vortex filament wa is constant where a is the cross-section of the filament and w the resultant spin (W. K. Clifford, Kinematic, book iii.).

So far these theorems on vortex motion are kinematical; but introducing the equations of motion of § 22, Du + dQ =o , Dv+dQ =o, Dw + dQ dt dx dt dy dt dz and taking dx, dy, dz in the direction of u, v, w, and dx: dy: dz=u: v: w, (udx + vdy + wdz) = Du dx +u 1+.. = -dQ+1dg2, and integrating round a closed curve (udx+vdy+wdz) =0, and the circulation in any circuit composed of the same fluid particles is constant; and if the motion is differential irrotational and due to a velocity function, the circulation is zero round all reconcilable paths. Interpreted dynamically the normal pressure of the surrounding fluid on a tube cannot create any circulation in the tube.

The circulation being always zero round a small plane curve passing through the axis of spin in vortical motion, it follows conversely that a vortex filament is composed always of the same fluid particles; and since the circulation round a cross-section of a vortex filament is constant, not changing with the time, it follows from the previous kinematical theorem that aw is constant for all time, and the same for every cross-section of the vortex filament.

A vortex filament must close on itself, or end on a bounding surface, as seen when the tip of a spoon is drawn through the surface of water.

Denoting the cross-section a of a filament by dS and its mass by dm, the quantity wdS/dm is called the vorticity; this is the same at all points of a filament, and it does not change during the motion; and the vorticity is given by w cos edS/dm, if dS is the oblique section of which the normal makes an angle e with the filament, while the aggregate vorticity of a mass M inside a surface S is M - l fw cos edS.

Employing the equation of continuity when the liquid is homogeneous, 2 (cly - d z) ? 2,. .., .. .,' d x 2 dy e d z2' (10) which is expressed . by p 2 (u, v, w) =2 curl (E, 1 7, 0, (, q, = 2 curl (u, v, w). (II) 38. Moving Axes in Hydrodynamics. - In many problems, such as the motion of a solid in liquid, it is convenient to take coordinate axes fixed to the solid and moving with it as the movable trihedron frame of reference. The components of velocity of the moving origin are denoted by U, V, W, and the components of angular velocity of the frame of reference by P, Q, R; and then if u, v, w denote the components of fluid velocity in space, and u', v', w' the components relative to the axes at a point (x, y, z) fixed to the frame of reference, we have u =U +u' - yR +zQ, v =V +v -zP +xR, w=W +w -xQ +yP.

Now if k denotes the component of absolute velocity in a direction fixed in space whose direction cosines are 1, m, n, k=lu+mv+nw; (2) and in the infinitesimal element of time dt, the coordinates of the fluid particle at (x, y, z) will have changed by (u', v', w')dt; so that Dk dl, do dt dt dt dt + dtw +1 (?t +u, dx +v, dy +w, dz) +m (d +u dx + v dy +w' dz) dw ,dw +n (dt ?dx+v?dy +w dz ) But as 1, m, n are the direction cosines of a line fixed in space, dl= m R-n Q, d m = nP-lR an =1Q-mPdt dt ' dt ' so that Dk __ du, du ,du ,du dt l (dt -vR+ wQ+u + v dy + w dz ) +m(.. .) +n(.. .) = (X-jfidp) +m (Y-dy) +n (Zp 2), for all values of 1, m, n, leading to the equations of motion with moving axes.

When the motion is such that d? d? dcp d11, do d? u '= - dx -md x ,' - dy -m dy , w = - dz-mdz' as in § 25 (I), a first integral of the equations in (5) may be written dp V + 2q 2 - d - n dt +14-14) (dx + m dz) +(v-v') (+m) +(w - w ) (+m) =F(t), (7) in which d4, do, d? d? dt-(u)dy- (w-w) dz = d - (U-yR+zQ) dy - (V-zP+xR)d -(W-xQ+yP) d z (8) is the time-rate of change of 49 at a point fixed in space, which is left behind with velocity components u-u', v-v', w-w'.

In the case of a steady motion of homogeneous liquid symmetrical about Ox, where 0 is advancing with velocity U, the equation (5) of § 34 p/p +V +Zq'2-f(,P') = constant becomes transformed into P +V + 2- dy + 2U 2 -f(t +2Uy 2 ) = constant, = 1,t+4Uy2, subject to the condition, from (4) § 34, Y -2 V = - f ' (Y', y 2 2 +2Uy2). Thus, for example, with = 4Uy 2 (r 2 a 2 -I), r2 = x2 +y 2, (13) for the space inside the sphere r=a, compared with the value of ,i' in § 34 (13) for the space outside, there is no discontinuity of the velocity in crossing the surface.

Inside the sphere d I d'd rd? I ,' 2 dx (y dx) +dy U dy) so that § 34 (4) is satisfied, with f' (W') =1.0 a2, f (Y") = 2 U'a2; and (ro) reduces to `)(() P +v-3 U j _ S = constant; (16) this gives the state of motion in M. J. M. Hill's spherical vortex, advancing through the surrounding liquid with uniform velocity.

39. As an application of moving axes, consider the motion of liquid filling the ellipsoidal case 2 y 2 z2 Ti + b1 +- 2 = I; (1) and first suppose the liquid be frozen, and the ellipsoid l3 (4) (I) (6) (9) (I o) (II) (12) (14) = 2 U ¢ 2, (15) rotating about the centre with components of angular velocity, 7 7, f'; then u= - y i +z'i, v = w = -x7 7 +y (2) Now suppose the liquid to be melted, and additional components of angular velocity S21, 522, S23 communicated to the ellipsoidal case; the additional velocity communicated to the liquid will be due to a velocity-function 2224_ - S2 b c 6 a 5 x b2xy, as may be verified by considering one term at a time.

If u', v', w' denote the components of the velocity of relative to the axes, = u +yR - zQ =a2+ b2S23y - c2 a2 ? Z, v' -V --zP - x R b2 2b2 c2 S2 1 z a2+ 2b2S23x, w' = w -{-xQ-yP2c2 cl? 2c2 c2. +a2 522x - b2 + c2121y, P =S21+, Q = 522+ 7 ], R =S23-H".

Thus a 2 vb 2 ? 2 = o, (8) so that a liquid particle remains always on a similar ellipsoid.

The hydrodynamical equations with moving axes, taking into account the mutual gravitation of the liquid, become dp +4 p Ax+ du - vR {-wQ? u ,du +v ,du +w ,du =o,...,..., (9) ax '? d t dy dz where _ oo abcdA A, B ' C ' - (a 2 +A, b 2 ±x, A, c 2 +A) P P 2 = 4(a 2 -F-A) (b 2 ±A) (c2+A).

With the values above of u, v, w, u', v', w', the equations become of the form p x + 4 7rpAx -Fax -{-hy-}-gz =o, - p - dy+ 4?pBy + hx+ay+fz =o, P d p + TpCZ +f y + yz = o, and integrating p p 1+27rp(Ax2+By2+CZ2) +z ('ax e +ay e + yz2 2 f yz + 2gzx + 2 hx y) = const., (14) so that the surfaces of equal pressure are similar quadric surfaces, which, symmetry and dynamical considerations show, must be coaxial surfaces; and f, g, h vanish, as follows also by algebraical reduction; and 4c2 (c 2 - a2) ? /c - a2, 1 z a = (c2 +a2)2 "2 +a21 J z z 2 ? z 2 a - 1 / 4 (a 2 +b 2) 2) ? 3 C a ' 2 +b 2 ? 3 /, with similar equations for 13 and If we can make (41 rpA + a) x2' = (4 7 P B +a) b2 = (4 7 P C +7)0, (16) the surfaces of equal pressure are similar to the external case, which can then be removed without affecting the motion, provided remain constant.

This is so when the axis of revolution is a principal axis, say Oz; when S21=0, t 2 =0, =o, o=0. (17) ellipsoid of liquid of three unequal axes, rotating bodily about the least axis;. and putting a=b, Maclaurin's solution is obtained of If 1 - =o or 8 3 =1" in addition, we obtain the solution of Jacobi's the rotating spheroid.

In the general motion again of the liquid filling a case, when a = b, 52 3 may be replaced by zero, and the equations, hydrodynamical and dynamical, reduce to d =- 2+ 2 J, = 2 x22111, d = 2 2`2 (+/'15-Om) (1 yy y n`t dt a +c dt a +c dt a +c) dc2, a2-1-c2 d122 a2 c2 dt ="2) +a2= G2y 71' dt = 121 1 - a 2 -c 2SJ, (19) of which three integrals are e +777 r z y 2= L -?2J2, (20) (a2 + c2) 2 2 121+14 =M+ 2c2(a2-c2)1 ' (21) 121+522hN = + x24 2,2 and then (dt / 2 = (a2 + c 2) 2(° v 2 - 12171) 2 4C4 2 2 - (+ c2)2?(E+77) (? 1 a +') - (121+52277)2] 4 - (a2+0)2 [L M -N2 { L 2c (a c 2 c 2) ae - N az+c2 l Y2 T L + 2 z 2 M (a2+c2) (9a2 - c2) 4 J 16c4 (a2-c2) = Z, where Z is a quadratic in ? 2, so that is an elliptic function of t, except when c =a, or 3a.

Put S2 1 =12 cos 4, 12 2 = -12 sin 4, d4 d52 1 dS22 Y a2+c2 122 7Ti = 71 22 CL2- c2(121+5221)J, a2 +c2 do a2+c2 + 4c2 z dt a'-c2 (a2+,c2)2 M+2c2(a2-c2 N-{-a2+c2 2 Ý_a 2 +c 2 (' 4c2 .?"d za 2 -c 2 c2)2 2'J Z M+ -c2) which, as Z is a quadratic function of i 2, are non-elliptic so also for ;G, where =co cos ,G, 7 7 = - sin 4.

In a state of steady motion d4- 121 _S22 Tit °' - fl 4=1G = nt, suppose, S21 -F9,277 = S2co, d4 a2+c2 WI- 1 a2-c2S21' _ 2a 2 SZ dt a2+c2cos' a 2 + c 2 a, 2 a 2 S2 I- a2_c22--a2+C2,0, 1a2 c2)2 (a 2 -c 2) (9a2-c2) ? 2 2a2 c2 4 (a2+L2) ' and a state of steady motion is impossible when 3a c >a. An experiment was devised by Lord Kelvin for demonstrating this, in which the difference of steadiness was shown of a copper shell filled with liquid and spun gyroscopically, according as the shell was slightly oblate or prolate. According to the theory above the stability is regained when the length is more than three diameters, so that a modern projectile with a cavity more than three diameters long should fly steadily when filled with water; while the old-fashioned type, not so elongated, would be highly unsteady; and for the same reason the gas bags of a dirigible balloon should be over rather than under three diameters long.

40. A Liquid Jet. - By the use of the complex variable and its conjugate functions, an attempt can be made to give a mathematical interpretation of problems such as the efflux of water in a jet or of smoke from a chimney, the discharge through a weir, the flow of water through the piers of a bridge, or past the side of a ship, the wind blowing on a sail or aeroplane, or against a wall, or impinging jets of gas or water; cases where a surface of discontinuity is observable, more or less distinct, which separates the running stream from the dead water or air.

Uniplanar motion alone is so far amenable to analysis; the velocity function 4 and stream function 1G are given as conjugate functions of the coordinates x, y by w=f(z), where z= x +yi, w=4-Plg, and then dw dod,y az = dx + i ax - -u+vi; so that, with u = q cos B, v = q sin B, the function - Q dw u_vi=g22(u-}-vi) = Q(cos 8+i sin 8), gives f' as a vector representing the reciprocal of the velocity in direction and magnitude, in terms of some standard velocity Q.

To determine the motion of a jet which issues from a vessel with plane walls, the vector I must be constructed so as to have a constant (to) (II) the liquid (15) 2, integrals; ,(29) (30) (I) direction 0 along a plane boundary, and to give a constant skin velocity over the surface of a jet, where the pressure is constant. It is convenient to introduce the function =log =log (Q / q) + 01 (4) so that the polygon representing 12 conformally has a boundary given by straight lines parallel to the coordinate axes; and then to determine S2 and w as functions of a variable u (not to be confused with the velocity component of q), such that in the conformal representation the boundary of the and w polygon is made to coincide with the real axis of u.

It will be sufficient to give a few illustrations.

Consider the motion where the liquid is coming from an infinite distance between two parallel walls at a distance xx' (fig. 4), and issues in a jet between two edges A and A'; the wall xA being bent at a corner B, with the external angle (3= 2Wr/n.

The theory of conformal representation shows that the motion is given by (b-a'.u -a) +? (b -a.u-a')1 m. [ - N -J (a -a'. u - b), u =ae ?, (5) where u =a, a' at the edge A, A l.; u at a corner B; u=o across xx' where 4 = oo; and u = oo, 43-= co across the end J J' of the jet, bounded by the curved lines APJ, A'P'J', over which the skin velocity is Q. The stream lines xBAJ, xA'J' are given by = 0, m; so that if c denotes the ultimate breadth JJ' of the jet, where the velocity may be supposed uniform and equal to the skin velocity Q, m=Qc, c=m/Q. If there are more B corners than one, either on xA or x'A', the expression for i is the product of corresponding factors, such as in (5) Restricting the attention to a single corner B, a = n(cos no +i sin 110) _ (b-a'.0-a) +1! (b-a.0-a) (6) (Q) AI (a-a .0b) ch n2= ch log (Q) n cos 114+i sh log (9) re sin n9 = 2(r+ fi n) = b - a ' ju -a (7) a-a' l u-b nf2 = sh log (cos nO +1 ch log (" sin 110 =2(?" - ? - n)= l b - au - ' (8) a - a u - b (9) dS2 I A I (b-a.b-a') dw m du = 21/(U - b)- ‘ 1 (u-a.0-a')' du - ,r u' Io) the formulas by which the conformal representation is obtained.

For the 2 polygon has a right angle at a=a, a', and a zero angle at u = b, where 0 changes from o to 27/n and 1 - 2 increases by 21-rr/n; so that dSt A (b -a.b -a') a - (u -b),/ (u - a.u - a where A= tn (II) And the w polygon has a zero angle at u =o, oo, where 4, changes from o to m and back again, so that w changes by im, and du =B, where B=-. (12) Along the stream line xBAPJ, t ' =0, u=ae-" c bl ,n; (13) and over the jet surface JPA, where the skin velocity is Q, - q = - Q, u = ae rs Q /m = ae rs lc, (14) ds denoting the arc AP by s, starting at u = a; a ' ch nS2=cos nB= -a' u u - - a b' (15) a l a - b l u - a' a-a' u-b' co > u = ae'" S " c > a, and this gives the intrinsic equation of the jet, and of curvature ds '&1) _ i dw i dw dS2 P= - dO = Q a0 - Q as2 = Q c u-b d (u -a.u -a') _ ? . 2n u (a -b.b -a')' not requiring the integration of (II) and (12) If 0=a across the end J J' of the jet, where u = oo, q= Q, b-a' a-b ch 7/2,=cos na = I ,, sh 162 = i sin na = i,, a - a a-a Then a-a'+(a+a) cos 2na-[a+a'+(a-a) cos 2najcos 2110 (a-a') sin' 2na X cos 211a - cos 2110 Along the wall AB, cos nO =o, sin n0= I, a> u> b, ch n2= i sh log (q) n= i? b - a ? `ia a - a u-b sh nf2= i ch log (Q) n =i a-b ? Iju -a' a-u a-b ds ds d? m _ c Q du dO dt - rrqu 2r qu AB _ Q du L bq (u -a') +V (b -a')1,l (a -u)11/ndu) L 1,1 (a-a/),,1 (u-b) ,f u Along the wall Bx, cos n0 =I, sin n0 =o, b >u>o ch nSt = ch log () n= ,, fb-a?, ? Ja - u  ?I a -a b -u' sh nS2=sh log (Q)= ?a - b a - a' b - u' At x where = co, u = o, and q= go, (O n b - a ' a + a -b a' cio) - ?a-a'?b a-a' q In crossing to the line of flow x'A'P'J', b changes from o to m, so that with q = Q across JJ', while across xx the velocity is qo, so that i n = go. xx' = Q.JJ' (31) JJ _ g o _ lb-a' la_I a-b a' 11n xx' Q L Va - a' b Va-a' b j ' (32) giving the contraction of the jet compared with the initial breadth of the stream.

Along the line of flow x'A'P'J', >L =m, u=a'e-rk/"n, and from x' to 41. The first problem of this kind, worked out by H. v. Helmholtz, of the efflux of a jet between two edges A and A 1 in an infinite wall, is obtained by the symmetrical duplication of the above, with n = 1, b = o, a' = - oo, as in fig. 5, I ch S2 = u a, sh C2= ' u (I) (I and along the jet APJ, oo > u=aerslc>a, sh S2=i sin 0 =iu=ie zrs/o, (2) PM sin 0 ds = f e ds = 1 = 1 sin 0 , (3) cos 272a - cos 2n0 = 2Q - ?ib L a b2 s i n' 27ta u-b A (a- (u -a.u -b') sin 2110 - 2 a-a .u-b  ?l (u -a.u -a') = s in 2na u-b 2n b) A (ab.ba') p l u -bJ (u -a.u -a') sh nS2=i sin 110=i then the radius sin 2170 (30) A', cos nO= i, sin n0=o, n 1 ' b-a' ch nS2= ch log (9) = Va -a' n shnS2= shlog (Q) q _ o> u>a'. Along the jet surface A'J', q = Q, b-a' ch nSl= cos 110= a-a la - b sh nft=i sin nO=i a'>u=a'erl"> -oo, giving the intrinsic equation.

rf > a> b> o> a'> -Do; and then so that PT =c/Zir, and the curve AP is the tractrix; and the coefficient of contraction, or breadth of the jet breadth of the orifice - +i' A change of S2 and 0 into nS2 and nO will give the solution for two walls converging symmetrically to the orifice AA 1 at an angle zr/n. With n=1, the re-entrant walls are given of Borda's mouthpiece, and the coefficient of contraction becomes 2. Generally, by making a' = -oo, the line x'A' may be taken as a straight stream line of infinite length, forming an axis of symmetry; and then by duplica tion the result can be ob A tained, with assigned n, a, and b, of the efflux from a symmetrical converging FIG. 5. FIG. 6. mouthpiece, or of the flow of water through the arches of a bridge, with wedge-shaped piers to divide the stream.

42. Other arrangements of the constants n, a, b, a' will give the results of special problems considered by J. M. Michell, Phil. Trans. 1890.

Thus with a' =o, a stream is split symmetrically by a wedge of angle ' zr/n as in Bobyleff's problem; and, by making a = oo, the wedge extends to infinity; then chnS2= u ,sh nS2= b n u. (I) Over the jet surface 4'=m, q=Q, u=-e rr,lm= -berslc, ch SZ=cos n0= e>rsle+I, shS2 =i sin ins =tan ds 2n (3) e2 =tan nO, - c dB sin 2,10' For a jet impinging normally on an infinite plane, as in fig. 6, n ei„s /c =tan 0, ch (z7rs/c) sin 20 = I, (4) sh 27rx/c = cot 0, sh 27ry /c = tan B, sh 27rx/c sh irry/c = I, ei r (x+Y)l ei rxlc +ei.Y/c +I. (5) With n =1, the jet is reversed in direction, and the profile is the atenary of equal strength.

', In Bobyleff's problem of the wedge of finite breadth, ch nS2 = ?a b' s h n S 2 = V b a a 1 u u b , (6) ? b a-b (7) cos = a, sin na = 1j a, nd along the free surface APJ, q =Q, 4) =o, u =e-. 4l m =aeffsle, e rs /° - I u =ae-"' - (a -b)w', (io) w'+ ch nS2 =, sh nS2 = w,, (II) I n which we may write w' =4)+41i. (12) Along the stream line xABPJ, 4) =o; and along the jet surface PJ, -1 >49> - oo; and putting 4 = -irs/c - I, the intrinsic quation is irs/c =cot 2 nO, (13) hich for n =I is the evolute of a catenary.

43. When the barrier AA' is held oblique to the current, the stream line xB is curved to the branch point B on AA' (fig. 7), and so must be excluded from the boundary of u; the conformal re presentation is made now with du= (b-a.b-a') du - (u-b) A l (u-a.0-a) (I) dw m I m' du = 7r u-j - u -j' _ m+m' u-b it u' j.0-j" b = mj i m'j m+m', taking u = co at the source where FIG.7.

0-00, u = b at the branch point B, u = j, j at the end of the two diverging streams where = -oo; while ¢=0 along the stream line which divides at B and passes through A, A'; and 4 ' =m, -m' along the outside boundaries, so that m/Q, m'/Q is the final breadth of the jets, and (m+m')/Q is the initial breadth, c, of the impinging stream. Then b - a' a -a ,b - a u-a' ch 252= a -a'? u -b' s h 2S2 =1 1a-a ? u -b' Along a jet surface, q=Q, and ch S2= cos 0 =cos a-i sin2a(a-a')/(u-b), (5) if 0 =-a at the source x of the jet xB, where u = co; and supposing 0=0,13 at the end of the streams where u =j, j', u-b i sin 2 a u - j cos 0-cos /3 i a -a cos a sin a -cos 0' aa' - 2 (cos a -cos (3) (cos a-cos 0)' u-j' 1 2 cos 0-cos, (6) a -a' - 2 S i n a (cos a -cos (3') (cos a -cos B)' and 4' being constant along a stream line d4 - dw ds _d8 d4 _ dw du du du' d- -dud0' 7rQ ds_ it ds (cos a-cos /3) (cos a -cos (3') sin 0 m+m' dB c d0 - (cos a-cos B) (cos 0-cos /3) (cos 0 -cos /3')' _ sin 0 cos a-cos 13 sin 0 - cos a-cos B + cos 0-cos (3' cos 0-cos 13 cos a -cos $ sin 6 cos (3-cos /3' cos 0-cos 0" giving the intrinsic equation of the surface of a jet, with proper attention to the sign.

From A to B, a>u >b, 0=0, ch S2= ch log Q=cos a-i sin 2a a-b I sh S2= sh log Q= I (a u-b-a/) s i n a Q = (u-b) cos a-2(a-a') sin 2 a+1,/ (a-u.u- a')sin a (8) u-b ds _ ds d4 _ Q dw Q du - Q d 4) du q du (u-b) cos a-2(a- a') sin 2 a (a-u.0 - a') sin a (9) it j- -j' AB _f a(2b - a - a')(u-b)-2(a-b)(b-a')+2V (a - b. b - a'. a - u. u - a') du, () - a - a' .j - u.0 - j' IO a with a similar expression for BA'.

The motion of a jet impinging on an infinite barrier is obtained by putting j = a, j' = a'; duplicated on the other side of the barrier, the motion reversed will represent the direct collision of two jets of unequal breadth and equal velocity. When the barrier is small compared with the jet, a=0=0', =a', and G. Kirchhoff's solution is obtained of a barrier placed obliquely in an infinite stream.

Two corners B 1 and in the wall xA, with a' = -00, and n =I, will give the solution, by duplication, of a jet issuing by a reentrant mouthpiece placed symmetrically in the end wall of the channel; or else of the channel blocked partially by a diaphragm across the middle, with edges turned back symmetrically, problems discussed by J. H. Michell, A. E. H. Love and M. Rethy.

ch S2 = 2b - a ?a' N a - a u - b ' sh S2 =1/ N V (2. au. u-al) u-b a -a' N =2 a -a' ' ? 8 (2) (4) e ,rs /c e ns/c + I' (2) cos n0= cos na-N e' 31 ' - cos'na' cos 2 na sin2n0 (8) sin 2 n0 - sin2na' he intrinsic equation, the other free surface A'P'J' being given by e m /? - cos 2 na sin2ng s (9) sin 2 na - sin2n0 Putting n =I gives the case of a stream of finite breadth disturbed y a transverse plane, a particular case of Fig. 7.

When a = b, a = o, and the stream is very broad compared with he wedge or lamina; so, putting w=w' (a-b)la in the penultimate ase, and A B A (7) When the polygon is closed by the walls joining, instead of reaching back to infinity at xx', the liquid motion must be due to a source, and this modification has been worked out by B. Hopkinson in the Proc. Lond. Math. Soc., 1898.

Michell has discussed also the hollow vortex stationary inside a polygon (Phil. Trans., 1890); the solution is given by ch nS2=sn w, shnS2=i cn w (II) so that, round the boundary of the polygon, ik = K', sin n8 =o; and on the surface of the vortex 1P= o, q = Q, and cos n8=sn4p,nB= Zit -am sic, (12) the intrinsic equation of the curve.

This is a closed Sumner line for n =I, when the boundary consists of two parallel walls; and n= z gives an Elastica.

44. The Motion of a Solid through a Liquid. - An important problem in the motion of a liquid is the determination of the state of velocity set up by the passage of a solid through it; and thence of the pressure and reaction of the liquid on the surface of the solid, by which its motion is influenced when it is free.

Beginning with a single body in liquid extending to infinity, and denoting by U, V, W, P, Q, R the components of linear and angular velocity with respect to axes fixed in the body, the velocity function takes the form = Ucb1+V42+W43+ P xi+Qx2+Rx3, (I) where the 0's and x's are functions of x, y, z depending on the shape of the body; interpreted dynamically, C -p0 represents the impulsive pressure required to stop the motion, or C +p4) to start it again from rest.

The terms of 0 may be determined one at a time, and this problem is purely kinematical; thus to determine 4)1, the component U alone is taken to exist, and then 1, m, n, denoting the direction cosines of the normal of the surface drawn into the exterior liquid, the function 01 must be determined to satisfy the conditions v 2 0 1 = o, throughout the liquid; (ii.) ' = -1, the gradient of 0 down the normal at the surface of the moving solid; 1 =0, over a fixed boundary, or at infinity; similarly for 02 and 03.

To determine x i the angular velocity P alone is introduced, and the conditions to be satisfied are (i.) 0 2 x1 = o, throughout the liquid; y l =mz - ny, at the surface of the moving body, but zero over a fixed surface, and at :infinity; the same for x 2 and x3. For a cavity filled with liquid in the interior of the body, since the liquid inside moves bodily for a motion of translation only, 41 = - x, 42 = - , 43 = - z; (2) but a rotation will stir up the liquid in the cavity, so that the'x's depend on the shape of the surface.

The ellipsoid was the shape first worked out, by George Green, in his Research on the Vibration of a Pendulum in a Fluid Medium (2833); the extension to any other surface will form an important step in this subject.

A system of confocal ellipsoids is taken y2 (3) a 2 +X b 2 +X c2 + A= I, and a velocity function of the form = x1 P, (4) where 4' is a function of X only, so that 4) is constant over an ellipsoid; and we seek to determine the motion set up, and the form of >G which will satisfy the equation of continuity.

Over the ellipsoid, p denoting the length of the perpendicular from the centre on a tangent plane, px _ pv _ _ pz 1= a2+X' b +A' n c2+A p2x2 + p2y2 p2z2 I (a2 - + X)2 (b 2 +x)2 + (0+X)2, p 2 = (a2+A)12+(b2+X)m2+(c2+X)n2, = a 2 1 2 +b 2 m 2 +c 2 n 2 +X, 2p d = ds; (8) Thence d? = dx ?+xd%y ds ds ds ds +2 l dd, so that the velocity of the liquid may be resolved into a component -41 parallel to Ox, and -2(a 2 +X)ld4/dX along the normal of the ellipsoid; and the liquid flows over an ellipsoid along a line of slope with respect to Ox, treated as the vertical.

Along the normal itself d= +2(a2+X) d ? l ' so that over the surface of an ellipsoid where X and ¢ are constant, the normal velocity is the same as that of the ellipsoid itself, moving as a solid with velocity parallel to Ox U = -q, - 2 (a2+X) dtP, and so the boundary condition is satisfied; moreover, any ellipsoidal surface X may be supposed moving as if rigid with the velocity in (I I), without disturbing the liquid motion for the moment.

The continuity is secured if the liquid between two ellipsoids X and X 11 moving with the velocity U and 15 1 of equation (II), is squeezed out or sucked in across the plane x=o at a rate equal to the integral flow of the velocity I across the annular area a l. - a of the two ellipsoids made by x=o; or if aU - a ? ' d?dA, a= 7r- % 1 (b2+a.c2+A). Expressed as a differential relation, with the value of c-TKL ai,G+2 (a 2 +A) a d ? J -1k di - d 3a dX +2(a2+X)d (a -) =o, and integrating ( a 2 + X) 3 /2ad? = a constant, so that we may put MdX (17) (a2+X)P' P2= 4(a2 + X)(b2 +X)(c 2 +X), (18) where M denotes a constant; so that 4) is an elliptic integral of th second kind.

The quiescent ellipsoidal surface, over which the motion is entirely tangential, is the one for which (a2+X)d? +4) =0, (19) and this is the infinite boundary ellipsoid if we make the upper limi =co.

The velocity of the ellipsoid defined by X =o is then U= - 2 __ M ((ro b J o (a2 =ab (i -A0), (20) with the notation A or A a a= a (a2bc+ = - 2abc d -- so that in (4) xA x 'UxA x A' 4)' 1 -Ao' (22) in (I) for an ellipsoid.

The impulse required to set up the motion in liquid of density p i the resultant of an impulsive pressure p4) over the surface S of th ellipsoid, and is therefore ffp4ldS = p4GoffxldS =p 40 (volume of the ellipsoid) =4)oW', (23) where W' denotes the weight of liquid displaced.

Denoting the effective inertia of the liquid parallel to Ox by aW' the momentum aW'U = 4)0W' (24) _ U i -AO' 25) in this way the air drag was calculated by Green for an ellipsoida pendulum.

Similarly, the inertia parallel to Oy and Oz is NW' - 1 B W', B C (b2 +-X, c 2 ab and A +C abc/ZP, Ao For a sphere a=b=c, Ao= Bo=Co =, 'a' = Q = = z, (9) U from (II), (16) so that the effective inertia of a sphere is increased by half the weight of liquid displaced; and in frictionless air or liquid the sphere, of weight W, will describe a parabola with vertical acceleration W - W' ,g (30) W+ aW Thus a spherical air bubble, in which W/W' is insensible, will begin to rise in water with acceleration 2g.

45. When the liquid is bounded externally by the fixed ellipsoid A = A I, a slight extension will give the velocity function 4 of the liquid in the interspace as the ellipsoid A=o is passing with velocity U through the confocal position; 4 must now take the formx(1'+N), and will satisfy the conditions in the shape CM abcdX ¢ = Ux - Ux a b x 2+X)P Bo+CoB I - C 1 (A 1 abcdX, I a1b1cl - J o (a2+ A)P and any'confocal ellipsoid defined by A, internal or external to A=A 1, may be supposed to swim with the liquid for an instant, without distortion or rotation, with velocity along Ox BA+CA-B 1 -C1 W'.

ZUy2BB0 Bll; reducing, when the liquid extends to infinity and B 3 =0, to = xA o' _ - zUy 2B o so that in the relative motion past the body, as when fixed in the current U parallel to xO, A 4)'=ZUx(I+Bo), 4)'= zUy2(I-B o)(6) Changing the origin from the centre to the focus of a prolate spheroid, then putting b 2 =pa, A = A'a, and proceeding to the limit where a = oo, we find for a paraboloid of revolution P B - p (7) B = 2p +A/' Bo p+A y2 i =p+A'- 2x, (8) p+?

with A' =0 over the surface of the paraboloid; and then' = ZU[y 2 - pJ (x2 + y2) + px ]; (9) =-2U p [1/ ( x2 + y2)-x]; (io) 4, = - ZUp log [J(x2+y2)+x](II) The relative path of a liquid particle is along a stream line 1,L'= 2Uc 2, a constant, (12) = /,2 3 ,2 _ (y 2 _ C 2) 2 2 2 2' - C2 2 x 2p(y2 - c2) /' J(x2 +y 2)= py ` 2p(y2_c2)) (13) a C4; while the absolute path of a particle in space will be given by dy_ r - x _ y 2 - c2 dx_ - y - 2py y 2 - c 2 = a 2 e -x 1 46. Between two concentric spheres, with a 2 +A = r 2, a2+A1=a12, A=B =C =a3/3r3, a 3 a 3 a3 _ a3 Cb _1U x r 3 a13 Y'=2 Uy2 r3 3 a13 .

2 V I - a /al ' Y' I-a /al ' and the effective inertia of the liquid in the interspace Ao+2A1 W, =1 a13 +2a3W'. 2A 0 - 2AI 2 a 1 3 - a3 When the spheres are not concentric, an expression for the effective inertia can be found by the method of images (W. M. Hicks, Phil. Trans., 1880).

The image of a source of strength p at S outside a sphere of radius a is a source of strength pa/f at H, where 'OS' =f, OH =a2/f, and a line sink reaching from the image H to the centre 0 of line strength - A la; this combination will be found to produce no flow across the surface of the sphere.

Taking Ox along OS, the Stokes' function at P for the source S is p cos PSx, and of the source H and line sink OH is p(a/f) cos PHx and - (p/a) (PO - PH); so that = p (cos PSx+f cos PHx PO a PH), (q) and Ili = -p, a constant, over the surface of the sphere, so that there is no flow across.

When the source S is inside the sphere and H outside, the line sink must extend from H to infinity in the image system; to realize physically the condition of zero flow across the sphere, an equal sink must be introduced at some other internal point S'.

When S and S' lie on the same radius, taken along Ox, the Stokes' function can; be written down; and when S and S' coalesce a doublet is produced, with a doublet image at H.

For a doublet at S, of moment m, the Stokes' function is M f cos PSx = - m p s3; and for its image at H the Stokes' function is m f cos PHx =m f 3 PH" (6) so that for the comnation _ a3 I I 2 4)-myb12 (f 3 PH PS 3) =m f 3 (pa ll 3 P53)' 3 and this vanishes over the surface of the sphere.

There is no Stokes' function when the axis of the doublet at S does not pass through 0; the image system will consist of an inclined doublet at H, making an equal angle with OS as the doublet S, and of a parallel negative line doublet, extending from H to 0, of moment varying as the distance from O.

A distribution of sources and doublets over a moving surface will enable an expression to be obtained for the velocity function of a body moving in the presence of a fixed sphere, or inside it.

The method of electrical images will enable the stream function ,)' to be inferred from a distribution of doublets, finite in number when the surface is composed of two spheres intersecting at an angle 7r/m, where m is an integer (R. A. Herman, Quart. Jour. of Math. xxii.).

Thus for m =2, the spheres are orthogonal, and it can be verified that a13 a2 3 aY3 i f /' = ZU (I - 13 - 7.2 3 + 3) ' (8) where a l , a2, a =a l a 2 /J (a 1 2 +a 2 2) is the radius of the spheres and their circle of intersection, and r 1, r 2, r the distances of a point from their centres.

The corresponding expression for two orthogonal cylinders will be With a 2 = co, these reduce to / y /, = Uy (I ra 2 p22 +-C24).. (9) a 5 ` = 2 (I -) y a x, or Uy (1- - 2 Y a 4 4) a, (io) for a sphere or cylinder, and a diametral plane.

Two equal spheres, intersecting at 120°, will require - I U j x _ a 3 a4(a 7 2 x) a3 a4(a+2x)] (II) 2 - _ 2 y a 271 3 271 +2Y2 3 2720 ' with a similar expression for cylinders; so that the plane x=o may be introduced as a boundary, cutting the surface at 60°. The motion of these cylinders across the line of centres is the equivalent of a line doublet along each axis.

47. The extension of Green's solution to a rotation of the ellipsoid was made by A. Clebsch, by taking a velocity function 4,=xyx (I) for a rotation R about Oz; and a similar procedure shows that an ellipsoidal surface A may be in rotation about Oz without disturbing the motion if I I dx + _ a2'-A) x 2 a R t i/(b2+A)- i/(a2+A) and that the continuity of the liquid is secured if (a 2 _ I -A) 3/2 (b 2 4 A)3f2(C2 -+- A) 2 ?? = constant, _ ff 00 NdA N BA-AA X - JA (a' +X) (b 2 +A)P - abc' a2 -b2 ' and at the surface A = o, I I N Bo-A 0 N I R - (a2+b2) abc a 2 -b 2 abc a2b2 I /b 2 N = R I /b2 - I /a2 abc I 1 I Bo - AO' a 2 b 2 - a2 b2 a 2 b2 = R (a 2 - b 2) /(a 22 + /b2) 2 - r (B o - Ao) U Bo+Co - B I - CI' Since - Ux is the velocity function for the liquid W' filling the ellipsoid A = o, and moving bodily with it, the effective inertia of the liquid in the interspace is Ao+B1+C1 Bo+Co - B1 - C, If the ellipsoid is of revolution, with b=c, - 2 XBo - - C BI' and the Stokes' current function 4, can be written down (I) is (5) (7) (6) The velocity function of the liquid inside the ellipsoid A=o due to the same angular velocity will be = Rxy (a2 - b2)/(a2 + b2), (7) and on the surface outside _ N Bo -Ao c1)0xy abc 2 62' so that the ratio of the exterior and interior value of at the surface is ?o= Bo-Ao (9) 4)1 (a 2 -6 2)/(a2 + b) - (Bo - Ao)' and this is the ratio of the effective angular inertia of the liquid, outside and inside the ellipsoid X = o.

The extension to the case where the liquid is bounded externally by a fixed ellipsoid X= X is made in a similar manner, by putting 4 = x y (x+ 11), (io) and the ratio of the effective angular inertia in (9) is changed to 2 (B0-A0) (B 1A1) +.a12 - a 2 +b 2 a b1c1 a -b -b12 abc (Bo-Ao)+(B1-A1) a 2 + b 2 a1 2 + b1 2 alblcl Make c= CO for confocal elliptic cylinders; and then _, 2 A? ? (a2 + A)? (4 a2 +A.b'-'-f-A) - a 2a b (I - V a2+A), (12) ab (. C - a 2 - 2 b2 +A , = and then as above in § 31, with a= c ch a, b=c sh a, a =-1 (a 2 +X) =c ch al, b1= c sh a (13) the ratio in (II) agrees with § 31 (6).

As before in § 31, the rotation may be resolved into a shear-pair, in planes perpendicular to Ox and Oy.

A torsion of the ellipsoidal surface will give rise to a velocity function of the form 4)--- where SZ can be expressed by the elliptic integrals in a similar manner, since dX/P3.

48. The determination of the O's and x's is a kinematical problem, solved as yet only for a few cases, such as those discussed above.

But supposing them determined for the motion of a body through a liquid, the kinetic energy T of the system, liquid and body, is expressible as a quadratic function of the components U, V, W, P, Q, R. The partial differential coefficient of T with respect to a component of velocity, linear or angular, will be the component of momentum, linear or angular, which corresponds.

Conversely, if the kinetic energy T is expressed as a quadratic function of x, x x3, y1, y2, y3, the components of momentum, the partial differential coefficient with respect to a momentum component will give the component of velocity to correspond.

These theorems, which hold for the motion of a single rigid body, are true generally for a flexible system, such as considered here for a liquid, with one or more rigid bodies swimming in it; and they express the statement that the work done by an impulse is the product of the impulse and the arithmetic mean of the initial and final velocity; so that the kinetic energy is the work done by the impulse in starting the motion from rest.

Thus if T is expressed as a quadratic function of U, V, W, P, Q, R, the components of momentum corresponding are dT dT dT (I) = dU + x2=dV, x3 =dW, dT dT dT Yi dp' dQ' y3=dR; but when it is expressed as a quadratic function of xi, 'x2, x3, yi, Y2, Y3, U = d, V= dx ,' w= ax dT Q_ dT dT dy 1 dy2 dy The second system of expression was chosen by Clebsch and adopted by Halphen in his Fonctions elliptiques; and thence the dynamical equations follow X = dt x2 dy +x3 d Y = ..., Z ..., (3) = dt1 -y2 ?y - '2dx3+x3 ' M =.. ., N =

, (4) where X, Y, Z, L, M, N denote components of external applied force on the body.

These equations are proved by taking a line fixed in space, whose direction cosines are 1, then dt=mR-nQ,' d'-t = nP =lQ-mP. (5) If P denotes the resultant linear impulse or momentum in this direction P =lxl+mx2+nx3, ' dP dt xl+, d y t x2' x3 +1 dtl dt 2 +n dt3, =1 ('+m (dt2-x3P+x1R) ' +n ('-x1Q-{-x2P) ' '= IX +mY+nZ, / (7) for all values of 1, Next, taking a fixed origin and axes parallel to Ox, Oy, Oz through 0, and denoting by x, y, z the coordinates of 0, and by G the component angular momentum about 1"2 in the direction (1, G =1(yi-x2z+x3y) m 2-+xlz) n(y(y 3x 1 x3x y + x 2 x)(8) Differentiating with respect to t, and afterwards moving the fixed. origin up to the moving origin 0, so that dy x=y=z=o, but dt U, dt= ' dG _ dyl =l (- yi y3Q x2w+xiv) +m (dY2yP+Yrxu+xw) +n (? 3 - y1Q+y2P-x1V+x2U =1L+mM+nN, (9) for all values of 1, m, n.

When no external force acts, the case which we shall consider, there are three integrals of the equations of motion (i.) T =constant, x 2 +x 2 +x 2 =F 2, a constant, (iii.) x1y1 +x2y2+x3y3 =n = GF, a constant; and the dynamical equations in (3) express the fact that x, x, xs. are the components of a constant vector having a fixed direction; while (4) shows that the vector resultant of y, y, y moves as if subject to a couple of components x Wx V, x Ux W, x V-x U, (Io) and the resultant couple is therefore perpendicular to F, the resultant of x, x, x, so that the component along OF is constant, as expressed by (iii).

If a fourth integral is obtainable, the solution is reducible to quadrature, but this is not possible except in a limited series of cases, investigated by H. Weber, F. Kotter, R. Liouville, Caspary, Jukovsky, Liapounoff, Kolosoff and others, chiefly Russian mathematicians; and the general solution requires the double-theta hyperelliptic function.

49. In the motion which can be solved by the elliptic function, the most general expression of the kinetic energy was shown by A. Clebsch to take the form T= 2p(x12 +x22)+2p'x32 + q (xiyi +x2y2) +q'x3y3 +2r(y12+y22)+2r'y32 so that a fourth integral is given by dy 3 /dt = o, y = constant; dx3 (4 y) (q + y) _ (y y) dt - xl 'x2 xl Y Y x l 2 - 1, y2 () = (x12 +x22) (y12 + y22) = (X 1 2 + X 2) +y22)-(FG-x3y3)2 = (x 1 y32-G2)-(Gx3-Fy3) 2, in which 2 = F 2 -x3 2, x l y l +x2y2 = FG-x3y3, Y(y1 2 +y2 2) = T -p(x12 +x22) -p'x32 -2q(xiyi 'x2y2)- 2 q ' x = (p -p') x 2 + 2 ( - q ') x 3 y 3+ m 1, (6) m1 = T 2 i y 3 2 (7) so that dt3) 2 =X3, (8) where X3 is a quartic function of x3, and thus t is given by an elliptic (8) (6) (I) integral of the first kind; and by inversion x 3 is in elliptic function of the time t. Now (x1 - x21) (y 1 +y21) = xl l +x2y2 + - (' r 1 2 - x2y1) = FG-x3y3+iV X3, yi+3 7 21_FG-x3y3+2V X3 xl+x21 X12 +X22 (x 1 +x 2 i) = - i{(q' - q)x3+r'y3]+irx3(y1+y21), = FG - x3y3 +ZJ X3 dt2log(x1+x22) - - (q g) x 3- r y3+rx3 F2x32 (12) d dl2 log V x1 ± x2 2 (q'-q)x3-(r'-r) y3FrFF2-x 2 3 ' (13) requiring the elliptic integral of the third kind; thence the expression of x1-f -x21 and yl-}-y21.

Introducing Euler's angles 0, c15, x1= F sin 0 sin 0, x 2 =F sin 0 cos 0, xl+x 2 i =iF sin 0e_ , x 3 = F cos 0; sin o t=P sin 4+Q cos 0, dT F sin 2 0d l - dy l + dy 2x = (qx1+ryi)xl +(qx2+ry2)x2 = q (x1 2 +x2 2) +r (xiyi +x2y2) = qF 2 sin 2 0-Fr (FG - x 3 y 3), (16) _Ft (FG _x 323 Frdx3 (17) F x3 X3 elliptic integrals of the third kind.

Employing G. Kirchhoff's expressions for X, Y, Z, the coordinates of the centre of the body, FX=y 1 cos xY--y 2 cos yY-{-y 3 cos zY, (18) FY = -y l cos xX -Hy2 cos yX+y 3 cos zX, (Ig) G=y 1 cos xZ+y 2 cos yZ+y 3 cos zZ, (20) (21) F(X+Yi) = Fy3-Gx3+i /) X 3epi. (22) Y (F2 x2) Suppose x 3 -F is a repeated factor of X3, then y 3 = G, and X 3 = (x 3 -F)2 [P' _ P(X3+F)2+2' _ G(X +F) -G 2 ], (23) nd putting x3-F=y, (y) 2= 7'3'2- [41' r 1' F 2 -{-4 g r qFG - G2 +2 (2P'r 19F+9 r q G) y+ r y (24) o that the stability of this axial movement is secured if A = 4 P' r ?'F 2 + 4 Y q FG - G 2 (25) s negative, and then the axis makes r J l (-A)/7r nutations per second. therwise, if A is positive rt= J y-s1 (A+2By+Cy') dy sh1 A'/ (A+2By+Cy 2) I ch1 A+By (26) -V A ch1 31, (B2--AC) - A sh - 1 (B2-AC)' nd the axis falls away ultimately from its original direction.

A number of cases are worked out in the American Journal of athematics (1907), in which the motion is made algebraical by the se of the pseudo-elliptic integral. To give a simple instance, hanging to the stereographic projection by putting tan 20=x, ill give a possible state of motion of the axis of the body; and the otion of the centre may then be inferred from (22).

50. The theory preceding is of practical application in the vestigation of the stability of the axial motion of a submarine oat, of the elongated gas bag of an airship, or of a spinning rifled rojectile. In the steady motion under no force of such a body in medium, the centre of gravity describes a helix, while the axis escribes a cone round the direction of motion of the centre of ravity, and the couple causing precession is due to the dislacement of the medium.

In the absence of a medium the inertia of the body to transtion is the same in all directions, and is measured by the (3) But the change of the resultant momentum F of the medium as. well as of the body from the vector OF to O'F' requires an impulse couple, tending to increase the angle F00', of magnitude, in sec. foot-pounds F.00'.sin FOO'=FVt sin (0-0), (4) equivalent to an incessant couple N=FV sin (0-0) = (F sin 0 cos 0-F cos 0 sin ¢)V = (c 2 -c i) (V /g) sin 0 cos 4) =W'(13-a)uv/g (5) This N is the couple in foot-pounds changing the momentum of the medium, the momentum of the body alone remaining the same; the medium reacts on the body with the same couple N in the opposite direction, tending when c 2 -c 1 is positive to set the body broadside to the advance.

An oblate flattened body, like a disk or plate, has c 2 -c 1 negative, so that the medium steers the body axially; this may be verified by a plate dropped in water, and a leaf or disk or rocket-stick or piece of paper falling in air. A card will show the influence of the couple N if projected with a spin in its plane, when it will be found to change its aspect in the air.

An elongated body like a ship has c 2 -c 1 positive, and the couple N tends to disturb the axial movement and makes it unstable, so that a steamer requires to be steered by constant attention at the helm.

Consider a submarine boat or airship moving freely with the direction of the resultant momentum horizontal, and the axis at a slight inclination 0. With no reserve of buoyancy W =W', and the couple N, tending to increase 0, has the effect of diminishing the metacentric height by h ft. vertical, where Wh tan 0 = N = (c 2 - cl) c2 g2 tan 0, (6) (7) in which we have put k' 2 = ek 2, where E is a numerical factor depending on the shape.

(Nxe, fri) 3 / 2 =(x+1) 1 / X2, X1 _ =ax 4 +2ax 3 =3(a+b)x 2 +2bx =b, X2 N3= - 8(a+b), weight W, and under no force the C.G. proceeds in a straight line, and the axis of rotation through the C.G. preserves its original direction, if a principal axis of the body; otherwise the axis describes a cone, right circular if the body has uniaxial symmetry, and a Poinsot cone in the general case.

But the presence of the medium makes the effective inertia depend on the direction of motion with respect to the external shape of the body, and on W' the weight of fluid medium displaced.

Consider, for example, a submarine boat under water; the inertia is different for axial and broadside motion, and may be represented by (1) c 1 =W+W'a, c2=W+W'/3' where a, R are numerical factors depending on the external shape; and if the C.G is moving with velocity V at an angle 4) with the axis, so that the axial and broadside component of velocity is u = V cos 0, v =V sin 4), the total momentum F of the medium, represented by the vector OF at an angle 0 with the axis, will have components, expressed in sec. Ib, F cos 0 =c 1 - = (W +W'a) V cos 43, F sin 0 = c 2.11 = (W +W'/3) V sin 4) . (2) g g a g Suppose the body is kept from turning as it advances; after t seconds the C.G. will have moved from 0 to 0', where 00' = Vt; and at 0' the momentum is the same in magnitude as before, but its vector is displaced from OF to O'F'.

For the body alone the resultant of the components of momentum W V -cos andW V sin 0 is W V -sec. lb, acting along 00', and so is unaltered.

(9) (io) (II) F 2 (X 2 +Y 2) = y12+y22 +y32G2, c 2 - c 1 +a u2 h W c2 g = (? - +13 g 51. An elongated shot is made to preserve its axial flight through the air by giving it the spin sufficient for stability, without which it would turn broadside to its advance; a top in the same way is made to stand upright on the point in the position of equilibrium, unstable statically but dynamically stable if the spin is sufficient; and the investigation proceeds in the same way for the two problems (see Gyroscope).

The effective angular inertia of the body in the medium is now required; denote it by C 1 about the axis of the figure, and by C2 about a diameter of the mean section. A rotation about the axis of a figure of revolution does not set the medium in motion, so that C 1 is. the moment of inertia of the body about the axis, denoted by But if is the moment of inertia of the body about a mean diameter, and w the angular velocity about it generated by an impluse couple M, and M' is the couple required to set the surrounding medium in motion, supposed of effective radius of gyration k', If the shot is spinning about its axis with angular velocity p, and is precessing steadily at a rate about a line parallel to the resultant momentum F at an angle 0, the velocity of the vector of angular momentum, as in the case of a top, is C i pµ sin 0- C2µ 2 sin 0 cos 0; (4) and equating this to the impressed couple (multiplied by g), that is, to gN = (c 1 -c 2)c2u 2 tan 0, (5) and dividing out sin 0, which equated to zero would imply perfect centring, we obtain C21 2 cos 0- (c 2 -c 1)c2u 2 sec 0 =o. (6) The least admissible value of p is that which makes the roots equal of this quadratic in µ, and then ICI s ec 0, ,u= z - p (7) the roots would be imaginary for a value of p smaller than given by Cip 2 - 4(c 2 -c i)c2C 2 u 2 =o, (8) p2 = 4(c 2 -c l)cl C2.

(9) c 2 Ci If the shot is moving as if fired from a gun of calibre d inches, in which the rifling makes one turn in a pitch of n calibres or nd inches, so that the angle S of the rifling is given by tan S = ird/nd = 2 d p/u, (10) '°If a denotes the density of the metal, and if the shell has a cavity homothetic with the external ellipsoidal shape, a fraction f of the linear scale; then the volume of a round shot being sird 3 , and sird 3 x of a shot x calibres long W =*ird 3 x(I -f 3)v, (20) 2 Wki 2= 61rd 3 xo(I-f 5)Q, (21) Wk22=67rd3x 2 2+0 2(I - f5)Q. If p denotes the density of the air or medium W' = sird 3 xp, (23) W' I p __ W I -1 3 k12 I k22 x2 ±i a 2= 101-1 3 '111 2= 2 tan g S = Q (l - a) x 2+ I (26) in which a/p may be replaced by 800 times the S.G. of the metal, taking water as 800 times denser than air on the average, in round numbers, and formula (to) may be written n tan 6=ir, or n6=180, when 6 is a small angle, and given in degrees.

From this formula (26) the table following has been calculated by A. G. Hadcock, and the results are in agreement with practical experience.

Table of Rifling for Stability of an Elongated Projectile, x Calibres long, giving S the Angle of Rifling, and n the Pitch of Rifling in Calibres. which is the ratio of the linear velocity of rotation 2dp to u, the velocity of advance, -T2 d2 C 22 tans = n 2 = 4 = (c 2 - Ct) cg C12 2 W! I +W a W a) ,' (k) 4 (I I) I+ w- R For a shot in air the ratio W'/W is so small that the square may be neglected, and formula (II) can be replaced for practical purpose in artillery by tan26= n2 = W i (0 - a) ( k ð)7()4 , (12) if then we can calculate /3, a, or (3-a for the external shape of the shot, this equation will give the value of 6 and n required for stability of flight in the air.

The ellipsoid is the only shape for which a and (3 have so far been determined analytically, as shown already in § 44, so we must restrict our calculation to an egg-shaped bullet, bounded by a prolate ellipsoid of revolution, in which, with b =c, Ao= fo (a2 + X)V [4(a2+X)(b +X)2]-J0 2(a2 +X)3/2(b2+X), (13) Ao+2Bo = I, (t4) _ B 0 t - A 0 I a?I-A0' Q I - Bo I-{- A o I- ? 2a

(15) The length of the shot being denoted by land the calibre by d, and the length in calibres by x l i d = 2a/2b = x, A 0 - (x2--St) 3 /2 I, 2B0= ch-tx+ x 2, (x -i)3/2 x2 -I x2Ao+2BO _ x s 2) (x 2 1) 1/ (x21)log[x-}-,1 (x 2 - I)]. 1r (52. In the steady motion the centre of the shot describes a helix, with axial velocity u cos 0+v sin 0 = (I+22tan 2 0) u cos 0?u sec 0, (I) and transverse velocity / u sin g -v cos 0= (1-? 2)a sin 0ti(a- a)u sin 0; (2) and the time of completing a turn of the spiral is 21r/µ. When µ has the critical value in (7), 2Tr _4 -C2 cos0 = 2 (x 2 +I) cos°, (3) -p Ct  which makes the circumference of the cylinder on which the helix is wrapped (u sin 0-v cos0)= 2 p (/ 3- a)(x2+I) sin 2 0 COs  0 = nd (s - a) (x 2 -1-1) sin 0 cos 0, (4) and the length of one turn of the helix -(u cos 0+v sin 0) = nd(x 2 + I); (5) thus for'x=3, the length is to times the pitch of the rifling.

53. The Motion of a Perforated Solid Liquid. the precedin investigation, the liquid stops dead when the body is brought to rest and when the body is in motion the surrounding liquid moves in uniform manner with respect to axes fixed in the body, and the force experienced by the body from the pressure of the liquid on it surface is the opposite of that required to change the motion of the liquid; this has been expressed by the dynamical equations give above. But if the body is perforated, the liquid can circulate throug a hole, in reentrant stream lines linked with the body, even whil: the body is at rest; and no reaction from the surface can influenc: this circulation, which may be supposed started in the ideal manne described in § 29, by the application of impulsive pressure across a ideal membrane closing the hole, by means of ideal mechanis connected with the body. The body is held fixed, and the reactio of the mechanism and the resultant of the impulsive pressure on th surface are a measure of the impulse, linear, ,, and angula A, µ, v, required to start the circulation.

00 ab2dX °0 ab2dX (22) L (I JS)' (18) (t9) This impulse will remain of constant magnitude, and fixed relatively to the body, which thus experiences an additional reaction from the circulation which is the opposite of the force required to change the position in space of the circulation impulse; and these extra forces must be taken into account in the dynamical equations.

An article may be consulted in the Phil. Mag., April 1893, by G. H. Bryan, in which the analytical equations of motion are deduced of a perforated solid in liquid, from considerations purely hydrodynamical.

The effect of an external circulation of vortex motion on the motion of a cylinder has been investigated in § 29; a similar procedure will show the influence of circulation through a hole in a solid, taking as the simplest illustration a ring-shaped figure, with uniplanar motion, and denoting by the resultant axial linear momentum of the circulation.

As the ring is moved from 0 to 0' in time t, with velocity Q, and angular velocity R, the components of liquid momentum change from aM'U +E and SM'V along Ox and Oy to aM'U'+ and /3M'V' along O'x' and O'y', (I) the axis of the ring changing from Ox to O'x'; and U = Q cos 0, V = Q sin 0, U' =Q cos (o - Rt), V' =Q sin (0 - Rt), (2) so that the increase of the components of momentum, X 1, Y 1, and N1, linear and angular, are X 1 = (aM'U'+ 0 cos Rt - aM'U - - 1 3M'V' sin Rt =(a - (3)M'Q sin_(0 - Rt) sin Rt - ver Rt (3) Y 1 = (aM'U'+) sin Rt-[-13M'V' cos Rt - (3M'V = (a - (3) M'Q cos (0 - Rt) sin Rt +t sin RT, N1=[ - (aM'U'+E) sin (0 - Rt)+ 1 3M'V' cos (o - Rt)]OO' = [ - (a - 1 3) M'Q cos (o - Rt) sin (o - Rt) - sin (o - Rt) ]Qt.

The components of force, X, Y, and N, acting on the liquid at 0, and reacting on the body, are then X=It. X i /t=(a - (3)M'QR sin 0= (a - (3)M'VR, (6) Y=It. Y I /t= (a - (3)M'QR cos o+ER= (a - (3)M'UR+0R, (7) Z = lt. ZI /t = - (a - s) M'Q 2 sine cos ° - EQ sin() =[ - (a - (3)M'U+E]V(8) Now suppose the cylinder is free; the additional forces acting on the body are the components of kinetic reaction of the liquid - aM' (Ç_vR), - (3M' (-- E -FUR), - EC' dR, (9) so that its equations of motion are M (Ç - vR) _ - aM' (_vR) - (a - $) M'VR, (io) M (Ç+uR) = - OM' (dV+U R) - (a - ()M'UR - R, '(II) C dR = dR + (a - Q)M'UV+0V; (12) and putting as before M+aM'=ci, M+13M' = c2, C+EC'=C3, ci dU - c2VR=o, dV +(c1U+E)R=o, c 3 dR - (c 1 U+ - c 2 U)V =o; showing the modification of the equations of plane motion, due to the component E of the circulation.

The integral of (14) and (15) may be written ciU+E=Fcoso, c 2 V= - Fsino, dx F cost o F sinz o 71 = U cos o - V sin o = cl + c c ic os o, chi = U sine +V coso= (F - F) sin cos o - l sino, (19) c i 2 2 2 sin o cos o - l ? sin o= F dl, (20) C3 do F2 h _ F2 cos 2 o F 2 sin z o F dt y - V C G c +2 c1 coso+H]; (21) 1 z so that cos 0 and y is an elliptic function of the time.

When is absent, dx/dt is always positive, and the centre of the body cannot describe loops; but with E, the influence may be great enough to make /dt change sign, and so loops occur, as shown in A. B. Basset's Hydrodynamics, i. 192, resembling the trochoidal curves, which can be looped, investigated in § 29 for the motion of a cylinder under gravity, when surrounded by a vortex.

The branch of hydrodynamics which discusses wave motion in a liquid or gas is given now in the articles Sound and Wave; while the influence of viscosity is considered under Hydraulics.

REFERENCES. - FOr the history and references to the original memoirs see Report to the British Association, by G. G. Stokes (1846), and W. M. Hicks (1882). See also the Fortschritte der Mathematik, and A. E. H. Love, " Hydrodynamik " in the Encyklopiidie der mathematischen Wissenschaften (1901). (A. G. G.)