Some questions about the functions ddtPhiCorr(rUA, U, phi) and adjustPhi(phi, U, p)
I am reading the icoFoam codes, where I find two functions ddtPhiCorr(rUA, U, phi) and adjustPhi(phi, U, p):
volScalarField rUA = 1.0/UEqn.A();
U = rUA*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf())
+ fvc::ddtPhiCorr(rUA, U, phi);
adjustPhi(phi, U, p);
Could you tell me why the two functions are used here? And where could I find the explanation of these two functions in official documentation? Thank you very much!
I'm afraid there aren't any official documents regrading this topic, but you can refer to the source code.
1) adjustPhi(phi, U, p)
2)ddtPhiCorr(rUA, U, phi)
You might want to read the following thread
Why is adjustPhi taken away in, e.g. compressibleInterFoam? Shouldn't the flux still be adjusted to satisfy continuity at the boundary?
adjustPhi is taken away from compressible solvers because you are no longer solving Poisson equation for getting the pressure. look at post #4 in http://www.cfd-online.com/Forums/ope...adjustphi.html
so there are no singularities to be dealt with. for example lets look at pEqn in laminarSMOKE :
- fvm::laplacian(rho*rUA, p)
|All times are GMT -4. The time now is 21:31.|