rhoPimpleDyMFoam/correctPhi.H: Add compressibility term in pcorr equation

This commit is contained in:
Henry 2014-05-18 21:19:26 +01:00
parent 08cff2ac41
commit 86498d3dbe

View file

@ -46,14 +46,18 @@ if (mesh.changing())
{
fvScalarMatrix pcorrEqn
(
fvm::laplacian(rAUf, pcorr) == fvc::div(phi) - divrhoU
fvm::ddt(psi, pcorr)
+ fvc::div(phi)
- fvm::laplacian(rAUf, pcorr)
==
divrhoU
);
pcorrEqn.solve();
if (pimple.finalNonOrthogonalIter())
{
phi -= pcorrEqn.flux();
phi += pcorrEqn.flux();
}
}
}