Entropy calculation in OpenFoam

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

 September 8, 2013, 15:32 #21 Super Moderator   Bruno Santos Join Date: Mar 2009 Location: Lisbon, Portugal Posts: 9,507 Blog Entries: 36 Rep Power: 96 OpenFOAM hasn't implemented all possible combinations of thermodynamic options. You have to first check the list of valid combinations, as implied here: http://www.openfoam.org/version2.2.0....php#x10-40000 - section "Dictionary-based Thermodynamics Selection". In other words: did you first run the solver with those options in the dictionary file, to confirm if it works? __________________ OpenFOAM: FAQ | Getting started Forum: How to get help and to post code/output What am I doing/planning: blog/wiki Read this before sending me PM ___ I'll be at OFW11 in Portugal

 September 8, 2013, 16:09 #22 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,209 Rep Power: 18 Hi Bruno yes that case is finished 10 days before and I'm engaging with making more than 150 plots,and its the thermodynamic dictionary: Code: ```thermoType { type hePsiThermo; mixture pureMixture; transport const; thermo hConst;//janaf;// equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } mixture { specie { nMoles 1; molWeight 28.96; } thermodynamics { /*Tlow 200; Thigh 6000; Tcommon 1000; highCpCoeffs (3.129672277499967 1.211766752750042e-3 -4.04761109000018e-7 6.409571300000324e-11 -3.780447620000211e-15 -.996955811e3 5.257756413); lowCpCoeffs (3.593868050000014 -8.419292810001117e-4 2.08457592775033e-6 -5.938441825004294e-10 -2.456771699998039e-13 -1.0512181e3 3.140021718);*/ Cp 1100; Hf 2.544e+06; } transport { /*As 1.4584e-6; Ts 110.4;*/ mu 0;//1.84e-05 Pr 0.75; } }``` maybe I have converted it by some mistakes. __________________ Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked.

September 8, 2013, 16:37
#23
Super Moderator

Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,507
Blog Entries: 36
Rep Power: 96
Hi Ehsan,

Unfortunately OpenFOAM doesn't have the calculation of entropy implemented for the "hConstThermo":
Quote:
 From file "\$FOAM_SRC/thermophysicalModels/specie/thermo/hConst/hConstThermoI.H": Code: ```template inline Foam::scalar Foam::hConstThermo::s ( const scalar p, const scalar T ) const { notImplemented ( "scalar hConstThermo::" "s(const scalar p, const scalar T) const" ); return T; }```
If you run this command:
Code:
`ls -l \$FOAM_SRC/thermophysicalModels/specie/thermo`
You should see some of the possible options for the "thermo" option.

After looking at the source codes, only the following options have an implemented entropy calculation:
Code:
```janaf
hPolynomial```
Best regards,
Bruno
__________________
___
I'll be at OFW11 in Portugal

 September 8, 2013, 17:18 #24 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,209 Rep Power: 18 many thanks Bruno for consideration and investigation, when I wanted to run the "inviscid case" because its inviscid,so viscosity had to be 0(constant) then Code: `transport const;` had to be set and I wanted to set Code: `thermo janaf;// hConst;//` but it didn't accept it according to possible combinations: Code: ```Valid psiThermo types are: type mixture transport thermo equationOfState specie energy hePsiThermo homogeneousMixture const hConst perfectGas specie sensibleEnthalpy hePsiThermo homogeneousMixture sutherland hConst perfectGas specie sensibleEnthalpy hePsiThermo homogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy hePsiThermo inhomogeneousMixture const hConst perfectGas specie sensibleEnthalpy hePsiThermo inhomogeneousMixture sutherland hConst perfectGas specie sensibleEnthalpy hePsiThermo inhomogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy hePsiThermo multiComponentMixture const hConst perfectGas specie sensibleEnthalpy hePsiThermo multiComponentMixture sutherland janaf perfectGas specie sensibleEnthalpy hePsiThermo pureMixture const eConst perfectGas specie sensibleInternalEnergy hePsiThermo pureMixture const hConst perfectGas specie sensibleEnthalpy hePsiThermo pureMixture const hConst perfectGas specie sensibleInternalEnergy hePsiThermo pureMixture sutherland eConst perfectGas specie sensibleInternalEnergy hePsiThermo pureMixture sutherland hConst perfectGas specie sensibleEnthalpy hePsiThermo pureMixture sutherland hConst perfectGas specie sensibleInternalEnergy hePsiThermo pureMixture sutherland janaf perfectGas specie sensibleEnthalpy hePsiThermo pureMixture sutherland janaf perfectGas specie sensibleInternalEnergy hePsiThermo reactingMixture const hConst perfectGas specie sensibleEnthalpy hePsiThermo reactingMixture sutherland janaf perfectGas specie sensibleEnthalpy hePsiThermo singleStepReactingMixture sutherland janaf perfectGas specie sensibleEnthalpy hePsiThermo veryInhomogeneousMixture const hConst perfectGas specie sensibleEnthalpy hePsiThermo veryInhomogeneousMixture sutherland hConst perfectGas specie sensibleEnthalpy hePsiThermo veryInhomogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy heheuPsiThermo egrMixture const hConst perfectGas specie absoluteEnthalpy heheuPsiThermo egrMixture sutherland janaf perfectGas specie absoluteEnthalpy heheuPsiThermo egrMixture sutherland janaf perfectGas specie absoluteInternalEnergy heheuPsiThermo homogeneousMixture const hConst perfectGas specie absoluteEnthalpy heheuPsiThermo homogeneousMixture sutherland janaf perfectGas specie absoluteEnthalpy heheuPsiThermo homogeneousMixture sutherland janaf perfectGas specie absoluteInternalEnergy heheuPsiThermo inhomogeneousMixture const hConst perfectGas specie absoluteEnthalpy heheuPsiThermo inhomogeneousMixture sutherland janaf perfectGas specie absoluteEnthalpy heheuPsiThermo inhomogeneousMixture sutherland janaf perfectGas specie absoluteInternalEnergy heheuPsiThermo veryInhomogeneousMixture const hConst perfectGas specie absoluteEnthalpy heheuPsiThermo veryInhomogeneousMixture sutherland janaf perfectGas specie absoluteEnthalpy heheuPsiThermo veryInhomogeneousMixture sutherland janaf perfectGas specie absoluteInternalEnergy``` it verifies what you mentioned,it doesn't calculate entropy because its not Janaf and it can't be set for constant transport parameters in ordinary thermodynamical combinations. __________________ Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36 wyldckat OpenFOAM Bugs 18 October 21, 2010 05:51 hansel OpenFOAM Installation 62 March 19, 2010 15:43 Valle OpenFOAM Running, Solving & CFD 4 August 19, 2009 05:53 Cormac Main CFD Forum 4 August 13, 2000 23:26

All times are GMT -4. The time now is 13:29.