|
[Sponsors] |
rhoCentralFoam not starting after changing thermophysicalProperties to polynomial |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 11, 2021, 06:53 |
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 |
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 } } 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 } } 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 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. |
|
August 11, 2021, 10:29 |
|
#2 |
Senior Member
|
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? |
|
August 11, 2021, 10:53 |
|
#3 |
New Member
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5 |
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. |
|
August 11, 2021, 11:12 |
|
#4 |
Senior Member
|
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 ?
|
|
August 11, 2021, 11:33 |
|
#5 |
New Member
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5 |
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; 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. |
|
August 16, 2021, 06:32 |
|
#6 |
New Member
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5 |
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 |
|
Tags |
rhocentralfoam, thermophysical properties |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |