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

reactingEulerFoam::dragModels::segregated: Corrected ReI

Browse Source
This commit is contained in:
Henry Weller
2019-08-09 16:56:17 +01:00
parent 61f9131389
commit 52cb55a27c
3 changed files with 35 additions and 26 deletions
Download Patch File Download Diff File Expand all files Collapse all files

36
applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/segregated/segregated.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2014-2018 OpenFOAM Foundation \\ / A nd | Copyright (C) 2014-2019 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -107,7 +107,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
L.primitiveFieldRef() = cbrt(mesh.V()); L.primitiveFieldRef() = cbrt(mesh.V());
L.correctBoundaryConditions(); L.correctBoundaryConditions();
volScalarField I const volScalarField I
( (
alpha1 alpha1
/max /max
@ -116,7 +116,8 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
pair_.phase1().residualAlpha() + pair_.phase2().residualAlpha() pair_.phase1().residualAlpha() + pair_.phase2().residualAlpha()
) )
); );
volScalarField magGradI
const volScalarField magGradI
( (
max max
( (
@ -125,28 +126,35 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
) )
); );
volScalarField muI const volScalarField muI
( (
rho1*nu1*rho2*nu2 rho1*nu1*rho2*nu2
/(rho1*nu1 + rho2*nu2) /(rho1*nu1 + rho2*nu2)
); );
volScalarField muAlphaI
const volScalarField limitedAlpha1
(
max(alpha1, pair_.phase1().residualAlpha())
);
const volScalarField limitedAlpha2
(
max(alpha2, pair_.phase2().residualAlpha())
);
const volScalarField muAlphaI
( (
alpha1*rho1*nu1*alpha2*rho2*nu2 alpha1*rho1*nu1*alpha2*rho2*nu2
/( /(limitedAlpha1*rho1*nu1 + limitedAlpha2*rho2*nu2)
max(alpha1, pair_.phase1().residualAlpha())*rho1*nu1
+ max(alpha2, pair_.phase2().residualAlpha())*rho2*nu2
)
); );
volScalarField ReI const volScalarField ReI
( (
pair_.rho() pair_.rho()*pair_.magUr()
*pair_.magUr() /(magGradI*limitedAlpha1*limitedAlpha2*muI)
/(magGradI*muI)
); );
volScalarField lambda(m_*ReI + n_*muAlphaI/muI); const volScalarField lambda(m_*ReI + n_*muAlphaI/muI);
return lambda*sqr(magGradI)*muI; return lambda*sqr(magGradI)*muI;
} }

21
applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/segregated/segregated.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2014-2018 OpenFOAM Foundation \\ / A nd | Copyright (C) 2014-2019 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -107,7 +107,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
L.primitiveFieldRef() = cbrt(mesh.V()); L.primitiveFieldRef() = cbrt(mesh.V());
L.correctBoundaryConditions(); L.correctBoundaryConditions();
volScalarField I const volScalarField I
( (
alpha1 alpha1
/max /max
@ -116,7 +116,8 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
pair_.phase1().residualAlpha() + pair_.phase2().residualAlpha() pair_.phase1().residualAlpha() + pair_.phase2().residualAlpha()
) )
); );
volScalarField magGradI
const volScalarField magGradI
( (
max max
( (
@ -125,25 +126,25 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
) )
); );
volScalarField muI const volScalarField muI
( (
rho1*nu1*rho2*nu2 rho1*nu1*rho2*nu2
/(rho1*nu1 + rho2*nu2) /(rho1*nu1 + rho2*nu2)
); );
volScalarField muAlphaI
const volScalarField muAlphaI
( (
alpha1*rho1*nu1*alpha2*rho2*nu2 alpha1*rho1*nu1*alpha2*rho2*nu2
/(alpha1*rho1*nu1 + alpha2*rho2*nu2) /(alpha1*rho1*nu1 + alpha2*rho2*nu2)
); );
volScalarField ReI const volScalarField ReI
( (
pair_.rho() pair_.rho()*pair_.magUr()
*pair_.magUr() /(magGradI*alpha1*alpha2*muI)
/(magGradI*muI)
); );
volScalarField lambda(m_*ReI + n_*muAlphaI/muI); const volScalarField lambda(m_*ReI + n_*muAlphaI/muI);
return lambda*sqr(magGradI)*muI; return lambda*sqr(magGradI)*muI;
} }

4
tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnIATE/system/setFieldsDict
View File

@ -17,7 +17,7 @@ FoamFile
defaultFieldValues defaultFieldValues
( (
volScalarFieldValue alpha1 1 volScalarFieldValue alpha.air 1
); );
regions regions
@ -27,7 +27,7 @@ regions
box (0 0 -0.1) (0.15 0.701 0.1); box (0 0 -0.1) (0.15 0.701 0.1);
fieldValues fieldValues
( (
volScalarFieldValue alphaair 0 volScalarFieldValue alpha.air 0
); );
} }
); );