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

Peng-Robinson roots problems for H2 at normal temperature

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 31, 2022, 09:56
Default Peng-Robinson roots problems for H2 at normal temperature
  #1
New Member
 
Luc Cena
Join Date: Oct 2022
Posts: 1
Rep Power: 0
cenaluc is on a distinguished road
Hi All

I have problems to work with Peng Robinson EoS for temperatures way above the Tcr of H2.

I am using the cardano method base to find the roots with the hydrogen properties: Tc=33K ; Pc=1.32MPa ; omega=-0.22
Target is to find the molar volume for normal pressure temperature range, like T1=20°C / P1=10MPa, so temperature and pressure are way above the critical values.

Well, I was expecting only one root because both P/T are above the critical values, but intermediate parameters analysis ( b^2 / 4 ) + ( m^3 / 27 ) is negative so I should expect 3 real roots with the one in the middle meaningless as usual. I can solve the cubic with the following solution:
Z³ - 0.933 * Z² - 0.105 * Z + 0.002 = 0
Z1= -0.117
Z2= +1.034
Z3= +0.017

All three are correct, so it looks ok, three solutions real, not complex as expected. But here comes the problems:
First solution is negative (!!!) : Theoretically it should provide the bubble point with but Vm= Z1 * R * T / p is negative if I use a negative Z, it does not make sense.
Second solution is above 1 : Theoretically all solutions should be in the range [0;1] if I understand the theory well.
Third solution is the middle one : I should not consider it.

Can you tell me why this math behaviour for H2? Am I too far from the critical point for this gas and I cannot use PR EOS or am I doing anything wrong?

I am obtaining the attached by plotting the data but I do not think that going negative with the volume is somethig good lol. Any suggestion??
Thanks
Attached Images
File Type: jpg plot.jpg (144.8 KB, 3 views)
cenaluc is offline   Reply With Quote

Reply

Tags
hydrogen, peng-robinson eoe, real fluids equation

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45
Problems referencing to temperature andreas CFX 7 October 30, 2006 08:26
Problems with temperature jump at the boundary tsjb00 OpenFOAM Running, Solving & CFD 4 August 25, 2005 14:33
Temperature calculation problems Julie Polyakh Siemens 1 January 15, 2002 21:50


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