Skip to content
Snippets Groups Projects
Commit cc6dbfa4 authored by andy's avatar andy
Browse files

BUG: Corrected temporary field usage - mantis #775

parent b45a4486
Branches
Tags
No related merge requests found
...@@ -406,9 +406,11 @@ void v2f::correct() ...@@ -406,9 +406,11 @@ void v2f::correct()
"v2f::alpha", "v2f::alpha",
1.0/T*((C1_ - N)*v2_ - 2.0/3.0*k_*(C1_ - 1.0)) 1.0/T*((C1_ - N)*v2_ - 2.0/3.0*k_*(C1_ - 1.0))
); );
const volScalarField Ceps1
tmp<volScalarField> Ceps1 = (
1.4*(1.0 + 0.05*min(sqrt(k_/v2_), scalar(100.0))); "Ceps1",
1.4*(1.0 + 0.05*min(sqrt(k_/v2_), scalar(100.0)))
);
// Update epsilon (and possibly G) at the wall // Update epsilon (and possibly G) at the wall
epsilon_.boundaryField().updateCoeffs(); epsilon_.boundaryField().updateCoeffs();
...@@ -420,7 +422,7 @@ void v2f::correct() ...@@ -420,7 +422,7 @@ void v2f::correct()
+ fvm::div(phi_, epsilon_) + fvm::div(phi_, epsilon_)
- fvm::laplacian(DepsilonEff(), epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_)
== ==
Ceps1()*G/T Ceps1*G/T
- fvm::SuSp(((2.0/3.0)*Ceps1 + Ceps3_)*rho_*divU, epsilon_) - fvm::SuSp(((2.0/3.0)*Ceps1 + Ceps3_)*rho_*divU, epsilon_)
- fvm::Sp(Ceps2_*rho_/T, epsilon_) - fvm::Sp(Ceps2_*rho_/T, epsilon_)
); );
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment