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

Fixed GGI transform

Browse source
This commit is contained in:
Ivor Clifford 2011-10-24 13:26:58 -06:00
parent 58e2ce9500
commit d0a6150dc3
6 changed files with 5 additions and 34 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

22
src/OpenFOAM/matrices/lduMatrix/lduAddressing/lduInterfaceFields/ggiLduInterfaceField/ggiLduInterfaceField.H
View file

@ -95,24 +95,8 @@ public:
scalarField& psiInternal, scalarField& psiInternal,
const direction cmpt const direction cmpt
) const; ) const;
//- Transform given patch internal field
template<class Type>
void transformCoupleField
(
Field<Type>& psiInternal
) const;
}; };
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//- Transform given patch internal field
template<>
void ggiLduInterfaceField::transformCoupleField
(
scalarField& psiInternal
) const;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -120,12 +104,6 @@ void ggiLduInterfaceField::transformCoupleField
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#ifdef NoRepository
# include "ggiLduInterfaceFieldTemplates.C"
#endif
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif #endif
// ************************************************************************* // // ************************************************************************* //

2
src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/cyclicGgi/cyclicGgiPolyPatch.C
View file

@ -241,7 +241,7 @@ const Foam::cyclicGgiPolyPatch& Foam::cyclicGgiPolyPatch::cyclicShadow() const
} }
void Foam::cyclicGgiPolyPatch::calcTransforms() void Foam::cyclicGgiPolyPatch::calcTransforms() const
{ {
if (active() && debug) if (active() && debug)
{ {

2
src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/cyclicGgi/cyclicGgiPolyPatch.H
View file

@ -75,7 +75,7 @@ class cyclicGgiPolyPatch
//- Calculate cyclic transforms (rotation and translation) //- Calculate cyclic transforms (rotation and translation)
// Virtual over-ride for base GGI patch. HJ, 14/Jan/2009 // Virtual over-ride for base GGI patch. HJ, 14/Jan/2009
virtual void calcTransforms(); virtual void calcTransforms() const;
//- Check definition: angles and offsets //- Check definition: angles and offsets
void checkDefinition() const; void checkDefinition() const;

2
src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C
View file

@ -859,7 +859,7 @@ void Foam::ggiPolyPatch::updateMesh()
} }
void Foam::ggiPolyPatch::calcTransforms() void Foam::ggiPolyPatch::calcTransforms() const
{ {
// Simplest GGI: no transform or separation. HJ, 24/Oct/2008 // Simplest GGI: no transform or separation. HJ, 24/Oct/2008
forwardT_.setSize(0); forwardT_.setSize(0);

2
src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.H
View file

@ -127,7 +127,7 @@ class ggiPolyPatch
void calcReconFaceCellCentres() const; void calcReconFaceCellCentres() const;
//- Force calculation of transformation tensors //- Force calculation of transformation tensors
virtual void calcTransforms(); virtual void calcTransforms() const;
// Parallel communication optimisation, stored on master processor // Parallel communication optimisation, stored on master processor

9
src/finiteVolume/fields/fvPatchFields/constraint/cyclicGgi/cyclicGgiFvPatchFields.C
View file

@ -106,14 +106,7 @@ void cyclicGgiFvPatchField<vector>::initInterfaceMatrixUpdate
sField[i] = psiInternal[sfc[i]]; sField[i] = psiInternal[sfc[i]];
} }
// Transform according to the transformation tensor, using the slave // Transformation is handled in interpolation. HJ, 7/Jan/2009
// side transform. Warning: forwardT() corresponds to the slave
// patch. HJ, 12/Jan/2009
transformCoupleField(sField);
// Note: scalar interpolate does not get a transform, so this is safe
// HJ, 12/Jan/2009
// TODO - IC, Is this safe for vectors?
Field<vector> pnf = cyclicGgiPatch_.interpolate(sField); Field<vector> pnf = cyclicGgiPatch_.interpolate(sField);
if (coeffs.activeType() == blockCoeffBase::SCALAR) if (coeffs.activeType() == blockCoeffBase::SCALAR)