This preview shows pages 1–2. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: 2 2 ( ) 2 m m m m RT a f V P V b V bV b =+ =+Where, r c T T T = ; R= 83.14 cm3bar/(molK) (a) Create a MATLAB function file called Peng_Robinson.m that when returns the value of f, when sent the Volume (Vm), Pressure (P), Temperature, critical temperature (Tc), critical pressure (Pc), and . Hint: In the function file, the values of a, b, alpha, and Tr will be first computed so that f can be calculated. (25 points) (b) Create a MATLAB script file to plot f(V) versus V for Butane with proper annotations (25 points). Use the plot to estimate the root(s) that makes the function zero at P=100 bar and T= 500 K. Use the following fluid properties for Butane: w = 0.200 Tc = 425.1 K Pc = 37.96 bar 3. Determine the root of f(x)=x 2 ex using false position. Start with x L =0 and x U =1. Carry out 4 iterations. (25 points). PLEASE FILL OUT THE COURSE SURVEY AT h ttps://web.austin.utexas.edu/diia/oca/ . Your feedback is very important to me....
View Full
Document
 Spring '06
 Klaus

Click to edit the document details