我需要实现一个 3D 版本的 Corner Transport Upwind (CTU) 有限体积方法(在 python 中);所以我一直在阅读 Leveque,“双曲问题的有限体积方法”,我认为这非常好。下面显示的是关于如何将线性平流的 CTU 推广到任意速度场(空间中的非常数)的关键部分。
我不了解通量更新和. 认为是负数并且是正的,但它们的大小相同。然后在两者之间都有横向通量和,因此这些通量更新将非零。进一步,将是正数(符号由于速度的符号而翻转,对吗?)和将是负面的。
然后将这些单元格用于单元格的通量差分公式作为:. 因为这两个项之前有相反的符号,所以减号使它们加在一起,而不是取消。
这对我来说似乎不对。关于细胞,我认为这些横向术语是为了解释不流入细胞的部分波,但在这种情况下,供体细胞波从进入在顶部和底部均等地伸展,所以不应该流入就好像横向通量根本不存在一样,就好像和都是零?也就是说,供体细胞流入扫出同一个区域,仿佛和都是零。