Mc Cabe Valentino

046f1006-22c4-473a-a942-f1efa5e3a38f 108 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 97 % * :

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 / 140

Flujo Molar en los Fondos B kmol hr / 210

Propiedades

Nobre de los Componentes

]]>

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

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

Constantes de Antoine

]]>

Introduce las constantes de Antoine para cada uno de los componentes

que conforman la mezcla binaria. En mmHg y K

]]>

P e A B T C + / -^≡

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

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

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

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

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

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

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 C.EMV B.emv B.EMV A.emv A.EMV D B R q F z x.B x.D CA C.emv C.EMV B.emv B.EMV A.emv A.EMV D B R q F z x.B x.D C.emv C.EMV B.emv B.EMV A.emv A.EMV D B R q F z x.B x.D

CA
Dato	A	B	C
Benzene	15.9007	2788.5100	-52.3600
Water	18.3036	3816.4400	-46.1300
#
#
#
#
#

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 373.1521 K *

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

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 :

L X 1 el X 2 el X 3 el X 4 el X 5 el X 6 el X 7 el X 8 el X 9 el X 10 el X 11 el 111 mat : 0 0.092 0.184 0.278 0.373 0.47 0.57 0.672 0.778 0.887 1111 mat

V Y 1 el Y 2 el Y 3 el Y 4 el Y 5 el Y 6 el Y 7 el Y 8 el Y 9 el Y 10 el Y 11 el 111 mat : 0 0.154 0.293 0.419 0.531 0.633 0.724 0.805 0.878 0.942 1111 mat

Tabla 1.2

]]> RSB C.emv C.EMV B.emv B.EMV A.emv A.EMV D R q F z x.B x.D V L RSB Y X P°.2 P°.1 T Temperature Temp T.K C A C.emv C.EMV B.emv B.EMV A.emv A.EMV D B R q F z x.B x.D C.emv C.EMV B.emv B.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
373.1521	1350.3554	760.0000	0.0000	0.0000
371.1624	1279.0480	707.5996	0.0917	0.1543
369.1727	1210.6808	658.2326	0.1842	0.2934
367.1830	1145.1725	611.7611	0.2779	0.4188
365.1933	1082.4426	568.0516	0.3732	0.5315
363.2036	1022.4112	526.9746	0.4703	0.6327
361.2139	964.9994	488.4047	0.5699	0.7236
359.2242	910.1292	452.2207	0.6721	0.8049
357.2345	857.7233	418.3054	0.7776	0.8776
355.2448	807.7053	386.5455	0.8867	0.9424
353.2551	760.0000	356.8317	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 2.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 x.pinch y.pinch augment 3 1 sys 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 01 solve (: 0.481

y.pinch X Y x.pinch cinterp : 0.64309597913494

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

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.2557

x I z x.pinch 21 mat z y.pinch 21 mat x linterp :

x II x.B 121 mat x.B Y.LA 21 mat x linterp :

x.LO x I x II ≡ x solve (: 0.457

y.LO x.LO I : 0.57099999829618

y.Interpolado X Y x.LO linterp : 0.6188

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 G L V x cinterp :

McCabe Codigo

x.D 0.97

x.B 0.02

q 1.5

R 3.5

LO G McCabe.ER y.LO y.Interpolado ≥ No es posible trzar el gráfico. Cambia q ó diminuye la presión x.propuesta 0 : x 1 el x.D : y 1 el x.D : x 2 el x.propuesta G y 1 el - x.propuesta 01 solve : 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 LO : x i el x.propuesta G y i 1 - el - x.propuesta 01 solve : 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 if 11 line :

xer LO G McCabe.ER 1 col :

yer LO G McCabe.ER 2 col :

Etapas xer rows 2 / Floor 1 - :

Pues LO G McCabe.ER :

Titulo Coordenada X Coordenada Y 12 mat :

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 19

Nota:

Si el punto rojo sale de la curva del equilibrio líquido - vapor

no existe McCabe, debes cambiar "q" o disminuir la presión

]]> 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 x.LO y.LO augment xer yer augment 10 1 sys Titulo Pues Etapas yer xer McCabe.ER G LO LOSA LOSE LA II I Y.LA y.LR x.1 y.0 x.0 y.1 XYPlot.YTick XYPlot.XTick XYPlot.YLimMax XYPlot.YLimMin XYPlot.XLimMax XYPlot.XLimMin C.emv C.EMV B.emv B.EMV A.emv A.EMV D R q F z x.B x.D Titulo Pues Etapas yer xer McCabe.ER G LO LOSA LOSE LA II I Y.LA y.LR x.1 y.0 x.0 y.1 Plot p2y p2x p1y p1x XYPlot.YTick XYPlot.XTick XYPlot.YLimMax XYPlot.YLimMin XYPlot.XLimMax XYPlot.XLimMin V L Y X P°.2 P°.1 T Temperature Temp T.K C A C.emv C.EMV B.emv B.EMV A.emv A.EMV D B R q F z x.B x.D Etapas LOSA LOSE LA Y.LA y.LR x.1 y.0 x.0 y.1 C.emv C.EMV B.emv B.EMV A.emv A.EMV D R q F z x.B x.D

Pues
0.970	0.970
0.970	0.970
0.939	0.970
0.939	0.946
0.893	0.946
0.893	0.910
0.830	0.910
0.830	0.861
0.752	0.861
0.752	0.801
0.666	0.801
0.666	0.734
0.582	0.734
0.582	0.668
0.508	0.668
0.508	0.610
0.448	0.610
0.448	0.559
0.399	0.559
0.399	0.497
0.343	0.497
0.343	0.427
0.285	0.427
0.285	0.354
0.228	0.354
0.228	0.282
0.176	0.282
0.176	0.217
0.132	0.217
0.132	0.161
0.096	0.161
0.096	0.116
0.068	0.116
0.068	0.080
0.047	0.080
0.047	0.054
0.031	0.054
0.031	0.034
0.019	0.034
0.019	0.019
#
Coordenada X	Coordenada Y
#
#
#