我知道你必须反映细胞穿过墙壁的速度并将反射的速度存储在幽灵细胞中（然后将用于通量/残差计算），但这是我理解的范围。我已经尝试在我的欧拉流求解器中实现这一点，但我似乎仍然遇到问题，这让我认为实现壁边界条件必须有更多。
如何在简单的 2d 欧拉流求解器中实现壁边界（使用幻像单元）？
计算科学
有限体积
2021-12-06 04:00:40
1个回答
指定速度边界条件是不够的。如果您确实在求解 2D Euler 方程，您还需要指定密度和温度的边界条件（如果您使用的是我假设的原始变量）。Maxim 是正确的，除了无滑移条件外
，您还需要强制垂直于壁的零速度。对于质量密度方程，您可以使用其中 0 是壁边界上内部单元的索引。
能量方程的墙 BC 需要更多信息。你的墙是绝热的吗？如果是这样，那么对于质量密度来说，鬼胞的实现类似于上面的那个。此外，您可以选择幻影单元中的温度值来表示通过壁的规定热通量。如果墙壁的温度由用户规定，则：
如果您使用一层幽灵单元，这一切都是有效的，但可以很容易地扩展到多层以用于高阶求解器。
其它你可能感兴趣的问题