[jeroendemto]

Dear OpenFOAMers

I'm trying to implement conjugate heat transfer into my phase change model, however I have a problem when implementing a volScalarField rho, which depends on the volScalarField T (temperature).

This is the code I use to define the temperature field:

PtrList<volScalarField> TFluid(fluidRegions.seize());

TFluid.set
(
    i,
    new volScalarField
    (
        IOobject
        (
                "T",
                runTime.timeName(),
                fluidRegions[i],
                IOobject::MUST_READ,
                IOobject::AUTO_WRITE
        ),
        fluidRegions[i]
    )
);

When I try to calculate the rho (refRhoLiquid*(scalar(1)-beta*(T-T0))) it gives the following error:

error: no match for 'opertor-' (operand types are '<unresolved overloaded function type>' and 'FOAM::dimesnionedScalar'

This means I did not define the temperature (T) correctly, but I can't see what went wrong with that?

If you need additional information feel free to ask.

regards,
Jeroen