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

Bugfix: volPointInterpolation across coupled boundaries. Vuko Vukcevic

Browse source
This commit is contained in:
Hrvoje Jasak 2016-06-22 17:06:40 +01:00
parent bab8ef6814
commit 0c113bfaf5
1 changed files with 23 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

24
src/finiteVolume/interpolation/volPointInterpolation/pointPatchInterpolation/pointPatchInterpolate.C
View file

@ -64,7 +64,7 @@ void pointPatchInterpolation::interpolate
forAll(bm, patchi)
{
if (!isA<emptyFvPatch>(bm[patchi]) && bm[patchi].coupled())
if (!isA<emptyFvPatch>(bm[patchi]) && !bm[patchi].coupled())
{
pointPatchField<Type>& ppf = pf.boundaryField()[patchi];
@ -110,6 +110,28 @@ void pointPatchInterpolation::interpolate
}
}
// Add missing contributions across the coupled boundaries. The weights are
// constructed and normalised by correctly taking into account their
// contributions from the other side, so we need to add the field
// contributions from the other side. VV, 21/June/2016.
forAll(pf.boundaryField(), patchi)
{
if (pf.boundaryField()[patchi].coupled())
{
pf.boundaryField()[patchi].initAddField();
}
}
forAll(pf.boundaryField(), patchi)
{
if (pf.boundaryField()[patchi].coupled())
{
pf.boundaryField()[patchi].addField
(
pf.internalField()
);
}
}
// Correct patch-patch boundary points by interpolation "around" corners
const labelListList& PointFaces = fvMesh_.pointFaces();