Handle bridge overlap in region couple

This commit is contained in:
Hrvoje Jasak 2012-01-29 11:45:46 +00:00
parent 7dd740bf24
commit b5d1c30299

View file

@ -353,6 +353,19 @@ void regionCouplingFvPatchField<Type>::initEvaluate
// Do interpolation
Field<Type>::operator=(weights*fOwn + (1.0 - weights)*fNei);
if (regionCouplePatch_.bridgeOverlap())
{
// Symmetry treatment used for overlap
vectorField nHat = this->patch().nf();
Field<Type> pif = this->patchInternalField();
Field<Type> bridgeField =
0.5*(pif + transform(I - 2.0*sqr(nHat), pif));
regionCouplePatch_.bridge(bridgeField, *this);
}
}