Multi component VLE

Multi component VLE - Matlab Functions multicomponent,...

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

View Full Document Right Arrow Icon
Matlab Functions multicomponent, non-ideal VLE Several Matlab functions have been written for help with non-ideal multicomponent simulations. The following table lists the functions that have been written and tested to date: function Description Example bubptg Determine bubble point of non-ideal liquid mixture >> [t,ys]=bubptg(p,xs) buberg Compare predicted and experimental t & ys for isobaric data >> buberg(3,amipa) bubtpg Determine bubble pressure of non-ideal liquid mixture >> [p,ys]=bubtpg(t,xs) buberg2 Compare predicted and experimental p & ys for isothermal data >> buberg2(3,amipa2) marg1 Margules 2 suffix activity determination >> gams=marg1(xs) minvle Finds sum of absolute errors between predicted and experimental gammas >> minvle(Amarg) terf1 VLE ternary conversion program >> XG=terf1(amipa) The ternary sytem: acetone, methanol, isopropyl alcohol will be the example that we use to demonstrate the programs. Data for this system was reported by Freshwater and Pike in the J. Chemical Engineering Data 12 page 179 1967. Two tables of ternary data were presented for this system. In the first set, the pressure was fixed at 760 mmHg and in the second set the temperature was fixed at 55C. The data for each set stored in a text file that was then read into a Matlab session. The isobaric data was stored in a file called amipa . It looked like the following when listed in Matlab: x1 x2 y1 y2 T in C p in mmHg 0.0360 0.5460 0.1050 0.6750 69.3000 760.0000 0.0930 0.5230 0.2390 0.5810 66.8400 760.0000 0.1360 0.2980 0.3350 0.3690 68.8300 760.0000 0.0370 0.1120 0.1210 0.1699 76.9700 760.0000 0.8708 0.0761 0.8866 0.0912 56.3200 760.0000 0.5459 0.2693 0.6680 0.2569 58.4500 760.0000 0.3555 0.1917 0.5780 0.2096 63.6500 760.0000 0.0520 0.0287 0.1606 0.0485 77.4100 760.0000 0.1324 0.2635 0.3218 0.3391 69.2300 760.0000 0.0950 0.4453 0.2306 0.5135 68.1800 760.0000 0.0601 0.8993 0.1300 0.8556 63.1000 760.0000 0.1231 0.6965 0.2469 0.6776 63.2700 760.0000 0.3948 0.4841 0.5346 0.4237 58.4300 760.0000 0.2635 0.3920 0.4576 0.3861 62.9500 760.0000 0.3089 0.1895 0.5465 0.2040 64.6500 760.0000 0.0820 0.5602 0.1955 0.6259 66.7900 760.0000 0.3167 0.4201 0.4923 0.4023 61.2900 760.0000 0.3732 0.5910 0.4982 0.4923 57.5400 760.0000 0.4062 0.5355 0.5315 0.4490 57.5600 760.0000 0.4089 0.0351 0.6596 0.0552 64.7300 760.0000
Background image of page 1

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

View Full DocumentRight Arrow Icon
0.4854 0.1379 0.6768 0.1480 61.7700 760.0000 The program terf1 was used to convert this data into an array containg the liquid mols
Background image of page 2
Image of page 3
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 09/02/2011 for the course CHEM 102 taught by Professor Kmn during the Spring '11 term at Babson College.

Page1 / 9

Multi component VLE - Matlab Functions multicomponent,...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online