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

BUG: ID 0000028 Inconsistency in nuEffa declaration

Browse Source
This commit is contained in:
sergio
2010-09-07 09:38:11 +01:00
parent 4e7fa4989b
commit 08a9c9da8f
1 changed files with 7 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
View File

@ -4,11 +4,6 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
{ {
{ {
volTensorField gradUaT = fvc::grad(Ua)().T(); volTensorField gradUaT = fvc::grad(Ua)().T();
volTensorField Rca
(
"Rca",
((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
);
if (kineticTheory.on()) if (kineticTheory.on())
{ {
@ -20,6 +15,13 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
nuEffa = sqr(Ct)*nutb + nua; nuEffa = sqr(Ct)*nutb + nua;
} }
volTensorField Rca
(
"Rca",
((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
);
if (kineticTheory.on()) if (kineticTheory.on())
{ {
Rca -= ((kineticTheory.lambda()/rhoa)*tr(gradUaT))*tensor(I); Rca -= ((kineticTheory.lambda()/rhoa)*tr(gradUaT))*tensor(I);