Understanding Cubic Equations of State (Pressure-Volume Relationship)

  • Thread starter ben_ob1
  • Start date
  • #1
2
0
I am trying to teach myself some chemical thermodynamics, in particular some methods for fluid property evaluation. I am reading an excellent book which has helped a lot ("Introduction to Chemical Engineering Thermodynamics", Smith-Van Ness-Abbott, Ed. 7, ISBN: 007-124708-4).

One thing that the book described, which I am now having some problems with, is the pressure-volume relationship of the cubic equations of state. Figure 3.11 in this book (page 92) shows a nice graph of this relationship with 3 isotherm lines. I have tried to make my own version of this figure (hopefully it is attached ok):

attachment.php?attachmentid=37873&stc=1&d=1312914978.jpg


The diagram shows that at temperatures less than the critical temperature (e.g. the T3 line) and for low pressures there are 3 vapour solutions (intersections with a horizontal pressure line).

Here is the equation:
[itex]P=\frac{RT}{V-b}-\frac{a}{\left(V+\epsilon b\right)\left(V+\sigma b\right)}[/itex]

Where [itex]P[/itex] is pressure, [itex]T[/itex] is temperature and [itex]V[/itex] is molar volume; [itex]a[/itex] is a function of temperature (but is always positive) and [itex]b[/itex], [itex]\epsilon[/itex] and [itex]\sigma[/itex] are positive constants; [itex]R[/itex] is the universal gas constant.

Differentiating this with respect to Volume while keeping T constant (an isotherm line) gives:
[itex]\left(\frac{dP}{dV}\right)_{T}=-\frac{RT}{\left(b-V\right)^{2}}-\frac{a\left(\epsilon b + \sigma b + 2V\right)}{\left(V+\epsilon b\right)^{2}\left(V+\sigma b\right)^{2}}[/itex]

Did I do this right? If this differential is correct then I don't see how there can be any extrema (maxima or minima) for any isotherm lines. The 2 terms in the equation would have to cancel each other out in order for the differential to equal zero. As far as I can see both terms will always have the same sign (for positive values of V) meaning they can never cancel out. Is my logic wrong? How does the isotherm T3 get its local minimum and local maximum?

Hope you can help explain this for me.
Ben.
 

Attachments

  • CubicEoS_PV_Relationship.jpg
    CubicEoS_PV_Relationship.jpg
    12.6 KB · Views: 1,036

Answers and Replies

  • #2
2
0
Differentiating this with respect to Volume while keeping T constant (an isotherm line) gives:
[itex]\left(\frac{dP}{dV}\right)_{T}=-\frac{RT}{\left(b-V\right)^{2}}-\frac{a\left(\epsilon b + \sigma b + 2V\right)}{\left(V+\epsilon b\right)^{2}\left(V+\sigma b\right)^{2}}[/itex]

Oops found the mistake.

[itex]\left(\frac{dP}{dV}\right)_{T}=-\frac{RT}{\left(b-V\right)^{2}}+\frac{a\left(\epsilon b + \sigma b + 2V\right)}{\left(V+\epsilon b\right)^{2}\left(V+\sigma b\right)^{2}}[/itex]
 

Related Threads on Understanding Cubic Equations of State (Pressure-Volume Relationship)

Replies
1
Views
1K
Replies
4
Views
2K
Replies
13
Views
4K
  • Last Post
Replies
1
Views
2K
Replies
14
Views
5K
  • Last Post
Replies
0
Views
2K
  • Last Post
Replies
3
Views
946
Replies
5
Views
71K
Replies
3
Views
10K
Top