[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

[mAlletto]

The problem is you try to subtract from a scalar 1 a dimensionised scalar beta. The operation is not defined anywhere. You need either to write beta.value() or assign the correct dimension to 1

[jeroendemto]

Thank you for your response, but I think the issue is not the beta function, as this worked before like this:

volScalarField rho
(
    IOobject
    (
        "rho",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    rhoLiquid*(scalar(1) - beta*(T - T0))
);

I also added an image to show the exact error, it shows that the problem is with T-T_0. T should be a volScalarField but is not.

[jeroendemto]

Found the solution it should be TFluid[i] - T0 instead of T-T0

[mAlletto]

Can you post the code snippet which is failing

[jeroendemto]

This is the original

rhoFluid.set
    (
    	    i,
	    new volScalarField
	    (
		    IOobject
		    (
			"rhoLiquid",
			runTime.timeName(),
			fluidRegions[i],
			IOobject::READ_IF_PRESENT,
			IOobject::AUTO_WRITE
		    ),
		    refRhoLiquid*(scalar(1) - beta*(T - T0))
	    )
    );

It should be:

rhoFluid.set
    (
    	    i,
	    new volScalarField
	    (
		    IOobject
		    (
			"rhoLiquid",
			runTime.timeName(),
			fluidRegions[i],
			IOobject::READ_IF_PRESENT,
			IOobject::AUTO_WRITE
		    ),
		    refRhoLiquid*(scalar(1) - beta*(TFluid[i] - T0))
	    )
    );