Appendix a Curvilinear Coordinates
Total Page:16
File Type:pdf, Size:1020Kb
Appendix A Curvilinear coordinates A.1 Lam´e coefficients Consider set of equations ξi = ξi(x1, x2, x3) , i = 1, 2, 3 where ξ1, ξ2, ξ3 independent, single-valued and continuous (x1, x2, x3) : coordinates of point P in Cartesian system with radius-vector x (ξ1, ξ2, ξ3) : coordinates of point P in curvilinear system e1, e2, e3: base unit vectors in Cartesian coordinates eξ1 , eξ2 , eξ3 : base unit vectors in curvilinear system we restrict ourselves to orthogonal curvilinear coordinates systems e e = δ ξi · ξj ij Consider x = x(xi(ξj)) and take the total differential (summation over double appearing indices understood) ∂x dx = dξj ∂ξj ∂x partial derivative keeping ξ2, ξ3 fixed: vector tangential to coordinate line ξ1 ∂ξ1 ∂x = h1 eξ1 −→ ∂ξ1 in general dx = h1 dξ1 eξ1 + h2 dξ2 eξ2 + h3 dξ3 eξ3 hi Lam´e coefficients Those coefficients (might) depend on the coordinates ξi hi = hi(ξ1, ξ2, ξ3) 157 hi dξi (no summation) components of length element along eξi Without restriction to orthogonal systems the length element ds = dx is defined as | | 2 ∂xj ∂xj ds = dx dx = dxj dxj = dξl dξm = glm dξl dξm · ∂ξl ∂ξm Coefficients ∂xj ∂xj glm = ∂ξl ∂ξm are called the elements of the metric tensor For an orthogonal curvilinear system the base vectors eξ1 , eξ2 , eξ3 are mutually perpendic- ular and form a right-handed system We get the length element squared 2 2 2 2 2 2 2 ds = h1 (dξ1) + h2 (dξ2) + h3 (dξ3) the Lam´e coefficients 2 2 2 g11 = h1 , g22 = h2 , g33 = h3 , g = 0 for i = k ik 6 the volume element (elementary rectangular box) 3 d x dV = h1h2h3 dξ1dξ2dξ3 ≡ 158 A.2 Some special orthogonal coordinates Cartesian coordinates ( < x1 < , < x2 < , < x3 < ) • −∞ ∞ −∞ ∞ −∞ ∞ h1 = h2 = h3 = 1 2 2 2 2 ds = dx1 + dx2 + dx3 3 d x = dx1 dx2 dx3 Cylindrical coordinates (0 ρ < , 0 ϕ < 2π , < z < ) • ≤ ∞ ≤ −∞ ∞ x1 = ρ cos ϕ , x2 = ρ sin ϕ , x3 = z hρ = hz = 1 , hϕ = ρ ds2 = dρ2 + ρ2dϕ2 + dz2 d3x = ρ dρ dϕ dz Spherical coordinates (0 r < , 0 θ π , 0 ϕ < 2π) • ≤ ∞ ≤ ≤ ≤ x1 = r sin θ cos ϕ , x2 = r sin θ sin ϕ , x3 = r cos θ hr = 1 , hθ = r , hϕ = r sin θ ds2 = dr2 + r2dθ2 + r2 sin θ2 dϕ2 d3x = r2 dr sin θ dθ dϕ = r2 dr d cos θ dϕ r2dr dΩ ≡ Parabolic cylindrical coordinates • parametrization in Mathematica (0 u < , 0 v < , < z < ) ≤ ∞ ≤ ∞ −∞ ∞ 1 2 2 x1 = (u v ) , x2 = u v , x3 = z 2 − 2 2 hu = hv = √u + v , hz = 1 ds2 = (u2 + v2)(du2 + dv2) + dz2 d3x = (u2 + v2) du dv dz parametrization of Arfken (0 ξ < , 0 η < , < z < ) ≤ ∞ ≤ ∞ −∞ ∞ x1 x2 , u η , v ξ ↔ → → Parabolic coordinates • parametrization of Arfken, Landau/Lifshitz (0 ξ < , 0 η < , 0 ϕ < 2π) ≤ ∞ ≤ ∞ ≤ 1 x1 = ξ η cos ϕ , x2 = ξ η sin ϕ , x3 = (ξ η) 2 − p p 159 1 ξ + η 1 ξ + η hξ = , hη = , hϕ = ξ η 2s ξ 2s η 1 ξ + η 1 ξ + η p ds2 = dξ2 + dη2 + ξη dϕ2 4 ξ 4 η 1 d3x = (ξ + η) dξ dη dϕ 4 parametrization in Mathematica (0 u < , 0 v < , 0 ϕ < 2π) ≤ ∞ ≤ ∞ ≤ 2 2 1 2 2 ξ = u , η = v , x1 = u v cos ϕ , x2 = u v sin ϕ , x3 = (u v ) 2 − 2 2 hu = hv = √u + v , hϕ = uv ds2 = (u2 + v2)(du2 + dv2) + uv dϕ2 d3x = uv(u2 + v2) du dv dϕ 160 A.3 Vector operations in orthogonal coordinates Gradient The length element along coordinate line ξ1 is h1 dξ1 1 ∂ψ eξ1 gradψ = eξ1 ∇ψ = −→ · · h1 ∂ξ1 ∇ 1 ∂ψ 1 ∂ψ 1 ∂ψ ψ = eξ1 + eξ2 + eξ3 h1 ∂ξ1 h2 ∂ξ2 h3 ∂ξ3 form of the nabla operator ⇒ 3 1 ∂ ∇ = e ξi h ∂ξ i=1 i i X Divergence and Laplacian Use the divergence definition 1 ∇ A = lim A n da · ∆V →0 ∆V · IS with A = A1(ξ1, ξ2, ξ3) eξ1 + A2(ξ1, ξ2, ξ2) eξ1 + A3(ξ1, ξ2, ξ3) eξ3 and choose as volume ∆V an elementary rectangular box of sides h1∆ξ1h2∆ξ2h3∆ξ3 Analogously to the derivation in Cartesian coordinates (compare Chapter 1.2) we calculate the net outward flux in direction of coordinate line ξ1 divided by the volume: lim A1(ξ1 + ∆ξ1, ξ2, ξ3) h2(ξ1 + ∆ξ1, ξ2, ξ3)∆ξ2 h3(ξ1 + ∆ξ1, ξ2, ξ3)∆ξ3 ∆V →0 − A1(ξ1, ξ2, ξ3) h2(ξ1, ξ2, ξ3)∆ξ2 h3(ξ1, ξ2, ξ3)∆ξ3 /(h1∆ξ1h2∆ξ2h3∆ξ3) 1 ∂ (h2h3A1) → h1h2h3 ∂ξ1 we get for the divergence ⇒ 1 ∂ ∂ ∂ ∇ A = (h2h3A1) + (h3h1A2) + (h1h2A3) · h1h2h3 ∂ξ1 ∂ξ2 ∂ξ3 The Laplacian operator follows (∇2 = 2 = ∆): ∇ 1 ∂ h2h3 ∂ ∂ h3h1 ∂ ∂ h1h2 ∂ ∆ = + + h1h2h3 ∂ξ1 h1 ∂ξ1 ∂ξ2 h2 ∂ξ2 ∂ξ3 h3 ∂ξ3 161 Curl Using a similar derivation as in Cartesian coordinates we get for the curl of vector A 1 ∂ ∂ A = (h3A3) (h2A2) eξ1 + cyclic permutations ∇ × h2h3 ∂ξ2 − ∂ξ3 Representation as determinant: h1 eξ1 h2 eξ2 h3 eξ3 ∇ 1 ∂ ∂ ∂ A = × h1h2h3 ∂ξ1 ∂ξ2 ∂ξ3 h1A1 h2A2 h3A3 162 A.4 Some explicit forms of vector operations Cartesian coordinates (x1, x2, x3) with base vectors e1, e2, e3 ∂ψ ∂ψ ∂ψ ∇ψ = e1 + e2 + e3 ∂x1 ∂x2 ∂x3 ∂A1 ∂A2 ∂A3 ∇ A = + + · ∂x1 ∂x2 ∂x3 ∂A3 ∂A2 ∂A1 ∂A3 ∂A2 ∂A1 ∇ A = e1 + e2 + e3 × ∂x2 − ∂x3 ∂x3 − ∂x1 ∂x1 − ∂x2 ∂2ψ ∂2ψ ∂2ψ ∆ ψ = 2 + 2 + 2 ∂x1 ∂x2 ∂x3 Cylindrical coordinates (ρ, ϕ, z) with base vectors eρ, eϕ, ez ∂ψ 1 ∂ψ ∂ψ ∇ψ = e + e + e ∂ρ ρ ρ ∂ϕ ϕ ∂z z 1 ∂ 1 ∂A ∂A ∇ A = (ρ A ) + ϕ + z · ρ ∂ρ ρ ρ ∂ϕ ∂z 1 ∂A ∂A ∂A ∂A 1 ∂ ∂A ∇ A = z ϕ e + ρ z e + (ρ A ) ρ e × ρ ∂ϕ − ∂z ρ ∂z − ∂ρ ϕ ρ ∂ρ ϕ − ∂ϕ z 1 ∂ ∂ψ 1 ∂2ψ ∂2ψ ∆ ψ = ρ + + ρ ∂ρ ∂ρ ρ2 ∂ϕ2 ∂z2 Spherical coordinates (r, θ, ϕ) with base vectors er, eθ, eϕ ∂ψ 1 ∂ψ 1 ∂ψ ∇ψ = e + e + e ∂r r r ∂θ θ r sin θ ∂ϕ ϕ 1 ∂ 1 ∂ 1 ∂A ∇ A = (r2A ) + (sin θ A ) + ϕ · r2 ∂r r r sin θ ∂θ θ r sin θ ∂ϕ 1 ∂ ∂A ∇ A = (sin θ A ) θ e × r sin θ ∂θ ϕ − ∂ϕ r 1 ∂A 1 ∂ 1 ∂ ∂A + r (r A ) e + (r A ) r e r sin θ ∂ϕ − r ∂r ϕ θ r ∂r θ − ∂θ ϕ 1 ∂ ∂ψ 1 ∂ ∂ψ 1 ∂2ψ ∆ ψ = r2 + sin θ + r2 ∂r ∂r r2 sin θ ∂θ ∂θ r2 sin2 θ ∂ϕ2 Note 1 ∂ ∂ψ 1 ∂2 r2 (r ψ) r2 ∂r ∂r ≡ r ∂r2 163 A.5 Relation between unit base vectors and their time derivatives Relation between the local unit base vectors of cylindrical coordinates eρ(t), eϕ(t), ez and their time derivatives and the global constant unit base vectors of Cartesian coordinates e1, e2, e3 e = cos ϕ e1 + sin ϕ e2 e1 = cos ϕ e sin ϕ e ρ ρ − ϕ e = sin ϕ e1 + cos ϕ e2 e2 = sin ϕ e + cos ϕ e ϕ − ρ ϕ ez = e3 e3 = ez e˙ = sin ϕ ϕ˙ e1 + cos ϕ ϕ˙ e2 = ϕ˙ e e˙ 1 = 0 ρ − ϕ e˙ = ϕ˙ e e˙ 2 = 0 ϕ − ρ e˙ z = 0 e˙ 3 = 0 Same for the local unit base vectors of spherical coordinates er(t), eθ(t), eϕ(t) and their time derivatives er = sin θ cos ϕ e1 + sin θ sin ϕ e2 + cos θ e3 e = cos θ cos ϕ e1 + cos θ sin ϕ e2 sin θ e3 θ − e = sin ϕ e1 + cos ϕ e2 ϕ − e1 = sin θ cos ϕ e + cos θ cos ϕ e sin ϕ e r θ − ϕ e2 = sin θ sin ϕ er + cos θ sin ϕ eθ + cos ϕ eϕ e3 = cos θ e sin θ e r − θ e˙ r = θ˙ eθ + sin θ ϕ˙ eϕ e˙ = θ˙ e + cos θ ϕ˙ e θ − r ϕ e˙ = ϕ˙ (sin θ e + cos θ e ) ϕ − r θ Relation between the local unit base vectors eρ(t), eϕ(t), ez and er(t), eθ(t), eϕ(t) eρ = sin θ er + cos θ eθ er = sin θ eρ + cos θ ez eϕ = eϕ eθ = cos θ eρ sin θ ez e = cos θ e sin θ e e = e − z r − θ ϕ ϕ 164.