IJM International Journal of Mechanics Research 2325-498X Scientific Research Publishing 10.12677/IJM.2021.102014 IJM-43509 IJM20210200000_32619500.pdf 数学与物理 曲线正交坐标系下 A 与∇ 2 A 的一种推导 A Derivation of A and ∇ 2 A in Curvilinear Orthogonal Coordinates 葛锐 2 1 梦冬 2 1 小艳 2 1 中南林业科技大学理学院,湖南 长沙 null 12 04 2021 10 02 143 152 © Copyright 2014 by authors and Scientific Research Publishing Inc. 2014 This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

电动力学教材一般仅仅在附录部分给出曲线正交坐标系下标量的梯度Ψ、矢量的散度A和旋度×A,以及标量的拉普拉斯算符∇2Ψ的一般表达式。在教学中我们发现,A和∇2A在曲线正交坐标系下的表达式也是很重要的。本文从微分几何的角度推导出三维欧氏空间中曲线正交坐标系下A和∇2A的一般表达式,以及它们在直角坐标系、柱坐标系和球坐标系的具体表达式,供感兴趣的学生和教师参考。 In textbooks of electrodynamics, the formulas of the gradientΨ, divergenceA, curl×A,and Laplace operator ∇2Ψ in curvilinear orthogonal coordinates are often given in the appendix. However, in the teaching process, we find that the expressions ofAand ∇2Ain curvilinear orthogonal coordinates are also very important. In this paper, the general expressions ofAand ∇2Ain curvilinear orthogonal coordinates in three-dimensional Euclidean space are derived based on the method of differential geometry. And the formulas ofAand ∇2Ain the Cartesian coordinates, cylindrical coordinates and spherical coordinates are also given for convenience of interested students and teachers.

曲线正交坐标系,微分几何, A 2 A,Curvilinear Orthogonal Coordinates Differential Geometry A 2 A
摘要

电动力学教材一般仅仅在附录部分给出曲线正交坐标系下标量的梯度 ∇ → ψ 、矢量的散度 ∇ → ⋅ A → 和旋度 ∇ → × A → ,以及标量的拉普拉斯算符 ∇ 2 ψ 的一般表达式。在教学中我们发现, ∇ → A → 和 ∇ 2 A → 在曲线正交坐标系下的表达式也是很重要的。本文从微分几何的角度推导出三维欧氏空间中曲线正交坐标系下 ∇ → A → 和 ∇ 2 A → 的一般表达式,以及它们在直角坐标系、柱坐标系和球坐标系的具体表达式,供感兴趣的学生和教师参考。

关键词

曲线正交坐标系,微分几何, ∇ → A → , ∇ 2 A →

A Derivation of ∇ → A → and ∇ 2 A → in Curvilinear Orthogonal Coordinates<sup> </sup>

Gerui Chen*, Mengdong He, Xiaoyan Chen

College of Science, Central South University of Forestry and Technology, Changsha Hunan

Received: May 13th, 2021; accepted: Jun. 18th, 2021; published: Jun. 28th, 2021

ABSTRACT

In textbooks of electrodynamics, the formulas of the gradient ∇ → ψ , divergence ∇ → ⋅ A → , curl ∇ → × A → ,and Laplace operator ∇ 2 ψ in curvilinear orthogonal coordinates are often given in the appendix. However, in the teaching process, we find that the expressions of ∇ → A → and ∇ 2 A → in curvilinear orthogonal coordinates are also very important. In this paper, the general expressions of ∇ → A → and ∇ 2 A → in curvilinear orthogonal coordinates in three-dimensional Euclidean space are derived based on the method of differential geometry. And the formulas of ∇ → A → and ∇ 2 A → in the Cartesian coordinates, cylindrical coordinates and spherical coordinates are also given for convenience of interested students and teachers.

Keywords:Curvilinear Orthogonal Coordinates, Differential Geometry, ∇ → A → , ∇ 2 A →

Copyright © 2021 by author(s) and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

在电动力学课程中,一般在讲授正式内容以前要先讲授教材附录部分中如下8个 ∇ → 算符运算公式的证明,因为它们在后面的课程中要频繁用到 [ 1 ]。

∇ → ( φ ψ ) = φ ∇ → ψ + ψ ∇ → φ (1)

∇ → ⋅ ( φ A → ) = ( ∇ → φ ) ⋅ A → + φ ∇ → ⋅ A → (2)

∇ → × ( φ A → ) = ( ∇ → φ ) × A → + φ ∇ → × A → (3)

∇ → ⋅ ( A → × B → ) = ( ∇ → × A → ) ⋅ B → − A → ⋅ ( ∇ → × B → ) (4)

∇ → × ( A → × B → ) = ( B → ⋅ ∇ → ) A → + ( ∇ → ⋅ B → ) A → − ( A → ⋅ ∇ → ) B → − ( ∇ → ⋅ A → ) B → (5)

∇ → ( A → ⋅ B → ) = A → × ( ∇ → × B → ) + ( A → ⋅ ∇ → ) B → + B → × ( ∇ → × A → ) + ( B → ⋅ ∇ → ) A → (6)

∇ → ⋅ ∇ → φ = ∇ 2 φ (7)

∇ → × ( ∇ → × A → ) = ∇ → ( ∇ → ⋅ A → ) − ∇ 2 A → (8)

电动力学课程中,除了直角坐标系,还经常用到柱坐标系和球坐标系等曲线正交坐标系。教材一般在附录部分不加证明地给出曲线正交坐标系下标量的梯度 ∇ → ψ 、矢量的散度 ∇ → ⋅ A → 和旋度 ∇ → × A → ,以及标量的拉普拉斯算符 ∇ 2 ψ 的一般表达式。然而仅仅利用这四个式子,上面的第(5)、(6)、(8)式却无法求出,因为里面涉及到 ∇ → A → 与 ∇ 2 A → 的运算。上面第(5)、(6)式需要利用 ∇ → A → ,而第(8)式需要利用 ∇ 2 A → 。电动力学课程中,电磁波的传播是重要的一章。时谐电磁波满足亥姆赫兹方程 ∇ 2 E → + k E → = 0 ,解此偏微分方程需要根据问题的对称性选择合适的坐标系。教材中讨论矩形谐振腔和波导,选择直角坐标系方便,但是如果谐振腔和波导具有球对称性和柱对称性,则需要 ∇ 2 E → 在球坐标系和柱坐标系的表达式。由此可见 ∇ → A → 与 ∇ 2 A → 在曲线正交坐标系下的表达式是重要的。本文从微分几何的角度推导三维欧氏空间中曲线正交坐标系下 ∇ → A → 与 ∇ 2 A → 的一般表达式。微分几何虽然是数学系本科生学习的内容,但是对学过高等数学和线性代数的二年级本科生在教师的指点下也可以大致掌握基本的思想和方法 [ 2 ]。作者讲授应用物理专业电动力学课程,课程最后要讲授狭义相对论,学生学完后就具有了一定的几何思想。有不少同学对狭义相对论之后的广义相对论兴趣十分浓厚,往往希望教师能够给出一些学习建议。近代物理的重要内容广义相对论和规范场论都是建立在整体微分几何的基础上的。鉴于近代微分几何对物理工作者的重要性,我们认为在本科教学中适当指导学生了解一点微分几何是极有帮助的,这对激发学生探索物理的兴趣起到很大的作用。文中详细展现每一步的推导计算过程,这样有助于读者初步了解微分几何的思想和方法。本文的计算采用著作 [ 2 ] 的符号约定。

2. 准备工作

三维欧氏空间在笛卡尔坐标系下的线元为:

d l 2 = d x d x + d y d y + d z d z (9)

设曲线正交坐标系的坐标为 u i ( i = 1 , 2 , 3 ) ,有变换关系:

x = x ( u 1 , u 2 , u 3 ) , y = y ( u 1 , u 2 , u 3 ) , z = z ( u 1 , u 2 , u 3 ) (10)

对x、y、z求全微分,带入线元(9)得到三维欧氏空间在曲线正交坐标系下的线元:

d l 2 = h 1 2 d u 1 d u 1 + h 2 2 d u 2 d u 2 + h 3 2 d u 3 d u 3 (11)

其中:

h i 2 = ( ∂ x ∂ u i ) 2 + ( ∂ y ∂ u i ) 2 + ( ∂ z ∂ u i ) 2 ,其中 i = 1 , 2 , 3 (12)

为标度因子。

由曲线正交坐标系下的线元表达式(11),可得度规:

g i j = ( h 1 2 0 0 0 h 2 2 0 0 0 h 3 2 ) ,其中 i , j = 1 , 2 , 3 (13)

及度规的逆:

g i j = ( 1 h 1 2 0 0 0 1 h 2 2 0 0 0 1 h 3 2 ) ,其中 i , j = 1 , 2 , 3 (14)

我们采用罗杰·彭罗斯1的抽象指标记号,矢量 A → 用记号 A a 表示,它在曲线正交坐标系的正交归一基底 ( e i ) a (即 e → i )下可表示为:

A a = A 1 ( e 1 ) a + A 2 ( e 2 ) a + A 3 ( e 3 ) a (15)

其中 A i ( i = 1 , 2 , 3 ) 是 A a 在正交归一基底的分量。

矢量 A a 在坐标基底 ( ∂ ∂ u i ) a 下可表示为:

A a = A u 1 ( ∂ ∂ u 1 ) a + A u 2 ( ∂ ∂ u 2 ) a + A u 3 ( ∂ ∂ u 3 ) a (16)

其中 A u i ( i = 1 , 2 , 3 ) 是矢量 A a 在坐标基底的分量。

注意正交归一基底和坐标基底有如下的关系:

( e i ) a = 1 h i ( ∂ ∂ u i ) a (17)

矢量 A a 在正交归一基底的分量 A i 与坐标基底的分量 A u i 有关系:

A u i = A i h i (18)

3. ∇ → A → 的计算

∇ → A → 用抽象指标可表示为 ∇ → A → = ∇ a A b ,则:

∇ → A → = ∇ a A b = g a c ∇ c A b = g a c ( ∂ c A b + Γ c e b A e ) = g a c ∂ c A b + g a c Γ c e b A e (19)

∇ a 是跟三维欧氏度规相适配的导数算符,与曲线正交坐标系的导数算符 ∂ a 有下面的关系:

∇ c A b = ∂ c A b + Γ c e b A e (20)

其中 Γ a b d 是反映 ∇ a 与 ∂ a 区别的克氏符,注意克氏符下面的两个指标是对称的,即 Γ a b d = Γ b a d 。

根据克氏符的定义 Γ μ ν σ = 1 2 g σ ρ ( g ρ μ , ν + g ν ρ , μ − g μ ν , ρ ) ,利用度规(13)以及度规的逆(14)可以算出非零克氏符,如下2

Γ 11 1 = ∂ 1 h 1 h 1 , Γ 12 1 = Γ 21 1 = ∂ 2 h 1 h 1 , Γ 13 1 = ∂ 3 h 1 h 1 , Γ 22 1 = − h 2 ∂ 1 h 2 h 1 2 , Γ 31 1 = ∂ 3 h 1 h 1 , Γ 33 1 = − h 3 ∂ 1 h 3 h 1 2 ,

Γ 11 2 = − h 1 ∂ 2 h 1 h 2 2 , Γ 12 2 = Γ 21 2 = ∂ 1 h 2 h 2 , Γ 22 2 = ∂ 2 h 2 h 2 , Γ 32 2 = ∂ 3 h 2 h 2 , Γ 33 2 = − h 3 ∂ 2 h 3 h 2 2 , Γ 11 3 = − h 1 ∂ 3 h 1 h 3 2 (21)

Γ 13 3 = Γ 31 3 = ∂ 1 h 3 h 3 , Γ 22 3 = − h 2 ∂ 3 h 2 h 3 2 , Γ 23 3 = Γ 32 3 = ∂ 2 h 3 h 3 , Γ 33 3 = ∂ 3 h 3 h 3

注意: ∂ i = ∂ u i = ∂ ∂ u i ,( i = 1 , 2 , 3 )。

把(19)式用曲线正交坐标基底 ( ∂ ∂ u i ) a 展开为:

g a c ∂ c A b + g a c Γ c e b A e = C 11 ( ∂ ∂ u 1 ) a ( ∂ ∂ u 1 ) b + C 12 ( ∂ ∂ u 1 ) a ( ∂ ∂ u 2 ) b + C 13 ( ∂ ∂ u 1 ) a ( ∂ ∂ u 3 ) b + C 21 ( ∂ ∂ u 2 ) a ( ∂ ∂ u 1 ) b + C 22 ( ∂ ∂ u 2 ) a ( ∂ ∂ u 2 ) b + C 23 ( ∂ ∂ u 2 ) a ( ∂ ∂ u 3 ) b + C 31 ( ∂ ∂ u 3 ) a ( ∂ ∂ u 1 ) b + C 32 ( ∂ ∂ u 3 ) a ( ∂ ∂ u 2 ) b + C 33 ( ∂ ∂ u 3 ) a ( ∂ ∂ u 3 ) b (22)

其中:

C 11 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 1 ) a ( d u 1 ) b , C 12 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 1 ) a ( d u 2 ) b

C 13 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 1 ) a ( d u 3 ) b , C 21 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 2 ) a ( d u 1 ) b

C 22 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 2 ) a ( d u 2 ) b , C 23 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 2 ) a ( d u 3 ) b (23)

C 31 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 3 ) a ( d u 1 ) b , C 32 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 3 ) a ( d u 2 ) b

C 33 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 3 ) a ( d u 3 ) b

其中 ( d u i ) a 是 ( ∂ ∂ u i ) a 的对偶坐标基底。

利用正交归一基底和坐标基底的关系(17)式,则(19)式用曲线正交坐标系的正交归一基底展开为:

g a c ∂ c A b + g a c Γ c e b A e = C 11 h 1 h 1 ( e 1 ) a ( e 1 ) b + C 12 h 1 h 2 ( e 1 ) a ( e 2 ) b + C 13 h 1 h 3 ( e 1 ) a ( e 3 ) b + C 21 h 2 h 1 ( e 2 ) a ( e 1 ) b + C 22 h 2 h 2 ( e 2 ) a ( e 2 ) b + C 23 h 2 h 3 ( e 2 ) a ( e 3 ) b + C 31 h 3 h 1 ( e 3 ) a ( e 1 ) b + C 32 h 3 h 2 ( e 3 ) a ( e 2 ) b + C 33 h 3 h 3 ( e 3 ) a ( e 3 ) b (24)

我们这里计算(24)式右边第一项。(23)式第一项为:

C 11 = ( g a c ∂ c A b + g a c Γ c e b A e ) ( d u 1 ) a ( d u 1 ) b = g 1 i ∂ i A u 1 + g 1 i Γ i e 1 A e = g 11 ∂ 1 A u 1 + g 11 Γ 1 e 1 A e = g 11 ∂ 1 A u 1 + g 11 Γ 11 1 A u 1 + g 11 Γ 12 1 A u 2 + g 11 Γ 13 1 A u 3 = g 11 ∂ 1 ( A 1 h 1 ) + g 11 Γ 11 1 A 1 h 1 + g 11 Γ 12 1 A 2 h 2 + g 11 Γ 13 1 A 3 h 3 (25)

上面最后一个等号我们利用了关系(18)。则得到:

C 11 h 1 h 1 = 1 h 1 ∂ 1 A 1 + 1 h 1 h 2 ( ∂ 2 h 1 ) A 2 + 1 h 1 h 3 ( ∂ 3 h 1 ) A 3 (26)

(24)式中其它的系数 C 12 h 1 h 2 ⋅ ⋅ ⋅ 可用类似的计算得到。最后把所求得的系数带入(24)式得到:

∇ → A → = [ 1 h 1 ∂ 1 A 1 + 1 h 1 h 2 ( ∂ 2 h 1 ) A 2 + 1 h 1 h 3 ( ∂ 3 h 1 ) A 3 ] e → 1 e → 1     + [ 1 h 1 ∂ 1 A 2 − 1 h 1 h 2 ( ∂ 2 h 1 ) A 1 ] e → 1 e → 2   + [ 1 h 1 ∂ 1 A 3 − 1 h 1 h 3 ( ∂ 3 h 1 ) A 1 ] e → 1 e → 3     + [ 1 h 2 ∂ 2 A 1 − 1 h 1 h 2 ( ∂ 1 h 2 ) A 2 ] e → 2 e → 1   + [ 1 h 2 ∂ 2 A 2 + 1 h 1 h 2 ( ∂ 1 h 2 ) A 1 + 1 h 2 h 3 ( ∂ 3 h 2 ) A 3 ] e → 2 e → 2   + [ 1 h 2 ∂ 2 A 3 − 1 h 2 h 3 ( ∂ 3 h 2 ) A 2 ] e → 2 e → 3 + [ 1 h 3 ∂ 3 A 1 − 1 h 1 h 3 ( ∂ 1 h 3 ) A 3 ] e → 3 e → 1   + [ 1 h 3 ∂ 3 A 2 − 1 h 2 h 3 ( ∂ 2 h 3 ) A 3 ] e → 3 e → 2 + [ 1 h 3 ∂ 3 A 3 + 1 h 1 h 3 ( ∂ 1 h 3 ) A 1 + 1 h 3 h 2 ( ∂ 2 h 3 ) A 2 ] e → 3 e → 3 (27)

这就是 ∇ → A → 在曲线正交坐标系的一般表达式,目前我们在通常的教科书中没有找到这个结果。常见的曲线正交坐标系有直角坐标系、柱坐标系和球坐标系,我们以这三个例子为例,得出具体的表达式。

1) 直角坐标系下

u 1 = x , u 2 = y , u 3 = z

h 1 = 1 , h 2 = 1 , h 3 = 1

∇ → A → = ∂ x A 1 e → x e → x + ∂ x A 2 e → x e → y + ∂ x A 3 e → x e → z + ∂ y A 1 e → y e → x + ∂ y A 2 e → y e → y + ∂ y A 3 e → y e → z + ∂ z A 1 e → z e → x + ∂ z A 2 e → z e → y + ∂ z A 3 e → z e → z (28)

2) 柱坐标系下

u 1 = r , u 2 = ϕ , u 3 = z

利用变换关系

x = r cos ϕ , y = r sin ϕ , z = z

带入(12)式可算出:

h 1 = 1 , h 2 = r , h 3 = 1

则得到:

∇ → A → = ∂ r A 1 e → r e → r + ∂ r A 2 e → r e → ϕ + ∂ r A 3 e → r e → z + [ 1 r ∂ ϕ A 1 − 1 r A 2 ] e → ϕ e → r + [ 1 r ∂ ϕ A 2 + 1 r A 1 ] e → ϕ e → ϕ + 1 r ∂ 2 A 3 e → ϕ e → z + ∂ z A 1 e → z e → r + ∂ z A 2 e → z e → ϕ + ∂ z A 3 e → z e → z (29)

3) 球坐标系下

u 1 = r , u 2 = θ , u 3 = ϕ

h 1 = 1 , h 2 = r , h 3 = r sin θ

∇ → A → = ∂ r A 1 e → r e → r + ∂ r A 2 e → r e → θ + ∂ r A 3 e → r e → ϕ + [ 1 r ∂ θ A 1 − 1 r A 2 ] e → θ e → r + [ 1 r ∂ θ A 2 + 1 r A 1 ] e → θ e → θ + 1 r ∂ θ A 3 e → θ e → ϕ + [ 1 r sin θ ∂ ϕ A 1 − 1 r A 3 ] e → ϕ e → r + [ 1 r sin θ ∂ ϕ A 2 − cos θ r sin θ A 3 ] e → ϕ e → θ + [ 1 r sin θ ∂ ϕ A 3 + 1 r A 1 + cos θ r sin θ A 2 ] e → ϕ e → ϕ (30)

4. 计算 ∇ 2 A →

∇ 2 A → 用抽象指标可表示为 ∇ a ∇ a A b ,则:

∇ 2 A → = ∇ a ∇ a A b = g c a ∇ c ∇ a A b (31)

= g c a [ ∂ c ( ∇ a A b ) + Γ c d b ∇ a A d − Γ c a d ∇ d A d ]

= g c a { ∂ c ( ∂ a A b + Γ a d b A d ) + Γ c d b ( ∂ a A d + Γ a e d A e ) − Γ c a d ( ∂ d A b + Γ d e b A e ) }

= g c a ∂ c ∂ a A b + g c a ∂ c ( Γ a d b ) A d + g c a Γ a d b ∂ c A d + g c a Γ c d b ∂ a A d + g c a Γ c d b Γ a e d A e − g c a Γ c a d ∂ d A b − g c a Γ c a d Γ d e b A e (32)

(31)式用曲线正交坐标基底 ( ∂ ∂ u i ) a 展开为:

∇ 2 A → = g c a ∇ c ∇ a A b = C 1 ( ∂ ∂ u 1 ) b + C 2 ( ∂ ∂ u 2 ) b + C 3 ( ∂ ∂ u 3 ) b (33)

其中:

C 1 = g c a ∇ c ∇ a A b ( d u 1 ) b C 2 = g c a ∇ c ∇ a A b ( d u 2 ) b C 3 = g c a ∇ c ∇ a A b ( d u 3 ) b (34)

利用正交归一基底和坐标基底的关系(17)式,(31)式可用曲线正交坐标系的正交归一基底展开为:

∇ 2 A → = g c a ∇ c ∇ a A b = C 1 h 1 ( e 1 ) b + C 2 h 2 ( e 2 ) b + C 3 h 3 ( e 3 ) b (35)

下面我们求 C 1 h 1 。

C 1 h 1 = g c a ∇ c ∇ a A b ( d u 1 ) b = h 1 g c a ∂ c ∂ a A b ( d u 1 ) b + h 1 g c a ∂ c ( Γ a d b ) A d ( d u 1 ) b + h 1 g c a Γ a d b ∂ c A d ( d u 1 ) b     + h 1 g c a Γ c d b ∂ a A d ( d u 1 ) b + h 1 g c a Γ c d b Γ a e d A e ( d u 1 ) b − h 1 g c a Γ c a d ∂ d A b ( d u 1 ) b     − h 1 g c a Γ c a d Γ d e b A e ( d u 1 ) b (36)

上式第二个等号用了(32)式。我们以(36)式第二个等号后的第二项的计算为例。

h 1 g c a ∂ c ( Γ a d b ) A d ( d u 1 ) b = h 1 g i j ∂ i ( Γ j d 1 ) A d = h 1 g 11 ∂ 1 ( Γ 1 d 1 ) A d + h 1 g 22 ∂ 2 ( Γ 2 d 1 ) A d + h 1 g 33 ∂ 3 ( Γ 3 d 1 ) A d = h 1 g 11 ∂ 1 ( Γ 11 1 ) A 1 h 1 + h 1 g 11 ∂ 1 ( Γ 12 1 ) A 2 h 2 + h 1 g 11 ∂ 1 ( Γ 13 1 ) A 3 h 3 + h 1 g 22 ∂ 2 ( Γ 21 1 ) A 1 h 1 + h 1 g 22 ∂ 2 ( Γ 22 1 ) A 2 h 2 + h 1 g 33 ∂ 2 ( Γ 23 1 ) A 3 h 3 + h 1 g 33 ∂ 3 ( Γ 31 1 ) A 1 h 1 + h 1 g 33 ∂ 3 ( Γ 32 1 ) A 2 h 2 + h 1 g 33 ∂ 3 ( Γ 33 1 ) A 3 h 3 (37)

同样的办法求出(36)式中剩下的6项,即可求得 (35)式中 ( e 1 ) b 前的系数:

C 1 h 1 = A 1 ( − 1 ) { ( ∂ 3 h 1 ) 2 h 1 2 h 3 2 + ( ∂ 2 h 1 ) 2 h 1 2 h 2 2 + ( ∂ 1 h 2 ) 2 h 1 2 h 2 2 + ( ∂ 1 h 3 ) 2 h 1 2 h 3 2 }   + ∂ 1 A 1 [ − ∂ 1 h 1 h 1 3 + ∂ 1 h 2 h 1 2 h 2 + ∂ 1 h 3 h 1 2 h 3 ] + ∂ 2 A 1 [ ∂ 2 h 1 h 1 h 2 2 − ∂ 2 h 2 h 2 3 + ∂ 2 h 3 h 2 2 h 3 ]   + ∂ 3 A 1 [ ∂ 3 h 1 h 1 h 3 2 − ∂ 3 h 3 h 3 3 + ∂ 3 h 2 h 2 h 3 2 ] + ∂ 1 ∂ 1 A 1 1 h 1 2 + ∂ 2 ∂ 2 A 1 1 h 2 2 + ∂ 3 ∂ 3 A 1 1 h 3 2   + A 2 [ − ∂ 2 h 1 ∂ 1 h 1 h 1 3 h 2 + ∂ 2 h 2 ∂ 1 h 2 h 1 h 2 3 − ∂ 2 h 3 ∂ 1 h 2 h 1 h 2 2 h 3 + ∂ 2 h 1 ∂ 1 h 3 h 1 2 h 2 h 3 − ∂ 2 h 3 ∂ 1 h 3 h 1 h 2 h 3 2 + ∂ 1 ∂ 2 h 1 h 1 2 h 2 − ∂ 1 ∂ 2 h 2 h 1 h 2 2 ]   + ∂ 1 A 2 2 ∂ 2 h 1 h 1 2 h 2 − ∂ 2 A 2 2 ∂ 1 h 2 h 1 h 2 2 + ∂ 1 A 3 2 ∂ 3 h 1 h 1 2 h 3 − ∂ 3 A 3 2 ∂ 1 h 3 h 1 h 3 2   + A 3 [ − ∂ 3 h 1 ∂ 1 h 1 h 1 3 h 3 + ∂ 3 h 1 ∂ 1 h 2 h 1 2 h 2 h 3 − ∂ 3 h 2 ∂ 1 h 2 h 1 h 2 2 h 3 − ∂ 3 h 2 ∂ 1 h 3 h 1 h 2 h 3 2 + ∂ 3 h 3 ∂ 1 h 3 h 1 h 3 3 + ∂ 1 ∂ 3 h 1 h 1 2 h 3 − ∂ 1 ∂ 3 h 3 h 1 h 3 2 ] (38)

根据对称性,把(38)式利用指标轮换( 1 → 2 ,   2 → 3 ,   3 → 1 )即可直接写出 ( e 2 ) b 前的系数:

C 2 h 2 = A 2 ( − 1 ) { ( ∂ 1 h 2 ) 2 h 2 2 h 1 2 + ( ∂ 3 h 2 ) 2 h 2 2 h 3 2 + ( ∂ 2 h 3 ) 2 h 2 2 h 3 2 + ( ∂ 2 h 1 ) 2 h 2 2 h 1 2 }   + ∂ 2 A 2 [ − ∂ 2 h 2 h 2 3 + ∂ 2 h 3 h 2 2 h 3 + ∂ 2 h 1 h 2 2 h 1 ] + ∂ 3 A 2 [ ∂ 3 h 2 h 2 h 3 2 − ∂ 3 h 3 h 3 3 + ∂ 3 h 1 h 3 2 h 1 ]   + ∂ 1 A 2 [ ∂ 1 h 2 h 2 h 1 2 − ∂ 1 h 1 h 1 3 + ∂ 1 h 3 h 3 h 1 2 ] + ∂ 2 ∂ 2 A 2 1 h 2 2 + ∂ 3 ∂ 3 A 2 1 h 3 2 + ∂ 1 ∂ 1 A 2 1 h 1 2   + A 3 [ − ∂ 3 h 2 ∂ 2 h 2 h 2 3 h 3 + ∂ 3 h 3 ∂ 2 h 3 h 2 h 3 3 − ∂ 3 h 1 ∂ 2 h 3 h 2 h 3 2 h 1 + ∂ 3 h 2 ∂ 2 h 1 h 2 2 h 3 h 1 − ∂ 3 h 1 ∂ 2 h 1 h 2 h 3 h 1 2 + ∂ 2 ∂ 3 h 2 h 2 2 h 3 − ∂ 2 ∂ 3 h 3 h 2 h 3 2 ]   + ∂ 2 A 3 2 ∂ 3 h 2 h 2 2 h 3 − ∂ 3 A 3 2 ∂ 2 h 3 h 2 h 3 2 + ∂ 2 A 1 2 ∂ 1 h 2 h 2 2 h 1 − ∂ 1 A 1 2 ∂ 2 h 1 h 2 h 1 2   + A 1 [ − ∂ 1 h 2 ∂ 2 h 2 h 2 3 h 1 + ∂ 1 h 2 ∂ 2 h 3 h 2 2 h 3 h 1 − ∂ 1 h 3 ∂ 2 h 3 h 2 h 3 2 h 1 − ∂ 1 h 3 ∂ 2 h 1 h 2 h 3 h 1 2 + ∂ 1 h 1 ∂ 2 h 1 h 2 h 1 3 + ∂ 2 ∂ 1 h 2 h 2 2 h 1 − ∂ 2 ∂ 1 h 1 h 2 h 1 2 ] (39)

同理把(39)式利用指标轮换( 1 → 2 ,   2 → 3 ,   3 → 1 )可直接写出 ( e 3 ) b 前的系数:

C 3 h 3 = A 3 ( − 1 ) { ( ∂ 2 h 3 ) 2 h 3 2 h 2 2 + ( ∂ 1 h 3 ) 2 h 3 2 h 1 2 + ( ∂ 3 h 1 ) 2 h 3 2 h 1 2 + ( ∂ 3 h 2 ) 2 h 3 2 h 2 2 }   + ∂ 3 A 3 [ − ∂ 3 h 3 h 3 3 + ∂ 3 h 1 h 3 2 h 1 + ∂ 3 h 2 h 3 2 h 2 ] + ∂ 1 A 3 [ ∂ 1 h 3 h 3 h 1 2 − ∂ 1 h 1 h 1 3 + ∂ 1 h 2 h 1 2 h 2 ]   + ∂ 2 A 3 [ ∂ 2 h 3 h 3 h 2 2 − ∂ 2 h 2 h 2 3 + ∂ 2 h 1 h 1 h 2 2 ] + ∂ 3 ∂ 3 A 3 1 h 3 2 + ∂ 1 ∂ 1 A 3 1 h 1 2 + ∂ 2 ∂ 2 A 3 1 h 2 2   + A 1 [ − ∂ 1 h 3 ∂ 3 h 3 h 3 3 h 1 + ∂ 1 h 1 ∂ 3 h 1 h 3 h 1 3 − ∂ 1 h 2 ∂ 3 h 1 h 3 h 1 2 h 2 + ∂ 1 h 3 ∂ 3 h 2 h 3 2 h 1 h 2 − ∂ 1 h 2 ∂ 3 h 2 h 3 h 1 h 2 2 + ∂ 3 ∂ 1 h 3 h 3 2 h 1 − ∂ 3 ∂ 1 h 1 h 3 h 1 2 ]   + ∂ 3 A 1 2 ∂ 1 h 3 h 3 2 h 1 − ∂ 1 A 1 2 ∂ 3 h 1 h 3 h 1 2 + ∂ 3 A 2 2 ∂ 2 h 3 h 3 2 h 2 − ∂ 2 A 2 2 ∂ 3 h 2 h 3 h 2 2   + A 2 [ − ∂ 2 h 3 ∂ 3 h 3 h 3 3 h 2 + ∂ 2 h 3 ∂ 3 h 1 h 3 2 h 1 h 2 − ∂ 2 h 1 ∂ 3 h 1 h 3 h 1 2 h 2 − ∂ 2 h 1 ∂ 3 h 2 h 3 h 1 h 2 2 + ∂ 2 h 2 ∂ 3 h 2 h 3 h 2 3 + ∂ 3 ∂ 2 h 3 h 3 2 h 2 − ∂ 3 ∂ 2 h 2 h 3 h 2 2 ] (40)

最后把(38)式、(39)式、(40)式带入(35)式就得到了曲线正交坐标系下 ∇ 2 A → 的一般表达式。需要说明的是, ∇ 2 A → 在正交曲线坐标系的一般表达式在题解 [ 3 ] 中附录部分不加证明地直接给出。题解 [ 3 ] 给出的是相对紧凑的形式,而我们给出的是完全展开的形式,简单计算容易发现两个结果是一样的。

下面我们给出 ∇ 2 A → 在直角坐标系、柱坐标系和球坐标系的具体表达式。

1) 直角坐标系

u 1 = x , u 2 = y , u 3 = z

h 1 = 1 , h 2 = 1 , h 3 = 1

∇ 2 A → = ∇ 2 A 1 e → x + ∇ 2 A 2 e → y + ∇ 2 A 3 e → z (41)

2) 柱坐标系

u 1 = r , u 2 = ϕ , u 3 = z

h 1 = 1 , h 2 = r , h 3 = 1

∇ 2 A → = ( − A 1 r 2 + ∂ 1 A 1 1 r + ∂ 1 ∂ 1 A 1 + ∂ 2 ∂ 2 A 1 1 r 2 + ∂ 3 ∂ 3 A 1 − ∂ 2 A 2 2 r 2 ) e → r + ( 1 r ∂ 1 A 2 + ∂ 1 ∂ 1 A 2 + ∂ 2 ∂ 2 A 2 1 r 2 + ∂ 3 ∂ 3 A 2 − A 2 r 2 + 2 r 2 ∂ 2 A 1 ) e → ϕ + ( 1 r ∂ 1 A 3 + ∂ 1 ∂ 1 A 3 + ∂ 2 ∂ 2 A 3 1 r 2 + ∂ 3 ∂ 3 A 3 ) e → z (42)

3) 球坐标系下

u 1 = r , u 2 = θ , u 3 = ϕ

h 1 = 1 , h 2 = r , h 3 = r sin θ

∇ 2 A → = [ − 2 r 2 A 1 + ∂ 1 A 1 2 r + ∂ 2 A 1 cos θ r 2 sin θ + ∂ 1 ∂ 1 A 1 + ∂ 2 ∂ 2 A 1 1 r 2 + ∂ 3 ∂ 3 A 1 1 r 2 ( sin θ ) 2     − A 2 2 cos θ r 2 sin θ − ∂ 2 A 2 2 r 2 − ∂ 3 A 3 2 r 2 sin θ ] e → r     + [ − 1 r 2 ( sin θ ) 2 A 2 + ∂ 2 A 2 cos θ r 2 sin θ + ∂ 1 A 2 2 r + ∂ 2 ∂ 2 A 2 1 r 2 + ∂ 3 ∂ 3 A 2 1 r 2 ( sin θ ) 2     + ∂ 1 ∂ 1 A 2 − ∂ 3 A 3 2 cos θ r 2 ( sin θ ) 2 + ∂ 2 A 1 2 r 2 ] e → θ     + [ − 1 r 2 ( sin θ ) 2 A 3 + ∂ 1 A 3 2 r + ∂ 2 A 3 cos θ r 2 sin θ + ∂ 3 ∂ 3 A 3 1 r 2 ( sin θ ) 2 + ∂ 1 ∂ 1 A 3     + ∂ 2 ∂ 2 A 3 1 r 2 + ∂ 3 A 1 2 r 2 sin θ + ∂ 3 A 2 2 cos θ r 2 ( sin θ ) 2 ] e → ϕ (43)

这样我们就得到了曲线正交坐标系下 ∇ → A → 与 ∇ 2 A → 在直角坐标系、柱坐标系和球坐标系的具体表达式。有兴趣的读者可以照着上面的步骤得出其它曲线正交坐标系下 ∇ → A → 与 ∇ 2 A → 的具体表达式。下面我们给出椭球坐标系与直角坐标系的变换关系,

x = r 2 + a 2 sin θ cos ϕ , y = r 2 + a 2 sin θ sin ϕ , z = r cos θ (44)

请读者朋友自己完成。

5. 结论

本文从微分几何的角度推导了曲线正交坐标系下 ∇ → A → 与 ∇ 2 A → 的一般表达式,以及它们在三维直角坐标系、柱坐标系和球坐标系的具体表达式,这些表达式在相对复杂的物理问题中会涉及到。近代物理的重要内容广义相对论和规范场论都是建立在整体微分几何的基础上的,因此在本科物理教学阶段适当指导学生了解一些微分几何知识对学生未来的发展是极有帮助的。文中详细展示了推导计算的每一步过程,这样有助于读者初步了解微分几何的基本思想和方法。

基金项目

湖南省教育厅一般项目19C1895,中南林业科技大学校级教改项目。

文章引用

陈葛锐,贺梦冬,陈小艳. 曲线正交坐标系下A与∇2A的一种推导A Derivation ofAand ∇2Ain Curvilinear Orthogonal Coordinates[J]. 力学研究, 2021, 10(02): 143-152. https://doi.org/10.12677/IJM.2021.102014

参考文献 References 郭硕鸿. 电动力学(第三版) [M]. 北京: 高等教育出版社, 2008. 梁灿彬, 周彬. 微分几何入门与广义相对论(上册) [M]. 北京: 科学出版社, 2006. 林璇英, 张之翔. 电动力学题解(第三版) [M]. 北京: 科学出版社, 2018.
Baidu
map