{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

# Example - T=380%K Tc=370%K Pc=42 bar R=8314 for i=1:50 if...

This preview shows pages 1–3. Sign up to view the full content.

Example No 13 Develop a Matlab code for predicting the state property of a fluid based on Peng-Robinson equation of state. You may refer to your note (pg 28-39) or any thermodynamics textbooks for documentation on Peng-Robinson equation of state. To test your Matlab code, determine the PV relation for propane at 360 K, 370 K and 380K. Plot the P-V diagram for propane for the above temperatures at a range of pressure between 30 – 80 bar Table 1. Data Material Propane Molecular Weight, MW 44.09 kg/kgmol Critical Pressure, Pc 42 bar Critical Temperature, Tc 370 K Acentric factor, ϖ 0.153 Solution Peng-Robinson Equation of State C r T T T = c r P P P = 2 26992 . 0 54226 . 1 37464 . 0 ϖ - ϖ + = κ ( 29 [ ] 2 2 1 r 2 r r T 1 1 T P 45724 . 0 A - κ + = r r T P 07780 . 0 B = ( 29 ( 29 ( 29 0 AB B B Z B 2 B 3 A Z 1 B Z 2 3 2 2 3 = - + + - - + - + specific volume P ZRT = ν Density ν = ρ MW

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Matlab Code: % Equation of State: Peng-Robinson % M=44.09; %kg/kgmol w=0.153; %P=120; % bar T=380;