CIVE1465: Water Engineering I

An Introduction to Engineering Fluid Mechanics

EXAMPLE SHEETS: WORKED SOLUTIONS

Prof Andrew Sleigh
School of Civil Engineering
University of Leeds

March 2026

This version was created on:
March 10, 2026

E6.1 Fluid Properties

  1. E6.1.1

    Explain why the viscosity of a liquid decreases while that of a gas increases with a temperature rise. The following is a table of measurements for a fluid at constant temperature.

    Determine the dynamic viscosity of the fluid.

    dudy(s1)\frac{du}{dy}(s^{-1}) 0.0 0.2 0.4 0.6 0.8
    τ(Nm2)\tau(Nm^{-2}) 0.0 1.0 1.9 3.1 4.0

    using Newton’s law of viscosity

    τ=μuy\tau=\mu\frac{\partial u}{\partial y}

    where μ\mu is the viscosity. So viscosity is the gradient of a graph of shear stress against the velocity gradient of the above data, or

    μ=τu/y\mu=\frac{\tau}{\partial u/\partial y}

    Plot the data as a graph

    Refer to caption
    Figure 1: τ\tau vs u/y\partial u/\partial y

    Calculate the gradient for each section of the line

    dudy(s1)\frac{du}{dy}(s^{-1}) 0.0 0.2 0.4 0.6 0.8
    τ(Nm2)\tau(Nm^{-2}) 0.0 1.0 1.9 3.1 4.0
    Gradient - 5.0 4.75 5.17 5.0

    Thus the mean gradient = viscosity = μ=4.98Ns/m2\mu=4.98Ns/m^{2}.

  2. E6.1.2

    The density of an oil is 850kg/m3850kg/m^{3}. Find its relative density and Kinematic viscosity if the dynamic viscosity is 5×103kg/ms5\times 10^{-3}kg/ms.

    ρoil=850kg/m3\rho_{oil}=850kg/m^{3}   ρwater=1000kg/m3\rho_{water}=1000kg/m^{3}

    γoil=850/1000=0.85\gamma_{oil}=850/1000=0.85

    Dynamic viscosity =μ=5×103kg/ms=\mu=5\times 10^{-3}kg/ms

    Kinematic viscosity =ν=μ/ρ=\nu=\mu/\rho

    ν=μρ=5×1031000=5×103m2s1\nu=\frac{\mu}{\rho}=\frac{5\times 10^{-3}}{1000}=5\times 10^{-3}\;m^{2}s^{-1}
  3. E6.1.3

    The velocity distribution of a viscous liquid (dynamic viscosity μ=0.9Ns/m2\mu=0.9Ns/m^{2}) flowing over a fixed plate is given by u=0.68yy2u=0.68y-y^{2} (uu is the velocity in m/sm/s and yy is the distance from the plate in mm).

    What are the shear stresses at the plate surface and at y=0.34my=0.34m?

    u=0.68yy2u=0.68y-y^{2}
    uy=0.682y\frac{\partial u}{\partial y}=0.68-2y

    At the plate surface y=0my=0m

    uy=0.68\frac{\partial u}{\partial y}=0.68

    Calculate the shear stress at the plate face

    τ=μuy=9×0.68=0.612N/m2\tau=\mu\frac{\partial u}{\partial y}=9\times 0.68=0.612N/m^{2}

    At y=0.34my=0.34m,

    uy=0.682×0.34=0.0\frac{\partial u}{\partial y}=0.68-2\times 0.34=0.0

    As the velocity gradient is zero at y=0.34my=0.34m then the shear stress must also be zero.

  4. E6.1.4

    5.6m35.6m^{3} of oil weighs 46800N46800\;N. Find its mass density, ρ\rho and relative density, γ\gamma.

    Weight 46800=mg46800=mg

    Mass m=46800/9.81=4770.6kgm=46800/9.81=4770.6kg

    Mass density ρ=Mass/volume=4770.6/5.6=852kg/m3\rho=Mass/volume=4770.6/5.6=852kg/m^{3}

    Relative density

    γ=ρρwater=8521000=0.852\gamma=\frac{\rho}{\rho_{water}}=\frac{852}{1000}=0.852
  5. E6.1.5

    From the table of fluid properties the viscosity of water is given as 0.01008poises0.01008\;poises.

    What is this value in Ns/m2Ns/m^{2} and PasPa\;s units? μ=0.01008poise\mu=0.01008\;poise

    1poise=0.1Pas=0.1Ns/m21\;poise=0.1\;Pa\;s=0.1Ns/m^{2}

    μ=0.001008Pas=0.001008Ns/m2\mu=0.001008\;Pa\;s=0.001008Ns/m^{2}

  6. E6.1.6

    In a fluid the velocity measured at a distance of 75mm from the boundary is 1.125m/s1.125m/s. The fluid has absolute viscosity 0.048Pas0.048\;Pa\;s and a relative density of 0.913. What is the velocity gradient and shear stress at the boundary assuming a linear velocity distribution?

    μ=0.048Pas\mu=0.048\;Pa\;s and γ=0.913\gamma=0.913

    uy=1.1250.075=15s1\frac{\partial u}{\partial y}=\frac{1.125}{0.075}=15\;s^{-1}
    τ=μuy=0.048×15=0.720Pas\tau=\mu\frac{\partial u}{\partial y}=0.048\times 15=0.720\;Pa\;s

E6.2 Fluids Statics

Pressure and Manometers

  1. E6.2.1

    What will be (a) the gauge pressure and (b) the absolute pressure of water at a depth 12m below the surface?
    ρwater=1000kg/m3\rho_{\text{water}}=1000kg/m^{3}, and patmospheric=101kN/m2p_{\text{atmospheric}}=101kN/m^{2}.
    [Ans: 117.72kN/m2117.72kN/m^{2}, 218.72kN/m2218.72kN/m^{2}] a)

    pgauge\displaystyle p_{gauge} =ρgh\displaystyle=\rho gh
    =1000×9.81×12\displaystyle=1000\times 9.81\times 12
    =117720N/m2\displaystyle=117720N/m^{2}
    =118kN/m2\displaystyle=118kN/m^{2}

    b)

    pabsolute\displaystyle p_{absolute} =pgauge+patmospheric\displaystyle=p_{gauge}+p_{atmospheric}
    =(117720+101000)N/m2\displaystyle=(117720+101000)N/m^{2}
    =219kN/m2\displaystyle=219kN/m^{2}
  2. E6.2.2

    At what depth below the surface of oil, relative density 0.8, will produce a pressure of 120kN/m2120kN/m^{2}? What depth of water is this equivalent to?
    [Ans: 15.3m15.3m, 12.2m12.2m]

    a)

    ρ\displaystyle\rho =γρwater\displaystyle=\gamma\rho_{water}
    =0.8×kg/m3\displaystyle=0.8\times kg/m^{3}
    p\displaystyle p =ρgh\displaystyle=\rho gh
    h\displaystyle h =pρg=120×103800×9.81=15.29m of oil\displaystyle=\frac{p}{\rho g}=\frac{120\times 10^{3}}{800\times 9.81}=15.29m% \text{ of oil}

    b)

    ρ\displaystyle\rho =1000kg/m3\displaystyle=1000kg/m^{3}
    h\displaystyle h =pρg=120×1031000×9.81=12.23m of water\displaystyle=\frac{p}{\rho g}=\frac{120\times 10^{3}}{1000\times 9.81}=12.23m% \text{ of water}
  3. E6.2.3

    What would the pressure in kN/m2kN/m^{2} be if the equivalent head is measured as 400mm400mm of (a) mercury γ=13.6\gamma=13.6 (b) water ( c) oil specific weight 7.9kN/m37.9kN/m^{3} (d) a liquid of density 520kg/m3520kg/m^{3}?
    [Ans: 53.4kN/m253.4kN/m^{2}, 3.92kN/m23.92kN/m^{2}, 3.16kN/m23.16kN/m^{2}, 2.04kN/m22.04kN/m^{2}]

    a)

    ρ\displaystyle\rho =γρwater\displaystyle=\gamma\rho_{\text{water }}
    =13.6×1000kg/m3\displaystyle=13.6\times 1000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}
    p\displaystyle p =ρgh\displaystyle=\rho gh
    =(13.6×103)×9.81×0.4=53366N/m2\displaystyle=\left(13.6\times 10^{3}\right)\times 9.81\times 0.4=53366\mathrm% {\leavevmode\nobreak\ N}/\mathrm{m}^{2}

    b)

    p\displaystyle p =ρgh\displaystyle=\rho gh
    =(103)×9.81×0.4=3924N/m2\displaystyle=\left(10^{3}\right)\times 9.81\times 0.4=3924\mathrm{\leavevmode% \nobreak\ N}/\mathrm{m}^{2}

    c)

    ω\displaystyle\omega =ρg\displaystyle=\rho g
    p\displaystyle p =ρgh\displaystyle=\rho gh
    =(7.9×103)×0.4=3160N/m2\displaystyle=\left(7.9\times 10^{3}\right)\times 0.4=3160\mathrm{\leavevmode% \nobreak\ N}/\mathrm{m}^{2}

    d)

    p\displaystyle p =ρgh\displaystyle=\rho gh
    =520×9.81×0.4=2040N/m2\displaystyle=520\times 9.81\times 0.4=2040\mathrm{\leavevmode\nobreak\ N}/% \mathrm{m}^{2}
  4. E6.2.4

    A manometer connected to a pipe indicates a negative gauge pressure of 50 mm of mercury. What is the absolute pressure in the pipe in Newtons per square metre is the atmospheric pressure is 1 bar?
    [ 93.3kN/m293.3\mathrm{kN}/\mathrm{m}^{2} ]

    patmosphere\displaystyle p_{\text{atmosphere }} =1 bar =1×105N/m2\displaystyle=1\text{ bar }=1\times 10^{5}\mathrm{\leavevmode\nobreak\ N}/% \mathrm{m}^{2}
    pabsolute\displaystyle p_{\text{absolute }} =pgauge +patmospheric\displaystyle=p_{\text{gauge }}+p_{\text{atmospheric }}
    =ρgh+patmospheric\displaystyle=\rho gh+p_{\text{atmospheric }}
    =13.6×103×9.81×0.05+105N/m2,(Pa)\displaystyle=-13.6\times 10^{3}\times 9.81\times 0.05+10^{5}\mathrm{% \leavevmode\nobreak\ N}/\mathrm{m}^{2},(\mathrm{\leavevmode\nobreak\ Pa})
    =93.33kN/m2,(kPa)\displaystyle=93.33\mathrm{kN}/\mathrm{m}^{2},(\mathrm{kPa})
  5. E6.2.5

    What height would a water barometer need to be to measure atmospheric pressure of 1 bar?
    [Ans: >10m>10m]

    patmosphere\displaystyle p_{\text{atmosphere }} 1bar=1×105N/m2\displaystyle\approx 1\mathrm{bar}=1\times 10^{5}\mathrm{\leavevmode\nobreak\ % N}/\mathrm{m}^{2}
    105\displaystyle 10^{5} =ρgh\displaystyle=\rho gh
    h\displaystyle h =1051000×9.81=10.19m of water\displaystyle=\frac{10^{5}}{1000\times 9.81}=10.19\mathrm{\leavevmode\nobreak% \ m}\text{ of water }
    h\displaystyle h =105(13.6×103)×9.81=0.75m of mercury\displaystyle=\frac{10^{5}}{\left(13.6\times 10^{3}\right)\times 9.81}=0.75% \mathrm{\leavevmode\nobreak\ m}\text{ of mercury }
  6. E6.2.6

    An inclined manometer is required to measure an air pressure of 3mm of water to an accuracy of ±3%\pm 3\%. The inclined arm is 8mm8mm in diameter, and the larger arm has a diameter of 24mm24mm. The manometric fluid has a density 740kg/m3,740kg/m^{3}, and the scale may be read to ±0.5mm\pm 0.5mm.
    What is the angle required to ensure the desired accuracy may be achieved?
    [Ans: 7.67.6^{\circ}]

    p1p2=ρmangh=ρmang(z1+z2)p_{1}-p_{2}=\rho_{man}gh=\rho_{man}g\left(z_{1}+z_{2}\right)

    Volume moved from left to right =z1A1=z2sinθA2=xA2\quad=z_{1}A_{1}=\frac{z_{2}}{\sin\theta}A_{2}=xA_{2}

    =z1\displaystyle=z_{1} πD24=z2sinθπd24=xπd24\displaystyle\frac{\pi D^{2}}{4}=\frac{z_{2}}{\sin\theta}\frac{\pi d^{2}}{4}=x% \frac{\pi d^{2}}{4}
    z1\displaystyle z_{1} =z2sinθd2D2=xd2D2\displaystyle=\frac{z_{2}}{\sin\theta}\frac{d^{2}}{D^{2}}=x\frac{d^{2}}{D^{2}}
    p1p2\displaystyle p_{1}-p_{2} =ρmangx(sinθ+d2D2)\displaystyle=\rho_{\operatorname{man}}gx\left(\sin\theta+\frac{d^{2}}{D^{2}}\right)
    ρwater gh\displaystyle\rho_{\text{water }}gh =ρmangx(sinθ+d2D2)\displaystyle=\rho_{\operatorname{man}}gx\left(\sin\theta+\frac{d^{2}}{D^{2}}\right)
    ρwater gh\displaystyle\rho_{\text{water }}gh =(0.74×ρwater )gx(sinθ+0.00820.0242)\displaystyle=\left(0.74\times\rho_{\text{water }}\right)gx\left(\sin\theta+% \frac{0.008^{2}}{0.024^{2}}\right)
    h\displaystyle h =0.74x(sinθ+0.1111)\displaystyle=0.74x(\sin\theta+0.1111)

    The head being measured is 3%3\% of 3mm=0.003×0.03=0.00009m3\mathrm{\leavevmode\nobreak\ mm}=0.003\times 0.03=0.00009\mathrm{\leavevmode% \nobreak\ m}
    This 3%3\% represents the smallest measurement possible on the manometer, 0.5mm=0.0005m0.5\mathrm{\leavevmode\nobreak\ mm}=0.0005\mathrm{\leavevmode\nobreak\ m}, giving

    0.00009\displaystyle 0.00009 =0.74×0.0005(sinθ+0.1111)\displaystyle=0.74\times 0.0005(\sin\theta+0.1111)
    sinθ\displaystyle\sin\theta =0.132\displaystyle=0.132
    θ\displaystyle\theta =7.6\displaystyle=7.6^{\circ}

    [This is not the same as the answer given on the question sheet.]

  7. E6.2.7

    Determine the resultant force due to the water acting on the 1m by 2m rectangular area AB shown in the figure 2 below.
    [43560N43560N, 2.37m2.37m from O]

    Refer to caption
    Figure 2: Tank with sloping side and gates.

    The magnitude of the resultant force on a submerged plane is:

    R\displaystyle\mathrm{R} = pressure at centroid × area of surface\displaystyle=\text{ pressure at centroid }\text{ {$\times$} }\text{ area of % surface }
    R\displaystyle R =ρgz¯A\displaystyle=\rho g\bar{z}A
    =1000×9.81×(1.22+1)×(1×2)\displaystyle=1000\times 9.81\times(1.22+1)\times(1\times 2)
    =43556N/m2\displaystyle=43556\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}

    This acts at right angles to the surface through the centre of pressure.

    Sc=IOOAx¯=2 nd moment of area about a line through O1 st moment of area about a line through OSc=\frac{I_{OO}}{A\bar{x}}=\frac{2\text{ nd moment of area about a line % through }\mathrm{O}}{1\text{ st moment of area about a line through }\mathrm{O}}

    By the parallel axis theorem (which will be given in an exam), Ioo=IGG+Ax¯2I_{oo}=I_{GG}+A\bar{x}^{2}, where IGGI_{GG} is the 2nd 2^{\text{nd }} moment of area about a line through the centroid and can be found in tables.

    Sc=IGGAx¯+x¯Sc=\frac{I_{GG}}{A\bar{x}}+\bar{x}
    [Uncaptioned image]

    As the wall is vertical, Sc=DSc=D and x¯=z¯\bar{x}=\bar{z},

    Sc\displaystyle Sc =1×2312(1×2)(1.22+1)+(1.22+1)\displaystyle=\frac{1\times 2^{3}}{12(1\times 2)(1.22+1)}+(1.22+1)
    =2.37m from O\displaystyle=2.37\mathrm{\leavevmode\nobreak\ m}\text{ from }\mathrm{O}
  8. E6.2.8

    Determine the resultant force due to the water acting on the 1.25m by 2.0m triangular area CD shown in the figure 2 above. The apex of the triangle is at C.
    [Ans: 23.8×103N23.8\times 10^{3}\;N, 2.821m2.821m from P]

Depth to centre of gravity is z¯=1.0+223cos45=1.943m\bar{z}=1.0+2\frac{2}{3}\cos 45=1.943m.

R\displaystyle R =ρgz¯A\displaystyle=\rho g\bar{z}A
=1000×9.81×1.943×(2.0×1.252.0)\displaystyle=1000\times 9.81\times 1.943\times\left(\frac{2.0\times 1.25}{2.0% }\right)
=23826N/m2\displaystyle=23826\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}

Distance from P is x¯=z¯/cos45=2.748m\bar{x}=\bar{z}/\cos 45=2.748m
Distance from P to the centre of pressure is

Sc\displaystyle Sc =IooAx¯\displaystyle=\frac{I_{oo}}{A\bar{x}}
Ioo\displaystyle I_{oo} =IGG+Ax¯2\displaystyle=I_{GG}+A\bar{x}^{2}
Sc\displaystyle Sc =IGGAx¯+x¯=1.25×2336(1.25)(2.748)+(2.748)\displaystyle=\frac{I_{GG}}{A\bar{x}}+\bar{x}=\frac{1.25\times 2^{3}}{36(1.25)% (2.748)}+(2.748)
=2.829m\displaystyle=2.829\mathrm{\leavevmode\nobreak\ m}

Forces on submerged surfaces

  1. E6.9

    Obtain an expression for the depth of the centre of pressure of a plane surface wholly submerged in a fluid and inclined at an angle to the free surface of the liquid.
    A horizontal circular pipe, 1.25m1.25m diameter, is closed by a butterfly disk which rotates about a horizontal axis through its centre. Determine the torque which would have to be applied to the disk spindle to keep the disk closed in a vertical position when there is a 3m head of fresh water above the axis.
    [Ans: 1176 Nm]

    Diagram of the forces on the disc valve, based on an imaginary water surface.
    h¯=3m\bar{h}=3m, the depth to the centroid of the disc
    h=h^{\prime}= depth to the centre of pressure (or line of action of the force)
    Calculate the force:

    F\displaystyle F =ρgh¯A\displaystyle=\rho g\bar{h}A
    =1000×9.81×3×π(1.252)2\displaystyle=1000\times 9.81\times 3\times\pi\left(\frac{1.25}{2}\right)^{2}
    =36.116kN\displaystyle=36.116\mathrm{kN}

    Calculate the line of action of the force, hh^{\prime}.

    h\displaystyle h^{\prime} =2 nd moment of area about water surface 1 st moment of area about water surface\displaystyle=\frac{2\text{ nd moment of area about water surface }}{1\text{ % st moment of area about water surface }}
    =IooAh¯\displaystyle=\frac{I_{oo}}{A\bar{h}}

    By the parallel axis theorem 2nd 2^{\text{nd }} moment of area about O (in the surface) Ioo=IGG+Ah¯2I_{oo}=I_{GG}+A\bar{h}^{2} where IGGI_{GG} is the 2nd moment of area about a line through the centroid of the disc and IGG=πr4/4I_{GG}=\pi r^{4}/4.

    h\displaystyle h^{\prime} =IGGAh¯+h¯\displaystyle=\frac{I_{GG}}{A\bar{h}}+\bar{h}
    =πr44(πr2)3+3\displaystyle=\frac{\pi r^{4}}{4\left(\pi r^{2}\right)3}+3
    =r212+3=3.0326m\displaystyle=\frac{r^{2}}{12}+3=3.0326\mathrm{\leavevmode\nobreak\ m}

    So the distance from the spindle to the line of action of the force is

    x=hh¯=3.03263=0.0326mx=h^{\prime}-\bar{h}=3.0326-3=0.0326m

    And the moment required to keep the gate shut is

     moment =Fx=36.116×0.0326=1.176kNm\text{ moment }=Fx=36.116\times 0.0326=1.176\mathrm{kN}\mathrm{\leavevmode% \nobreak\ m}
  2. E6.10

    (HARDER) A dock gate is to be reinforced with three horizontal beams. If the water acts on one side only, to a depth of 6m6m, find the positions of the beams measured from the water surface so that each will carry an equal load. Give the load per meter.
    [Ans: 58860N/m58860N/m, 2.31m2.31m, 4.22m4.22m, 5.47m5.47m]

    The resultant force per unit length of gate is the area of the pressure diagram. So the total resultant force is

    R=12ρgh2=0.5×1000×9.81×62=176580N( per m length )R=\frac{1}{2}\rho gh^{2}=0.5\times 1000\times 9.81\times 6^{2}=176580N(\text{ % per }m\text{ length })

    Alternatively the resultant force is, R=\mathrm{R}= Pressure at centroid ×\times Area , (take width of gate as 1 m to give force per m)

    R=ρgh2×(h×1)=176580N( per m length )R=\rho g\frac{h}{2}\times(\mathrm{h}\times 1)=176580N(\text{ per }m\text{ % length })

    This is the resultant force exerted by the gate on the water.
    The three beams should carry an equal load, so each beam carries the load ff, where

    f=R3=58860Nf=\frac{R}{3}=58860N

    If we take moments from the surface,

    DR=fd1+fd2+fd3\displaystyle DR=fd_{1}+fd_{2}+fd_{3}
    D(3f)=f(d1+d2+d3)\displaystyle D(3f)=f\left(d_{1}+d_{2}+d_{3}\right)
    12=d1+d2+d3\displaystyle 12=d_{1}+d_{2}+d_{3}

    Taking the first beam, we can draw a pressure diagram for this, (ignoring what is below),

    We know that the resultant force, F=12ρgH2F=\frac{1}{2}\rho gH^{2}, so H=2FρgH=\sqrt{\frac{2F}{\rho g}}

    H=2Fρg=2×588601000×9.81=3.46mH=\sqrt{\frac{2F}{\rho g}}=\sqrt{\frac{2\times 58860}{1000\times 9.81}}=3.46% \mathrm{\leavevmode\nobreak\ m}

    And the force acts at 2H/32\mathrm{H}/3, so this is the position of the 1st 1^{\text{st }} beam,

     position of 1st beam =23H=2.31m\text{ position of 1st beam }=\frac{2}{3}H=2.31\mathrm{\leavevmode\nobreak\ m}

    Taking the second beam into consideration, we can draw the following pressure diagram,

    The reaction force is equal to the sum of the forces on each beam, so as before

    H=2Fρg=2×(2×58860)1000×9.81=4.9mH=\sqrt{\frac{2F}{\rho g}}=\sqrt{\frac{2\times(2\times 58860)}{1000\times 9.81% }}=4.9\mathrm{\leavevmode\nobreak\ m}

    The reaction force acts at 2H/32\mathrm{H}/3, so H=3.27m\mathrm{H}=3.27\mathrm{\leavevmode\nobreak\ m}. Taking moments from the surface,

    (2×58860)×3.27\displaystyle(2\times 58860)\times 3.27 =58860×2.31+58860×d2\displaystyle=58860\times 2.31+58860\times d_{2}
    depth to second beam d2\displaystyle\text{ depth to second beam }d_{2} =4.22m\displaystyle=4.22\mathrm{\leavevmode\nobreak\ m}

    For the third beam, from before we have,

    12\displaystyle 12 =d1+d2+d3\displaystyle=d_{1}+d_{2}+d_{3}
    depth to third beam d3\displaystyle\text{ depth to third beam }d_{3} =122.314.22=5.47m\displaystyle=12-2.31-4.22=5.47\mathrm{\leavevmode\nobreak\ m}
  3. E6.11

    The profile of a masonry dam is an arc of a circle, the arc having a radius of 30m30m and subtending an angle of 6060^{\circ} at the centre of curvature, which lies in the water surface. Determine (a) the load on the dam in N/mN/m length, (b) the position of the line of action to this pressure.
    [Ans: 4.28×106N/m4.28\times 10^{6}N/m length at depth 19.0m19.0m]

    h=30sin60=25.98m\displaystyle h=30\sin 60=25.98\mathrm{\leavevmode\nobreak\ m}
    a=30cos60=15.0m\displaystyle a=30\cos 60=15.0\mathrm{\leavevmode\nobreak\ m}

    Calculate Fv=\mathrm{F}_{\mathrm{v}}= total weight of fluid above the curved surface (per m length)

    Fv\displaystyle F_{v} =ρg( area of sector - area of triangle )\displaystyle=\rho g(\text{ area of sector - area of triangle })
    =1000×9.81×[(π302×60360)(25.98×152)]\displaystyle=1000\times 9.81\times\left[\left(\pi 30^{2}\times\frac{60}{360}% \right)-\left(\frac{25.98\times 15}{2}\right)\right]
    =2711.375kN/m\displaystyle=2711.375\mathrm{kN}/\mathrm{m}

    Calculate Fh=\mathrm{F}_{\mathrm{h}}= force on projection of curved surface onto a vertical plane

    Fh\displaystyle F_{h} =12ρgh2\displaystyle=\frac{1}{2}\rho gh^{2}
    =0.5×1000×9.81×25.982=3310.681kN/m\displaystyle=0.5\times 1000\times 9.81\times 25.98^{2}=3310.681\mathrm{kN}/% \mathrm{m}

    The resultant,

    FR\displaystyle F_{R} =Fv2+Fh2=3310.6812+2711.3752\displaystyle=\sqrt{F_{v}^{2}+F_{h}^{2}}=\sqrt{3310.681^{2}+2711.375^{2}}
    =4279.27kN/m\displaystyle=4279.27\mathrm{kN}/\mathrm{m}

    acting at the angle

    tanθ\displaystyle\tan\theta =FvFh=0.819\displaystyle=\frac{F_{v}}{F_{h}}=0.819
    θ\displaystyle\theta =39.32\displaystyle=39.32^{\circ}

    As this force act normal to the surface, it must act through the centre of radius of the dam wall. So the depth to the point where the force acts is,

    y=30sin39.31=19my=30\sin 39.31^{\circ}=19m
  4. E6.12

    (HARDER) The arch of a bridge over a stream is in the form of a semi-circle of radius 2m2m. The bridge width is 4m4m. Due to a flood, the water level is now 1.25m1.25m above the crest of the arch. Calculate (a) the upward force on the underside of the arch, (b) the horizontal thrust on one half of the arch.
    [Ans: 263.6kN263.6kN, 176.6kN176.6kN]


    a) The upward force on the arch == weight of (imaginary) water above the arch.

    Rv\displaystyle R_{v} =ρg× volume of water\displaystyle=\rho g\times\text{ volume of water }
    volume =((1.25+2)×4π222)×4=26.867m3\displaystyle=\left((1.25+2)\times 4-\frac{\pi 2^{2}}{2}\right)\times 4=26.867% \mathrm{\leavevmode\nobreak\ m}^{3}
    Rv\displaystyle R_{v} =1000×9.81×26.867=263.568kN\displaystyle=1000\times 9.81\times 26.867=263.568\mathrm{kN}

    b)

    The horizontal force on half of the arch, is equal to the force on the projection of the curved surface onto a vertical plane.

    Fh\displaystyle F_{h} = pressure at centroid × area\displaystyle=\text{ pressure at centroid {$\times$} area }
    =ρg(1.25+1)×(2×4)\displaystyle=\rho g(1.25+1)\times(2\times 4)
    =176.58kN\displaystyle=176.58\mathrm{kN}
  5. E6.13

    The face of a dam is vertical to a depth of 7.5m7.5m below the water surface then slopes at 3030^{\circ} to the vertical. If the depth of water is 17m17m what is the resultant force per metre acting on the whole face?
    [Ans: 1563.29kN1563.29kN]


    h2=17.0m\mathrm{h}_{2}=17.0\mathrm{\leavevmode\nobreak\ m}, so h1=17.07.5=9.5.x=9.5/tan60=5.485m\mathrm{h}_{1}=17.0-7.5=9.5.\quad\mathrm{x}=9.5/\tan 60=5.485\mathrm{% \leavevmode\nobreak\ m}.
    Vertical force == weight of water above the surface,

    Fv\displaystyle F_{v} =ρg(h2×x+0.5h1×x)\displaystyle=\rho g\left(h_{2}\times x+0.5h_{1}\times x\right)
    =9810×(7.5×5.485+0.5×9.5×5.485)\displaystyle=9810\times(7.5\times 5.485+0.5\times 9.5\times 5.485)
    =659.123kN/m\displaystyle=659.123\mathrm{kN}/\mathrm{m}

    The horizontal force == force on the projection of the surface on to a vertical plane.

    Fh\displaystyle F_{h} =12ρgh2\displaystyle=\frac{1}{2}\rho gh^{2}
    =0.5×1000×9.81×172\displaystyle=0.5\times 1000\times 9.81\times 17^{2}
    =1417.545kN/m\displaystyle=1417.545\mathrm{kN}/\mathrm{m}

    The resultant force is

    FR\displaystyle F_{R} =Fv2+Fh2=659.1232+1417.5452\displaystyle=\sqrt{F_{v}^{2}+F_{h}^{2}}=\sqrt{659.123^{2}+1417.545^{2}}
    =1563.29kN/m\displaystyle=1563.29\mathrm{kN}/\mathrm{m}

    And acts at the angle

    tanθ\displaystyle\tan\theta =FvFh=0.465\displaystyle=\frac{F_{v}}{F_{h}}=0.465
    θ\displaystyle\theta =24.94\displaystyle=24.94^{\circ}
  6. E6.14

    A tank with vertical sides is square in plan with 3m3m long sides. The tank contains oil of relative density 0.9 to a depth of 2.0m2.0m which is floating on water a depth of 1.5m1.5m. Calculate the force on the walls and the height of the centre of pressure from the bottom of the tank.
    [Ans: 165.54kN165.54kN, 1.15m1.15m]

Consider one wall of the tank. Draw the pressure diagram:


density of oil ρoil =0.9ρwater =900kg/m3\rho_{\text{oil }}=0.9\rho_{\text{water }}=900\mathrm{\leavevmode\nobreak\ kg}% /\mathrm{m}^{3}.
Force per unit length, F=\mathrm{F}= area under the graph == sum of the three areas =f1+f2+f3=f_{1}+f_{2}+f_{3}

f1=(900×9.81×2)×22×3=52974N\displaystyle f_{1}=\frac{(900\times 9.81\times 2)\times 2}{2}\times 3=52974% \mathrm{\leavevmode\nobreak\ N}
f2=(900×9.81×2)×1.5×3=79461N\displaystyle f_{2}=(900\times 9.81\times 2)\times 1.5\times 3=79461\mathrm{% \leavevmode\nobreak\ N}
f3=(1000×9.81×1.5)×1.52×3=33109N\displaystyle f_{3}=\frac{(1000\times 9.81\times 1.5)\times 1.5}{2}\times 3=33% 109\mathrm{\leavevmode\nobreak\ N}
F=f1+f2+f3=165544N\displaystyle F=f_{1}+f_{2}+f_{3}=165544\mathrm{\leavevmode\nobreak\ N}

To find the position of the resultant force F , we take moments from any point. We will take moments about the surface.

DF\displaystyle DF =f1d1+f2d2+f3d3\displaystyle=f_{1}d_{1}+f_{2}d_{2}+f_{3}d_{3}
165544D\displaystyle 165544D =52974×232+79461×(2+1.52)+33109×(2+231.5)\displaystyle=52974\times\frac{2}{3}2+79461\times\left(2+\frac{1.5}{2}\right)+% 33109\times\left(2+\frac{2}{3}1.5\right)
D\displaystyle D =2.347m( from surface )\displaystyle=2.347m(\text{ from surface })
=1.153m( from base of wall )\displaystyle=1.153m(\text{ from base of wall })

E6.3 Fluid Dynamics: Bernoulli Equation

Application of the Bernoulli Equation

  1. E6.3.1

    In a vertical pipe carrying water, pressure gauges are inserted at points A and B where the pipe diameters are 0.15 m and 0.075 m respectively. The point BB is 2.5 m below AA and when the flow rate down the pipe is 0.02 cumecs, the pressure at B is 14715N/m214715\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2} greater than that at A .

    Assuming the losses in the pipe between A and B can be expressed as kv22gk\frac{v^{2}}{2g} where vv is the velocity at A , find the value of kk.
    If the gauges at A and B are replaced by tubes filled with water and connected to a U -tube containing mercury of relative density 13.6, give a sketch showing how the levels in the two limbs of the U-tube differ and calculate the value of this difference in metres.
    [k=0.319,0.0794m][k=0.319,0.0794\mathrm{\leavevmode\nobreak\ m}]

    Part i)

    dA=0.15mdB=0.075mQ=0.02m3/s\displaystyle d_{A}=0.15m\quad d_{B}=0.075m\quad Q=0.02\mathrm{\leavevmode% \nobreak\ m}^{3}/\mathrm{s}
    pBpA=14715N/m2\displaystyle p_{B}-p_{A}=14715\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}
    hf=kv22g\displaystyle h_{f}=\frac{kv^{2}}{2g}

    Taking the datum at B , the Bernoulli equation becomes:

    pAρg+uA22g+zA=pBρg+uB22g+zB+kuA22g\displaystyle\frac{p_{A}}{\rho g}+\frac{u_{A}^{2}}{2g}+z_{A}=\frac{p_{B}}{\rho g% }+\frac{u_{B}^{2}}{2g}+z_{B}+k\frac{u_{A}^{2}}{2g}
    zA=2.5zB=0\displaystyle z_{A}=2.5\quad\quad z_{B}=0

    By continuity: Q=uAAA=uBAB\mathrm{Q}=\mathrm{u}_{\mathrm{A}}\mathrm{A}_{\mathrm{A}}=\mathrm{u}_{\mathrm{% B}}\mathrm{A}_{\mathrm{B}}

    uA=0.02/(π0.0752)=1.132m/s\displaystyle u_{A}=0.02/\left(\pi 0.075^{2}\right)=1.132\mathrm{\leavevmode% \nobreak\ m}/\mathrm{s}
    uB=0.02/(π0.03752)=4.527m/s\displaystyle u_{B}=0.02/\left(\pi 0.0375^{2}\right)=4.527\mathrm{\leavevmode% \nobreak\ m}/\mathrm{s}

    giving

    pBpA1000gzA+uB2uA22g\displaystyle\frac{p_{B}-p_{A}}{1000g}-z_{A}+\frac{u_{B}^{2}-u_{A}^{2}}{2g} =kuA22g\displaystyle=-k\frac{u_{A}^{2}}{2g}
    1.52.5+1.0450.065\displaystyle 1.5-2.5+1.045-0.065 =0.065k\displaystyle=-0.065k
    k\displaystyle k =0.319\displaystyle=0.319

    Part ii)

    pxxL\displaystyle p_{xxL} =ρwgzB+pB\displaystyle=\rho_{w}gz_{B}+p_{B}
    pxxR\displaystyle p_{xxR} =ρmgRp+ρwgzAρwgRp+pA\displaystyle=\rho_{m}gR_{p}+\rho_{w}gz_{A}-\rho_{w}gR_{p}+p_{A}
    pxxL\displaystyle p_{xxL} =pxxR\displaystyle=p_{xxR}
    ρwgzB+pB\displaystyle\rho_{w}gz_{B}+p_{B} =ρmgRp+ρwgzAρwgRp+pA\displaystyle=\rho_{m}gR_{p}+\rho_{w}gz_{A}-\rho_{w}gR_{p}+p_{A}
    pBpA\displaystyle p_{B}-p_{A} =ρwg(zAzB)+gRP(ρmρw)\displaystyle=\rho_{w}g\left(z_{A}-z_{B}\right)+gR_{P}\left(\rho_{m}-\rho_{w}\right)
    14715\displaystyle 14715 =1000×9.81×2.5+9.81Rp(136001000)\displaystyle=1000\times 9.81\times 2.5+9.81R_{p}(13600-1000)
    Rp\displaystyle R_{p} =0.079m\displaystyle=-0.079m
  2. E6.3.2

    A Venturimeter with an entrance diameter of 0.3 m and a throat diameter of 0.2 m is used to measure the volume of gas flowing through a pipe. The discharge coefficient of the meter is 0.96 .
    Assuming the specific weight of the gas to be constant at 19.62N/m319.62\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{3}, calculate the volume flowing when the pressure difference between the entrance and the throat is measured as 0.06 m on a water U-tube manometer.
    [ 0.816m3/s0.816\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s} ]

    What we know from the question:

    ρgg\displaystyle\rho_{g}g =19.62N/m2\displaystyle=19.62\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}
    Cd\displaystyle C_{d} =0.96\displaystyle=0.96
    d1\displaystyle d_{1} =0.3m\displaystyle=0.3\mathrm{\leavevmode\nobreak\ m}
    d2\displaystyle d_{2} =0.2m\displaystyle=0.2\mathrm{\leavevmode\nobreak\ m}

    Calculate Q.

    u1=Q/0.0707u_{1}=Q/0.0707
    u2=Q/0.0314u_{2}=Q/0.0314

    For the manometer:

    p1+ρggz\displaystyle p_{1}+\rho_{g}gz =p2+ρgg(z2Rp)+ρwgRp\displaystyle=p_{2}+\rho_{g}g\left(z_{2}-R_{p}\right)+\rho_{w}gR_{p}
    p1p2\displaystyle p_{1}-p_{2} =19.62(z2z1)+587.423\displaystyle=19.62\left(z_{2}-z_{1}\right)+587.423 (1)

    For the Venturimeter

    p1ρgg+u122g+z1\displaystyle\frac{p_{1}}{\rho_{g}g}+\frac{u_{1}^{2}}{2g}+z_{1} =p2ρgg+u222g+z2\displaystyle=\frac{p_{2}}{\rho_{g}g}+\frac{u_{2}^{2}}{2g}+z_{2}
    p1p2\displaystyle p_{1}-p_{2} =19.62(z2z1)+0.803u22\displaystyle=19.62\left(z_{2}-z_{1}\right)+0.803u_{2}^{2} (2)

    Combining (1) and (2)

    0.803u22=587.423\displaystyle 0.803u_{2}^{2}=587.423
    u2ideal =27.047m/s\displaystyle u_{2_{\text{ideal }}}=27.047\mathrm{\leavevmode\nobreak\ m}/% \mathrm{s}
    Qideal =27.047×π(0.22)2=0.85m3/s\displaystyle Q_{\text{ideal }}=27.047\times\pi\left(\frac{0.2}{2}\right)^{2}=% 0.85\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}
    Q=CdQidea =0.96×0.85=0.816m3/s\displaystyle Q=C_{d}Q_{\text{idea }}=0.96\times 0.85=0.816\mathrm{\leavevmode% \nobreak\ m}^{3}/\mathrm{s}
  3. E6.3.3

    A Venturimeter is used for measuring flow of water along a pipe. The diameter of the Venturi throat is two fifths the diameter of the pipe. The inlet and throat are connected by water filled tubes to a mercury U-tube manometer. The velocity of flow along the pipe is found to be 2.5Hm/s2.5\sqrt{H}\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}, where HH is the manometer reading in metres of mercury. Determine the loss of head between inlet and throat of the Venturi when HH is 0.49 m . (Relative density of mercury is 13.6).
    [0.23m of water]

    For the manometer:

    p1+ρwgz1\displaystyle p_{1}+\rho_{w}gz_{1} =p2+ρwg(z2H)+ρmgH\displaystyle=p_{2}+\rho_{w}g\left(z_{2}-H\right)+\rho_{m}gH
    p1p2\displaystyle p_{1}-p_{2} =ρwgz2ρwgH+ρmgHρwgz1\displaystyle=\rho_{w}gz_{2}-\rho_{w}gH+\rho_{m}gH-\rho_{w}gz_{1} (1)

    For the Venturimeter

    p1ρwg+u122g+z1\displaystyle\frac{p_{1}}{\rho_{w}g}+\frac{u_{1}^{2}}{2g}+z_{1} =p2ρwg+u222g+z2+ Losses\displaystyle=\frac{p_{2}}{\rho_{w}g}+\frac{u_{2}^{2}}{2g}+z_{2}+\text{ Losses }
    p1p2\displaystyle p_{1}-p_{2} =ρwu222+ρwgz2ρwu122ρwgz1+Lρwg\displaystyle=\frac{\rho_{w}u_{2}^{2}}{2}+\rho_{w}gz_{2}-\frac{\rho_{w}u_{1}^{% 2}}{2}-\rho_{w}gz_{1}+L\rho_{w}g (2)

    Combining (1) and (2)

    p1ρwg+u122g+z1\displaystyle\frac{p_{1}}{\rho_{w}g}+\frac{u_{1}^{2}}{2g}+z_{1} =p2ρwg+u222g+z2+ Losses\displaystyle=\frac{p_{2}}{\rho_{w}g}+\frac{u_{2}^{2}}{2g}+z_{2}+\text{ Losses }
    Lρwg\displaystyle L\rho_{w}g =Hg(ρmρw)ρw2(u22u12)\displaystyle=Hg\left(\rho_{m}-\rho_{w}\right)-\frac{\rho_{w}}{2}\left(u_{2}^{% 2}-u_{1}^{2}\right) (3)

    but at 1 . From the question

    u1\displaystyle u_{1} =2.5H=1.75m/s\displaystyle=2.5\sqrt{H}=1.75\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    u1A1\displaystyle u_{1}A_{1} =u2A2\displaystyle=u_{2}A_{2}
    1.75×πd24\displaystyle 1.75\times\pi\frac{d^{2}}{4} =u2π(2d10)2\displaystyle=u_{2}\pi\left(\frac{2d}{10}\right)^{2}
    u2\displaystyle u_{2} =10.937m/s\displaystyle=10.937\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}

    Substitute in (3)

    Losses =L=0.49×9.81(136001000)(1000/2)(10.93721.752)9.81×1000\displaystyle=L=\frac{0.49\times 9.81(13600-1000)-(1000/2)\left(10.937^{2}-1.7% 5^{2}\right)}{9.81\times 1000}
    =0.233m\displaystyle=0.233\mathrm{\leavevmode\nobreak\ m}
  4. E6.3.4

    Water is discharging from a tank through a convergent-divergent mouthpiece. The exit from the tank is rounded so that losses there may be neglected and the minimum diameter is 0.05 m .
    If the head in the tank above the centre-line of the mouthpiece is 1.83 m . a) What is the discharge?
    b) What must be the diameter at the exit if the absolute pressure at the minimum area is to be 2.44 m of water? c) What would the discharge be if the divergent part of the mouthpiece were removed. (Assume atmospheric pressure is 10 m of water).
    [0.0752m,0.0266m3/s,0.0118m3/s]\left[0.0752\mathrm{\leavevmode\nobreak\ m},0.0266\mathrm{\leavevmode\nobreak% \ m}^{3}/\mathrm{s},0.0118\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}\right]

    From the question:

    d2\displaystyle d_{2} =0.05m\displaystyle=0.05\mathrm{\leavevmode\nobreak\ m}
    minimum pressure =p2ρg=2.44m\displaystyle=\frac{p_{2}}{\rho g}=2.44\mathrm{\leavevmode\nobreak\ m}
    p1ρg\displaystyle\frac{p_{1}}{\rho g} =10m=p3ρg\displaystyle=10\mathrm{\leavevmode\nobreak\ m}=\frac{p_{3}}{\rho g}

    Apply Bernoulli:

    p1ρg+u122g+z1=p2ρg+u222g+z2=p3ρg+u322g+z3\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}=\frac{p_{3}}{\rho g}+\frac{u_{3}^{2}}{2g}+z_{3}

    If we take the datum through the orifice:

    z1=1.83mz2=z3=0u1= negligible z_{1}=1.83m\quad z_{2}=z_{3}=0\quad u_{1}=\text{ negligible }

    Between 1 and 2

    10+1.83\displaystyle 10+1.83 =2.44+u222g\displaystyle=2.44+\frac{u_{2}^{2}}{2g}
    u2\displaystyle u_{2} =13.57m/s\displaystyle=13.57\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =u2A2=13.57×π(0.052)2=0.02665m3/s\displaystyle=u_{2}A_{2}=13.57\times\pi\left(\frac{0.05}{2}\right)^{2}=0.02665% \mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    Between 1 and 3p1=p33p_{1}=p_{3}

    1.83\displaystyle 1.83 =u322g\displaystyle=\frac{u_{3}^{2}}{2g}
    u3\displaystyle u_{3} =5.99m/s\displaystyle=5.99\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =u3A3\displaystyle=u_{3}A_{3}
    0.02665\displaystyle 0.02665 =5.99×πd324\displaystyle=5.99\times\pi\frac{d_{3}^{2}}{4}
    d3\displaystyle d_{3} =0.0752m\displaystyle=0.0752\mathrm{\leavevmode\nobreak\ m}

    If the mouth piece has been removed, p1=p2p_{1}=p_{2}

    p1ρg+z1\displaystyle\frac{p_{1}}{\rho g}+z_{1} =p2ρg+u222g\displaystyle=\frac{p_{2}}{\rho g}+\frac{u_{2}^{2}}{2g}
    u2\displaystyle u_{2} =2gz1=5.99m/s\displaystyle=\sqrt{2gz_{1}}=5.99\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =5.99π0.0524=0.0118m3/s\displaystyle=5.99\pi\frac{0.05^{2}}{4}=0.0118\mathrm{\leavevmode\nobreak\ m}^% {3}/\mathrm{s}
  5. E6.3.5

    A closed tank has an orifice 0.025 m diameter in one of its vertical sides. The tank contains oil to a depth of 0.61 m above the centre of the orifice and the pressure in the air space above the oil is maintained at 13780N/m213780\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2} above atmospheric. Determine the discharge from the orifice.
    (Coefficient of discharge of the orifice is 0.61 , relative density of oil is 0.9 ).
    [ 0.00195m3/s0.00195\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s} ]

    From the question

    σ\displaystyle\sigma =0.9=ρoρw\displaystyle=0.9=\frac{\rho_{o}}{\rho_{w}}
    ρo\displaystyle\rho_{o} =900\displaystyle=900
    Cd\displaystyle C_{d} =0.61\displaystyle=0.61

    Apply Bernoulli,

    p1ρg+u122g+z1=p2ρg+u222g+z2\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}

    Take atmospheric pressure as 0 ,

    13780ρog+0.61\displaystyle\frac{13780}{\rho_{o}g}+0.61 =u222g\displaystyle=\frac{u_{2}^{2}}{2g}
    u2\displaystyle u_{2} =6.53m/s\displaystyle=6.53\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =0.61×6.53×π(0.0252)2=0.00195m3/s\displaystyle=0.61\times 6.53\times\pi\left(\frac{0.025}{2}\right)^{2}=0.00195% \mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}
  6. E6.3.6

    The discharge coefficient of a Venturimeter was found to be constant for rates of flow exceeding a certain value. Show that for this condition the loss of head due to friction in the convergent parts of the meter can be expressed as KQ2mKQ^{2}m where KK is a constant and QQ is the rate of flow in cumecs.
    Obtain the value of KK if the inlet and throat diameter of the Venturimeter are 0.102 m and 0.05 m respectively and the discharge coefficient is 0.96 .
    [ K=1060K=1060 ]

  7. E6.3.7

    A Venturimeter is to fitted in a horizontal pipe of 0.15 m diameter to measure a flow of water which may be anything up to 240m3/240\mathrm{\leavevmode\nobreak\ m}^{3}/ hour. The pressure head at the inlet for this flow is 18 m above atmospheric and the pressure head at the throat must not be lower than 7 m below atmospheric. Between the inlet and the throat there is an estimated frictional loss of 10%10\% of the difference in pressure head between these points. Calculate the minimum allowable diameter for the throat.
    [0.063m]

    From the question:

    d1=0.15mQ=240m3/hr=0.667m3/su1=Q/A=3.77m/sp1ρg=18mp2ρg=7m\begin{array}[]{ll}d_{1}=0.15\mathrm{\leavevmode\nobreak\ m}&Q=240\mathrm{% \leavevmode\nobreak\ m}^{3}/\mathrm{hr}=0.667\mathrm{\leavevmode\nobreak\ m}^{% 3}/\mathrm{s}\\ u_{1}=Q/A=3.77\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}&\\ \frac{p_{1}}{\rho g}=18\mathrm{\leavevmode\nobreak\ m}&\frac{p_{2}}{\rho g}=-7% \mathrm{\leavevmode\nobreak\ m}\end{array}

    Friction loss, from the question:

    hf=0.1(p1p2)ρgh_{f}=0.1\frac{\left(p_{1}-p_{2}\right)}{\rho g}

    Apply Bernoulli:

    p1ρg+u122g+\displaystyle\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+ =p2ρg+u222g+hf\displaystyle=\frac{p_{2}}{\rho g}+\frac{u_{2}^{2}}{2g}+h_{f}
    p1ρgp2ρg+u122ghf\displaystyle\frac{p_{1}}{\rho g}-\frac{p_{2}}{\rho g}+\frac{u_{1}^{2}}{2g}-h_% {f} =u222g\displaystyle=\frac{u_{2}^{2}}{2g}
    253.7722g2.5\displaystyle 25-\frac{3.77^{2}}{2g}-2.5 =u222g\displaystyle=\frac{u_{2}^{2}}{2g}
    u2\displaystyle u_{2} =21.346m/s\displaystyle=21.346\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =u2A2\displaystyle=u_{2}A_{2}
    0.0667\displaystyle 0.0667 =21.346×πd224\displaystyle=21.346\times\pi\frac{d_{2}^{2}}{4}
    d2\displaystyle d_{2} =0.063m\displaystyle=0.063\mathrm{\leavevmode\nobreak\ m}
  8. E6.3.8

    A Venturimeter of throat diameter 0.076 m is fitted in a 0.152 m diameter vertical pipe in which liquid of relative density 0.8 flows downwards. Pressure gauges are fitted to the inlet and to the throat sections. The throat being 0.914 m below the inlet. Taking the coefficient of the meter as 0.97 find the discharge a) when the pressure gauges read the same b)when the inlet gauge reads 15170N/m215170\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2} higher than the throat gauge.
    [0.0192m3/s,0.034m3/s]\left[0.0192\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s},0.034\mathrm{% \leavevmode\nobreak\ m}^{3}/\mathrm{s}\right]

    From the question:

    d1=0.152mA1=0.01814md2=0.076mA2=0.00454mρ=800kg/m3Cd=0.97\begin{array}[]{ll}d_{1}=0.152m&A_{1}=0.01814m\\ d_{2}=0.076m&A_{2}=0.00454m\\ \rho=800\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}&\\ C_{d}=0.97&\end{array}

    Apply Bernoulli:

    p1ρg+u122g+z1=p2ρg+u222g+z2\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}

    a) p1=p2p_{1}=p_{2}

    u122g+z1=u222g+z2\frac{u_{1}^{2}}{2g}+z_{1}=\frac{u_{2}^{2}}{2g}+z_{2}

    By continuity:

    Q=u1A1=u2A2u2=u1A1A2=u14u122g+0.914=16u122gu1=0.914×2×9.8115=1.0934m/sQ=CdA1u1Q=0.96×0.01814×1.0934=0.019m3/s\begin{gathered}Q=u_{1}A_{1}=u_{2}A_{2}\\ u_{2}=u_{1}\frac{A_{1}}{A_{2}}=u_{1}4\\ \frac{u_{1}^{2}}{2g}+0.914=\frac{16u_{1}^{2}}{2g}\\ u_{1}=\sqrt{\frac{0.914\times 2\times 9.81}{15}}=1.0934\mathrm{\leavevmode% \nobreak\ m}/\mathrm{s}\\ Q=C_{d}A_{1}u_{1}\\ Q=0.96\times 0.01814\times 1.0934=0.019\mathrm{\leavevmode\nobreak\ m}^{3}/% \mathrm{s}\end{gathered}

    b)

    p1p2=15170\displaystyle p_{1}-p_{2}=15170
    p1p2ρg=\displaystyle\frac{p_{1}-p_{2}}{\rho g}= u22u122g0.914\displaystyle\frac{u_{2}^{2}-u_{1}^{2}}{2g}-0.914
    15170ρg=\displaystyle\frac{15170}{\rho g}= Q2(220.43255.112)2g0.914\displaystyle\frac{Q^{2}\left(220.43^{2}-55.11^{2}\right)}{2g}-0.914
    55.8577=\displaystyle 55.8577= Q2(220.43255.112)\displaystyle Q^{2}\left(220.43^{2}-55.11^{2}\right)
    Q=\displaystyle Q= 0.035m3/s\displaystyle 0.035\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}
  9. E6.3.9

    A vertical cylindrical tank 2 m diameter has, at the bottom, a 0.05 m diameter sharp edged orifice for which the discharge coefficient is 0.6 .
    a) If water enters the tank at a constant rate of 0.0095 cumecs find the depth of water above the orifice when the level in the tank becomes stable.
    b) Find the time for the level to fall from 3 m to 1 m above the orifice when the inflow is turned off.
    c) If water now runs into the tank at 0.02 cumecs, the orifice remaining open, find the rate of rise in water level when the level has reached a depth of 1.7 m above the orifice.
    [a) 3.314 m , b) 881 seconds, c) 0.252m/min0.252\mathrm{\leavevmode\nobreak\ m}/\mathrm{min} ]

    From the question: Qin=0.0095m3/s,do=0.05m,Cd=0.6Q_{in}=0.0095\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s},d_{o}=0.05m,C_{d}=% 0.6
    Apply Bernoulli from the water surface (1) to the orifice (2),

    p1ρg+u122g+z1=p2ρg+u222g+z2\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}

    p1=p2p_{1}=p_{2} and u1=0.u2=2ghu_{1}=0.u_{2}=\sqrt{2gh}.
    With the datum the bottom of the cylinder, z1=h,z2=0z_{1}=h,z_{2}=0
    We need QQ in terms of the height hh measured above the orifice.

    Qout\displaystyle Q_{\text{out }} =Cdaou2=Cdao2gh\displaystyle=C_{d}a_{o}u_{2}=C_{d}a_{o}\sqrt{2gh}
    =0.6π(0.052)22×9.81h\displaystyle=0.6\pi\left(\frac{0.05}{2}\right)^{2}\sqrt{2\times 9.81}\sqrt{h}
    =0.00522h\displaystyle=0.00522\sqrt{h} (1)

    For the level in the tank to remain constant:

    inflow = out flow\displaystyle\text{ inflow }=\text{ out flow }
    Qin =Qout\displaystyle Q_{\text{in }}=Q_{\text{out }}
    0.0095=0.00522hh=3.314m\displaystyle\begin{aligned} 0.0095&=0.00522\sqrt{h}\\ h&=3.314\mathrm{\leavevmode\nobreak\ m}\end{aligned}

    (b) Write the equation for the discharge in terms of the surface height change:

    Qδt\displaystyle Q\delta t =Aδh\displaystyle=-A\delta h
    δt\displaystyle\delta t =AQδh\displaystyle=-\frac{A}{Q}\delta h

    Integrating between h1h_{1} and h2h_{2}, to give the time to change surface level

    T\displaystyle T =h1h2AQ𝑑h\displaystyle=-\int_{h_{1}}^{h_{2}}\frac{A}{Q}dh
    =601.8h1h2h1/2𝑑h\displaystyle=-601.8\int_{h_{1}}^{h_{2}}h^{-1/2}dh
    =1203.6[h1/2]h1h2\displaystyle=-1203.6\left[h^{1/2}\right]_{h_{1}}^{h_{2}}
    =1203.6[h21/2h11/2]\displaystyle=-1203.6\left[h_{2}^{1/2}-h_{1}^{1/2}\right]

    h1=3h_{1}=3 and h2=1h_{2}=1 so

    T=881secT=881\mathrm{sec}

    c) Qin Q_{\text{in }} changed to Qin =0.02m3/sQ_{\text{in }}=0.02\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    From (1) we have Qout =0.00522hQ_{\text{out }}=0.00522\sqrt{h}. The question asks for the rate of surface rise when h=1.7mh=1.7\mathrm{\leavevmode\nobreak\ m}.
    i.e.

    Qout=0.005221.7=0.0068m3/sQ_{out}=0.00522\sqrt{1.7}=0.0068\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    The rate of increase in volume is:

    Q=Qin Qout =0.020.0068=0.0132m3/sQ=Q_{\text{in }}-Q_{\text{out }}=0.02-0.0068=0.0132\mathrm{\leavevmode\nobreak% \ m}^{3}/\mathrm{s}

    As Q=Q= Area ×\times Velocity, the rate of rise in surface is

    Q=Au\displaystyle Q=Au
    u=QA=0.0132(π224)=0.0042m/s=0.252m/min\displaystyle u=\frac{Q}{A}=\frac{0.0132}{\left(\frac{\pi 2^{2}}{4}\right)}=0.% 0042\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}=0.252\mathrm{\leavevmode% \nobreak\ m}/\mathrm{min}
  10. E6.3.10

    A horizontal boiler shell (i.e. a horizontal cylinder) 2 m diameter and 10 m long is half full of water. Find the time of emptying the shell through a short vertical pipe, diameter 0.08 m , attached to the bottom of the shell. Take the coefficient of discharge to be 0.8 .
    [1370 seconds]

    d=2md=2\mathrm{\leavevmode\nobreak\ m}

    From the question W=10m,D=10mdo=0.08mCd=0.8W=10m,D=10md_{o}=0.08mC_{d}=0.8

    Apply Bernoulli from the water surface (1) to the orifice (2),

    p1ρg+u122g+z1=p2ρg+u222g+z2\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}

    p1=p2p_{1}=p_{2} and u1=0.u2=2ghu_{1}=0.u_{2}=\sqrt{2gh}.
    With the datum the bottom of the cylinder, z1=h,z2=0z_{1}=h,z_{2}=0
    We need QQ in terms of the height hh measured above the orifice.

    Qout\displaystyle Q_{\text{out }} =Cdaou2=Cdao2gh\displaystyle=C_{d}a_{o}u_{2}=C_{d}a_{o}\sqrt{2gh}
    =0.8π(0.082)22×9.81h\displaystyle=0.8\pi\left(\frac{0.08}{2}\right)^{2}\sqrt{2\times 9.81}\sqrt{h}
    =0.0178h\displaystyle=0.0178\sqrt{h}

    Write the equation for the discharge in terms of the surface height change:

    Qδt\displaystyle Q\delta t =Aδh\displaystyle=-A\delta h
    δt\displaystyle\delta t =AQδh\displaystyle=-\frac{A}{Q}\delta h

    Integrating between h1h_{1} and h2h_{2}, to give the time to change surface level

    T=h1h2AQ𝑑hT=-\int_{h_{1}}^{h_{2}}\frac{A}{Q}dh

    But we need AA in terms of hh

    Surface area A=10LA=10L, so need LL in terms of hh

    12\displaystyle 1^{2} =a2+(L2)2\displaystyle=a^{2}+\left(\frac{L}{2}\right)^{2}
    a\displaystyle a =(1h)\displaystyle=(1-h)
    12\displaystyle 1^{2} =(1h)2+(L2)2\displaystyle=(1-h)^{2}+\left(\frac{L}{2}\right)^{2}
    L\displaystyle L =2(2hh2)\displaystyle=2\sqrt{\left(2h-h^{2}\right)}
    A\displaystyle A =202hh2\displaystyle=20\sqrt{2h-h^{2}}

    Substitute this into the integral term,

    T\displaystyle T =h1h2202hh20.1078h𝑑h\displaystyle=-\int_{h_{1}}^{h_{2}}\frac{20\sqrt{2h-h^{2}}}{0.1078\sqrt{h}}dh
    =1123.6h1h22hh2h𝑑h\displaystyle=-1123.6\int_{h_{1}}^{h_{2}}\frac{\sqrt{2h-h^{2}}}{\sqrt{h}}dh
    =1123.6h1h22hh2h𝑑h\displaystyle=-1123.6\int_{h_{1}}^{h_{2}}\sqrt{\frac{2h-h^{2}}{h}}dh
    =1123.6h1h22h𝑑h\displaystyle=-1123.6\int_{h_{1}}^{h_{2}}\sqrt{2-h}dh
    =1123.6(23)[(2h)3/2]h1h2\displaystyle=-1123.6\left(-\frac{2}{3}\right)\left[(2-h)^{3/2}\right]_{h_{1}}% ^{h_{2}}
    =749.07[2.8281]=1369.6sec\displaystyle=749.07[2.828-1]=1369.6\mathrm{sec}
  11. E6.3.11

    Two cylinders standing upright contain liquid and are connected by a submerged orifice. The diameters of the cylinders are 1.75 m and 1.0 m and of the orifice, 0.08 m . The difference in levels of the liquid is initially 1.35 m . Find how long it will take for this difference to be reduced to 0.66 m if the coefficient of discharge for the orifice is 0.605 . (Work from first principles.)
    [30.7 seconds]

    do=0.108m\mathrm{d}_{\mathrm{o}}=0.108\mathrm{\leavevmode\nobreak\ m}
    A1=π(1.752)2=2.4m2A2=π(12)2=0.785m2\displaystyle A_{1}=\pi\left(\frac{1.75}{2}\right)^{2}=2.4m^{2}\quad A_{2}=\pi% \left(\frac{1}{2}\right)^{2}=0.785m^{2}
    do=0.08m,ao=π(0.082)2=0.00503m2Cd=0.605\displaystyle d_{o}=0.08m,\quad a_{o}=\pi\left(\frac{0.08}{2}\right)^{2}=0.005% 03m^{2}\quad C_{d}=0.605

    by continuity,

    A1δh1=A2δh2=Qδt-A_{1}\delta h_{1}=A_{2}\delta h_{2}=Q\delta t (1)

    defining, h=h1h2h=h_{1}-h_{2}

    δh=δh1+δh2-\delta h=-\delta h_{1}+\delta h_{2}

    Substituting this in (1) to eliminate δh2\delta h_{2}

    A1δh1\displaystyle-A_{1}\delta h_{1} =A2(δh1δh)=A2δh1A2δh\displaystyle=A_{2}\left(\delta h_{1}-\delta h\right)=A_{2}\delta h_{1}-A_{2}\delta h
    δh1\displaystyle\delta h_{1} =A2δhA1+A2\displaystyle=\frac{A_{2}\delta h}{A_{1}+A_{2}}
    A1A2δhA1+A2\displaystyle-A_{1}\frac{A_{2}\delta h}{A_{1}+A_{2}} =Qδt\displaystyle=Q\delta t (2)

    From the Bernoulli equation we can derive this expression for discharge through the submerged orifice:

    Q=Cdao2ghQ=C_{d}a_{o}\sqrt{2gh}

    So

    A1A2δhA1+A2=Cdao2ghδt\displaystyle-A_{1}\frac{A_{2}\delta h}{A_{1}+A_{2}}=C_{d}a_{o}\sqrt{2gh}\delta t
    δt=A1A2(A1+A2)Cdao2g1hδh\displaystyle\delta t=-\frac{A_{1}A_{2}}{\left(A_{1}+A_{2}\right)C_{d}a_{o}% \sqrt{2g}}\frac{1}{\sqrt{h}}\delta h

    Integrating

    T\displaystyle T =A1A2(A1+A2)Cdao2gh1h21h𝑑h\displaystyle=-\frac{A_{1}A_{2}}{\left(A_{1}+A_{2}\right)C_{d}a_{o}\sqrt{2g}}% \int_{h_{1}}^{h_{2}}\frac{1}{\sqrt{h}}dh
    =2A1A2(A1+A2)Cdao2g(h2h1)\displaystyle=-\frac{2A_{1}A_{2}}{\left(A_{1}+A_{2}\right)C_{d}a_{o}\sqrt{2g}}% \left(\sqrt{h_{2}}-\sqrt{h_{1}}\right)
    =2×2.4×0.785(2.4+0.785)×0.605×0.005032×9.81(0.81241.1619)\displaystyle=\frac{-2\times 2.4\times 0.785}{(2.4+0.785)\times 0.605\times 0.% 00503\sqrt{2\times 9.81}}(0.8124-1.1619)
    =30.7sec\displaystyle=30.7\mathrm{sec}
  12. E6.3.12

    A rectangular reservoir with vertical walls has a plan area of 60000m260000\mathrm{\leavevmode\nobreak\ m}^{2}. Discharge from the reservoir take place over a rectangular weir. The flow characteristics of the weir is Q=0.678H3/2Q=0.678H^{3/2} cumecs where HH is the depth of water above the weir crest. The sill of the weir is 3.4 m above the bottom of the reservoir. Starting with a depth of water of 4 m in the reservoir and no inflow, what will be the depth of water after one hour?
    [3.98m][3.98\mathrm{\leavevmode\nobreak\ m}]

    From the question A=60000m2,Q=0.678h3/2A=60000m^{2},Q=0.678h^{3/2}
    Write the equation for the discharge in terms of the surface height change:

    Qδt\displaystyle Q\delta t =Aδh\displaystyle=-A\delta h
    δt\displaystyle\delta t =AQδh\displaystyle=-\frac{A}{Q}\delta h

    Integrating between h1h_{1} and h2h_{2}, to give the time to change surface level

    T\displaystyle T =h1h2AQ𝑑h\displaystyle=-\int_{h_{1}}^{h_{2}}\frac{A}{Q}dh
    =600000.678h1h21h3/2𝑑h\displaystyle=-\frac{60000}{0.678}\int_{h_{1}}^{h_{2}}\frac{1}{h^{3/2}}dh
    =2×88495.58[h1/2]h1h2\displaystyle=2\times 88495.58\left[h^{-1/2}\right]_{h_{1}}^{h_{2}}

    From the question T=3600secT=3600\mathrm{sec} and hl=0.6mh_{l}=0.6\mathrm{\leavevmode\nobreak\ m}

    3600=176991.15[h21/20.61/2]\displaystyle 3600=176991.15\left[h_{2}^{-1/2}-0.6^{-1/2}\right]
    h2=0.5815m\displaystyle h_{2}=0.5815m

    Total depth =3.4+0.58=3.98m=3.4+0.58=3.98\mathrm{\leavevmode\nobreak\ m}

    Notches and weirs

  13. E6.3.13

    Deduce an expression for the discharge of water over a right-angled sharp edged V-notch, given that the coefficient of discharge is 0.61 .
    A rectangular tank 16 m by 6 m has the same notch in one of its short vertical sides. Determine the time taken for the head, measured from the bottom of the notch, to fall from 15 cm to 7.5 cm .
    [1399 seconds]

    From your notes you can derive:

    Q=815Cdtanθ22gH5/2Q=\frac{8}{15}C_{d}\tan\frac{\theta}{2}\sqrt{2g}H^{5/2}

    For this weir the equation simplifies to

    Q=1.44H5/2Q=1.44H^{5/2}

    Write the equation for the discharge in terms of the surface height change:

    Qδt\displaystyle Q\delta t =Aδh\displaystyle=-A\delta h
    δt\displaystyle\delta t =AQδh\displaystyle=-\frac{A}{Q}\delta h

    Integrating between h1h_{1} and h2h_{2}, to give the time to change surface level

    T\displaystyle T =h1h2AQ𝑑h\displaystyle=-\int_{h_{1}}^{h_{2}}\frac{A}{Q}dh
    =16×61.44h1h21h5/2𝑑h\displaystyle=-\frac{16\times 6}{1.44}\int_{h_{1}}^{h_{2}}\frac{1}{h^{5/2}}dh
    =23×66.67[h3/2]h1h2\displaystyle=\frac{2}{3}\times 66.67\left[h^{-3/2}\right]_{h_{1}}^{h_{2}}

    h1=0.15m,h2=0.075mh_{1}=0.15m,h_{2}=0.075m

    T\displaystyle T =44.44[0.0753/20.153/2]\displaystyle=44.44\left[0.075^{-3/2}-0.15^{-3/2}\right]
    =1399sec\displaystyle=1399\mathrm{sec}
  14. E6.3.14

    Derive an expression for the discharge over a sharp crested rectangular weir. A sharp edged weir is to be constructed across a stream in which the normal flow is 200 litres/sec. If the maximum flow likely to occur in the stream is 5 times the normal flow then determine the length of weir necessary to limit the rise in water level to 38.4 cm above that for normal flow. Cd=0.61C_{d}=0.61.
    [1.24m][1.24\mathrm{\leavevmode\nobreak\ m}]

    From your notes you can derive:

    Q=23Cdb2gh3/2Q=\frac{2}{3}C_{d}b\sqrt{2g}h^{3/2}

    From the question:

    Q1=0.2m3/s,h1=xQ2=1.0m3/s,h2=x+0.384\begin{array}[]{ll}Q_{1}=0.2\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s},&h_% {1}=x\\ Q_{2}=1.0\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s},&h_{2}=x+0.384\end{array}

    where xx is the height above the weir at normal flow.
    So we have two situations:

    0.2=23Cdb2gx3/2=1.801bx3/2\displaystyle 0.2=\frac{2}{3}C_{d}b\sqrt{2g}x^{3/2}=1.801bx^{3/2} (1)
    1.0=23Cdb2g(x+0.384)3/2=1.801b(x+0.384)3/2\displaystyle 1.0=\frac{2}{3}C_{d}b\sqrt{2g}(x+0.384)^{3/2}=1.801b(x+0.384)^{3% /2} (2)

    From (1) we get an expression for bb in terms of xx

    b=0.111x3/2b=0.111x^{-3/2}

    Substituting this in (2) gives,

    1.0\displaystyle 1.0 =1.801×0.111(x+0.384x)3/2\displaystyle=1.801\times 0.111\left(\frac{x+0.384}{x}\right)^{3/2}
    52/3\displaystyle 5^{2/3} =(x+0.384x)\displaystyle=\left(\frac{x+0.384}{x}\right)
    x\displaystyle x =0.1996m\displaystyle=0.1996\mathrm{\leavevmode\nobreak\ m}

    So the weir breadth is

    b\displaystyle b =0.111(0.1996)3/2\displaystyle=0.111(0.1996)^{-3/2}
    =1.24m\displaystyle=1.24\mathrm{\leavevmode\nobreak\ m}
  15. E6.3.15

    Show that the rate of flow across a triangular notch is given by Q=CdKH5/2Q=C_{d}KH^{5/2} cumecs, where CdC_{d} is an experimental coefficient, KK depends on the angle of the notch, and HH is the height of the undisturbed water level above the bottom of the notch in metres. State the reasons for the introduction of the coefficient.
    Water from a tank having a surface area of 10m210\mathrm{\leavevmode\nobreak\ m}^{2} flows over a 9090^{\circ} notch. It is found that the time taken to lower the level from 8 cm to 7 cm above the bottom of the notch is 43.5 seconds . Determine the coefficient CdC_{d} assuming that it remains constant during his period.
    [0.635]
    The proof for Q=815Cdtanθ22gH5/2=CdKH5/2Q=\frac{8}{15}C_{d}\tan\frac{\theta}{2}\sqrt{2g}H^{5/2}=C_{d}KH^{5/2} is in the notes.
    From the question:

    A=10m2θ=90h1=0.08mh2=0.07mT=43.5secA=10\mathrm{\leavevmode\nobreak\ m}^{2}\quad\theta=90^{\circ}\quad h_{1}=0.08% \mathrm{\leavevmode\nobreak\ m}\quad h_{2}=0.07\mathrm{\leavevmode\nobreak\ m}% \quad T=43.5\mathrm{sec}

    So

    Q=2.36Cdh5/2Q=2.36C_{d}h^{5/2}

    Write the equation for the discharge in terms of the surface height change:

    Qδt\displaystyle Q\delta t =Aδh\displaystyle=-A\delta h
    δt\displaystyle\delta t =AQδh\displaystyle=-\frac{A}{Q}\delta h

    Integrating between h1h_{1} and h2h_{2}, to give the time to change surface level

    T\displaystyle T =h1h2AQ𝑑h\displaystyle=-\int_{h_{1}}^{h_{2}}\frac{A}{Q}dh
    =102.36Cdh1h21h5/2𝑑h\displaystyle=-\frac{10}{2.36C_{d}}\int_{h_{1}}^{h_{2}}\frac{1}{h^{5/2}}dh
    =23×4.23Cd[h3/2]0.070.08\displaystyle=\frac{2}{3}\times\frac{4.23}{C_{d}}\left[h^{-3/2}\right]_{0.07}^% {0.08}
    43.5\displaystyle 43.5 =2.82Cd[0.073/20.083/2]\displaystyle=\frac{2.82}{C_{d}}\left[0.07^{-3/2}-0.08^{-3/2}\right]
    Cd\displaystyle C_{d} =0.635\displaystyle=0.635
  16. E6.3.16

    A reservoir with vertical sides has a plan area of 56000m256000\mathrm{\leavevmode\nobreak\ m}^{2}. Discharge from the reservoir takes place over a rectangular weir, the flow characteristic of which is Q=1.77BH3/2m3/sQ=1.77BH^{3/2}\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}. At times of maximum rainfall, water flows into the reservoir at the rate of 9m3/s9\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}. Find a) the length of weir required to discharge this quantity if head must not exceed 0.6 m ; b) the time necessary for the head to drop from 60 cm to 30 cm if the inflow suddenly stops.
    [10.94m, 3093seconds]
    From the question:

    A=56000m2Q=1.77BH3/2Qmax=9m3/sA=56000\mathrm{\leavevmode\nobreak\ m}^{2}\quad Q=1.77\mathrm{BH}^{3/2}\quad Q% _{\max}=9\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    a) Find B for H=0.6\mathrm{H}=0.6

    9=1.77B0.63/2\displaystyle 9=1.77B0.6^{3/2}
    B=10.94m\displaystyle B=10.94m

    b) Write the equation for the discharge in terms of the surface height change:

    Qδt\displaystyle Q\delta t =Aδh\displaystyle=-A\delta h
    δt\displaystyle\delta t =AQδh\displaystyle=-\frac{A}{Q}\delta h

    Integrating between h1h_{1} and h2h_{2}, to give the time to change surface level

    T\displaystyle T =h1h2AQ𝑑h\displaystyle=-\int_{h_{1}}^{h_{2}}\frac{A}{Q}dh
    =560001.77Bh1h21h3/2𝑑h\displaystyle=-\frac{56000}{1.77B}\int_{h_{1}}^{h_{2}}\frac{1}{h^{3/2}}dh
    =2×560001.77B[h1/2]0.60.3\displaystyle=\frac{2\times 56000}{1.77B}\left[h^{-1/2}\right]_{0.6}^{0.3}
    =5784[0.31/20.61/2]\displaystyle=5784\left[0.3^{-1/2}-0.6^{-1/2}\right]
    T\displaystyle T =3093sec\displaystyle=3093\mathrm{sec}
  17. E6.3.17

    Develop a formula for the discharge over a 90V90^{\circ}\mathrm{V}-notch weir in terms of head above the bottom of the V. A channel conveys 300 litres /sec/\mathrm{sec} of water. At the outlet end there is a 90V90^{\circ}\mathrm{V}-notch weir for which the coefficient of discharge is 0.58 . At what distance above the bottom of the channel should the weir be placed in order to make the depth in the channel 1.30 m ? With the weir in this position what is the depth of water in the channel when the flow is 200 litres/sec?
    [0.755m,1.218m][0.755\mathrm{\leavevmode\nobreak\ m},1.218\mathrm{\leavevmode\nobreak\ m}]
    Derive this formula from the notes: Q=815Cdtanθ22gH5/2Q=\frac{8}{15}C_{d}\tan\frac{\theta}{2}\sqrt{2g}H^{5/2}
    From the question:

    θ=90Cd0.58Q=0.3m3/s, depth of water, Z=0.3m\theta=90^{\circ}\quad C_{d}0.58\quad Q=0.3\mathrm{\leavevmode\nobreak\ m}^{3}% /\mathrm{s},\quad\text{ depth of water, }Z=0.3\mathrm{\leavevmode\nobreak\ m}

    giving the weir equation:

    Q=1.37H5/2Q=1.37H^{5/2}

    a) As HH is the height above the bottom of the V , the depth of water =Z=D+H=Z=D+H, where DD is the height of the bottom of the V from the base of the channel. So

    Q=1.37(ZD)5/2\displaystyle Q=1.37(Z-D)^{5/2}
    0.3=1.37(1.3D)5/2\displaystyle 0.3=1.37(1.3-D)^{5/2}
    D=0.755m\displaystyle D=0.755m

    b) Find ZZ when Q=0.2m3/sQ=0.2\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    0.2=1.37(Z0.755)5/2\displaystyle 0.2=1.37(Z-0.755)^{5/2}
    Z=1.218m\displaystyle Z=1.218m
  18. E6.3.18

    Show that the quantity of water flowing across a triangular V-notch of angle 2θ2\theta is
    Q=Cd815tanθ2gH5/2Q=C_{d}\frac{8}{15}\tan\theta\sqrt{2g}H^{5/2}. Find the flow if the measured head above the bottom of the V is 38 cm , when θ=45\theta=45^{\circ} and Cd=0.6C_{d}=0.6. If the flow is wanted within an accuracy of 2%2\%, what are the limiting values of the head.
    [0.126m3/s,0.377m,0.383m]\left[0.126\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s},0.377\mathrm{% \leavevmode\nobreak\ m},0.383\mathrm{\leavevmode\nobreak\ m}\right]
    Proof of the v -notch weir equation is in the notes.
    From the question:

    H=0.38mθ=45Cd=0.6H=0.38m\theta=45^{\circ}\quad C_{d}=0.6

    The weir equation becomes:

    Q\displaystyle Q =1.417H5/2\displaystyle=1.417H^{5/2}
    =1.417(0.38)5/2\displaystyle=1.417(0.38)^{5/2}
    =0.126m3/s\displaystyle=0.126\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}
    Q+2%=0.129m3/s0.129=1.417H5/2H=0.383m\begin{gathered}Q+2\%=0.129\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}\\ 0.129=1.417H^{5/2}\\ H=0.383\mathrm{\leavevmode\nobreak\ m}\end{gathered}
    Q2%\displaystyle Q-2\% =0.124m3/s\displaystyle=0.124\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}
    0.124\displaystyle 0.124 =1.417H5/2\displaystyle=1.417\mathrm{H}^{5/2}
    H\displaystyle H =0.377m\displaystyle=0.377\mathrm{\leavevmode\nobreak\ m}

E6.4 Fluid Dynamics: Momentum Equation

  1. E6.4.1

    The figure below shows a smooth curved vane attached to a rigid foundation. The jet of water, rectangular in section, 75 mm wide and 25 mm thick, strike the vane with a velocity of 25m/s25\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}. Calculate the vertical and horizontal components of the force exerted on the vane and indicate in which direction these components act.
    [Horizontal 233.4 N acting from right to left. Vertical 1324.6 N acting downwards]

    From the question:

    a1\displaystyle a_{1} =0.075×0.025=1.875×103m2\displaystyle=0.075\times 0.025=1.875\times 10^{-3}\mathrm{\leavevmode\nobreak% \ m}^{2}
    u1\displaystyle u_{1} =25m/s\displaystyle=25\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =1.875×103×25m3/s\displaystyle=1.875\times 10^{-3}\times 25\mathrm{\leavevmode\nobreak\ m}^{3}/% \mathrm{s}
    a1\displaystyle a_{1} =a2, so u1=u2\displaystyle=a_{2},\quad\text{ so }\quad u_{1}=u_{2}

    Calculate the total force using the momentum equation:

    FTx\displaystyle F_{T_{x}} =ρQ(u2cos25u1cos45)\displaystyle=\rho Q\left(u_{2}\cos 25-u_{1}\cos 45\right)
    =1000×0.0469(25cos2525cos45)\displaystyle=1000\times 0.0469(25\cos 25-25\cos 45)
    =233.44N\displaystyle=233.44\mathrm{\leavevmode\nobreak\ N}

    Ftx=ρau1(u2cosθ2u1cosθ1)F_{t_{x}}=\rho au_{1}\left(u_{2}\cos\theta_{2}-u_{1}\cos\theta_{1}\right)

    FTy\displaystyle F_{T_{y}} =ρQ(u2sin25u1sin45)\displaystyle=\rho Q\left(u_{2}\sin 25-u_{1}\sin 45\right)
    =1000×0.0469(25sin2525sin45)\displaystyle=1000\times 0.0469(25\sin 25-25\sin 45)
    =1324.6N\displaystyle=1324.6\mathrm{\leavevmode\nobreak\ N}

    Body force and pressure force are 0 .
    So force on vane:

    Rx=Ftx=233.44N\displaystyle R_{x}=-F_{t_{x}}=-233.44\mathrm{\leavevmode\nobreak\ N}
    Ry=Fty=1324.6N\displaystyle R_{y}=-F_{t_{y}}=-1324.6\mathrm{\leavevmode\nobreak\ N}
  2. E6.4.2

    A 600 mm diameter pipeline carries water under a head of 30 m with a velocity of 3m/s3\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}. This water main is fitted with a horizontal bend which turns the axis of the pipeline through 7575^{\circ} (i.e. the internal angle at the bend is 105105^{\circ} ). Calculate the resultant force on the bend and its angle to the horizontal.
    [104.044kN,5229]\left[104.044\mathrm{kN},52^{\circ}29^{\prime}\right]

    From the question:

    a=π(0.62)2=0.283m2d=0.6mh=30m\displaystyle a=\pi\left(\frac{0.6}{2}\right)^{2}=0.283\mathrm{\leavevmode% \nobreak\ m}^{2}\quad d=0.6\mathrm{\leavevmode\nobreak\ m}\quad h=30\mathrm{% \leavevmode\nobreak\ m}
    u1=u2=3m/sQ=0.848m3/s\displaystyle u_{1}=u_{2}=3\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}\quad Q=0% .848\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    Calculate total force.

    FTx=ρQ(u2xu1x)=FRx+FPx+FBx\displaystyle F_{Tx}=\rho Q\left(u_{2x}-u_{1x}\right)=F_{Rx}+F_{Px}+F_{Bx}
    FTx=1000×0.848(3cos753)=1.886kN\displaystyle F_{Tx}=1000\times 0.848(3\cos 75-3)=-1.886\mathrm{kN}
    FTy=ρQ(u2yu1y)=FRy+FPy+FBy\displaystyle F_{Ty}=\rho Q\left(u_{2y}-u_{1y}\right)=F_{Ry}+F_{Py}+F_{By}
    FTy=1000×0.848(3sin750)=2.457kN\displaystyle F_{Ty}=1000\times 0.848(3\sin 75-0)=2.457\mathrm{kN}

    Calculate the pressure force

    p1\displaystyle p_{1} =p2=p=hρg=30×1000×9.81=294.3kN/m2\displaystyle=p_{2}=p=h\rho g=30\times 1000\times 9.81=294.3\mathrm{kN}/% \mathrm{m}^{2}
    FTx\displaystyle F_{Tx} =p1a1cosθ1p2a2cosθ2\displaystyle=p_{1}a_{1}\cos\theta_{1}-p_{2}a_{2}\cos\theta_{2}
    =294300×0.283(1cos75)\displaystyle=294300\times 0.283(1-\cos 75)
    =61.73kN\displaystyle=61.73\mathrm{kN}
    FTy\displaystyle F_{Ty} =p1a1sinθ1p2a2sinθ2\displaystyle=p_{1}a_{1}\sin\theta_{1}-p_{2}a_{2}\sin\theta_{2}
    =294300×0.283(0sin75)\displaystyle=294300\times 0.283(0-\sin 75)
    =80.376kN\displaystyle=-80.376\mathrm{kN}

    There is no body force in the x or y directions.

    FRx\displaystyle F_{Rx} =FTxFPxFBx\displaystyle=F_{Tx}-F_{Px}-F_{Bx}
    =1.88661.730=63.616kN\displaystyle=-1.886-61.73-0=-63.616\mathrm{kN}
    FRy\displaystyle F_{Ry} =FTyFPyFBy\displaystyle=F_{Ty}-F_{Py}-F_{By}
    =2.457+80.3760=82.833kN\displaystyle=2.457+80.376-0=-82.833\mathrm{kN}

    These forces act on the fluid
    The resultant force on the fluid is

    FR=FRx+FRy=104.44kN\displaystyle F_{R}=\sqrt{F_{Rx}+F_{Ry}}=104.44\mathrm{kN}
    θ=tan1(FRyFRx)=5229\displaystyle\theta=\tan^{-1}\left(\frac{F_{Ry}}{F_{Rx}}\right)=52^{\circ}29^{\prime}
  3. E6.4.3

    A horizontal jet of water 2×103mm22\times 10^{3}\mathrm{\leavevmode\nobreak\ mm}^{2} cross-section and flowing at a velocity of 15m/s15\mathrm{\leavevmode\nobreak\ m}/\mathrm{s} hits a flat plate at 6060^{\circ} to the axis (of the jet) and to the horizontal. The jet is such that there is no side spread. If the plate is stationary, calculate a) the force exerted on the plate in the direction of the jet and b) the ratio between the quantity of fluid that is deflected upwards and that downwards. (Assume that there is no friction and therefore no shear force.)
    [338N, 3:1]

    From the question a2=a3=2×103m2u=15m/s\quad a_{2}=a_{3}=2\times 10^{-3}\mathrm{\leavevmode\nobreak\ m}^{2}\quad u=15% \mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Apply Bernoulli,

    p1ρg+u122g+z1=p2ρg+u222g+z2=p3ρg+u322g+z3\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}=\frac{p_{3}}{\rho g}+\frac{u_{3}^{2}}{2g}+z_{3}

    Change in height is negligible so z1=z2=z3z_{1}=z_{2}=z_{3} and pressure is always atmospheric p1=p2=p3=0p_{1}=p_{2}=p_{3}=0. So

    u1=u2=u3=15m/su_{1}=u_{2}=u_{3}=15\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}

    By continuity Q1=Q2+Q3Q_{1}=Q_{2}+Q_{3}
    so

    u1a1=u2a2+u3a3\displaystyle u_{1}a_{1}=u_{2}a_{2}+u_{3}a_{3}
    a1=a2+a3\displaystyle a_{1}=a_{2}+a_{3}

    Put the axes normal to the plate, as we know that the resultant force is normal to the plate.

    Q1=a1u=2×103×15=0.03\displaystyle Q_{1}=a_{1}u=2\times 10^{-3}\times 15=0.03
    Q1=(a2+a3)u\displaystyle Q_{1}=\left(a_{2}+a_{3}\right)u
    Q2=a2u\displaystyle Q_{2}=a_{2}u
    Q3=(a1a2)u\displaystyle Q_{3}=\left(a_{1}-a_{2}\right)u

    Calculate total force.

    FTx=ρQ(u2xu1x)=FRx+FPx+FBx\displaystyle F_{Tx}=\rho Q\left(u_{2x}-u_{1x}\right)=F_{Rx}+F_{Px}+F_{Bx}
    FTx=1000×0.03(015sin60)=390N\displaystyle F_{Tx}=1000\times 0.03(0-15\sin 60)=390\mathrm{\leavevmode% \nobreak\ N}

    Component in direction of jet =390sin60=338N=390\sin 60=338N

    As there is no force parallel to the plate Fty =0F_{\text{ty }}=0

    FTy=ρu22a2ρu32a3ρu12a1cosθ=0\displaystyle F_{Ty}=\rho u_{2}^{2}a_{2}-\rho u_{3}^{2}a_{3}-\rho u_{1}^{2}a_{% 1}\cos\theta=0
    a2a3a1cosθ=0\displaystyle a_{2}-a_{3}-a_{1}\cos\theta=0
    a1=a2+a3\displaystyle a_{1}=a_{2}+a_{3}
    a3+a1cosθ=a1a3\displaystyle a_{3}+a_{1}\cos\theta=a_{1}-a_{3}
    4a3=a1=43a2\displaystyle 4a_{3}=a_{1}=\frac{4}{3}a_{2}
    a3=13a2\displaystyle a_{3}=\frac{1}{3}a_{2}

    Thus 3/43/4 of the jet goes up, 1/41/4 down

  4. E6.4.4

    A 75 mm diameter jet of water having a velocity of 25m/s25\mathrm{\leavevmode\nobreak\ m}/\mathrm{s} strikes a flat plate, the normal of which is inclined at 3030^{\circ} to the jet. Find the force normal to the surface of the plate.
    [2.39kN][2.39\mathrm{kN}]

    From the question, djet =0.075mu1=25m/sQ=25π(0.075/2)2=0.11m3/s\quad d_{\text{jet }}=0.075m\quad u_{1}=25m/s\quad Q=25\pi(0.075/2)^{2}=0.11% \mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}
    Force normal to plate is

    FTx=ρQ(0ulx)\displaystyle F_{Tx}=\rho Q\left(0-u_{lx}\right)
    FTx=1000×0.11(025cos30)=2.39kN\displaystyle F_{Tx}=1000\times 0.11(0-25\cos 30)=2.39\mathrm{kN}
  5. E6.4.5

    The outlet pipe from a pump is a bend of 4545^{\circ} rising in the vertical plane (i.e. and internal angle of 135135^{\circ} ). The bend is 150 mm diameter at its inlet and 300 mm diameter at its outlet. The pipe axis at the inlet is horizontal and at the outlet it is 1 m higher. By neglecting friction, calculate the force and its direction if the inlet pressure is 100kN/m2100\mathrm{kN}/\mathrm{m}^{2} and the flow of water through the pipe is 0.3m3/s0.3\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}. The volume of the pipe is 0.075m30.075\mathrm{\leavevmode\nobreak\ m}^{3}.
    [13.94kN at 674067^{\circ}40^{\prime} to the horizontal]

    1&2 Draw the control volume and the axis system

    p1=100kN/m2,Q=0.3m3/sθ=45d1=0.15md2=0.3mA1=0.177m2A2=0.0707m2\begin{array}[]{lll}\mathrm{p}_{1}=100\mathrm{kN}/\mathrm{m}^{2},&\mathrm{Q}=0% .3\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}&\theta=45^{\circ}\\ \mathrm{d}_{1}=0.15\mathrm{\leavevmode\nobreak\ m}&\mathrm{\leavevmode\nobreak% \ d}_{2}=0.3\mathrm{\leavevmode\nobreak\ m}&\\ \mathrm{\leavevmode\nobreak\ A}_{1}=0.177\mathrm{\leavevmode\nobreak\ m}^{2}&% \mathrm{\leavevmode\nobreak\ A}_{2}=0.0707\mathrm{\leavevmode\nobreak\ m}^{2}&% \end{array}

    Calculate the total force in the x direction

    FTx\displaystyle F_{T_{x}} =ρQ(u2xu1x)\displaystyle=\rho Q\left(u_{2_{x}}-u_{1_{x}}\right)
    =ρQ(u2cosθu1)\displaystyle=\rho Q\left(u_{2}\cos\theta-u_{1}\right)

    by continuity A1u1=A2u2=QA_{1}u_{1}=A_{2}u_{2}=Q, so

    u1=0.3π(0.152/4)=16.98m/s\displaystyle u_{1}=\frac{0.3}{\pi\left(0.15^{2}/4\right)}=16.98\mathrm{% \leavevmode\nobreak\ m}/\mathrm{s}
    u2=0.30.0707=4.24m/s\displaystyle u_{2}=\frac{0.3}{0.0707}=4.24\mathrm{\leavevmode\nobreak\ m}/% \mathrm{s}
    FTx\displaystyle F_{T_{x}} =1000×0.3(4.24cos4516.98)\displaystyle=1000\times 0.3(4.24\cos 45-16.98)
    =1493.68N\displaystyle=-1493.68\mathrm{\leavevmode\nobreak\ N}

    and in the y -direction

    FTy\displaystyle F_{T_{y}} =ρQ(u2yu1y)\displaystyle=\rho Q\left(u_{2_{y}}-u_{1_{y}}\right)
    =ρQ(u2sinθ0)\displaystyle=\rho Q\left(u_{2}\sin\theta-0\right)
    =1000×0.3(4.24sin45)\displaystyle=1000\times 0.3(4.24\sin 45)
    =899.44N\displaystyle=899.44\mathrm{\leavevmode\nobreak\ N}

    4 Calculate the pressure force.

    FP= pressure force at 1 - pressure force at 2\displaystyle F_{P}=\text{ pressure force at }1\text{ - pressure force at }2
    FPx=p1A1cos0p2A2cosθ=p1A1p2A2cosθ\displaystyle F_{P_{x}}=p_{1}A_{1}\cos 0-p_{2}A_{2}\cos\theta=p_{1}A_{1}-p_{2}% A_{2}\cos\theta
    FPy=p1A1sin0p2A2sinθ=p2A2sinθ\displaystyle F_{P_{y}}=p_{1}A_{1}\sin 0-p_{2}A_{2}\sin\theta=-p_{2}A_{2}\sin\theta

    We know pressure at the inlet but not at the outlet. we can use Bernoulli to calculate this unknown pressure.

    p1ρg+u122g+z1=p2ρg+u222g+z2+hf\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}+h_{f}

    where hfh_{f} is the friction loss
    In the question it says this can be ignored, hf=0h_{f}=0
    The height of the pipe at the outlet is 1 m above the inlet.
    Taking the inlet level as the datum:

    z1=0z2=1m\mathrm{z}_{1}=0\quad\mathrm{z}_{2}=1\mathrm{\leavevmode\nobreak\ m}

    So the Bernoulli equation becomes:

    1000001000×9.81+16.9822×9.81+0\displaystyle\frac{100000}{1000\times 9.81}+\frac{16.98^{2}}{2\times 9.81}+0 =p21000×9.81+4.2422×9.81+1.0\displaystyle=\frac{p_{2}}{1000\times 9.81}+\frac{4.24^{2}}{2\times 9.81}+1.0
    p2\displaystyle p_{2} =225361.4N/m2\displaystyle=225361.4\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}
    FPx\displaystyle F_{P_{x}} =100000×0.0177225361.4cos45×0.0707\displaystyle=100000\times 0.0177-225361.4\cos 45\times 0.0707
    =177011266.34=9496.37kN\displaystyle=1770-11266.34=-9496.37\mathrm{kN}
    FPy\displaystyle F_{P_{y}} =225361.4sin45×0.0707\displaystyle=-225361.4\sin 45\times 0.0707
    =11266.37\displaystyle=-11266.37

    Calculate the body force The only body force is the force due to gravity. That is the weight acting in the y direction.

    FBy\displaystyle F_{B_{y}} =ρg× volume\displaystyle=-\rho g\times\text{ volume }
    =1000×9.81×0.075\displaystyle=-1000\times 9.81\times 0.075
    =12901.56N\displaystyle=-12901.56\mathrm{\leavevmode\nobreak\ N}

    There are no body forces in the x direction,

    FBx=0F_{B_{x}}=0

    6 Calculate the resultant force

    FTx=FRx+FPx+FBx\displaystyle F_{T_{x}}=F_{R_{x}}+F_{P_{x}}+F_{B_{x}}
    FTy=FRy+FPy+FBy\displaystyle F_{T_{y}}=F_{Ry}+F_{Py}+F_{B_{y}}
    FRx\displaystyle F_{R_{x}} =FTxFPxFBx\displaystyle=F_{T_{x}}-F_{P_{x}}-F_{B_{x}}
    =4193.6+9496.37\displaystyle=-4193.6+9496.37
    =5302.7N\displaystyle=5302.7\mathrm{\leavevmode\nobreak\ N}
    FRy\displaystyle F_{R_{y}} =FTyFPyFBy\displaystyle=F_{T_{y}}-F_{P_{y}}-F_{B_{y}}
    =899.44+11266.37+735.75\displaystyle=899.44+11266.37+735.75
    =12901.56N\displaystyle=12901.56\mathrm{\leavevmode\nobreak\ N}

    And the resultant force on the fluid is given by

    FR\displaystyle F_{R} =FRx2FRy2\displaystyle=\sqrt{F_{R_{x}}^{2}-F_{R_{y}}^{2}}
    =5302.72+12901.562\displaystyle=\sqrt{5302.7^{2}+12901.56^{2}}
    =13.95kN\displaystyle=13.95\mathrm{kN}

    And the direction of application is

    ϕ=tan1(FRyFRx)=tan1(12901.565302.7)=67.66\phi=\tan^{-1}\left(\frac{F_{R_{y}}}{F_{R_{x}}}\right)=\tan^{-1}\left(\frac{12% 901.56}{5302.7}\right)=67.66^{\circ}

    The force on the bend is the same magnitude but in the opposite direction

    R=FRR=-F_{R}
  6. E6.4.6

    The force exerted by a 25 mm diameter jet against a flat plate normal to the axis of the jet is 650 N . What is the flow in m3/s\mathrm{m}^{3}/\mathrm{s} ?
    [ 0.018m3/s0.018\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s} ]

    From the question, djet =0.025mFTx=650N\quad d_{\text{jet }}=0.025m\quad F_{Tx}=650N
    Force normal to plate is

    FTx=ρQ(0ulx)\displaystyle F_{Tx}=\rho Q\left(0-u_{lx}\right)
    650=1000×Q(0u)\displaystyle 650=1000\times Q(0-u)

    Q=au=(πd2/4)uQ=au=\left(\pi d^{2}/4\right)u

    650=1000au2=1000Q2/a\displaystyle 650=-1000au^{2}=-1000Q^{2}/a
    650=1000Q2/(π0.0252/4)\displaystyle 650=-1000Q^{2}/\left(\pi 0.025^{2}/4\right)
    Q=0.018m3/s\displaystyle Q=0.018m^{3}/s
  7. E6.4.7

    A curved plate deflects a 75 mm diameter jet through an angle of 4545^{\circ}. For a velocity in the jet of 40m/s40\mathrm{\leavevmode\nobreak\ m}/\mathrm{s} to the right, compute the components of the force developed against the curved plate. (Assume no friction). [ Rx=2070N,Ry=5000N\mathrm{R}_{\mathrm{x}}=2070\mathrm{\leavevmode\nobreak\ N},\mathrm{R}_{% \mathrm{y}}=5000\mathrm{\leavevmode\nobreak\ N} down]

    From the question:

    a1\displaystyle a_{1} =π0.0752/4=4.42×103m2\displaystyle=\pi 0.075^{2}/4=4.42\times 10^{-3}\mathrm{\leavevmode\nobreak\ m% }^{2}
    u1\displaystyle u_{1} =40m/s\displaystyle=40\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    Q\displaystyle Q =4.42×103×40=0.1767m3/s\displaystyle=4.42\times 10^{-3}\times 40=0.1767\mathrm{\leavevmode\nobreak\ m% }^{3}/\mathrm{s}
    a1\displaystyle a_{1} =a2, so u1=u2\displaystyle=a_{2},\quad\text{ so }\quad u_{1}=u_{2}

    Calculate the total force using the momentum equation:

    FTx\displaystyle F_{T_{x}} =ρQ(u2cos45u1)\displaystyle=\rho Q\left(u_{2}\cos 45-u_{1}\right)
    =1000×0.1767(40cos4540)\displaystyle=1000\times 0.1767(40\cos 45-40)
    =2070.17N\displaystyle=-2070.17\mathrm{\leavevmode\nobreak\ N}
    FTy\displaystyle F_{T_{y}} =ρQ(u2sin450)\displaystyle=\rho Q\left(u_{2}\sin 45-0\right)
    =1000×0.1767(40sin45)\displaystyle=1000\times 0.1767(40\sin 45)
    =4998N\displaystyle=4998\mathrm{\leavevmode\nobreak\ N}

    Body force and pressure force are 0 .
    So force on vane:

    Rx=Ftx=2070N\displaystyle R_{x}=-F_{t_{x}}=2070\mathrm{\leavevmode\nobreak\ N}
    Ry=Fty=4998N\displaystyle R_{y}=-F_{t_{y}}=-4998\mathrm{\leavevmode\nobreak\ N}
  8. E6.4.8

    A 4545^{\circ} reducing bend, 0.6 m diameter upstream, 0.3 m diameter downstream, has water flowing through it at the rate of 0.45m3/s0.45\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s} under a pressure of 1.45 bar. Neglecting any loss is head for friction, calculate the force exerted by the water on the bend, and its direction of application.
    [R=34400N\left[\mathrm{R}=34400\mathrm{\leavevmode\nobreak\ N}\right. to the right and down, θ=14]\left.\theta=14^{\circ}\right]

    1&2 Draw the control volume and the axis system

    p1=1.45×105N/m2,Q=0.45m3/sθ=45d1=0.6md2=0.3mA1=0.283m2A2=0.0707m2\begin{array}[]{lll}\mathrm{p}_{1}=1.45\times 10^{5}\mathrm{\leavevmode% \nobreak\ N}/\mathrm{m}^{2},&\mathrm{Q}=0.45\mathrm{\leavevmode\nobreak\ m}^{3% }/\mathrm{s}&\theta=45^{\circ}\\ \mathrm{d}_{1}=0.6\mathrm{\leavevmode\nobreak\ m}&\mathrm{\leavevmode\nobreak% \ d}_{2}=0.3\mathrm{\leavevmode\nobreak\ m}&\\ \mathrm{\leavevmode\nobreak\ A}_{1}=0.283\mathrm{\leavevmode\nobreak\ m}^{2}&% \mathrm{\leavevmode\nobreak\ A}_{2}=0.0707\mathrm{\leavevmode\nobreak\ m}^{2}&% \end{array}

    Calculate the total force in the x direction

    FTx\displaystyle F_{T_{x}} =ρQ(u2xu1x)\displaystyle=\rho Q\left(u_{2_{x}}-u_{1_{x}}\right)
    =ρQ(u2cosθu1)\displaystyle=\rho Q\left(u_{2}\cos\theta-u_{1}\right)

    by continuity A1u1=A2u2=QA_{1}u_{1}=A_{2}u_{2}=Q, so

    u1=0.45π(0.62/4)=1.59m/s\displaystyle u_{1}=\frac{0.45}{\pi\left(0.6^{2}/4\right)}=1.59\mathrm{% \leavevmode\nobreak\ m}/\mathrm{s}
    u2=0.450.0707=6.365m/s\displaystyle u_{2}=\frac{0.45}{0.0707}=6.365\mathrm{\leavevmode\nobreak\ m}/% \mathrm{s}
    FTx\displaystyle F_{T_{x}} =1000×0.45(6.365cos451.59)\displaystyle=1000\times 0.45(6.365\cos 45-1.59)
    =1310N\displaystyle=1310\mathrm{\leavevmode\nobreak\ N}

    and in the y -direction

    FTy\displaystyle F_{T_{y}} =ρQ(u2yu1y)\displaystyle=\rho Q\left(u_{2y}-u_{1y}\right)
    =ρQ(u2sinθ0)\displaystyle=\rho Q\left(u_{2}\sin\theta-0\right)
    =1000×0.45(6.365sin45)\displaystyle=1000\times 0.45(6.365\sin 45)
    =1800N\displaystyle=1800\mathrm{\leavevmode\nobreak\ N}

    4 Calculate the pressure force.

    FP= pressure force at 1 - pressure force at 2\displaystyle F_{P}=\text{ pressure force at }1\text{ - pressure force at }2
    FPx=p1A1cos0p2A2cosθ=p1A1p2A2cosθ\displaystyle F_{P_{x}}=p_{1}A_{1}\cos 0-p_{2}A_{2}\cos\theta=p_{1}A_{1}-p_{2}% A_{2}\cos\theta
    FPy=p1A1sin0p2A2sinθ=p2A2sinθ\displaystyle F_{P_{y}}=p_{1}A_{1}\sin 0-p_{2}A_{2}\sin\theta=-p_{2}A_{2}\sin\theta

    We know pressure at the inlet but not at the outlet.
    we can use Bernoulli to calculate this unknown pressure.

    p1ρg+u122g+z1=p2ρg+u222g+z2+hf\frac{p_{1}}{\rho g}+\frac{u_{1}^{2}}{2g}+z_{1}=\frac{p_{2}}{\rho g}+\frac{u_{% 2}^{2}}{2g}+z_{2}+h_{f}

    where hfh_{f} is the friction loss
    In the question it says this can be ignored, hf=0h_{f}=0
    Assume the pipe to be horizontal

    z1=z2\mathrm{z}_{1}=\mathrm{z}_{2}

    So the Bernoulli equation becomes:

    1450001000×9.81+1.5922×9.81\displaystyle\frac{145000}{1000\times 9.81}+\frac{1.59^{2}}{2\times 9.81} =p21000×9.81+6.36522×9.81\displaystyle=\frac{p_{2}}{1000\times 9.81}+\frac{6.365^{2}}{2\times 9.81}
    p2\displaystyle p_{2} =126007N/m2\displaystyle=126007\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}
    FPx\displaystyle F_{P_{x}} =145000×0.283126000cos45×0.0707\displaystyle=145000\times 0.283-126000\cos 45\times 0.0707
    =410356300=34735N\displaystyle=41035-6300=34735\mathrm{\leavevmode\nobreak\ N}
    FPy\displaystyle F_{P_{y}} =126000sin45×0.0707\displaystyle=-126000\sin 45\times 0.0707
    =6300N\displaystyle=-6300\mathrm{\leavevmode\nobreak\ N}

    The only body force is the force due to gravity.
    There are no body forces in the x or y directions,

    FBx=FBy=0F_{B_{x}}=F_{B_{y}}=0

    6 Calculate the resultant force

    FTx=FRx+FPx+FBx\displaystyle F_{T_{x}}=F_{R_{x}}+F_{P_{x}}+F_{B_{x}}
    FTy=FRy+FPy+FBy\displaystyle F_{T_{y}}=F_{R_{y}}+F_{P_{y}}+F_{B_{y}}
    FRx\displaystyle F_{R_{x}} =FTxFPxFBx\displaystyle=F_{T_{x}}-F_{P_{x}}-F_{B_{x}}
    =131034735\displaystyle=1310-34735
    =33425N\displaystyle=-33425\mathrm{\leavevmode\nobreak\ N}
    FRy\displaystyle F_{R_{y}} =FTyFPyFBy\displaystyle=F_{T_{y}}-F_{P_{y}}-F_{B_{y}}
    =1800+6300\displaystyle=1800+6300
    =8100N\displaystyle=8100\mathrm{\leavevmode\nobreak\ N}

    And the resultant force on the fluid is given by

    FR\displaystyle F_{R} =FRx2FRy2\displaystyle=\sqrt{F_{R_{x}}^{2}-F_{R_{y}}^{2}}
    =334252+81002\displaystyle=\sqrt{33425^{2}+8100^{2}}
    =34392kN\displaystyle=34392\mathrm{kN}

    And the direction of application is

    ϕ=tan1(FRyFRx)=tan1(810033425)=13.62\phi=\tan^{-1}\left(\frac{F_{R_{y}}}{F_{R_{x}}}\right)=\tan^{-1}\left(\frac{81% 00}{-33425}\right)=13.62^{\circ}

    The force on the bend is the same magnitude but in the opposite direction

    R=FRR=-F_{R}

E6.5 Boundary Effects: Laminar Flow

  1. E6.5.1

    The distribution of velocity, uu, in metres/sec with radius rr in metres in a smooth bore tube of 0.025 m bore follows the law, u=2.5kr2\mathrm{u}=2.5-\mathrm{kr}^{2}. Where k is a constant. The flow is laminar and the velocity at the pipe surface is zero. The fluid has a coefficient of viscosity of 0.00027kg/ms0.00027\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}\mathrm{s}. Determine (a) the rate of flow in m3/s\mathrm{m}^{3}/\mathrm{s} (b) the shearing force between the fluid and the pipe wall per metre length of pipe.
    [6.14×104m3/s,8.49×103N]\left[6.14\times 10^{-4}\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s},8.49% \times 10^{-3}\mathrm{\leavevmode\nobreak\ N}\right]
    The velocity at distance r from the centre is given in the question:

    u=2.5kr2u=2.5-kr^{2}

    Also we know: μ=0.00027kg/ms2r=0.025m\quad\mu=0.00027\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}\quad 2\mathrm{r}=% 0.025\mathrm{\leavevmode\nobreak\ m}
    We can find kk from the boundary conditions:
    when r=0.0125,u=0.0r=0.0125,\quad u=0.0\quad (boundary of the pipe)

    0.0=2.5k0.01252\displaystyle 0.0=2.5-k0.0125^{2}
    k=16000\displaystyle k=16000
    u=2.51600r2\displaystyle u=2.5-1600r^{2}

    a)

    Following along similar lines to the derivation seen in the lecture notes, we can calculate the flow δQ\delta Q through a small annulus δr\delta r :

    δQ\displaystyle\delta Q =urAannulus\displaystyle=u_{r}A_{\text{annulus }}
    Aannulus\displaystyle A_{\text{annulus }} =π(r+δr)2πr22πrδr\displaystyle=\pi(r+\delta r)^{2}-\pi r^{2}\approx 2\pi r\delta r
    δQ\displaystyle\delta Q =(2.516000r2)2πrδr\displaystyle=\left(2.5-16000r^{2}\right)2\pi r\delta r
    Q\displaystyle Q =2π00.0125(2.5r16000r3)𝑑r\displaystyle=2\pi\int_{0}^{0.0125}\left(2.5r-16000r^{3}\right)dr
    =2π[2.5r22160004r4]00.0125\displaystyle=2\pi\left[\frac{2.5r^{2}}{2}-\frac{16000}{4}r^{4}\right]_{0}^{0.% 0125}
    =6.14m3/s\displaystyle=6.14\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    b)

    The shear force is given by F=τ×(2πr)F=\tau\times(2\pi r)
    From Newtons law of viscosity

    τ\displaystyle\tau =μdudr\displaystyle=\mu\frac{du}{dr}
    dudr\displaystyle\frac{du}{dr} =2×16000r=32000r\displaystyle=-2\times 16000\mathrm{r}=-32000\mathrm{r}
    F\displaystyle F =0.00027×32000×0.0125×(2×π×0.0125)\displaystyle=-0.00027\times 32000\times 0.0125\times(2\times\pi\times 0.0125)
    =8.48×103N\displaystyle=8.48\times 10^{-3}\mathrm{\leavevmode\nobreak\ N}
  2. E6.5.2

    A liquid whose coefficient of viscosity is m flows below the critical velocity for laminar flow in a circular pipe of diameter d and with mean velocity u . Show that the pressure loss in a length of pipe is 32um/d232\mathrm{um}/\mathrm{d}^{2}.
    Oil of viscosity 0.05kg/ms0.05\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms} flows through a pipe of diameter 0.1 m with a velocity of 0.6m/s0.6\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}. Calculate the loss of pressure in a length of 120 m .
    [ 11520N/m211520\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2} ]
    See the proof in the lecture notes for
    Consider a cylinder of fluid, length L , radius r , flowing steadily in the centre of a pipe

    The fluid is in equilibrium, shearing forces equal the pressure forces.

    τ2πrL\displaystyle\tau 2\pi rL =ΔpA=Δpπr2\displaystyle=\Delta pA=\Delta p\pi r^{2}
    τ\displaystyle\tau =ΔpLr2\displaystyle=\frac{\Delta p}{L}\frac{r}{2}

    Newtons law of viscosity τ=μdudy\tau=\mu\frac{du}{dy},
    We are measuring from the pipe centre, so τ=μdudr\tau=-\mu\frac{du}{dr}
    Giving:

    ΔpLr2=μdudr\displaystyle\frac{\Delta p}{L}\frac{r}{2}=-\mu\frac{du}{dr}
    dudr=ΔpLr2μ\displaystyle\frac{du}{dr}=-\frac{\Delta p}{L}\frac{r}{2\mu}

    In an integral form this gives an expression for velocity,

    u=ΔpL12μr𝑑ru=-\frac{\Delta p}{L}\frac{1}{2\mu}\int rdr

    The value of velocity at a point distance r from the centre

    ur=ΔpLr24μ+Cu_{r}=-\frac{\Delta p}{L}\frac{r^{2}}{4\mu}+C

    At r=0r=0, (the centre of the pipe), u=umax u=u_{\text{max }}, at r=Rr=R (the pipe wall) u=0u=0;

    C=ΔpLR24μC=\frac{\Delta p}{L}\frac{R^{2}}{4\mu}

    At a point rr from the pipe centre when the flow is laminar:

    ur=ΔpL14μ(R2r2)u_{r}=\frac{\Delta p}{L}\frac{1}{4\mu}\left(R^{2}-r^{2}\right)

    The flow in an annulus of thickness δr\delta r

    δQ\displaystyle\delta Q =urAannulus\displaystyle=u_{r}A_{\text{annulus }}
    Aannulus\displaystyle A_{\text{annulus }} =π(r+δr)2πr22πrδr\displaystyle=\pi(r+\delta r)^{2}-\pi r^{2}\approx 2\pi r\delta r
    δQ\displaystyle\delta Q =ΔpL14μ(R2r2)2πrδr\displaystyle=\frac{\Delta p}{L}\frac{1}{4\mu}\left(R^{2}-r^{2}\right)2\pi r\delta r
    Q\displaystyle Q =ΔpLπ2μ0R(R2rr3)𝑑r\displaystyle=\frac{\Delta p}{L}\frac{\pi}{2\mu}\int_{0}^{R}\left(R^{2}r-r^{3}% \right)dr
    =ΔpLπR48μ=Δpπd4L128μ\displaystyle=\frac{\Delta p}{L}\frac{\pi R^{4}}{8\mu}=\frac{\Delta p\pi d^{4}% }{L128\mu}

    So the discharge can be written

    Q=ΔpLπd4128μQ=\frac{\Delta p}{L}\frac{\pi d^{4}}{128\mu}

    To get pressure loss in terms of the velocity of the flow, use the mean velocity:

    u\displaystyle u =Q/A\displaystyle=Q/A
    u\displaystyle u =Δpd232μL\displaystyle=\frac{\Delta pd^{2}}{32\mu L}
    Δp\displaystyle\Delta p =32μLud2\displaystyle=\frac{32\mu Lu}{d^{2}}
    Δp\displaystyle\Delta p =32μud2 per unit length\displaystyle=\frac{32\mu u}{d^{2}}\quad\text{ per unit length }

    b) From the question

    μ=0.05kg/msd=0.1mu=0.6m/sL=120.0m\begin{array}[]{ll}\mu=0.05\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}&d=0.1% \mathrm{\leavevmode\nobreak\ m}\\ u=0.6\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}&L=120.0\mathrm{\leavevmode% \nobreak\ m}\end{array}
    Δp=32×0.05×120×0.60.12=11520N/m2\Delta p=\frac{32\times 0.05\times 120\times 0.6}{0.1^{2}}=11520\mathrm{% \leavevmode\nobreak\ N}/\mathrm{m}^{2}
  3. E6.5.3

    A plunger of 0.08 m diameter and length 0.13 m has four small holes of diameter 5/1600m5/1600\mathrm{\leavevmode\nobreak\ m} drilled through in the direction of its length. The plunger is a close fit inside a cylinder, containing oil, such that no oil is assumed to pass between the plunger and the cylinder. If the plunger is subjected to a vertical downward force of 45 N (including its own weight) and it is assumed that the upward flow through the four small holes is laminar, determine the speed of the fall of the plunger. The coefficient of velocity of the oil is 0.2 kg/ms\mathrm{kg}/\mathrm{ms}.
    [0.00064m/s][0.00064\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}]

    Flow through each tube given by Hagen-Poiseuille equation

    Q=ΔpLπd4128μQ=\frac{\Delta p}{L}\frac{\pi d^{4}}{128\mu}

    There are 4 of these so total flow is

    Q=4ΔpLπd4128μ=Δp4π(5/1600)40.13×128×0.2=Δp3.601×1010\displaystyle Q=4\frac{\Delta p}{L}\frac{\pi d^{4}}{128\mu}=\Delta p\frac{4\pi% (5/1600)^{4}}{0.13\times 128\times 0.2}=\Delta p3.601\times 10^{-10}
    Force = pressure × area\displaystyle\text{ Force }=\text{ pressure {$\times$} area }
    F=45=Δp(π(0.082)24π(5/16002)2)\displaystyle F=45=\Delta p\left(\pi\left(\frac{0.08}{2}\right)^{2}-4\pi\left(% \frac{5/1600}{2}\right)^{2}\right)
    Δp=9007.206N/m2\displaystyle\Delta p=9007.206\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}
    Q=3.24×106m3/s\displaystyle Q=3.24\times 10^{-6}\mathrm{\leavevmode\nobreak\ m}^{3}/\mathrm{s}

    Flow up through piston = flow displaced by moving piston

    Q=Avpiston\displaystyle\mathrm{Q}=\mathrm{Av}_{\text{piston }}
    3.24×106=π×0.042×Vpiston\displaystyle 3.24\times 10^{-6}=\pi\times 0.04^{2}\times\mathrm{V}_{\text{% piston }}
    Vpiston =0.00064m/s\mathrm{V}_{\text{piston }}=0.00064\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
  4. E6.5.4

    A vertical cylinder of 0.075 metres diameter is mounted concentrically in a drum of 0.076 metres internal diameter. Oil fills the space between them to a depth of 0.2 m . The rotque required to rotate the cylinder in the drum is 4 Nm when the speed of rotation is 7.5revs/sec7.5\mathrm{revs}/\mathrm{sec}. Assuming that the end effects are negligible, calculate the coefficient of viscosity of the oil.
    [0.638kg/ms][0.638\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}]
    From the question r1=0.076/2r2=0.075/2\quad r_{1}=0.076/2\quad r_{2}=0.075/2\quad Torque =4Nm,L=0.2m=4Nm,L=0.2m
    The velocity of the edge of the cylinder is:

    ucyl=7.5×2πr=7.5×2×π×0.0375=1.767m/s\displaystyle u_{cyl}=7.5\times 2\pi r=7.5\times 2\times\pi\times 0.0375=1.767% \mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    udrum=0.0\displaystyle u_{drum}=0.0

    Torque needed to rotate cylinder

    T\displaystyle T =τ× surface area\displaystyle=\tau\times\text{ surface area }
    4\displaystyle 4 =τ(2πr2×L)\displaystyle=\tau\left(2\pi r_{2}\times L\right)
    τ\displaystyle\tau =2263.54N/m2\displaystyle=2263.54\mathrm{\leavevmode\nobreak\ N}/\mathrm{m}^{2}

    Distance between cylinder and drum =r1r2=0.0380.0375=0.005m=r_{1}-r_{2}=0.038-0.0375=0.005\mathrm{\leavevmode\nobreak\ m}
    Using Newtons law of viscosity:

    τ\displaystyle\tau =μdudr\displaystyle=\mu\frac{du}{dr}
    dudr\displaystyle\frac{du}{dr} =1.76700.0005\displaystyle=\frac{1.767-0}{0.0005}
    τ\displaystyle\tau =2263.5=μ3534\displaystyle=2263.5=\mu 3534
    μ\displaystyle\mu =0.64kg/ms\displaystyle=0.64\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}

E6.6 Dimensional Analysis: Physical Modelling

NOT STUDIED ANY MORE ON THIS COURSE

  1. E6.6.1

    A stationary sphere in water moving at a velocity of 1.6m/s1.6\mathrm{\leavevmode\nobreak\ m}/\mathrm{s} experiences a drag of 4 N . Another sphere of twice the diameter is placed in a wind tunnel. Find the velocity of the air and the drag which will give dynamically similar conditions. The ratio of kinematic viscosities of air and water is 13, and the density of air 1.28kg/m31.28\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}.
    [10.4m/s0.865N][10.4\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}0.865\mathrm{\leavevmode% \nobreak\ N}]
    Draw up the table of values you have for each variable:

    variable water air
    u 1.6m/s1.6\mathrm{\leavevmode\nobreak\ m}/\mathrm{s} uair\mathrm{u}_{\text{air }}
    Drag 4 N Dair\mathrm{D}_{\text{air }}
    vv vv 13v13v
    ρ\rho 1000kg/m31000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3} 1.28kg/m31.28\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}
    d d 2 d

    Kinematic viscosity is dynamic viscosity over density =ν=μ/ρ=\nu=\mu/\rho.
    The Reynolds number =Re=ρudμ=udv=\operatorname{Re}=\frac{\rho ud}{\mu}=\frac{ud}{v}
    Choose the three recurring (governing) variables; u,d,ρ\mathrm{u},\mathrm{d},\rho.
    From Buckinghams π\pi theorem we have mn=53=2\mathrm{m}-\mathrm{n}=5-3=2 non-dimensional groups.

    ϕ(u,d,ρ,D,v)\displaystyle\phi(u,d,\rho,D,v) =0\displaystyle=0
    ϕ(π1,π2)\displaystyle\phi\left(\pi_{1},\pi_{2}\right) =0\displaystyle=0
    π1\displaystyle\pi_{1} =ua1db1ρc1D\displaystyle=u^{a_{1}}d^{b_{1}}\rho^{c_{1}}D
    π2\displaystyle\pi_{2} =ua2db2ρc2v\displaystyle=u^{a_{2}}d^{b_{2}}\rho^{c_{2}}v

    As each π\pi group is dimensionless then considering the dimensions, for the first group, π1\pi_{1} : (note D is a force with dimensions MLT2\mathrm{MLT}^{-2} )

    M0L0T0=(LT1)a1(L)b1(ML3)c1MLT2M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right)^{c% _{1}}MLT^{-2}

    M]

    0=c1+10=\mathrm{c}_{1}+1

    c1=1\mathrm{c}_{1}=-1
    L]

    0=a1+b13c1+10=a_{1}+b_{1}-3c_{1}+1

    T]

    0=a120=-a_{1}-2
    π1=u2d2ρ1D\pi_{1}=u^{-2}d^{-2}\rho^{-1}D
    =Dρu2d2=\frac{D}{\rho u^{2}d^{2}}

    And the second group π2\pi_{2} :

    M0L0T0=(LT1)a2(L)b2(ML3)c2L2T1\displaystyle M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{2}}(L)^{b_{2}}\left(ML^% {-3}\right)^{c_{2}}L^{2}T^{-1}
    0=c2\displaystyle 0=\mathrm{c}_{2}
    0=a2+b23c2+2\displaystyle 0=\mathrm{a}_{2}+\mathrm{b}_{2}-3\mathrm{c}_{2}+2
    2=a2+b2\displaystyle-2=\mathrm{a}_{2}+\mathrm{b}_{2}
    0=a21\displaystyle 0=-\mathrm{a}_{2}-1
    a2=1\displaystyle\mathrm{a}_{2}=-1
    b2=1\displaystyle\mathrm{\leavevmode\nobreak\ b}_{2}=-1
    π2=u1d1ρ0v\displaystyle\pi_{2}=u^{-1}d^{-1}\rho^{0}v
    =vud\displaystyle=\frac{v}{ud}

    T]

    So the physical situation is described by this function of nondimensional numbers,

    ϕ(π1,π2)=ϕ(Dρu2d2,vud)=0\phi\left(\pi_{1},\pi_{2}\right)=\phi\left(\frac{D}{\rho u^{2}d^{2}},\frac{v}{% ud}\right)=0

    For dynamic similarity these non-dimensional numbers are the same for the both the sphere in water and in the wind tunnel i.e.

    π1air =π1water\displaystyle\pi_{1_{\text{air }}}=\pi_{1_{\text{water }}}
    π2air =π2water\displaystyle\pi_{2_{\text{air }}}=\pi_{2_{\text{water }}}

    For π1\pi_{1}

    (Dρu2d2)air\displaystyle\left(\frac{D}{\rho u^{2}d^{2}}\right)_{\text{air }} =(Dρu2d2)water\displaystyle=\left(\frac{D}{\rho u^{2}d^{2}}\right)_{\text{water }}
    Dair 1.28×10.42×(2d)2\displaystyle\frac{D_{\text{air }}}{1.28\times 10.4^{2}\times(2d)^{2}} =41000×1.62×d2\displaystyle=\frac{4}{1000\times 1.6^{2}\times d^{2}}
    Dair\displaystyle D_{\text{air }} =0.865N\displaystyle=0.865\mathrm{\leavevmode\nobreak\ N}

    For π2\pi_{2}

    (vud)air\displaystyle\left(\frac{v}{ud}\right)_{\text{air }} =(vud)water\displaystyle=\left(\frac{v}{ud}\right)_{\text{water }}
    13vuair ×2d\displaystyle\frac{13v}{u_{\text{air }}\times 2d} =v1.6×d\displaystyle=\frac{v}{1.6\times d}
    uair\displaystyle u_{\text{air }} =10.4m/s\displaystyle=10.4\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
  2. E6.6.2

    Explain briefly the use of the Reynolds number in the interpretation of tests on the flow of liquid in pipes. Water flows through a 2 cm diameter pipe at 1.6m/s1.6\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}. Calculate the Reynolds number and find also the velocity required to give the same Reynolds number when the pipe is transporting air. Obtain the ratio of pressure drops in the same length of pipe for both cases. For the water the kinematic viscosity was 1.31×106m2/s1.31\times 10^{-6}\mathrm{\leavevmode\nobreak\ m}^{2}/\mathrm{s} and the density was 1000kg/m31000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}. For air those quantities were 15.1×106m2/s15.1\times 10^{-6}\mathrm{\leavevmode\nobreak\ m}^{2}/\mathrm{s} and 1.19kg/m31.19\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}.
    [24427, 18.4m/s,0.15718.4\mathrm{\leavevmode\nobreak\ m}/\mathrm{s},0.157 ]
    Draw up the table of values you have for each variable:

    variable water air
    u 1.6m/s1.6\mathrm{\leavevmode\nobreak\ m}/\mathrm{s} uair\mathrm{u}_{\text{air }}
    p pwater\mathrm{p}_{\text{water }} pair\mathrm{p}_{\text{air }}
    ρ\rho 1000kg/m31000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3} 1.19kg/m31.19\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}
    ν\nu 1.31×106m2/s1.31\times 10^{-6}\mathrm{\leavevmode\nobreak\ m}^{2}/\mathrm{s} 15.1×106m2/s15.1\times 10^{-6}\mathrm{\leavevmode\nobreak\ m}^{2}/\mathrm{s}
    ρ\rho 1000kg/m31000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3} 1.28kg/m31.28\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}
    d 0.02 m 0.02 m

    Kinematic viscosity is dynamic viscosity over density =ν=μ/ρ=\nu=\mu/\rho.
    The Reynolds number =Re=ρudμ=udv=\operatorname{Re}=\frac{\rho ud}{\mu}=\frac{ud}{v}
    Reynolds number when carrying water:

    Rewater =udv=1.6×0.021.31×106=24427\mathrm{Re}_{\text{water }}=\frac{ud}{v}=\frac{1.6\times 0.02}{1.31\times 10^{% -6}}=24427

    To calculate Reair \mathrm{Re}_{\text{air }} we know,

    Rewater\displaystyle\operatorname{Re}_{\text{water }} =Reair\displaystyle=\operatorname{Re}_{\text{air }}
    24427\displaystyle 24427 =uair 0.0215×106\displaystyle=\frac{u_{\text{air }}0.02}{15\times 10^{-6}}
    uair\displaystyle u_{\text{air }} =18.44m/s\displaystyle=18.44\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}

    To obtain the ratio of pressure drops we must obtain an expression for the pressure drop in terms of governing variables.
    Choose the three recurring (governing) variables; u,d,ρ\mathrm{u},\mathrm{d},\rho.
    From Buckinghams π\pi theorem we have mn=53=2\mathrm{m}-\mathrm{n}=5-3=2 non-dimensional groups.

    ϕ(u,d,ρ,v,p)\displaystyle\phi(u,d,\rho,v,p) =0\displaystyle=0
    ϕ(π1,π2)\displaystyle\phi\left(\pi_{1},\pi_{2}\right) =0\displaystyle=0
    π1\displaystyle\pi_{1} =ua1db1ρc1v\displaystyle=u^{a_{1}}d^{b_{1}}\rho^{c_{1}}v
    π2\displaystyle\pi_{2} =ua2db2ρc2p\displaystyle=u^{a_{2}}d^{b_{2}}\rho^{c_{2}}p

    As each π\pi group is dimensionless then considering the dimensions, for the first group, π1\pi_{1} :

    M0L0T0=(LT1)a1(L)b1(ML3)c1L2T1M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right)^{c% _{1}}L^{2}T^{-1}

    M] 0=c1\quad 0=\mathrm{c}_{1}
    L] 0=a1+b13c1+2\quad 0=\mathrm{a}_{1}+\mathrm{b}_{1}-3\mathrm{c}_{1}+2

    T]

    2=a1+10=a11a1=1b1=1π1=u1d1ρ0v=vud\begin{gathered}\begin{array}[]{c}-2=a_{1}+1\\ 0=-a_{1}-1\\ a_{1}=-1\\ b_{1}=-1\\ \pi_{1}=u^{-1}d^{-1}\rho^{0}v\\ =\frac{v}{ud}\end{array}\end{gathered}

    And the second group π2\pi_{2} :
    (note p is a pressure (force/area) with dimensions ML1T2\mathrm{ML}^{-1}\mathrm{\leavevmode\nobreak\ T}^{-2} )

    M0L0T0=(LT1)a1(L)b1(ML3)c1MT2L1M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right)^{c% _{1}}MT^{-2}L^{-1}

    M]

    0=c2+10=\mathrm{c}_{2}+1

    c2=1\mathrm{c}_{2}=-1
    L]

    T]

    0=a2+b23c212=a2+b20=a22a2=2b2=0π2=u2ρ1p=pρu2\begin{gathered}0=a_{2}+b_{2}-3c_{2}-1\\ -2=a_{2}+b_{2}\\ 0=-a_{2}-2\\ a_{2}=-2\\ b_{2}=0\\ \pi_{2}=u^{-2}\rho^{-1}p\\ =\frac{p}{\rho u^{2}}\end{gathered}

    So the physical situation is described by this function of nondimensional numbers,

    ϕ(π1,π2)=ϕ(vud,pρu2)=0\phi\left(\pi_{1},\pi_{2}\right)=\phi\left(\frac{v}{ud},\frac{p}{\rho u^{2}}% \right)=0

    For dynamic similarity these non-dimensional numbers are the same for the both water and air in the pipe.

    π1air =π1water\displaystyle\pi_{1_{\text{air }}}=\pi_{1_{\text{water }}}
    π2air =π2water\displaystyle\pi_{2_{\text{air }}}=\pi_{2_{\text{water }}}

    We are interested in the relationship involving the pressure i.e. π2\pi_{2}

    (pρu2)air\displaystyle\left(\frac{p}{\rho u^{2}}\right)_{\text{air }} =(pρu2)water\displaystyle=\left(\frac{p}{\rho u^{2}}\right)_{\text{water }}
    pwater pair\displaystyle\frac{p_{\text{water }}}{p_{\text{air }}} =ρwater uwater 2ρair uair 2\displaystyle=\frac{\rho_{\text{water }}u_{\text{water }}^{2}}{\rho_{\text{air% }}u_{\text{air }}^{2}}
    =1000×1.621.19×18.442=10.158=6.327\displaystyle=\frac{1000\times 1.6^{2}}{1.19\times 18.44^{2}}=\frac{1}{0.158}=% 6.327
  3. E6.6.3

    Show that Reynold number, ρud/μ\rho ud/\mu, is non-dimensional. If the discharge Q through an orifice is a function of the diameter d , the pressure difference p , the density ρ\rho, and the viscosity μ\mu, show that Q=Cp1/2d2/ρ1/2\mathrm{Q}=\mathrm{Cp}^{1/2}\mathrm{\leavevmode\nobreak\ d}^{2}/\rho^{1/2} where CC is some function of the non-dimensional group (dρ1/2d1/2/μ)\left(d\rho^{1/2}d^{1/2}/\mu\right).
    Draw up the table of values you have for each variable:
    The dimensions of these following variables are

    ρML3uLT1dLμML1T1Re=ML3LT1L(ML1T1)1=ML3LT1LM1LT=1\begin{array}[]{ll}\rho&\mathrm{ML}^{-3}\\ \mathrm{u}&\mathrm{LT}^{-1}\\ \mathrm{\leavevmode\nobreak\ d}&\mathrm{\leavevmode\nobreak\ L}\\ \mu&\mathrm{ML}^{-1}\mathrm{\leavevmode\nobreak\ T}^{-1}\\ \operatorname{Re}=&\mathrm{ML}^{-3}\mathrm{LT}^{-1}\mathrm{\leavevmode\nobreak% \ L}\left(\mathrm{ML}^{-1}\mathrm{\leavevmode\nobreak\ T}^{-1}\right)^{-1}=% \mathrm{ML}^{-3}\mathrm{LT}^{-1}\mathrm{\leavevmode\nobreak\ L}\mathrm{M}^{-1}% \mathrm{LT}=1\end{array}

    i.e. Re is dimensionless.

    We are told from the question that there are 5 variables involved in the problem: d,p,ρ,μ\mathrm{d},\mathrm{p},\rho,\mu and Q .
    Choose the three recurring (governing) variables; Q,d,ρ\mathrm{Q},\mathrm{d},\rho.
    From Buckinghams π\pi theorem we have mn=53=2\mathrm{m}-\mathrm{n}=5-3=2 non-dimensional groups.

    ϕ(Q,d,ρ,μ,p)\displaystyle\phi(Q,d,\rho,\mu,p) =0\displaystyle=0
    ϕ(π1,π2)\displaystyle\phi\left(\pi_{1},\pi_{2}\right) =0\displaystyle=0
    π1\displaystyle\pi_{1} =Qa1db1ρc1μ\displaystyle=Q^{a_{1}}d^{b_{1}}\rho^{c_{1}}\mu
    π2\displaystyle\pi_{2} =Qa2db2ρc2p\displaystyle=Q^{a_{2}}d^{b_{2}}\rho^{c_{2}}p

    As each π\pi group is dimensionless, then considering the dimensions, for the first group, π1\pi_{1} :

    M0L0T0=(L3T1)a1(L)b1(ML3)c1ML1T1M^{0}L^{0}T^{0}=\left(L^{3}T^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right% )^{c_{1}}ML^{-1}T^{-1}

    M]

    0=c1+1\displaystyle 0=c_{1}+1
    c1=1\displaystyle c_{1}=-1

    L]

    0=3a1+b13c110=3a_{1}+b_{1}-3c_{1}-1

    T]

    0=a110=-a_{1}-1

    a1=1\mathrm{a}_{1}=-1
    b1=1\mathrm{b}_{1}=1

    π1\displaystyle\pi_{1} =Q1d1ρ1μ\displaystyle=Q^{-1}d^{1}\rho^{-1}\mu
    =dμρQ\displaystyle=\frac{d\mu}{\rho Q}

    And the second group π2\pi_{2} :
    (note p is a pressure (force/area) with dimensions ML1T2\mathrm{ML}^{-1}\mathrm{\leavevmode\nobreak\ T}^{-2} )

    M0L0T0=(L3T1)a1(L)b1(ML3)c1MT2L1M^{0}L^{0}T^{0}=\left(L^{3}T^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right% )^{c_{1}}MT^{-2}L^{-1}

    M] 0=c2+1\quad 0=\mathrm{c}_{2}+1
    c2=1\mathrm{c}_{2}=-1
    L]

    0=3a2+b23c21\displaystyle 0=3\mathrm{a}_{2}+\mathrm{b}_{2}-3\mathrm{c}_{2}-1
    2=3a2+b2\displaystyle-2=3\mathrm{a}_{2}+\mathrm{b}_{2}

    T]

    0=a220=-a_{2}-2

    a2=2\mathrm{a}_{2}=-2

    b2=4\mathrm{b}_{2}=4
    π2\displaystyle\pi_{2} =Q2d4ρ1p\displaystyle=Q^{-2}d^{4}\rho^{-1}p
    =d4pρQ2\displaystyle=\frac{d^{4}p}{\rho Q^{2}}

    So the physical situation is described by this function of non-dimensional numbers,

    ϕ(π1,π2)=ϕ(dμQρ,d4pρQ2)=0\phi\left(\pi_{1},\pi_{2}\right)=\phi\left(\frac{d\mu}{Q\rho},\frac{d^{4}p}{% \rho Q^{2}}\right)=0

    or

    dμQρ=ϕ1(d4pρQ2)\frac{d\mu}{Q\rho}=\phi_{1}\left(\frac{d^{4}p}{\rho Q^{2}}\right)

    The question wants us to show : Q=f(dρ1/2p1/2μ)(d2p1/2ρ)Q=f\left(\frac{d\rho^{1/2}p^{1/2}}{\mu}\right)\left(\frac{d^{2}p^{1/2}}{\rho}\right)
    Take the reciprocal of square root of π2:1π2=ρ1/2Qd2p1/2=π2a\pi_{2}:\frac{1}{\sqrt{\pi_{2}}}=\frac{\rho^{1/2}Q}{d^{2}p^{1/2}}=\pi_{2a},
    Convert π1\pi_{1} by multiplying by this number

    π1a=π1π2a=dμQρρ1/2Qd2p1/2=μdρ1/2p1/2\pi_{1a}=\pi_{1}\pi_{2a}=\frac{d\mu}{Q\rho}\frac{\rho^{1/2}Q}{d^{2}p^{1/2}}=% \frac{\mu}{d\rho^{1/2}p^{1/2}}

    then we can say

    ϕ(1/π1a,π2a)=ϕ(p1/2ρ1/2dμ,d2p1/2ρ1/2)=0\phi\left(1/\pi_{1a},\pi_{2a}\right)=\phi\left(\frac{p^{1/2}\rho^{1/2}d}{\mu},% \frac{d^{2}p^{1/2}}{\rho^{1/2}}\right)=0

    or

    Q=ϕ(p1/2ρ1/2dμ)d2p1/2ρ1/2Q=\phi\left(\frac{p^{1/2}\rho^{1/2}d}{\mu}\right)\frac{d^{2}p^{1/2}}{\rho^{1/2}}
  4. E6.6.4

    A cylinder 0.16 m in diameter is to be mounted in a stream of water in order to estimate the force on a tall chimney of 1 m diameter which is subject to wind of 33m/s33\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}. Calculate (A) the speed of the stream necessary to give dynamic similarity between the model and chimney, (b) the ratio of forces.
    Chimney: ρ=1.12kg/m3μ=16×106kg/ms\quad\rho=1.12\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}\quad\mu=16\times 1% 0^{-6}\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}
    Model: ρ=1000kg/m3μ=8×104kg/ms\quad\rho=1000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}\quad\mu=8\times 1% 0^{-4}\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}
    [11.55m/s,0.057][11.55\mathrm{\leavevmode\nobreak\ m}/\mathrm{s},0.057]
    Draw up the table of values you have for each variable:

    variable water air
    u uwater\mathrm{u}_{\text{water }} 33m/s33\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    F Fwater\mathrm{F}_{\text{water }} Fair\mathrm{F}_{\text{air }}
    ρ\rho 1000kg/m31000\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3} 1.12kg/m31.12\mathrm{\leavevmode\nobreak\ kg}/\mathrm{m}^{3}
    μ\mu 8×104kg/ms8\times 10^{-4}\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms} 16×106kg/ms16\times 10^{-6}\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}
    d 0.16 m 1 m

    Kinematic viscosity is dynamic viscosity over density =ν=μ/ρ=\nu=\mu/\rho.
    The Reynolds number =Re=ρudμ=udv=\operatorname{Re}=\frac{\rho ud}{\mu}=\frac{ud}{v}
    For dynamic similarity:

    Rewater\displaystyle\mathrm{Re}_{\text{water }} =Reair\displaystyle=\mathrm{Re}_{\text{air }}
    1000uwater 0.168×104\displaystyle\frac{1000u_{\text{water }}0.16}{8\times 10^{-4}} =1.12×33×116×106\displaystyle=\frac{1.12\times 33\times 1}{16\times 10^{-6}}
    uwater\displaystyle u_{\text{water }} =11.55m/s\displaystyle=11.55\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}

    To obtain the ratio of forces we must obtain an expression for the force in terms of governing variables.
    Choose the three recurring (governing) variables; u,d,ρ,F,μ\mathrm{u},\mathrm{d},\rho,\mathrm{F},\mu.
    From Buckinghams π\pi theorem we have mn=53=2\mathrm{m}-\mathrm{n}=5-3=2 non-dimensional groups.

    ϕ(u,d,ρ,μ,F)\displaystyle\phi(u,d,\rho,\mu,F) =0\displaystyle=0
    ϕ(π1,π2)\displaystyle\phi\left(\pi_{1},\pi_{2}\right) =0\displaystyle=0
    π1\displaystyle\pi_{1} =ua1db1ρc1μ\displaystyle=u^{a_{1}}d^{b_{1}}\rho^{c_{1}}\mu
    π2\displaystyle\pi_{2} =ua2db2ρc2F\displaystyle=u^{a_{2}}d^{b_{2}}\rho^{c_{2}}F

    As each π\pi group is dimensionless then considering the dimensions, for the first group, π1\pi_{1} :

    M0L0T0=(LT1)a1(L)b1(ML3)c1ML1T1M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right)^{c% _{1}}ML^{-1}T^{-1}

    M] 0=c1+1\quad 0=\mathrm{c}_{1}+1

    c1=1c_{1}=-1

    L] 0=a1+b13c11\quad 0=\mathrm{a}_{1}+\mathrm{b}_{1}-3\mathrm{c}_{1}-1

    2=a1+b1-2=a_{1}+b_{1}

    T]

    0=a110=-a_{1}-1

    a1=1\mathrm{a}_{1}=-1
    b1=1\mathrm{b}_{1}=-1

    π1=u1d1ρ1μ\pi_{1}=u^{-1}d^{-1}\rho^{-1}\mu
    =μρud=\frac{\mu}{\rho ud}

    i.e. the (inverse of) Reynolds number
    And the second group π2\pi_{2} :

    M0L0T0=(LT1)a2(L)b2(ML3)c2ML1T2M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{2}}(L)^{b_{2}}\left(ML^{-3}\right)^{c% _{2}}ML^{-1}T^{-2}

    M]

    0=c2+10=\mathrm{c}_{2}+1

    c2=1\mathrm{c}_{2}=-1
    L]

    0=a2+b23c21\displaystyle 0=a_{2}+b_{2}-3c_{2}-1
    3=a2+b2\displaystyle-3=a_{2}+b_{2}

    T]

    0=a220=-a_{2}-2
    π2=u2d1ρ1F\pi_{2}=u^{-2}d^{-1}\rho^{-1}F
    =Fu2dρ=\frac{F}{u^{2}d\rho}

    So the physical situation is described by this function of nondimensional numbers,

    ϕ(π1,π2)=ϕ(μρud,Fρdu2)=0\phi\left(\pi_{1},\pi_{2}\right)=\phi\left(\frac{\mu}{\rho ud},\frac{F}{\rho du% ^{2}}\right)=0

    For dynamic similarity these non-dimensional numbers are the same for the both water and air in the pipe.

    π1air =π1water\displaystyle\pi_{1_{\text{air }}}=\pi_{1_{\text{water }}}
    π2air =π2water\displaystyle\pi_{2_{\text{air }}}=\pi_{2_{\text{water }}}

    To find the ratio of forces for the different fluids use π2\pi_{2}

    π2air\displaystyle\pi_{2_{\text{air }}} =π2water\displaystyle=\pi_{2_{\text{water }}}
    (Fρu2d)air\displaystyle\left(\frac{F}{\rho u^{2}d}\right)_{\text{air }} =(Fρu2d)water\displaystyle=\left(\frac{F}{\rho u^{2}d}\right)_{\text{water }}
    (Fρu2d)air\displaystyle\left(\frac{F}{\rho u^{2}d}\right)_{\text{air }} =(Fρu2d)water\displaystyle=\left(\frac{F}{\rho u^{2}d}\right)_{\text{water }}
    Fair Fwater\displaystyle\frac{F_{\text{air }}}{F_{\text{water }}} =1.12×332×11000×11.552×0.16=0.057\displaystyle=\frac{1.12\times 33^{2}\times 1}{1000\times 11.55^{2}\times 0.16% }=0.057
  5. E6.6.5

    If the resistance to motion, RR, of a sphere through a fluid is a function of the density ρ\rho and viscosity μ\mu of the fluid, and the radius rr and velocity uu of the sphere, show that RR is given by

    R=μ2ρf(ρurμ)R=\frac{\mu^{2}}{\rho}f\left(\frac{\rho ur}{\mu}\right)

    Hence show that if at very low velocities the resistance R is proportional to the velocity u , then R=kμru\mathrm{R}=\mathrm{k}\mu\mathrm{ru} where k is a dimensionless constant.
    A fine granular material of specific gravity 2.5 is in uniform suspension in still water of depth 3.3 m .
    Regarding the particles as spheres of diameter 0.002 cm find how long it will take for the water to clear.
    Take k=6π\mathrm{k}=6\pi and μ=0.0013kg/ms\mu=0.0013\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}.
    [218mins 39.3 sec ]
    Choose the three recurring (governing) variables; u,r,ρ,R,μu,r,\rho,R,\mu.
    From Buckinghams π\pi theorem we have mn=53=2\mathrm{m}-\mathrm{n}=5-3=2 non-dimensional groups.

    ϕ(u,r,ρ,μ,R)\displaystyle\phi(u,r,\rho,\mu,R) =0\displaystyle=0
    ϕ(π1,π2)\displaystyle\phi\left(\pi_{1},\pi_{2}\right) =0\displaystyle=0
    π1\displaystyle\pi_{1} =ua1rb1ρc1μ\displaystyle=u^{a_{1}}r^{b_{1}}\rho^{c_{1}}\mu
    π2\displaystyle\pi_{2} =ua2rb2ρc2R\displaystyle=u^{a_{2}}r^{b_{2}}\rho^{c_{2}}R

    As each π\pi group is dimensionless then considering the dimensions, for the first group, π1\pi_{1} :

    M0L0T0=(LT1)a1(L)b1(ML3)c1ML1T1M^{0}L^{0}T^{0}=\left(LT^{-1}\right)^{a_{1}}(L)^{b_{1}}\left(ML^{-3}\right)^{c% _{1}}ML^{-1}T^{-1}

    M]

    0=c1+1\displaystyle 0=c_{1}+1
    c1=1\displaystyle c_{1}=-1

    L]

    0=a1+b13c110=\mathrm{a}_{1}+\mathrm{b}_{1}-3\mathrm{c}_{1}-1
    2=a1+b1-2=a_{1}+b_{1}

    T]

    0=a110=-a_{1}-1

    a1=1\mathrm{a}_{1}=-1

    b1=1\mathrm{b}_{1}=-1
    π1\displaystyle\pi_{1} =u1r1ρ1μ\displaystyle=u^{-1}r^{-1}\rho^{-1}\mu
    =μρur\displaystyle=\frac{\mu}{\rho ur}

    i.e. the (inverse of) Reynolds number

    And the second group π2\pi_{2} :

    M0L0T0\displaystyle M^{0}L^{0}T^{0} =(LT1)a2(L)b2(ML3)c2ML1T2\displaystyle=\left(LT^{-1}\right)^{a_{2}}(L)^{b_{2}}\left(ML^{-3}\right)^{c_{% 2}}ML^{-1}T^{-2}
    0\displaystyle 0 =c2+1\displaystyle=\mathrm{c}_{2}+1

    L] 0=a2+b23c21\quad 0=\mathrm{a}_{2}+\mathrm{b}_{2}-3\mathrm{c}_{2}-1

    T]

    3=a2+b20=a22a2=2b2=1π2=u2r1ρ1R=Ru2rρ\begin{gathered}-3=a_{2}+b_{2}\\ 0=-a_{2}-2\\ a_{2}=-2\\ b_{2}=-1\\ \pi_{2}=u^{-2}r^{-1}\rho^{-1}R\\ =\frac{R}{u^{2}r\rho}\end{gathered}

    So the physical situation is described by this function of nondimensional numbers,

    ϕ(π1,π2)=ϕ(μρur,Rρru2)=0\phi\left(\pi_{1},\pi_{2}\right)=\phi\left(\frac{\mu}{\rho ur},\frac{R}{\rho ru% ^{2}}\right)=0

    or

    Rρru2=ϕ1(μρur)\frac{R}{\rho ru^{2}}=\phi_{1}\left(\frac{\mu}{\rho ur}\right)

    he question asks us to show R=μ2ρf(ρurμ)R=\frac{\mu^{2}}{\rho}f\left(\frac{\rho ur}{\mu}\right) or Rρμ2=f(ρurμ)\frac{R\rho}{\mu^{2}}=f\left(\frac{\rho ur}{\mu}\right)
    Multiply the LHS by the square of the RHS: (i.e. π2×(1/π1)2\pi_{2}\times\left(1/\pi_{1}{}^{2}\right) )

    Rρru2×ρ2u2r2μ2=Rρμ2\frac{R}{\rho ru^{2}}\times\frac{\rho^{2}u^{2}r^{2}}{\mu^{2}}=\frac{R\rho}{\mu% ^{2}}

    So

    Rρμ2=f(ρurμ)\frac{R\rho}{\mu^{2}}=f\left(\frac{\rho ur}{\mu}\right)

    The question tells us that R is proportional to u so the function ff must be a constant, kk

    Rρμ2\displaystyle\frac{R\rho}{\mu^{2}} =kρurμ\displaystyle=k\frac{\rho ur}{\mu}
    R\displaystyle R =μkru\displaystyle=\mu kru

    The water will clear when the particle moving from the water surface reaches the bottom. At terminal velocity there is no acceleration - the force R=mg\mathrm{R}=\mathrm{mg} - upthrust.
    From the question:

    σ=2.5 so ρ=2500kg/m3μ=0.0013kg/msk=6π\sigma=2.5\quad\text{ so }\quad\rho=2500\mathrm{\leavevmode\nobreak\ kg}/% \mathrm{m}^{3}\quad\mu=0.0013\mathrm{\leavevmode\nobreak\ kg}/\mathrm{ms}\quad% \mathrm{k}=6\pi
    r=0.00001m depth =3.3m\mathrm{r}=0.00001\mathrm{\leavevmode\nobreak\ m}\quad\text{ depth }=3.3% \mathrm{\leavevmode\nobreak\ m}
    mg=43π0.000013×9.81×(25001000)\displaystyle mg=\frac{4}{3}\pi 0.00001^{3}\times 9.81\times(2500-1000)
    =6.16×1011\displaystyle=6.16\times 10^{-11}
    μkru\displaystyle\mu kru =0.0013×6π×0.00001u=6.16×1011\displaystyle=0.0013\times 6\pi\times 0.00001u=6.16\times 10^{-11}
    u\displaystyle u =2.52×104m/s\displaystyle=2.52\times 10^{-4}\mathrm{\leavevmode\nobreak\ m}/\mathrm{s}
    t=\displaystyle t= 3.32.52×104=218min39.3sec\displaystyle\frac{3.3}{2.52\times 10^{-4}}=218\mathrm{\leavevmode\nobreak\ % min}39.3\mathrm{sec}