三维不可压流体求解taichi代码

各位同学老师好,我用taichi写了一个三维的,时间空间二阶精度,基于交错网格, 采用两步映射法求解NS方程的流体求解器,
https://github.com/MsureCFD/FVM_Flow3D_taichi
想随着taichi数值计算兴趣小组的不断壮大,将其优化为一个良好的代码从而应用于工程领域的数值模拟。因为业余时间有限,当前该代码的边界条件还没设置太多,周期性边界试验了一下有点小问题,但是无滑移和自由滑边界都是正常的。先用标准的cavity flow验证了一下,结果没问题。周期性边界修正好后可以拿来直接数值模拟 (quasi-DNS) 算一些紊流case比如channel flow去解析紊流。当应用于紊流直接数值模拟时,case对并行的要求极高,这也许是taichi可大展拳脚的部分。目前这个代码写的比较丑陋且有极大可优化空间,希望随着对taichi的不断学习可以将其逐渐优化改进。后期希望添加vof和浸入边界法,进行一些波流与结构物相互作用的数值模拟。

5 Likes

赞~~期待优化后的成果~