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 cm3-bar/(mol-K) (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- e-x 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
- Numerical Analysis, Binary numeral system, meaningful variable names, relative error respectively.