CFD中文网

前辈，麻烦问一下，您的问题解决了吗？我刚接触这块知识，也是不会设置颗粒初始位置。（已知初始浓度）不知您有什么好的解决办法？

@youv 好的，我会去看看书上的内容，谢谢您的指点！

可以使用cyclicAMI边界条件

首先你需要做两套网格，可以要么用第三方网格软件来做，要么用openfoam自带的topoSet之类的来做，两个网格最重要的是，上图中的交界面要有边界。可以是wall，也可以是patch，反正后来要改。

导入到OpenFOAM之后，要把这个交界面设置为cyclicAMI，这个可以参考pimpleDyMFoam下面的propeller算例的设置过程

滑移网格是滑移和插值的过程。就不存在网格负体积了。你的运动发散是因为变形过大产生了负体积。

0_1528241930740_无标题.png

if (potentialFlow.finalNonOrthogonalIter()) { phi -= PhiEqn.flux(); }

这段我觉得倒不是因为phi被初始化了，好像上面求解的Phi方程就不是 http://dyfluid.com/potentialFoam.html 里头的方程（6），而应该是
▽(▽P)=▽(U-▽P)
也就是代码中的phi=(U-▽P)，这样子的话也能解释速度边界条件是怎么添加到Phi 方程中的了。但有一点不很明白，就是代码中好像表达的并不完全是phi=(U-▽P)的意思，所以还是有些糊涂。

另，

MRF.makeAbsolute(phi);//此处的makeAbsolute按字面意思是取绝对值，那这个对下文中重构速度场不会产生影响么？

这句话应该不是取绝对值，可能用了多重坐标系，之前有个变换，现在变回来了，但不是很确定是不是这个意思
@东岳 @张某人
上面是我的理解，如果有不对的地方还请直接指出来，一起讨论讨论