EXTENSEON OF FRE$YQN'S SHEAR MEASUREMENT TECHNEQUE TO EQUGH BQUE‘QDARE‘ES‘ Thesis {or ”19 Degree of M. 5. MICHIGAN STATE UNWERSITY Li- San Hwang 1962 2n U) THE This is to certify that the thesis entitled Emmanuel: of bestow a Shem: Mutant Technique to W Boundaries” presented by Li-San Huang has been accepted towards fulfillment of the requirements for __!1_~S.__degree mwneering Major professor Date February 23. 1962 0-169 —___ LIBRARY Michigan State University EXTENSION OF PRESTON'S SHEAR MEASUREMENT TECHNIQUE TO ROUGH BOUNDARIES By Li - San Hwang A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Civil Engineering 1962 ii ACKNOW LEDGEMENTS The author is indebted to Dr. E. M. Laursen, without whose guidance and encouragement, this thesis would not have been possible; he also wishes to express his gratitude to Dr. H. R. Henry for his review of the manuscript, to Dr. C. P. Wells for his review of the analytical result, and to his colleague, Mr. J. R. Adams for his assistance. The research upon which this thesis is based was part of a larger study on "Pressure and Shear Distribution on Dune-shaped Boundaries, " conducted under a grant from the National Science Foundation. iii TABLE OF CONTENTS Page NOTATION ............................................. iv I INTRODUCTION .............................. 1 II. ANALYSIS ................................... 4 III. EXPERIMENTS .............................. 10 IV. RESULTS AND DISCUSSION ................... 13 V. CONCLUSIONS ............................... l7 BIBLIOGRAPHY ....................................... . Z7 NOTATION 3o h/k 5 Inner radius of the stagnation tube 30 a/k 5 Outer radius of the stagnation tube Coefficient in Karman-Prandtl velocity equation Coefficient in Karmén-Prandtl velocity equation Resistance coefficient Height of center of stagnation tube from zero datum Sand roughness Dummy integration parameter Stagnation pressure Wall static pressure Reynolds number Inner radius of the tubing Velocity at distance y from wall Shear velocity = To p Distance from zero datum Height of bottom of stagnation tube from zero datum Gamma function Thickness of laminar sublayer Integration parameter iv U Kinematic viscosity p Density of fluid 0' Area of stagnation tube opening To Wall shear stress P—P ( O )a Analytical pressure-shear ratio, in completely rough regime o P-P ( o ) Experimental pressure-shear ratio, in completely rough regime e o ABSTRACT EXTENSION OF PRESTON'S SHEAR MEASUREMENT TECHNIQUE TO ROUGH BOUNDARIES by Li-San Hwang Preston's shear measurement technique consists of placing a pitot tube in contact with a wall and interpreting the dynamic pressure reading obtained as a measure of the local shear on the wall. The rationale of the technique is that the velocity distribution near the wall is a function only of conditions at the wall, and the dynamic pres- sure reading in the tube is determined by the velocity distribution and the tube size. Assuming the Karman-Prandtl velocity distribution, an analysis was performed for the fully rough flow regime which gives the ratio of the dynamic pressure reading to the wall shear as a function of the tube size and? of the groughness element. _ Experimental measurements provided correction factors for the particular roughness used and for the transition regime offlow. Satisfactory results were obtained indicating that Preston's technique can be used for rough boundaries. I. INTRODUCTION Detailed knowledge of the local boundary shear is needed in many problems of fluid mechanics. Since analysis can seldom provide the desired answer if the flow is turbulent, one must resort to experi- mental measurement techniques. The empirical data from measure- ments may be immediately helpful in practical problems, or may be useful in the building of theory. Direct measurement of boundary shear has been attempted, but has not been entirely successful. To obtain a direct measurement of the shear on the boundary, one must measure the shear force on a small isolated floating element of the surface.1 Maintaining a small gap around the floating element, keeping it in the correct position, and recording the small shear force acting on the element are all difficult problems. An inherent disadvantage to this type of instrument is that it is restricted to a predetermined fixed point, thus inhibiting the free exploration of the shear distribution on the boundary surfaces. A standard indirect technique in boundary layer‘studies has been the measurement of the velocity and pressure profile normal to the boundary at successive sections and solving for the shear by the momentum principle. 2 If the shear is small or varies in a manner which is not simple, the results are less satisfactory, since the shear is obtained as the small difference between large values. Ludweig developed an indirect method which relates shear force to the heat loss from a hot-spot on the boundary to the flow medium. The insulation between the hot-spot and its surroundings must be excellent or considerable error may result. This method is not adapted to the evaluation of unknown shear distributions as it is restricted to measurements at pre-chosen fixed points. Moreover, the method can only be used when there is a laminar sublayer. In 1954 J. A. Preston successfully developed a simple technique for measuring the local shear on smooth boundaries using a pitot tube in contact with the surfaces. 4 This method is based upon the assumption of an inner law relating the local shear to the velocity distribution near the wall. Using the pressure drop in a pipe to cali- brate the instrument, Preston obtained equations relating the shear to the pitot tube reading both for the case of the laminar sublayer enveloping the tube, and for the case of the tube in the turbulent boundary layer on the smooth surface. E. Y. Hsu used the velocity distribution equations, u/u"< = k y1/7 (turbulent boundary layer) and u/u* = u"< y/V (laminar sublayer), to establish analytical relationships between the dynamic pressure read- ing'and the local boundary shear. 5 Hsu's analysis agrees well with Preston's experimental results and also appears to give good results for the boundary layer on a flat plate with ambient pressure gradients. In the present investigation an experiment was performed with a pipe to check equipment and technique. The results agreed, as expected, with Preston's and Hsu's. In practice the boundaries of the flow are more likely to be rough than smooth. Thus, an extension of Preston's relatively simple technique to rough boundaries would be very useful. An analytical relationship has been developed between the dynamic pressures acting on the pitot tube in contact with the rough surface and the local boundary shear. The relation is a function of the inside diameter of the tube, the size of the roughness, and the position of the tube in relation to the zero datum for flow in the completely rough regime. For the transition and hydraulically smooth regimes one more parameter, the ratio of the roughness to the laminar sublayer is needed. The effect of this parameter has been assessed experimentally. II. ANALYSIS The concept of the inner law is that the velocity distribution near a rough boundary depends only upon the viscosity and the density of the fluid, the roughness of the boundary, and the shear stress at the wall. The Karman-Prandtl velocity distribution equation which is generally considered to be a satisfactory approximation for pipes and channels and is often used in other situations, can be expressed as follow 5 u/u,< = C log y/kS + D ................ (l) The classic systematic experiments of J. Nikuradse utilized pipes covered tightly on the inside with ordinary building sand. 6 For better adherence of the sand grains, the pipe was filled and emptied a second time with Japanese lacquer after the inside had been coated with sand grains. The values usually quoted for C and D are based on Nikuradse's data, C = 5. 75 and D = 8. 5, except in the transition regime where D is a function of u* kS/U or kS/o'. Although the adequacy of this logarithmic equation can be debated, especially near the boundary, and the indeterminancy of the zero datum presents difficulties, Eq. (1) has been adopted for this analysis. Further assumptions which have been used in the analysis are: (l) The disturbance of the flow caused by the presence of the tube on the rough boundary may be neglected. (Z) The dynamic pressure acting on surface pitot tube, and therefore the reading obtained, is the average of puZ/Z over the open area. With the above assumptions, the dynamic pressure acting on the tube, as shown in Fig. 1, can be correlated to the local boundary shear by - 2 _ l f 2 (P—Poy-n'a, — 2 p 0' L1 dO' ........... (2) 01' P - 1Do 1 u 2 = (—) d U ‘r 2 a u* _ o ZTra In the completely rough regime, i. e. , where u* kS/V > 70, Eq. (1) can be written in the form u/uak = 5. 75 log 30 y/ks .............. (3) I do = “/ a2--(y-h)Z dy ‘1 Li t | JV ‘7? Y . ”h / ,.--2° G} 6 6366””; line of zero datum sand grains Fig. l. A stagnation‘tube resting on a rough boundary Substituting Eq. (3) into Eq. (2) results in h+a P-P o .1 3O 2 2 Z Z = E [5.75 u* 10g —k—Y] ° 2 ° \/a - (y-h) dy (4) Tra 3 h-a This integral relation can be evaluated in the following way: Let y = h + a sin 9, A: 30 h/k , s B = 30 a/ks- (where A > B) Making these substitutions in Eq.‘ (4), and changing integration limits, one obtains, (" Tr P-P 2 '— . 2 ___i :(_5_7_5_)__ 2 [log(A+B sin 9] coszede T TT 0 E Z 2 2n . Z 2 : (2—2???— [log(A+B sin 9)] cos 9 d9 ...... (5) 0 Zn m B . m 2 Now let f(m) = A (1+; sm 6) cos 6 d6 ----------- (6) 0 Upon expanding the integrand, there is obtained, 217 2 B Z f(m)=Arn (1+m— sine + (111) -B—- sin 9 + (r3) . A 2 A2 3 O K. sinad + . . .)cosZGd8 bulbs The odd powers of sin 9 disappear and the even powers can be integrated through the relation between the gamma and trigonometric functions, 1 1 5-1 1_ r<3s1r<§n1 ' 2 1 1 O F(ZS+En) then 3 3 2 I‘(—)1“(-) _ m B} E}. _.__Z _2__ f ;> The gamma functions can be evaluated by the relation I‘(s+l) = sF(s), and rm =1 and 11—5) 2W. Then ml m B2 1 m 13-1A" f(m):ZTrA [2+(2): 2.4 +(4)2 4-63" I (7) 2 4 1 B ml 3-1 B Let¢(m)=—+(m)— H) . — + 2 2 4A2 4 2 4 6 A4 Then f(m) = 21TAm ¢(m) .................... (8) Differentiating Eqs. (6) and (8) with respect to m twice and putting m = 0, one obtains 2n 2 2 f" (0) =f [log (A+B sin 8)] cos 9 d9 0 — 2w log A¢(O) +4Tl' log A¢'(O) + 2n¢"(0) Substituting Eq. (9) and the values of A and B into Eq. (5) there results P-P o _ 30h Z —16.53l<:[log k ] O S 30h a 2 a4 a6 -log ks [0. 25 (E) +0.0833(B-) +0.0704(h) +. . .] a 2 a4 a 6 +[o.25(-}—1-) +0.1146(-fi) +0.0586(-£) +...]} (10) This series converges quite rapidly, especially when the ratio a to h is small. Given an inside tube diameter and a roughness size a/ks, the pressure-shear ratio (P-Po)/'ro, varies with the relative position of the tube as shown in Fig. 4. The relative position of the tube (h - a)/k.S can also be written as zO/kS + t/ks, where 20 is the distance from the zero datum to the outside bottom of the tube and t is the thickness of the tube wall. If the tube is large compared to the sand particle, 20 equals the distance from the datum to the top of the particles, or something less than 0. 5 ks. If the tube is small it may rest below the top of the particles but zo/ks would probably be greater than zero. Intuitively, one might estimate that zo/ks would fall between 0.1 and 0.4. Examination of Fig. 4 shows that the larger the tube opening, the greater the pressure-shear ratio and the less the effect of unknowable zo/ks values. Therefore, the tube should be chosen as large as possible, but within the region in which the velocity distribution is primarily determined by the boundary shear and in which the general overall boundary configuration has little influence. The above analysis assumes that the flow is in the completely rough regime. In the hydraulically smooth and transitional regimes, D is no longer constant, but a function of u>z1/ I 0'301 z - 25 ogflwon £03353» 9% HOW mHOuodm cofiuoouuoU .w ondwwh a a s. m NS \ as 2 Ha: o2 Let \9 4| 4 as .2 o # WW” W 4\ .P .7 we as. . w or? \ \ \ A 0.2 on as or o.~ card ed «5 N. xe.:\ has: n...w\oa. Z6 own—mu Hgnmiogmuohm do 0::on 3°C.. Ho accumu— .@ mark A\ ax}: oo ow 2. co cm a on ON A: m l 0 . L b . .s a. rial r Ml. at/ . of /+. h/ww /o . 0 i4 6:. Romanian ¢ v o m o N d. a 0:3 0 m v n N o A ne.:\-uo~¢: u . «\dx o.N 27 BIBLIOGRAPHY SMITH, D. W. , and WALKER, J. H., "Skin-Friction Meas- urement in Incompressible Flow, " NACA Technical Note 4231, 1958. ABARBANEL, S. S., HAKKINEN, R. J., and TRILLING, L., "Use of a Stanton Tube for Skin-Friction Measurements, " NASA Memorandum 2-17-59W, 1959. LUDWEIG, H. , “Instrument for Measuring the Wall Shearing Stress of Turbulent Boundary Layer, " NACA Technical Memoran- dum 1284, 1950. PRESTON, J. H. , "The Determination of Turbulent Skin Friction by Means of Pitot Tubes, " Journal of the Royal Aeronautical Society, Vol. 54, February, 1954. HSU, E.’ Y.‘, "'The Measurement of Local Skin Friction by Means of Surface Pitot Tubes, " David W. Taylor Model Basin Research and Development Report No. 957, August, 1955. NIKURADSE J., "Laws of Flow in Rough Pipes, ” translated as NACA Technical Memorandum 1292, 1950. WHITTAKER, E. T., and WATSON, G. N., A Course of Modern Analysis, Cambridge, New York, 1927. DIXON, W. J. , and MASSEY, F. J. , Introduction to Statistical Analysis, McGraw-Hill Book Company, New York, 1957. ROBERTSON, J. M. , "The Turbulent Velocity Distribution in Rough Pipe, " Proceedings of the Fifth Midwestern Conference on Fluid Mechanics. I It I A ' USE only ECG '"'717111111le11[1111317111 11171111“ 83 0578