zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

multiphaseEuler/interfacialModels: Reverted handling of the continuous phase-fraction

Browse Source 
back to the form in the original multiphaseEulerFoam, i.e. using the continuous
phase-fraction directly rather than 1 - dispersed phase fraction.
This commit is contained in:
Henry Weller
2023-04-04 14:59:20 +01:00
parent b414884142
commit d33eafe7df
7 changed files with 19 additions and 28 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/modules/multiphaseEuler/interfacialModels/dragModels/Beetstra/Beetstra.C
View File

@ -63,12 +63,12 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::Beetstra::CdRe() const
{
const volScalarField alpha1
(
max(interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.dispersed(), interface_.dispersed().residualAlpha())
);
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
const volScalarField Res(max(alpha2*interface_.Re(), scalar(0)));

7
applications/solvers/modules/multiphaseEuler/interfacialModels/dragModels/Gibilaro/Gibilaro.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2023 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -63,13 +63,12 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::Gibilaro::CdRe() const
{
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
return
(4.0/3.0)
*(17.3/alpha2 + 0.336*interface_.Re())
*max(interface_.continuous(), interface_.continuous().residualAlpha())
*(17.3/alpha2 + 0.336*interface_.Re())*alpha2
*pow(alpha2, -2.8);
}

15
applications/solvers/modules/multiphaseEuler/interfacialModels/dragModels/GidaspowSchillerNaumann/GidaspowSchillerNaumann.C
View File

@ -64,21 +64,18 @@ Foam::dragModels::GidaspowSchillerNaumann::CdRe() const
{
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
const volScalarField Re(alpha2*interface_.Re());
const volScalarField Res(alpha2*interface_.Re());
const volScalarField CdsRe
const volScalarField CdsRes
(
neg(Re - 1000)*24*(1.0 + 0.15*pow(Re, 0.687))/alpha2
+ pos0(Re - 1000)*0.44*Re
neg(Res - 1000)*24*(1.0 + 0.15*pow(Res, 0.687))/alpha2
+ pos0(Res - 1000)*0.44*Res
);
return
CdsRe
*pow(alpha2, -2.65)
*max(interface_.continuous(), interface_.continuous().residualAlpha());
return CdsRes*pow(alpha2, -1.65);
}

9
applications/solvers/modules/multiphaseEuler/interfacialModels/dragModels/SyamlalOBrien/SyamlalOBrien.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2023 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -63,7 +63,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::SyamlalOBrien::CdRe() const
{
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
const volScalarField A(pow(alpha2, 4.14));
@ -82,10 +82,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::SyamlalOBrien::CdRe() const
);
volScalarField CdsRe(sqr(0.63*sqrt(Re) + 4.8*sqrt(Vr)));
return
CdsRe
*max(interface_.continuous(), interface_.continuous().residualAlpha())
/sqr(Vr);
return CdsRe*alpha2/sqr(Vr);
}

2
applications/solvers/modules/multiphaseEuler/interfacialModels/dragModels/Tenneti/Tenneti.C
View File

@ -69,7 +69,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::Tenneti::CdRe() const
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
const volScalarField Res(alpha2*interface_.Re());

8
applications/solvers/modules/multiphaseEuler/interfacialModels/dragModels/WenYu/WenYu.C
View File

@ -63,20 +63,18 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::WenYu::CdRe() const
{
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
const volScalarField Res(alpha2*interface_.Re());
const volScalarField CdsRes
(
neg(Res - 1000)*24*(1.0 + 0.15*pow(Res, 0.687))
+ pos0(Res - 1000)*0.44*Res
);
return
CdsRes
*pow(alpha2, -3.65)
*max(interface_.continuous(), interface_.continuous().residualAlpha());
return CdsRes*pow(alpha2, -2.65);
}

2
applications/solvers/modules/multiphaseEuler/interfacialModels/heatTransferModels/Gunn/Gunn.C
View File

@ -68,7 +68,7 @@ Foam::heatTransferModels::Gunn::K(const scalar residualAlpha) const
{
const volScalarField alpha2
(
max(1 - interface_.dispersed(), interface_.continuous().residualAlpha())
max(interface_.continuous(), interface_.continuous().residualAlpha())
);
const volScalarField sqrAlpha2(sqr(alpha2));