src/finiteVolume/fvMesh/wallDist/wallDistData.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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 "wallDistData.H"
  28 #include "patchDataWave.H"
  29 #include "wallPolyPatch.H"
  30 #include "emptyFvPatchFields.H"
  31
  32 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  33
  34 // Construct from components
  35 template<class TransferType>
  36 Foam::wallDistData<TransferType>::wallDistData
  37 (
  38     const Foam::fvMesh& mesh,
  39     GeometricField<Type, fvPatchField, volMesh>& field,
  40     const bool correctWalls
  41 )
  42 :
  43     volScalarField
  44     (
  45         IOobject
  46         (
  47             "y",
  48             mesh.time().timeName(),
  49             mesh
  50         ),
  51         mesh,
  52         dimensionedScalar("y", dimLength, GREAT)
  53     ),
  54     cellDistFuncs(mesh),
  55     field_(field),
  56     correctWalls_(correctWalls),
  57     nUnset_(0)
  58 {
  59     correct();
  60 }
  61
  62
  63 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  64
  65 template<class TransferType>
  66 Foam::wallDistData<TransferType>::~wallDistData()
  67 {}
  68
  69
  70 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  71
  72 // Correct for mesh geom/topo changes
  73 template<class TransferType>
  74 void Foam::wallDistData<TransferType>::correct()
  75 {
  76   Info<< "wallDistData.correct() called" << endl;
  77     const polyMesh& mesh = cellDistFuncs::mesh();
  78
  79     //
  80     // Fill data on wall patches with initial values
  81     //
  82
  83     const polyBoundaryMesh& bMesh = mesh.boundaryMesh();
  84     labelHashSet wallPatchIDs(bMesh.size());
  85     forAll(bMesh, patchI)
  86     {
  87         if (bMesh[patchI].isWall())
  88         {
  89             wallPatchIDs.insert(patchI);
  90         }
  91     }
  92
  93     // Get patchids of walls
  94     // labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
  95
  96     // Collect pointers to data on patches
  97     UPtrList<Field<Type> > patchData(mesh.boundaryMesh().size());
  98
  99     forAll(field_.boundaryField(), patchI)
 100     {
 101         patchData.set(patchI, &field_.boundaryField()[patchI]);
 102     }
 103
 104     // Do mesh wave
 105     patchDataWave<TransferType> wave
 106     (
 107         mesh,
 108         wallPatchIDs,
 109         patchData,
 110         correctWalls_
 111     );
 112
 113     // Transfer cell values from wave into *this and field_
 114     transfer(wave.distance());
 115
 116     field_.transfer(wave.cellData());
 117
 118     // Transfer values on patches into boundaryField of *this and field_
 119     forAll(boundaryField(), patchI)
 120     {
 121         scalarField& waveFld = wave.patchDistance()[patchI];
 122
 123         if (!isA<emptyFvPatchScalarField>(boundaryField()[patchI]))
 124         {
 125             boundaryField()[patchI].transfer(waveFld);
 126
 127             Field<Type>& wavePatchData = wave.patchData()[patchI];
 128
 129             field_.boundaryField()[patchI].transfer(wavePatchData);
 130         }
 131     }
 132
 133     // Transfer number of unset values
 134     nUnset_ = wave.nUnset();
 135 }
 136
 137
 138 // ************************************************************************* //