Mc Cabe Valentino

046f1006-22c4-473a-a942-f1efa5e3a38f 105 Mc Cabe Valentino Giovanni Valentino Español Este es un programa que calcula el No. de Platos de una mezcla binaria Giovanni Valentino 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

Protada

-------------------- Mc Cabe – Thiele --------------------

----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

PARA MEZCLAS BINARIA

Creado por Valentino S. Erick Giovanni

]]> Primeros Datos

-------------------- Mc Cabe – Thiele --------------------

PARA MEZCLAS BINARIA

Creado por Valentino S. Erick Giovanni

]]> DATOS NECESARIOS

Porcentaje de la de especie más volátil de la mezcla binaria en el Destilado "XD" x.D 99 % * :

Porcentaje de la especie más volátil de la mezcla binaria en los Fondos "XB" x.B 2 % * :

Porcentaje de las especie más volátil de la mezcla que entra a la columna "z" z 40 % * :

Flujo molar que ingresa a la columna "F" F 350 kmol hr / * :

Ingresa el valor de la calidad de alimentación "q" q 1.5 :

q 1 ≡ q 1.001 : q q : if

Ingresa el valor del Reflujo dado R 3.5 :

BALANCE DE MATERIA

F D B + ≡

F z * D x.D * B x.B * + ≡

B F F z * x.D / - 1 x.B x.D / - / :

D F B - :

Flujo Molar del Destilado D kmol hr / 137.1134

Flujo Molar en los Fondos B kmol hr / 212.8866

Propiedades

Nobre de los Componentes

]]>

Nombre del componente más volátil Componente.1 Benceno :

Nombre del componente menos volátil Componente.2 Tolueno :

Constantes de Antoine

]]>

Introduce las constantes de Antoine para cada uno de los componentes

que conforman la mezcla binaria. En mmHg y K

]]>

Constante "A" del componente más volátil A.EMV 15.9007 :

Constante "A" del menos volátil A.emv 16.0137 :

Constante "C" del componente más volátil C.EMV 52.36 - :

Constante "B" del menos volátil B.emv 3096.52 :

Constante "B" del componente más volátil B.EMV 2788.51 :

Constante "C" del menos volátil C.emv 53.67 - :

CA Dato A B C Componente.1 A.EMV B.EMV C.EMV Componente.2 A.emv B.emv C.emv 34 mat :

Tabla 1.1

]]> CA C.emv B.EMV B.emv C.EMV A.emv A.EMV D B R q F z x.B x.D CA C.emv B.EMV B.emv C.EMV A.emv A.EMV D B R q F z x.B x.D C.emv B.EMV B.emv C.EMV A.emv A.EMV D B R q F z x.B x.D

CA
Dato	A	B	C
Benceno	15.9007	2788.5100	-52.3600
Tolueno	16.0137	3096.5200	-53.6700
#
#
#
#
#

ELV

Presión del Sistema

]]>

Ingresa la presión en mmHg P 760 :

A A.EMV A.emv 21 mat :

B B.EMV B.emv 21 mat :

C C.EMV C.emv 21 mat :

Temperaturas de Sasaturación

]]>

T.K B A P ln - / C - vectorize :

Temp T.K K * :

Temperatura de Sasaturación 1 en K Temp 1 el 353.2551 K *

Temperatura de Saturación 1 en °C Temp 1 el °C 80.1051

]]>

Temperatura de Sasaturación 2 en K Temp 2 el 383.776 K *

Temperatura de Saturación 2 en °C Temp 2 el °C 110.626

Temperature T.K 2 el T.K 2 el T.K 2 el T.K 1 el - 10 / - T.K 2 el T.K 2 el T.K 1 el - 10 / 2 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 3 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 4 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 5 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 6 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 7 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 8 * (- T.K 2 el T.K 2 el T.K 1 el - 10 / 9 * (- T.K 1 el 111 mat :

T Temperature :

P°.1 e A.EMV B.EMV Temperature C.EMV + / - (^vectorize :

P°.2 e A.emv B.emv Temperature C.emv + / - (^vectorize :

X P P°.2 - P°.1 P°.2 - / vectorize :

Y X P°.1 * P / vectorize :

RSB Temperatra K P°1 mmHg P°2 mmHg X Y T.K 2 el P°.1 1 el P°.2 1 el X 1 el Y 1 el T.K 2 el T.K 2 el T.K 1 el - 10 / - P°.1 2 el P°.2 2 el X 2 el Y 2 el T.K 2 el T.K 2 el T.K 1 el - 10 / 2 * (- P°.1 3 el P°.2 3 el X 3 el Y 3 el T.K 2 el T.K 2 el T.K 1 el - 10 / 3 * (- P°.1 4 el P°.2 4 el X 4 el Y 4 el T.K 2 el T.K 2 el T.K 1 el - 10 / 4 * (- P°.1 5 el P°.2 5 el X 5 el Y 5 el T.K 2 el T.K 2 el T.K 1 el - 10 / 5 * (- P°.1 6 el P°.2 6 el X 6 el Y 6 el T.K 2 el T.K 2 el T.K 1 el - 10 / 6 * (- P°.1 7 el P°.2 7 el X 7 el Y 7 el T.K 2 el T.K 2 el T.K 1 el - 10 / 7 * (- P°.1 8 el P°.2 8 el X 8 el Y 8 el T.K 2 el T.K 2 el T.K 1 el - 10 / 8 * (- P°.1 9 el P°.2 9 el X 9 el Y 9 el T.K 2 el T.K 2 el T.K 1 el - 10 / 9 * (- P°.1 10 el P°.2 10 el X 10 el Y 10 el T.K 1 el P°.1 11 el P°.2 11 el X 11 el Y 11 el 125 mat :

Tabla 1.2

]]> RSB C.emv B.EMV B.emv C.EMV A.emv A.EMV D R q F z x.B x.D RSB Y X P°.2 P°.1 T Temperature Temp T.K C A C.emv B.EMV B.emv C.EMV A.emv A.EMV D B R q F z x.B x.D C.emv B.EMV B.emv C.EMV A.emv A.EMV D R q F z x.B x.D

RSB
Temperatra K	P°1 mmHg	P°2 mmHg	X	Y
383.7760	1784.2839	760.0000	0.0000	0.0000
380.7239	1650.0593	696.2997	0.0668	0.1450
377.6718	1523.6942	636.8873	0.1388	0.2783
374.6197	1404.8846	581.5571	0.2167	0.4006
371.5676	1293.3298	530.1077	0.3012	0.5126
368.5156	1188.7326	482.3431	0.3931	0.6148
365.4635	1090.7998	438.0716	0.4932	0.7079
362.4114	999.2422	397.1069	0.6027	0.7924
359.3593	913.7754	359.2679	0.7227	0.8689
356.3072	834.1196	324.3786	0.8546	0.9379
353.2551	760.0000	292.2685	1.0000	1.0000
#
#
#
#
#

Primeros Gráficos

XYPlot.XLimMin 1 el X 1 el :

XYPlot.XLimMax 1 el X X rows el :

XYPlot.YLimMin 1 el T.K 1 el :

XYPlot.YLimMax 1 el T.K 2 el :

XYPlot.XTick 1 el 0.1 :

XYPlot.YTick 1 el 5 :

X T augment Y T augment 2 1 sys

XYPlot.XLimMin 1 el X 1 el :

XYPlot.XLimMax 1 el X X rows el :

Plot p1x X P°.1 augment : p1y Y P°.1 augment : p2x X P°.2 augment : p2y Y P°.2 augment : 41 sys :

XYPlot.YLimMin 1 el P°.1 X rows el :

XYPlot.YLimMax 1 el P°.1 1 el :

XYPlot.XTick 1 el 0.1 :

XYPlot.YTick 1 el 100 :

Plot

XYPlot.XLimMin 1 el X 1 el :

XYPlot.XLimMax 1 el X X rows el :

Plot p1x X P°.1 augment : p1y Y P°.1 augment : p2x X P°.2 augment : p2y Y P°.2 augment : 41 sys :

XYPlot.YLimMin 1 el P°.2 X rows el :

XYPlot.YLimMax 1 el P°.2 1 el :

XYPlot.XTick 1 el 0.1 :

XYPlot.YTick 1 el 100 :

Plot

XYPlot.XLimMin 1 el X 1 el :

XYPlot.XLimMax 1 el X X rows el :

XYPlot.YLimMin 1 el Y 1 el :

XYPlot.YLimMax 1 el Y Y rows el :

XYPlot.XTick 1 el 0.1 :

XYPlot.YTick 1 el 0.1 :

X Y augment X X augment 2 1 sys Funciones

x LA q q 1 - / x * z q 1 - / - :

x LOSE x R 1 R + / * x.D 1 R + / + :

x LOSA x R z x.B - (* q x.D x.B - (* + q x.D x.B - (* R z x.B - (* + x.D - z + / * x.B z x.D - (* q x.D x.B - (* R z x.B - (* + x.D - z + / + :

x LO x x.D q 1 - (* z R 1 + (* + q R + / < x LOSA x LOSE if 11 line :

X 0 0.0668 0.1388 0.2167 0.3012 0.3931 0.4932 0.6027 0.7227 0.8546 1111 mat

x g X Y x cinterp eval :

x f x x.D q 1 - (* z R 1 + (* + q R + / < x R z x.B - (* q x.D x.B - (* + q x.D x.B - (* R z x.B - (* + x.D - z + / * x.B z x.D - (* q x.D x.B - (* R z x.B - (* + x.D - z + / + x R 1 R + / * x.D 1 R + / + if :

McCabe Codigo

McCabe.ER x.propuesta Clear x 1 el x.D : y 1 el x.D : x 2 el x.propuesta g y 1 el - x 1 el Broyden : y 2 el y 1 el : i 3 : x i 2 - el x.B > x i 1 - el x i 2 - el : y i 1 - el x i 1 - el f : x i el x.propuesta g y i 1 - el - x i 1 - el Broyden : y i el y i 1 - el : i i 2 + : 51 line while x i 1 - el x i 2 - el : y i 1 - el x i 2 - el : x y augment 101 line :

McCabe.ER McCabe.ER eval :

xer McCabe.ER 1 col :

yer McCabe.ER 2 col :

Etapas xer rows 2 / Floor 1 - :

Point Pinch

y.1 q q 1 - / x.1 * z q 1 - / - :

x.int 1 :

x.0 z :

y.0 x.0 :

x.1 1 z - (:

x.pinch y.0 y.1 y.0 - x.1 x.0 - / x.int x.0 - (* + (X Y x.int cinterp - x.int solve (: 0.5063

y.0 y.1 y.0 - x.1 x.0 - / xxx.int x.0 - (* + (X Y xxx.int cinterp - 1 Broyden 0.5063

y.pinch X Y x.pinch cinterp : 0.7188

y.LR 0 R 1 R + / * x.D 1 R + / + : 0.22

Y.LA 1 R z x.B - (* q x.D x.B - (* + q x.D x.B - (* R z x.B - (* + x.D - z + / * x.B z x.D - (* q x.D x.B - (* R z x.B - (* + x.D - z + / + : 1.2634

XYPlot.XLimMin 1 el X 1 el :

XYPlot.XLimMax 1 el X X rows el :

XYPlot.YLimMin 1 el Y 1 el :

XYPlot.YLimMax 1 el Y Y rows el :

XYPlot.XTick 1 el 0.1 :

XYPlot.YTick 1 el 0.1 :

Etapas 13

X Y augment X X augment x.D 0 2 1 mat x.D y.LR 2 1 mat augment x.B 1.0 2 1 mat x.B Y.LA 2 1 mat augment z x.pinch 2 1 mat z y.pinch 2 1 mat augment x.D x.D 2 1 mat x.D 0 2 1 mat augment x.B x.B 2 1 mat x.B 0 2 1 mat augment z z 2 1 mat z 0 2 1 mat augment xer yer augment 9 1 sys