research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

continuousGasKEpsilon: Omega now consistent with Lahey 2005 paper

Browse source
This commit is contained in:
Henry 2014-03-26 12:48:37 +00:00
parent 7cc1ed27b3
commit c957f2450d
1 changed files with 11 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/TurbulenceModels/phaseIncompressible/RAS/continuousGasKEpsilon/continuousGasKEpsilon.C
View file

@ -124,9 +124,17 @@ void continuousGasKEpsilon<BasicTurbulenceModel>::correctNut()
const transportModel& liquid = fluid.otherPhase(gas);
volScalarField thetal(liquidTurbulence.k()/liquidTurbulence.epsilon());
volScalarField thetag((1.0/(18*liquid.nu()))*sqr(gas.d()));
volScalarField expThetar(exp(min(thetal/thetag, scalar(50))));
volScalarField omega(sqr(expThetar - 1)/(sqr(expThetar) - 1));
volScalarField rhodv(gas.rho() + fluid.virtualMass(gas).Cvm()*liquid.rho());
volScalarField thetag((rhodv/(18*liquid.rho()*liquid.nu()))*sqr(gas.d()));
volScalarField expThetar
(
min
(
exp(min(thetal/thetag, scalar(50))),
scalar(1)
)
);
volScalarField omega((1 - expThetar)/(1 + expThetar));
nutEff_ = omega*liquidTurbulence.nut();
}