Mini Project 1: Matlab Program for Multicomponent Flash Distillation with Rachford-Rice Equation Problem Statement A feed that is 50 mol% methane, 10 mol% n-butane, 15 mol% n-pentane, and 25 mol% n-hexane is flash distilled. F = 150 kmol/h. Drum pressure = 250 kPa, drum temperature = 10. Figure 1 shows the diagram illustration of the flash distillation. Use the DePriester chart (Figure 2) to find K values and Rachford-Rice equation to find V/F, Xi, yi, V, L. Write a Matlab script for solving the Rachford- Rice equation, using V/F=0.1 as the initial guess. V, y F = 150 kmol/h Tdrum = 10C Parum = 250kPa L, X PRESSURE, PO 5000 4000 3500 METHANE METHANE ETHYLENE ETHANE PROPYLENE PROPANE ETHYLENE ETHANE PROPYLENE PROPANE Figure 2. DePriester Chart. 68 ISOBUTANE N-BUTANE ISOPENTANE N-PENTANE but TUOTTOTETTOREN Step UTIN apg ISOBUTANE N-BUTANE ISOPENTANE N-PENTANE N-HEXANE N-HEPTANE N-OCTANE N-NONANE N-DECANE N-HEXANE . N-HEPTANE N-OCTANE ww 88 N- NONANE N-DECANE "" TEMPERATURE, C Information About the Matlab Script A Matlab function is expected to be submitted. The inputs argument to the function include feed flow rate F, mol fractions of each species in the feed steam z and discovered K values from De Priester chart, and initial guess of V/F. The output arguments include the computed V/F, Xi, yi, V and L. Inside the function, Newton's iteration will be used to find the root of the Rachford-Rice equation. Please use a loop (with a stopping criterion) to conduct the Newton's iteration, instead of doing the iteration manually one by one. The components in the source code can be referred to Table 1. Results to be Submitted in the Report and as Attachments (1) Demonstrate the discovered K values from the DePriester chart. (2) Paste your source code of the function to the report. (3) Show the computed final V/F, Xi, yi, V, L values. Plot the function value f[(V/F)k] over iterations. (4) Attach your.m or .py" file as a separate file during submission so we can run your scripts. Table 1. Grading Rubric for the Project Points Note 10 Total: 70 10 z and K are vectors . 10 . 10 Items (1) K values (2) Souce code for Rachford-Rice equation Input and output arguments Stopping criterion for the loop for Newton's method Calculating the function value at (V/F) Calculating the function derivative at (V/F)k Update the value of V/F Calculating Xi, yi Calculating V, and L (3) Correctness of final V/F, Xi, yi, V, L and the plot of f[(V/F)k] over iterations 10 . 10 10 10 20 Total: 100