×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

calculating enthalpy of mixture from Peng-Robinson EOS

calculating enthalpy of mixture from Peng-Robinson EOS

calculating enthalpy of mixture from Peng-Robinson EOS

(OP)
I know the equation for calculating the enthalpy of a mixture using Peng-Robinson.  However, there is a derivative, d(am(T))/dT,
where am(T) = xi*xj*(a(T,i)*a(T,j))^(1/2)
       a(T) = ac*(1+k*(1-(T/Tc)^(1/2))^2
         k  = 0.37464+1.54226*om - 0.26992*om^2
that I am having a difficult time solving.  I used my TI-89 and it gave me an extremely long result.  Does any one know the solution to this derivative or know where I may find it to verify my solution (which is too long to reproduce here)
Thanks!

RE: calculating enthalpy of mixture from Peng-Robinson EOS

Can you not make life easy for yourself by just using a numerical derivative? Just a thought, sshep

RE: calculating enthalpy of mixture from Peng-Robinson EOS

I haven't bothere to expand the equation, but it shouldn't be too hard to solve using a combination of chain and product rules .... to verify you have the correct solution, as sshep said, do it numerically. You could always use Simpson's Rule to solve the derivative. Set it up in excel, pretty simple to do!

Have a read of this
http://www.mactech.com/articles/mactech/Vol.09/09.10/SimpsonsRule/

or do a search on the net

RE: calculating enthalpy of mixture from Peng-Robinson EOS

(OP)
How do I use Simpson's rule to calculate a derivative?  I thought it was for integrals.

RE: calculating enthalpy of mixture from Peng-Robinson EOS

Sorry, you are right ... was having a blonde moment!!!! For a derivative, just do it from first principles using df/dx = lim (DeltaX ---> 0) (f(x+DeltaX)-f(x))/DeltaX ...

ut I still reckon you should solve the derivative! maybe two pages of work max?

RE: calculating enthalpy of mixture from Peng-Robinson EOS

Ok ... I must be bored today, or you are lazy ;)

d(am)/dT = (xi.xj/2)(a(T,i).d(a(T,j))/dT + a(T,j).d(a(T,i))/dT)^(-0.5)

Now (and I did this very quickly so it is probably wrong) ...
da/dT = ac(k^2/Tc-ac.(T.Tc)^(-0.5).(1-k^2))

Also where is the dependance on i or j in the equation for a?

RE: calculating enthalpy of mixture from Peng-Robinson EOS

(OP)
am(T) = (sumproduct of) xi*xj*(a(T,i)*a(T,j))^(1/2)
 a(T) = ac(i)*(1+k*(1-(T/Tc(i))^(1/2))^2
   k  = 0.37464+1.54226*om(i) - 0.26992*(om(i))^2

ac, om, and Tc are component specific, therefore a(T) is different for each component and am(T) is the sumproduct of the following equation, where i, j are different components:

am(T) = (sumproduct of) xi*xj*(a(T,i)*a(T,j))^(1/2)

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources