CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Newer Janaf Model (http://www.cfd-online.com/Forums/openfoam/76123-newer-janaf-model.html)

iyer.arvind May 14, 2010 13:50

Newer Janaf Model
 
1 Attachment(s)
The older janaf thermodynamic model is implemented in OpenFOAM. For my use i have implemented the newer Janaf model, which models Cp0 H0 S0 for most species from 200K to 20000K.

The current code is for OpenFOAM-1.5 ( :( ) but must not be too difficult to port to 1.6.x. (I will do that if i get some time for it). However if someone wants to use it, i am attaching the patch for thermophysicalModels.

I would like to know if people find it useful.

NOTE This code is hardly tested! I coded it just yesterday, and it seemed to work fine. Please test it for your purpose, and kindly let me know of any bugs.

NOTE 2: use bzip2 to unzip the patch file

iyer.arvind May 18, 2010 01:03

any problems?
 
If people are finding any difficulty in applying the patch, I would be glad to mention the procedure to do that. 70 views and no replies is slightly discouraging ! :(

Chris Lucas May 18, 2010 04:33

Hi,

Thank you for posting your work.

I tried to run the script but it didn't work. I'm using OpenFoam1.6.x.
Would it be possible to post the main files of the new Janaf Tables?

Regards
Christian

iyer.arvind May 18, 2010 07:00

Arvind Iyer
 
Thanks for reponding. The patch is for OpenFOAM-1.5 (sorry to say) but I guess if i could give you the files you should ot find it difficult in implementing in OpenFOAM-1.6. Also i will try to spare some time this weekend and try and get it up for 1.6. However, I am unable to upload all the files here because they are beyond 100kb. Could you please let me know a mode of transfer? I would be glad to share it with you.

Chris Lucas May 19, 2010 07:04

Openfoam Wiki
 
Hi,

you could upload the new Janaf Tables onto the Openfoam Wiki homepage. That way, everybody who wants to use them can download them.

Regards,
Christian

iyer.arvind May 19, 2010 07:19

I think it would be better if i could get the 1.6.x version, and push to the OpenFOAM-extend project. However i will see if i can put the files in the wiki for now.

CedricVH May 19, 2010 10:44

Thanks for your efforts Iyer! I have never used any Janaf model, but I have read on the forums that there are some problems with the current Janaf model in OpenFOAM so your contribution is appreciated!

chandramurthy May 19, 2010 11:09

good work arvind...when i go through this file, i found some things,
1) I guess it is trying to replace existing Janaf, can both(janaf-I and janf-II) simultaneously work?
2) any validation being done?

iyer.arvind May 20, 2010 00:14

The problem is that janaf (reader especially) is quite deeply coded. the reader appears in the same file which reads the chemkin data.

_____/thermophysicalModels/combustion/chemistryReaders/chemkinReader/chemkinLexer.L

Also,

_____/thermophysicalModels/specie/reaction/reactions/makeReactionJanaf.H:
typedef sutherlandTransport<specieThermo<janafThermo<perfe ctGas> > > reactionThermo;

As you can see janaf is quite deep embedded in the code. If it were abstracted to a generic thermophysical model interface (abstract class), and janaf being a child class with its own reader, it would have been much better, and easily we could have either models to work with.

universez May 28, 2010 19:05

Great work.
 
Good work !

thekay May 29, 2010 11:56

Thank you Iyer, good piece of work. I will probably need to test it in a couple of weeks. Will let you know of my results.

skarnani July 16, 2010 01:09

Hello All.
I was curious about this new Janaf model. I'm currently running OpenFOAM 1.6.x. It would be an interesting exercise to try and port it over but, in all honesty, I wouldn't really be sure where to begin. To start, how would I go about applying a patch?

Thanks.
Sunny

iyer.arvind July 16, 2010 07:34

Patch is for 1.5
 
Sorry to say but the patch is only for 1.5. I am not getting the breathing space to implement for 1.7. I cannot comment on how easy it will be with 1.6x. If you are unlucky, you can at best take inspiration from it :);

chbenz October 6, 2010 07:12

Can you give us the source code files for the new janaf model. I'm not able to extract them from the patch-file.

regards,
christof

bfisch90 April 2, 2011 06:10

Thank you for your work. How does the thermo.dat file look like? I tried to implement your new janaf model to OF-1.7.1 an would like to test it. Maybe you can post it. If it works I would post the updated files here.

Thomas

iyer.arvind April 2, 2011 06:42

Link
 
http://www.galcit.caltech.edu/EDL/pu...rmo/thermo.inp

iyer.arvind April 2, 2011 06:43

If you require the source code, please mail me, i will send you the source code too

bfisch90 April 14, 2011 21:36

Updated new janaf model
 
Hallo,

I have implemented Iyers newer janaf model to OpenFOAM 1.7.1. You can find and download it in the OpenFOAM wiki:
http://openfoamwiki.net/index.php/Janaf.
There you can also find a modificated test case, which should run if the installation was correct.

Thomas

yingkun June 24, 2011 04:32

Hi,Thomas
I try as you say,then Allwmake,and run my case,it seems doesn't work, it shows on the terminal:
--> FOAM FATAL ERROR:
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = 5014.46
why the temperature range doesn't change?
may be there're something I missing?

bfisch90 June 24, 2011 10:47

Quote:

Originally Posted by yingkun (Post 313370)
Hi,Thomas
I try as you say,then Allwmake,and run my case,it seems doesn't work, it shows on the terminal:
--> FOAM FATAL ERROR:
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = 5014.46
why the temperature range doesn't change?
may be there're something I missing?

Hallo Yingkun,

there must be at least one species in the janaf table, which has only a temperature range between 0 and 5000 K. A quick and dirty solution would be to find these species in the table and set the temperature limit from 5000 K to 6000 K and look what happens. Maybe (or better hopefully) these species have such a low concentration, that it does not have a negative effect on your calculation.

Hope this helps!
Thomas


All times are GMT -4. The time now is 23:48.