From 0f3d393ec528eb63b3741c11c4cc890d58caffe6 Mon Sep 17 00:00:00 2001 From: Henry <Henry> Date: Mon, 21 Mar 2011 10:49:20 +0000 Subject: [PATCH] bubbleFoam, twoPhaseEulerFoam: Updated handling of the transpose terms --- applications/solvers/multiphase/bubbleFoam/UEqns.H | 4 ++-- applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/applications/solvers/multiphase/bubbleFoam/UEqns.H b/applications/solvers/multiphase/bubbleFoam/UEqns.H index 82d53d84b00..3239df21cb3 100644 --- a/applications/solvers/multiphase/bubbleFoam/UEqns.H +++ b/applications/solvers/multiphase/bubbleFoam/UEqns.H @@ -2,7 +2,7 @@ fvVectorMatrix UaEqn(Ua, Ua.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); surfaceScalarField phiRa @@ -36,7 +36,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime); 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); surfaceScalarField phiRb diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H index b07f3e70dd3..0a590e6bf93 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H @@ -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()) { @@ -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 ( "Rcb", -- GitLab