[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

[Tobermory]

The error message tells you that rhoTabulated<Specie>::E() is not implemented. Sure enough, if you look in the class inline source file (https://www.cfd-online.com/Forums/ne...ote=1&p=846511) you'll see:

template<class Specie>
inline Foam::scalar Foam::rhoTabulated<Specie>::E
 (
     scalar p,
     scalar T
 ) const
 {
     NotImplemented;
     return 0;
 }

... i.e. E is not implemented in this class. Infact you'll see that a whole bunch of thernodynamic functions are also not implemented, so this is probably not a good EOS for you to use.

[sreehahaha]

Thanks Tobermory for the quick reply.

Could you please suggest how I rectify this issue? I'm new to using OF. I think I have to use rhoTabulated as the rPolynomial I've been using isn't very accurate.

SP

[Tobermory]

Sreehari

The short answer is that you can't use rhoTabulated at the moment - it doesn't do what you need it to do, i.e. solve the energy equation. I think your best bet is to try and improve the polynomial fit as best you can, and continue using rPolynomial.

Otherwise, you could try and implement your own version of the rhoTabulated class, and implement the energy functions ... however, there may be a good reason why they have not been implemented so far (I suspect that it's not just laziness on the Developer's side), and so this may be a dead end - I don't have enough knowledge of these particular EOS's to be able to comment.

Good luck!

[sreehahaha]

Thank you very much, Tobermory.

SP