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?
谢谢
@Cp_Zhao
好的，感谢您及时的回复哈！
我看controlDict 里用的pimpleFoam, 于是就试了下，但是给出：
--> FOAM FATAL ERROR:
Residual data for p must be specified as a dictionary
From function virtual bool Foam::solutionControl::read(bool)
in file cfdTools/general/solutionControl/solutionControl/solutionControl.C at line 85.
FOAM exiting
确认下您用的是pimpleFoam求解器吗？感觉这个error应该是版本问题造成的吧，我再去装一个OF7，试试
非常感谢！
