ENGI 5432 Advanced Calculus

    Faculty of Engineering and Applied Science
    2009 Winter

    Problem Set 5   -   Solutions

    1. Find the work done when an object is moved along the curve of intersection   C   of the circular paraboloid   z = x 2 + y 2   and the plane   2x + y = 2   from (1, 0, 1) to (0, 2, 4) by a force F  =  < 2/x, 1/x, 1 >.


      2x + y = 2       Þ     y = 2(1 - x)
      Let   x = t   then   y = 2(1 - t)     (from the equation of the plane)
      and   z = t2 + 4(1 - t) 2     (from the equation of the paraboloid)

      The value of   t   at   (1, 0, 1)   is   1.
      The value of   t   at   (0, 2, 4)   is   0.

      Therefore, on the curve   C,
      vector r  =  < t, 2 - 2t, 5t^2 - 8t + 4 >
     ==>  dr/dt = < 1, -2, 10t - 8 >
      vector F  =  > 2/t, 1/t, 1 <
      F dot dr/dt  =  10t - 8
      work done = Integ_1^0 {10t - 8} dt = 3
      Note the order of the limits of integration.
      The path is from the point where t = 1 to the point where t = 0.
      Therefore

        W = 3 Nm  

      [Note:
      curl F  not=  0
      Therefore   F   is not conservative and the value of   Integral_C {F dot dr}   depends on the path taken.]


    1. Find the work done in travelling in   real 2 space   once anti-clockwise around the unit circle   C, centered on the origin, in the presence of the force   F  =  < x - y,  x + y >, using both of the following methods:

      1. by direct evaluation of the line integral Integral_C {F dot dr} and

        A parameterization for the unit circle   C   is
        vector r = < cos t, sin t >
        [values of dr/dt  and F]
        F dot dr/dt = 1  for all t
        line integral = integ_0^2pi 1 dt

        Line integral  =  2 pi


      2. by using Green’s theorem and evaluating Integ Integ {df2/dx - df1/dy} dA.


        df2/dx - df1/dy  =  ...  =  2
        Double integral  =  ...  =  2 A  =  2 pi
        F   has no singularities anywhere on   real 2 space, so, by Green’s theorem,

        Line integral  =  2 pi


      3. Is the vector field   F   conservative?

        Line integral not= 0     Therefore

          NO,     F is not conservative.  


    1. In class we saw that the line integral   Integral_C {F dot dr}   has the value   2p, where   F  =  < -y / (x^2+y^2), x / (x^2+y^2) >   and the path   C   is one anti-clockwise circuit of the unit circle centered on the origin.

      1. Show that the value of the line integral remains   2p   when the vector field is replaced by
        F = < -y / sqrt{x^2+y^2}^n, x / sqrt{x^2+y^2}^n >

        F  =  < -y/r^n, +x/r^n >
        A parameterization for the unit circle   C   is
        r  =  < cos t, sin t >
        r = 1   for all t   and
         r^n = 1  for all n
        dr/dt  =  < -s, c > = F
        F . dr/dt  =  1   for all t
        Line integral  =  2 pi


      2. Show that the value of   Integ Integ {df2/dx - df1/dy} dA     (where   D   is the circular area enclosed by   C) is not 2p   unless   n < 2.

        Hint:   transform the area integral into plane polar coordinates, with
        x = r cos q,       y = r sin q       and     dA = dx dy = r dr dq.


        r = sqrt{x^2 + y^2}   ==>  dr/dx = x/r
        By symmetry,   dr/dy  =  y/r
        df2/dx  =  (r^2 - n*x^2)/r^(n+2)
        By symmetry,   df1/dy  =  -(r^2 - n*y^2)/r^(n+2)
        Integrand  =  (2-n)*r^(-n)
        Double integral = ...
        Integ 1 dtheta × Integ (2-n)r^(1-n) dr
        Integ 1 dtheta = 2 pi
        Integ (2-n)r^(1-n) dr = 1 for n < 2 only
        Integral  =  { 2 pi (n<2);  0 (n=2);  undefined else}

        Therefore the integral ¹ 2p   unless n < 2.


      3. Parts (a) and (b) clearly show that Green’s theorem does not apply when n > 2.
        Why should Green’s theorem not be used when   n > 0 ?

        Green’s theorem should not be used if   C   or   D   contain any singularities.
        There is a singularity at the origin (r = 0) if   n > 0.
        There are no singularities anywhere in real 2 space otherwise.
        Therefore do not use Green’s theorem for   n > 0.


    1. Find a potential function for the vector field   F  =  < 2xy cos z, x^2 cos z, -x^2y sin z >
      Use this potential function to evaluate the line integral Integral_C {F dot dr} for any piecewise smooth simple curve   C   from the point (1, 1, -1) to the point (2, 0, 5).


      F  =  < 2xy cos z, x^2 cos z, -x^2y sin z >
      Match this with   grad phi  =  < d phi/dx, d phi/dy, d phi/dz >
      One of several valid solutions is:
      d phi/dx = 2 xy cos z    
    ==>  phi = x^2 y cos z  +  g(y,z)
      where   g(y, z)   is an arbitrary function of   y   and   z.
      This expression for   f   leads to   d phi/dy  =  x^2 cos z  +  dg/dy
      But  d phi/dy  =  x^2 cos z   ==>  dg/dy = 0
      Þ     g   is a function of   z   alone,     g(z).
      d phi/dz  =  -x^2 y sin z  +  dg/dz
      But  d phi/dz  =  -x^2 y sin z  ==>  dg/dz = 0
      Þ     g = A, an arbitrary constant.     Therefore the potential function is

        f (x, y, z)   =   x2y cos z   +   A  

      [The fact that we were able to find a potential function     Þ     F   is conservative.]

      C   is any path from   (1, 1, -1)   to   (2, 0, 5).
      Line integral = phi(2, 0, 5) - phi(1, 1, -1)
      = (22×0×cos 5 + A)   -   (12×1×cos(-1) + A)       Therefore

      Line integral  =  - cos 1  approx= .540


    1. In the evaluation of an area integral   Double integral f(A) dA   or a volume integral   Triple integral f(V) dV,   the Cartesian differentials are related to the differentials of a different parameterization (u, v, w) of the surface or volume by the Jacobian determinant:

      dA = dx dy = |partial(x,y) / partial(u,v)| du dv
      and   dV = dx dy dz = |partial(x,y,z) / partial(u,v,w)| du dv dw,
      where the Jacobian is
      det
dx/du  dy/du
dx/dv  dy/dv       
or  det 
dx/du  dy/du  dz/du
dx/dv  dy/dv  dz/dv
dx/dw  dy/dw  dz/dw.

      Evaluate the Jacobian in order to determine the conversion formula from Cartesian differentials   dx dy   (or   dx dy dz) to

      1. Plane polar area differentials   dr dq

        x = r cos q ,       y = r sin q
        Jacobian = r

        dA  =  r dr dtheta


      2. Cylindrical polar volume differentials   dr   df   dz

        x = r cos f,     y = r sin f,     z = z
        Jacobian = rho

        dV  =  rho drho dphi dz


      3. Spherical polar volume differentials   dr dq   df

        x = r sin q cos f,     y = r sin q sin f,     z = r cos q
        [Jacobian determinant]
        =   | cos q (r2 sin q cos q (cos2f + sin2f) ) + r sin q (r sin2q (cos2f + sin2f) ) + 0 |
        =   | r2 sin q (cos2q + sin2q) |   =   r2 sin q
        [Note that   sin q ³ 0   because   0 £ q £ p]

        dV  =  r^2sin theta dr dtheta dphi


      4. Parabolic area differentials   du dv, where   x = u v,   y = v 2 - u 2.


        x = uv,     y = v 2 - u 2
        Jacobian  =  2(u^2 + v^2)

        dA  =  2(u^2 + v^2) du dv


      5. In part (d), sketch on the same x-y plane any three members of each of the two families of coordinate curves   u = constant   and   v = constant.


        y = v^2 - (x/v)^2  and  y = (x/u)^2 - u^2
        Choosing (arbitrarily) the values   1/2, 1 and 2   for each of   u   and   v:
        Parabolas, up, vertices (0, -u^2)
        Parabolas, down, vertices  (0, +v^2)
        [Note that   u = 0   is the positive y-axis and   v = 0   is the negative y-axis.
        u = c   is the same curve as   u = -c.
        v = c   is the same curve as   v = -c.
        However,   u   and   v   must be of opposite sign wherever   x < 0
        and   u   and   v   must be the same sign wherever   x > 0.]

        [Graph of parabolic coordinate grid]


    1. A thin wire of line density   r   =   (ax + b) e z   is laid out along a circle, centre the origin, radius   r, in the   x-y plane, (where   a   and   b   are any real constants and   r   is any positive real constant).

      1. Show that a requirement that all parts of the wire have positive mass leads to the constraint   b > | a | r.

        [graph of unit circle] r   =   (ax + b) e z
        a, b are real numbers,
         r is a positive real number
        All of the wire is in the plane   z = 0
        On C :
        x = r cos t,     y = r sin t,     z = 0
            with     0 < t < 2p
        Þ     r   =   (ar cos t + b) e 0
        But -1 < cos t < 1
        The minimum value of the line density
        anywhere on C is therefore
        rmin   =   (- | a | r + b)
        Positive mass     Þ   r > 0   everywhere on C
        rmin > 0     Þ     b   >   | a | r


      2. Find the mass of the wire (as a function of   a, b   and   r ).

        vector r  =  < rc, rs, 0 >  ==>
         (d/dt)(vector r) = < -rs, rc, 0 >
        ==> ds/dt = sqrt(r^2 s^2 + r^2 c^2 + 0) = r
        ==> m = (Int over c)(rho)ds/dt dt = 
         [ ar^2 s + brt]| evaluated at 2pi and 0
                (note that   a, b, r   are all constant )
        = (-0 + 2pibr) - (-0 + 0)

        m = 2 pi b r


      3. Show that the centre of mass is at the point ( a*r^2/(2*b),  0,  0 ).

        M = Int. over c) (rho)(r-arrow)ds/dt dt =
         (Int 0 - 2pi)(arc+b) < rc,rs,0 > r dt
        r^2(Int 0 - 2pi) < arc^2 + bc, arcs + bs, 0 > dt
        Evaluating each of the trigonometric integrals,
        (Int)cos^2tdt = pi
        (Int 0 - 2pi)cost dt = [sin t](eval 2pi and 0) = 0
        (Int 0 - 2pi) cost sint dt =  0
        (Int 0 - 2pi) sint dt = 0
        vector M = r^2 < ar pi, 0, 0 > = pi ar^3 iHat
        ==>  < r-arrow > = (vector M)/m
         = ( pi ar^3 < 1, 0, 0 > )/(2 pi br)
         = (ar^2)/(2b)i
        Therefore the centre of mass is at

        ( (ar^2)/(2b), 0, 0 )


      4. With the requirement of part (a) in place, what are the maximum and minimum possible distances of the centre of mass from the origin?


        b > |a| r ==> 
         | (ar^2)/(2b) | < |a|r^2 /(2|a|r) = r/2
        Therefore the maximum distance of the centre of mass from O is   r/2
        [This happens when   rho goes to 0   at one of the x-intercepts.]
        When   a = 0   the centre of mass is at the origin.   Therefore

        dmin = 0, dmax = r/2


      5. The moment of inertia   I   of a body about an axis of rotation   L   is defined by
                DI   =   r 2 Dm
        (where   r   is the distance of the element of mass   Dm   from the axis   L).
        [The moment of inertia is thus a second moment.]
        Find the moment of inertia of the thin wire about the   z-axis.


        delta I = r^2 delta m = r^2 rho delta s 
         = r^2 rho (delta s)/(delta t) delta t
        But   r   is a constant   ==> I = r^2 (Int on C) rho ds/dt dt 
         = r^2 m = 2 pi b r^3
        (from part (b)).   Therefore

          I   =   2pbr3  


      6. The kinetic energy   E   of a rigid body rotating with angular velocity   w   about the axis   L   is given by E   =   ½ Iw   2.   Find the kinetic energy of the thin wire when it rotates about the   z-axis with an angular velocity   w.


        E = Iw^2 / 2 E = pi b r^3 w^2

        [Note that this result is consistent with the kinetic energy formula
        K.E. = (1/2) m v^2


      7. The angular momentum (or moment of momentum) of a rigid body rotating with angular velocity   w   about the axis   L   is   Iw.   [In the absence of any friction or externally applied torque, the angular momentum is conserved.]   Find the angular momentum of the thin wire when it rotates about the   z-axis with an angular velocity   w.


        I w = 2 pi b r^3 w

        The angular momentum may also be calculated from
        L = mv r


    1. Find the mass and centre of mass of a thin wire that is stretched along a straight line between the origin and the point (6, 6, 6), given that the line density at   (x, y, z)   is   (x + y + z ) / 100   kg m-1.


      [graph of line through origin and (6,6,6)] Clearly a vector parallel to   Vector OQ   is   v = < 1, 1, 1 >
      The line passes through the origin O.
      Therefore the equation of the line segment is:
      vector r = t < 1, 1, 1 > ;
     (0 <= t <= 6)
      ==> (d/dt)r = < 1, 1, 1 > for all t
      Line density:
      rho( x, y, z) = 3t/100
      ds/dt = | dr/dt | = sqrt{3}
      mass = (Int over C) rho ds = 27sqrt{3}/50

      mass = 27sqrt{3}/50 = 0.935 kg

      Both the line density   r   and the geometry of the wire are completely symmetric with respect to any pair of x, y, z
      xBar = yBar = zBar
      Therefore one need take moments about any one coordinate plane only.

      OR   (ignoring this shortcut),

      Finding the vector moment M:
      vector M = (Int over C) (rho)(r)(ds/dt)dt
      Moment = 54sqrt{3}/25 < 1,1,1 >
      The position vector of the centre of mass is
      vector r = M/m = 4 < 1,1,1 >

      centre of mass at (4,4,4)


    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 06 and most recently modified 2008 12 27 by Dr. G.H. George