“其中数值通量在有界集上连续并一致:”
形式陈述 ​
对标量守恒律
用来近似界面处沿正向法向流出的物理通量。最基本的一致性条件是
它保证常状态不会被离散更新改变,并使光滑极限中的通量回到原方程。收敛分析通常还要求
标量两点通量若对左状态
常见耗散形式为
其中
Godunov 通量定义为局部 Riemann 熵解在界面射线
线性通量
边界数值通量的两个参数不总是“两个内部单元平均”。一个参数可由内部迹给出,另一个由流入数据、反射条件或边界 Riemann 问题构造。法向反转时通量符号和左右状态约定必须一致变换。
直觉
数值通量是相邻单元之间唯一的交易凭证。左右单元可以各自拥有不同重构值,但界面最终只能确认一个穿越量;若双方各算一个近似物理通量,两个账本不会精确相消。通量函数把连续波在一个时间步内的复杂界面行为压缩成这个共同数字。
不同通量主要在“对尚未解析的波做多少假设”上分岔。Godunov 通量解局部波问题,尽量按真实传播方向取值;中心加耗散通量只需波速上界,成本低而较模糊;迎风通量在线性单速问题中恰好精确表达方向。高阶重构提升送入通量函数的左右迹质量,并不会自动把一个不稳定通量变稳定。
例子与边界
对 Burgers 通量,取
若使用耗散参数
两者都一致,却对非相等状态给出不同通量。Rusanov 式中的居中项为
再取线性通量
若误取左态会得到
非凸标量通量不能把上述全区间极值偷换成只比较端点:状态区间跨过多个极值时,熵 Riemann 解可能是复合波,但全区间极值公式仍给出精确 Godunov 通量。系统通量是向量,Roe、HLL 与 HLLC 等近似求解器保留不同波族,不能由一个标量
边界通量也会决定全局总量。周期网格最后一个单元与第一个单元必须共享同一通量;流出边界若用固定外部状态,可能人为送回反射波。守恒检查通过只说明内部成对抵消,边界物理是否正确仍需独立验证。
推论与应用
一个通量可从四个互不替代的角度审查:一致性检查常状态,守恒装配检查共享界面,单调性或熵稳定性控制非物理振荡,精度分析检查光滑状态下的展开。数值耗散常改善稳定与熵选择,却会降低间断分辨率;减少耗散则可能暴露高频振荡。
Godunov、Lax–Friedrichs 与迎风格式的差别可以统一归结为所选界面通量及时间更新。MUSCL 等高阶方法仍调用同一两态或近似 Riemann 通量,只是两态来自受限重构。实现时适合把通量当作纯函数测试:核对
参考资料
- Randall J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002, Chs. 11–12.
- E. Godlewski and P.-A. Raviart, Numerical Approximation of Hyperbolic Systems of Conservation Laws, Springer, 1996, Ch. 3.
- Eleuterio F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics, 3rd ed., Springer, 2009, Chs. 10–12.