我试图了解当使用欧拉方程的非线性 LU-SGS 算法时边界条件是如何实现的。大多数论文描述了对网格单元的 Gauss-Seidel 扫描，但没有解释边界条件是如何应用的。我可以考虑与显式情况下相同的 BC 实现（即扫描单元，然后在边界上添加通量校正），或者将由于边界条件导致的雅可比贡献添加到元素雅可比。然而，后者导致更复杂的实现，其中必须检查每个元素是否应该接收边界雅可比贡献。您是否知道任何更详细地讨论可压缩 Euler/Navier-Stokes 方程的 LU-SGS和边界条件的参考资料？

你认为有可能用更简单的方法代替“边界条件雅可比”方法吗？这会降低求解器的收敛速度吗？