[sreehahaha]

Hi Foamers,

First of all, thank you for all the support I've been receiving so far. Many suggestions from the forum have been incredibly helpful for me.

To the issue: I'm simulating a compressible system using a solver modified from rhoPimpleFoam on OF9. It ran successfully when I used a "thermophysicalProperties" file like:

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          eConst;
    equationOfState rPolynomial;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture                                                                        
{                                                                              
    specie                                                                     
    {                                                                          
        molWeight       18.0;                                                  
    }                                                                          
    thermodynamics                                                             
    {                              
	Cv		4233.7894;                                            
        Sf              0;                                                     
        Hf              0;                                                     
    }                                                                          
    equationOfState                                                            
    {                                                                          
        C (0.001866 -6.186e-06 1.105e-08 5.697e-13 0);               
    }                                                                          
    transport                                                                  
    {                                                                       
	mu	5.0e-03;
	Pr	26.4;   
    }                                                                          
}

But when I change the equationOfState to rhoTabulated and give rho as a table, it gives me the following error:

Using dynamicCode for codedBase heRhoThermo_pureMixture_const_eConst_rhoTabulated_specie___sensibleInternalEnergy___ at line 17 in "/desktop/2018/software/OpenFOAM/9-foss-2021b/OpenFOAM-9/etc/codeTemplates/dynamicCode/fluidThermo"
[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Not implemented
[1] 
[1]     From function Foam::scalar Foam::rhoTabulated<Specie>::E(Foam::scalar, Foam::scalar) const [with Specie = Foam::specie; Foam::scalar = double]
[1]     in file /desktop/2018/software/OpenFOAM/9-foss-2021b/OpenFOAM-9/src/thermophysicalModels/specie/lnInclude/rhoTabulatedI.H at line 122.
[1] 
FOAM parallel run aborting

Can someone please help me understand what's going on? I'm confused, because rPolynomial gives rho as a function of (p,T), and so does rhoTabulated.

Many thanks,
SP