compressibleInterFilmFoam::VoFPatchTransfer: corrected parallel handling of all patch option

This commit is contained in:
Henry Weller
2017-06-28 15:03:40 +01:00
parent 0c24327b32
commit 2d866e281b
5 changed files with 25 additions and 11 deletions

View File

@ -1,4 +1,4 @@
EXE_INC = -ggdb3 \ EXE_INC = \
-I. \ -I. \
-I.. \ -I.. \
-I../../VoF \ -I../../VoF \

View File

@ -72,7 +72,10 @@ VoFPatchTransfer::VoFPatchTransfer
) )
{ {
const polyBoundaryMesh& pbm = film.regionMesh().boundaryMesh(); const polyBoundaryMesh& pbm = film.regionMesh().boundaryMesh();
patchIDs_.setSize(pbm.size()); patchIDs_.setSize
(
pbm.size() - film.regionMesh().globalData().processorPatches().size()
);
if (coeffDict_.found("patches")) if (coeffDict_.found("patches"))
{ {
@ -100,7 +103,7 @@ VoFPatchTransfer::VoFPatchTransfer
patchIDs_[patchi] = patchi; patchIDs_[patchi] = patchi;
} }
patchTransferredMasses_.setSize(pbm.size(), 0); patchTransferredMasses_.setSize(patchIDs_.size(), 0);
} }
if (!patchIDs_.size()) if (!patchIDs_.size())

View File

@ -130,7 +130,11 @@ void injectionModelList::info(Ostream& os)
const polyBoundaryMesh& pbm = film().regionMesh().boundaryMesh(); const polyBoundaryMesh& pbm = film().regionMesh().boundaryMesh();
scalar injectedMass = 0; scalar injectedMass = 0;
scalarField patchInjectedMasses(pbm.size(), 0); scalarField patchInjectedMasses
(
pbm.size() - film().regionMesh().globalData().processorPatches().size(),
0
);
forAll(*this, i) forAll(*this, i)
{ {
@ -141,7 +145,7 @@ void injectionModelList::info(Ostream& os)
os << indent << "injected mass = " << injectedMass << nl; os << indent << "injected mass = " << injectedMass << nl;
forAll(pbm, patchi) forAll(patchInjectedMasses, patchi)
{ {
if (mag(patchInjectedMasses[patchi]) > VSMALL) if (mag(patchInjectedMasses[patchi]) > VSMALL)
{ {
@ -150,7 +154,7 @@ void injectionModelList::info(Ostream& os)
} }
} }
scalarField mass0(massInjected_.size(), 0.0); scalarField mass0(massInjected_.size(), 0);
this->getBaseProperty("massInjected", mass0); this->getBaseProperty("massInjected", mass0);
scalarField mass(massInjected_); scalarField mass(massInjected_);

View File

@ -52,7 +52,10 @@ patchInjection::patchInjection
deltaStable_(coeffDict_.lookupOrDefault<scalar>("deltaStable", 0.0)) deltaStable_(coeffDict_.lookupOrDefault<scalar>("deltaStable", 0.0))
{ {
const polyBoundaryMesh& pbm = film.regionMesh().boundaryMesh(); const polyBoundaryMesh& pbm = film.regionMesh().boundaryMesh();
patchIDs_.setSize(pbm.size()); patchIDs_.setSize
(
pbm.size() - film.regionMesh().globalData().processorPatches().size()
);
if (coeffDict_.found("patches")) if (coeffDict_.found("patches"))
{ {
@ -80,7 +83,7 @@ patchInjection::patchInjection
patchIDs_[patchi] = patchi; patchIDs_[patchi] = patchi;
} }
patchInjectedMasses_.setSize(pbm.size(), 0); patchInjectedMasses_.setSize(patchIDs_.size(), 0);
} }
if (!patchIDs_.size()) if (!patchIDs_.size())

View File

@ -158,7 +158,11 @@ void transferModelList::info(Ostream& os)
const polyBoundaryMesh& pbm = film().regionMesh().boundaryMesh(); const polyBoundaryMesh& pbm = film().regionMesh().boundaryMesh();
scalar transferredMass = 0; scalar transferredMass = 0;
scalarField patchTransferredMasses(pbm.size(), 0); scalarField patchTransferredMasses
(
pbm.size() - film().regionMesh().globalData().processorPatches().size(),
0
);
forAll(*this, i) forAll(*this, i)
{ {
@ -169,7 +173,7 @@ void transferModelList::info(Ostream& os)
os << indent << "transferred mass = " << transferredMass << nl; os << indent << "transferred mass = " << transferredMass << nl;
forAll(pbm, patchi) forAll(patchTransferredMasses, patchi)
{ {
if (mag(patchTransferredMasses[patchi]) > VSMALL) if (mag(patchTransferredMasses[patchi]) > VSMALL)
{ {
@ -178,7 +182,7 @@ void transferModelList::info(Ostream& os)
} }
} }
scalarField mass0(massTransferred_.size(), 0.0); scalarField mass0(massTransferred_.size(), 0);
this->getBaseProperty("massTransferred", mass0); this->getBaseProperty("massTransferred", mass0);
scalarField mass(massTransferred_); scalarField mass(massTransferred_);