[james.t]

Hello,

in the k-epsilon turbulence model the total stress term is calculated as volSymmTensorField

tmp<volSymmTensorField> kEps::devRhoReff() const
{
    return tmp<volSymmTensorField>
    (
        new volSymmTensorField
        (
            IOobject
            (
                "devRhoReff",
                runTime_.timeName(),
                mesh_,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
           -muEff()*dev(twoSymm(fvc::grad(U_)))
        )
    );
}

and the source term for the momentum equation as fvVectorMatrix

tmp<fvVectorMatrix> kEps::divDevRhoReff(volVectorField& U) const
{
    return
    (
      - fvm::laplacian(muEff(), U)
      - fvc::div(muEff()*dev2(fvc::grad(U)().T()))
    );
}

Now I'd like to modify the total stress tensor devRhoReff and use it for the calculation of the source term of the momentum equation divDevRhoReff like this:

fvVectorMatrix divDevRhoReff = fvc::div(devRhoReff);

However the following error occurs:

error: conversion from ‘Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >’ to non-scalar type ‘Foam::fvVectorMatrix’ requested

Is there another way to convert the total stress tensor devRhoReff to the source term divDevRhoReff?

Thanks alot!

James