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: ************************ //