src/finiteVolume/fvMesh/fvPatches/constraint/cyclic/cyclicFvPatch.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "cyclicFvPatch.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvMesh.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  37
  38 defineTypeNameAndDebug(cyclicFvPatch, 0);
  39 addToRunTimeSelectionTable(fvPatch, cyclicFvPatch, polyPatch);
  40
  41
  42 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  43
  44 // Make patch weighting factors
  45 void cyclicFvPatch::makeWeights(scalarField& w) const
  46 {
  47     const scalarField& magFa = magSf();
  48
  49     scalarField deltas = nf() & fvPatch::delta();
  50     label sizeby2 = deltas.size()/2;
  51
  52     for (label facei = 0; facei < sizeby2; facei++)
  53     {
  54         scalar avFa = (magFa[facei] + magFa[facei + sizeby2])/2.0;
  55
  56         if (mag(magFa[facei] - magFa[facei + sizeby2])/avFa > 1e-4)
  57         {
  58             FatalErrorIn("cyclicFvPatch::makeWeights(scalarField& w) const")
  59                 << "face " << facei << " and " << facei + sizeby2
  60                 <<  " areas do not match by "
  61                 << 100*mag(magFa[facei] - magFa[facei + sizeby2])/avFa
  62                 << "% -- possible face ordering problem"
  63                 << abort(FatalError);
  64         }
  65
  66         scalar di = deltas[facei];
  67         scalar dni = deltas[facei + sizeby2];
  68
  69         w[facei] = dni/(di + dni);
  70         w[facei + sizeby2] = 1 - w[facei];
  71     }
  72 }
  73
  74
  75 // Make patch face - neighbour cell distances
  76 void cyclicFvPatch::makeDeltaCoeffs(scalarField& dc) const
  77 {
  78     scalarField deltas = nf() & fvPatch::delta();
  79     label sizeby2 = deltas.size()/2;
  80
  81     for (label facei = 0; facei < sizeby2; facei++)
  82     {
  83         scalar di = deltas[facei];
  84         scalar dni = deltas[facei + sizeby2];
  85
  86         dc[facei] = 1.0/(di + dni);
  87         dc[facei + sizeby2] = dc[facei];
  88     }
  89 }
  90
  91
  92 // Return delta (P to N) vectors across coupled patch
  93 tmp<vectorField> cyclicFvPatch::delta() const
  94 {
  95     vectorField patchD = fvPatch::delta();
  96     label sizeby2 = patchD.size()/2;
  97
  98     tmp<vectorField> tpdv(new vectorField(patchD.size()));
  99     vectorField& pdv = tpdv();
 100
 101     // To the transformation if necessary
 102     if (parallel())
 103     {
 104         for (label facei = 0; facei < sizeby2; facei++)
 105         {
 106             vector ddi = patchD[facei];
 107             vector dni = patchD[facei + sizeby2];
 108
 109             pdv[facei] = ddi - dni;
 110             pdv[facei + sizeby2] = -pdv[facei];
 111         }
 112     }
 113     else
 114     {
 115         for (label facei = 0; facei < sizeby2; facei++)
 116         {
 117             vector ddi = patchD[facei];
 118             vector dni = patchD[facei + sizeby2];
 119
 120             pdv[facei] = ddi - transform(forwardT()[0], dni);
 121             pdv[facei + sizeby2] = -transform(reverseT()[0], pdv[facei]);
 122         }
 123     }
 124
 125     return tpdv;
 126 }
 127
 128
 129 tmp<labelField> cyclicFvPatch::interfaceInternalField
 130 (
 131     const unallocLabelList& internalData
 132 ) const
 133 {
 134     return patchInternalField(internalData);
 135 }
 136
 137
 138 tmp<labelField> cyclicFvPatch::transfer
 139 (
 140     const Pstream::commsTypes,
 141     const unallocLabelList& interfaceData
 142 ) const
 143 {
 144     tmp<labelField> tpnf(new labelField(this->size()));
 145     labelField& pnf = tpnf();
 146
 147     label sizeby2 = this->size()/2;
 148
 149     for (label facei=0; facei<sizeby2; facei++)
 150     {
 151         pnf[facei] = interfaceData[facei + sizeby2];
 152         pnf[facei + sizeby2] = interfaceData[facei];
 153     }
 154
 155     return tpnf;
 156 }
 157
 158
 159 tmp<labelField> cyclicFvPatch::internalFieldTransfer
 160 (
 161     const Pstream::commsTypes commsType,
 162     const unallocLabelList& iF
 163 ) const
 164 {
 165     const unallocLabelList& faceCells = this->patch().faceCells();
 166
 167     tmp<labelField> tpnf(new labelField(this->size()));
 168     labelField& pnf = tpnf();
 169
 170     label sizeby2 = this->size()/2;
 171
 172     for (label facei=0; facei<sizeby2; facei++)
 173     {
 174         pnf[facei] = iF[faceCells[facei + sizeby2]];
 175         pnf[facei + sizeby2] = iF[faceCells[facei]];
 176     }
 177
 178     return tpnf;
 179 }
 180
 181
 182 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 183
 184 } // End namespace Foam
 185
 186 // ************************************************************************* //