EXCHANGE H!l 18 I'M' THE MEASUREMENT OF THE COEFFICIENT OF VISCOSITY BY MEANS OF THE FORCED VIBRATIONS OF A SPHERE BY GEORGE FRANCIS McEWEN A.B., Stanford University, 1908 A THESIS Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Physics at The Leland Stanford Junior University, 1911 mss OP THE NEW ERA PR'KTtNO COMPANY LAftCASTCR. PA THE MEASUREMENT OF THE COEFFICIENT OF VISCOSITY BY MEANS OF THE FORCED VIBRATIONS OF A SPHERE BY GEORGE FRANCIS McEWEN A.B., Stanford University, 1908 A THESIS Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Physics at The Leland Stanford Junior University, 1911 PRCSS OF THE NEW ERA PP'NTINO COMPANY LANCASTER. PA IQTI (Reprinted from the PHYSICAL REVIEW, Vol. XXXIII., No. 6, December, 1911.1 THE MEASUREMENT OF THE FRICTIONAL FORCE EX- ERTED ON A SPHERE BY A VISCOUS FLUID, WHEN THE CENTER OF THE SPHERE PERFORMS SMALL PERIODIC OSCILLATIONS ALONG A STRAIGHT LINE. BY GEORGE F. McEwEN. ^HE purpose of the following investigation was to devise a method by which Stokes's law of the frictional resistance of a fluid to the motion of a sphere whose center performs small periodic oscilla- tions could be applied to the measurement of the coefficient of viscosity of fluids. The two main objects were: first, to obtain as close an agree- ment as possible between the actual working conditions and those de- manded by theory; second, to devise a process of measuring the force acting on the sphere, even for fluids having a very large coefficient of viscosity. The contents of this paper fall under the following five heads: I. The effect of the internal friction of fluids on the motion of pendu- lums, from Sir G. G. Stokes's Math, and Phys. Papers, Cambridge, 1880 and 1901, Vols. I. and III. II. An account of the method adopted in the present investigation to overcome the difficulties mentioned by Stokes, and to more nearly realize in the experimental work the ideal condition assumed in the theory from which Stokes's law was deduced. III. Experimental tests of the present method. IV. Suggestions for future research. V. Summary of the paper. I. THE EFFECT OF THE INTERNAL FRICTION OF FLUIDS ON THE MOTION OF PENDULUMS. I. Observations on the Motion of Pendulums. An account of the experiments made by Bessel, Baily, Dubuat, and Sabine, and the theoretical results obtained by Poisson, Challis and Plana is given in Stokes's Math, and Phys. Papers, Vol. III., pp. 1-7. The effect of the surrounding fluid on the time of vibration of a pendu lum was computed by Poisson, Challis, Green, and Plana from the hydrodynamical theory of a frictionless fluid. A fair agreement with the observations was found in some cases, but in many cases, especially 493 GEORGE F. McEWEN. [VOL. XXXIII. where the dimensions were small, the theory failed entirely to account for the experimental results. Because of this failure, Stokes was led to apply the equations of motion 1 of a viscous fluid to pendulum problems. 2. Stokes' s Deduction of the Law of Resistance. 2 In 1850 Stokes completed the solution of the following problem: The center of a sphere performs small periodic oscillations along a straight line; the sphere having a motion of translation only; it is re- quired to determine from the ordinary hydrodynamic equations of motion of a viscous fluid, the motion of the surrounding fluid, and the force exerted on the sphere. He assumed the velocities to be so small that their squares could be neglected, that there was no slipping of the fluid along the surface of the solid in contact with it. and that the amplitude of vibration of the sphere was very small and remained constant. His result for the force is p A d " y R dy F - A dt*- B dt> where y is the displacement and t is the time, and the coefficients A and B have the following values : where Mi = the mass of the fluid displaced by the sphere, r = the radius of the sphere, a = , where T is the period of oscillation, X = ^1 , , and fji f = coefficient of viscosity -5- density. From the same assumptions and equations, the same result has been obtained by a different method. 3 For a cylinder oscillating in a direction perpendicular to its axis, Stokes deduced the following expression : where F = the force acting upon unit length of the cylinder, m 3 equals the displacement of the center of mass of m-l from the position of equilibrium. Mi equals the mass of the sphere and its suspending wire. 6 equals the angular displacement of FD from its equilibrium position which is assumed to be horizontal. Mi equals the mass of the fluid displaced by the sphere. g equals the acceleration of gravity 2 5 = s + s 3 ; h = -- s 2 . do equals the angular displacement of the heavy pendulum used to maintain the vibrations. /i0 equals yo, the linear displacement of the point A. |8 equals the change in torque due to the elasticity of one suspending strip when bent through unit angle. &\y equals the change in the force on the wire holding the sphere, due to the elasticity of the surface film of the liquid, and to the varying length of the wire immersed. The variable force exerted on the sphere by the fluid equals Denote the variable force exerted on the wire by the fluid by: < equals the angular displacement of the lower bar from the equilibrium position. d kz ~r = the torque on the lower bar due to the friction of the air and the suspending wires. ki -,- = the torque on the upper bar due to the friction of the air, the suspending wires, and the knife-edge. 497 GEORGE F. McEWEN. [VOL. XXXIII. Fg equals the horizontal force at E acting on the knife-edge. F K equals the vertical force at E acting on the knife-edge. F A equals the vertical force at A acting on the lower bar. F B equals the vertical force at B acting on the lower bar. F c equals the vertical force at C acting on the lower bar. F D equals the vertical force at D acting on the upper bar. F equals the vertical force at F acting on the upper p*. m s equals the mass of each scale-pan and contents, and wires between C and D. m equals the mass of each scale pan and contents, and wires attached at F. The meaning of the other symbols used is indicated on the diagram of the apparatus. Fundamental Dynamical Relations: Mi'g -(A+A') 2 -(B + B') -B ~ = - F B + F A + F c - m, f g t Io"^ = - Fj>Si + FS* - F E 'L cos 6 - F E L sin 6 - k^ - 2/30 + FL sin 6 + /y sin 0, M" ^ = K" I (cos 0) (^ )V (sin 0) ^ J / = F E - F - ^ - M"g, c , . - ^ . , , de 01 1 r c {?2>3 (cos0)^-(sin f r / do \ 2 d*d i 1 /(sin e I M"l [ (cos 0) ( d - ) + (sin 0) ^ - J + M"g | No. 6.] MEASUREMENTS OF FRICTIONAL FORCE. 498 Fundamental Kinemalical Relations: x and y are the coordinates of C\, x = I sin 6, y I cos 9, dx d9 dy ^0 dt =l(c s6} dt< dt = l(sm ^dt' d$ o fl dy*_ o^ dt ~ d? ' 6l The above dynamical and kinematical relations can be so combined as to give the following ordinary differential equation : 499 GEORGE F. McEWEN. [VOL. XXXIII. in which the coefficients are constants and have the following values: P = [/*" + m*S? + w 3 Si 2 + Sl ^ (M, + A + A') + S ~m, f + ^'~], Assume jo = ce~ a ' sin at = lido, and substitute in the second member of the differential equation; the result will be : c s where P " = (P 2 _ [ a s _ o2]p" _ ap') f and P ' = o(P' - Assume that the following ordinary equation is a particular integral of the differential equation : 6 = Atf~ at sin (at - 61) + B,e~ at cos (at - 2 ), substitute in the equation and solve for the quantities: Ai, 0i, jBi, and 2 - The following values of the constants will make the assumed form for 0, a solution: P 2 + P[ 2 - a 2 ] - Pi + (P\a - 2aaP) tan t ' P 2 + P[a 2 - a 2 ] - aPi + (Pia - 2ao:P) tan 2 ' tan 0i = - No. 6.] MEASUREMENTS OF FRICTION A L FORCE. 5OO To obtain the complete integral of the equation, the complementary function, Aae~ Plt/zp sin (&/ + 3 ) would have to be added to the above value of 0. A$ and 3 are arbitrary constants, depending for their values on how the system is started, and = 2P In the present method of measurement, no observations are taken until the complementary function has a negligible value, so only the first or particular integral will be retained. In practical applications, the ratio of the maximum value of to the maximum value of 0o is determined by experiment. The relation of the coefficient B to that ratio, which is required for determining B will now be derived. The following equation for 0o is assumed: 0o = r e~ at sin at. l\ Then the expression for can be reduced to the form = A 2 e~ at sin (at - 0i + 3 ) where + B i 2 , and tan 3 = - A l P"' Let t\ be the time at which has a maximum value, and k be the time at which 0o first reaches its maximum value after t\. Let k ti = A/. By taking the first derivative of and 0o, the maximum values, 0' and 0o' given below are obtained : 0' = A 2 e Dividing the first equation by the second gives - 9J This ratio is obtained experimentally by first observing the maximum value of 0, then observing the next maximum value of 0o, and dividing the first by the second. Substituting the value of A and B in the expres- sion for A gives 2 + (A") 2 a 2 ] 5 denote (aPi-2aP) 2 by zi, denote (P 2 + P[ 2 - a 2 ] - aPi) by H^ZZ - z', 1 and let i - = Wm Then A / .. . (P 2 + P[<* 2 - fl 2 ] - c and the equation for R takes the following simplified form : 2 z 8 Now by varying n, the number of complete turns of W s , R n will vary but 22, z 3 , and z\ will remain constant, and are positive, w is a con- stant, positive or negative, and depends on the position of PF 3 when n is assumed to be zero. The maximum value of R, denoted by R w is : R w = therefore R w is a function of z\ and z 3 only, and corresponds to the special case when there is resonance. The above equation involving R n can be written in the following forms: (n - w) = WF - I = -71 (n - where V = . 22 An inspection of the above equations shows that if R n and R w are both multiplied by an arbitrary constant, the roots of the equation will be 1 Assume that A/j corresponds to one turn of Ws, and denote by n the number of turns from a given position near the top of the screw. Then n and P increase as Wi is lowered. 503 GEORGE F. McEWEN. [VOL. XXXIII. unaltered. Therefore the ratios used need only be proportional to the true values of the ratios. Also, if (A) has a constant amplitude, only the readings for 6 need to be taken. From three observations, three correspondeing values of n and R n can be obtained, and the two quantities R w and w, which will be the same in all cases can be eliminated, thus giving the value of ^ V. Now, V^ = y/~V / = 18.1 n = .495 cm. RJ = 8.30 w = 1.384. n *n *n 1*\ i+A"=39S.7^i 1 2.612 6.822 2.973 7 1.5513 2.4065 2.950 11 .9375 .8789 2.965 17 .5709 .3259 2.954 2.960 1,171 = 18.05 r 2 = .792 cm. R w * = 10.12 w = -1.700. n * *n