场论初步:梯度、散度、旋度与Hamilton算子完整知识体系
📐 核心知识架构思维导图
场论基础 (Field Theory Fundamentals)
|
________________|_________________
| | |
标量场与向量场 微分算子 场的性质与分类
| | |
___________|_______ _____|_____ ______|______
| | | | | |
标量场(φ) 向量场(F) ∇ ∇·∇× 保守场 无源场
温度/势能 力场/速度 梯度 散度旋度 ∇×F=0 ∇·F=0
| | | | | |
等值面 场线 Hamilton Laplace 势函数 管量守恒
方向导数 通量 算子 算子 存在性 不可压缩
| | | | | |
最速上升 环流 ∇φ ∇²φ 路径无关 源密度=0
第一章:场论的基本概念
1.1 场的定义
1.1.1 标量场
定义 1.1(标量场)
若对空间区域 内的每一点 ,都对应一个确定的标量 (数值),则称在 内确定了一个标量场,记作:
经典例子:
-
温度场: 表示空间各点在时刻 的温度
- 稳定温度场: 不随时间变化
- 非稳定温度场: 随时间变化
-
电势场: 表示静电场的电势分布
-
压强场: 表示流体内部各点的压强
-
密度场: 表示物质密度分布
-
高度场: 表示地形的海拔高度
1.1.2 向量场
定义 1.2(向量场)
若对空间区域 内的每一点 ,都对应一个确定的向量 ,则称在 内确定了一个向量场,记作:
或分量形式:
经典例子:
-
速度场: 表示流体在各点的流速
- 稳定流动: 不随时间变化
- 非稳定流动: 随时间变化
-
力场: 表示空间各点受到的力
- 引力场:
- 电场:
-
磁场: 表示磁感应强度分布
-
梯度场:由标量场的梯度形成的向量场
1.2 等值面与场线
1.2.1 等值面(等位面)
定义 1.3(等值面)
标量场 (常数)所确定的曲面称为该标量场的等值面或等位面。
性质:
- 不同常数值对应不同的等值面
- 等值面族充满整个场所在的空间
- 不同的等值面不会相交(除了奇点)
物理意义:
- 等温面:温度场的等值面
- 等压面:压强场的等值面
- 等势面:电势场的等值面
例子:
点电荷产生的电势场:
等势面方程:
即:
这是一族以原点为中心的同心球面。
1.2.2 场线(力线、流线)
定义 1.4(场线)
向量场 的场线是这样的曲线:曲线上每一点的切线方向都与该点的向量 的方向一致。
微分方程:
场线满足:
这是场线的对称形式微分方程。
物理意义:
- 电力线:电场的场线,表示电场方向
- 磁力线:磁场的场线,表示磁场方向
- 流线:速度场的场线,表示流体质点的瞬时运动方向
注意:
- 流线 ≠ 轨迹线(非定常流动时)
- 定常流动时,流线 = 轨迹线
例题 1(原文例题)
**问题:**求向量场 的场线方程。
解:
场线方程为:
从第一、二项:
积分:
即:
这是圆柱面。
从第二、三项:
由 ,可设 ,
积分:
场线方程:
或参数形式:
这是圆柱面上的螺旋线。
几何图像: 场线是绕 轴的圆柱螺旋线,随 增大,按右手螺旋向上旋转。
第二章:方向导数与梯度
2.1 方向导数
2.1.1 定义与概念
定义 2.1(方向导数)
设函数 在点 的某邻域内有定义, 是从 出发的射线, 是 上的点,以 表示 到 的距离。
若极限
存在,则称此极限为函数 在点 沿方向 的方向导数,记作:
其中 是方向 的方向余弦。
2.1.2 计算公式
定理 2.1(方向导数的计算)
若函数 在点 可微,则 在 点沿任一方向 的方向导数存在,且:
向量形式:
设 是方向 的单位向量,则:
例题 2
**问题:**求函数 在点 沿方向 的方向导数。
解:
第一步:计算偏导数
在点 :
第二步:单位化方向向量
第三步:计算方向导数
答案:
2.2 梯度
2.2.1 定义
定义 2.2(梯度)
设函数 在点 具有一阶偏导数,则称向量
为函数 在点 的梯度(gradient)。
记号:
其中 (nabla)称为 Hamilton 算子或梯度算子:
2.2.2 梯度的几何与物理意义
定理 2.2(梯度的性质)
-
方向导数与梯度的关系:
其中 是 与 的夹角。
-
最大增长方向:
- 当 (即 )时,方向导数达到最大值
- 梯度方向是函数增长最快的方向
- 梯度的模 表示最大变化率
-
与等值面的关系:
- 梯度垂直于等值面
- 梯度指向函数值增大的方向
几何解释:
考虑标量场 的等值面族。在点 处:
- 等值面的法向量方向是函数值变化最快的方向
- 这个方向正是梯度 的方向
- 在等值面内(切平面上)移动,函数值不变,方向导数为零
物理解释:
-
温度场: 指向温度上升最快的方向,热量沿 方向流动
-
电势场:电场强度 ,电场指向电势降低的方向
-
高度场: 指向最陡上坡方向,水流沿 方向流动
例题 3(原文例题)
**问题:**求函数 在点 处的梯度,并求出其在该点沿哪个方向方向导数最大?最大值是多少?
解:
计算梯度:
在点 :
最大方向导数方向:
梯度方向,即 或单位化后
最大值:
物理意义:
函数 的等值面是同心球面。在点 处,梯度指向径向外,沿径向离开原点时, 增长最快。
2.2.3 梯度的运算法则
定理 2.3(梯度的运算性质)
设 是可微函数, 是常数,则:
-
线性性:
-
乘积法则:
-
商法则:
-
复合函数: 若 ,则:
-
链式法则: 若 ,,则:
例题 4
**问题:**设 ,其中 ,求 。
解:
方法一:直接计算
类似地:
因此:
方法二:利用链式法则
其中 (径向单位向量)
物理意义:
这正是点电荷或点质量产生的场的形式:
- 电场:()
- 引力场:()
第三章:散度(Divergence)
3.1 通量的概念
3.1.1 物理背景
通量(flux)是向量场穿过曲面的总量的度量。
定义 3.1(通量)
设 是向量场 中的定向曲面,则向量场 通过曲面 的通量定义为:
其中 是 的单位法向量。
物理意义:
-
流体力学:流速场 通过曲面 的通量表示单位时间内穿过 的流体体积
-
电磁学:电场 通过曲面 的通量(电通量)
-
热传导:热流密度 通过曲面 的通量表示热流量
3.2 散度的定义
3.2.1 直观引入
考虑点 处的一个小闭合曲面 ,体积为 。向量场 通过 的通量为:
单位体积的通量:
当 时,这个比值的极限反映了点 处向量场的"发散程度"。
定义 3.2(散度)
设向量场 在点 的某邻域内有定义,若极限
存在,则称此极限为向量场 在点 的散度(divergence),记作:
3.2.2 散度的计算公式
定理 3.1(散度的坐标表示)
若向量场 的分量函数有一阶连续偏导数,则:
记号:
形式上看作 Hamilton 算子 与向量 的点积。
证明思路:
利用高斯公式:
应用积分中值定理,存在 ,使得:
因此:
当 时,,由散度的连续性:
3.3 散度的物理意义
3.3.1 源与汇
物理解释:
散度 表示向量场 在某点的"源密度"或"流出强度":
- :该点是源(source),有净流出
- :该点是汇(sink),有净流入
- :该点无源无汇,流入 = 流出
3.3.2 应用实例
1. 流体力学
连续性方程:
设流体密度为 ,速度场为 ,则质量守恒要求:
不可压缩流体( 常数):
2. 静电学
高斯定律(微分形式):
其中 是电荷密度, 是真空介电常数。
**物理意义:**电场的散度等于电荷密度(电荷是电场的"源")。
3. 热传导
热传导方程:
设温度场为 ,热流密度为 (傅里叶定律),则:
其中 是拉普拉斯算子。
3.4 散度的运算性质
定理 3.2(散度的运算法则)
设 是向量场, 是标量场, 是常数,则:
-
线性性:
-
数乘法则:
-
梯度的散度:
(拉普拉斯算子)
-
旋度的散度恒为零:
例题 5(原文例题)
**问题:**计算向量场 的散度。
解:
在点 :
例题 6
**问题:**证明向量场 (,)在 处散度为零。
解:
类似地:
相加:
物理意义:
点电荷或点质量产生的场在源点以外处处无散度(无其他源)。但在原点处,散度是 函数(狄拉克函数),表示点源。
第四章:旋度(Curl)
4.1 环流量的概念
4.1.1 物理背景
环流量(circulation)度量向量场沿闭曲线的"旋转强度"。
定义 4.1(环流量)
设 是向量场 中的定向闭曲线,则向量场 沿 的环流量定义为:
物理意义:
-
流体力学:速度场 沿闭曲线 的环流量表示流体绕 的总"旋转"
-
电磁学:磁场 沿闭曲线 的环流量与穿过 的电流有关(安培定律)
4.2 旋度的定义
4.2.1 直观引入
考虑点 处垂直于方向 的小闭曲线 ,所围面积为 。向量场 沿 的环流量为:
单位面积的环流量:
当 时,这个比值的极限反映了点 处向量场绕方向 的"旋转强度"。
定义 4.2(旋度)
设向量场 在点 的某邻域内有定义,若存在一个向量 ,使得对任意单位向量 ,有:
则称 为向量场 在点 的旋度(curl 或 rotation),记作:
4.2.2 旋度的计算公式
定理 4.1(旋度的坐标表示)
若向量场 的分量函数有一阶连续偏导数,则:
$$\boxed{\text{curl}, \mathbf{F} = \nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ P & Q & R \end{vmatrix}}$$
展开得:
或分量形式:
记忆技巧:
将 形式上看作 Hamilton 算子 与向量 的叉积(向量积)。
4.3 旋度的物理意义
4.3.1 涡量
物理解释:
旋度 表示向量场 在某点的"旋转强度"或"涡量"(vorticity):
- :旋转的强度(角速度的两倍)
- 的方向:旋转轴的方向(右手螺旋)
- :无旋场(irrotational field)
4.3.2 应用实例
1. 流体力学
涡量:
**物理意义:**涡量是流体微团旋转角速度的两倍
无旋流动:(势流)
2. 电磁学
法拉第电磁感应定律(微分形式):
安培-麦克斯韦定律(微分形式):
3. 静电场
静电场是保守场:
因此:
**推论:**任何梯度场都是无旋场。
4.4 旋度的运算性质
定理 4.2(旋度的运算法则)
设 是向量场, 是标量场, 是常数,则:
-
线性性:
-
数乘法则:
-
梯度的旋度恒为零:
-
旋度的散度恒为零:
-
旋度的旋度:
其中
例题 7(原文例题)
**问题:**计算向量场 的旋度。
解:
$$\nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ y & -x & 0 \end{vmatrix}$$
物理意义:
这是绕 轴逆时针旋转的速度场(角速度 ),旋度为 ,指向 轴负方向(从俯视角度看,旋度向下表示逆时针旋转)。
例题 8
**问题:**证明 。
证明:
设 ,其中:
计算旋度的 分量:
(由混合偏导数相等)
类似地, 和 分量也为零。
因此 。
**物理意义:**梯度场(保守场)必为无旋场。
第五章:Hamilton 算子
5.1 Hamilton 算子的定义
定义 5.1(Hamilton 算子)
符号
称为 Hamilton 算子(Hamilton operator)或劈形算子(nabla operator)。
**历史:**由爱尔兰数学家 William Rowan Hamilton(1805-1865)引入, 的名称"nabla"源于古希腊的一种竖琴。
5.2 作用于标量场
5.2.1 梯度
性质:
- 输入:标量场
- 输出:向量场
5.2.2 拉普拉斯算子
记号:
应用:
-
拉普拉斯方程:(调和函数)
-
泊松方程:
-
热传导方程:
-
波动方程:
-
薛定谔方程:
5.3 作用于向量场
5.3.1 散度
性质:
- 输入:向量场
- 输出:标量场
5.3.2 旋度
$$\nabla \times \mathbf{F} = \text{curl}, \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ P & Q & R \end{vmatrix}$$
性质:
- 输入:向量场
- 输出:向量场
5.3.3 向量拉普拉斯
5.4 的重要恒等式
5.4.1 基本恒等式
定理 5.1( 运算恒等式)
设 是标量场, 是向量场,则:
-
梯度的旋度:
-
旋度的散度:
-
散度的梯度(不存在): 没有定义(标量无梯度)
但可以定义: 作为向量
-
旋度的旋度:
-
数乘的梯度:
-
数乘的散度:
-
数乘的旋度:
-
点积的散度:
-
叉积的旋度:
5.4.2 推论
推论 5.1
设 是调和函数(),则:
推论 5.2
若 (保守场),则:
推论 5.3
若 (无散场),则:
5.5 的运算规则总结
| 运算 | 输入 | 输出 | 记号 | 名称 |
|---|---|---|---|---|
| 标量 | 向量 | grad | 梯度 | |
| 向量 | 标量 | div | 散度 | |
| 向量 | 向量 | curl | 旋度 | |
| 标量 | 标量 | 拉普拉斯 | ||
| 向量 | 向量 | 向量拉普拉斯 |
第六章:保守场与势函数
6.1 保守场的定义
定义 6.1(保守场)
若向量场 可以表示为某个标量场 的梯度,即:
则称 为保守场(conservative field)或有势场, 称为 的势函数(potential function)。
**注:**物理中常用 ,其中 称为势能。
6.2 保守场的充要条件
定理 6.1(单连通区域中保守场的判别)
设 是单连通区域,向量场 在 内有一阶连续偏导数,则以下条件等价:
(i) 是保守场,即
(ii) 是无旋场,即
(iii) 对 内任意闭曲线 ,有
(iv) 在 内,曲线积分 与路径无关
无旋条件(显式):
6.3 势函数的求法
方法一:路径积分法
固定起点 ,定义:
选择简单路径(如折线)计算。
方法二:不定积分法
从 对 积分:
利用 和 确定 。
方法三:观察法(凑微分)
直接识别 是某个函数的全微分。
例题 9(原文例题)
**问题:**验证 是保守场,并求势函数。
解:
验证无旋:
所以 , 是保守场。
求势函数:
从 :
从 :
因此 ,即 。
从 :
因此 ,即 。
势函数:
6.4 物理中的保守场
6.4.1 引力场
**万有引力定律:
第六章:保守场与势函数(续)
6.4 物理中的保守场(续)
6.4.1 引力场
万有引力定律:
质量为 的质点在原点产生的引力场为:
其中 ,。
引力势能:
验证:
(利用例题 4 的结果)
**物理意义:**引力场是保守场,做功与路径无关,只与起点和终点位置有关。
6.4.2 静电场
库仑定律:
电荷为 的点电荷在原点产生的电场为:
其中 。
电势:
关系:
**推论:**静电场是保守场:
6.4.3 弹性力场
胡克定律:
弹性势能:
验证:
6.5 非保守场的例子
例 1:摩擦力场
摩擦力 始终与运动方向相反,做功依赖于路径长度,是非保守力。
沿闭路一周,。
例 2:涡旋场
旋度:
这是有旋场,不是保守场。
环流量:
沿单位圆 :
第七章:不同坐标系下的场论公式
7.1 柱坐标系
7.1.1 坐标变换
柱坐标 与直角坐标 的关系:
单位向量:
注意: 和 随位置变化!
7.1.2 梯度
7.1.3 散度
设 ,则:
7.1.4 旋度
7.1.5 拉普拉斯算子
7.2 球坐标系
7.2.1 坐标变换
球坐标 与直角坐标 的关系:
其中:
- :径向距离
- :方位角
- :极角(与 轴的夹角)
单位向量:
7.2.2 梯度
7.2.3 散度
7.2.4 旋度
7.2.5 拉普拉斯算子
例题 10
**问题:**在球坐标下,验证 满足拉普拉斯方程 ()。
解:
由于 不依赖于 和 :
因此 是调和函数()。
**注:**在 处,(狄拉克 函数)。
第八章:积分定理与场论的关系
8.1 高斯公式(散度定理)
定理 8.1(高斯公式)
设 是由分片光滑闭曲面 围成的空间区域,向量场 在 上有连续一阶偏导数,则:
其中 是 的外法向量。
微分形式:
物理意义:
区域内的总源强度(散度积分)= 穿过边界的总通量
推论:散度的物理意义验证
对于小区域 :
因此:
这正是散度的定义。
8.2 斯托克斯公式(旋度定理)
定理 8.2(斯托克斯公式)
设 是光滑曲面,其边界 是按段光滑的连续曲线,向量场 在 (含 )上有连续一阶偏导数,则:
其中 的方向与 的方向按右手法则确定。
微分形式:
物理意义:
穿过曲面的总涡通量(旋度通量)= 沿边界的环流量
推论:旋度的物理意义验证
对于小曲面 ,面积为 :
因此:
这正是旋度在法向 上的投影。
8.3 格林公式(平面情形)
定理 8.3(格林公式)
设 是平面区域,其边界 为按段光滑的闭曲线, 在 上有连续一阶偏导数,则:
其中 取逆时针方向。
与斯托克斯公式的关系:
格林公式是斯托克斯公式在平面 的特例。设 ,则:
8.4 三大积分定理的统一
| 定理 | 维度 | 关系 | 公式 |
|---|---|---|---|
| 牛顿-莱布尼茨 | 1D → 0D | 区间 → 端点 | |
| 格林公式 | 2D → 1D | 区域 → 边界 | |
| 斯托克斯公式 | 2D → 1D | 曲面 → 边界 | |
| 高斯公式 | 3D → 2D | 体积 → 表面 |
广义斯托克斯定理(微分形式版本):
其中 是微分形式, 是外微分算子, 是 的边界。
第九章:Helmholtz 分解定理
9.1 定理陈述
定理 9.1(Helmholtz 分解定理)
设 是定义在全空间 上的充分光滑的向量场,且在无穷远处充分快地趋于零,则 可以唯一分解为:
其中:
- :标量势(scalar potential)
- :向量势(vector potential)
术语:
- :无旋部分(irrotational part)
- :无散部分(solenoidal part)
9.2 构造方法
标量势:
向量势:
9.3 物理应用
9.3.1 电磁场
Maxwell 方程组:
势的引入:
由 ,引入向量势 :
由 :
因此引入标量势 :
即:
9.3.2 流体力学
速度场分解:
其中:
- :速度势(velocity potential)
- :流函数(stream function)
无旋流动(势流):,则
不可压缩流动:,则 (拉普拉斯方程)
第十章:经典场论问题
10.1 泊松方程与格林函数
10.1.1 泊松方程
方程:
特殊情形:
- :拉普拉斯方程
- :静电势的泊松方程
- :引力势的泊松方程
10.1.2 格林函数
**定义:**格林函数 满足:
加上适当的边界条件。
无界空间的格林函数:
泊松方程的解:
10.2 调和函数的性质
定义 10.1(调和函数)
满足拉普拉斯方程 的函数称为调和函数(harmonic function)。
定理 10.1(平均值定理)
若 在球 内调和,则:
即中心点的值等于球面上的平均值。
定理 10.2(最大值原理)
调和函数在有界区域内部不能取到最大值或最小值(除非为常数)。
**推论:**调和函数的极值必在边界上取得。
定理 10.3(唯一性)
Dirichlet 问题(给定边界值的拉普拉斯方程)的解唯一。
10.3 位势理论
10.3.1 牛顿位势
定义:
性质:
特殊情况:
-
点质量:,则
-
球对称分布: 其中 是总质量, 是球半径。
10.3.2 对数位势(二维)
定义:
性质:
第十一章:综合例题精讲
11.1 梯度与方向导数
例题 11
**问题:**设温度场为 ,求:
- 点 处温度变化最快的方向
- 该方向上的变化率
解:
1. 最快变化方向 = 梯度方向
在 :
单位化:
2. 最大变化率
物理意义:
- 温度在点 沿方向 下降最快
- 下降率为 °C/单位长度
- 这个方向指向原点(温度最高点)
11.2 散度与通量
例题 12
**问题:**计算向量场 通过单位球面 的通量。
解:
方法一:直接计算曲面积分(困难)
方法二:使用高斯公式
应用高斯公式:
使用球坐标:
$$= 3 \cdot 2\pi \cdot 2 \cdot \frac{1}{5}
第十一章:综合例题精讲(续)
11.2 散度与通量(续)
$$= 3 \cdot 2\pi \cdot 2 \cdot \frac{1}{5} = \frac{12\pi}{5}$$
答案:$\frac{12\pi}{5}$
**优势:**使用高斯公式将复杂的曲面积分转化为简单的三重积分,大大简化了计算。
$\square$
11.3 旋度与环流
例题 13
**问题:**计算向量场 $\mathbf{F} = (y - z, z - x, x - y)$ 沿圆周 $L: x^2 + y^2 = 1, z = 0$(逆时针)的环流量。
解:
方法一:直接计算曲线积分
参数化:$x = \cos t, y = \sin t, z = 0$,$0 \leq t \leq 2\pi$
$$\mathbf{F}|_L = (\sin t, -\cos t, \cos t - \sin t)d\mathbf{r} = (-\sin t, \cos t, 0) , dt\mathbf{F} \cdot d\mathbf{r} = (\sin t)(-\sin t) + (-\cos t)(\cos t) = -\sin^2 t - \cos^2 t = -1\Gamma = \oint_L \mathbf{F} \cdot d\mathbf{r} = \int_0^{2\pi} (-1) , dt = -2\pi$$
方法二:使用斯托克斯公式
取曲面 $S$ 为圆盘 $x^2 + y^2 \leq 1, z = 0$,法向量 $\mathbf{n} = (0, 0, 1)$(上侧)。
计算旋度:
$$\nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ y - z & z - x & x - y \end{vmatrix}= \mathbf{i}\left[\frac{\partial(x-y)}{\partial y} - \frac{\partial(z-x)}{\partial z}\right] - \mathbf{j}\left[\frac{\partial(x-y)}{\partial x} - \frac{\partial(y-z)}{\partial z}\right] + \mathbf{k}\left[\frac{\partial(z-x)}{\partial x} - \frac{\partial(y-z)}{\partial y}\right]= \mathbf{i}(-1 - 1) - \mathbf{j}(1 - (-1)) + \mathbf{k}(-1 - 1)= (-2, -2, -2)\Gamma = \iint_S (\nabla \times \mathbf{F}) \cdot \mathbf{n} , dS = \iint_S (-2, -2, -2) \cdot (0, 0, 1) , dS= \iint_S (-2) , dS = -2 \cdot \pi \cdot 1^2 = -2\pi$$
答案:$-2\pi$
**验证:**两种方法结果一致!✓
$\square$
11.4 保守场判断与势函数
例题 14
**问题:**判断向量场 $\mathbf{F} = (yz, xz, xy)$ 是否为保守场,若是,求势函数。
解:
第一步:检验无旋条件
$$\frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} = x - x = 0\frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} = y - y = 0\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = z - z = 0$$
因此 $\nabla \times \mathbf{F} = \mathbf{0}$,$\mathbf{F}$ 是保守场。
第二步:求势函数
设 $\phi$ 满足 $\nabla \phi = \mathbf{F}$,即:
$$\frac{\partial \phi}{\partial x} = yz, \quad \frac{\partial \phi}{\partial y} = xz, \quad \frac{\partial \phi}{\partial z} = xy$$
从第一式对 $x$ 积分:
$$\phi = \int yz , dx = xyz + f(y, z)$$
对 $y$ 求偏导:
$$\frac{\partial \phi}{\partial y} = xz + \frac{\partial f}{\partial y} = xz$$
因此 $\frac{\partial f}{\partial y} = 0$,即 $f = g(z)$。
$$\phi = xyz + g(z)$$
对 $z$ 求偏导:
$$\frac{\partial \phi}{\partial z} = xy + g'(z) = xy$$
因此 $g'(z) = 0$,即 $g(z) = C$。
势函数:
$$\phi(x, y, z) = xyz + C\nabla \phi = (yz, xz, xy) = \mathbf{F} \quad \checkmark$$
$\square$
11.5 拉普拉斯方程
例题 15
**问题:**证明函数 $u(x, y, z) = \frac{1}{\sqrt{x^2 + y^2 + z^2}}$ 在原点以外满足拉普拉斯方程 $\nabla^2 u = 0$。
证明:
设 $r = \sqrt{x^2 + y^2 + z^2}$,则 $u = \frac{1}{r}$。
方法一:直角坐标
$$\frac{\partial u}{\partial x} = -\frac{x}{r^3}\frac{\partial^2 u}{\partial x^2} = \frac{\partial}{\partial x}\left(-\frac{x}{r^3}\right) = -\frac{1}{r^3} + \frac{3x^2}{r^5}\frac{\partial^2 u}{\partial y^2} = -\frac{1}{r^3} + \frac{3y^2}{r^5}\frac{\partial^2 u}{\partial z^2} = -\frac{1}{r^3} + \frac{3z^2}{r^5}\nabla^2 u = -\frac{3}{r^3} + \frac{3(x^2 + y^2 + z^2)}{r^5} = -\frac{3}{r^3} + \frac{3r^2}{r^5} = -\frac{3}{r^3} + \frac{3}{r^3} = 0$$
方法二:球坐标
在球坐标中,对于只依赖于 $r$ 的函数:
$$\nabla^2 u = \frac{1}{r^2}\frac{d}{dr}\left(r^2 \frac{du}{dr}\right)u = \frac{1}{r}, \quad \frac{du}{dr} = -\frac{1}{r^2}r^2 \frac{du}{dr} = r^2 \cdot \left(-\frac{1}{r^2}\right) = -1\frac{d}{dr}(-1) = 0$$
因此 $\nabla^2 u = 0$。
**注:**在原点 $r = 0$ 处,$\nabla^2 u = -4\pi \delta(\mathbf{r})$(广义函数意义下)。
$\square$
11.6 场的分解
例题 16
**问题:**将向量场 $\mathbf{F} = (x + y, y + z, z + x)$ 分解为无旋部分和无散部分。
解:
根据 Helmholtz 分解定理:$\mathbf{F} = \nabla \phi + \nabla \times \mathbf{A}$
计算散度:
$$\nabla \cdot \mathbf{F} = \frac{\partial}{\partial x}(x + y) + \frac{\partial}{\partial y}(y + z) + \frac{\partial}{\partial z}(z + x) = 1 + 1 + 1 = 3\nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ x + y & y + z & z + x \end{vmatrix}= \mathbf{i}(0 - 1) - \mathbf{j}(1 - 0) + \mathbf{k}(0 - 1) = (-1, -1, -1)$$
无旋部分(纵向分量):
寻找 $\phi$ 使得 $\nabla^2 \phi = \nabla \cdot \mathbf{F} = 3$
最简单的选择:$\phi = \frac{3}{2}(x^2 + y^2 + z^2)$
验证:$\nabla^2 \phi = 3 + 3 + 3 = 9 \neq 3$(错误)
正确的:$\phi = \frac{3}{2}r^2$ 满足 $\nabla^2 \phi = 3 \cdot 2 = 6 \neq 3$
实际上,$\phi = \frac{1}{2}(x^2 + y^2 + z^2)$ 给出 $\nabla^2 \phi = 3$ ✓
$$\nabla \phi = (x, y, z)\mathbf{F}_{\text{solenoidal}} = \mathbf{F} - \nabla \phi = (x + y, y + z, z + x) - (x, y, z) = (y, z, x)\nabla \cdot (y, z, x) = 0 + 0 + 0 = 0 \quad \checkmark\nabla \times (y, z, x) = (0 - 1, 1 - 0, 0 - 1) = (-1, 1, -1)$$
**注:**这个场实际上不能完全用简单形式分解,需要更仔细的构造向量势 $\mathbf{A}$。
重新分析:
直接观察:
$$\mathbf{F} = (x + y, y + z, z + x) = (x, y, z) + (y, z, x)$$
- $(x, y, z) = \nabla\left(\frac{x^2 + y^2 + z^2}{2}\right)$:无旋部分
- $(y, z, x)$:需要找 $\mathbf{A}$ 使得 $\nabla \times \mathbf{A} = (y, z, x)$
尝试 $\mathbf{A} = (0, xz, -xy)$:
$$\nabla \times \mathbf{A} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ 0 & xz & -xy \end{vmatrix}= \mathbf{i}(-x - 0) - \mathbf{j}(0 - 0) + \mathbf{k}(z - 0) = (-x, 0, z) \neq (y, z, x)$$
(此问题较复杂,实际需要解偏微分方程系统)
$\square$
第十二章:物理应用专题
12.1 静电学中的场论
12.1.1 电场与电势
基本关系:
$$\mathbf{E} = -\nabla \phi\oiint_S \mathbf{E} \cdot \mathbf{n} , dS = \frac{Q_{\text{enc}}}{\varepsilon_0}\nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0}$$
结合 $\mathbf{E} = -\nabla \phi$:
$$\nabla \cdot (-\nabla \phi) = \frac{\rho}{\varepsilon_0}\nabla^2 \phi = -\frac{\rho}{\varepsilon_0}$$
这是泊松方程。
无电荷区域($\rho = 0$):
$$\nabla^2 \phi = 0$$
这是拉普拉斯方程,$\phi$ 是调和函数。
12.1.2 电偶极子
电偶极矩:$\mathbf{p} = q\mathbf{d}$(正电荷到负电荷的向量,乘以电荷量)
电势(远场近似):
$$\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \frac{\mathbf{p} \cdot \mathbf{r}}{r^3} = \frac{1}{4\pi\varepsilon_0} \frac{\mathbf{p} \cdot \hat{\mathbf{r}}}{r^2}\mathbf{E} = -\nabla \phi = \frac{1}{4\pi\varepsilon_0} \left[\frac{3(\mathbf{p} \cdot \hat{\mathbf{r}})\hat{\mathbf{r}} - \mathbf{p}}{r^3}\right]$$
12.2 磁场与安培定律
12.2.1 毕奥-萨伐尔定律
电流元产生的磁场:
$$d\mathbf{B} = \frac{\mu_0}{4\pi} \frac{I , d\mathbf{l} \times \mathbf{r}}{r^3}\mathbf{B}(\mathbf{r}) = \frac{\mu_0 I}{4\pi} \oint \frac{d\mathbf{l}' \times (\mathbf{r} - \mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|^3}$$
12.2.2 安培定律
积分形式:
$$\oint_L \mathbf{B} \cdot d\mathbf{l} = \mu_0 I_{\text{enc}}\nabla \times \mathbf{B} = \mu_0 \mathbf{J}\nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}}{\partial t}$$
12.2.3 磁场的性质
无散性:
$$\nabla \cdot \mathbf{B} = 0\mathbf{B} = \nabla \times \mathbf{A}$$
其中 $\mathbf{A}$ 是向量势(vector potential)。
12.3 流体力学中的场论
12.3.1 连续性方程
质量守恒:
$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$
不可压缩流体($\rho =$ 常数):
$$\nabla \cdot \mathbf{v} = 0$$
12.3.2 涡量动力学
涡量:$\boldsymbol{\omega} = \nabla \times \mathbf{v}$
涡量方程(理想流体):
$$\frac{D\boldsymbol{\omega}}{Dt} = (\boldsymbol{\omega} \cdot \nabla)\mathbf{v}$$
其中 $\frac{D}{Dt} = \frac{\partial}{\partial t} + \mathbf{v} \cdot \nabla$ 是物质导数。
12.3.3 势流理论
无旋流动:$\nabla \times \mathbf{v} = \mathbf{0}$
引入速度势 $\Phi$:
$$\mathbf{v} = \nabla \Phi\nabla \cdot \mathbf{v} = \nabla \cdot (\nabla \Phi) = \nabla^2 \Phi = 0$$
速度势满足拉普拉斯方程。
12.4 热传导
12.4.1 傅里叶定律
热流密度:
$$\mathbf{q} = -k \nabla T$$
其中 $k$ 是热导率。
**物理意义:**热量从高温流向低温,流速与温度梯度成正比。
12.4.2 热传导方程
能量守恒:
$$\rho c \frac{\partial T}{\partial t} = -\nabla \cdot \mathbf{q}\rho c \frac{\partial T}{\partial t} = -\nabla \cdot (-k\nabla T) = k \nabla^2 T\frac{\partial T}{\partial t} = \alpha \nabla^2 T$$
其中 $\alpha = \frac{k}{\rho c}$ 是热扩散系数(thermal diffusivity)。
12.4.3 稳态热传导
稳态条件:$\frac{\partial T}{\partial t} = 0$
$$\nabla^2 T = 0$$
温度分布满足拉普拉斯方程,是调和函数。
第十三章:场论与微分形式
13.1 微分形式简介
13.1.1 外代数
0-形式(标量):$f$
1-形式(余切向量):$\omega^1 = P , dx + Q , dy + R , dz$
2-形式(面元):$\omega^2 = A , dy \wedge dz + B , dz \wedge dx + C , dx \wedge dy$
3-形式(体元):$\omega^3 = F , dx \wedge dy \wedge dz$
外积(wedge product):
$$dx \wedge dy = -dy \wedge dxdx \wedge dx = 0$$
**性质:**反对称、反交换
13.1.2 外微分算子
定义:$d: \Omega^k \to \Omega^{k+1}$
0-形式的外微分:
$$df = \frac{\partial f}{\partial x} dx + \frac{\partial f}{\partial y} dy + \frac{\partial f}{\partial z} dzd(P , dx + Q , dy + R , dz) = dP \wedge dx + dQ \wedge dy + dR \wedge dz= \left(\frac{\partial P}{\partial y} dy + \frac{\partial P}{\partial z} dz\right) \wedge dx + \left(\frac{\partial Q}{\partial x} dx + \frac{\partial Q}{\partial z} dz\right) \wedge dy + \left(\frac{\partial R}{\partial x} dx + \frac{\partial R}{\partial y} dy\right) \wedge dz= \left(\frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z}\right) dy \wedge dz + \left(\frac{\partial P}{\partial z} - \frac{\partial R}{\partial x}\right) dz \wedge dx + \left(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right) dx \wedge dyd(A , dy \wedge dz + B , dz \wedge dx + C , dx \wedge dy)= \left(\frac{\partial A}{\partial x} + \frac{\partial B}{\partial y} + \frac{\partial C}{\partial z}\right) dx \wedge dy \wedge dz$$
(对应散度)
13.1.3 基本性质
定理 13.1
-
$d^2 = 0$:$d(d\omega) = 0$ 对任意 $k$-形式 $\omega$ 成立
-
Leibniz 法则:$d(\omega \wedge \eta) = d\omega \wedge \eta + (-1)^k \omega \wedge d\eta$
其中 $\omega$ 是 $k$-形式
物理对应:
- $d^2 = 0$ 对应于:
- $\nabla \times (\nabla f) = \mathbf{0}$(梯度无旋)
- $\nabla \cdot (\nabla \times \mathbf{F}) = 0$(旋度无散)
13.2 广义斯托克斯定理
定理 13.2(广义斯托克斯定理)
设 $\Omega$ 是 $n$ 维带边流形,$\omega$ 是 $\Omega$ 上的 $(n-1)$-形式,则:
$$\boxed{\int_\Omega d\omega = \int_{\partial \Omega} \omega}$$
其中 $\partial \Omega$ 是 $\Omega$ 的边界,定向由右手法则确定。
特殊情况:
| $n$ | $\Omega$ | 定理 | 形式 |
|---|---|---|---|
| 1 | 区间 $[a, b]$ | 微积分基本定理 | $\int_a^b f'(x) , dx = f(b) - f(a)$ |
| 2 | 平面区域 $D$ | 格林公式 | $\iint_D d\omega = \oint_{\partial D} \omega$ |
| 2 | 曲面 $S$ | 斯托克斯公式 | $\iint_S d\omega = \oint_{\partial S} \omega$ |
| 3 | 立体 $V$ | 高斯公式 | $\iiint_V d\omega = \oiint_{\partial V} \omega$ |
13.3 德拉姆上同调
13.3.1 定义
闭形式(closed form):$d\omega = 0$
恰当形式(exact form):$\omega = d\eta$ 对某个 $\eta$
**关系:**恰当 $\Rightarrow$ 闭(因为 $d^2 = 0$)
**反之:**闭 $\not\Rightarrow$ 恰当(在非单连通区域)
德拉姆上同调群:
$$H^k_{dR}(M) = \frac{{\text{闭 } k\text{-形式}}}{{\text{恰当 } k\text{-形式}}} = \frac{\ker d}{\text{im } d}$$
物理意义:$H^k_{dR}(M)$ 测量空间 $M$ 的 "$k$ 维洞"的数量。
13.3.2 例子
1. $\mathbb{R}^3$(单连通):
$$H^0_{dR}(\mathbb{R}^3) = \mathbb{R}, \quad H^1_{dR}(\mathbb{R}^3) = 0, \quad H^2_{dR}(\mathbb{R}^3) = 0$$
**物理意义:**所有无旋场都是保守场(可表示为梯度)。
2. $\mathbb{R}^3 \setminus {0}$(去掉原点,非单连通):
$$H^1_{dR}(\mathbb{R}^3 \setminus {0}) \neq 0$$
**物理意义:**存在闭但不恰当的 1-形式,对应于围绕原点有非零环流的无旋场。
例子:$\omega = \frac{-y , dx + x , dy}{x^2 + y^2}$
- $d\omega = 0$(闭)
- $\oint_C \omega = 2\pi \neq 0$(不恰当)
其中 $C$ 是绕原点一周的闭曲线。
第十四章:高级专题
14.1 Green 恒等式
14.1.1 第一 Green 恒等式
定理 14.1
$$\iiint_V (\phi \nabla^2 \psi + \nabla \phi \cdot \nabla \psi) , dV = \oiint_S \phi \frac{\partial \psi}{\partial n} , dS$$
证明:
考虑向量场 $\mathbf{F} = \phi \nabla \psi$
$$\nabla \cdot \mathbf{F} = \nabla \cdot (\phi \nabla \psi) = \phi \nabla^2 \psi + \nabla \phi \cdot \nabla \psi\iiint_V \nabla \cdot \mathbf{F} , dV = \oiint_S \mathbf{F} \cdot \mathbf{n} , dS = \oiint_S \phi \nabla \psi \cdot \mathbf{n} , dS = \oiint_S \phi \frac{\partial \psi}{\partial n} , dS$$
$\square$
14.1.2 第二 Green 恒等式
定理 14.2
$$\iiint_V (\phi \nabla^2 \psi - \psi \nabla^2 \phi) , dV = \oiint_S \left(\phi \frac{\partial \psi}{\partial n} - \psi \frac{\partial \phi}{\partial n}\right) dS$$
证明:
应用第一 Green 恒等式两次,交换 $\phi$ 和 $\psi$,然后相减。
$\square$
应用:
- 证明调和函数的唯一性
- 推导 Green 函数
- 证明平均值定理
14.2 向量场的分类
14.2.1 Helmholtz 分类
任何充分光滑且在无穷远衰减的向量场 $\mathbf{F}$ 可以唯一分解为:
$$\mathbf{F} = \mathbf{F}_L + \mathbf{F}_T$$
其中:
- $\mathbf{F}_L = \nabla \phi$:纵向场(longitudinal field),无旋($\nabla \times \mathbf{F}_L = \mathbf{0}$)
- $\mathbf{F}_T = \nabla \times \mathbf{A}$:横向场(transverse field),无散($\nabla \cdot \mathbf{F}_T = 0$)
14.2.2 判别方法
| 场类型 | 条件 | 名称 | 例子 |
|---|---|---|---|
| 无旋场 | $\nabla \times \mathbf{F} = \mathbf{0}$ | 保守场、有势场 | 静电场 $\mathbf{E}$ |
| 无散 |
第十四章:高级专题(续)
14.2 向量场的分类(续)
14.2.2 判别方法(续)
| 场类型 | 条件 | 名称 | 例子 |
|---|---|---|---|
| 无旋场 | $\nabla \times \mathbf{F} = \mathbf{0}$ | 保守场、有势场 | 静电场 $\mathbf{E}$ |
| 无散场 | $\nabla \cdot \mathbf{F} = 0$ | 螺线场、管形场 | 磁场 $\mathbf{B}$ |
| 调和场 | $\nabla \times \mathbf{F} = \mathbf{0}$ 且 $\nabla \cdot \mathbf{F} = 0$ | 无源无旋场 | 真空中的静电场 |
| 拉普拉斯场 | $\mathbf{F} = \nabla \phi$,$\nabla^2 \phi = 0$ | 调和势场 | 无电荷区域的电场 |
14.3 Hodge 分解定理
定理 14.3(Hodge 分解定理)
设 $M$ 是紧致定向 Riemannian 流形,则任何光滑微分形式 $\omega$ 可以唯一分解为:
$$\omega = d\alpha + \delta\beta + \gamma$$
其中:
- $d\alpha$:恰当部分(exact part)
- $\delta\beta$:余恰当部分(coexact part),$\delta$ 是 $d$ 的伴随算子
- $\gamma$:调和部分(harmonic part),$d\gamma = 0$ 且 $\delta\gamma = 0$
三维向量场的特殊情况:
$$\mathbf{F} = \nabla \phi + \nabla \times \mathbf{A} + \mathbf{H}$$
其中 $\mathbf{H}$ 是调和向量场($\nabla \cdot \mathbf{H} = 0$,$\nabla \times \mathbf{H} = \mathbf{0}$)。
14.4 张量场论初步
14.4.1 张量的定义
二阶张量:
$$T_{ij}, \quad i, j = 1, 2, 3$$
对称张量:$T_{ij} = T_{ji}$
反对称张量:$T_{ij} = -T_{ji}$
14.4.2 张量的散度
定义:
$$(\nabla \cdot \mathbf{T})i = \sum{j=1}^3 \frac{\partial T_{ij}}{\partial x_j}\nabla \cdot \boldsymbol{\sigma} = \mathbf{f}$$
其中 $\boldsymbol{\sigma}$ 是应力张量,$\mathbf{f}$ 是体力密度。
14.4.3 应变张量
位移场:$\mathbf{u} = (u_1, u_2, u_3)$
应变张量:
$$\varepsilon_{ij} = \frac{1}{2}\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right)\theta = \varepsilon_{11} + \varepsilon_{22} + \varepsilon_{33} = \nabla \cdot \mathbf{u}$$
14.5 广义坐标下的场论
14.5.1 度规张量
在广义坐标 $(q^1, q^2, q^3)$ 中,度规张量:
$$g_{ij} = \frac{\partial \mathbf{r}}{\partial q^i} \cdot \frac{\partial \mathbf{r}}{\partial q^j}$$
标度因子:$h_i = \sqrt{g_{ii}}$(无求和)
14.5.2 正交曲线坐标
梯度:
$$\nabla \phi = \sum_{i=1}^3 \frac{1}{h_i} \frac{\partial \phi}{\partial q^i} \hat{\mathbf{e}}i\nabla \cdot \mathbf{F} = \frac{1}{h_1 h_2 h_3} \sum{i=1}^3 \frac{\partial}{\partial q^i}(h_1 h_2 h_3 \frac{F_i}{h_i})$$
(其中 $h_1 h_2 h_3 / h_i$ 表示去掉 $h_i$ 的乘积)
旋度:
$$\nabla \times \mathbf{F} = \frac{1}{h_1 h_2 h_3} \begin{vmatrix} h_1 \hat{\mathbf{e}}_1 & h_2 \hat{\mathbf{e}}_2 & h_3 \hat{\mathbf{e}}3 \ \frac{\partial}{\partial q^1} & \frac{\partial}{\partial q^2} & \frac{\partial}{\partial q^3} \ h_1 F_1 & h_2 F_2 & h_3 F_3 \end{vmatrix}\nabla^2 \phi = \frac{1}{h_1 h_2 h_3} \sum{i=1}^3 \frac{\partial}{\partial q^i}\left(\frac{h_1 h_2 h_3}{h_i^2} \frac{\partial \phi}{\partial q^i}\right)$$
14.5.3 常见坐标系的标度因子
柱坐标 $(r, \theta, z)$:
$$h_r = 1, \quad h_\theta = r, \quad h_z = 1$$
球坐标 $(r, \theta, \varphi)$:
$$h_r = 1, \quad h_\theta = r\sin\varphi, \quad h_\varphi = rh_u = h_v = a\sqrt{\sinh^2 u + \sin^2 v}, \quad h_z = 1$$
第十五章:场论在数学物理方程中的应用
15.1 波动方程
15.1.1 标量波动方程
$$\frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u$$
物理例子:
- 弦振动:$u(x, t)$ 是位移
- 声波:$u(\mathbf{r}, t)$ 是压强扰动
- 电磁波(标量近似)
15.1.2 向量波动方程
$$\frac{\partial^2 \mathbf{E}}{\partial t^2} = c^2 \nabla^2 \mathbf{E}\nabla^2 \mathbf{E} - \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} = \mathbf{0}\nabla^2 \mathbf{B} - \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{B}}{\partial t^2} = \mathbf{0}$$
其中 $c = \frac{1}{\sqrt{\mu_0 \varepsilon_0}}$ 是光速。
15.1.3 达朗贝尔算子
定义:
$$\Box = \nabla^2 - \frac{1}{c^2}\frac{\partial^2}{\partial t^2}\Box u = 0$$
15.2 热传导方程
15.2.1 抛物型方程
$$\frac{\partial u}{\partial t} = \alpha \nabla^2 u$$
性质:
- 不可逆(时间箭头)
- 无限传播速度
- 平滑解(瞬间平滑化初始不连续)
15.2.2 基本解(热核)
三维热核:
$$G(x, y, z, t) = \frac{1}{(4\pi \alpha t)^{3/2}} e^{-\frac{x^2 + y^2 + z^2}{4\alpha t}}u(\mathbf{r}, t) = \iiint_{\mathbb{R}^3} G(\mathbf{r} - \mathbf{r}', t) u_0(\mathbf{r}') , d\mathbf{r}'$$
其中 $u_0(\mathbf{r})$ 是初始温度分布。
15.3 Laplace 方程与调和函数
15.3.1 边值问题类型
1. Dirichlet 问题:
$$\begin{cases} \nabla^2 u = 0 & \text{在 } \Omega \text{ 内} \ u = f & \text{在 } \partial\Omega \text{ 上} \end{cases}\begin{cases} \nabla^2 u = 0 & \text{在 } \Omega \text{ 内} \ \frac{\partial u}{\partial n} = g & \text{在 } \partial\Omega \text{ 上} \end{cases}\alpha u + \beta \frac{\partial u}{\partial n} = h$$
15.3.2 分离变量法
例:矩形区域上的 Laplace 方程
$$\nabla^2 u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0, \quad 0 < x < a, 0 < y < bu(0, y) = u(a, y) = u(x, 0) = 0, \quad u(x, b) = f(x)$$
分离变量:$u(x, y) = X(x)Y(y)$
$$\frac{X''}{X} + \frac{Y''}{Y} = 0\frac{X''}{X} = -\lambda, \quad \frac{Y''}{Y} = \lambdaX'' + \lambda X = 0, \quad X(0) = X(a) = 0$$
本征值:$\lambda_n = \left(\frac{n\pi}{a}\right)^2$,$n = 1, 2, 3, \ldots$
本征函数:$X_n(x) = \sin\frac{n\pi x}{a}$
对应的 $Y$:
$$Y_n'' - \left(\frac{n\pi}{a}\right)^2 Y_n = 0Y_n(y) = \sinh\frac{n\pi y}{a}$$
(利用 $Y_n(0) = 0$)
通解:
$$u(x, y) = \sum_{n=1}^\infty A_n \sin\frac{n\pi x}{a} \sinh\frac{n\pi y}{a}$$
利用 $u(x, b) = f(x)$ 确定 $A_n$(Fourier 级数)。
15.4 Helmholtz 方程
15.4.1 方程形式
$$\nabla^2 u + k^2 u = 0$$
其中 $k$ 是波数(wave number)。
**来源:**时谐波动方程的空间部分($u(\mathbf{r}, t) = \phi(\mathbf{r})e^{-i\omega t}$)
15.4.2 球对称解
设 $u = u(r)$(球对称),则:
$$\frac{1}{r^2}\frac{d}{dr}\left(r^2 \frac{du}{dr}\right) + k^2 u = 0$$
令 $u = \frac{v(r)}{r}$:
$$v'' + k^2 v = 0v(r) = A\cos(kr) + B\sin(kr)u(r) = \frac{A\cos(kr) + B\sin(kr)}{r}u(r) = \frac{C e^{ikr} + D e^{-ikr}}{r}$$
**物理意义:**球面波(外传或内收)
15.5 Maxwell 方程组的完整形式
15.5.1 微分形式
$$\begin{cases} \nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0} & \text{(高斯定律)} \ \nabla \cdot \mathbf{B} = 0 & \text{(无磁单极)} \ \nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t} & \text{(法拉第定律)} \ \nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0\varepsilon_0 \frac{\partial \mathbf{E}}{\partial t} & \text{(安培-麦克斯韦定律)} \end{cases}$$
15.5.2 势的表示
引入标量势 $\phi$ 和向量势 $\mathbf{A}$:
$$\mathbf{B} = \nabla \times \mathbf{A}\mathbf{E} = -\nabla \phi - \frac{\partial \mathbf{A}}{\partial t}\mathbf{A}' = \mathbf{A} + \nabla \chi, \quad \phi' = \phi - \frac{\partial \chi}{\partial t}$$
不改变 $\mathbf{E}$ 和 $\mathbf{B}$。
15.5.3 Lorenz 规范
选择:
$$\nabla \cdot \mathbf{A} + \mu_0\varepsilon_0 \frac{\partial \phi}{\partial t} = 0\Box \phi = -\frac{\rho}{\varepsilon_0}\Box \mathbf{A} = -\mu_0 \mathbf{J}$$
其中 $\Box = \nabla^2 - \frac{1}{c^2}\frac{\partial^2}{\partial t^2}$ 是达朗贝尔算子。
第十六章:计算技巧与实用公式
16.1 常用恒等式汇总
16.1.1 一级运算
$$\nabla(c) = \mathbf{0}\nabla \cdot \mathbf{c} = 0 \quad (\mathbf{c} \text{ 是常向量})\nabla \times \mathbf{c} = \mathbf{0}$$
16.1.2 二级运算
$$\nabla \times (\nabla \phi) = \mathbf{0}\nabla \cdot (\nabla \times \mathbf{F}) = 0\nabla \times (\nabla \times \mathbf{F}) = \nabla(\nabla \cdot \mathbf{F}) - \nabla^2 \mathbf{F}\nabla \cdot (\nabla \phi) = \nabla^2 \phi$$
16.1.3 乘积法则
$$\nabla(\phi \psi) = \phi \nabla \psi + \psi \nabla \phi\nabla \cdot (\phi \mathbf{F}) = \phi(\nabla \cdot \mathbf{F}) + \mathbf{F} \cdot \nabla \phi\nabla \times (\phi \mathbf{F}) = \phi(\nabla \times \mathbf{F}) + (\nabla \phi) \times \mathbf{F}\nabla(\mathbf{F} \cdot \mathbf{G}) = (\mathbf{F} \cdot \nabla)\mathbf{G} + (\mathbf{G} \cdot \nabla)\mathbf{F} + \mathbf{F} \times (\nabla \times \mathbf{G}) + \mathbf{G} \times (\nabla \times \mathbf{F})\nabla \cdot (\mathbf{F} \times \mathbf{G}) = \mathbf{G} \cdot (\nabla \times \mathbf{F}) - \mathbf{F} \cdot (\nabla \times \mathbf{G})\nabla \times (\mathbf{F} \times \mathbf{G}) = \mathbf{F}(\nabla \cdot \mathbf{G}) - \mathbf{G}(\nabla \cdot \mathbf{F}) + (\mathbf{G} \cdot \nabla)\mathbf{F} - (\mathbf{F} \cdot \nabla)\mathbf{G}$$
16.2 常见函数的场量
16.2.1 径向函数 $f(r)$
其中 $r = |\mathbf{r}| = \sqrt{x^2 + y^2 + z^2}$,$\hat{\mathbf{r}} = \frac{\mathbf{r}}{r}$
梯度:
$$\nabla f(r) = f'(r) \hat{\mathbf{r}}\nabla^2 f(r) = f''(r) + \frac{2}{r}f'(r) = \frac{1}{r^2}\frac{d}{dr}\left(r^2 \frac{df}{dr}\right)\nabla\left(\frac{1}{r}\right) = -\frac{\hat{\mathbf{r}}}{r^2} = -\frac{\mathbf{r}}{r^3}\nabla^2\left(\frac{1}{r}\right) = 0 \quad (r \neq 0)\nabla(r^n) = nr^{n-2}\mathbf{r}\nabla^2(r^n) = n(n+1)r^{n-2}$$
16.2.2 径向向量场
$$\mathbf{F} = f(r)\hat{\mathbf{r}} = \frac{f(r)}{r}\mathbf{r}\nabla \cdot \mathbf{F} = \frac{1}{r^2}\frac{d}{dr}(r^2 f(r)) = 2f(r) + rf'(r)\nabla \times \mathbf{F} = \mathbf{0}\mathbf{F} = \frac{\mathbf{r}}{r^3}\nabla \cdot \mathbf{F} = \frac{1}{r^2}\frac{d}{dr}\left(r^2 \cdot \frac{1}{r^3}\right) = \frac{1}{r^2}\frac{d}{dr}\left(\frac{1}{r}\right) = \frac{1}{r^2} \cdot \left(-\frac{1}{r^2}\right) \cdot (-1) = 0 \quad (r \neq 0)\nabla \times \mathbf{F} = \mathbf{0}$$
16.3 积分计算技巧
16.3.1 利用对称性
**例:**计算 $\iiint_V x^2 , dV$,其中 $V$ 是单位球。
由对称性:
$$\iiint_V x^2 , dV = \iiint_V y^2 , dV = \iiint_V z^2 , dV\iiint_V (x^2 + y^2 + z^2) , dV = \iiint_V r^2 , dV= \int_0^{2\pi} \int_0^\pi \int_0^1 r^2 \cdot r^2 \sin\varphi , dr d\varphi d\theta = 2\pi \cdot 2 \cdot \frac{1}{5} = \frac{4\pi}{5}\iiint_V x^2 , dV = \frac{1}{3} \cdot \frac{4\pi}{5} = \frac{4\pi}{15}$$
16.3.2 利用散度定理简化
**例:**计算 $\oiint_S \mathbf{r} \cdot \mathbf{n} , dS$,其中 $S$ 是任意闭曲面,包围体积 $V$。
$$\nabla \cdot \mathbf{r} = \nabla \cdot (x, y, z) = 1 + 1 + 1 = 3\oiint_S \mathbf{r} \cdot \mathbf{n} , dS = \iiint_V 3 , dV = 3V$$
**结果:**等于体积的三倍。
16.3.3 利用 Stokes 定理简化
**例:**计算 $\oint_L (y , dx + z , dy + x , dz)$,其中 $L$ 是单位圆 $x^2 + y^2 = 1, z = 0$(逆时针)。
设 $\mathbf{F} = (y, z, x)$
$$\nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ y & z & x \end{vmatrix} = (-1, -1, -1)$$
取曲面 $S$ 为圆盘 $x^2 + y^2 \leq 1, z = 0$,$\mathbf{n} = (0, 0, 1)$
$$\oint_L \mathbf{F} \cdot d\mathbf{r} = \iint_S (\nabla \times \mathbf{F}) \cdot \mathbf{n} , dS = \iint_S (-1, -1, -1) \cdot (0, 0, 1) , dS= \iint_S (-1) , dS = -\pi$$
16.4 数值方法简介
16.4.1 有限差分法
梯度的差分近似:
$$\frac{\partial u}{\partial x}\bigg|{(x_i, y_j, z_k)} \approx \frac{u{i+1,j,k} - u_{i-1,j,k}}{2\Delta x}\nabla^2 u\big|{(x_i, y_j)} \approx \frac{u{i+1,j} + u_{i-1,j} + u_{i,j+1} + u_{i,j-1} - 4u_{i,j}}{\Delta x^2}$$
16.4.2 有限元法
将区域离散化为有限个单元(三角形、四面体等),在每个单元上用简单函数(多项式)近似解,通过变分原理得到离散方程组。
**优点:**适用于复杂几何形状和边界条件。
16.4.3 谱方法
用正交函数系(Fourier 级数、Legendre 多项式等)展开解,将偏微分方程转化为代数方程。
**优点:**高精度(指数收敛)。
**缺点:**对复杂边界不友好。
第十七章:习题精选与解答
17.1 基础题
习题 1
**问题:**求函数 $f(x, y, z) = x^2y + y^2z + z^2x$ 在点 $(1, 1, 1)$ 处沿方向 $(1, 1, 1)$ 的方向导数。
解答:
$$\nabla f = (2xy + z^2, x^2 + 2yz, y^2 + 2zx)$$
在 $(1, 1, 1)$:
$$\nabla f\big|_{(1,1,1)} = (2 + 1, 1 + 2, 1 + 2) = (3, 3, 3)\mathbf{e} = \frac{1}{\sqrt{3}}(1, 1, 1)\frac{\partial f}{\partial \mathbf{l}} = \nabla f \cdot \mathbf{e} = (3, 3, 3) \cdot \frac{1}{\sqrt{3}}(1, 1, 1) = \frac{9}{\sqrt{3}} = 3\sqrt{3}$$
$\square$
习题 2
**问题:**证明向量场 $\mathbf{F} = (2xy + z, x^2 + 2yz, y^2 + x)$ 是保守场,并求势函数。
解答:
检验无旋:
$$\frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} = 2y - 2y = 0\frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} = 1 - 1 = 0\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = 2x - 2x = 0\frac{\partial \phi}{\partial x} = 2xy + z \Rightarrow \phi = x^2y + xz + f(y, z)\frac{\partial \phi}{\partial y} = x^2 + \frac{\partial f}{\partial y} = x^2 + 2yz \Rightarrow \frac{\partial f}{\partial y} = 2yz \Rightarrow f = y^2z + g(z)\phi = x^2y + xz + y^2z + g(z)\frac{\partial \phi}{\partial z} = x + y^2 + g'(z) = y^2 + x \Rightarrow g'(z) = 0 \Rightarrow g = C\phi(x, y, z) = x^2y + xz + y^2z + C$$
$\square$
17.2 中等题
习题 3
**问题:**设 $\mathbf{F} = \nabla \times \mathbf{A}$,其中 $\mathbf{A} = (xy^2, 2x^2y, 0)$。求 $\mathbf{F}$ 通过上半球面 $x^2 + y^2 + z^2 = 1, z \geq 0$ 的通量。
解答:
由于 $\mathbf{F} = \nabla \times \mathbf{A}$,有:
$$\nabla \cdot \mathbf{F} = \nabla \cdot (\nabla \times \mathbf{A}) = 0$$
设 $S_1$ 是上半球面,$S_2$ 是圆盘 $x^2 + y^2 \leq 1, z = 0$(向下),它们组成闭曲面 $S = S_1 \cup S_2$。
由高斯公式:
$$\oiint_S \mathbf{F} \cdot \mathbf{n} , dS = \iiint_V \nabla \cdot \mathbf{F} , dV = 0\iint_{S_1} \mathbf{F} \cdot \mathbf{n}1 , dS + \iint{S_2} \mathbf{F} \cdot \mathbf{n}_2 , dS = 0$$
其中 $\mathbf{n}_1$ 是向外(向上),$\mathbf{n}_2 = (0, 0, -1)$。
$$\iint_{S_1} \mathbf{F} \cdot \mathbf{n}1 , dS = -\iint{S_2} \mathbf{F} \cdot (0, 0, -1) , dS = \iint_{S_2} F_z , dS$$
计算 $\mathbf{F}$:
$$\mathbf{F} = \nabla \times \mathbf{A} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \ xy^2 & 2x^2y & 0 \end{vmatrix}= \mathbf{i}(0 - 0) - \mathbf{j}(0 - 0) + \mathbf{k}(4xy - 2xy) = (0, 0, 2xy)$$
在 $S_2$($z = 0$ 上):
$$\iint_{S_2} F_z , dS = \iint_{x^2+y^2 \leq 1} 2xy , dx dy$$
由对称性(积分区域关于 $x$ 轴和 $y$ 轴对称,被积函数是奇函数):
$$\iint_{x^2+y^2 \leq 1} 2xy , dx dy = 0$$
**答案:**通量为 $\boxed{0}$。
$\square$
习题 4
**问题:**证明:若 $\mathbf{F}$ 和 $\mathbf{G}$ 都是调和向量场(即 $\nabla^2 \mathbf{F} = \mathbf{0}$,$\nabla^2 \mathbf{G} = \mathbf{0}$),则 $\mathbf{F} \cdot \mathbf{G}$ 满足:
$$\nabla^2(\mathbf{F} \cdot \mathbf{G}) = 2[(\nabla \mathbf{F}) : (\nabla \mathbf{G})] + 2\mathbf{F} \cdot \nabla^2 \mathbf{G} + 2\mathbf{G} \cdot \nabla^2 \mathbf{F}\nabla(\mathbf{F} \cdot \mathbf{G}) = (\mathbf{F} \cdot \nabla)\mathbf{G} + (\mathbf{G} \cdot \nabla)\mathbf{F} + \mathbf{F} \times (\nabla \times \mathbf{G}) + \mathbf{G} \times (\nabla \times \mathbf{F})$$
再取散度,过程较复杂,直接给出结果:
在调和条件 $\nabla^2 \mathbf{F} = \mathbf{0}$,$\nabla^2 \mathbf{G} = \mathbf{0}$ 下:
$$\nabla^2(\mathbf{F} \cdot \mathbf{G}) = 2\sum_{i,j} \frac{\partial F_i}{\partial x_j} \frac{\partial G_i}{\partial x_j}$$
(这是张量缩并 $(\nabla \mathbf{F}) : (\nabla \mathbf{G})$)
$\square$
17.3 挑战题
习题 5
**问题:**设 $u$ 和 $v$ 都满足拉普拉斯方程 $\nabla^2 u = 0$,$\nabla^2 v = 0$。证明:
$$\iiint_V |\nabla u|^2 , dV = \oiint_S u \frac{\partial u}{\partial n} , dS$$
并说明物理意义。
证明:
利用第一 Green 恒等式,取 $\phi = u$,$\psi = u$:
$$\iiint_V (u \nabla^2 u + \nabla u \cdot \nabla u) , dV = \oiint_S u \frac{\partial u}{\partial n} , dS$$
由于 $\nabla^2 u = 0$:
$$\iiint_V |\nabla u|^2 , dV = \oiint_S u \frac{\partial u}{\partial n} , dS$$
物理意义:
- 左边:$\iiint_V |\nabla u|^2 , dV$ 代表场的"能量"(如静电能量)
- 右边:边界上的能量流入/流出
- 恒等式表明:调和函数的内部能量完全由边界条件决定
$\square$
总结与展望
核心知识框架回顾
场论体系
│
├─ 标量场
│ ├─ 梯度 ∇φ:方向导数、等值面
│ └─ 拉普拉斯算子 ∇²φ:调和函数
│
├─ 向量场
│ ├─ 散度 ∇·F:通量、源
│ ├─ 旋度 ∇×F:环流、涡
│ └─ 分类
│ ├─ 保守场:∇×F = 0,F = ∇φ
│ ├─ 无散场:∇·F = 0,F = ∇×A
│ └─ Helmholtz 分解:F = ∇φ + ∇×A
│
├─ Hamilton 算子 ∇
│ ├─ 作用于标量:∇φ(梯度)
│ ├─ 点积向量:∇·F(散度)
│ ├─ 叉积向量:∇×F(旋度)
│ └─ 重要恒等式
│
├─ 积分定理
│ ├─ 高斯公式:∭(∇·F)dV = ∬F·ndS
│ ├─ Stokes 公式:∬(∇×F)·ndS = ∮F·dr
│ └─ Green 恒等式
│
└─ 物理应用
├─ 电磁学:Maxwell 方程组
├─ 流体力学:Navier-Stokes 方程
├─ 热传导:热方程
└─ 量子力学:薛定谔方程
学习建议
- 理解物理意义:每个数学概念都有直观的物理图像
- 掌握计算技巧:熟练运用恒等式和积分定理
- 重视几何直觉:场线、等值面、通量、环流的几何意义
- 练习不同坐标系:直角、柱、球坐标的灵活运用
- 联系实际问题:通过物理应用加深理解
进一步学习方向
- 微分几何:流形上的微分算子、联络、曲率
- 微分拓扑:德拉姆上同调、Morse 理论
- 数学物理方程:偏微分方程理论与方法
- 张量分析:Einstein 求和约定、协变微分
- 广义相对论:Riemann 几何、Einstein 场方程
- 规范场论:Yang-Mills 理论、纤维丛
本讲义完
希望这份完整的场论知识体系能帮助您系统掌握梯度、散度、旋度与 Hamilton 算子的理论与应用!