Original RC interplation
No time dependency, No relaxation of the velocity field:

// Rhie-Chow interplation phi = (U_avg_f & mesh.Sf()) - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() );

Time dependency, relaxation of the velocity field:

// Rhie-Chow interplation phi = (U_avg_f & mesh.Sf()) - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() ) + (scalar(1) - URFU)*(phi.prevIter() - (U_avg_prevIter_f & mesh.Sf())) + DTf*( phi_old - (U_old_f& mesh.Sf()));

OpenFOAM
No time dependency, with or without relaxation of the velocity field (check simpleFoam solver):

surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());

Questions: Where is the second term of the original RC? Where is the third term of the original RC if taking the relaxation of the velocity field into consideration?

Time dependency, relaxation of the velocity field (check pisoFoam solver):

surfaceScalarField phiHbyA ( "phiHbyA", (fvc::interpolate(HbyA) & mesh.Sf()) + fvc::ddtPhiCorr(rAU, U, phi) );

Questions: ddtPhiCorr takes the function of (the third term+the fourth term of the original RC)? I guess NO. Let's recall the definition of ddtPhiCorr:

tmp<fluxFieldType> ddtPhiCorr ( new fluxFieldType ( ddtIOobject, rDeltaT * this->fvcDdtPhiCoeff(U.oldTime(), phiAbs.oldTime()) * ( fvc::interpolate(rA*rho.oldTime())*phiAbs.oldTime() - ( fvc::interpolate(rA*rho.oldTime()*U.oldTime()) & mesh().Sf() ) ) ) );

Look like ddtPhiCorr returns the following codes:

fvcDdtPhiCoeff*DTf*( phi_old - (U_old_f& mesh.Sf()));

Last question: How to take the effect of relaxation of the velocity field into ddtPhiCorr?
谢谢

meshToMeshInterpolate.C 可以看看这个内容，对，就是插值

@东岳 :xinlei:

@东岳
老师这是我目前算出来的一些结果：
链接: https://pan.baidu.com/s/163fY1uXF_sCCcN9nt1MiCQ
提取码: w3gp
网格量很少的，我现在只是想通过调试把路给走通，方便以后的使用

@东岳 谢谢老师，这个问题当时思维定式了，把U_p提出来，立马就懂了。

@D-Benjamin 求交集的意思

@小考拉 你好这个问题我和我师兄讨论过，你不可能输出瞬时的脉动速度场，因为湍流模型根本就没有求解脉动速度，只能输出二阶矩这种统计值，建模的话可以假设三个方向的速度大小相同，给出脉动速度的大小，但是不能给出速度的方向

@wwzhao

是的,OpenFOAM中面的方向定义永远是从小编号指向大编号，Gauss定理中面的方向指向外。也许把这些面分为大面和小面区分一下更容易理解。

对于编号大于$P$的单元$U$,他们之间面上的量为:

$$\phi_ f = \varpi \phi_{P} + (1 - \varpi) \phi_{U}$$

$$\dot { m } _ { f } = ( \rho \mathbf { u } ) _ { f } \cdot \mathbf { S } _ { f }$$

对于编号小于$P$的单元$L$,他们之间面上的量为:

$$\phi_ f = (1 - \varpi) \phi_{P} + \varpi\phi_{L}$$

$$\dot { m } _ { f } =- ( \rho \mathbf { u } ) _ { f } \cdot \mathbf { S } _ { f }$$

$$
\sum _ { f \sim n b ( P ) } ( \rho \mathbf { u } \phi ) _ { f } \cdot \mathbf { S } _ { f }
= \sum _ { N \in L ( P ) } \left( - \dot { m } _ { f } \phi _ { f } \right) + \sum _ { N \in U ( P ) } \dot { m } _ { f } \phi _ { f }
$$
$$
= \sum _ { N \in L ( P ) } - \dot { m } _ { f } \left[ ( 1 - \varpi ) \phi _ { P } + \varpi \phi _ { N } \right] + \sum _ { N \in U ( P ) } \dot { m } _ { f } \left[ \varpi \phi _ { P } + ( 1 - \varpi ) \phi _ { N } \right]
$$

$$
= \left( \sum _ { N \in L ( P) } - \dot { m } _ { f } ( 1 - \varpi ) + \sum _ { N \in U ( P ) } \dot { m } _ { f } \varpi \right) \phi _ { P } + \sum _ { N \in L ( P ) } - \dot { m } _ { f } \varpi \phi _ { N } + \sum _ { N \in U ( P ) } \dot { m } _ { f } ( 1 - \varpi ) \phi _ { N }
$$

这样可以看出lower存的是 $- \varpi _ { f } \dot { m } _ { f } $, upper存的是$\dot { m } _ { f } ( 1 - \varpi ) _ { f } $.并且diag存的是
$$a _ { P } = -[\sum _ { N \in L ( P ) } \dot { m } _ { f } ( 1 - \varpi ) + \sum _ { N\in U ( P ) } (-\dot { m } _ { f } \varpi )]$$

写个类调用下fftw的函数接口就好了

@东岳 您好，东岳老师。我最近才看到您之前发的这个动图。这个是用mppicinterfoam做出来的吗？我尝试了很多次都不能堆积出您的这种沙床的结果。您能简单说一下您怎么做的吗？等您回复！

http://www.cfd-china.com/topic/3182

有一些介绍，你试试二楼的

@赵一铭

感谢分享！最近安装低版本的OpenFOAM做原有的程序二次开发，用到了这个语句~~，

在此分享一个基于不同Linux系统的OpenFOAM不同版本安装总结，亲测OF211，OF231 再 Ubuntu 16.04 安装成功无Error！

https://openfoamwiki.net/index.php/Installation/Linux

再次感谢！

@东岳 谢谢老师，太久没看，给忘了

@Cp_Zhao 明白了，感谢指点！！

@gyzhangqm 是否有足够的数据支撑？有的话可以。

这是div(phi,T)与-div(phi,T)等的解释

微信图片_20200320103052.jpg

没有simpleR，http://dyfluid.com/simplefoam.html

可以直接在of中sample提取数据

@Samuel-Tu Lib那改到Foam_USER_Libbin

1.JPG
2.JPG