FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

wallDist.C

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 \*---------------------------------------------------------------------------*/
00025 
00026 #include "wallDist.H"
00027 #include <meshTools/patchWave.H>
00028 #include <finiteVolume/fvMesh.H>
00029 #include <OpenFOAM/wallPolyPatch.H>
00030 #include <finiteVolume/fvPatchField.H>
00031 #include <OpenFOAM/Field.H>
00032 #include <finiteVolume/emptyFvPatchFields.H>
00033 
00034 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
00035 
00036 Foam::wallDist::wallDist(const fvMesh& mesh, const bool correctWalls)
00037 :
00038     volScalarField
00039     (
00040         IOobject
00041         (
00042             "y",
00043             mesh.time().timeName(),
00044             mesh
00045         ),
00046         mesh,
00047         dimensionedScalar("y", dimLength, GREAT)
00048     ),
00049     cellDistFuncs(mesh),
00050     correctWalls_(correctWalls),
00051     nUnset_(0)
00052 {
00053     wallDist::correct();
00054 }
00055 
00056 
00057 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
00058 
00059 Foam::wallDist::~wallDist()
00060 {}
00061 
00062 
00063 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
00064 
00065 // Correct for mesh geom/topo changes. Might be more intelligent in the
00066 // future (if only small topology change)
00067 void Foam::wallDist::correct()
00068 {
00069     // Get patchids of walls
00070     labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
00071 
00072     // Calculate distance starting from wallPatch faces.
00073     patchWave wave(cellDistFuncs::mesh(), wallPatchIDs, correctWalls_);
00074 
00075     // Transfer cell values from wave into *this 
00076     transfer(wave.distance());
00077 
00078     // Transfer values on patches into boundaryField of *this
00079     forAll(boundaryField(), patchI)
00080     {
00081         if (!isA<emptyFvPatchScalarField>(boundaryField()[patchI]))
00082         {
00083             scalarField& waveFld = wave.patchDistance()[patchI];
00084 
00085             boundaryField()[patchI].transfer(waveFld);
00086         }
00087     }
00088 
00089     // Transfer number of unset values
00090     nUnset_ = wave.nUnset();
00091 }
00092 
00093 
00094 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines