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 "reflectionVectors.H" 00027 #include <finiteVolume/wallFvPatch.H> 00028 #include <finiteVolume/surfaceFields.H> 00029 00030 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // 00031 00032 // Construct from components 00033 Foam::reflectionVectors::reflectionVectors(const Foam::fvMesh& mesh) 00034 : 00035 n_ 00036 ( 00037 IOobject 00038 ( 00039 "reflectionVectors", 00040 mesh.time().timeName(), 00041 mesh 00042 ), 00043 mesh, 00044 dimensionedVector("n", dimless, vector::zero) 00045 ) 00046 { 00047 correct(); 00048 } 00049 00050 00051 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // 00052 00053 00054 void Foam::reflectionVectors::correct() 00055 { 00056 const fvMesh& mesh = n_.mesh(); 00057 const fvPatchList& patches = mesh.boundary(); 00058 00059 forAll(patches, patchi) 00060 { 00061 // find the nearest face for every cell 00062 if (isA<wallFvPatch>(patches[patchi])) 00063 { 00064 n_.boundaryField()[patchi] = 00065 mesh.Sf().boundaryField()[patchi] 00066 /mesh.magSf().boundaryField()[patchi]; 00067 } 00068 } 00069 } 00070 00071 00072 // ************************ vim: set sw=4 sts=4 et: ************************ //
