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 <finiteVolume/volFields.H> 00027 00028 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00029 00030 namespace Foam 00031 { 00032 00033 namespace limitFuncs 00034 { 00035 00036 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00037 00038 template<class Type> 00039 inline tmp<volScalarField> magSqr<Type>::operator() 00040 ( 00041 const GeometricField<Type, fvPatchField, volMesh>& phi 00042 ) const 00043 { 00044 return Foam::magSqr(phi); 00045 } 00046 00047 00048 template<> 00049 inline tmp<volScalarField> magSqr<scalar>::operator() 00050 ( 00051 const volScalarField& phi 00052 ) const 00053 { 00054 return phi; 00055 } 00056 00057 00058 template<> 00059 inline tmp<volScalarField> magSqr<tensor>::operator() 00060 ( 00061 const volTensorField& phi 00062 ) const 00063 { 00064 return Foam::tr(phi); 00065 } 00066 00067 00068 template<class Type> 00069 inline tmp<volScalarField> rhoMagSqr<Type>::operator() 00070 ( 00071 const GeometricField<Type, fvPatchField, volMesh>& phi 00072 ) const 00073 { 00074 const volScalarField& rho = 00075 phi.db().objectRegistry::lookupObject<volScalarField>("rho"); 00076 return Foam::magSqr(phi/rho); 00077 } 00078 00079 00080 template<> 00081 inline tmp<volScalarField> rhoMagSqr<scalar>::operator() 00082 ( 00083 const volScalarField& phi 00084 ) const 00085 { 00086 const volScalarField& rho = 00087 phi.db().objectRegistry::lookupObject<volScalarField>("rho"); 00088 return phi/rho; 00089 } 00090 00091 00092 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00093 00094 } // End namespace limitFuncs 00095 } // End namespace Foam 00096 00097 // ************************ vim: set sw=4 sts=4 et: ************************ //
