Fix mesh volume handling in mesh flux updates

This commit is contained in:
Hrvoje Jasak 2018-05-16 15:11:12 +01:00
parent 4a16503359
commit 25b3940175

View file

@ -336,6 +336,7 @@ void fvMesh::updatePhi(const scalarField& sweptVols) const
// Make sure V and V0 are constructed before the correction
// HJ, 22/Dec/2017
V();
if (!V0Ptr_)
{
const_cast<fvMesh&>(*this).setV0();
@ -363,7 +364,7 @@ void fvMesh::updatePhi(const scalarField& sweptVols) const
phi.boundaryField()[patchI] *= rDeltaT;
}
// Boundary update. Used in complex geometries, eg. immersed boundary
// Boundary update. Used in complex geometry updates, eg. immersed boundary
// HJ, 29/Nov/2017
forAll (phi.boundaryField(), patchI)
{