songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bubbleFoam, twoPhaseEulerFoam: Updated handling of the transpose terms

Browse Source
This commit is contained in:
Henry 2011-03-21 10:49:20 +00:00
parent f663562ca0
commit 0f3d393ec5
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/multiphase/bubbleFoam/UEqns.H
View File

@ -2,7 +2,7 @@ fvVectorMatrix UaEqn(Ua, Ua.dimensions()*dimVol/dimTime);
fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime); fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
{ {
volTensorField Rca(-nuEffa*(fvc::grad(Ua)().T())); volTensorField Rca(-nuEffa*(T(fvc::grad(Ua))));
Rca = Rca + (2.0/3.0)*sqr(Ct)*I*k - (2.0/3.0)*I*tr(Rca); Rca = Rca + (2.0/3.0)*sqr(Ct)*I*k - (2.0/3.0)*I*tr(Rca);
surfaceScalarField phiRa surfaceScalarField phiRa
@ -36,7 +36,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
UaEqn.relax(); UaEqn.relax();
volTensorField Rcb(-nuEffb*fvc::grad(Ub)().T()); volTensorField Rcb(-nuEffb*T(fvc::grad(Ub)));
Rcb = Rcb + (2.0/3.0)*I*k - (2.0/3.0)*I*tr(Rcb); Rcb = Rcb + (2.0/3.0)*I*k - (2.0/3.0)*I*tr(Rcb);
surfaceScalarField phiRb surfaceScalarField phiRb

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

@ -3,7 +3,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
{ {
{ {
volTensorField gradUaT(fvc::grad(Ua)().T()); volTensorField gradUaT(T(fvc::grad(Ua)));
if (kineticTheory.on()) if (kineticTheory.on())
{ {
@ -58,7 +58,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
} }
{ {
volTensorField gradUbT(fvc::grad(Ub)().T()); volTensorField gradUbT(T(fvc::grad(Ub)));
volTensorField Rcb volTensorField Rcb
( (
"Rcb", "Rcb",