方程是：

$\partial f/\partial t + \nabla \cdot (v f) = 0$ $, \;\; f \in [0,1]$

是每个网格单元已知的速度。一个更精确的约束是要么是 0 要么是 1，但我也可以忍受它在域 [0,1] 中。几何是 2D 圆柱 RZ 坐标。$v$$f$

对于初始条件，域中的所有区域都是已知的，并且不是一个简单的函数。请注意，在时间和空间中不断变化，即实际上是。我不希望解决方案具有扩散性。$f(R, Z)$$R-Z$$v$$v$$v(R, Z, t)$