1. 工程问题中的偏微分方程建模

偏微分方程（PDE）数值求解的第一步，不是写代码，也不是选格式，而是把工程问题翻译成一个可计算的数学模型。这个模型由三部分组成：控制方程描述物理守恒律，边界条件描述求解域边界上的约束，初始条件描述时间演化问题的起点。三者缺一不可，且必须与所选的坐标系、符号约定保持一致。本节从热传导、结构振动、流体流动三个典型场景出发，说明如何完成这一翻译过程。

一、从物理守恒到控制方程

工程中的PDE几乎都来自守恒律。以一块厚度为 L 的均质平板一维导热为例，取能量守恒的微元分析：流入与流出微元的热流之差，等于微元内能的增量。设温度 T(x,t)，导热系数 k，密度 ρ，比热容 c，无内热源时得到

ρc ∂T/∂t = ∂/∂x ( k ∂T/∂x )

若 k 为常数，即化为经典的一维热传导方程 ∂T/∂t = α ∂²T/∂x²，其中热扩散率 α = k/(ρc)。这一步的关键是明确未知量是什么、自变量有哪些、系数是否随空间或温度变化。工程上很多错误并非来自数值格式，而是建模阶段把变系数问题当成了常系数问题。

结构振动则来自牛顿第二定律与材料本构。一根均匀细长杆的纵向振动，设轴向位移 u(x,t)，弹性模量 E，密度 ρ，截面积 A，可得

ρA ∂²u/∂t² = ∂/∂x ( EA ∂u/∂x )

与热传导方程相比，它关于时间是二阶的，因此需要两个初始条件（初始位移与初始速度），这是建模时必须提前识别的差异。流体流动的控制方程来自质量与动量守恒，一维不可压缩流动可写成连续性方程 ∂u/∂x = 0 与动量方程 ρ(∂u/∂t + u ∂u/∂x) = −∂p/∂x + μ ∂²u/∂x²。这里的非线性对流项 u ∂u/∂x 是后续数值求解困难的重要根源，建模阶段就应记录其存在。

二、边界条件与初始条件的工程识别

控制方程只描述域内规律，边界条件才把具体工程对象“钉”成唯一解。以平板导热为例，若左表面 x=0 与高温流体接触，右表面 x=L 保持恒定温度，则左端为第三类边界 −k ∂T/∂x|₀ = h(T_∞ − T(0,t))，右端为第一类边界 T(L,t) = T_L。注意法向方向的定义：若取外法向为正，左端外法向指向 −x，符号必须相应调整。工程建模中大量符号错误源于没有先声明法向约定。

结构振动的边界条件同样需要从实际约束翻译：固定端为 u=0（位移边界），自由端为 EA ∂u/∂x=0（力边界），若端部连接弹簧则写成 EA ∂u/∂x + k_su = 0。初始条件则来自静止释放或给定初速度，例如 u(x,0)=0、∂u/∂t(x,0)=v₀。流体问题中，壁面通常取无滑移 u=0，入口给定速度或流量，出口常取压力参考值或充分发展条件。这些选择直接决定数值求解的稳定性与物理合理性，必须在建模阶段写清楚，而不是留给求解器“默认处理”。

三、建模交付物：一页式模型说明

为了让建模结果可复核、可交接，建议每个问题形成一份简短的模型说明。以下是一个平板非稳态导热的示例（参数为教学假设值，非实测）：

有了这张表，后续选择空间离散格式、时间推进方法和验证方案才有依据。例如由 α 与 L 可估算特征时间 L²/α ≈ 2.9×10³ s，从而判断需要模拟多长时间、时间步长应取多大。建模阶段的这些量纲分析与特征尺度估计，往往比盲目试算更能提前暴露问题。

最后要强调，建模不是一次完成的。数值结果若出现非物理振荡、能量不守恒或对网格异常敏感，应回到模型层检查边界条件是否过约束、系数是否突变、初始条件是否与边界条件相容（例如初始温度与边界温度在角点处矛盾会引入人为奇性）。把建模说明与数值设置分开记录，是排查这类问题的有效习惯。