1. 课程导入与守恒定律数值校验总览

在计算流体力学、结构动力学、电磁场仿真以及多物理场耦合分析中，数值模拟已经成为与理论分析、实验观测并列的第三种研究手段。然而，任何离散化方法都会引入截断误差、舍入误差和时间推进误差，使得数值解偏离真实解。守恒定律数值校验的核心任务，就是利用质量、动量、能量等物理量在离散系统中的收支关系，判断计算结果是否"物理上可信"。本课程将围绕这一目标，系统讲解从概念到实现的完整方法体系。

为什么守恒性是最可靠的校验标尺

守恒定律是自然界最普遍的约束。无论采用有限体积、有限元还是有限差分，只要离散格式对守恒量满足"局部通量平衡、全局收支相消"，数值解就应当保持相应的守恒性质。与逐点误差比较不同，守恒校验不依赖精确解，也不需要实验数据，仅凭计算域边界上的通量记录和内部源项即可完成判断，因此特别适合工程中缺乏参考解的复杂问题。

以一个二维不可压缩流动的有限体积求解为例：若每个控制体的质量通量满足 Σ(ρu·A)_out − Σ(ρu·A)_in = 0，则全域质量守恒残差应仅由边界条件与舍入误差贡献。实际计算中，若发现某时刻总质量相对漂移超过 1e-6，往往意味着通量重构或边界处理存在缺陷。这类判断正是本课程要训练的核心能力。

校验方法的分类框架

守恒校验方法可按"校验对象"和"校验尺度"两个维度划分。按对象分，包括质量守恒、动量守恒、能量守恒、组分守恒以及电荷守恒等；按尺度分，包括局部（单元级）校验、区域（子域级）校验和全局（全域级）校验。局部校验灵敏度高，能定位问题单元；全局校验计算量小，适合作为每次运行后的快速健康检查。实际工程中通常采用"全局粗筛 + 局部精查"的两级策略。

此外，还可按实现方式分为解析校验（构造满足守恒的解析解）、收支审计（记录边界通量与源项）、残差监控（跟踪离散方程残差范数）和交叉验证（不同格式或网格对比）。本课程将逐一展开这些方法的实现细节与适用边界。

守恒校验在数值模拟工作流中的位置

一个完整的仿真工作流通常包括几何建模、网格划分、物理设置、求解、后处理五个阶段。守恒校验不应只在后处理阶段"看一眼总量"，而应贯穿求解全过程：在求解器内部按时间步记录边界通量累计值，形成"守恒审计日志"；在自适应网格加密时，检查粗细网格交界面的通量是否满足守恒约束；在并行计算中，核对分区边界上的通量交换是否对称。把校验嵌入工作流，才能及早发现格式缺陷，而不是等到结果异常再回头排查。

需要强调的是，守恒性成立并不等于数值解精确。一个一阶迎风格式可能严格守恒，但数值耗散很大；反之，某些高阶格式局部精度高，却可能因边界处理不当破坏全局守恒。因此守恒校验是"必要条件"而非"充分条件"，应与网格收敛性分析、解析解对比配合使用。

实战项目预告与学习建议

本课程的实战项目将围绕二维不可压缩方腔流动展开。学员需要自行实现或调用一个有限体积求解器，在顶盖驱动条件下运行若干时间步，并按以下思路完成校验：首先在求解器中插入通量累加器，记录每个时间步四个边界上的质量通量与动量通量；其次计算总质量相对漂移 |M(t)−M(0)|/M(0) 与总动量变化率，与理论收支对比；然后设定阈值（例如质量相对漂移 1e-6、动量收支相对误差 1e-4），标记超阈值时刻；最后回溯该时刻的边界单元，检查是否存在通量符号错误或边界条件不一致。

建议学习时准备一份"校验日志"模板，字段包括：时间步、边界通量累计、源项累计、理论收支、相对误差、判定结论。每完成一个知识点，就把对应方法填入日志并实际运行验证。守恒校验是一项动手性很强的工作，只有亲手记录数据、观察漂移、定位异常，才能真正建立对数值结果可信度的判断力。