src/lagrangian/intermediate/submodels/Kinematic/DragModel/SphereDrag/SphereDrag.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "SphereDrag.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 template <class CloudType>
  32 Foam::SphereDrag<CloudType>::SphereDrag
  33 (
  34     const dictionary& dict,
  35     CloudType& owner
  36 )
  37 :
  38     DragModel<CloudType>(dict, owner)
  39 {}
  40
  41
  42 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  43
  44 template <class CloudType>
  45 Foam::SphereDrag<CloudType>::~SphereDrag()
  46 {}
  47
  48
  49 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  50
  51 template <class CloudType>
  52 bool Foam::SphereDrag<CloudType>::active() const
  53 {
  54     return true;
  55 }
  56
  57
  58 template <class CloudType>
  59 Foam::scalar Foam::SphereDrag<CloudType>::Cd(const scalar Re) const
  60 {
  61     scalar Cd;
  62     if (Re < SMALL)
  63     {
  64         Cd = GREAT;
  65     }
  66     else if (Re > 1000.0)
  67     {
  68         Cd =  0.424;
  69     }
  70     else
  71     {
  72         Cd = 24.0/Re*(1.0 + 1.0/6.0*pow(Re, 2.0/3.0));
  73     }
  74
  75     return Cd;
  76 }
  77
  78
  79 // ************************************************************************* //