zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: jump should be negated on slave patch.

Browse Source
This commit is contained in:
mattijs
2010-04-16 21:52:39 +01:00
parent 0712a68639
commit 55d6460445
1 changed files with 8 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/finiteVolume/fields/fvPatchFields/constraint/jumpCyclic/jumpCyclicFvPatchField.C
View File

@ -106,6 +106,10 @@ tmp<Field<Type> > jumpCyclicFvPatchField<Type>::patchNeighbourField() const
Field<Type>& pnf = tpnf(); Field<Type>& pnf = tpnf();
tmp<Field<scalar> > tjf = jump(); tmp<Field<scalar> > tjf = jump();
if (!this->cyclicPatch().owner())
{
tjf = -tjf;
}
const Field<scalar>& jf = tjf(); const Field<scalar>& jf = tjf();
if (this->doTransform()) if (this->doTransform())
@ -149,6 +153,10 @@ void jumpCyclicFvPatchField<Type>::updateInterfaceMatrix
if (&psiInternal == &this->internalField()) if (&psiInternal == &this->internalField())
{ {
tmp<Field<scalar> > tjf = jump(); tmp<Field<scalar> > tjf = jump();
if (!this->cyclicPatch().owner())
{
tjf = -tjf;
}
const Field<scalar>& jf = tjf(); const Field<scalar>& jf = tjf();
forAll(*this, facei) forAll(*this, facei)