CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

rhoCentralFoam not starting after changing thermophysicalProperties to polynomial

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 11, 2021, 06:53
Default rhoCentralFoam not starting after changing thermophysicalProperties to polynomial
  #1
New Member
 
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5
janibo is on a distinguished road
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:
Code:
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:
Code:
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:
Code:
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

Last edited by janibo; August 11, 2021 at 10:58.
janibo is offline   Reply With Quote

Old   August 11, 2021, 10:29
Default
  #2
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 781
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
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?
dlahaye is offline   Reply With Quote

Old   August 11, 2021, 10:53
Default
  #3
New Member
 
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5
janibo is on a distinguished road
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.
janibo is offline   Reply With Quote

Old   August 11, 2021, 11:12
Default
  #4
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 781
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
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 ?
dlahaye is offline   Reply With Quote

Old   August 11, 2021, 11:33
Default
  #5
New Member
 
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5
janibo is on a distinguished road
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
Code:
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 is offline   Reply With Quote

Old   August 16, 2021, 06:32
Default
  #6
New Member
 
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5
janibo is on a distinguished road
So I played a little bit around with the DebugSwitches, namely those
for the thermoTypes.

log-file:
Code:
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.
janibo is offline   Reply With Quote

Reply

Tags
rhocentralfoam, thermophysical properties

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to apply a turbulence model to the rhoCentralFoam solver? arussell92 OpenFOAM Pre-Processing 18 July 25, 2022 05:26
Modify rhoCentralFoam: other equations of state fivos OpenFOAM Programming & Development 5 July 29, 2020 14:17
big difference between clockTime and executionTime LM4112 OpenFOAM Running, Solving & CFD 21 February 15, 2019 04:05
[OpenFOAM] XYZ starting points of streamtraces in Paraview deniggo ParaView 3 May 22, 2013 04:02
dynamic mesh refinement and rhoCentralFoam ChrisA OpenFOAM Running, Solving & CFD 1 March 21, 2013 09:00


All times are GMT -4. The time now is 07:19.