matlab优化微分方程组代码3D
Navier
Stokes方程求解器
使用有限差分法和均匀网格并行求解3D不可压缩Navier
Stokes方程。
不可压缩性是通过使用压力芯实现的方案和线性(泊松)求解器是使用多网格v周期实现的。
该代码几乎是基于此的。
该文档位于[此处](se
/
codes
/
mit18086_navierstokes.pdf)。
问题陈述
不可压缩的Navier
Stokes方程:
不可压缩条件:
3D立方域[0,1]
^
3。
实现了速度的Dicirhlet边界条件和压力的Neumann边界条件。
边界上的速度可以用双bcs
[3]
[6]数组指定(第一个维度指定x-y-z-速度,第二个维数指定立方体的面)。
执行
立方域通过不规则的交错网格离散化。
偏微分方程采用有限差分法离散化。
并行化由OpenMP实现。
压力校正方案用于增强不可压缩性。
用法
./multigrid
[线程数]
[v周期的最大级别]
[x网格大小]
[y网格大小]
[z网格大小]
2021-07-24 17:28:26
29KB
系统开源
1