4.5 Curve integral of the second kind

Definition 4.5.1 Directional curve (path). Given x:[a,b]→Ω (a kind of motion), x(t) represents the position (t is the time), γ={x(t)|t∈[a,b]}.

Note x(t) is not a parametric representative of γ! Curve γ is a static geometric object, while path x(t) is a dynamic motion (meaning it could turn back along γ).

Example 4.5.1 Work done by force. F(x):Ω→Rm (continuous) is a vector field where x represents the position. Work

      ∫          ∫  b          ′        ∫           x ′(t)           ′        ∫
W  =    F ⋅ dx =    ⟨F (x(t)),x (t)⟩dt =    ∖lef t⟨F, ---′---∖right ⟩∥x (t)∥dt =   ⟨F, T ⟩dl
       γ           a                     γ         ∥x (t)∥                    γ

where T(x) represents a forward-direction unit tangent vector field of γ.

Flow velocity field F (continuous). Assuming T(x) represents a forward-direction unit tangent vector field of closed curve γ, then the circulation is

∫
  ⟨F(x ),T(x )⟩dl
 γ

Assuming

F (x) =,     x(t) =,a ≤  t ≤ b

Then we have

∫             ∫          ∫  b                                        ∫
                                      ′                     ′
  γ⟨F, T ⟩dl =  γ F ⋅ dx =  a [F1(x (t))x1(t) + ⋅⋅⋅ + Fm (x(t))x m(t)]dt = γ ω

where ω=F1(x)dx1+⋯+Fm(x)dxm has the form of first-order derivative.

Example 4.5.2 Assuming γ is a simple closed curve in plain (boundary of plain), N(x) represents a forward-direction unit outer-norm vector field of γ, then the flux is

∫

 γ⟨F(x ),N (x )⟩dl

The relation between circulation and flux. Assuming the norm vector k of the plain (direction determined by the forward direction of γ and right-hand rule), then construct F~=k×F where F~ is derived by rotating F by 90∘ counterclockwise. Since N=T×k, naturally we have

∫             ∫                 ∫                ∫
   ⟨F, N ⟩dl =   ⟨F, T × k ⟩dl =   ⟨k × F, T⟩dl =   ⟨F˜, T ⟩dl
  γ            γ                 γ                γ

Assuming

      (x′,y ′)          (y′,− x′)
T =  ---′--′--⇒  N  = ---′--′--,    F  =
     ∥(x,y )∥         ∥(x ,y )∥

then we have

∫             ∫                                   ∫
                b          ′             ′
   ⟨F,N ⟩dl =    [X (x(t))y(t) − Y(x (t))x (t)]dt =   (− Y dx + Xdy  )
  γ            a                                   γ

It also has a form of first-order derivative.

Example 4.5.3 Assuming

γ :

and its direction is determined by right-hand rule circling around +z axis. Seek ∫γzdx+xdy+ydz.

Firstly, seek the parametric equation of γ. Eliminate z by z=−x−y to get x2+y2+(−x−y)2=R2, so

          y        2  3  2   R2
∖left(x + -∖right ) + --y  = ---
          2           4       2

Let (x+y2,32y)=R2(cos⁡θ,sin⁡θ), then z=−x−y=⋯.

  • The range of θ. It is trivial that θ∈[0,2π].
  • Is the increasing direction of θ the same as the direction of γ? Take θ=0, we can get that (x,y)=(R2,0) and x′<0,y′>0. Then we can verify that these directions are the same.

Then on γ, we have ω=(−x−y)dx+xdy+yd(−x−y)=−xdx−2ydx+xdy−ydy=d(−x22−y22+xy)−3ydx∫ω=∫γd(−x22−y22+xy)⏟0 due to the closure of γ−∫γ3ydx=−3∫γydx

Then it’s left as exercise. :)

Example 4.5.4 Assuming γ is a simple closed (Jordan) oriented (counterclockwise, natural positive direction) plain curve of C1, then

∫           ∫          ∫
   − ydx =    xdy  = 1-   xdy − ydx =  A(γ )
  γ          γ       2  γ

where A(γ) represents the area of closed domain bounded with γ.

Assuming (x(t),y(t))=(r(t)cos⁡θ(t),r(t)sin⁡θ(t)),t∈[a,b] where (x(a),y(a))=(x(b),y(b)), then we have

∫        ∫        ∫
                                     b
  xdy +    ydx  =    d(xy) = x(t)y(t)|a = 0
 γ        γ        γ

Notice that 12(xdy−ydx)=12r(t)cos⁡θ(t)d(r(t)sin⁡θ(t))−12r(t)sin⁡θ(t)d(r(t)cos⁡θ(t))=r2(t)θ′(t)∫γ12(xdy−ydx)=12∫abr2(t)θ′(t)dt

So the former 2 equalities are proved.

Assuming θ′(t)≠0, then θ(t) has inverse function t(θ), we have

∫                    ∫                 ∫             ∫   ∫
   1               1   b 2           1   b 2           2π  r(𝜃)
   -(xdy − ydx ) = --   r (t(𝜃))d𝜃 = --   r (𝜃)d𝜃 =            rdrd𝜃 = A (γ)
 γ 2               2  a              2  a             0   0

Theorem is proved.

Definition 4.5.2 F:Ω→Rm is a field with potential if ∃f:Ω→R such that ∀x∈Ω, F(x)=∇f(x). f is called a potential function of F.

Theorem 4.5.3 For any path γ (starting from A, ending at B) of C1 in Ω, we have

∫
   ∇f  ⋅ dx = f (B) − f(A )
 γ

Assuming

F ==  ∇f  =

then we have

                                   ∂f              ∂f
ω = F1 (x )dx1 + ⋅⋅⋅ + Fm (x)dxm =  ---dx1 +  ⋅⋅⋅ +----dxm  =  df
                                   ∂x1            ∂xm

where ω is called the total differential of f. Then ∫γdf=f|AB (Newton-Leibniz formula).

Definition 4.5.4 If the integral of a vector field on any path only relates to the starting point and the destination, i.e. it has nothing to do with the path itself, then the vector field is conservative.

Corollary 4.5.5 Field with potential is conservative.

Theorem 4.5.6 Assuming ω⊆Rm is a connected open set (domain), then any conservative field on Ω has potential function.

Proof Mark F as the conservative field, assuming γ is a path from A0 to B(x). Mark f(x)=∫γF⋅dx, it is to be proven that ∂f∂x1=F1(x).

Notice that f(x+te1)−f(x)=∫0tF(x+se1)⋅dds(x+se1)ds=∫0tF(x+se1)⋅e1ds=∫0tF1(x+se1)ds∂f∂x1(x)=ddtf(x+te1)|t=0=F1(x+te1)|t=0=F1(x)

Repeat this process, we can get that ∇f=F. ◻

Example 4.5.5 Assuming F(x)=−x∥x∥3, seek its potential function.

Notice that

                                           2
ω = − --x1-dx1 − ⋅⋅⋅ − -xm--dxm =  − d(∥x∥--)=  − d∥x∥-=  d--1--
      ∥x ∥3            ∥x ∥3          2∥x∥3       ∥x∥2     ∥x ∥2

hence f(x)=1∥x∥+C.

Example 4.5.6 Seek

∫
     y                y
   (e  + sinx )dx + (xe −  cosy)dy
  γ

where γ is part of (x−π)2+y2=π2, circling from (0,0) to (π,π) in counterclockwise direction.

Notice that

ω =  d(− cosx − sin y + xey)

hence

∫
   ω = (−  cos x − sin y + xey)|(π,π)=  ⋅⋅⋅
  γ                           (0,0)

Example 4.5.7 Seek

∫
  (x2 − yz)dx + (y2 − zx )dy +  (z2 − xy )dz
 γ

where γ:x=acos⁡t,y=asin⁡t,z=bt, t∈[0,2π], a,b>0.

Notice that

            x3-+-y3-+-z3
ω =  d∖left(     3       − xyz ∖right)

hence

∫
                   x3-+-y3 +-z3                    2π
   ω = ∖left.∖left(     3      −  xyz∖right)∖right |t=0 = ⋅⋅⋅
 γ
Assuming ω=F1(x)dx1+⋯+Fm(x)dxm=df (F=∇f), then Fi=∂f∂i. Assuming f∈C2, we have ∂2f∂xj∂xi=∂2f∂xi∂xj.

Definition 4.5.7 Assuming F is a vector field. If

∂Fi-   ∂Fj-
∂x   = ∂x  ,     ∀1 ≤ i < j ≤ m
   j      i

then F is irrotational.

Corollary 4.5.8 Conservative field is irrotational.

Example 4.5.8 Irrotational may not be conservative. Assuming F=(−yx2+y2,xx2+y2)T is defined on R2∖{0,0}, it is trivial that it is irrotational since

-∂-      --−-y---          ∂--      --x-----
∂y ∖left(x2 + y2∖right ) = ∂x∖lef t( x2 + y2∖right)

Yet after transformation (x,y)→(r,θ),

∫           ∫ 2π
       ω =      d𝜃 =  2π ⁄= 0
 ∂B(0,a)      0

In fact, the integral reflects the geometric property of γ – how many turns does γ circle around the origin by.

Assuming θ=arctan⁡yx where −π2<θ<π2, then

         y-   − ydx-+-xdy-                        y-
d arctan x =    x2 + y2   =  ω,     ∇ ∖lef t(arctan x ∖right) = F

Since we can’t define a singular-value continuous function θ:R2→R, so the vector field above is not conservative! If the domain of definition Ω:R2∖{(x,0)|x≥0}, then the vector field is conservative. It relates to the topological structure of the domain.

Plain vector field.

Definition 4.5.9 Assuming F=(X,Y)T is a vector field of C1, mark

        ∂Y     ∂X               ∂X    ∂Y
rotF  = ----−  ---,     div F =  ----+ ----
        ∂x     ∂y               ∂x     ∂x

respectively as the rotation and divergence of F.

Theorem 4.5.10 (Green formula in physical form) Assuming Ω⊆R2 is a bounded closed domain, ∂Ω is piecewisely of C1, F:Ω→R2 is a vector field of C1, then

1.

∫∂ΩF⋅dx=∫ΩrotFdxdy.

2.

∫∂Ω⟨F,n⟩dl=∫ΩdivFdxdy.

where the forward direction of ∂Ω is the natual positive direction.

Mark ∇=(∂∂x,∂∂y,∂∂z), we could mark
rot F = ∇  × F,     div F = ∇  ⋅ F

Notice that

divF =  tr ∂(X,-Y-)
           ∂(x,y)

indicating that the divergence has nothing to do with the selection of the coordinate system, since the similarity transformation of a matrix doesn’t change its trace.

Definition 4.5.11 Wedge (exterior) product. Define a formal product dx∧dy satisfying

  • Binary linearity.

    (adx  +  bdx ) ∧ dy = adx  ∧ dy + bdx  ∧ dy
     1      2             1           2
  • Anti-symmetry.

    dx ∧ dy =  − dy ∧ dx

It is trivial that dx∧dx=0. Actually, dx∧dy represents a oriented area, yet dxdy represents a non-oriented area. When moving forward along ∂Ω, the outside of Ω is on the right, i.e. the forward direction of ∂Ω is the natual positive direction, we have dx∧dy=dxdy.

Definition 4.5.12 Exterior differential. Assuming a form of first-order derivative ω=Xdx+Ydy, define

dω = dX  ∧ dx +  dY ∧ dy =  ∖lef t(∂Y--− ∂X--∖right)dx ∧ dy
                                  ∂x    ∂y

Then the third form of Green formula is

Theorem 4.5.13 (Green formula in mathematical form) ∫∂Ωω=∫Ωdω.

Notice that
⟨F, n⟩dl = − Y dx + XdY   = ˜ω

so we have

            ∂X    ∂Y
d˜ω = ∖lef t(∂x--+ -∂y-∖right)dx ∧ dy

verifying the second form of Green formula.

Note The physical meaning of divergence. Assuming F is a vector field of C1, Ω(ϵ)=B(P0,ϵ), then we have

                   ∫                       ∫
divF (P0) = lim  -1--  div Fdxdy  = lim -1--   ⟨F, n⟩dl
            𝜖→0 π𝜖2  Ω             𝜖→0 π𝜖2  ∂Ω

Note

Corollary 4.5.14 If Ω⊆R2 is a single-connected domain, i.e. any continuous closed loop in Ω could become a single point after continuous transformation in Ω, then any irrotational vector field of C1 on Ω is conservative.

Proof Assuming D⊆Ω, γ=∂D is a closed loop, then we have

∫      ∫           ∫

   ω =    F ⋅ dx =    rotFdxdy  =  0
 γ      γ           D

meaning F is conservative. ◻

Example 4.5.9 Prove

∫           ∫        1 ∫
   − ydx =    xdy  = --   xdy − ydx =  A(γ )
  γ          γ       2  γ

with Green formula.

Assuming ω1=−ydx, then we have

∫      ∫         ∫                 ∫             ∫
  ω  =     dω  =    (− 1)dy ∧ dx =    dx ∧ dy =     dxdy =  A(γ )
 γ  1    D   1    D                 D             D

Then it’s left as exercise. :)

Example 4.5.10 Assuming γ:x2a2+y2b2=1(y≥0) from A(a,0) to B(−a,0), seek

∫
  (1 + yex)dx + (x + ex)dy
 γ

Let γ′:{(x,0)|−a≤x≤a} from B to A, γ1=γ∪γ′ is a closed loop. Notice that

    ∫          x   ∫                x B   ∫         ∫                π
I =   d (x + ye ) +   xdy =  (x + ye )A +    xdy −     xdy =  − 2a + -ab
     γ               γ                     γ1       ◟γ′◝◜--◞          2
                                                    γ′:dy=0