exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bridging overlap bugfix in ggiPolyPatch.C

Browse source
This commit is contained in:
Vuko Vukcevic 2018-03-02 13:14:36 +01:00
parent 28e6b2de20
commit 9a3b72961c
1 changed files with 6 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
src/foam/meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C
View file

@ -294,17 +294,16 @@ void Foam::ggiPolyPatch::calcReconFaceCellCentres() const
- boundaryMesh()[shadowID].faceCentres()
);
// Get face centres on master side
const vectorField::subField cf = faceCentres();
if (bridgeOverlap_)
{
// Get necessary mesh data from polyPatch on this (master) side
const vectorField::subField cf = faceCentres();
const vectorField::subField Sf = faceAreas();
// Get face cell centres on master side
const vectorField ccf = faceCellCentres();
const vectorField nf = Sf/mag(Sf);
// Deltas for fully uncovered faces
const vectorField uncoveredDeltas(2.0*(cf - ccf));
const vectorField uncoveredDeltas(cf - ccf);
// Scale partially overlapping faces and set uncovered deltas to
// fully uncovered faces
@ -312,7 +311,7 @@ void Foam::ggiPolyPatch::calcReconFaceCellCentres() const
}
// Calculate the reconstructed cell centres
reconFaceCellCentresPtr_ = new vectorField(tdf() + faceCentres());
reconFaceCellCentresPtr_ = new vectorField(tdf() + cf);
}
else
{