[shock77]

Hi,


I would like to calculate the specific heat capacity cp, specific heat ratio gamma and Ma number during the computation.


I have compiled a new solver based on rhoCentralFoam for this purpose.
In createFields.H I have added:

volScalarField cp
    (
    IOobject
        (
            "cp",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
        ),
         thermo.Cp()
     );
     
volScalarField gamma
    (
    IOobject
        (
            "gamma",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
        ),
         thermo.gamma()
     );
     
volScalarField Ma
    (
    IOobject
        (
            "Ma",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
        ),
         mag(U)/sqrt(thermo.gamma()*thermo.p()/thermo.rho())
     );

I get the right values at time t = 0. Unfortunately they do not change over time, which is strange. Pressure, velocity and density do change, so the Mach Number should change aswell, but I cant see this. It works perfectly, if I use the functionObjects for the Mach Number.


I wanted to add also the total temperature and total pressure for compressible fluids to be calculated during runtime, but I would like to calculate them using the Mach Number, which is not calculated right at the moment.




Can anyone tell what I did wrong?






Kind regards,
shock77