From 4a46fab1ba825ef7b8bc7a887f8b2384c292062e Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 14 Jun 2018 20:09:58 +0100 Subject: [PATCH] *EulerFoam/interfacialModels/dragModels/IshiiZuber: Set constants to correspond to original paper "Drag Coefficient and relative velocity in bubbly, droplet and particulate flows", Ishii, M., Zuber, N., AIChE Journal 5, Vol. 25, 1979, pp. 843-855. --- .../dragModels/IshiiZuber/IshiiZuber.C | 25 +++++++++-------- .../dragModels/IshiiZuber/IshiiZuber.C | 27 +++++++++---------- 2 files changed, 25 insertions(+), 27 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C index bb2eca183..1efa90073 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C @@ -63,33 +63,32 @@ Foam::dragModels::IshiiZuber::~IshiiZuber() Foam::tmp Foam::dragModels::IshiiZuber::CdRe() const { - volScalarField Re(pair_.Re()); - volScalarField Eo(pair_.Eo()); + const volScalarField Re(pair_.Re()); + const volScalarField Eo(pair_.Eo()); - volScalarField mud(pair_.dispersed().mu()); - volScalarField muc(pair_.continuous().mu()); + const volScalarField mud(pair_.dispersed().mu()); + const volScalarField muc(pair_.continuous().mu()); - volScalarField muStar((mud + 0.4*muc)/(mud + muc)); + const volScalarField muStar((mud + 0.4*muc)/(mud + muc)); - volScalarField muMix + const volScalarField muMix ( - muc - *pow(max(1 - pair_.dispersed(), scalar(1e-3)), -2.5*muStar) + muc*pow(max(1 - pair_.dispersed(), scalar(1e-3)), -2.5*muStar) ); - volScalarField ReM(Re*muc/muMix); - volScalarField CdRe + const volScalarField ReM(Re*muc/muMix); + const volScalarField CdRe ( - pos0(1000 - ReM)*24*(scalar(1) + 0.15*pow(ReM, 0.687)) + pos0(1000 - ReM)*24*(scalar(1) + 0.1*pow(ReM, 0.75)) + neg(1000 - ReM)*0.44*ReM ); volScalarField F((muc/muMix)*sqrt(1 - pair_.dispersed())); F.max(1e-3); - volScalarField Ealpha((1 + 17.67*pow(F, 0.8571428))/(18.67*F)); + const volScalarField Ealpha((1 + 17.67*pow(F, 0.8571428))/(18.67*F)); - volScalarField CdReEllipse(Ealpha*0.6666*sqrt(Eo)*Re); + const volScalarField CdReEllipse(Ealpha*0.6666*sqrt(Eo)*Re); return pos0(CdReEllipse - CdRe) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C index 74adfb84a..1efa90073 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/IshiiZuber/IshiiZuber.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2014-2017 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2014-2018 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -63,33 +63,32 @@ Foam::dragModels::IshiiZuber::~IshiiZuber() Foam::tmp Foam::dragModels::IshiiZuber::CdRe() const { - volScalarField Re(pair_.Re()); - volScalarField Eo(pair_.Eo()); + const volScalarField Re(pair_.Re()); + const volScalarField Eo(pair_.Eo()); - volScalarField mud(pair_.dispersed().mu()); - volScalarField muc(pair_.continuous().mu()); + const volScalarField mud(pair_.dispersed().mu()); + const volScalarField muc(pair_.continuous().mu()); - volScalarField muStar((mud + 0.4*muc)/(mud + muc)); + const volScalarField muStar((mud + 0.4*muc)/(mud + muc)); - volScalarField muMix + const volScalarField muMix ( - muc - *pow(max(1 - pair_.dispersed(), scalar(1e-3)), -2.5*muStar) + muc*pow(max(1 - pair_.dispersed(), scalar(1e-3)), -2.5*muStar) ); - volScalarField ReM(Re*muc/muMix); - volScalarField CdRe + const volScalarField ReM(Re*muc/muMix); + const volScalarField CdRe ( - pos0(1000 - ReM)*24.0*(scalar(1) + 0.15*pow(ReM, 0.687)) + pos0(1000 - ReM)*24*(scalar(1) + 0.1*pow(ReM, 0.75)) + neg(1000 - ReM)*0.44*ReM ); volScalarField F((muc/muMix)*sqrt(1 - pair_.dispersed())); F.max(1e-3); - volScalarField Ealpha((1 + 17.67*pow(F, 0.8571428))/(18.67*F)); + const volScalarField Ealpha((1 + 17.67*pow(F, 0.8571428))/(18.67*F)); - volScalarField CdReEllipse(Ealpha*0.6666*sqrt(Eo)*Re); + const volScalarField CdReEllipse(Ealpha*0.6666*sqrt(Eo)*Re); return pos0(CdReEllipse - CdRe)