1. 课程导入与PDE定解问题概览

偏微分方程（PDE）本身描述的是某一类物理或工程现象在区域内部必须满足的普遍规律，例如热传导中的能量守恒、弦振动中的牛顿第二定律、静电场中的高斯定律。但仅有方程，解并不唯一——同一根杆，两端温度固定为0℃和100℃，与两端都绝热，温度分布完全不同。因此，一个可求解、有物理意义的定解问题，必须由三部分共同构成：偏微分方程（控制方程）、初始条件与边界条件。本课程的核心，就是系统掌握后两者的设计方法，并理解它们如何与方程一起决定唯一解。

一、定解问题的三要素与角色分工

可以把定解问题类比为“剧本+开场状态+舞台边界”：PDE是剧情发展的规则，初始条件规定t=0时刻整个系统的状态，边界条件规定区域边界上每个时刻发生什么。三者缺一不可，且必须相容——初始条件在边界上的取值要与边界条件在t=0时刻的取值一致，否则会出现物理上不合理的跳跃，导致数值求解时在角点附近产生剧烈振荡甚至发散。

从数学上看，初始条件的数量由方程对时间的阶数决定：热传导方程（一阶时间导数）只需一个初始条件 u(x,0)=f(x)；波动方程（二阶时间导数）需要两个，即初始位移 u(x,0)=f(x) 与初始速度 u_t(x,0)=g(x)。边界条件的数量与类型则由空间算子的阶数和边界上的物理约束决定，常见的有Dirichlet（给定函数值）、Neumann（给定法向导数，常对应通量或绝热）、Robin（给定值与导数的线性组合，对应对流换热）三类。

二、初始条件与边界条件的典型设计错误

初学者最常见的问题不是不会解方程，而是条件设计不当。以下三类错误在工程中反复出现：

• 条件数量不匹配：对二阶时间导数的波动方程只给初始位移、漏掉初始速度，导致解中含任意函数，数值上表现为结果依赖求解器内部默认值。
• 物理类型错配：把“绝热”边界误设为Dirichlet条件（固定温度），相当于在边界上强行加了一个恒温源，热量被凭空抽走或注入，能量不再守恒。
• 角点不相容：矩形区域左下角同时受 x 方向与 y 方向边界条件约束，若两者在角点给出不同值，离散后该点方程矛盾，迭代残差无法下降。

三、本课程的整体框架

课程按“概念—分类—设计—验证—实战”推进。前段建立定解问题的数学语言，明确三类边界条件的物理含义与数学表达；中段进入设计方法，讨论非齐次边界条件的齐次化、移动边界与自由边界、周期性边界与对称性简化；后段聚焦数值实现中的条件离散、稳定性与收敛性验证，最后以一个完整的实战项目收束。

贯穿全课的实战项目设定为：一维金属杆的瞬态热传导分析。杆长 L=1 m，材料导热系数 k=45 W/(m·K)，密度 ρ=7800 kg/m³，比热容 c=460 J/(kg·K)，热扩散率 α=k/(ρc)≈1.25×10⁻⁵ m²/s。初始温度场设为 u(x,0)=20+80·sin(πx/L) ℃，左端 x=0 保持 100℃（Dirichlet），右端 x=L 为对流换热边界（Robin），环境温度 20℃、对流系数 h=50 W/(m²·K)。学员将亲手写出该问题的完整定解描述，用有限差分离散，观察不同时刻温度曲线的演化，并通过能量平衡校验结果合理性。

四、学习目标与自检标准

完成本章及后续内容后，你应当能够：给定一个物理场景，独立写出包含方程、初始条件、边界条件的完整定解问题；识别并修正条件数量、类型、相容性方面的错误；对数值结果给出基于物理守恒的判断，而不是只看曲线是否“光滑”。

一个可操作的自检清单如下，建议在每次建模后逐项核对：