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

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2011 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 "GidaspowErgunWenYu.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     defineTypeNameAndDebug(GidaspowErgunWenYu, 0);
  34
  35     addToRunTimeSelectionTable
  36     (
  37         dragModel,
  38         GidaspowErgunWenYu,
  39         dictionary
  40     );
  41 }
  42
  43
  44 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  45
  46 Foam::GidaspowErgunWenYu::GidaspowErgunWenYu
  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::GidaspowErgunWenYu::~GidaspowErgunWenYu()
  61 {}
  62
  63
  64 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  65
  66 Foam::tmp<Foam::volScalarField> Foam::GidaspowErgunWenYu::K
  67 (
  68     const volScalarField& Ur
  69 ) const
  70 {
  71     volScalarField beta = max(scalar(1) - alpha_, scalar(1.0e-6));
  72
  73     volScalarField bp = pow(beta, -2.65);
  74     volScalarField Re = max(Ur*phasea_.d()/phaseb_.nu(), scalar(1.0e-3));
  75
  76     volScalarField Cds =
  77         neg(Re - 1000)*(24.0*(1.0 + 0.15*pow(Re, 0.687))/Re)
  78       + pos(Re - 1000)*0.44;
  79
  80     // Wen and Yu (1966)
  81     return
  82         pos(beta - 0.8)
  83        *(0.75*Cds*phaseb_.rho()*Ur*bp/phasea_.d())
  84       + neg(beta - 0.8)
  85        *(
  86            150.0*alpha_*phaseb_.nu()*phaseb_.rho()/(sqr(beta*phasea_.d()))
  87          + 1.75*phaseb_.rho()*Ur/(beta*phasea_.d())
  88         );
  89 }
  90
  91
  92 // ************************************************************************* //