Del in cylindrical and spherical coordinates

From testwiki
Jump to navigation Jump to search

Template:Short description This is a list of some vector calculus formulae for working with common curvilinear coordinate systems.

Notes

  • This article uses the standard notation ISO 80000-2, which supersedes ISO 31-11, for spherical coordinates (other sources may reverse the definitions of θ and φ):
    • The polar angle is denoted by θ∈[0,π]: it is the angle between the z-axis and the radial vector connecting the origin to the point in question.
    • The azimuthal angle is denoted by φ∈[0,2π]: it is the angle between the x-axis and the projection of the radial vector onto the xy-plane.
  • The function Template:Nowrap can be used instead of the mathematical function Template:Nowrap owing to its domain and image. The classical arctan function has an image of Template:Nowrap, whereas atan2 is defined to have an image of Template:Nowrap.

Coordinate conversions

Conversion between Cartesian, cylindrical, and spherical coordinates[1]
From
Cartesian Cylindrical Spherical
To Cartesian x=xy=yz=z x=ρcos⁡φy=ρsin⁡φz=z x=rsin⁡θcos⁡φy=rsin⁡θsin⁡φz=rcos⁡θ
Cylindrical ρ=x2+y2φ=arctan⁡(yx)z=z ρ=ρφ=φz=z ρ=rsin⁡θφ=φz=rcos⁡θ
Spherical r=x2+y2+z2θ=arccos⁡(zx2+y2+z2)φ=arctan⁡(yx) r=ρ2+z2θ=arctan⁡(ρz)φ=φ r=rθ=θφ=φ

Note that the operation arctan⁡(AB) must be interpreted as the two-argument inverse tangent, atan2.

Unit vector conversions

Conversion between unit vectors in Cartesian, cylindrical, and spherical coordinate systems in terms of destination coordinates[1]
Cartesian Cylindrical Spherical
Cartesian 𝐱̂=𝐱̂𝐲̂=𝐲̂𝐳̂=𝐳̂ 𝐱̂=cos⁡φ𝝆̂−sin⁡φ𝝋̂𝐲̂=sin⁡φ𝝆̂+cos⁡φ𝝋̂𝐳̂=𝐳̂ 𝐱̂=sin⁡θcos⁡φ𝐫̂+cos⁡θcos⁡φ𝜽̂−sin⁡φ𝝋̂𝐲̂=sin⁡θsin⁡φ𝐫̂+cos⁡θsin⁡φ𝜽̂+cos⁡φ𝝋̂𝐳̂=cos⁡θ𝐫̂−sin⁡θ𝜽̂
Cylindrical 𝝆̂=x𝐱̂+y𝐲̂x2+y2𝝋̂=−y𝐱̂+x𝐲̂x2+y2𝐳̂=𝐳̂ 𝝆̂=𝝆̂𝝋̂=𝝋̂𝐳̂=𝐳̂ 𝝆̂=sin⁡θ𝐫̂+cos⁡θ𝜽̂𝝋̂=𝝋̂𝐳̂=cos⁡θ𝐫̂−sin⁡θ𝜽̂
Spherical 𝐫̂=x𝐱̂+y𝐲̂+z𝐳̂x2+y2+z2𝜽̂=(x𝐱̂+y𝐲̂)z−(x2+y2)𝐳̂x2+y2+z2x2+y2𝝋̂=−y𝐱̂+x𝐲̂x2+y2 𝐫̂=ρ𝝆̂+z𝐳̂ρ2+z2𝜽̂=z𝝆̂−ρ𝐳̂ρ2+z2𝝋̂=𝝋̂ 𝐫̂=𝐫̂𝜽̂=𝜽̂𝝋̂=𝝋̂
Conversion between unit vectors in Cartesian, cylindrical, and spherical coordinate systems in terms of source coordinates
Cartesian Cylindrical Spherical
Cartesian 𝐱̂=𝐱̂𝐲̂=𝐲̂𝐳̂=𝐳̂ 𝐱̂=x𝝆̂−y𝝋̂x2+y2𝐲̂=y𝝆̂+x𝝋̂x2+y2𝐳̂=𝐳̂ 𝐱̂=x(x2+y2𝐫̂+z𝜽̂)−yx2+y2+z2𝝋̂x2+y2x2+y2+z2𝐲̂=y(x2+y2𝐫̂+z𝜽̂)+xx2+y2+z2𝝋̂x2+y2x2+y2+z2𝐳̂=z𝐫̂−x2+y2𝜽̂x2+y2+z2
Cylindrical 𝝆̂=cos⁡φ𝐱̂+sin⁡φ𝐲̂𝝋̂=−sin⁡φ𝐱̂+cos⁡φ𝐲̂𝐳̂=𝐳̂ 𝝆̂=𝝆̂𝝋̂=𝝋̂𝐳̂=𝐳̂ 𝝆̂=ρ𝐫̂+z𝜽̂ρ2+z2𝝋̂=𝝋̂𝐳̂=z𝐫̂−ρ𝜽̂ρ2+z2
Spherical 𝐫̂=sin⁡θ(cos⁡φ𝐱̂+sin⁡φ𝐲̂)+cos⁡θ𝐳̂𝜽̂=cos⁡θ(cos⁡φ𝐱̂+sin⁡φ𝐲̂)−sin⁡θ𝐳̂𝝋̂=−sin⁡φ𝐱̂+cos⁡φ𝐲̂ 𝐫̂=sin⁡θ𝝆̂+cos⁡θ𝐳̂𝜽̂=cos⁡θ𝝆̂−sin⁡θ𝐳̂𝝋̂=𝝋̂ 𝐫̂=𝐫̂𝜽̂=𝜽̂𝝋̂=𝝋̂

Del formula

Template:Sticky table start

Table with the del operator in cartesian, cylindrical and spherical coordinates
Operation Cartesian coordinates Template:Math Cylindrical coordinates Template:Math Spherical coordinates Template:Math,
where Template:Math is the polar angle and Template:Math is the azimuthal angleTemplate:Ref
Vector field Template:Math Ax𝐱̂+Ay𝐲̂+Az𝐳̂ Aρ𝝆̂+Aφ𝝋̂+Az𝐳̂ Ar𝐫̂+Aθ𝜽̂+Aφ𝝋̂
Gradient Template:Math[1] ∂f∂x𝐱̂+∂f∂y𝐲̂+∂f∂z𝐳̂ ∂f∂ρ𝝆̂+1ρ∂f∂φ𝝋̂+∂f∂z𝐳̂ ∂f∂r𝐫̂+1r∂f∂θ𝜽̂+1rsin⁡θ∂f∂φ𝝋̂
Divergence Template:Math[1] ∂Ax∂x+∂Ay∂y+∂Az∂z 1ρ∂(ρAρ)∂ρ+1ρ∂Aφ∂φ+∂Az∂z 1r2∂(r2Ar)∂r+1rsin⁡θ∂∂θ(Aθsin⁡θ)+1rsin⁡θ∂Aφ∂φ
Curl Template:Math[1] (∂Az∂y−∂Ay∂z)𝐱̂+(∂Ax∂z−∂Az∂x)𝐲̂+(∂Ay∂x−∂Ax∂y)𝐳̂ (1ρ∂Az∂φ−∂Aφ∂z)𝝆̂+(∂Aρ∂z−∂Az∂ρ)𝝋̂+1ρ(∂(ρAφ)∂ρ−∂Aρ∂φ)𝐳̂ 1rsin⁡θ(∂∂θ(Aφsin⁡θ)−∂Aθ∂φ)𝐫̂+1r(1sin⁡θ∂Ar∂φ−∂∂r(rAφ))𝜽̂+1r(∂∂r(rAθ)−∂Ar∂θ)𝝋̂
Laplace operator Template:Math[1] ∂2f∂x2+∂2f∂y2+∂2f∂z2 1ρ∂∂ρ(ρ∂f∂ρ)+1ρ2∂2f∂φ2+∂2f∂z2 1r2∂∂r(r2∂f∂r)+1r2sin⁡θ∂∂θ(sin⁡θ∂f∂θ)+1r2sin2θ∂2f∂φ2
Vector gradient Template:MathTemplate:Ref ∂Ax∂x𝐱̂⊗𝐱̂+∂Ax∂y𝐱̂⊗𝐲̂+∂Ax∂z𝐱̂⊗𝐳̂+∂Ay∂x𝐲̂⊗𝐱̂+∂Ay∂y𝐲̂⊗𝐲̂+∂Ay∂z𝐲̂⊗𝐳̂+∂Az∂x𝐳̂⊗𝐱̂+∂Az∂y𝐳̂⊗𝐲̂+∂Az∂z𝐳̂⊗𝐳̂ ∂Aρ∂ρ𝝆̂⊗𝝆̂+(1ρ∂Aρ∂φ−Aφρ)𝝆̂⊗𝝋̂+∂Aρ∂z𝝆̂⊗𝐳̂+∂Aφ∂ρ𝝋̂⊗𝝆̂+(1ρ∂Aφ∂φ+Aρρ)𝝋̂⊗𝝋̂+∂Aφ∂z𝝋̂⊗𝐳̂+∂Az∂ρ𝐳̂⊗𝝆̂+1ρ∂Az∂φ𝐳̂⊗𝝋̂+∂Az∂z𝐳̂⊗𝐳̂ ∂Ar∂r𝐫̂⊗𝐫̂+(1r∂Ar∂θ−Aθr)𝐫̂⊗𝜽̂+(1rsin⁡θ∂Ar∂φ−Aφr)𝐫̂⊗𝝋̂+∂Aθ∂r𝜽̂⊗𝐫̂+(1r∂Aθ∂θ+Arr)𝜽̂⊗𝜽̂+(1rsin⁡θ∂Aθ∂φ−cot⁡θAφr)𝜽̂⊗𝝋̂+∂Aφ∂r𝝋̂⊗𝐫̂+1r∂Aφ∂θ𝝋̂⊗𝜽̂+(1rsin⁡θ∂Aφ∂φ+cot⁡θAθr+Arr)𝝋̂⊗𝝋̂
Vector Laplacian Template:Math[2] ∇2Ax𝐱̂+∇2Ay𝐲̂+∇2Az𝐳̂

(∇2Aρ−Aρρ2−2ρ2∂Aφ∂φ)𝝆̂+(∇2Aφ−Aφρ2+2ρ2∂Aρ∂φ)𝝋̂+∇2Az𝐳̂

(∇2Ar−2Arr2−2r2sin⁡θ∂(Aθsin⁡θ)∂θ−2r2sin⁡θ∂Aφ∂φ)𝐫̂+(∇2Aθ−Aθr2sin2θ+2r2∂Ar∂θ−2cos⁡θr2sin2θ∂Aφ∂φ)𝜽̂+(∇2Aφ−Aφr2sin2θ+2r2sin⁡θ∂Ar∂φ+2cos⁡θr2sin2θ∂Aθ∂φ)𝝋̂

Directional derivative Template:Math[3] 𝐀⋅∇Bx𝐱̂+𝐀⋅∇By𝐲̂+𝐀⋅∇Bz𝐳̂ (Aρ∂Bρ∂ρ+Aφρ∂Bρ∂φ+Az∂Bρ∂z−AφBφρ)𝝆̂+(Aρ∂Bφ∂ρ+Aφρ∂Bφ∂φ+Az∂Bφ∂z+AφBρρ)𝝋̂+(Aρ∂Bz∂ρ+Aφρ∂Bz∂φ+Az∂Bz∂z)𝐳̂

(Ar∂Br∂r+Aθr∂Br∂θ+Aφrsin⁡θ∂Br∂φ−AθBθ+AφBφr)𝐫̂+(Ar∂Bθ∂r+Aθr∂Bθ∂θ+Aφrsin⁡θ∂Bθ∂φ+AθBrr−AφBφcot⁡θr)𝜽̂+(Ar∂Bφ∂r+Aθr∂Bφ∂θ+Aφrsin⁡θ∂Bφ∂φ+AφBrr+AφBθcot⁡θr)𝝋̂

Tensor divergence Template:MathTemplate:Ref

(∂Txx∂x+∂Tyx∂y+∂Tzx∂z)𝐱̂+(∂Txy∂x+∂Tyy∂y+∂Tzy∂z)𝐲̂+(∂Txz∂x+∂Tyz∂y+∂Tzz∂z)𝐳̂

[∂Tρρ∂ρ+1ρ∂Tφρ∂φ+∂Tzρ∂z+1ρ(Tρρ−Tφφ)]𝝆̂+[∂Tρφ∂ρ+1ρ∂Tφφ∂φ+∂Tzφ∂z+1ρ(Tρφ+Tφρ)]𝝋̂+[∂Tρz∂ρ+1ρ∂Tφz∂φ+∂Tzz∂z+Tρzρ]𝐳̂

[∂Trr∂r+2Trrr+1r∂Tθr∂θ+cot⁡θrTθr+1rsin⁡θ∂Tφr∂φ−1r(Tθθ+Tφφ)]𝐫̂+[∂Trθ∂r+2Trθr+1r∂Tθθ∂θ+cot⁡θrTθθ+1rsin⁡θ∂Tφθ∂φ+Tθrr−cot⁡θrTφφ]𝜽̂+[∂Trφ∂r+2Trφr+1r∂Tθφ∂θ+1rsin⁡θ∂Tφφ∂φ+Tφrr+cot⁡θr(Tθφ+Tφθ)]𝝋̂

Template:Sticky table end Template:Unbulleted list

Differential elements

Operation Cartesian coordinates
Template:Math
Cylindrical coordinates
Template:Math
Spherical coordinatesTemplate:Ref
Template:Math
Differential displacement Template:Math[1] dx𝐱̂+dy𝐲̂+dz𝐳̂ dρ𝝆̂+ρdφ𝝋̂+dz𝐳̂ dr𝐫̂+rdθ𝜽̂+rsin⁡θdφ𝝋̂
Differential normal area Template:Math dydz𝐱̂+dxdz𝐲̂+dxdy𝐳̂ ρdφdz𝝆̂+dρdz𝝋̂+ρdρdφ𝐳̂ r2sin⁡θdθdφ𝐫̂+rsin⁡θdrdφ𝜽̂+rdrdθ𝝋̂
Differential volume Template:Math[1] dxdydz ρdρdφdz r2sin⁡θdrdθdφ

Calculation rules

  1. div⁡grad⁡f≡∇⋅∇f≡∇2f
  2. curl⁡grad⁡f≡∇×∇f=𝟎
  3. div⁡curl⁡𝐀≡∇⋅(∇×𝐀)=0
  4. curl⁡curl⁡𝐀≡∇×(∇×𝐀)=∇(∇⋅𝐀)−∇2𝐀 (Lagrange's formula for del)
  5. ∇2(fg)=f∇2g+2∇f⋅∇g+g∇2f
  6. ∇2(𝐏⋅𝐐)=𝐐⋅∇2𝐏−𝐏⋅∇2𝐐+2∇⋅[(𝐏⋅∇)𝐐+𝐏×∇×𝐐] (From [4] )

Cartesian derivation

div⁡𝐀=limV→0∬∂V𝐀⋅d𝐒∭VdV=[Ax(x+dx)−Ax(x)]dydz+[Ay(y+dy)−Ay(y)]dxdz+[Az(z+dz)−Az(z)]dxdydxdydz=∂Ax∂x+∂Ay∂y+∂Az∂z

(curl⁡𝐀)x=limS⊥𝐱̂→0∫∂S𝐀⋅dℓ∬SdS=[Az(y+dy)−Az(y)]dz−[Ay(z+dz)−Ay(z)]dydydz=∂Az∂y−∂Ay∂z

The expressions for (curl⁡𝐀)y and (curl⁡𝐀)z are found in the same way.

Cylindrical derivation

div⁡𝐀=limV→0∬∂V𝐀⋅d𝐒∭VdV=[Aρ(ρ+dρ)(ρ+dρ)−Aρ(ρ)ρ]dϕdz+[Aϕ(ϕ+dϕ)−Aϕ(ϕ)]dρdz+[Az(z+dz)−Az(z)]dρ(ρ+dρ/2)dϕρdϕdρdz=1ρ∂(ρAρ)∂ρ+1ρ∂Aϕ∂ϕ+∂Az∂z

(curl⁡𝐀)ρ=limS⊥𝝆̂→0∫∂S𝐀⋅dℓ∬SdS=Aϕ(z)(ρ+dρ)dϕ−Aϕ(z+dz)(ρ+dρ)dϕ+Az(ϕ+dϕ)dz−Az(ϕ)dz(ρ+dρ)dϕdz=−∂Aϕ∂z+1ρ∂Az∂ϕ

(curl⁡𝐀)ϕ=limS⊥𝝓̂→0∫∂S𝐀⋅dℓ∬SdS=Az(ρ)dz−Az(ρ+dρ)dz+Aρ(z+dz)dρ−Aρ(z)dρdρdz=−∂Az∂ρ+∂Aρ∂z

(curl⁡𝐀)z=limS⊥𝒛̂→0∫∂S𝐀⋅dℓ∬SdS=Aρ(ϕ)dρ−Aρ(ϕ+dϕ)dρ+Aϕ(ρ+dρ)(ρ+dρ)dϕ−Aϕ(ρ)ρdϕρdρdϕ=−1ρ∂Aρ∂ϕ+1ρ∂(ρAϕ)∂ρ

curl⁡𝐀=(curl⁡𝐀)ρ𝝆̂+(curl⁡𝐀)ϕ𝝓̂+(curl⁡𝐀)z𝒛̂=(1ρ∂Az∂ϕ−∂Aϕ∂z)𝝆̂+(∂Aρ∂z−∂Az∂ρ)𝝓̂+1ρ(∂(ρAϕ)∂ρ−∂Aρ∂ϕ)𝒛̂

Spherical derivation

div⁡𝐀=limV→0∬∂V𝐀⋅d𝐒∭VdV=[Ar(r+dr)(r+dr)2−Ar(r)r2]sin⁡θdθdϕ+[Aθ(θ+dθ)sin⁡(θ+dθ)−Aθ(θ)sin⁡θ]rdrdϕ+[Aϕ(ϕ+dϕ)−Aϕ(ϕ)]rdrdθdrrdθrsin⁡θdϕ=1r2∂(r2Ar)∂r+1rsin⁡θ∂(Aθsin⁡θ)∂θ+1rsin⁡θ∂Aϕ∂ϕ

(curl⁡𝐀)r=limS⊥𝒓̂→0∫∂S𝐀⋅dℓ∬SdS=Aθ(ϕ)rdθ+Aϕ(θ+dθ)rsin⁡(θ+dθ)dϕ−Aθ(ϕ+dϕ)rdθ−Aϕ(θ)rsin⁡(θ)dϕrdθrsin⁡θdϕ=1rsin⁡θ∂(Aϕsin⁡θ)∂θ−1rsin⁡θ∂Aθ∂ϕ

(curl⁡𝐀)θ=limS⊥𝜽̂→0∫∂S𝐀⋅dℓ∬SdS=Aϕ(r)rsin⁡θdϕ+Ar(ϕ+dϕ)dr−Aϕ(r+dr)(r+dr)sin⁡θdϕ−Ar(ϕ)drdrrsin⁡θdϕ=1rsin⁡θ∂Ar∂ϕ−1r∂(rAϕ)∂r

(curl⁡𝐀)ϕ=limS⊥𝝓̂→0∫∂S𝐀⋅dℓ∬SdS=Ar(θ)dr+Aθ(r+dr)(r+dr)dθ−Ar(θ+dθ)dr−Aθ(r)rdθrdrdθ=1r∂(rAθ)∂r−1r∂Ar∂θ

curl⁡𝐀=(curl⁡𝐀)r𝒓̂+(curl⁡𝐀)θ𝜽̂+(curl⁡𝐀)ϕ𝝓̂=1rsin⁡θ(∂(Aϕsin⁡θ)∂θ−∂Aθ∂ϕ)𝒓̂+1r(1sin⁡θ∂Ar∂ϕ−∂(rAϕ)∂r)𝜽̂+1r(∂(rAθ)∂r−∂Ar∂θ)𝝓̂

Unit vector conversion formula

The unit vector of a coordinate parameter u is defined in such a way that a small positive change in u causes the position vector 𝐫 to change in 𝐮 direction.

Therefore, ∂𝐫∂u=∂s∂u𝐮 where Template:Mvar is the arc length parameter.

For two sets of coordinate systems ui and vj, according to chain rule, d𝐫=∑i∂𝐫∂uidui=∑i∂s∂ui𝐮̂idui=∑j∂s∂vj𝐯̂jdvj=∑j∂s∂vj𝐯̂j∑i∂vj∂uidui=∑i∑j∂s∂vj∂vj∂ui𝐯̂jdui.

Now, we isolate the ith component. For i≠k, let duk=0. Then divide on both sides by dui to get: ∂s∂ui𝐮̂i=∑j∂s∂vj∂vj∂ui𝐯̂j.

See also

References

Template:Reflist