经典物理偏微分方程

物理问题引入：当物理量在三维空间中弥散时

常微分方程研究的是孤立质点随时间演化的轨迹 $x (t)$ x(t)．然而，真实物理世界充满了 连续场（Continuous Fields）：

传播在空间中的光与声波 $u (x, y, z, t)$ u(x, y, z, t)；
传导在金属块各处的温度场 $T (x, y, z, t)$ T(x, y, z, t)；
充满全宇宙的静电势与引力势 $ϕ (x, y, z)$ phi(x, y, z)．

这些物理量同时依赖于时间 $t$ t 与空间坐标 $(x, y, z)$ (x, y, z)．描述连续介质动力学与场论的数学工具，就是包含了多元偏导数的 偏微分方程（Partial Differential Equations, PDE）．

1. 三类经典物理偏微分方程的物理推导与分类

大学物理主干理论几乎全部围绕以下三类二阶线性偏微分方程展开：

flowchart LR
    A["二阶线性偏微分方程"] --> B["双曲型 (Hyperbolic)<br/>波动方程"]
    A --> C["抛物型 (Parabolic)<br/>热传导 / 扩散方程"]
    A --> D["椭圆型 (Elliptic)<br/>拉普拉斯 / 泊松方程"]

    B --- B1["时间二阶导数 ∂²/∂t²<br/>有限波速传播、因果延迟、振荡"]
    C --- C1["时间一阶导数 ∂/∂t<br/>时间不可逆、能量耗散、光滑平滑"]
    D --- D1["无时间导数 (稳恒平衡态)<br/>全空间相互约束、调和势场"]

1.1 波动方程 (Wave Equation，双曲型）

考虑线密度为 $ρ$ rho、张力为 $T$ T 的一维绷紧弹性弦，取一小微元 $Δ x$ Delta x．微元两侧的垂直合力为 $T \sin θ_{2} - T \sin θ_{1} \approx T ({\frac{\partial u}{\partial x} |}_{x + Δ x} - {\frac{\partial u}{\partial x} |}_{x}) \approx T \frac{\partial^{2} u}{\partial x^{2}} Δ x$ T sintheta_2 - T sintheta_1 ~ Tleft(left.d u/d xright|_x+Delta x - left.d u/d xright|_xright) ~ T d^2 ud x^2Delta x．根据牛顿第二定律 $(ρ Δ x) \frac{\partial^{2} u}{\partial t^{2}} = T \frac{\partial^{2} u}{\partial x^{2}} Δ x$ (rhoDelta x)d^2 ud t^2 = T d^2 ud x^2Delta x，令波速 $v = \sqrt{T / ρ}$ v = T/rho：

\nabla^{2} u = \frac{1}{v^{2}} \frac{\partial^{2} u}{\partial t^{2}}

nabla^2 u = 1v^2d^2 ud t^2

物理特性：以确定速度 $v$ v 向外传播，波形在无耗散介质中保持不变，服从有限波速的因果约束．

1.2 热传导与扩散方程 (Heat Equation，抛物型）

考虑各向同性介质中的热流动：

傅里叶定律：热流密度与温度负梯度成正比： $q = - κ \nabla T$ q = -kappa nabla T；
能量守恒连续性：微元净吸收热量等于内能增加率： $- \nabla \cdot q = c ρ \frac{\partial T}{\partial t}$ -nabla*q = crho d Td t．

将两式联立，令热扩散率 $a^{2} = \frac{κ}{c ρ}$ a^2 = crho：

\nabla^{2} T = \frac{1}{a^{2}} \frac{\partial T}{\partial t}

nabla^2 T = 1a^2d Td t

物理特性：时间具有单向不可逆性（热力学第二定律熵增方向）；初始尖锐分布会被迅速平滑耗散．

1.3 泊松方程与拉普拉斯方程 (Poisson & Laplace Equation，椭圆型）

在静电场中，根据麦克斯韦高斯通量定理 $\nabla \cdot E = \frac{ρ}{ε_{0}}$ nabla*E = varepsilon_0，电场由标量电势负梯度给出 $E = - \nabla ϕ$ E = -nablaphi：

\nabla^{2} ϕ = - \frac{ρ (r)}{ε_{0}} (泊松方程)

nabla^2 phi = -rho(r)varepsilon_0 quad (泊松方程)

在无电荷的真空区域（ $ρ = 0$ rho = 0），退化为著名的 拉普拉斯方程：

\nabla^{2} ϕ = 0

nabla^2 phi = 0

物理特性：描述稳恒终态与平衡场，区域内任意一点的势等于其周围球面上势的平均值（平均值定理），不存在孤立极大/极小值．

2. 定解问题的完整构成

偏微分方程的通解通常包含任意未知函数，无法单独确定物理系统的具体行为．一个良定的 定解问题（Well-posed Problem） 必须包含三要素：

定解问题 = 偏微分方程 + 空间边界条件 + 初始条件

定解问题 = 偏微分方程 + 空间边界条件 + 初始条件

2.1 初始条件（Initial Conditions, IC）

给出时间起算时刻（ $t = 0$ t=0）全空间的物理状态：

波动方程（时间二阶）：需要初始位移 $u (x, 0) = ϕ (x)$ u(x,0) = phi(x) 和初始速度 ${\frac{\partial u}{\partial t} |}_{t = 0} = ψ (x)$ left.d u/d tright|_t=0 = psi(x)；
热传导方程（时间一阶）：只需初始温度分布 $T (x, 0) = f (x)$ T(x,0) = f(x)；
泊松方程（静态平衡）：无需初始条件．

2.2 三类典型边界条件（Boundary Conditions, BC）

在空间区域边界 $\partial Ω$ dOmega 上施加的物理约束：

边界类型	数学表述	典型物理场景
第一类边界（狄利克雷 Dirichlet）	${u \|}_{\partial Ω} = f (r)$ left.urightrvert_dOmega=f(boldsymbol r)	固定弦端位移为零；接地导体电势为零
第二类边界（诺伊曼 Neumann）	${\frac{\partial u}{\partial n} \|}_{\partial Ω} = g (r)$ left.d ud nrightrvert_dOmega=g(boldsymbol r)	绝热边界的法向温度梯度为零
第三类边界（罗宾 Robin）	${(α u + β \frac{\partial u}{\partial n}) \|}_{\partial Ω} = h (r)$ left.(alpha u+betad ud n)rightrvert_dOmega=h(boldsymbol r)	牛顿冷却对流换热边界

3. 线性偏微分方程的叠加原理

若 $u_{1} (r, t)$ u_1(r, t) 与 $u_{2} (r, t)$ u_2(r, t) 是齐次偏微分方程及齐次边界条件的两个特解，则它们的线性组合：

u (r, t) = C_{1} u_{1} (r, t) + C_{2} u_{2} (r, t)

u(r, t) = C_1 u_1(r, t) + C_2 u_2(r, t)

依然满足方程与齐次边界条件．求解偏微分方程的通用纲领：通过某种结构化方法寻找一组正交完备的「简单特解基底」，再通过积分或级数无限叠加，去逐一满足给定的复杂初始条件或非齐次边界！

4. 学习衔接

下一节：学习将偏微分降维求解的最强方法，进入偏微分方程分离变量法；
物理应用：在电磁学：静电场中掌握拉普拉斯方程边值求解．

本页面最近更新：2026/10/1 10:49:01，更新历史
发现错误？想一起完善？在 GitHub 上编辑此页！
本页面贡献者：Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用