Separación

Basadas en 2024

el Equilibrio

ECUACIONES MESH

Diego Octavio Perez Molina

20208250
6to semestre
Ingeniería Química
 UNIVERSIDAD AUTONOMA DE TLAXCALA
FACULTYAD DE CIENCIAS BASICAS INGENIERIA Y TECNOLOGIA
LIC. INGENIERIA QUIMICA

ECUACIONES MESH.
1. Condiciones térmicas del caso de estudio
ComponentesZi F lbmol/hr Pf psia Etapas (N)
propano 0.3 30 100 R 5
n-butano 0.3 30 P mmhg 1.5 N alimentacion
N-Pentano 0.4 40 5171.6 3
1 100

2. Ecuaciones para usar.

Balance de Masa, Equilibrio

Formulas y Restricciones
Masa
1. 𝐹𝑗 𝑧𝑖𝑗 𝐿𝑗 −1 𝑋𝑖,𝑗 −1 + 𝑉𝑗 +1 𝑦𝑖𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )𝑋𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )𝑌𝑖𝑗 = 0
Equilibrio
𝑦 𝑖𝑗
2. 𝐾𝑖𝑗 − =0
𝑋 𝑖𝑗

Formulas Sumatoria
Masa 𝑐
3. 𝑋𝑖𝑗 − 1 = 0
1. 𝐹𝑗 𝑧𝑖𝑗 𝐿𝑗 −1 𝑋𝑖,𝑗 −1 + 𝑉𝑗 +1 𝑦𝑖𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )𝑋𝑖𝑗 − (𝑉𝑗 + 𝑤𝑖=1
𝑗 )𝑌𝑖𝑗 = 0
𝑐
𝑖=1 𝑌𝑖𝑗 −1 =0
Equilibrio 4.
𝑦 𝑖𝑗
2. 𝐾𝑖𝑗 − =0
𝑋 𝑖𝑗 H Entalpia
Sumatoria
5. 𝐹𝑗 𝑧𝑖𝑗𝑐 𝐿𝑗 −1 ℎ𝐿,𝑗 −1 + 𝑉𝑗 +1 ℎ𝑣𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )ℎ𝐿𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )ℎ𝑣𝑗 − 𝑄𝑗 = 0
Balance de energía 3. 𝑖=1 𝑋𝑖𝑗 − 1 = 0
𝑐 𝐾𝑖𝑗 = (𝐾𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗 , 𝑦𝑖𝑗 )
4. 𝑖=1 𝑌𝑖𝑗 −1 =0
H Entalpia ℎ𝑖𝑗 = ℎ𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗
5. 𝐹𝑗 𝑧𝑖𝑗 𝐿𝑗 −1 ℎ𝐿,𝑗 −1 + 𝑉𝑗 +1 ℎ𝑣𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )ℎ𝐿𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )ℎ𝑣𝑗 − 𝑄𝑗 = 0
ℎ𝑣𝑗 = (ℎ𝑣𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑦𝑖𝑗 )
𝐾𝑖𝑗 = (𝐾𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗 , 𝑦𝑖𝑗 )
Especificaciones
ℎ𝑖𝑗 = ℎ𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗
Formulas 6. 𝐹𝑗 , 𝑧𝑖𝑗 , 𝑇𝐹𝑗 , 𝑃𝐹 , 𝑃𝑗 , 𝑈𝑗 , 𝑊𝑗 , 𝑄𝑗
ℎ𝑣𝑗 = (ℎ𝑣𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑦𝑖𝑗 )
Masa
Especificaciones Ecuaciones
1. 𝐹𝑗 𝑧𝑖𝑗 𝐿𝑗 −1 𝑋𝑖,𝑗 −1 + 𝑉𝑗 +1 𝑦𝑖𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )𝑋𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )𝑌𝑖𝑗 = 0
𝑇𝐹𝑗 , 𝑃𝐹 , 𝑃𝑗 , 𝑈𝑗 , 𝑊𝑗 , 7.
6. 𝐹𝑗 , 𝑧𝑖𝑗 ,Equilibrio 𝑄𝑗 2𝐶 + 3 𝑁(2𝐶 + 3)
2.Ecuaciones
𝑦 𝑖𝑗
𝐾𝑖𝑗 −
𝑋
=0
𝑖𝑗
Incógnitas
7. 2𝐶 + Sumatoria
3 𝑁(2𝐶 + 3)
8. 𝑋𝑖,𝑗 , 𝑦𝑖,𝑗 , 𝐿,𝑗 , 𝑉,𝑗 𝑇,𝑗
3. Incógnitas
𝑐
𝑖=1 𝑋𝑖𝑗 − 1 = 0
8. 𝑋𝑐𝑖,𝑗 , 𝑦𝑖,𝑗 , 𝐿,𝑗 , 𝑉,𝑗 𝑇,𝑗
Matriz tridiagonal
4. 𝑖=1 𝑌𝑖𝑗 − 1 = 0
Matriz tridiagonal 9. 𝐴𝑖 𝑋𝑗 −1 𝐵,𝑗 𝑋𝑖𝑗 + 𝐶𝑗 𝑋𝑖𝑗 +1 = 𝐷
H Entalpia

5. 𝐹𝑗 𝑧𝑖𝑗 𝐿𝑗 −1 ℎ𝐿,𝑗 −1 + 9.
𝑉𝑗 +1𝐴ℎ𝑖 𝑋𝑗 −1 𝐵,𝑗 𝑋𝑖𝑗 + 𝐶𝑗 𝑋𝑖𝑗 +1 = 𝐷
𝑣𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )ℎ𝐿𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )ℎ𝑣𝑗 − 𝑄𝑗 = 0
Balance de Masa, Equilibrio
Formulas y Restricciones
𝐾𝑖𝑗 = (𝐾𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗 ,10.
𝑦𝑖𝑗 𝐴
) 𝐽
𝐽 = 𝑉𝑗 + (𝐹 − 𝑊𝑚 − 𝑈𝑚 )𝑉1
𝐽 Masa ,𝑚
𝑀=1
10. 𝐴𝐽 = 𝑉𝑗 + 𝑀=1( 𝐹,𝑚 − 𝑊𝑚 − 𝑈𝑚 )𝑉1
ℎ𝑖𝑗 = ℎ 𝑇,𝑃,𝑋
1. 𝑖𝑗 𝐹 𝑗𝑧 𝑗𝐿 𝑖𝑗 𝑋 + 𝑉𝑗 +1 𝑦𝑖𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )𝑋𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )𝑌𝑖𝑗 = 0
𝑗 𝑖𝑗 𝑗 −1 𝑖,𝑗 −1
ℎ𝑣𝑗 = (ℎ𝑣𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑦𝑖𝑗 ) Equilibrio
Especificaciones 2. 𝐾𝑖𝑗 −
𝑦 𝑖𝑗
=0
𝑋 𝑖𝑗
6. 𝐹𝑗 , 𝑧𝑖𝑗 , 𝑇𝐹𝑗 , 𝑃𝐹 , 𝑃𝑗 , 𝑈𝑗 , 𝑊𝑗 , 𝑄𝑗
Sumatoria
Ecuaciones 𝑐
3. 𝑖=1 𝑋𝑖𝑗 − 1 =0
7. 2𝐶 + 3 𝑁(2𝐶 + 3)
𝑐
Incógnitas
4. 𝑖=1 𝑌𝑖𝑗 −1 =0
8. 𝑋𝑖,𝑗 , 𝑦𝑖,𝑗 , 𝐿,𝑗 , 𝑉,𝑗 𝑇,𝑗
H Entalpia

5. 𝐹𝑗 𝑧tridiagonal
Matriz 𝑖𝑗 𝐿𝑗 −1 ℎ𝐿,𝑗 −1 + 𝑉𝑗 +1 ℎ𝑣𝑗 +1 − (𝐿𝑗 + 𝑈𝑗 )ℎ𝐿𝑖𝑗 − (𝑉𝑗 + 𝑤𝑗 )ℎ𝑣𝑗 − 𝑄𝑗 = 0

𝐾𝑖𝑗
9. 𝐴𝑖 𝑋𝑗 −1 𝐵,𝑗 𝑋𝑖𝑗 + 𝐶𝑗 𝑋𝑖𝑗 +1 = 𝐷 = (𝐾𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗 , 𝑦𝑖𝑗 )

ℎ𝑖𝑗 = ℎ𝑖𝑗 𝑇𝑗 , 𝑃𝑗 , 𝑋𝑖𝑗

𝐽
10. 𝐴𝐽 = 𝑉𝑗 + 𝑀=1( 𝐹,𝑚 − 𝑊𝑚 − 𝑈𝑚 )𝑉1
Matriz a resolver

𝑁 𝑁−1

𝑄𝑁 = ∑(𝐹𝑗 ℎ𝐹,𝑗 − 𝑈𝑗 ℎ𝐿,𝑗 − 𝑊ℎ𝑉,𝑗 ) − ∑(𝑄𝑗 −𝑉1 ℎ𝑉,𝑗 − 𝐿𝑁 ℎ𝐿,𝑁 )

𝑗=1 𝑗=1
𝑖

𝐿𝑗 = 𝑉𝑗 +1 + ∑ (𝐹𝑚 − 𝑈𝑚 − 𝑊𝑚 ) − 𝑉1
𝑚 =1

3. Algoritmo de solución.
 4. Tabla con variables especificadas: Fj, Zij, TFj, PF, Pj, Uj, Wj, Qj con las variables a
resolver:Xij, Yij, Lj, Vj, Tj
ITERACIÓN 6 t 0.019073519 vapor liquido
Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 300.34373 300.2999181 100 0 0.858767423 0.136547134 0.004685443 0.593429489 0.362924152 0.043646358
2 320.8202482 320.7482668 89.12557308 150 0.59352331 0.362854773 0.043621911 0.256134771 0.540620336 0.203244893
3 337.779634 337.7153679 186.4746444 139.1255731 0.374278042 0.47709643 0.148625528 0.114519706 0.466571474 0.41890882
4 352.4826746 352.4126919 185.7810184 136.4746444 0.152835101 0.541186736 0.305978163 0.03573222 0.381104119 0.583163661
5 364.1141766 364.0599086 50 135.7810184 0.045627683 0.428585788 0.525786529 0.008763226 0.237546062 0.753690712

ITERACIÓN 1 t 558.5554848 vapor liquido

Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 286.1277579 300.5952046 100 0 0.853609144 0.142012102 0.004378754 0.585577014 0.374070801 0.040352185
2 311.4469253 321.6397294 88.234587 150 0.571782277 0.384442765 0.043774958 0.242096947 0.559559925 0.198343129
3 336.7660928 339.2899293 186.5983286 138.234587 0.335463576 0.508466168 0.156070256 0.099607025 0.479333366 0.421059609
4 358.426286 354.335807 187.2607971 136.5983286 0.124190818 0.542539594 0.33326958 0.028080763 0.366806031 0.605113206
5 380.0864793 365.1783122 50 137.2607971 0.035974081 0.413082644 0.550943275 0.006784732 0.223948046 0.769267222

ITERACIÓN 2 t 2.212652243 vapor liquido

Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 300.5952046 300.7849917 100 0 0.853522445 0.141334004 0.00514355 0.58278312 0.370143365 0.047073515
2 321.6397294 321.5883583 89.33360684 150 0.582850508 0.370094508 0.047054985 0.247053572 0.539399366 0.213547062
3 339.2899293 338.4627344 186.8870823 139.3336068 0.366132461 0.478036655 0.155830883 0.110437792 0.459397266 0.430164942
4 354.335807 353.2205296 186.4303527 136.8870823 0.147128758 0.533146774 0.319724455 0.033916842 0.369050399 0.597032759
5 365.1783122 364.6824302 50 136.4303527 0.043271209 0.415780885 0.540947906 0.008226878 0.227626696 0.764146426

ITERACIÓN 3 t 0.597608013 vapor liquido

Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 300.7849917 300.5599282 100 0 0.855987023 0.139076763 0.004936214 0.587719801 0.366732217 0.045547983
2 321.5883583 321.1756858 89.2199488 150 0.588258764 0.366337718 0.045403256 0.251551439 0.539721079 0.208727483
3 338.4627344 338.098413 186.6689624 139.2199488 0.370269306 0.477376062 0.152354632 0.112465403 0.462679838 0.424854759
4 353.2205296 352.8283402 186.0984112 136.6689624 0.149915357 0.537066104 0.313018539 0.034796124 0.374871507 0.590332369
5 364.6824302 364.3822416 50 136.0984112 0.044420437 0.421965466 0.533614098 0.008486691 0.232387937 0.759125372

ITERACIÓN 4 t 0.190991217 vapor liquido

Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 300.5599282 300.4214886 100 0 0.857464847 0.13773436 0.004800794 0.590751101 0.364726144 0.044522756
2 321.1756858 320.9480474 89.16910434 150 0.591048105 0.364507558 0.044444276 0.253977611 0.540240492 0.205781897
3 338.098413 337.8946789 186.5660767 139.1691043 0.372376445 0.477267945 0.15035561 0.113547169 0.464783927 0.421668904
4 352.8283402 352.6069199 185.9295409 136.5660767 0.151453558 0.539295307 0.309251135 0.035289249 0.378204466 0.586506285
5 364.3822416 364.2105677 50 135.9295409 0.04505737 0.425507467 0.529435163 0.008632473 0.235136629 0.756230898

ITERACIÓN 5 t 0.060342513 vapor liquido

Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 300.4214886 300.34373 100 0 0.85829762 0.136975708 0.004726672 0.592462563 0.363577439 0.043959999
2 320.9480474 320.8202482 89.14175115 150 0.592629198 0.363454424 0.043916359 0.255354795 0.540490447 0.204154758
3 337.8946789 337.779634 186.507199 139.1417512 0.373596495 0.477159205 0.149244299 0.114170477 0.465931032 0.419898492
4 352.6069199 352.4826746 185.8345768 136.507199 0.152333126 0.540513115 0.307153759 0.035571318 0.380062644 0.584366037
5 364.2105677 364.1141766 50 135.8345768 0.045420731 0.42748063 0.527098639 0.008715766 0.236678714 0.75460552

ITERACIÓN 6 t 0.019073519 vapor liquido

Etapa T°k sup T°k cal L lbmol/hr| V lbmol/h propano N-butano Npentano propano N-butano Npentano
1 300.34373 300.2999181 100 0 0.858767423 0.136547134 0.004685443 0.593429489 0.362924152 0.043646358
2 320.8202482 320.7482668 89.12557308 150 0.59352331 0.362854773 0.043621911 0.256134771 0.540620336 0.203244893
3 337.779634 337.7153679 186.4746444 139.1255731 0.374278042 0.47709643 0.148625528 0.114519706 0.466571474 0.41890882
4 352.4826746 352.4126919 185.7810184 136.4746444 0.152835101 0.541186736 0.305978163 0.03573222 0.381104119 0.583163661
5 364.1141766 364.0599086 50 135.7810184 0.045627683 0.428585788 0.525786529 0.008763226 0.237546062 0.753690712
5. Perfiles de composición (número de etapas vs x y número de etapas vs y) y perfil de
Temperatura

N vs T (°k
370

360

350

340
T°K

330

320

310

300
1 1.5 2 2.5 3 3.5 4 4.5 5 5.5
Etapas

N vs XY
1
0.9
0.8
0.7
0.6
0.5
Y

0.4
0.3
0.2
0.1
0
1 1.5 2 2.5 3 3.5 4 4.5 5 5.5
Etapas

6. Cálculo de las cargas térmicas en rehervidor y condensador (explicar cómo se calculan)

𝑁 𝑁−1

𝑄𝑁 = ∑(𝐹𝑗 ℎ𝐹,𝑗 − 𝑈𝑗 ℎ𝐿,𝑗 − 𝑊ℎ𝑉,𝑗 ) − ∑(𝑄𝑗 −𝑉1 ℎ𝑉,𝑗 − 𝐿𝑁 ℎ𝐿,𝑁 )

𝑗=1 𝑗=1

Para realizar el calculo del rehervidor y condensador se debe usar la ecuación de Qn y

realizar el cálculo de las incógnitas faltantes como hF,j; hL,j; hv,j.
𝑏 𝑐 𝑑
𝐻𝑖 = 𝑎 ∗ 𝑇 + ∗ 𝑇2 + ∗ 𝑇3 + ∗ 𝑇4
2 3 4
𝐻𝑣𝑖 = ℎ𝑖 ∗ 𝑦𝑖
 𝑅𝑇 2 𝐵
𝜆=
(𝑇 + 𝐶)2

Se debe calcular λvap con la ecuación de Antoine. Se procede a usar la eci=u para calcular
Qn tomando como 1 la etapa del rehervidor y Q5 del condensador siendo esta la ultima
etapa 5.
7. 4 capturas de pantalla, donde expliques qué es lo que estás calculando en cada una de ellas
Calculo de T supuesta para la matriz y calculo de Burbuja

Resolución de la matriz A y B para encontrar la composición en cada etapa.

Calculo de burbuja para cada etapa

 Calculo de Qn, entalpias.

8. Conclusiones
En el presente escrito se desarrolló el cálculo de ecuaciones mesh. Estas ecuaciones son útiles para
evaluar diseños previos de destilación en mezclas multicomponentes y permiten calcular los flujos
de vapor y líquido en cada etapa. Para resolver estas ecuaciones, se utilizan métodos de
convergencia y se deben tomar en cuenta las propiedades no lineales de la relación de equilibrio.
Este método proporciona una mayor precisión en los cálculos y permite verificar la adecuación del
diseño.

BIBLIOGRAFIAS
https://www.youtube.com/watch?v=tL718a0aWKc&t=799s
 𝑄1 = 𝐹1 ℎ𝐹,1 + 𝑉2 ℎ𝑉,2 − (𝐿1 + 𝑈1 )ℎ𝐿,1 − 𝑉1 + 𝑊𝑗 ℎ𝑉,1

𝑄1 = 𝐹1 ℎ𝐹,1 + 𝑉2 ℎ𝑉,2 − (𝐿1 + 𝑈1 )ℎ𝐿,1 − 𝑉1 + 𝑊𝑗 ℎ𝑉,1

𝑖

𝐿𝑗 = 𝑉𝑗 +1 + ∑ (𝐹𝑚 − 𝑈𝑚 − 𝑊𝑚 ) − 𝑉1
𝑚 =1
𝑁 𝑁−1

𝑄𝑁 = ∑(𝐹𝑗 ℎ𝐹,𝑗 − 𝑈𝑗 ℎ𝐿,𝑗 − 𝑊ℎ𝑉,𝑗 ) − ∑(𝑄𝑗 −𝑉1 ℎ𝑉,𝑗 − 𝐿𝑁 ℎ𝐿,𝑁 )

𝑗=1 𝑗=1