MULES: correct the handling of boundary conditions for limiting the sum

This commit is contained in:
Henry
2011-09-05 19:12:42 +01:00
parent 8cd7ca50a3
commit 94e86ddd97

View File

@ -670,17 +670,20 @@ void Foam::MULES::limitSum(SurfaceScalarFieldList& phiPsiCorrs)
forAll(phiPsiCorrs[0].boundaryField(), patchi)
{
UPtrList<scalarField> phiPsiCorrsPatch(phiPsiCorrs.size());
forAll(phiPsiCorrs, phasei)
if (phiPsiCorrs[0].boundaryField()[patchi].coupled())
{
phiPsiCorrsPatch.set
(
phasei,
&phiPsiCorrs[phasei].boundaryField()[patchi]
);
}
UPtrList<scalarField> phiPsiCorrsPatch(phiPsiCorrs.size());
forAll(phiPsiCorrs, phasei)
{
phiPsiCorrsPatch.set
(
phasei,
&phiPsiCorrs[phasei].boundaryField()[patchi]
);
}
limitSum(phiPsiCorrsPatch);
limitSum(phiPsiCorrsPatch);
}
}
}