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

If you had a problem inputting your own peng-rob fluid

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

Like Tree1Likes
  • 1 Post By Anthony_impeller

Reply
 
LinkBack Thread Tools Display Modes
Old   October 6, 2011, 01:15
Default If you had a problem inputting your own peng-rob fluid
  #1
New Member
 
Anthony
Join Date: Dec 2010
Location: Newcastle, Australia
Posts: 2
Rep Power: 0
Anthony_impeller is on a distinguished road
Hello Forum-goers,

I've been doing turbomachinery CFD for the last year and have ran into so many problems that I want no one else to hit. So rather than asking a question I will post my solution.

The most recent one was specifying a "user-specified" peng robinson fluid.

My fluid of interest was R245fa and I had obtained my critical pressure, density (volume) and temperature, boiling temperature, accentric factor and molecular mass. This is easy to grab off NIST Refprop if you are looking, I think it is free.

So, the problem was in specifying the five coefficients for the fourth order polynomial of specific heat at constant pressure. After researching ANSYS help files and thermodynamic brain racking I found that is was in fact a zero pressure polynomial. I found the expected order of magnitude by looking at ANSYS' peng-rob data files for r134a.

To determine those specific to R245fa I plotted a bunch of points (C_p) from REFPROP for zero pressure and varying temperature in EXCEL. There was a lower and higher limit to where REFPROP was useful but the middle ground was enough to fit a polynomial to. When the fitted polynomial coefficients were input into CFX my solution would fail before any iteration. The error was:

"Error detected by routine MAKDAT" Blah blah blah.

That means I was making the specific heat capacity at constant volume C_v go negative.

The coefficient values obtained were much smaller than those ANSYS had for R134a. On second look I noted that I had not divided by the fluids specific gas constant (R with a bar in kJ/kgK). On rectifying this I had values on the same order of magnitude.

These values worked fine in ANSYS CFX and my solution solved like a charm (aka I had to work through convergence problems too but that is covered in FAQs and things).

So if you are having this problem now you can fix it.

Good Show,

Anthony.
Medison likes this.
Anthony_impeller is offline   Reply With Quote

Old   June 15, 2013, 06:57
Default R245fa
  #2
New Member
 
Gato Lopez
Join Date: Jun 2013
Location: Stuttgart
Posts: 8
Rep Power: 4
enuano is on a distinguished road
Hi!
I'm writting my master's thesis and have some problems when modeling an expantion for R245fa. Which kind of pure substance do you use? Peng Robinson WET refrigerant? because my problem is that pressure at outlet is ok, but temperature is too high and enthalpy and entropy make no sense. Also, do you still have the values for the coefficients you used in your simulation? THank you very much
Quote:
Originally Posted by Anthony_impeller View Post
Hello Forum-goers,

I've been doing turbomachinery CFD for the last year and have ran into so many problems that I want no one else to hit. So rather than asking a question I will post my solution.

The most recent one was specifying a "user-specified" peng robinson fluid.

My fluid of interest was R245fa and I had obtained my critical pressure, density (volume) and temperature, boiling temperature, accentric factor and molecular mass. This is easy to grab off NIST Refprop if you are looking, I think it is free.

So, the problem was in specifying the five coefficients for the fourth order polynomial of specific heat at constant pressure. After researching ANSYS help files and thermodynamic brain racking I found that is was in fact a zero pressure polynomial. I found the expected order of magnitude by looking at ANSYS' peng-rob data files for r134a.

To determine those specific to R245fa I plotted a bunch of points (C_p) from REFPROP for zero pressure and varying temperature in EXCEL. There was a lower and higher limit to where REFPROP was useful but the middle ground was enough to fit a polynomial to. When the fitted polynomial coefficients were input into CFX my solution would fail before any iteration. The error was:

"Error detected by routine MAKDAT" Blah blah blah.

That means I was making the specific heat capacity at constant volume C_v go negative.

The coefficient values obtained were much smaller than those ANSYS had for R134a. On second look I noted that I had not divided by the fluids specific gas constant (R with a bar in kJ/kgK). On rectifying this I had values on the same order of magnitude.

These values worked fine in ANSYS CFX and my solution solved like a charm (aka I had to work through convergence problems too but that is covered in FAQs and things).

So if you are having this problem now you can fix it.

Good Show,

Anthony.
enuano 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
Fluid Structure Interaction Apollo Main CFD Forum 5 July 4, 2011 16:15
VOF - fluid property problem weechristo FLUENT 1 April 11, 2009 15:08
How to apply negtive pressure to outlet bioman66 CFX 5 June 3, 2006 01:40
Intl Conf Computational Methods in Fluid Power Jacek Stecki Main CFD Forum 0 November 10, 2002 06:49
Simulation of Two Phase Fluid Flow Problem Using Fluent 5.0 Mohammad Al-Shannag Main CFD Forum 1 July 16, 1999 11:28


All times are GMT -4. The time now is 11:08.