![]() |
|
00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 1991-2005 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 the 00013 Free Software Foundation; either version 2 of the License, or (at your 00014 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, write to the Free Software Foundation, 00023 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA 00024 00025 Class 00026 midPoint 00027 00028 Description 00029 Mid-point interpolation (weighting factors = 0.5) scheme class. 00030 00031 SourceFiles 00032 midPoint.C 00033 00034 \*---------------------------------------------------------------------------*/ 00035 00036 #ifndef midPoint_H 00037 #define midPoint_H 00038 00039 #include "surfaceInterpolationScheme.H" 00040 #include "volFields.H" 00041 00042 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00043 00044 namespace Foam 00045 { 00046 00047 /*---------------------------------------------------------------------------*\ 00048 Class midPoint Declaration 00049 \*---------------------------------------------------------------------------*/ 00050 00051 template<class Type> 00052 class midPoint 00053 : 00054 public surfaceInterpolationScheme<Type> 00055 { 00056 // Private Member Functions 00057 00058 //- Disallow default bitwise assignment 00059 void operator=(const midPoint&); 00060 00061 00062 public: 00063 00064 //- Runtime type information 00065 TypeName("midPoint"); 00066 00067 00068 // Constructors 00069 00070 //- Construct from mesh 00071 midPoint(const fvMesh& mesh) 00072 : 00073 surfaceInterpolationScheme<Type>(mesh) 00074 {} 00075 00076 //- Construct from Istream 00077 midPoint(const fvMesh& mesh, Istream&) 00078 : 00079 surfaceInterpolationScheme<Type>(mesh) 00080 {} 00081 00082 //- Construct from faceFlux and Istream 00083 midPoint 00084 ( 00085 const fvMesh& mesh, 00086 const surfaceScalarField&, 00087 Istream& 00088 ) 00089 : 00090 surfaceInterpolationScheme<Type>(mesh) 00091 {} 00092 00093 00094 // Member Functions 00095 00096 //- Return the interpolation weighting factors 00097 tmp<surfaceScalarField> weights 00098 ( 00099 const GeometricField<Type, fvPatchField, volMesh>& 00100 ) const 00101 { 00102 tmp<surfaceScalarField> taw 00103 ( 00104 new surfaceScalarField 00105 ( 00106 IOobject 00107 ( 00108 "midPointWeights", 00109 this->mesh().time().timeName(), 00110 this->mesh() 00111 ), 00112 this->mesh(), 00113 dimensionedScalar("0.5", dimless, 0.5) 00114 ) 00115 ); 00116 00117 surfaceScalarField::GeometricBoundaryField& awbf = 00118 taw().boundaryField(); 00119 00120 forAll(awbf, patchi) 00121 { 00122 if (!awbf[patchi].coupled()) 00123 { 00124 awbf[patchi] = 1.0; 00125 } 00126 } 00127 00128 return taw; 00129 } 00130 }; 00131 00132 00133 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00134 00135 } // End namespace Foam 00136 00137 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00138 00139 #endif 00140 00141 // ************************************************************************* //