CFD Online Discussion Forums - double dot tensor product (double inner product) implementation

[RESOLVED]

Hello,

I'm trying to solve energy equation using a modified icoFoam solver.
I would like to take into account viscous dissipations in energy equation and to do that I need to calculate the double inner product between the viscous stress tensor and gradient of velocity.
I tried this:

Code:

//tensors computation

volTensorField gradU = fvc::grad(U);

volTensorField tau = nu * (gradU + gradU.T());



//solve the steady-state energy equation in temperature.

solve ( rho * Cp * fvm::div(phi, T) - fvm::laplacian(ka, T) - (tau && gradU));

but it doesn't compile leading to the error :

Code:

error: no match for ‘operator-’ in ‘Foam::operator-(const Foam::tmp<Foam::fvMatrix<Type> >&, const Foam::tmp<Foam::fvMatrix<Type> >&) [with Type = Foam::Vector<double>](((const Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > >&)((const Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > >*)(& Foam::fvm::laplacian(const Foam::dimensioned<Type2>&, Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = Foam::Vector<double>, GType = double](((Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&)(& T)))))))

- Foam::operator&&(const Foam::GeometricField<TypeR, PatchField, GeoMesh>&, const Foam::GeometricField<Type1, PatchField, GeoMesh>&) [with Type1 = Foam::Tensor<double>, Type2 = Foam::Tensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&)((const Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>*)(& gradU))))’

When I remove the term

Code:

- (tau && gradU)

of the formlation, the solver compiles perfectely...:confused:

So it's seems to be the double dot product implementation which doesn't work...

Any ideas ? (I use OpenFoam1.7.1)

Thanks in advance

Clément