CFD中文网

对粒子并行，每个cell生成concurrent vector，用粒子坐标计算出所在网格，pushback即可。

@东岳
之前忘了@了，麻烦东岳老师有空的时候再帮忙看下哈，谢谢！

感谢李老师的解答！我后来推出来了。关于同位网格的simple算法还有个小细节想请教您。

下图中的(6)中的计算界面速度修正值，我个人感觉这一步是不是没有必要？因为第五步已经计算出了压力修正方程，可以求出各个节点的压力修正值，进一步可以得出个节点的压力值，然后插值便有了界面压力值，可以通过6-72a和6-72b计算节点的速度修正值？最后更新节点的速度和压力值即可？

近日在读数值传热学和Versteeg的计算流体力学时发现一些疑问，想请教下各位。

在数值传热学247页的6-67式中写道

而在Versteeg的计算流体力学339页中有

感觉这两个式子怎么也不一样，请问这两个式子到底等价吗？为什么？

另外如果使用英文书中的这种格式，那么离散动量方程的对流项时并线性化后，会出现$$U^n S U^r $$，n为当前步，r为下一步。在这里是要将U^n和U^r的界面速度全都替换为11.83的格式吗？还是只替换U^r的？？

@东岳 !
下图是我从您说的Peric的书里截出来的
我的理解是7.55,7.56式是用了第n步的压力先求出一个中间的速度场$$u^{**}$$。再通过对7.57式求散度，利用连续性方程消除$$u^{n+1}$$,使之成为关于压力的泊松方程，再对压力进行求解，最后返回来再求$$u^{n+1}$$.
不过请问东岳老师为什么7.55式子中有第n步的压力梯度，而7.56式中又消掉了一半呢？请问这个是什么原理？

本人是19Fall机械工程还没入学的新生PhD，方向是流固耦合。
我看课题组里的论文多是用immersed-boundary method处理流体部分，用finite element method处理固体部分。据我目前上网查的资料所知，OpenFOAM似乎是主要用有限体积法处理流体问题的，想请教下各位使用OpenFOAM的前辈，这样研究方向是否适合使用Openfoam？