1.4 从物理到方程的建模之路

1.4 从物理到方程的建模之路

本节摘要：方程不是查表得来的。本节完整走一遍建模路径：从微元守恒律出发推出热传导方程，指出每一步近似藏在哪里；再讲清初始条件与边界条件为何是方程的"另一半"，以及三类边界条件各自的物理含义。学完本节，面对新场景应能自己写出方程与定解条件。

建模的三步会谈

把物理场景翻译成偏微分方程，通常经过三步。第一步找守恒量：质量、动量、能量，总有一个在所研究的过程中守恒或按已知速率变化。第二步选描述量：决定用哪个函数刻画系统状态，温度、浓度、位移或速度，这一步决定了未知函数的维数与阶数。第三步补本构关系：守恒律只给出通量的散度等于源，通量本身如何响应状态梯度（傅里叶定律、菲克定律、胡克定律），需要额外的经验定律封口。

三步会谈的产物总是同一个结构：一个含时间与空间导数的方程，加上一套初始条件与边界条件。缺少任何一件，问题都算不出来。本节以热传导为例把三步走全，因为它是推导最短、近似最透明的一条。

手推热传导方程：每一步都看得见

考虑一根均匀细杆，侧壁绝热，温度只沿杆轴变化：温度场是 u(x,t)。取杆上一小段，位于 x 到 x 加 Δx 之间。能量守恒说：这一小段内能的增加率，等于从左端流入的热功率减去从右端流出的热功率，加上段内热源的产生率。

左端流入功率记 q(x)，右端流出功率记 q(x + Δx)。段内内能正比于温度对位置的积分，增加率是密度与比热的乘积乘以温度时间导数的区间积分。写成等式：

∫[x, x+Δx] ρc · ∂u/∂t ds = q(x) − q(x + Δx) + ∫[x, x+Δx] Q ds，

其中 Q 是单位长度热源强度。两边除以 Δx 并令 Δx 趋于零，左边趋于 ρc·∂u/∂t，右边的差商趋于 −∂q/∂x，得

ρc·∂u/∂t = −∂q/∂x + Q。

到目前为止只用了一般的守恒律，对任何输运过程都成立。现在需要第三步：热流 q 与温度梯度怎么挂钩？傅里叶定律给出 q = −k·∂u/∂x——热往冷处流，流量正比于温度梯度。代入并整理，令 a² = k/(ρc)：

∂u/∂t = a²·∂²u/∂x² + Q/(ρc)。

这就是带源项的热传导方程。回头看，推导里藏了三个近似：杆无限细（截面温度均匀）、材料均匀（ρ、c、k 为常数）、傅里叶定律在极端梯度下失效（极低温或极高频时需要非局域修正）。第三章将看到非局域修正的一种形态——分数阶扩散在图像处理中的应用，那里梯度假设同样需要放松。

这条推导的价值不在结论（教科书都有），而在可迁移性：把"温度"换成"浓度"，傅里叶定律换成菲克定律，得到扩散方程；换成"电荷密度"，得到电荷弛豫方程。守恒律是骨架，本构是血肉——血肉一换，新的方程就有了。

定解条件：方程的另一半

方程只规定"相邻时刻、相邻位置之间如何协调"，要确定唯一解，还须告诉它"起点在哪、边上有谁管着"。这就是定解条件，分两类。

初始条件：给整个空间区域在初始时刻的完整状态。热传导方程是一阶时间方程，给一个初始温度分布 u(x,0) = φ(x) 即可；波动方程是二阶时间方程，必须给两个：初始位移 u(x,0) 与初始速度 u_t(x,0)。少给一个，解不唯一；这是"时间阶数决定初始条件个数"的一般规律。

边界条件：区域边界上的约束，物理上最常见的有三种，以传热为语境列表如下：

三类边界的选择不是数学趣味，而是工程描述能力：同一个杆，三种边界给出三种完全不同的温度场。尤其是绝热端对应诺伊曼零通量，这个条件在第四章的差分实验里会实际用到——处理不当，数值解会凭空丢失或注入热量。

图解：一根杆上的三类边界

三类边界条件物理场景示意

建模中的高频错误

错误一：条件给多。 有人在杆的两端都给温度、又给初始分布，还抱怨解不出来。超定的定解问题一般无解；条件是"钉死"自由度用的，给过头了反而互相矛盾。判断标准是数自由度：时间每高一阶多给一层初始条件，空间每延展一维多给一圈边界条件。

错误二：边界条件与物理不符。 最常见的把"离热源很远"设成固定温度（狄利克雷），实际那里更接近"不受扰动"（诺伊曼零梯度）。远场条件的选型直接影响解在全域的形态，第六章应用案例里会看到它对散热器设计结论的实质影响。

错误三：把源项忘了量纲。 Q 的单位是单位体积单位时间的能量，很容易跟总功率混淆。写完方程后做一次量纲检查——方程两边量纲不一致，后面的推导再漂亮也是空中楼阁。

⚠️ 常见坑：写程序前先把定解问题用纸写全——方程、区域、每一类边界、初始分布。数值实验里大量"解看起来不对"的案例，根源是边界条件写串了类型或下标错位，而不是求解器本身的问题。

迁移演练：同一条推导换两个肉身

热传导推导的价值在可迁移性，光说不练不算数，这里把"守恒律加本构"的模板再执行两遍。肉身一：污染物扩散。守恒量换成质量，描述量换成浓度 c(x,t)，微元质量平衡的推演与温度逐字平行：流入减流出加源等于存量变化率；本构换成菲克定律——通量 q = −D·∂c/∂x（物质从浓处向稀处扩散）。得到的 c_t = D·c_xx 与热方程只差记号名姓，连稳态、核函数、谱衰减全套家当都可整体搬用。地下水污染预测、药物缓释设计，用的就是这条方程。

肉身二：弦的横振动。守恒量换成动量，微元是长度 Δx 的一小段弦，牛顿第二定律的纵向形式换成横向：质量乘横向加速度等于两端张力的横向分量之差。小振动近似下斜角很小，横向合力约等于 T·(u_x 的差商) 乘 Δx，除以质量密度 ρ 后得到 u_tt = (T/ρ)·u_xx——波动方程带着波速 c = √(T/ρ) 落地。这个公式顺带回答了一个调音问题：弦的音高与张力平方根成正比、与密度平方根成反比——吉他手拧弦轴的行为，就是在调波速。

两条演练的共同复盘：守恒律的骨架一个字没改，改的只是本构与物理参数。新场景建模时的自查次序由此固定：先问什么量守恒，再问通量怎么响应梯度，最后问界面（边界）上发生什么——三问齐全，方程自然成形。

本节要点

• 建模三步：找守恒量、选描述量、补本构关系；守恒律给骨架，本构定血肉；
• 热传导推导：微元能量守恒加傅里叶定律，过程中藏有细杆、均匀材料、局域梯度三个近似；
• 定解条件：初始条件个数由时间阶数决定；边界条件三型的物理含义各有所指；
• 高频错误：超定、远场条件选型不当、源项量纲混乱，都属于"方程写错"而非"算错"。

方程与条件凑齐，一条定解问题就诞生了。但它一定有解吗？解唯一吗？对输入的小扰动敏感吗？下一章的分类与适定性理论，回答的正是这三个问题。