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

TurbulenceModels: Reorder stress component evaluation

Browse Source 
This ensures the U boundary conditions are updated before the explicit
parts of the stress are evaluated.

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1988
This commit is contained in:
Henry Weller
2016-02-04 20:23:22 +00:00
parent 7265d90d0d
commit 8e04a529d9
2 changed files with 17 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

28
src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2015-2016 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -238,43 +238,43 @@ Foam::ReynoldsStress<BasicTurbulenceModel>::divDevRhoReff
{
return
(
fvc::div
fvc::laplacian
(
(1.0 - couplingFactor_)*this->alpha_*this->rho_*this->nut(),
U,
"laplacian(nuEff,U)"
)
+ fvc::div
(
this->alpha_*this->rho_*R_
+ couplingFactor_
*this->alpha_*this->rho_*this->nut()*fvc::grad(U),
"div(devRhoReff)"
)
+ fvc::laplacian
(
(1.0 - couplingFactor_)*this->alpha_*this->rho_*this->nut(),
U,
"laplacian(nuEff,U)"
)
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
- fvc::div(this->alpha_*this->rho_*this->nu()*dev2(T(fvc::grad(U))))
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
);
}
else
{
return
(
fvc::div(this->alpha_*this->rho_*R_)
+ fvc::laplacian
fvc::laplacian
(
this->alpha_*this->rho_*this->nut(),
U,
"laplacian(nuEff,U)"
)
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
+ fvc::div(this->alpha_*this->rho_*R_)
- fvc::div(this->alpha_*this->rho_*this->nu()*dev2(T(fvc::grad(U))))
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
);
}
return
(
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
- fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U))))
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
);
}
@ -289,8 +289,8 @@ Foam::ReynoldsStress<BasicTurbulenceModel>::divDevRhoReff
{
return
(
- fvm::laplacian(this->alpha_*rho*this->nuEff(), U)
- fvc::div((this->alpha_*rho*this->nuEff())*dev2(T(fvc::grad(U))))
- fvm::laplacian(this->alpha_*rho*this->nuEff(), U)
);
}

6
src/TurbulenceModels/turbulenceModels/linearViscousStress/linearViscousStress.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2013-2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -97,8 +97,8 @@ Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff
{
return
(
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
- fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U))))
- fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
);
}
@ -113,8 +113,8 @@ Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff
{
return
(
- fvm::laplacian(this->alpha_*rho*this->nuEff(), U)
- fvc::div((this->alpha_*rho*this->nuEff())*dev2(T(fvc::grad(U))))
- fvm::laplacian(this->alpha_*rho*this->nuEff(), U)
);
}