1. 非线性PDE问题识别与建模
在进入离散格式与求解器设计之前，最容易被低估却决定后续全部工作成败的一步，是把物理现象翻译成正确的数学方程。线性PDE的建模几乎可以靠"套模板"完成，而非线性PDE的建模必须回答三个问题：非线性从哪里来？它属于哪一类？这个方程在数学上是否"讲得通"（适定性）？本章围绕这三个问题展开，并以Burgers方程和反应扩散系统作为贯穿案例。
一、非线性项的来源与分类
非线性项并非凭空出现，它总是对应某种具体的物理机制。工程与物理中最高频的三类非线性项可以这样识别：
- 对流非线性（输运型）：当输运速度本身依赖于被输运的量时出现，典型形式为 \(u\,u_x\) 或 \(\nabla\cdot(u\,\mathbf{u})\)。它来自流体惯性、激波形成、交通流拥堵波等。物理直觉是"快的地方推得更快"，因此容易自发陡化。
- 反应非线性（源汇型）：当局部生成/消耗速率是状态量的非线性函数时出现，如 \(u^2\)、\(u(1-u)\)、\(\mathrm{e}^{u}\)。它来自化学反应动力学、种群增长、燃烧放热。其危险在于可能有限时间爆破。
- 扩散耦合非线性：扩散系数依赖状态量，如 \(\nabla\cdot\big(D(u)\nabla u\big)\)，见于多孔介质渗流（\(D\propto u^{m}\)）、非线性热传导。它改变的是"传播速度随幅度变化"这一特性。
识别时的一个实用判据是：把方程中每一项单独拿出来，问它是否对 \(u\) 齐次线性。若某项含 \(u\) 的乘积、幂次或非线性函数，即为非线性项。进一步要判断它是"守恒型"（可写成散度形式）还是"非守恒型"，这直接决定有限体积离散能否直接使用。
二、从现象到方程的建模流程
建模不是写下一个方程就结束，而是一个可复现的推导链条。以Burgers方程为例：考虑一维无粘流体，取控制体，写出质量守恒的积分形式，令控制体收缩并利用散度定理，得到守恒律 \(u_t + \big(\tfrac{1}{2}u^2\big)_x = 0\)；若再叠加粘性扩散通量 \(-\nu u_x\)，则得到粘性Burgers方程
u_t + u u_x = ν u_xx , x ∈ (0,L), t > 0
u(x,0) = u0(x) （初值）
u(0,t) = u(L,t) = 0 （Dirichlet边界）注意这里 \(u u_x\) 与 \(\big(\tfrac12 u^2\big)_x\) 在光滑解下等价，但在激波处不等价——这正是后续需要熵条件与守恒型格式的根源。建模阶段就应明确采用哪一种写法，并记录在文档中。
反应扩散系统则来自"局部反应 + 空间扩散"的叠加。以经典的FitzHugh–Nagumo型双组分系统为例：
u_t = D_u u_xx + u(1-u)(u-a) - v
v_t = D_v v_xx + ε(b u - c v)
参数示例：D_u = 0.01, D_v = 0.05, a = 0.1, ε = 0.05, b = 0.5, c = 0.2其中 \(u(1-u)(u-a)\) 是三次反应非线性，\(-v\) 是耦合项。参数取值仅为教学示例，用于说明不同参数区间会带来单稳、双稳或振荡行为，并非实测数据。
三、适定性初步分析：建模阶段就要做的判断
适定性包含存在性、唯一性与对初值/参数的连续依赖。建模阶段不必做严格证明，但必须做定性判断，因为它决定离散策略。三个高频风险点值得逐一排查：
- 有限时间爆破：若反应项增长过快（如 \(u^2\) 且无耗散），解可能在有限时刻发散。判断方法是看源项与扩散项的量级竞争，或对常微分比较方程做估计。若存在爆破风险，时间步长必须自适应，且需监控解的范数。
- 激波与不唯一性：对流非线性使光滑解在有限时间陡化，弱解不唯一，必须附加熵条件。这意味着离散必须采用守恒型格式（如Godunov型、迎风通量），否则会收敛到非物理解。
- 退化扩散：当 \(D(u)\to 0\) 时方程类型退化，解的正则性下降，标准中心差分可能产生负值。此时需保证离散的保正性。
四、案例对照：Burgers方程与反应扩散系统
两个案例恰好代表两类不同的非线性主导机制，对照理解有助于建立"看到方程就知道该防什么"的直觉。
| 对比维度 | 粘性Burgers方程 | 反应扩散系统 |
|---|---|---|
| 主导非线性 | 对流型 \(u u_x\) | 反应型 \(f(u,v)\) + 耦合 |
| 典型现象 | 激波陡化、粘性层 | 斑图、行波、振荡 |
| 主要数学风险 | 弱解不唯一，需熵条件 | 刚性源项、可能爆破 |
| 离散关键 | 守恒型迎风通量 | 源项隐式/算子分裂 |
| 时间尺度 | 由 \(\nu\) 与网格决定 | 反应与扩散尺度差异大 |
以Burgers方程为例，若取 \(\nu=0.01\)、\(L=1\)、初值 \(u_0(x)=\sin(2\pi x)\)，在 \(t\approx 1/(2\pi)\) 附近解开始陡化，粘性会阻止真正间断的形成，但梯度极大。建模时就应预判：网格需在陡化区加密，时间推进宜用隐式或半隐式以放宽 \(\Delta t\) 限制。反应扩散系统中，若 \(\varepsilon\) 很小，反应项的特征时间远小于扩散时间，属于典型刚性问题，显式格式会被稳定性死死限制，必须对源项做隐式处理或采用IMEX分裂。
完成本章的识别与建模后，我们得到的不仅是一组方程，还包括一份关于非线性类型、守恒性质、适定性风险与参数来源的结构化记录。这份记录将直接指导后续的网格设计、格式选择与时间推进策略——建模的质量，决定了离散与求解的上限。