ENGI 5432 Advanced Calculus

    Faculty of Engineering and Applied Science
    2009 Winter

    Problem Set 6   -   Solutions

    1. Find the centroid (centre of mass when the surface density is constant) of the frustum of the circular cone   z2 = x2 + y2   between the origin and the plane   z = a, (where   a   is a positive constant).


      Surface density   r = constant
      Use parameters (q, z)
      (0 £ q £ 2p,     0 £ z £ a )
      vector r  =  < z cos q, z sin q, z >
      Normal  =  dr/dtheta cross dr/dz
      = det [...]

      [Cross-section of cone]

      = ...   = z < cos theta, sin theta, -1 >
      scalar N  =  z * sqrt{2}
      m  =  Integral rho dS  =  ...   = ...
      = ...   = rho sqrt{2} (a^2 / 2) (2 pi)
      mass = pi rho a^2 sqrt{2}

      Taking moments about the   x-y plane:
      DMz = z Dm   Moment z  =  Integral z rho dS   = ...
      =  rho sqrt{2} Integ z^2 dz Integ dtheta   = ...   = ...
      Mz  =  (2/3) pi rho a^3 sqrt{2}   zBar  =  Mz / m  =  2a/3
      By symmetry,  xBar = yBar = 0
      Therefore the centroid of the frustum of the cone is at   (0, 0, (2/3)a )


    1. Find the total flux that passes through the hemisphere   x2 + y2 + z2 = a2,   z > 0 ,   due to the electric flux density D = eE = Q / (4*pi*r^3) * vector-r   where  r  =  x i + y j + z k   using:

      1. the surface projection method   and


        S:     x 2 + y 2 + z 2 = a 2,     z > 0   z  =  +sqrt{a^2 - x^2 - y^2}
        Partial dz/dx  =  - x/z
        Partial dz/dy  =  - y/z     (by symmetry)
        Normal  =  ...   ... = (1/z) * vector r
        D dot N  =  Qr^2 / (4 pi r^3 z)
        But everywhere on the hemisphere   r = a
        D dot N  =  Q / (4 pi a z)
==>  Flux  =  Integral {Q / (4 pi a z)} dA
        Switch to plane polar coordinates, in which   dA = r dr df.
        x = r cos f,     y = r sin f,     Þ     x 2 + y 2 = r 2

        [Projection on x-y plane]

        z  =  sqrt{a^2 - rho^2}   everywhere on the hemisphere
        Flux  =  Integral ...    = ...
        ... = (Q / (4 pi a) * 2 pi * a

        Flux  =  Q/2


      2. the parametric surface net method.


        Use the parameters   (q, f)       (0 £ q £ p/2,     0 £ f < 2p).
        r  =  < x, y, z>   in terms of theta, phi

        [3D view of hemisphere]

        N  =  dr/dtheta cross dr/dphi
        [Note:   this vector   N   has a different magnitude from the vector   N   in part (a) of this question.]
        = det [...]
        = ...
        = ...
        N  =  a^2 sin theta * vector r-hat
        But everywhere on the hemisphere   r = a   D  =  (Q / (4 pi a^3) * vector r
        D dot N = ... = Q sin theta / (4 pi)
        Integral D dot dS = ...
        = ...   = ...   =  Q / 2


      3. What can you say about the relative orientation of the vector field D and the normal N to the hemisphere?


        From part (b) above, vectors   D   and   N   are both multiples of the displacement vector   r.   D   and   N   are therefore parallel to each other everywhere on the hemisphere.


    1. Find the total flux of the vector field F  =  (x/2) a-hat_x through that part of the circular paraboloid   z = 4 - x2 - y2   that lies in the first octant.


      Projection method:
      z = 4 - x 2 - y 2     (first octant only)
      N  =  < 2x, 2y, 1 >
      F dot N  =  x^2
      Use plane polar coordinates
            x = r cos f,     y = r sin f,       Þ     x 2 + y 2 = r 2
      Þ     z = 4 - r 2     and
      Integral F dot N dS = Integral x^2 dA = ...
      = ...
      = ...
      Therefore

      Net flux through   S = p.

      or

      Surface parametric net method:
      Use parameters   r   and   f:
            x = r cos f,     y = r sin f,     z = 4 - r 2
      N  =  det [...]   = < 2r^2 cos f,  2r^2 sin f, r >
      F dot N  =  ...  =  rho^3 cos^2 phi
      Flux = ... = pi  (as before)


    1. A thin sheet, in the shape of the paraboloid   z = 16 - (x2 + y2), where   x > 0   and   y > 0, has a surface density of   rho  =  x*y / sqrt{1 + 4(x^2 + y^2)}   It lies within the coaxial region between the circular cylinders   C1:  x^2 + y^2 = 1, z any real #   and   C2:  x^2 + y^2 = 9, z any real #

      1. Determine the total mass of the sheet, and

        In order to obtain a better visualization of the surface S, examine its projection onto any vertical plane containing the z-axis and its projection onto the x-y plane:

        [Vertical cross-section]   [Projection on x-y plane]

        Let   r 2 = x 2 + y 2,   then the paraboloid is
              z = 16 - r 2
        Cylinder   C1   is     r 2 = 1
        Cylinder   C2   is     r 2 = 9
        The sheet is therefore defined by
              z = 16 - r 2,     7 £ z £ 15,     x ³ 0 ,     y ³ 0.

        Surface parametric net method:

        Use parameters   (r, f):
        vector r  =  < r cos f, r sin f, 16 - r^2 >
        1 £ r £ 3 ,     0 £ f £ p/2
        The surface density then becomes
        rho  =  r^2 cos f sin f / sqrt{1 + 4r^2}
        N  =  dR/dr cross dR/dphi
        = det [...]
        = ...   = r < 2r cos phi, 2r sin phi, 1 >
        N  =  r sqrt{1 + 4r^2}
        Total mass:
        mass  =  Integral rho dS  =  ...
        = ... =  (1+1)/4 * (81-1)/4
        Therefore (assuming S.I. units)

        Total mass = 10 kg

        Projection Method:

        z = 16 - (x 2 + y 2)
        Partial dz/dx = -2x, Partial dz/dy = -2y
        vector N  =  < 2x, 2y, 1 >
        scalar N  =  sqrt{1 + 4(x^2 + y^2)}
        mass  =  Integral  rho N dA  =  Integral xy dA
        Convert to plane polar coordinates:
        m  =  Integral ...
        = ... =  (1+1)/4 * (81-1)/4
        Therefore (assuming S.I. units)

        Total mass = 10 kg


      2. Locate the centre of mass for the sheet.

        Elements of first moments:
        DMx   =   x . Dm   =   x r DS   =   x r N Dr Df
        = ...
        Mx  =  Integral r^4 sin f cos^2 f dr df
        = ... = 242/15
        xBar  =  121/75  =  1.613
        By symmetry,   yBar  =  xBar

        DMz   =   z . Dm   =   z r DS   =   z r N Dr Df
        = ...
        Mz  =  Integral (16r^3 - r^5) sin f cos f dr df
        = ...
        = ...   = ...
        Mz  =  298 / 3
        zBar  = ... =  149/15  =  9.933
        Therefore the centre of mass is at

        (121/75, 121/75, 149/15)

        Or, to 3 d.p.,

        (1.613,   1.613,   9.933)

    1. Fluid is flowing along a cylindrical pipe.   The circular cross section inside the pipe has a constant radius of a (m).   The Cartesian coordinate system is aligned with the x axis along the line of symmetry (centre line) of the pipe.   The other two axes are in the plane of one of the circular cross sections.   The velocity of fluid at any point (x, y, z) inside the pipe is
                    v = v0 / a * sqrt{a^2 - y^2 - z^2} iHat
      Find the flux   Q (in m3/s)   (the rate at which fluid is flowing through the pipe).


      The x-z and y-z planes are illustrated here, together with the velocity profile in the former case:
      [side and end views of the pipe]
      The speed increases from zero on the wall of the pipe to a maximum of   v 0   along the centre line.
      The flux   Q   is the flux of the velocity vector   v   of the fluid across a circular cross section in the direction of a unit normal to that cross section.   The unit normal is
      N^ = i^   ==>   
     v dot N^ = v0 / a * sqrt{a^2 - y^2 - z^2}
      Q = Integral over A of v0 / a * sqrt{a^2 - y^2 - z^2}
      Transform to plane polar coordinates:
      y = r cos theta,  z = r sin theta
      Q = v0 / a * Int Int sqrt{a^2 - r^2} r dr d_theta
      Q = 2 pi / 3 v0 a^2

      Q = 2 pi / 3 v0 a^2

      [We can also verify that this answer is dimensionally consistent:
      the units on the right side are [m s–1] × [m]2 = [m3 s–1]
      which are the units for volume per unit time.]


    1. Complete the evaluation (in Example 2.6.3 of the lecture notes) of the line integral   line integral F.dr   around the unit square in the x-z plane for F = < xyz, xz, e^(xy) >, (without using Stokes’ theorem).


      [unit square in the x-z plane] From example 2.6.3, the positive orientation around the square is   OGHJ (the y axis is directed into the page).
      In the x-z plane y = 0   vector F  =  < 0, xz, 1 >
      Compute the line integral around the four sides of the square, starting with OG:
      r = < 0, 0, t >  ==>  dr/dt = < 0, 0, 1 >
      F = < 0, 0, 1 >  ==>  F dot dr/dt = 1
      line integral(OG) = 1
      The above was presented in a lecture.
      Below are the calculations for the other three sides of the square:
      GH:  r = < t, 0, 1 >  ==>  dr/dt = < 1, 0, 0 >
      F = < 0, t, 1 >  ==>  F dot dr/dt = 0
      line integral(GH) = 0
      HJ:  r = < 1, 0, -t >  ==> dr/dt = < 0, 0, -1 >
      F = < 0, -t, 1 >  ==>  F dot dr/dt = -1
      line integral(HJ) = -1
      JO:  r = < -t, 0, 0 >  ==> dr/dt = < -1, 0, 0 >
      F = < 0, 0, 1 >  ==>  F dot dr/dt = 0
      line integral(JO) = 0
      Total line integral = 0
      Therefore

      Total line integral = 0


    1. Find the mass and the location of the centre of mass of the upper half of the ellipsoid, whose equation in Cartesian coordinates is
      x2 + y2 + 4z2   =   4 ,       z > 0
      and whose surface density is
      rho = 1 / sqrt{1 + 3z^2}
      1. using a parametric net method. and

        Use a parametric net (q, f) that is a close relative of the spherical polar coordinate system:
        x = 2 sin q cos f ,       y = 2 sin q sin f ,       z = cos q
        Þ     x2 + y2   =   4 sin2q
        and     x2 + y2 + 4z2   =   4 sin2q   + 4 cos2q   =   4
        The parameter ranges are     0 < q < p/2 ,     0 < f < 2p
        Employ the abbreviations     cq = cos q ,   sq = sin q ,   cf = cos f ,   sf = sin f
        N = det[
    i     j     k
  2ctcf 2ctsf  -st
 -2stsf 2stcf   0  ]
        = 2st < stcf, stsf, 2ct >
        N = 2st sqrt{1 + 3 ct^2}
        density = 1/sqrt{1 + 3z^2} = 1/sqrt{1 + 3 ct^2}
        m = Integral{rho dS} = ...
        = 2[phi]_0^2pi [-cos theta]_0^pi/2
        Therefore

        m = 4 pi

        By symmetry,   xBar = yBar = 0
        Taking moments about the xy-plane:
        Mz = Integral{rho z dS} = ...
        = 2[phi]_0^2pi [-(cos 2theta)/2]_0^pi/2
         = 2pi
        zBar = Mz / m = 1/2
        Therefore the centre of mass is at the point

        (0, 0, 1/2)


      2. using a projection method.

        The shadow of the ellipsoid on the xy-plane is the area bounded by the circle   x2 + y2 = 4.
        Everywhere on the ellipsoid,   z^2  =  (4 - (x^2 + y^2))/4
        ==>  dz/dx  =  -x/(4z)
        and  dz/dy  =  -y/(4z)
        N  =  sqrt{(dz/dx)^2 + (dz/dy)^2 + 1}
    =  sqrt{(x/(4z))^2 + (y/(4z))^2 + 1}
        =  sqrt{((x^2 + y^2 + 4z^2) + 12z^2)/(16z^2)}
        =  sqrt{1 + 3z^2}/(2z)
        m  =  Integral{rho dS}
         = Integral{1/sqrt{1+3z^2} * sqrt{1+3z^2}/(2z) dx dy}
        But  2z = sqrt{4 - x^2 - y^2}
         = sqrt{4 - r^2} ,  where  r^2 = x^2 + y^2
        Switch to plane polar coordinates (r, q):
        m = Integral {r / sqrt{4-r^2} dr dtheta}
     =  [ theta ]_0^2pi [-sqrt{4-r^2}]_0^2  =  2pi(0 + 2)
        Therefore

        m = 4 pi

        By symmetry,   xBar = yBar = 0
        Taking moments about the xy-plane:
        Mz = Integral{rho z dS} =
     = Integral{1/sqrt{1+3z^2} * z * sqrt{1+3z^2}/(2z) dA}
        = Integral{1/2 dA}  =  A/2
        But the projection is the region bounded by a circle of radius 2.
        Þ     A   =   p (2)2
        Mz = (4pi)/2 ==>
         zBar = Mz / m = 1/2
        Therefore the centre of mass is at the point

        (0, 0, 1/2)


    Return to Index of Problem Set Solutions   [Return to the Index of Solutions.]
    Back to previous page   [Return to your previous page]
    Created 2008 01 11 and most recently modified 2008 12 27 by Dr. G.H. George