[atoof]

Dears,

I am going to implement a new model for calculating viscosity that needs to heat capacity of each species. The new solver is based on reactingFoam.

Is it possible to get heat capacity of all species in all computational cells of the domain?

I used the following with no success.

volScalarField cpi = composition.Cp(i,T);

where
i is the label of species and

basicMultiComponentMixture& composition = thermo.composition();

Thank you in advance,

Hossein