11梯度、散度和旋度
▽算子不是一个 矢量,除非你把它 作用在一个函数上,否则它没啥意义。但是,它在各个方面的表现确实又像一个矢量,只要你把▽算子的“ 作用”看成矢量的“ 相乘”。 一个 矢量一般来说有3种“ 乘法”:1、矢量 A和一个标量a 相乘:a A 。比如我把一个矢量 A大小变为原来的2倍,方向不变,那么这时候就可以写成2 A。2、矢量 A和一个矢量 B进行 点乘: A·B。这个点乘我们上面介绍很多了, A·B=| A|| B|Cosθ,这里就不说了。3、矢量 A和一个矢量 B进行 叉乘: A×B。这个叉乘跟点乘类似,也是我们单独针对矢量定义的另外一种乘法, A×B=| A|| B|Sinθ。大家可以看到,这个叉乘跟点乘唯一的区别就是: 点乘是两个矢量的大小乘以它们的 余弦值Cosθ,叉乘是两个矢量的大小乘以它们的 正弦值Sinθ(在直角三角形里,角的 对边和斜边的比为正弦Sinθ, 邻边和斜边的比值为余弦Cosθ)。 那么,同样的,我们的 ▽算子也有3种 作用方式:1、▽算子作用在一个 标量函数 z上:▽z。这个▽z我们上面说过了,它表示函数z的 梯度,它表示这个函数z变化最快的方向。2、▽算子跟一个 矢量函数E 点乘: ▽·E。这就表示E的 散度,我们开篇讲的高斯电场定律的左边就是电场E的散度,它就是表示成 ▽·E这样。3、▽算子跟一个 矢量函数E 叉乘: ▽×E。它叫E的 旋度,这个我们后面会再详细说。
这样,我们就以一种很自然的方式引出了这三个非常重要的概念: 梯度(▽z )、 散度(▽·E)和 旋度(▽×E)。大家可以看到,▽算子的这三种作用跟矢量的三种乘法是非常相似的,只不过▽是一个算子,它必须作用在一个 函数上才行,所以我们把上面的标量和矢量换成了 标量函数和 矢量函数。 我们在描述山的 高度的函数 z=f(x,y)的时候,不同的点(x,y)对应不同的山的高度,而山的高度只有大小没有方向,所以这是个标量函数,我们可以求它的 梯度▽z。但是, 电场E既有大小又有方向,这是一个矢量,所以我们可以用一个矢量函数 E=f(x,y)表示空间中不同点(x,y)的电场E的分布情况。那么对这种 矢量函数,我们就不能去求它的 梯度了,我们只能去求它的 散度▽·E和 旋度▽×E。
为了让大家对这些能够有更直观的概念,我们接下来就来仔细看看电场的散度▽·E。
12电场的散度
当我们把电场的散度写成 ▽·E这样的时候,我们会觉得:啊,好简洁!但是我们也知道▽算子的定义是这样的:
也就是说,我们最开始从 无穷小曲面的通量定义来的 散度和我们上面通过 偏导数定义来的 散度▽·指的是同一个东西。即:
13为何这两种散度是等价的?
很多人可能觉得难以理解,这两个东西的表达形式和来源都完全不一样,它们怎么会是同一个东西呢?但是它们确实是同一个东西,那我们为什么要弄两套东西出来呢?在最开始我也说了,通过 无穷小曲面的通量定义的散度很容易理解,跟麦克斯韦方程组的积分形式的通量也有非常大的联系,但是这种定义 不好计算(上面的例2,你用这种方式去求它的散度试试?),所以我们需要找一种能方便计算、实际可用的方式,这样才出现了▽·形式的散度。 至于为什么这两种形式是等价的,我给大家提供一个 简单的思路。因为这毕竟是面向大众的科普性质的文章,具体的证明过程我就不细说了。真正感兴趣的朋友可以顺着这个思路去完成自己的证明,或者来我的 社群(回复“ 社群”即可)里讨论。 证明思路:我们假设有一个边长分别为Δx、Δy、Δz的小长方体,空间中的电场为E(x,y,z),然后假设在这个长方体的正中心有一个点(x,y,z),那么这个电场通过这个长方体 前面(沿着x轴 正方向)的 电场就可以表示为: Ex(x+Δx/2,y,z)。Ex表示电场在x方向上的分量(因为我们是考虑长方体上表面的通量,所以只用考虑电场的x分量),因为中心坐标为(x,y,z),那么沿着x轴移动到表面的坐标自然就是( x+Δx/2,y,z)。而这个面的 面积为 ΔyΔz,那么通过前面的电通量就可以写成: Ex(x+Δx/2,y,z)·ΔyΔz。 同样的,通过长方体 后面(沿着x轴的 负方向)的 电通量,就可以写成 Ex(x-Δx/2,y,z)·ΔyΔz。因为这两个面的方向是 相反的(前面后面,一个沿着x轴正方向,一个沿着负方向),所以,这两个沿着x轴方向的面的电通量之和 Φx就应该是两者相减: Φx=( Ex(x+Δx/2,y,z)·ΔyΔz- Ex (x-Δx/2,y,z)·ΔyΔz)。 如果我们两边都除以 Δv(其中, Δv=ΔxΔyΔz),那么就得到: Φx/Δv=(Ex(x+Δx/2,y,z)- Ex (x-Δx/2,y,z))/Δx,然后你会发现等式的右边刚好就是 偏导数的 定义(标准的极限定义 )。也就是说, 电场通过沿着x轴的两个面(前后两面)的通量之和就等于电场的x分量对x的偏导数: Φx/Δv=Ex/x。 同样的,我们发现电场沿着 y轴的两面(左右两面)和 z轴的两面(上下两面)的电通量之和分别就等于电场的y分量和z分量对y和z的 偏导: Φy/Δv=Ey/y, Φz/Δv=Ez/z。然后我们把这三个式子 加起来, 左边就是电场通过六个面的通量除以体积,也就是通过这个长方体的通量除以体积,右边就是我们▽·E的形式,这分别就是我们上面两种 散度的表示方式, 证明完成。 这个证明一时半会没看懂也没关系,感兴趣的可以后面慢慢去琢磨。我只是想通过这种方式让大家明白 通过某一方向的两个面的通量跟 这方向的偏导数之间是存在这种对应关系的,这样我们就容易接受 无穷小曲面的通量和 ▽·这两种散度的定义方式了。 这两种散度的定义方式各有所长,比如我们在判断某一点的散度是否为零的时候,我用第一个定义,去看看包含这个点的无穷小曲面的通量是不是为零就行了。 如果这一点有电荷,那么这个无穷小曲面的电通量肯定就不为零,它的散度也就不为零;如果这个无穷小曲面没有包含电荷,那这一点的散度一定为0,这就是 高斯电场定律的 微分方程想要告诉我们的东西。但是,如果你要计算这一点的散度是多少,那还是乖乖的拿起▽·去计算吧。
14散度的几何意义
此外,跟梯度一样,散度这个名字也是非常形象的。很多人会跟你说散度表示的是“ 散开的程度”,这种说法很容易让初学者误解或者迷惑,比如一个正电荷产生的产生的如下的电场线,它看起来是散开的,所以很多就会认为这里 所有的点的散度都是不为零的,都是正的。


15方程一: 高斯电场定律
说了这么多,又是证明不同散度形式( 无穷小曲面的通量和▽·)的 等价性,又是说明不同散度理解方式的 同一性( 无穷小曲面的通量和 散开的程度),都是为了让大家从更多的维度全方位的理解 散度的概念,尽量避开初学者学习散度会遇到的各种坑。理解了这个散度的概念之后,我们再来看 麦克斯韦方程组的 第一个方程—— 高斯电场定律的 微分形式就非常容易理解了:


16方程二: 高斯磁场定律
理解了高斯电场定律的微分形式,那么 高斯磁场定律的微分形式就能轻松写出来了。因为现在还没有找到磁单极子,磁感线都是闭合的曲线,所以 闭合曲面的磁通量一定恒为0,这就是 高斯磁场定律积分形式的思想:

17旋度
静电和静磁的微分形式我们已经说完了,那么接下来就是 磁如何生电的 法拉第定律了。关于 法拉第是如何通过实验一步一步发现法拉第定律的内容,我在积分篇里已经详细说了,这里就不再多说。对 法拉第定律的 基本思想和 积分形式的内容还不太熟悉的请先去看上一篇
积分篇
的内容。 法拉第定律是法拉第对 电磁感应现象的一个总结,他发现只要一个曲面的 磁通量(B·a)发生了改变,那么就会在曲面的边缘感生出一个 旋涡状的 电场E出来。这个旋涡状的 感生电场我们是用电场的 环流来描述的,也就是 电场沿着曲面边界进行的线积分。



18矢量的叉乘
因为旋度是 ▽算子以 叉乘×的方式作用在 矢量场上,所以这里我们来简单的看一下 叉乘。两个矢量 A和 B的 点乘被定义为: A·B=| A|| B|Cosθ,它们的 叉乘则被定义为 A×B=| A|| B|Sinθ,其中θ为它们的夹角。单从这样看,它们之间的差别好像很小,只不过一个是乘以 余弦Cosθ,另一个是乘以 正弦Sinθ。 从它们的 几何意义来说, 点乘表示的是 投影,因为|OA|Cosθ刚好就是OA在OB上的投影,也就是OC的长度。如下图:


19方程三: 法拉第定律
好,知道了矢量的叉乘,知道了 ▽×E可以表示 电场的旋度,而且知道旋度的定义是: 无穷小非闭合曲面的环流和这个曲面的面积之比。那我们再来回过头看一看 法拉第定律的 积分形式:

积分篇
里说过,磁通量( B·a)的变化可以有两种方式: 磁场(B)的变化和 通过曲面面积(S)的变化,我们上面这种方式是把这两种情况都算在内。但是,还有的学者认为只有磁场(B)的变化产生的电场才算法拉第定律,所以法拉第定律还有另外一个版本:



20旋度的几何意义
我们知道 旋度的定义是 无穷小曲面的环流和面积的比值,但是它既然取了旋度这个名字,那么它跟 旋转应该还是有点关系的。我们变化的磁场感生出来的电场也是一个 旋涡状的电场。那么,是不是只要看起来像漩涡状的矢量场,它就一定有旋度呢?

这个问题我们在讨论 散度的时候也遇到过,很多初学者认为只要看起来发散的东西就是有散度的,然后我们通过分析知道这是不对的。一个点电荷产生静电场,只要在电荷处散度不为零的,在其他地方,虽然看起来是散开的,其实它的 散度是 零。如果我们放一个非常轻的橡皮筋在上面,除了电荷所在处,其它地方这个橡皮筋是不会被撑开的(即便会被冲走),所以其他地方的散度都为零。
同样的,在旋度这里,一个变换的磁场会产生一个旋涡状的电场,在旋涡的中心,在磁场变化的这个中心点这里,它的旋度肯定是不为零的。但是,在其它地方呢?从公式上看,其它地方的旋度一定为零,为什么?因为其他地方 并没有变化的磁场啊,所以按照 法拉第定律的 微分形式, 没有变化的磁场的地方的电场的旋度肯定是0。 跟 散度一样,我们不能仅凭一个感生电场是不是旋转状的来判断这点旋度是否为0,我们也需要借助一个小道具: 小风车。我们把一个小风车放在某一点上,如 果这个风车能转起来,就说明这点的旋度不为0。你只要把风车放在 感生电场中心 以外的地方,就会发现如果外层的电场线让小风车 顺时针转,内层的电场线就会让小风车 逆时针转,这两股力刚好 抵消了。最终风车不会转,所以旋度为0。

如果大家能理解 静电场除了中心点以外的地方 散度处处为零,那么理解 感生电场除了中心点以外的地方 旋度处处为零就不是什么难事。在非中心点的地方, 散度的流入流出两股力量抵消了,旋度顺时针逆时针的两股力量抵消了,为什么刚好他们能抵消呢?本质原因还是因为 这两种电场都是随着距离的平方反比减弱。如果它们不遵守 平方反比定律,那么你去计算里外的散度和旋度,它们就不再为零。 关于 旋度的事情就先说这么多,大家如果理解了旋度,对比法拉第定律的积分方程,要理解它的微分方程是很容易的。我前面花了很大的篇幅给大家讲了矢量的 点乘和 散度,作为类比,理解矢量的 叉乘和 旋度也不是什么难事,它们确实太相似了。
21方程四: 安培-麦克斯韦定律
讲完了磁生电的法拉第定律,我们麦克斯韦方程组就只剩最后一个电生磁的 安培-麦克斯韦定律了。它描述的是 电流和 变化的电场如何产生 旋涡状的感生磁场的,因为它电的来源有电流和变化的电场两项,所以它的形式也是最复杂的。方程的 积分形式如下(具体过程见
积分篇
):


22麦克斯韦方程组
至此,麦克斯韦方程组的 四个方程:描述 静电的 高斯电场定律、描述 静磁的 高斯磁场定律、描述 磁生电的 法拉第定律和描述 电生磁的 安培-麦克斯韦定律的 微分形式就都说完了。把它们都写下来就是这样:


从思想上来讲, 微分形式和 积分形式表达的思想是一样的,毕竟它们都是麦克斯韦方程组。它们的差别仅仅在于 积分形式是从 宏观的角度描述问题,我们面对的 宏观上的曲面,所以要用 通量和 环流来描述电场、磁场;而 微分形式是从 微观的角度来描述问题,这时候曲面缩小都 无穷小,我们面对的东西就变成了一个 点,所以我们使用 散度和 旋度来描述电场、磁场。 这一点是特别要强调的: 通量和 环流是定义在 曲面上的,而 散度和 旋度是定义在一个 点上的。我们可以说通过通过一个曲面的 通量或者沿曲面边界的 环流,但是当我们在说 散度和 旋度的时候,我们都是在说 一个点的 散度和 旋度。
理解了这些,你再回过头去看看 麦克斯韦方程组的 积分形式:

23结语
到这里,麦克斯韦方程组的 积分篇和 微分篇就都说完了。 长尾君在这两篇文章里先从零开始引出了 通量,然后从通量的概念慢慢引出了 麦克斯韦方程组的 积分形式,再从积分形式用“ 把曲面压缩到无穷小”推出了对应的 微分形式。整个过程我都极力做到“ 通俗但不失准确”,所有新概念的引出都会先做层层铺垫,绝不从天而降的抛出一个新东西。目的就是为了让多的人能够更好的了解麦克斯韦方程组,特别是让 中学生也能看懂,能理解麦克斯韦方程组的美妙,同时也激发出他们对科学的好奇和热爱之心,打消他们对“高深”科学的 畏惧之心: 看,这么高大上的麦克斯韦方程组,年纪轻轻的我也能看懂,也能掌握~ 此外,麦克斯韦方程组是真的很美,你掌握的物理知识越多,就会越觉得它美。我也更希望 大家是因为它的美而喜欢这个方程组,而不仅仅是因为它的“重要性”。我们也都知道,麦克斯韦写出这套方程组以后,就从方程推导出了 电磁波,当他把相关的参数代入进去算出电磁波的速度的时候,他惊呆了!他发现这个 电磁波的速度跟人们实验测量的 光速极为接近,于是他给出了一个大胆的预测: 光就是一种电磁波。

可惜的是,英年早逝的麦克斯韦(48岁去世)并没能看到他的预言被证实,人类直到他去世9年后,也就是1888年才由赫兹首次证实了“光是一种电磁波”。那么, 麦克斯韦是怎么从方程组导出电磁波的呢?既然我们已经学完了 麦克斯韦方程组,想必大家也很知道如何从这套方程组推导出电磁波的方程,然后亲眼见证“ 电磁波的速度等于光速”这一奇迹时刻。这部分的内容, 长尾科技下篇文章再说。
最后,这篇文章主要参考了《 电动力学导论》(格里菲斯)和《 麦克斯韦方程直观》(Daniel Fleisch),大家想对麦克斯韦方程组做进一步了解的可以看看这两本书,需要电子档的可以在后台回复“ 麦克斯韦方程组”。 最美的方程,愿你能懂她的美~