applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2010 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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "Ergun.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     defineTypeNameAndDebug(Ergun, 0);
  34
  35     addToRunTimeSelectionTable
  36     (
  37         dragModel,
  38         Ergun,
  39         dictionary
  40     );
  41 }
  42
  43
  44 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  45
  46 Foam::Ergun::Ergun
  47 (
  48     const dictionary& interfaceDict,
  49     const volScalarField& alpha,
  50     const phaseModel& phasea,
  51     const phaseModel& phaseb
  52 )
  53 :
  54     dragModel(interfaceDict, alpha, phasea, phaseb)
  55 {}
  56
  57
  58 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  59
  60 Foam::Ergun::~Ergun()
  61 {}
  62
  63
  64 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  65
  66 Foam::tmp<Foam::volScalarField> Foam::Ergun::K
  67 (
  68     const volScalarField& Ur
  69 ) const
  70 {
  71     volScalarField beta(max(scalar(1) - alpha_, scalar(1.0e-6)));
  72
  73     return
  74         150.0*alpha_*phaseb_.nu()*phaseb_.rho()
  75        /sqr(beta*phasea_.d())
  76       + 1.75*phaseb_.rho()*Ur/(beta*phasea_.d());
  77 }
  78
  79
  80 // ************************************************************************* //