# nasa polynomials

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

 March 9, 2011, 05:32 nasa polynomials #1 New Member   anton.gardi Join Date: Aug 2010 Posts: 21 Rep Power: 9 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.

 March 9, 2011, 19:45 #2 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 13,808 Rep Power: 107 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?

 March 10, 2011, 05:53 nasa polynomials #3 New Member   anton.gardi Join Date: Aug 2010 Posts: 21 Rep Power: 9 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

 March 10, 2011, 07:29 #4 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 13,808 Rep Power: 107 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.

 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 NewKid OpenFOAM 18 January 20, 2011 17:55 mystix OpenFOAM 4 August 20, 2010 07:21 Tomislav Sencic Main CFD Forum 8 December 15, 2009 11:09 Stone CFX 6 March 20, 2006 03:03 Stone CFX 2 February 7, 2006 11:17

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