CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

nasa polynomials

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 9, 2011, 05:32
Default nasa polynomials
  #1
New Member
 
anton.gardi
Join Date: Aug 2010
Posts: 21
Rep Power: 7
gardian is on a distinguished road
Hello

I m trying to simulate a helium 4 gas in a closed tube 10 m long.
the ends are warm like 120K while the rest of the tube is at 2K.
I wanted to simulate the gas as real so using the peng-robinson method
and the nasa polynomials for the pressure i have inside the tube.
in the CFX-help it is written that in this case all i want is the specific heat , the enthalpy and entropy.
i found them from a software and I have all the above values.
for CFX Cp/R= a1+ a2*T + a3*T^2/2 + a4*T^3/3 +a5*T^4/4
similar equations are given for Ho/R and So/R.

I made a fit for Cp,T and extracted the coefficients . likewise i found a6 and a7 from enthalpy and entropy.

In my case because of the Helium 4 lambda point i divided the interval of Temperatures in the lambda point because there is a discontinuity .

that is around 2K .

So I have the upper coefficients for the range 1.8-130K and the lower one for the range 1.2-1.8

but trying to run this (while I have done it in the past for Helium 3 with success ) the solver crashes with exit code error 255.

can you see something wrong in my steps i followed ??

thank you


** i also receive this error +--------------------------------------------------------------------+
| ERROR #001100279 has occurred in subroutine ErrAction. |
| Message: |
| c_fpx_handler: Floating point exception: Invalid Operand |
| |
|

i m using a coefficient for the fit that is like 1.2e-19 !!!!!!
is it too much for cfx?

Last edited by gardian; March 9, 2011 at 06:54.
gardian is offline   Reply With Quote

Old   March 9, 2011, 19:45
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,929
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
The coefficient of 1.2e-19 just means it is almost zero. Should be fine.

Almost certainly there is a problem with your EOS. Does this run OK with an ideal gas with similar properties?
ghorrocks is offline   Reply With Quote

Old   March 10, 2011, 05:53
Default nasa polynomials
  #3
New Member
 
anton.gardi
Join Date: Aug 2010
Posts: 21
Rep Power: 7
gardian is on a distinguished road
yes . i have also tried the zero's in the very small numbers and is fine.
also the ideal gas case runs properly and gives nice results .
the problem seems to be in the lower temperature interval that is 1.2 -1.8 K.
Is it correct for ansys to divide the intervals like that (1.2-1.8) and (1.8-130) ?
in this low interval the enthalply as i got it from HEPAK software for Helium 4 gives enthalpies that are very small.
PRESSURE TEMP DENSITY CP ENTROPY ENTHALPY QUALITY
[Pa] [K] [Kg/m3] [J/Kg-K] [J/Kg-K] [J/Kg] [-]
1200. 1.200 145.2 317.7 50.96 59.47 -2.000

is it ok ?
thank you
gardian is offline   Reply With Quote

Old   March 10, 2011, 07:29
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,929
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Am I correct in reading you have one EOS for above 1.8K and another for below? This type of discontinuity is always difficult to converge. Not only do you need to get the function value smooth over the 1.8K transition, but you also need to minimise the gradient discontinuity (Or preferably eliminate it) at the transition point. Often this means you need to "bend" your EOS a bit at the transition point so the gradients are also continuous.
ghorrocks is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Segmentation fault in running alternateSteadyReactingFoam,why? NewKid OpenFOAM 18 January 20, 2011 17:55
therm.dat NASA polynomials mystix OpenFOAM 4 August 20, 2010 07:21
NASA polynomials, thermo.dat tables Tomislav Sencic Main CFD Forum 8 December 15, 2009 11:09
How to Read Nasa Format Stone CFX 6 March 20, 2006 03:03
NASA Format Stone CFX 2 February 7, 2006 11:17


All times are GMT -4. The time now is 06:05.