BUGFIX: Correct averaging in ProcessorPointPatchField. Author: Vanja Skuric. Merge: Dominik Christ.

This commit is contained in:
Dominik Christ 2015-04-23 12:08:10 +01:00
commit 8f829cc2e7
6 changed files with 13 additions and 26005 deletions

View file

@ -80,3 +80,4 @@ Contents:
Tian Tang Tian Tang
Zeljko Tukovic Zeljko Tukovic
Niklas Wikstrom Niklas Wikstrom
Vanja Skuric

View file

@ -476,7 +476,7 @@ initEvaluate
} }
// Field transfer // Initialise field transfer
template template
< <
template<class> class PatchField, template<class> class PatchField,
@ -511,18 +511,21 @@ evaluate
} }
// Average over two sides // Average over two sides
// VS, 2015-04-12 - This works in parallel!
tpn = 0.5*(this->patchInternalField(this->internalField()) + tpn);
// ZT, 22-07-2014 - point ordering is not same // ZT, 22-07-2014 - point ordering is not same
// at master and slave side after topology change // at master and slave side after topology change
const labelList& neiPoints = const labelList& neiPoints =
procPatch_.procPolyPatch().neighbPoints(); procPatch_.procPolyPatch().neighbPoints();
tpn = // VS, 2015-04-12 - This doesn't work in parallel!
0.5* // tpn =
( // 0.5*
this->patchInternalField(this->internalField()) // (
+ Field<Type>(tpn, neiPoints) // this->patchInternalField(this->internalField())
); // + Field<Type>(tpn, neiPoints)
// );
// Get internal field to insert values into // Get internal field to insert values into
Field<Type>& iF = const_cast<Field<Type>&>(this->internalField()); Field<Type>& iF = const_cast<Field<Type>&>(this->internalField());

View file

@ -5,3 +5,4 @@
cleanCase cleanCase
\rm -rf constant/polyMesh/* \rm -rf constant/polyMesh/*
\rm -f naca0012.cas

View file

@ -7,4 +7,5 @@ application="dbnsTurbFoam"
gunzip naca0012.cas.gz gunzip naca0012.cas.gz
runApplication fluentMeshToFoam naca0012.cas runApplication fluentMeshToFoam naca0012.cas
gzip naca0012.cas
runApplication $application runApplication $application

File diff suppressed because one or more lines are too long