1 离散守恒的基本概念
1.1 守恒量与守恒律的离散化
离散守恒是指把连续模型（通常是守恒律形式的偏微分方程）离散成可计算的网格方程或更新格式后，所对应的某种物理或几何量在计算演化过程中不发生“系统性漂移”，理想情况下严格不变，实际工程中在给定误差容限内维持守恒。守恒量可能对应质量、动量、能量，或更一般的“通过守恒律定义的量”（例如标量场的某种整体量、某类积分不变量等）。
离散化并不只是把方程换成代数方程；关键在于离散格式应保留连续方程所蕴含的结构，例如“通量进出抵消”“局部守恒导致整体守恒”等性质。若离散化过程引入了等价于“源项/汇项”的数值误差，就可能出现守恒破坏。
1.2 从连续到离散：关键映射思想
连续守恒律常写作“时间变化等于通量的散度（或其负值）”，其积分形式揭示了区域边界上的流入流出平衡。离散守恒的核心映射思想，是选择一种离散算子与离散几何/网格结构，使连续的积分恒等式在离散层面仍可成立。实践中通常体现为以下对应关系：
- 连续的体积分与离散网格单元的加权求和相对应；
- 连续的边界通量与离散面通量相对应；
- 连续的散度与离散散度（或差分散度）相对应；
- 连续的时间导数与离散时间推进格式相对应。
当这些映射建立了恰当的“结构同构”或“配对关系”时，离散系统才能在数值层面复现守恒律的抵消机制。
1.3 “保持不变”与“数值误差”边界
“保持不变”并不等同于零误差。离散守恒可被理解为两类目标：
- 结构性守恒：在格式的离散代数层面，守恒误差不会产生与时间单调增长相关的系统偏差，或其增长被控制在更高阶或与网格尺寸相关的范围内。
- 误差容限内守恒：考虑浮点误差、截断误差与迭代误差，实际计算结果在允许的误差界内满足守恒指标。
边界的划分通常取决于模型类别、网格拓扑与离散格式的设计。例如有些方法保证局部单元守恒，从而自动给出全局守恒；也有些方法保证能量或范数的“稳定性形式”，其意义相对偏向“结构保持”。在严格意义上，离散守恒往往与“格式满足离散积分恒等式”或“数值源项为零（或可严格界定）”相关。
2 数学表述与离散算子
2.1 通量守恒形式（离散积分恒等式）
许多守恒律的连续形式可转写为积分守恒：对任意控制体，量的时间变化等于边界通量的净流量。离散守恒通常以“单元级”的通量平衡实现：每个离散单元在一步更新中只通过其面（或边）的通量与邻接单元交换，从而使得内部面通量在全局求和中相互抵消。
因此，离散积分恒等式的典型形态是：把时间推进写成单元量的更新，并把通量通过单元界面上的离散通量相加。若该结构保持得当，全局守恒可由“内边界抵消 + 外边界给定”自然获得。
2.2 离散差分与离散散度/梯度的配对
离散算子之间的“配对”是离散守恒能否成立的关键。连续情形中存在诸如分部积分、对偶性与格林恒等式的结构；离散版本若能保留这些关系，则更可能实现守恒或能量一致性。
常见的做法是构造离散散度与离散梯度，使它们之间满足某种离散对偶关系（例如离散散度是离散梯度的伴随算子，或满足相应的离散分部积分公式）。这种配对不仅影响守恒，还影响稳定性：若离散算子破坏了连续结构，数值方案可能在每一步引入等价于源项的残差，导致守恒漂移。
2.3 离散守恒与离散连通性（如网格拓扑）
守恒是“局部交换、整体抵消”的结果，因此离散连通性和拓扑结构直接影响守恒实现。网格拓扑决定了单元如何共享面、面如何共享边，从而决定通量能否在离散层面正确抵消。
例如：
- 若网格与通量未知量放置位置不匹配，可能导致内部面不能精确抵消；
- 若存在不规则连接（例如非流形边界、粗细网格接口），需要额外的接口守恒策略；
- 对周期性边界，离散拓扑应正确处理“首尾相连”，否则守恒会被边界处理误差破坏。
因此，离散守恒并非仅是局部公式正确，更依赖全局网格拓扑与变量布置方式的一致性。
3 数值方法中的实现路径
3.1 有限体积法与通量守恒
有限体积法以控制体思想为出发点，将守恒律直接在单元上积分。它通常把未知量存于单元中心或单元平均值，把通量存于面，并通过通量通用性保证内部面抵消。由于有限体积天然对应“离散积分恒等式”，因此常被视为实现离散守恒的标准路径。
在实际算法中，关键环节包括：面通量的数值近似（由重构与黎曼求解或通量函数给出）、时间推进的稳定性选择，以及边界通量的正确施加。若面通量计算在离散层面满足一致的符号约定和共享关系，则守恒往往可在全局层面得到保证。
3.2 有限差分格式与守恒差分形式
有限差分方法可通过把方程改写成差分的守恒形式来实现离散守恒。具体地，把时间变化写成离散的通量差：单元（或网格点附近的控制体）量的变化由相邻界面通量之差决定。只要这种“通量差”结构被严格采用，并保证相邻格点之间通量使用一致，就能避免非物理源项。
有限差分的难点通常在于非均匀网格、边界处理与高阶格式下的离散算子配对：若高阶重构导致通量不一致或破坏对偶结构，守恒性质可能下降到更低阶或出现局部漂移。
3.3 有限元/谱方法中的结构保持
有限元与谱方法并非总以“控制体通量”形式呈现，但通过选择合适的弱形式、基函数空间与离散算子，仍可实现某类守恒或结构保持。例如，在弱形式中，守恒律对应的对偶结构可通过质量矩阵、刚度矩阵与数值积分方式得到体现。
在工程实践中常见的策略包括：
- 选择与守恒律相关的测试函数空间，使离散对偶关系成立；
- 使用适当的数值通量或边界项处理，以匹配连续的积分抵消；
- 在混合有限元或离散外微分结构中，通过几何一致性实现守恒。
需要注意的是，有限元/谱方法有时更强调“能量稳定”“范数保持”等替代指标；这些指标与守恒可能密切相关，也可能在特定情形下形成“近似守恒”。
3.4 守恒约束与投影技术（结构保真）
当原始离散格式未自动满足某个守恒律时，可以通过约束或投影方式修正数值解，使其回到守恒流形附近。这类方法通常把“守恒量”写成约束方程，随后对更新结果进行最小改动或投影，使守恒误差被消除或缩小。
投影技术常见形式包括：
- 在每个时间步后对全局守恒量进行校正；
- 使用拉格朗日乘子把守恒约束内嵌到离散系统；
- 对局部守恒量采用守恒投影，确保单元或面级别的平衡。
这类方法的优点在于实现直接、效果可控；其潜在代价是可能影响高阶精度或引入额外噪声，因此通常需与误差分析与稳定性评估共同考虑。
4 结构性条件与充分性准则
4.1 一致性条件：局部到全局的匹配
离散守恒往往需要“局部结构正确”并能传递到“全局守恒”。一致性条件强调：每个单元的离散更新必须与共享界面的离散通量一致，保证内部贡献抵消；同时边界条件必须以与连续模型一致的方式进入离散通量。
常用的充分性线索包括：
- 内部面通量的对称性（同一通量以相反符号贡献两侧单元）；
- 单元体积与界面面积的几何离散一致；
- 对源项为零的纯守恒模型，离散方程不应产生等价源项。
当这些匹配被满足时，局部平衡即可推出全局守恒，从而形成可验证的结构性准则。
4.2 稳定性与守恒的协同（能否同时做到）
守恒与稳定性经常是协同目标，但并非总能自动同时满足。许多守恒格式在理论上可保证守恒，却可能在某些时间步长或离散参数下出现振荡或数值不稳定；反之，一些追求强稳定性的格式可能会引入额外数值耗散，从而影响某种“严格守恒”的形式。
因此需要考察两类关系：
- 守恒对稳定的作用：守恒本身不等价于稳定，然而守恒往往限制了某些能量或整体量的无界增长。
- 稳定对守恒的影响：为抑制振荡采取的滤波或耗散可能会破坏守恒，除非耗散被设计为“守恒兼容”（例如通过结构保持或守恒修正实现）。
在更一般的框架下，可以把守恒与稳定性视为同一离散结构中的不同面向：通过合适的离散算子与时间推进策略，使二者同时满足或在误差阶上协调。
4.3 误差分析中的守恒主导项
误差分析通常包含截断误差、离散几何误差与数值求解误差。离散守恒强调一种“误差主导项”的观点：如果守恒结构被正确保留，守恒量的误差不再由一阶或低阶的数值源项主导，而可能由更高阶的截断项或与网格尺寸相关的项主导，从而表现为守恒误差随网格细化而衰减。
在分析中常出现以下现象：
- 非守恒源项会导致守恒量的误差在时间上累积，表现为漂移；
- 结构兼容的格式会使内部抵消成立，漂移被抑制到可控阶次；
- 即便有耗散，若其形式仍满足守恒律所要求的约束结构，则“能量/物理量的衰减机制”可以被解释为物理或数值中的合理效应，而非凭空生成。
5 收敛性与物理一致性
5.1 守恒性对收敛的影响
守恒性与收敛性存在关联，但关系并非一一对应。守恒保证了某些整体量不发生不物理增长，有助于建立先验界，从而为收敛证明提供基础条件。尤其在非线性或激波等情形中，全局约束可减少“离散解跑偏”的可能。
同时也要认识到：守恒并不必然保证收敛到正确的弱解；若离散格式在通量函数、重构或边界处理上不一致，可能仍出现偏差。因而，守恒往往是“必要的结构条件之一”，但不是唯一充分条件。
5.2 长时间积分下的误差增长控制
许多数值方法在短时间内误差可接受，但在长时间积分中出现累积效应。离散守恒的优势在于：当守恒量对应的误差不呈持续积累，长时间计算更可能保持物理一致。
常见机制包括：
- 守恒结构限制了整体量的漂移，使误差被“困在”稳定吸收域或受控范数范围；
- 若同时保持能量相关结构，还可减少相位误差导致的能量泄漏；
- 时间离散若满足相容性（例如与守恒结构一致），则可避免因时间离散造成的逐步源项。
5.3 避免非物理源项：数值“凭空产生/消失”的问题
离散守恒的直接物理意义在于避免“守恒律被数值破坏”。当格式不满足离散积分恒等式时，可能产生等效源项，使质量、能量或其他量在没有物理原因的情况下被生成或消除。此类问题在工程观测上常表现为：
- 质量不守恒：例如流体计算中密度或体积积分随时间变化；
- 能量不守恒：例如波动问题出现逐步衰减或增长；
- 电磁计算中约束违反：例如由守恒律导出的约束量偏离。
因此，守恒不仅是抽象数学性质，也是一种对物理一致性的底线约束。
6 典型应用场景（非对立性综述）
6.1 流体与守恒律（质量、动量、能量）
在可压缩或不可压缩流体模型中，连续方程常体现质量守恒（不可压缩的散度约束或可压缩情况下的守恒形式）、动量守恒与能量守恒。离散守恒方法通过通量平衡实现这些量的整体一致，特别在包含激波、接触间断或强对流时，保守型格式能降低非物理振荡导致的守恒破坏。
不同量的守恒可能需要不同的数值细节：质量和动量通常可通过面通量守恒框架实现，而能量还涉及方程组耦合、压力功与数值耗散的协调。
6.2 热传导与能量守恒的离散处理
热传导问题常以扩散型守恒律为主，热通量的散度决定温度的变化。离散守恒强调：在每个控制体上，热量的变化应由界面热通量之差产生。只要离散导数与通量的几何一致性匹配，整体热量变化会与边界条件相符。
在含有变系数介质或复杂边界的情形，守恒性的关键落在网格度量与通量离散的配对上，以避免“数值源项”造成的温度整体漂移。
6.3 计算电磁学中的离散守恒观念
计算电磁学涉及电荷守恒与麦克斯韦方程相关结构。离散守恒在此常表现为对电荷守恒约束或相关散度条件的保持。若离散格式不能维持这些约束，可能出现约束违反并向时间推进传播，导致数值解偏离物理行为。
因此，电磁离散守恒常与“约束保持”“结构保持”的思想紧密结合：通过合适的离散微分结构或网格上对应的对偶算子，使散度与旋度等算子之间的关系在离散层面保持，从而避免不符合守恒律的误差增长。
6.4 图模型/网络中的流量守恒类比
离散守恒的思想可迁移到图模型与网络流：在每个节点上，流入与流出满足平衡，形成守恒约束。虽然这与连续介质的几何意义不同，但形式上仍体现为“局部交换—全局平衡”的结构。
在这种类比中，“网格拓扑”被图的连通关系替代，“通量”由边上的流量变量或概率流表示。保持流量守恒可以用于解释稳定性、守恒损失在学习中的正则项，或用于确保网络上迭代过程不会无源增殖。
7 离散守恒的评价指标与诊断
7.1 守恒误差的定义与测量
评价离散守恒通常需要明确守恒量的离散定义与误差度量方式。常见指标包括：
- 全局守恒误差：守恒量在计算终止时刻与初值的差；
- 相对误差：误差除以初值或特征量，便于跨尺度比较；
- 随时间的轨迹：记录误差随步数的变化，用于判断漂移还是受控振荡；
- 局部守恒残差：每个单元或节点的通量不平衡，用于定位问题来源。
这些指标共同帮助区分“误差来自截断”“误差来自边界处理”“误差来自离散算子不配对”。
7.2 结构破坏的常见征兆（漂移、涡流式误差）
结构破坏常见表现并不只是一条曲线的误差变大，还包括空间分布上的异常：
- 单调漂移：守恒量随时间单方向改变，提示存在系统性的等效源项；
- 局部残差反复正负抵消不足：全局看似不大，但局部长期偏离会影响后续演化；
- 涡流式误差：在对流主导问题中，离散不一致可能导致数值“回流”或错误传播，形成类似涡动结构的误差组织方式；
- 边界相关的集中误差：在复杂边界或接口处残差显著，通常说明通量一致性或边界项离散不匹配。
7.3 回退测试：守恒回归与回归曲线（梗式：别让能量“跑路”）
回退测试是指构造容易验证守恒性的基准问题，用来检查离散格式的结构是否“回归正确”。典型做法包括：
- 选择纯守恒模型的解析解或构造满足守恒的初值与边界条件；
- 做网格加密实验，观察守恒误差是否随网格尺寸按预期衰减；
- 对比不同格式或不同参数下的守恒轨迹，以判断是否存在结构破坏。
梗式表达“别让能量‘跑路’”，对应的就是观察能量或其他关键守恒量是否出现不合理逃逸：例如在无耗散物理模型中能量却持续衰减或增长，往往意味着格式破坏了守恒结构或出现与时间推进不相容的误差累积。
8 相关概念与扩展方向
8.1 边界条件对守恒的影响
边界条件是离散守恒成败的重要因素。即使内部通量抵消正确，若边界通量与连续模型不一致，仍可能导致守恒误差。常见挑战包括：
- 边界通量的数值一致性（符号约定、法向方向等）；
- 复杂几何下的边界离散误差；
- 周期边界、滑移边界、吸收边界等类型的实现与守恒相容性。
因此，边界处理不仅是“补上信息”，还需要保证其与离散积分恒等式在同一框架下成立。
8.2 离散对称性、Noether思想的离散版本
Noether思想表明连续对称性与守恒律之间存在对应关系。离散守恒的一个扩展方向是研究：当连续系统的某种对称结构被离散化保留时，是否能得到相应的离散守恒量。与此相关的概念包括结构保持、变分离散与几何数值积分等，它们希望在离散层面保留“对称性—不变量”的对应，从而提升长期稳定性与物理一致性。
8.3 变网格与自适应中守恒的维护
自适应网格细化/粗化会改变控制体与界面几何，从而影响通量抵消。要在变网格中保持守恒，通常需要：
- 进行守恒插值或守恒重映射，使守恒量在网格转换时不被破坏；
- 在粗细接口处采用一致的通量守恒策略；
- 处理网格运动或动态边界问题，保证离散几何更新与通量计算相互兼容。
缺乏这些处理时，网格更新步骤可能成为守恒破坏的主要来源。
8.4 随机离散与期望意义下的守恒
8.4.1 统计守恒与波动幅度控制
在随机离散或随机驱动系统中，精确守恒往往不再适用，但可讨论“期望意义下”的守恒：例如守恒量在统计平均意义上保持，或其偏差在概率意义上受控。评价重点从确定性误差转向统计量控制，包括：
- 期望守恒误差随时间的增长；
- 方差或高阶矩随网格与时间步的变化；
- 随机噪声与离散结构之间的相容性，避免噪声被结构性地放大为系统偏差。
在这种框架下，“守恒”更多体现为统计平衡与波动受限，而非每次样本路径都严格守恒。