[Gian Maria]

Hi foamers,

I'm trying to improve a new solver in OpenFOAM, but I have some problems.
I need to calculate Cp (Heat Capacity at Constant Pressure), and I must add it in the standard solver "rhoPorousSimpleFoam".

I tried to add a line in "createFields.H", but it doesn't work.

// these lines are already present

autoPtr<basicPsiThermo> pThermo
(
basicPsiThermo::New(mesh)
);

basicPsiThermo& thermo = pThermo();

volScalarField& p = thermo.p();
volScalarField& h = thermo.h();
const volScalarField& psi = thermo.psi();

// this is my add
volScalarField& Cp = thermo.Cp();

In "basicPsiThermo" Cp is called with "T" and "patchi" as arguments. Must I do the same?

Or can I load properties of air by an external file?

[l_r_mcglashan]

What does the error say? Is it a compiler error or a runtime error? You will probably find the mistake is in the type you have given for the variable Cp. Check it against the return type of Cp in basicThermo.

[Gian Maria]

Hi Laurence, thanks for your reply.

So, I will explain you what is my problem.

I want to create a new solutor, and I need for my calculations some variables:
- enthalpy
- pressure
- heat costant Cp
- rho

But whenever h, p and rho are already gived, I must find the way to calculate Cp.

My "createFields.H" file start with these lines:

Info<< "Reading thermophysical properties\n" << endl;

autoPtr<basicPsiThermo> pThermo
(
basicPsiThermo::New(mesh)
);

basicPsiThermo& thermo = pThermo();

volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
thermo.rho()
);

volScalarField& p = thermo.p();
volScalarField& h = thermo.h();
const volScalarField& psi = thermo.psi();

so I tried to load Cp value using the same templates:

volScalarField& Cp = thermo.Cp();

But when I try to compile I have this error:

createFields.H:27: error: invalid initialization of non-const reference of type ‘Foam::volScalarField&’ from a temporary of type ‘Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’

Why compare this error?
How can I move to correct it?

PS: if I check in "basicPsiThermo" file, I find this

// Fields derived from thermodynamic state variables

//- Temperature [K]
virtual const volScalarField& T() const;

//- Heat capacity at constant pressure for patch [J/kg/K]
virtual tmp<scalarField> Cp
(
const scalarField& T,
const label patchi
) const;

Can anyone explain me how to use this information, plese?

Thanks in advance!

[l_r_mcglashan]

In short: Change 'volScalarField& Cp = thermo.Cp()' to 'volScalarField Cp = thermo.Cp();'

[alimea]

Hi foamers,

Why does cp(p,T) function return zero for ideal gas?

How can I calculate it for other equation of states?

Thanks

[labyrinth01]

Originally Posted by alimea

Hi foamers,

Why does cp(p,T) function return zero for ideal gas?

How can I calculate it for other equation of states?

Thanks

Hey, Cp=Cp.Janaf+Cp.EOS. It includes two parts: one is calculated using the Janaf relationship and the other one is calculated using the EOS. For ideal gas equation, the EOS part is zero. But when you are using the real-gas EOS, it is not zero and there is a departure function from the reference state.