[janibo]

Dear Foamers,

I am able to run a cd-nozzle case with the fluid air with the following thermophysicalProperties
successfully and with good results:

original thermophysicalProperties-file:

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo         	hConst;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture	//Air
{
    specie
    {
	nMoles			1;		//mol
        molWeight       28.96; // g/mol // kg/kmol
    }
    thermodynamics
    {
	
        Cp              1004.5;	
        Hf              2.544e+06;
    }
    transport
    {
	As  		1.4584e-6;
	Ts  		110.33;
        mu           1.8e-5;
        Pr             0.71;
    }
	equationOfState
	{
		Tc  132.52; 	//K
		Vc  0.08443;	//m^3/kmol
		Pc  37.56e05;	//Pa
		omega  0.038;	//dimensionless
	}
}

I changed the thermodynamics and transport to
the polynomial options and everything else stayed the same:

thermophysicalProperties-file changed to polynomial:

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo         	hPolynomial;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture	//Air
{
    specie
    {
	nMoles	      1;		//mol
        molWeight       28.96; // g/mol // kg/kmol
    }
    thermodynamics
    {
	//Cv				718;	//	J/g*K
        //Cp                         1004.5;	//	J/g*K
        Hf                             0; //2.544e+06;
	Sf				0.0;
	CpCoeffs<8>		(1.3211 -0.0047 3e-5 -8e-8 1e-10 -1e-13 4e-17 0); //based on [J/kg*K]
    }
    transport
    {
	//As  			1.4584e-6;
	//Ts  			110.33;
        //mu                         1.8e-5;
        //Pr                           0.71;
	muCoeffs<8>     ( 6e-7 7e-8 -4e-11 8e-15 8e-18 0 0 0 ); //based on [Pa s]
        kappaCoeffs<8>  ( -0.0018 0.0001 -1e-7 9e-11 -3e-14 0 0 0 ); // based on [W/mk]
    }
	equationOfState
	{
		Tc  132.52; 	//K
		Vc  0.08443;	//m^3/kmol
		Pc  37.56e05;	//Pa
		omega  0.038;	//dimensionless
	}
}

After the change I have the issue that rhoCentralFoam does not start at all,
although all information are read and no error occurs:

log-file:

Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    RASModel        kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
    decayControl    false;
    kInf            0;
    omegaInf        0;
}

fluxScheme: Kurganov

Starting time loop

I took the handbook as guidance for the declaration of variables. The coefficients are taken from the VDI heat atlas.

Does anyone have a clue or might point me in the right dircetion?

Thanks in advance. Niklas

[dlahaye]

1/ Does rhoCentralFoam use the constant/thermoPhysicalProperties file at all? That is, do you observe any any changes in the convergence and/or final solution when making changes in the constant/thermoPhysicalProperties file?

2/ What exactly do you mean by rhoCentralFoam not starting. Does rhoCentralFoam print message that we can learn from?

[janibo]

Hi Domenico,

thanks for the reply.

1) Yes, rhoCentralFoam uses the thermophysicalProperty file as shown
in the first file code, that one worls like a charm with nice results.

2) No unfortunetly not. The posted log-file is all I get. rhoCentralFoam
just exits the simulation after "Starting time loop", but does not yield
an error message.

[dlahaye]

Does it help to increase the verbose by increasing the value of the debug switches? See e.g. line 87 to line 96 of https://github.com/OpenFOAM/OpenFOAM...tc/controlDict ?

[janibo]

So interestingly my etc\controlDict file yields around 1000 DebugSwitches,
i'm using v2012, unlike the controlDict from the docs you linked.

I am able to find entries in the etc\controlDict like

pureMixture<constTransport<thermo<hConstThermo<perfectGas>>>> 0;

but i cannot find an entry fitting the selection I have.

If I would have choosen a wrong combination of options,
OpenFOAM would have yield an error accordingly.

I picked my options from the list of valid thermoTypes.

If I had declared a varaible wrong OpenFOAM would yield an error accordingly, too.

So that is what really confuses me.

[janibo]

So I played a little bit around with the DebugSwitches, namely those
for the thermoTypes.

log-file:

Create time

Overriding DebugSwitches according to controlDict
    dimensionSet    2;

    fileName (unregistered)
    GAMGAgglomeration 2;

    level (unregistered)
    lduMatrix       2;

    SolverPerformance 2;

    vtkUnstructuredReader 2;

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    RASModel        kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
    decayControl    false;
    kInf            0;
    omegaInf        0;
}

fluxScheme: Kurganov

Starting time loop

The only thing that changed is the overwriting of the DebugSwitches according to the controlDict, but other than that nothing changed.