zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

driftFluxFoam: Changed the laminar viscosity used in the k-epsilon model that of the mixture

Browse Source
This commit is contained in:
Henry
2014-03-26 12:49:18 +00:00
committed by Andrew Heather
parent 5ba7f46a20
commit 7640ec88b0
4 changed files with 16 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
View File

@ -87,7 +87,6 @@ int main(int argc, char *argv[])
#include "alphaEqnSubCycle.H"
twoPhaseProperties.correct();
Info<< average(twoPhaseProperties.mu()) << endl;
#include "UEqn.H"

20
applications/solvers/multiphase/driftFluxFoam/kEpsilon.H
View File

@ -10,8 +10,6 @@ if (turbulence)
dimensionedScalar epsilon0("epsilon0", epsilon.dimensions(), 0);
dimensionedScalar epsilonMin("epsilonMin", epsilon.dimensions(), SMALL);
volScalarField divU(fvc::div(phi));
tmp<volTensorField> tgradU = fvc::grad(U);
volScalarField G(mut*(tgradU() && dev(twoSymm(tgradU()))));
tgradU.clear();
@ -21,7 +19,7 @@ if (turbulence)
Cmu*k/sigmak*(g & fvc::grad(rho))/(epsilon + epsilonMin)
);
volScalarField muc(twoPhaseProperties.nucModel().nu()*rho2);
volScalarField mul(twoPhaseProperties.mu());
#include "wallFunctions.H"
@ -32,12 +30,12 @@ if (turbulence)
+ fvm::div(rhoPhi, epsilon)
- fvm::laplacian
(
mut/sigmaEps + muc, epsilon,
mut/sigmaEps + mul, epsilon,
"laplacian(DepsilonEff,epsilon)"
)
==
C1*G*epsilon/(k + kMin)
- fvm::SuSp(C1*(1.0 - C3)*Gcoef + (2.0/3.0*C1)*rho*divU, epsilon)
- fvm::SuSp(C1*(1.0 - C3)*Gcoef, epsilon)
- fvm::Sp(C2*rho*epsilon/(k + kMin), epsilon)
);
@ -56,12 +54,12 @@ if (turbulence)
+ fvm::div(rhoPhi, k)
- fvm::laplacian
(
mut/sigmak + muc, k,
mut/sigmak + mul, k,
"laplacian(DkEff,k)"
)
==
G
- fvm::SuSp(Gcoef + 2.0/3.0*rho*divU, k)
- fvm::SuSp(Gcoef, k)
- fvm::Sp(rho*epsilon/(k + kMin), k)
);
@ -75,6 +73,10 @@ if (turbulence)
mut = rho*Cmu*sqr(k)/(epsilon + epsilonMin);
#include "wallViscosity.H"
}
muEff = mut + twoPhaseProperties.mu();
muEff = mut + mul;
}
else
{
muEff = mut + twoPhaseProperties.mu();
}

4
applications/solvers/multiphase/driftFluxFoam/wallFunctions.H
View File

@ -33,7 +33,7 @@
if (isA<wallFvPatch>(curPatch))
{
const scalarField& mutw = mut.boundaryField()[patchi];
const scalarField& mucw = muc.boundaryField()[patchi];
const scalarField& mulw = mul.boundaryField()[patchi];
scalarField magFaceGradU
(
@ -55,7 +55,7 @@
/(kappa_*y[patchi][facei]);
G[faceCelli] +=
(mutw[facei] + mucw[facei])
(mutw[facei] + mulw[facei])
*magFaceGradU[facei]
*Cmu25*::sqrt(k[faceCelli])
/(kappa_*y[patchi][facei]);

6
applications/solvers/multiphase/driftFluxFoam/wallViscosity.H
View File

@ -12,7 +12,7 @@
if (isA<wallFvPatch>(curPatch))
{
scalarField& mutw = mut.boundaryField()[patchi];
const scalarField& mucw = muc.boundaryField()[patchi];
const scalarField& mulw = mul.boundaryField()[patchi];
forAll(curPatch, facei)
{
@ -20,12 +20,12 @@
scalar yPlus =
Cmu25*y[patchi][facei]*::sqrt(k[faceCelli])
/(mucw[facei]/rho2.value());
/(mulw[facei]/rho2.value());
if (yPlus > 11.6)
{
mutw[facei] =
mucw[facei]*(yPlus*kappa_/::log(E_*yPlus) - 1);
mulw[facei]*(yPlus*kappa_/::log(E_*yPlus) - 1);
}
else
{