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 <lagrangianIntermediate/SphereDrag.H> 00027 00028 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // 00029 00030 template <class CloudType> 00031 Foam::SphereDrag<CloudType>::SphereDrag 00032 ( 00033 const dictionary& dict, 00034 CloudType& owner 00035 ) 00036 : 00037 DragModel<CloudType>(dict, owner) 00038 {} 00039 00040 00041 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * // 00042 00043 template <class CloudType> 00044 Foam::SphereDrag<CloudType>::~SphereDrag() 00045 {} 00046 00047 00048 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // 00049 00050 template <class CloudType> 00051 bool Foam::SphereDrag<CloudType>::active() const 00052 { 00053 return true; 00054 } 00055 00056 00057 template <class CloudType> 00058 Foam::scalar Foam::SphereDrag<CloudType>::Cd(const scalar Re) const 00059 { 00060 scalar Cd; 00061 if (Re < SMALL) 00062 { 00063 Cd = GREAT; 00064 } 00065 else if (Re > 1000.0) 00066 { 00067 Cd = 0.424; 00068 } 00069 else 00070 { 00071 Cd = 24.0/Re*(1.0 + 1.0/6.0*pow(Re, 2.0/3.0)); 00072 } 00073 00074 return Cd; 00075 } 00076 00077 00078 // ************************ vim: set sw=4 sts=4 et: ************************ //
