请问您能帮我解决使用标记和单元方法（MAC）编写的斯托克斯流的问题吗？我只需要解决eq。连续性+动量方程。对于给定的条件（稳态）。我不明白如何在边界上做事情（假设我想在边界上设置 v=0），如何填充矩阵 - 我将把它解决为 Ax=b 其中 A 是矩阵 3*N*Nx3*如果我的盒子是 NxN 并且我有 3 个未知数 (vx,vy,p) 和三个 eq，则 N*N 和 b 右手边向量的昏暗 3*N*N。（连续性当量 + 2 动量当量）或者有什么代码我可以看一下以了解它的味道吗？非常感谢
标记和单元方法 (MAC) - STOKES FLOW - 边界？
计算科学
边界条件
矩阵
2021-12-08 06:55:08
1个回答
Harlow 和 Welch 的原始出版物讨论了这一点。如果您在您最喜欢的搜索引擎中搜索“Harlow Welch 标记单元格”，您应该会获得一些副本链接。如果这没有帮助，我自己的 Int 文章。J.数字。冰毒。标题为“Stokes 方程和 MAC 方案的无发散不连续 Galerkin 方案”的流体，可在http://www.isc.tamu.edu/publications-reports/tr/0705.pdf的预印本中获得，可能会有所帮助。
对于条件 u dot n = f，您将边界上的自由度设置为 f 的值。
对于切向速度的条件，您可以应用惩罚公式，例如，对于源自边界的边缘上的所有自由度，该公式可以简化为 alpha (uf) = 0 形式的附加项，但不是边界。