我尝试对一阶 Godunov 求解器实施约束传输方法，当我在 2D 转子测试用例上运行它时，求解器变得不稳定并崩溃。求解器能够保持无发散磁场，但它是高度振荡的，这使其不稳定。我发现问题在于使用算术平均面心磁场作为细胞中心磁场来找到通过黎曼求解器的通量。是否应该实施约束传输方法，使得算术平均面心磁场仅用作最后产生的磁场，并且在时间演化过程中不应使用它来代替细胞中心磁场？