拉普拉斯算子在不同坐标系中的表达式推导

理解核心概念：拉普拉斯算子（记作 $\nabla^2$ 或 $\Delta$）是向量分析中的一个核心算子，可以简单理解为描述某个物理量（如温度、电势）在空间中某一点变化“激烈程度”的二阶微分工具。在直角坐标系中，它的形式我们很熟悉，但在其他坐标系（如柱坐标、球坐标）下，形式会变得复杂。掌握其推导方法是解决许多物理和工程问题的基础。

1. 从熟悉的起点开始：直角坐标系下的定义

在三维直角坐标系 $(x, y, z)$ 中，拉普拉斯算子作用于一个标量函数 $f$ 的定义最为直观。

1. 写出 一阶偏导数：首先计算函数 $f$ 对三个坐标变量 $x, y, z$ 的一阶偏导数，分别记为 $\frac{\partial f}{\partial x}$, $\frac{\partial f}{\partial y}$, $\frac{\partial f}{\partial z}$。
2. 计算 二阶偏导数：接着对每个一阶偏导数再次求导，得到三个二阶偏导数：$\frac{\partial^2 f}{\partial x^2}$, $\frac{\partial^2 f}{\partial y^2}$, $\frac{\partial^2 f}{\partial z^2}$。
3. 求和：将这三个二阶偏导数相加，就得到了直角坐标系下的拉普拉斯算子表达式。
   $$ \nabla^2 f = \Delta f = \frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2} + \frac{\partial^2 f}{\partial z^2} $$
   这就是我们进行坐标变换的基准。下一步，我们将看到如何在不同的“描述位置的方式”（坐标系）下，写出这个“求和”的等效形式。

2. 推导的核心工具：坐标变换与度规张量

推导的关键在于，我们不能简单地“套用”直角坐标下的求和公式。必须考虑新坐标系中坐标轴可能弯曲、不正交的特性。

1. 引入 新坐标系：设新坐标系为 $(u^1, u^2, u^3)$。例如，柱坐标系中 $(u^1, u^2, u^3) = (r, \theta, z)$，球坐标系中 $(u^1, u^2, u^3) = (r, \theta, \phi)$。
2. 定义 位置矢量：空间任一点的位置可以用新坐标表示为 $\mathbf{r}(u^1, u^2, u^3)$。
3. 计算 基向量：对每个坐标 $u^i$ 求偏导，得到一组基向量：
   $$ \mathbf{e}_i = \frac{\partial \mathbf{r}}{\partial u^i}, \quad i=1,2,3 $$
   这组基向量通常不是单位向量，也不一定相互垂直。例如在柱坐标中，$\mathbf{e}_r$ 和 $\mathbf{e}_\theta$ 的方向随点变化。
4. 计算 度规张量：为了量化新坐标系下的“长度”和“角度”，我们需要计算基向量的内积，构成一个对称矩阵，称为度规张量 $g_{ij}$：
   $$ g_{ij} = \mathbf{e}_i \cdot \mathbf{e}_j $$
   度规张量包含了坐标系的所有几何信息。例如，在直角坐标系中，$g_{ij}$ 是一个单位矩阵。

3. 推导拉普拉斯算子的通用公式

利用度规张量，我们可以建立一个适用于任意正交曲线坐标系（即基向量相互垂直的坐标系，柱坐标和球坐标都属于此类）的拉普拉斯算子通用公式。

1. 写出 度规张量的行列式：计算 $g = \det(g_{ij})$。对于正交坐标系，由于非对角元素为零，这个行列式就是对角线元素的乘积，即 $g = g_{11}g_{22}g_{33}$。
2. 定义 尺度因子：令 $h_i = \sqrt{g_{ii}} = |\mathbf{e}_i|$。$h_i$ 被称为尺度因子，它表示当坐标 $u^i$ 发生单位变化时，空间中移动的实际弧长。例如，柱坐标中 $h_r=1$, $h_\theta=r$, $h_z=1$。
3. 套用 通用公式：对于正交曲线坐标系，拉普拉斯算子的通用表达式为：
   $$ \nabla^2 f = \frac{1}{h_1 h_2 h_3} \left[ \frac{\partial}{\partial u^1} \left( \frac{h_2 h_3}{h_1} \frac{\partial f}{\partial u^1} \right) + \frac{\partial}{\partial u^2} \left( \frac{h_3 h_1}{h_2} \frac{\partial f}{\partial u^2} \right) + \frac{\partial}{\partial u^3} \left( \frac{h_1 h_2}{h_3} \frac{\partial f}{\partial u^3} \right) \right] $$
   这个公式是后续所有推导的“万能钥匙”。

4. 实战一：推导柱坐标下的表达式

柱坐标 $(r, \theta, z)$ 与直角坐标 $(x, y, z)$ 的关系为：$x = r\cos\theta$, $y = r\sin\theta$, $z = z$。

1. 计算 基向量和尺度因子：
   • $\mathbf{e}_r = (\cos\theta, \sin\theta, 0)$，故 $h_r = |\mathbf{e}_r| = 1$。
   • $\mathbf{e}_\theta = (-r\sin\theta, r\cos\theta, 0)$，故 $h_\theta = |\mathbf{e}_\theta| = r$。
   • $\mathbf{e}_z = (0, 0, 1)$，故 $h_z = 1$。
2. 代入 通用公式：将 $h_r=1$, $h_\theta=r$, $h_z=1$ 和坐标 $(u^1, u^2, u^3) = (r, \theta, z)$ 代入步骤 3 的公式。
   $$ \nabla^2 f = \frac{1}{1 \cdot r \cdot 1} \left[ \frac{\partial}{\partial r} \left( \frac{r \cdot 1}{1} \frac{\partial f}{\partial r} \right) + \frac{\partial}{\partial \theta} \left( \frac{1 \cdot 1}{r} \frac{\partial f}{\partial \theta} \right) + \frac{\partial}{\partial z} \left( \frac{1 \cdot r}{1} \frac{\partial f}{\partial z} \right) \right] $$
3. 化简 得到最终形式：
   $$ \nabla^2 f = \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial f}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 f}{\partial \theta^2} + \frac{\partial^2 f}{\partial z^2} $$
   请注意，与直角坐标系相比，第一项多了一个 $1/r$ 的系数，第二项的分母变成了 $r^2$。这反映了在极坐标系下，径向和角向的二阶导数需要特殊的修正。

5. 实战二：推导球坐标下的表达式

球坐标 $(r, \theta, \phi)$ 与直角坐标的关系为：$x = r\sin\theta\cos\phi$, $y = r\sin\theta\sin\phi$, $z = r\cos\theta$。这里 $\theta$ 是极角，$\phi$ 是方位角。

1. 计算 基向量和尺度因子：
   • $\mathbf{e}_r = (\sin\theta\cos\phi, \sin\theta\sin\phi, \cos\theta)$，故 $h_r = 1$。
   • $\mathbf{e}_\theta = (r\cos\theta\cos\phi, r\cos\theta\sin\phi, -r\sin\theta)$，故 $h_\theta = r$。
   • $\mathbf{e}_\phi = (-r\sin\theta\sin\phi, r\sin\theta\cos\phi, 0)$，故 $h_\phi = r\sin\theta$。
2. 代入 通用公式：将 $h_r=1$, $h_\theta=r$, $h_\phi=r\sin\theta$ 和坐标 $(u^1, u^2, u^3) = (r, \theta, \phi)$ 代入步骤 3 的公式。
   $$ \nabla^2 f = \frac{1}{1 \cdot r \cdot r\sin\theta} \left[ \frac{\partial}{\partial r} \left( \frac{r \cdot r\sin\theta}{1} \frac{\partial f}{\partial r} \right) + \frac{\partial}{\partial \theta} \left( \frac{r\sin\theta \cdot 1}{r} \frac{\partial f}{\partial \theta} \right) + \frac{\partial}{\partial \phi} \left( \frac{1 \cdot r}{r\sin\theta} \frac{\partial f}{\partial \phi} \right) \right] $$
3. 化简 得到最终形式：
   $$ \nabla^2 f = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial f}{\partial r} \right) + \frac{1}{r^2 \sin\theta} \frac{\partial}{\partial \theta} \left( \sin\theta \frac{\partial f}{\partial \theta} \right) + \frac{1}{r^2 \sin^2\theta} \frac{\partial^2 f}{\partial \phi^2} $$
   这个表达式比柱坐标更复杂，包含了 $\sin\theta$ 项，这源于球坐标中经线圈和纬线圈的弯曲特性。