Skip to content

数值通量

Numerical flux · Conservative interface flux · 数值界面通量

以左右离散状态给出共享界面通量,并用一致性、单调性与耗散量区分不同守恒格式。

条目类型
定义

形式陈述

标量守恒律 ut+f(u)x=0,两点数值通量是函数

f^:R×RR,(uL,uR)f^(uL,uR),

用来近似界面处沿正向法向流出的物理通量。最基本的一致性条件是

f^(u,u)=f(u).

它保证常状态不会被离散更新改变,并使光滑极限中的通量回到原方程。收敛分析通常还要求 f^ 在有界状态集上连续或 Lipschitz;“守恒”则来自相邻单元共享同一个 f^,不是来自函数本身的某个单独代数标签。

标量两点通量若对左状态 uL 单调不减、对右状态 uR 单调不增,称为单调通量。它是完整更新单调的重要条件,但还不够:时间步也必须使更新对中心状态的系数非负。

常见耗散形式为

f^(uL,uR)=f(uL)+f(uR)2α2(uRuL),

其中 α 至少覆盖相关状态区间内的最大波速时,得到局部 Lax–Friedrichs/Rusanov 型通量。第一项居中,第二项抑制跳跃;α 太小可能失去单调性,过大则把接触和激波涂抹得更宽。

Godunov 通量定义为局部 Riemann 熵解在界面射线 ξ=0 的物理通量。对标量通量可写成全区间极值公式

f^G(uL,uR)={minu[uL,uR]f(u),uLuR,maxu[uR,uL]f(u),uL>uR.

线性通量 f(u)=au 则给 f^=auLa>0,给 f^=auRa<0,即迎风通量。

边界数值通量的两个参数不总是“两个内部单元平均”。一个参数可由内部迹给出,另一个由流入数据、反射条件或边界 Riemann 问题构造。法向反转时通量符号和左右状态约定必须一致变换。

直觉

数值通量是相邻单元之间唯一的交易凭证。左右单元可以各自拥有不同重构值,但界面最终只能确认一个穿越量;若双方各算一个近似物理通量,两个账本不会精确相消。通量函数把连续波在一个时间步内的复杂界面行为压缩成这个共同数字。

不同通量主要在“对尚未解析的波做多少假设”上分岔。Godunov 通量解局部波问题,尽量按真实传播方向取值;中心加耗散通量只需波速上界,成本低而较模糊;迎风通量在线性单速问题中恰好精确表达方向。高阶重构提升送入通量函数的左右迹质量,并不会自动把一个不稳定通量变稳定。

例子与边界

对 Burgers 通量,取 uL=2,uR=0。Riemann 解是速度 s=1>0 的激波,界面 ξ=0 位于左状态侧,所以

f^G(2,0)=f(2)=2.

若使用耗散参数 α=2 的 Rusanov 通量,则

f^R(2,0)=2+0222(02)=3.

两者都一致,却对非相等状态给出不同通量。Rusanov 式中的居中项为 1、耗散项贡献为 2,所以总通量为 3;它比 Godunov 通量多 1,但不表示真实物理通量在任何状态上等于 3

再取线性通量 f(u)=3u,左右状态为 (1,4)。速度为负,信息从右侧进入界面,所以迎风通量是

f^up(1,4)=34=12.

若误取左态会得到 3,方向正好相反。这个错误在常状态测试 f^(u,u)=f(u) 中仍可能通过,说明一致性测试不能替代方向和稳定性测试。

非凸标量通量不能把上述全区间极值偷换成只比较端点:状态区间跨过多个极值时,熵 Riemann 解可能是复合波,但全区间极值公式仍给出精确 Godunov 通量。系统通量是向量,Roe、HLL 与 HLLC 等近似求解器保留不同波族,不能由一个标量 α 的优劣概括全部接触分辨率。

边界通量也会决定全局总量。周期网格最后一个单元与第一个单元必须共享同一通量;流出边界若用固定外部状态,可能人为送回反射波。守恒检查通过只说明内部成对抵消,边界物理是否正确仍需独立验证。

推论与应用

一个通量可从四个互不替代的角度审查:一致性检查常状态,守恒装配检查共享界面,单调性或熵稳定性控制非物理振荡,精度分析检查光滑状态下的展开。数值耗散常改善稳定与熵选择,却会降低间断分辨率;减少耗散则可能暴露高频振荡。

Godunov、Lax–Friedrichs 与迎风格式的差别可以统一归结为所选界面通量及时间更新。MUSCL 等高阶方法仍调用同一两态或近似 Riemann 通量,只是两态来自受限重构。实现时适合把通量当作纯函数测试:核对 f^(u,u)、交换法向后的符号、给定端态的手算值以及非有限输入的失败行为。

参考资料
  • 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.
关系图谱12 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系