1 概述与基本思想
守恒律（conservation laws）用来描述物理量在时空演化过程中不会凭空产生或消失的规律。所谓“不会凭空”，意味着对任意空间区域而言，区域内物理量的增加只可能来自边界处的净流入；区域内的减少则对应边界处的净流出或内部向外输运。把这种“区域—边界”关系写成方程，就得到守恒律的结构。
守恒律并不总是意味着物理过程完全不发生变化。物理量可能在不同形式之间转换（例如能量在动能与势能间转化），也可能在局部发生再分布，但只要所选取的物理量在定义范围内没有源与汇，就应满足“流入与流出守恒”的原则。
1.1 守恒律要解决的问题：从“全局不变”到“局域不变”
传统直觉往往从全局观点出发：例如封闭系统中的某些量总和是否保持不变。守恒律将这种直觉提升到局域层面，即不仅关注某个整体量在整个系统上的变化，还要描述任意小控制体内发生了什么。关键转折在于：一旦局域守恒成立，就可以通过边界通量把整体行为追溯到空间分布的原因，从而为连续介质建模与方程推导提供统一框架。
在工程与物理建模中，局域守恒还带来可验证的结构约束：当某个物理量守恒时，方程形式必然呈现“时间变化 + 空间输运 = 源项”的特征，使得理论推导、边界条件设定与数值离散都更有针对性。
1.2 通量（flux）与净流入/净流出的物理含义
通量（flux）刻画物理量沿空间方向的“穿过某个面”的速率。若考虑某种标量物理量，其通量可以理解为单位面积上、单位时间内通过界面的量。对于控制体，边界由若干面组成，通量方向决定了“流入”或“流出”。因此，边界上各面通量的合计（更准确地说是其几何加权的净效果）等于该控制体内物理量变化的驱动项之一。
从形式上看，通量守恒常写作“通量散度”出现在线性或非线性偏微分方程中；当通量在空间上有不均匀变化时，散度反映了净流入/净流出在局部的差别。
1.3 源项（source/sink）与“允许生成或湮灭”的情形
源项（source/sink）表示物理量在控制体内部由于某些机制被生成或被消耗。若源项为零，方程退化为纯守恒形式；若源项非零，则守恒律扩展为“守恒 + 生成/湮灭”。例如质量守恒在许多经典场景中可视为无源，但在带有化学反应或相变的模型中，质量可能以不同形式或组分之间的方式转化，从而在“组分层面”出现有效源项。
源项的存在并不破坏框架本身，而是把“局域变化来源”分解为两类：一类来自边界净输运，另一类来自内部源汇过程。这样既能描述生成/消耗，也保留通过通量理解输运机制的优势。
2 通量守恒的数学框架
守恒律最常见的表达方式是将“守恒量的变化”与“通量通过空间的变化”联系起来。数学上，这通常体现为连续方程的守恒形式，以及散度算子与坐标系几何因素之间的匹配关系。
2.1 连续方程的一般守恒形式
设某个物理量用场变量 \(q(x,t)\) 表示，其在时空中的演化满足守恒形式的通量方程。常见写法是 \[ \frac{\partial q}{\partial t}+\nabla\cdot \mathbf{F}=S, \] 其中 \(\mathbf{F}\) 为与 \(q\) 相关的通量，\(\nabla\cdot\mathbf{F}\) 是散度，\(S\) 为源项。若无源情形 \(S=0\)，则 \[ \frac{\partial q}{\partial t}+\nabla\cdot \mathbf{F}=0. \] 该结构的意义在于：局部的时间变化由通量在空间上的发散（净流入/净流出差异）以及内部生成或湮灭共同决定。
在更一般的情形中，\(q\) 可以是向量或多分量组（例如动量、能量、粒子分布），此时散度和源项都对应相应的分量结构，形成守恒型系统。
2.2 散度与通量在不同坐标系中的表达
散度算子刻画矢量场的“局部发散程度”，其具体计算依赖坐标系的度量与体元/面元表达。在直角坐标中，散度通常是简单的三方向偏导之和；在柱坐标、球坐标等曲线坐标中，散度包含额外的几何系数，从而保证守恒律在这些坐标系下依然与物理含义一致。
因此，“通量守恒”不是仅仅把 \(\nabla\cdot\mathbf{F}\) 写出来就结束，还要确保通量分量与度量因子配套。正确的几何表达保证了积分形式中边界通量与体内变化严格对应。
2.3 Reynolds 运输定理与从控制体到偏微分方程
控制体（control volume）方法是从工程与物理直觉直接出发的推导起点：考虑一个随时间可能运动或变形的体，分析体内量的变化与穿过边界的通量、以及体内源项之间关系。Reynolds 运输定理为这种“随动/控制体”的变化提供了数学桥梁。
当选择的控制体运动方式与物理过程相协调，Reynolds 运输定理可以把积分表达转化为局域的偏微分形式。最终得到的就是守恒律的一般形式：时间变化项 +（与通量有关的）空间发散 = 源项。这样推导的好处是逻辑闭环：从边界净流与内部源汇，严格导出局部方程。
2.4 积分形式与微分形式的等价性条件
守恒律可有积分形式与微分形式两种表达。积分形式直接反映控制体内量变化等于边界净通量加源项；微分形式则把上述关系“局部化”。两者等价通常需要一定条件，例如场变量足够光滑以便使用适当的微积分工具（如散度定理）把积分转为微分；或者在弱解框架下，等价关系按分布意义解释。
在数值与理论中，“弱化”要求往往是必要的，因为实际解可能出现间断（如激波）。此时，守恒律更自然地以积分形式保持其正确性，从而保证物理总量的平衡不因间断而丢失。
3 典型守恒律示例
守恒律在不同物理领域以不同变量与通量形式出现，但其共同点是：所选物理量在没有源汇时满足“局域净输运为零”的结构，或在有源时满足“源项补偿”的结构。
3.1 质量守恒（连续性方程）
质量守恒描述流体或介质中密度的演化。对密度 \(\rho\) 而言，连续性方程常写作 \[ \frac{\partial \rho}{\partial t}+\nabla\cdot(\rho \mathbf{u})=0, \] 其中 \(\mathbf{u}\) 是速度场。该形式意味着：密度在空间的变化来源于携带质量的流动穿过控制体边界的净效应。
在存在质量源或质量注入/消耗的模型中，右端会出现对应源项，从而体现“允许生成或湮灭”的一般扩展。
3.2 动量守恒（应力与体力的作用）
动量守恒通常与应力张量、压力梯度以及体力（如重力）联系。对不可压缩或更一般流体，动量方程可理解为：单位体积动量的变化等于应力带来的空间输运与体力的共同作用，形式上仍可写成守恒或准守恒结构：
- 通量项来自应力对面的“传递效应”；
- 源项（体力）表示内部受力导致的动量增减。
在固体力学中，应力与应变的关系决定了动量通量的细节；在流体中，黏性与压力共同构成应力贡献。尽管具体形式不同，其守恒的核心仍是“边界上通量的净作用对应体内动量变化”。
3.3 能量守恒（能量通量与功率项）
能量守恒讨论的是系统内部能量与其通过空间的输运之间的关系。能量的通量通常与速度场和压力/应力做功有关，并且可能包含额外的扩散或热传导效应（在扩散或黏性耗散模型中会以相应项体现）。一般守恒结构体现为：能量的时间变化由能量通量的发散和源项共同决定，其中源项可表示外部功率输入或体积内的能量转化机制。
当无外部功率且模型不含热源时，能量方程常呈现更接近纯守恒的形式；引入黏性耗散或热传导后，仍可以以扩展的守恒型方程描述，只是“耗散”往往表现为源项或扩展通量中的额外机制。
3.4 电荷守恒（连续方程的电磁对应）
电荷守恒是电磁理论中的基础要求。设电荷密度为 \(\rho_e\)，电流密度为 \(\mathbf{J}\)，则连续方程给出 \[ \frac{\partial \rho_e}{\partial t}+\nabla\cdot \mathbf{J}=0. \] 其物理含义是：电荷密度的增加只能来自电流造成的净流入；反之，电荷的减少对应净流出。
在更完整的电磁模型中，电荷守恒与麦克斯韦方程之间存在一致性关系：只要电磁场满足相应约束，电荷连续方程就应当在同一物理框架下成立。
4 边界条件与可积性/约束
守恒律不仅由体内方程决定，还强烈依赖边界如何处理。对守恒问题而言，边界条件决定了通量如何穿过边界，从而决定总量随时间的变化。对于数值计算与理论解，正确的边界设定是“守恒结构能否保持”的关键。
4.1 控制体边界上的通量指定
在控制体分析中，边界上通量可以采取不同形式：直接指定某种通量，或指定允许通过边界的状态量再由物理关系推导通量。对于守恒律，若边界上给出的条件导致通量净流入/流出被正确控制，那么控制体内量的变化将与方程保持一致。
因此，“边界条件”在守恒问题里不仅是补充信息，也是物理量守恒机制的一部分。错误的通量边界会造成总量误差，表现为非物理的增减。
4.2 初值问题：适定性与物理解读
对于以守恒形式表述的偏微分方程，初值条件通常给出初始时刻的场变量分布。适定性（解是否存在、是否唯一以及是否连续依赖初始数据）与方程的类型有关，例如线性/非线性、是否出现间断、是否为双曲型等。在许多守恒律场景下，方程具有双曲或双曲-抛物混合特征，初值问题的解读往往涉及特征传播、波动与间断的形成。
物理解读上，初值决定了初始量的空间分布，而通量与源项决定后续如何演化；若边界封闭且无源，初值对应的总量应在守恒意义下保持。
4.3 周期边界、封闭边界与开边界的差异
- 周期边界：边界“首尾相连”，净通量在全空间闭合意义上可为零，从而更容易保证整体守恒。
- 封闭边界：边界不允许物理量穿过或只允许特定形式的通量约束，闭合系统总量变化通常由内部源项决定。
- 开边界：允许通量通过系统边界进入或离开，系统总量会随外部环境输入输出而改变。此时“总量如何变化”更多依赖外部通量模型或边界数据。
在数值实现中，对不同边界类型的处理会直接影响守恒误差的累积方式。
4.4 分界面条件：跳跃与守恒的一致性
当解在空间上出现跳跃（例如接触间断、激波），分界面附近需要满足适当的跳跃条件（通常由积分形式推导得到）。在守恒律框架下，界面处的条件确保：尽管局部场量可能不连续，但控制体跨越界面时的净通量平衡仍成立，从而保证整体物理量不凭空改变。
因此，分界面条件本质上是对“局域守恒”在非光滑情形下的延续：用合适的弱形式把守恒结构继续保持。
5 从物理到“方程结构”：守恒形式的优势
把物理规律写成守恒形式，带来结构性的优势：它不仅反映物理含义，也在推导与计算中提供稳定性与可检验性。
5.1 为什么要用守恒型变量
守恒型变量指那些满足守恒结构直接出现在方程中的变量（或其等价表达）。使用这类变量的好处是：总量守恒往往在数学上更自然地体现为可对方程进行积分得到的约束，从而便于分析与验证模型。
当变量选择不当时，即使物理过程本质上是守恒的，方程可能被写成非守恒形式，导致在数值或分析中出现“守恒量被污染”的现象。使用守恒型变量能降低此类风险。
5.2 与非守恒形式（平流型）的关系与差别
非守恒形式通常表现为将某些项展开后，方程中不再直接是散度结构的形式。两者在光滑情况下可能通过代数恒等式相互转换，但在间断存在时，转换过程可能失效或引入额外项的歧义。由此，守恒形式更适合处理具有间断解的双曲守恒律问题。
从物理角度看，守恒形式强调“通过边界的通量平衡”，而非守恒形式更强调“局部场量的演化”。当解不光滑时，前者更能保证总量的物理正确性。
5.3 误差累积与守恒性质对数值稳定性的影响
数值离散不可避免地产生截断误差。若离散格式不满足守恒结构，误差可能在时间推进中系统性累积，表现为总量随时间漂移，即使其他局部指标看似合理。守恒型离散通过在离散层面控制通量，可以在一定程度上抑制这种漂移，从而提高长期计算的可靠性。
因此，守恒结构不仅是物理准确性的要求，也是数值稳定与误差管理的重要原则。
6 数值方法中的守恒型离散
守恒律的计算通常需要离散化。对于双曲守恒律而言，“守恒型离散”的思想尤为关键：离散格式应在离散控制体上保持通量平衡，从而尽可能复现连续层面的守恒。
6.1 有限体积法：天然匹配通量守恒
有限体积法（finite volume method, FVM）以控制体为基本单元，对每个网格单元积分守恒律。由于守恒律本质是“体内变化 = 边界净通量 + 源项”，这种积分直接对应离散方程中通量的计算与更新。
在这种框架下，同一条内部网格边界由相邻单元共享，因此只要数值通量在两侧取相反方向、保证一致性，整体上就能实现更严格的守恒。
6.2 通量计算与Riemann 问题的引入
在有限体积法中，更新单元内量需要知道单元界面处的通量。常见做法是把界面附近的局部状态近似为左右极限，并求解相应的局部初值问题，得到界面通量。局部初值问题的经典形式就是 Riemann 问题（左右两侧给定常值状态，界面初始为分段常值）。
Riemann 求解器或其近似方法能够提供界面通量，从而把物理传播（特征结构）嵌入到离散更新里。对于含间断的场景，这种做法通常比简单平均更符合真实传播机制。
6.3 保守性、单调性与通量修正
仅保证守恒并不足以避免数值振荡。对于非线性双曲守恒律，离散方案还需要处理高频误差，常用手段包括单调性保持策略、限制器（limiters）以及通量修正（例如避免非物理的新极值）。
从目标上看，守恒型离散负责“总量不漂移”，限制器或修正项负责“局部不产生虚假振荡”。两者结合通常能够在间断附近兼顾稳定与物理合理性。
6.4 离散误差如何影响总量守恒
即便采用守恒型格式，仍会存在与源项离散、时间推进误差、边界通量实现方式等有关的误差。尤其是：
- 源项离散若与通量离散不同阶或不一致，可能导致平衡误差；
- 时间积分若过于粗糙，可能破坏理想的守恒更新；
- 边界通量若处理不当，会将非物理输入输出引入系统。
因此，讨论数值守恒不能只看空间离散，还要关注全离散的一致性与误差阶数。
7 理论延伸与常见变体
守恒律的框架可延伸到含源模型、非线性激波问题、以及带扩散与混合机制的方程类型。同时，与弱解与测度等更一般的数学框架也存在直观联系：当解不光滑时，仍需以适当意义保持守恒结构。
7.1 含源守恒律与反应-输运类模型
含源守恒律在右端包含源项，源项常来自化学反应、相变、外场作用或粒子产生/湮灭等。反应-输运类模型通常把输运（带走）与反应（生成/消耗）耦合在一起。此时整体结构仍可以写成“守恒 + 源”，但反应项可能带来新的尺度与刚性问题。
在数值上，常见挑战是保持对平衡态的正确刻画（例如稳态解不被数值扰动破坏），这与离散源项和通量的匹配程度有关。
7.2 非线性守恒律与间断（激波）问题的基本图景
非线性守恒律方程中，通量取决于未知量本身。该耦合会导致特征线发生会聚，从而引发间断形成。物理上可对应激波或其他类型的波前。严格的解可能不再是处处光滑，而需要采用弱解概念来定义。
在弱解框架中，守恒律以积分意义表达，间断处通过合适的跳跃关系维持守恒一致性，从而确保物理总量不出现非物理偏差。
7.3 平流-扩散形式与“散度结构”保留的重要性
当加入扩散或黏性效应时，方程常呈现平流-扩散结构。扩散部分通常引入二阶导数项，但仍可通过适当重写保持散度结构（例如把扩散通量写成与梯度有关的矢量场，再取散度）。保留这种结构对边界条件与数值离散都很重要：它确保了扩散项对应的通量解释清晰，避免仅在代数层面展开导致守恒结构在数值上丢失。
7.4 与测度/弱解框架的联系（直观层面）
当解出现尖峰、浓度分布或其他分布性特征时，点值意义下的函数模型可能不够。直观上，弱解或测度框架允许未知量在更一般的数学对象中出现，从而仍能表达“任意测试函数意义下的守恒”。这种扩展使得守恒律能够覆盖更多非光滑现象，同时保持物理上“总量变化受通量与源项控制”的思想。
8 相关条目与阅读路径
这一部分给出进一步阅读的方向，帮助把守恒律与控制体方法、方程分类以及计算策略建立联系。
8.1 控制体方法、连续介质与流体力学
控制体方法贯穿守恒律推导，是理解“边界通量—体内变化”对应关系的直接途径。连续介质假设为把物理量表示为场提供基础，而流体力学中的质量、动量、能量守恒方程则是常见落地形式。
8.2 偏微分方程的分类：守恒律与相关方程
从偏微分方程角度，守恒律常与双曲型方程相关，非线性守恒律则与间断传播密切相连。理解方程类型有助于选择合适的初边值条件处理方式，以及判定间断、特征与数值稳定性的关键影响因素。
8.3 计算物理中的守恒型格式与误差分析
计算物理视角强调离散与误差：守恒型格式通过控制离散通量保持全局或局部守恒性质；误差分析则关注截断误差、源项匹配、时间积分与边界处理对守恒误差的影响。把物理结构与数值实现对应起来，能更系统地解释“为什么某些数值方案总量更稳”。