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

PDRkEpsilon: Removed second bounding of epsilon to avoid turbulence-combustion coupling instability

Browse Source
This commit is contained in:
Henry Weller
2024-06-26 15:59:20 +01:00
parent 8bf9afff11
commit b8ad3a9036
2 changed files with 7 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
applications/legacy/combustion/PDRFoam/PDRModels/turbulence/PDRkEpsilon/PDRkEpsilon.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) 2011-2023 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2024 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -46,16 +46,16 @@ addToRunTimeSelectionTable(RASModel, PDRkEpsilon, dictionary);
tmp<volScalarField> PDRkEpsilon::boundEpsilon() tmp<volScalarField> PDRkEpsilon::boundEpsilon()
{ {
tmp<volScalarField> tCmuk2(Cmu_*sqr(k_)); tmp<volScalarField> tCmuk2(Cmu_*sqr(k_));
epsilon_ = max(epsilon_, tCmuk2()/(this->nutMaxCoeff_*this->nu())); epsilon_ = max(epsilon_, tCmuk2()/(nutMaxCoeff_*nu()));
return tCmuk2; return tCmuk2;
} }
void PDRkEpsilon::correctNut() void PDRkEpsilon::correctNut()
{ {
this->nut_ = boundEpsilon()/epsilon_; nut_ = Cmu_*sqr(k_)/epsilon_;
this->nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
fvConstraints::New(this->mesh_).constrain(this->nut_); fvConstraints::New(mesh_).constrain(nut_);
} }

4
tutorials/legacy/combustion/PDRFoam/flamePropagationWithObstacles/system/controlDict
View File

@ -16,7 +16,7 @@ FoamFile
application PDRFoam; application PDRFoam;
startFrom latestTime; startFrom startTime;
startTime 0; startTime 0;
@ -32,7 +32,7 @@ writeInterval 0.01;
purgeWrite 0; purgeWrite 0;
writeFormat ascii; writeFormat binary;
writePrecision 8; writePrecision 8;