Introductory Chemical Engineering Thermodynamics By J.R. Elliott - - PowerPoint PPT Presentation

Introductory Chemical Engineering Thermodynamics

By J.R. Elliott and C.T. Lira

Chapter 6 - Engineering Equations of State

Chapter 6 - Engineering Equations of State Slide 1

• II. Generalized Fluid Properties

The principle of two-parameter corresponding states

• 10

10 20 30 40 50 0.1 0.2 0.3 0.4 Density (g/cc) T=133K T=191K T=286K

ρL ρV

• 10

10 20 30 40 50 0.1 0.2 0.3 0.4 0.5 Density (g/cc) T=329K T=470 T=705K

VdW Pressure (bars) in Methane VdW Pressure (bars) in Pentane

Critical Definitions: Tc - critical temperature - the temperature above which no liquid can exist. Pc - critical pressure - the pressure above which no vapor can exist. ω - acentric factor - a third parameter which helps to specify the vapor pressure curve which, in turn, affects the rest of the thermodynamic variables. Note: at the critical point,

c c T T

P and T P P at and

2 2

=         =         ∂ρ ∂ ∂ρ ∂

Chapter 6 - Engineering Equations of State Slide 2

• II. Generalized Fluid Properties

The van der Waals (1873) Equation Of State (vdW-EOS)

Based on some semi-empirical reasoning about the ways that temperature and density affect the pressure, van der Waals (1873) developed the equation below, which he considered to be fairly crude. We will discuss the reasoning at the end of the chapter, but it is useful to see what the equations are and how we use them before deriving the details. The vdW-EOS is:

( ) ( )

Z b b a RT b a RT = + − − = − − 1 1 1 1 ρ ρ ρ ρ ρ van der Waals’ trick for characterizing the difference between subcritical and supercritical fluids was to recognize that, at the critical point,

c c T T

, P T P P at and

2 2

=         =         ∂ρ ∂ ∂ρ ∂

Since there are only two “undetermined parameters” in his EOS (a and b), he has reduced the problem to one of two equations and two unknowns. Running the calculus gives: a = 0.475 R2Tc

2/Pc; b = 0.125 RTc/Pc

The capability of this simple approach to represent all of the properties and processes that we will discuss below is a tribute to the genius of van der Waals.

Chapter 6 - Engineering Equations of State Slide 3

• II. Generalized Fluid Properties

The principle of three-parameter corresponding states Reduced vapor pressure behavior:

1 / T r 1 1 .2 1 .4 1 .6 1 .8 2 2 .2 2 .4

To improve our accuracy over the VdW EOS, we can generate a different set of PVT curves for each family of compounds. We specify the family of compounds via the "third parameter" i.e. ω. Note: The specification of Tc , pc , and ω provides two points

• n the vapor pressure curve. The key to accurate characterization of the vapor liquid

behavior of mixtures of fluids is the accurate characterization of the vapor pressure of pure fluids. VLE was central to the development of distillation technology for the petrochemical industry and provided the basis for most of today’s process simulation technology.

Chapter 6 - Engineering Equations of State Slide 4

• II. Generalized Fluid Properties

The Peng-Robinson (1976) Equation of State

( ) ( )

2 2 2 2 2

2 1 1 1 = Z

• r

2 1 1 ρ ρ ρ ρ ρ ρ ρ ρ b b b bRT a

• b

b b a b RT p − + ⋅ − − + − − =

where ρ = molar density = n/V By fitting the critical point, where ∂P/∂V = 0 and ∂2P/∂V2 = 0, a = 0.457235528 αR2Tc

2/Pc; b = 0.0777960739 RTc/Pc

α = [1+ κ (1-√Tr)]2 ; κ = 0.37464 + 1.54226ω - 0.26992ω2 ω ≡ -1 - log10(Psat/Pc)Tr =0.7 ≡ “acentric factor” Tc , Pc , and ω are reducing constants according to the principle of corresponding states. By applying Maxwell's relations, we can calculate the rest of the thermodynamic properties (H,U,S) based on this single equation.

Chapter 6 - Engineering Equations of State Slide 5

• II. Generalized Fluid Properties

Solving the Equation of State for Z

( )

Z = 1 1- / Z B A B B Z B Z B Z − ⋅ + − / / ( / ) 1 2

2

bρ ≡ B/Z where B ≡ bP/RT Z ≡ P/ρRT A ≡ aP/R

2T 2

Rearranging yields a cubic function in Z Z3 -(1-B)Z2 + (A-3B2-2B)Z - (AB-B2-B3) = 0 Naming this function F(Z), we can plot F(Z) vs. Z to gain some understanding about its roots

Chapter 6 - Engineering Equations of State Slide 6

• II. Generalized Fluid Properties

Solving the Equation of State for Z

• 0.04
• 0.02

0.02 0.04 0.06 0.08 0.1 0.5 1 P= Pvap T< Tc F Z

• 0.04
• 0.02

0.02 0.04 0.06 0.08 0.1 0.5 1 P> > Pvap T< Tc F Z

• 0.04
• 0.02

0.02 0.04 0.06 0.08 0.1 0.5 1 P< < Pvap T< Tc F Z

• 0.06
• 0.04
• 0.02

0.02 0.04 0.06 0.08 0.1 0.5 1 T> Tc F Z

Chapter 6 - Engineering Equations of State Slide 7

• II. Generalized Fluid Properties

Solving the Equation of State for Z (cont.)

• 1. Guess Zold=1 or Zold=0 and compute Fold(Zold).
• 2. Compute Fnow(Z) at Z=1.0001 or Z=0.0001
• 3. "Interpolate" between these guesses to estimate where F(Z)=0.

∆Z = (0 - Fnow)*(Znow-Zold)/(Fnow - Fold).

• 4. Set Fold=Fnow, Zold=Znow, Znow=Zold+∆Z
• 5. If |∆Z/Znow| < 1.E-5, print the value of Znow and stop.
• 6. Compute Fnow(Znow) and return to step 3 until step 5 terminates.

Chapter 6 - Engineering Equations of State Slide 8

• II. Generalized Fluid Properties

An Introduction to the Radial Distribution Function

Nc g r

• R o

=

∫

ρ π ( ) 4 r dr;

2

where g(r) is our "weighting factor" henceforth referred to as the radial distribution function (rdf).

g r / 1 1 2

The body centered cubic unit cell

The radial distribution function for the bcc hard sphere fluid

Chapter 6 - Engineering Equations of State Slide 9 An Introduction to the Radial Distribution Function

dr r r g Nc

• R
• 2

4 ) ( π ρ∫ =

g r / 1

1 2 3 4 1 2 3 4 r/σ g

The low density hard-sphere fluid The high density hard sphere fluid

g r / 1 1 2 3 4 1 2 3 4 r/σ g

η=0

βε =1.0

η=0.4

The low density square-well fluid The high density square-well fluid

Chapter 6 - Engineering Equations of State Slide 10

• II. Generalized Fluid Properties

The connection from the molecular scale to the macroscopic scale The Energy Equation U U nRT N RT N u g 4

id A A

− =

∞

∫

ρ π 2 r dr

2

The Pressure Equation P RT N RT rN du dr

A A

ρ ρ π = −      

∞

∫

1 6 g 4 r dr

2

P RT N RT rN du dr N rN du dr g 4

A A A A

/ ρ ρ π ρ π

σ σ

= −            

∫ ∫

∞

1 6 g 4 r dr - 6RT r dr

2 2

Chapter 6 - Engineering Equations of State Slide 11

• II. Generalized Fluid Properties

The Van der Waals Equation of State

( )

−       ≈

∫

N RT rN du dr

A O A

ρ π ρ ρ

σ

6 g 4 r dr b 1-b

2

where b~close-packed volume N RT rN du dr g 4 N N RT x du dx g 4

A A A A

ρ π ρ σ ε π σ

σ

6 6

1 ∞ ∞

∫ ∫

            ≡ r dr = x dx where x r /

2 3 2

a N N x du dx g 4

A A

≡      

∞

∫

σ ε ε π

σ 3

6 / x dx

2

The resulting equation of state is:

( ) ( )

Z b b a RT b a RT = + − − = − − 1 1 1 1 ρ ρ ρ ρ ρ By fitting the critical point, where ∂P/∂V = 0 and ∂2P/∂V2 = 0, a = 0.475 R2Tc

2/Pc; b = 0.125 RTc/Pc

Chapter 6 - Engineering Equations of State Slide 12

• II. Generalized Fluid Properties
• Example. From the molecular scale to the continuum

Suppose that the radial distribution function can be reasonably represented by:

g b kT x ≈ + 1 1

4

ρε

where x ≡ r/σ and b ≡ π/6 NAσ3 Derive expressions for the compressibility factors of fluids that can be accurately represented by the square-well potential. Evaluate this expression at bρ = 0.2 and ε/kT = 1. Solution: First read Appendix C then note that, for the square-well potential, exp(-u/kT) = exp(ε/kT) H(r-σ) + [exp(ε/kT)-1] [1-H(r-1.5σ)] Taking the derivative of the Heaviside function gives the Dirac delta in two places: ( ) ( )

[ ] [ ]

∫

− − − − + = dr r kT r y r kT r y r N Z

A 3

4 1 / exp ) ( ) 5 . 1 ( / exp ) ( ) ( 6 1 π ε σ δ ε σ δ ρ

( ) ( ) [ ]

{ }

1 / exp ) 5 . 1 ( 5 . 1 / exp ) ( 6 4 1

3 3

− − + = kT y kT y N Z

A

ε σ ε σ ρσ π

Noting that y(r) ≡ g(r) exp(u/kT) and that exp(u/kT) is best evaluated inside the well: Z = 1+ 4bρ { g(σ+) - 1.53[1-exp(-ε/kT)] g(1.5σ-)} This is true for SW with any g(r). For the above expression: g(σ+) = 1 + bρ ε/kT and g(1.5σ-) = 1 + 0.198 bρ ε/kT Z = 1+ 4bρ {1 + bρ ε/kT - 1.53[1-exp(-ε/kT)]( 1 + 0.198 bρ ε/kT)} Z = 1+4(0.2){1+0.2-2.1333*1.0396} = 0.1858