diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/Make/files b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/Make/files index aff39a6ded..14ed51d6ba 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/Make/files +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/Make/files @@ -1,5 +1,6 @@ dragModels/dragModel/dragModel.C dragModels/dragModel/newDragModel.C +dragModels/Beetstra/Beetstra.C dragModels/segregated/segregated.C dragModels/Ergun/Ergun.C dragModels/Gibilaro/Gibilaro.C @@ -10,6 +11,7 @@ dragModels/SchillerNaumann/SchillerNaumann.C dragModels/SyamlalOBrien/SyamlalOBrien.C dragModels/TomiyamaCorrelated/TomiyamaCorrelated.C dragModels/TomiyamaAnalytic/TomiyamaAnalytic.C +dragModels/Tenneti/Tenneti.C dragModels/TomiyamaKataokaZunSakaguchi/TomiyamaKataokaZunSakaguchi.C dragModels/WenYu/WenYu.C dragModels/IshiiZuber/IshiiZuber.C diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Beetstra/Beetstra.C b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Beetstra/Beetstra.C new file mode 100755 index 0000000000..932bff3fb6 --- /dev/null +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Beetstra/Beetstra.C @@ -0,0 +1,101 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | + \\ / A nd | Copyright (C) 2016 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +\*---------------------------------------------------------------------------*/ + +#include "Beetstra.H" +#include "phasePair.H" +#include "addToRunTimeSelectionTable.H" + +// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // + +namespace Foam +{ +namespace dragModels +{ + defineTypeNameAndDebug(Beetstra, 0); + addToRunTimeSelectionTable(dragModel, Beetstra, dictionary); +} +} + + +// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // + +Foam::dragModels::Beetstra::Beetstra +( + const dictionary& dict, + const phasePair& pair, + const bool registerObject +) +: + dragModel(dict, pair, registerObject), + residualRe_("residualRe", dimless, dict.lookup("residualRe")) +{} + + +// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * // + +Foam::dragModels::Beetstra::~Beetstra() +{} + + +// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // + +Foam::tmp Foam::dragModels::Beetstra::CdRe() const +{ + volScalarField alpha1 + ( + max(pair_.dispersed(), pair_.continuous().residualAlpha()) + ); + + volScalarField alpha2 + ( + max(scalar(1) - pair_.dispersed(), pair_.continuous().residualAlpha()) + ); + + volScalarField Re(pair_.Re()); + volScalarField ReLim + ( + "ReLim", + max(Re, residualRe_) + ); + + volScalarField F0 + ( + "F0", + 10.0*alpha1/sqr(alpha2) + sqr(alpha2)*(1.0 + 1.5*sqrt(alpha1)) + ); + + volScalarField F1 + ( + "F1", + 0.413*Re/(24.0*sqr(alpha2))*(1.0/alpha2 + + 3.0*alpha1*alpha2 + 8.4*pow(ReLim, -0.343)) + /(1.0 + pow(10.0, 3*alpha1)*pow(ReLim, -(1.0 + 4.0*alpha1)/2.0)) + ); + + return 24.0*alpha2*(F0 + F1); +} + + +// ************************************************************************* // diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Beetstra/Beetstra.H b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Beetstra/Beetstra.H new file mode 100755 index 0000000000..4a44e944c7 --- /dev/null +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Beetstra/Beetstra.H @@ -0,0 +1,111 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | + \\ / A nd | Copyright (C) 2016 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +Class + Foam::dragModels::Beetstra + +Description + Drag model of Beetstra et al. for monodisperse gas-particle flows obtained + with direct numerical simulations with the Lattice-Boltzmann method and + accounting for the effect of particle ensembles. + + Reference: + \verbatim + Beetstra, R., van der Hoef, M. A., & Kuipers, J. a. M. (2007). + Drag force of intermediate Reynolds number flow past mono- and + bidisperse arrays of spheres. + AIChE Journal, 53(2), 489–501. + \endverbatim + +SourceFiles + Beetstra.C + +\*---------------------------------------------------------------------------*/ + +#ifndef Beetstra_H +#define Beetstra_H + +#include "dragModel.H" + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +namespace Foam +{ + +class phasePair; + +namespace dragModels +{ + +/*---------------------------------------------------------------------------*\ + Class Beetstra Declaration +\*---------------------------------------------------------------------------*/ + +class Beetstra +: + public dragModel +{ + // Private data + + //- Residual Reynolds Number + const dimensionedScalar residualRe_; + + +public: + + //- Runtime type information + TypeName("Beetstra"); + + + // Constructors + + //- Construct from a dictionary and a phase pair + Beetstra + ( + const dictionary& dict, + const phasePair& pair, + const bool registerObject + ); + + + //- Destructor + virtual ~Beetstra(); + + + // Member Functions + + //- Drag coefficient + virtual tmp CdRe() const; +}; + + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +} // End namespace dragModels +} // End namespace Foam + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +#endif + +// ************************************************************************* // diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.C b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.C new file mode 100755 index 0000000000..7fcda489a9 --- /dev/null +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.C @@ -0,0 +1,106 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | + \\ / A nd | Copyright (C) 2016 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +\*---------------------------------------------------------------------------*/ + +#include "Tenneti.H" +#include "phasePair.H" +#include "SchillerNaumann.H" +#include "addToRunTimeSelectionTable.H" + +// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // + +namespace Foam +{ +namespace dragModels +{ + defineTypeNameAndDebug(Tenneti, 0); + addToRunTimeSelectionTable(dragModel, Tenneti, dictionary); +} +} + + +// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // + +Foam::dragModels::Tenneti::Tenneti +( + const dictionary& dict, + const phasePair& pair, + const bool registerObject +) +: + dragModel(dict, pair, registerObject), + SchillerNaumann_ + ( + new SchillerNaumann + ( + dict, + pair, + false + ) + ), + residualRe_("residualRe", dimless, dict.lookup("residualRe")) +{} + + +// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * // + +Foam::dragModels::Tenneti::~Tenneti() +{} + + +// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // + +Foam::tmp Foam::dragModels::Tenneti::CdRe() const +{ + volScalarField alpha1 + ( + max(pair_.dispersed(), pair_.continuous().residualAlpha()) + ); + + volScalarField alpha2 + ( + max(scalar(1) - pair_.dispersed(), pair_.continuous().residualAlpha()) + ); + + volScalarField F0 + ( + 5.81*alpha1/pow3(alpha2) + 0.48*pow(alpha1, 1.0/3.0)/pow4(alpha2) + ); + + volScalarField F1 + ( + pow(alpha1, 3)*max(pair_.Re(), residualRe_) + *(0.95 + 0.61*pow3(alpha1)/sqr(alpha2)) + ); + + // Tenneti et al. correlation includes the mean pressure drag. + // This was removed here by multiplying F by alpha2 for consistency with + // the formulation used in OpenFOAM + return + SchillerNaumann_->CdRe()/(alpha2*max(pair_.Re(), residualRe_)) + + 24.0*sqr(alpha2)*(F0 + F1); +} + + +// ************************************************************************* // diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.H b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.H new file mode 100755 index 0000000000..3003dbc594 --- /dev/null +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.H @@ -0,0 +1,115 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | + \\ / A nd | Copyright (C) 2016 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +Class + Foam::dragModels::Tenneti + +Description + Drag model of Tenneti et al. for monodisperse gas-particle flows obtained + with particle-resolved direct numerical simulations and accounting for the + effect of particle ensembles. + + Reference: + \verbatim + Tenneti, S., Garg, R., & Subramaniam, S. (2011). + Drag law for monodisperse gas–solid systems using particle-resolved + direct numerical simulation of flow past fixed assemblies of spheres. + International Journal of Multiphase Flow, 37(9), 1072–1092. + \verbatim + +SourceFiles + Tenneti.C + +\*---------------------------------------------------------------------------*/ + +#ifndef Tenneti_H +#define Tenneti_H + +#include "dragModel.H" + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +namespace Foam +{ + +class phasePair; + +namespace dragModels +{ + +class SchillerNaumann; +/*---------------------------------------------------------------------------*\ + Class Tenneti Declaration +\*---------------------------------------------------------------------------*/ + +class Tenneti +: + public dragModel +{ + // Private data + + //- Ergun drag model + autoPtr SchillerNaumann_; + + //- Residual Reynolds Number + const dimensionedScalar residualRe_; + + +public: + + //- Runtime type information + TypeName("Tenneti"); + + + // Constructors + + //- Construct from a dictionary and a phase pair + Tenneti + ( + const dictionary& dict, + const phasePair& pair, + const bool registerObject + ); + + + //- Destructor + virtual ~Tenneti(); + + + // Member Functions + + //- Drag coefficient + virtual tmp CdRe() const; +}; + + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +} // End namespace dragModels +} // End namespace Foam + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +#endif + +// ************************************************************************* //