src/finiteVolume/fvMesh/wallDist/nearWallDist.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 "nearWallDist.H"
  28 #include "fvMesh.H"
  29 #include "cellDistFuncs.H"
  30 #include "wallFvPatch.H"
  31 #include "surfaceFields.H"
  32
  33 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  34
  35 void Foam::nearWallDist::doAll()
  36 {
  37     cellDistFuncs wallUtils(mesh_);
  38
  39     // Get patch ids of walls
  40     labelHashSet wallPatchIDs
  41     (
  42         wallUtils.getPatchIDs
  43         (
  44             wallPolyPatch::typeName
  45         )
  46     );
  47
  48     // Size neighbours array for maximum possible
  49
  50     labelList neighbours(wallUtils.maxPatchSize(wallPatchIDs));
  51
  52
  53     // Correct all cells with face on wall
  54
  55     const volVectorField& cellCentres = mesh_.C();
  56
  57     forAll(mesh_.boundary(), patchI)
  58     {
  59         fvPatchScalarField& ypatch = operator[](patchI);
  60
  61         const fvPatch& patch = mesh_.boundary()[patchI];
  62
  63         if (patch.type() == wallFvPatch::typeName)
  64         {
  65             const polyPatch& pPatch = patch.patch();
  66
  67             const unallocLabelList& faceCells = patch.faceCells();
  68
  69             // Check cells with face on wall
  70             forAll(patch, patchFaceI)
  71             {
  72                 label nNeighbours = wallUtils.getPointNeighbours
  73                 (
  74                     pPatch,
  75                     patchFaceI,
  76                     neighbours
  77                 );
  78
  79                 label minFaceI = -1;
  80
  81                 ypatch[patchFaceI] = wallUtils.smallestDist
  82                 (
  83                     cellCentres[faceCells[patchFaceI]],
  84                     pPatch,
  85                     nNeighbours,
  86                     neighbours,
  87                     minFaceI
  88                 );
  89             }
  90         }
  91         else
  92         {
  93             ypatch = 0.0;
  94         }
  95     }
  96 }
  97
  98
  99 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 100
 101 Foam::nearWallDist::nearWallDist(const Foam::fvMesh& mesh)
 102 :
 103     volScalarField::GeometricBoundaryField
 104     (
 105         mesh.boundary(),
 106         mesh.V(),           // Dummy internal field,
 107         calculatedFvPatchScalarField::typeName
 108     ),
 109     mesh_(mesh)
 110 {
 111     doAll();
 112 }
 113
 114
 115 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 116
 117 Foam::nearWallDist::~nearWallDist()
 118 {}
 119
 120
 121 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 122
 123 void Foam::nearWallDist::correct()
 124 {
 125     if (mesh_.changing())
 126     {
 127         // Update size of GeometricBoundaryField
 128         forAll(mesh_.boundary(), patchI)
 129         {
 130             operator[](patchI).setSize(mesh_.boundary()[patchI].size());
 131         }
 132     }
 133
 134     doAll();
 135 }
 136
 137
 138 // ************************************************************************* //