Punto de Burbuja, Punto de Rocío, y Destilación Flash

Objetivo del artículo

El objetivo de este artículo es mostrar cómo se debe calcular la composición de un flash adiabático, y exponer el papel que juegan el punto de burbuja y el punto de rocío en el cálculo del Flash.

El artículo está inspirado en unas hojas del libro del Profesor Skogestad, de la Politécnica de Noruega a juzgar por el sitio web (http://folk.ntnu.no/skoge/bok/mer/flash_english_edition_2009), que está escrito en un lenguaje sencillo y muy entendible.

Restricciones al modelo matemático

Se ha tratado a los componentes de la mezcla mediante la Ley de Raoult; se ha estimado sus presiones de vapor mediante el modelo de Antoine, .

Donde

Ecuacion1

Ecuacion2

Mezcla que se alimenta al tambor Flash

A continuación se indica la mezcla que se utilizó, así como los respectivos coeficientes de Antoine, que se encuentran en la fuente citada.

Datos del problema

Figura 1. Datos del Problema

Para estos datos se escribió un pequeño programa, que se transcribe más adelante, que maneja tres subrutinas, quq calculan el punto de burbuja, el punto de rocío , y la composición del Flash.

Resultados

Los resultados del programa se muestran a continuación.

Resultados

Figura 2. Resultados del programa

Cálculo del punto de burbuja

En este caso se calculó el punto de burbuja a una presión dada, para lo que se trabajó con la ecuación siguiente (Skogestad);

Ecuacion3

La solución del problema se reduce a encontrar la temperatura que satisfaga la ecuación, lo que se realiza mediante un For-Next  y una condición de convergencia. También se muestra en el programa, y en la Figura 1, la composición de la primera burbuja que se forma.

Cálculo del punto de rocío

Para este propósito se considera que la mezcla de alimentación está constituida de vapor, y se busca la temperatura a la que se condensa la primera gota de líquido, lo que significa que en la composición de la mezcla que está expresada en fracciones molares Zi, éstas deben expresarse como yi,

Lo anterior significa que, para fines de cálculo se utiliza la expresión (Skogestad) siguiente:

Ecuacion4

El cálculo, dada la presión, implica encontrar una temperatura que satisfaga la ecuación anterior de acuerdo a un criterio de convergencia que se muestra en el programa. En el programa también se muestran las fracciones molares de la primera gota que se condensa.

Cálculo del Flash

Para los fines consiguientes se utilizó la ecuación siguiente con un criterio de convergencia adecuado que se puede observar en el programa, y que obtiene el V/F. Los cálculos de la composición pueden observarse en el programa.

Ecuacion5

El programa

 

Option Base 1

Option Explicit

Dim ArrayZ(1 To 5), ArrayA(1 To 5), ArrayB(1 To 5), ArrayC(1 To 5), ArrayY(1 To 5), Arrayp(1 To 5), p As Variant

Dim ArrayK(1 To 5), ArrayX(1 To 5) As Variant

Dim Summx As Variant

Dim kk, Nc, Fila, N, PP As Integer

Dim SumFPB, SumFPPB, T, Tn, Tv, M, MM, MMM, MMMM, Sump, j, Psiv, Psin, FsumFPsi, F2Psi, Sumx, Psi, Sumy, SumPsi, SumFPsi, SumFLV, SumPFLV As Variant

 

Sub destilacion_flash()

ArrayZ(1) = Worksheets(“Hoja1”).Cells(2, 3)

ArrayZ(2) = Worksheets(“Hoja1”).Cells(3, 3)

ArrayZ(3) = Worksheets(“Hoja1”).Cells(4, 3)

ArrayZ(4) = Worksheets(“Hoja1”).Cells(5, 3)

ArrayA(1) = Worksheets(“Hoja1”).Cells(6, 3)

ArrayA(2) = Worksheets(“Hoja1”).Cells(7, 3)

ArrayA(3) = Worksheets(“Hoja1”).Cells(8, 3)

ArrayA(4) = Worksheets(“Hoja1”).Cells(9, 3)

ArrayB(1) = Worksheets(“Hoja1”).Cells(10, 3)

ArrayB(2) = Worksheets(“Hoja1”).Cells(11, 3)

ArrayB(3) = Worksheets(“Hoja1”).Cells(12, 3)

ArrayB(4) = Worksheets(“Hoja1”).Cells(13, 3)

ArrayC(1) = Worksheets(“Hoja1”).Cells(14, 3)

ArrayC(2) = Worksheets(“Hoja1”).Cells(15, 3)

ArrayC(3) = Worksheets(“Hoja1”).Cells(16, 3)

ArrayC(4) = Worksheets(“Hoja1”).Cells(17, 3)

Log (2) / Log(10)

p = Worksheets(“Hoja1”).Cells(18, 3)

Nc = Worksheets(“Hoja1”).Cells(19, 3)

T = Worksheets(“Hoja1”).Cells(20, 3)

Psiv = Worksheets(“Hoja1”).Cells(21, 3)

T = T + 273.2

Tv = T

N = 0: Sump = 0

‘Limpieza

Worksheets(“Hoja1”).Range(“A25:g33000”).Clear

Fila = 25

Call Punto_de_Burbuja(p, ArrayA(), ArrayC(), ArrayZ(), ArrayB(), Tv, Tn, Fila)

Call Dew_Point(31, ArrayA(), ArrayB(), ArrayC(), ArrayZ(), ArrayY(), 5)

Call Tanteo(ArrayZ(), ArrayK(), ArrayA(), ArrayB(), ArrayC(), ArrayX(), 4, 375, 37)

End Sub

Sub Punto_de_Burbuja(p, ArrayA(), ArrayC(), ArrayZ(), ArrayB(), Tv, Tn, Fila)

SumFPPB = 0

For T = 354 To 355 Step 0.1

SumFPB = 0:

For kk = 1 To Nc

MM = ArrayZ(kk) * (10 ^ (ArrayA(kk) – (ArrayB(kk) / (T + ArrayC(kk)))))

SumFPB = SumFPB + MM

Next kk

If Abs(SumFPB – p) <= 0.05 Then

N = N + 1

Worksheets(“Hoja1”).Cells(Fila, 1) = “El Punto de Burbuja es  ” & T & ” oK” & ” ó ” & “T= ” & (T – 273.2) & ” oC ”

Exit For

End If

Next T

Sump = 0

For kk = 1 To Nc

Arrayp(kk) = ArrayZ(kk) * (10 ^ (ArrayA(kk) – (ArrayB(kk) / (T + ArrayC(kk)))))

Worksheets(“Hoja1″).Cells(Fila + kk, 1) = ” y(” & kk & “)= “: Worksheets(“Hoja1”).Cells(Fila + kk, 2) = Arrayp(kk) / p

Sump = Sump + (Arrayp(kk) / p)

‘ArrayY(kk) = Arrayp(kk) / p

Next kk

Worksheets(“Hoja1”).Cells(Fila + kk, 2) = Sump

Fila = Fila + 6

Exit Sub

Fila = Fila + 1

End Sub

Sub Tanteo(ArrayZ(), ArrayK(), ArrayA(), ArrayB(), ArrayC(), ArraX(), p, T, Fila)

Worksheets(“Hoja1”).Cells(Fila, 1) = “Flash:”

For Psi = 0.02 To 1 Step 0.02

Sumx = 0: Sumy = 0

For kk = 1 To Nc

ArrayK(kk) = (10 ^ (ArrayA(kk) – ((ArrayB(kk) / (T + ArrayC(kk)))))) / p

ArrayX(kk) = ArrayZ(kk) / (1 + Psi * (ArrayK(kk) – 1))

ArrayY(kk) = ArrayX(kk) * ArrayK(kk)

Sumy = Sumy + ArrayY(kk)

Sumx = Sumx + ArrayX(kk)

Next kk

If Abs(Sumx – 1) < 0.001 Then Exit For

Next Psi

Worksheets(“Hoja1”).Cells(Fila + 1, 1) = “V/F=” & Psi: Worksheets(“Hoja1”).Cells(Fila + 1, 2) = “T= ” & T & ” oK o ” & T – 273.2 & “oC”

Worksheets(“Hoja1”).Cells(Fila + 1, 3) = “p=” & p

For kk = 1 To Nc

Worksheets(“Hoja1”).Range(“A43:B43”).HorizontalAlignment = xlLeft

Worksheets(“Hoja1”).Cells(Fila + 1 + kk, 1) = “x(” & kk & “)= ” & ArrayX(kk): Worksheets(“Hoja1”).Cells(Fila + 1 + kk, 2) = “y(” & kk & “)= ” & ArrayY(kk)

Next kk

Worksheets(“Hoja1”).Cells(Fila + 1 + kk, 1) = Sumx: Worksheets(“Hoja1”).Cells(Fila + 1 + kk, 2) = Sumy

End Sub

Sub Dew_Point(Fila, ArrayA(), ArrayB(), ArrayC(), ArrayZ(), Arrayp(), p)

For T = 354.3 To 410

Sump = 0: Summx = 0

For kk = 1 To Nc

Sump = Sump + ArrayZ(kk) / 10 ^ (ArrayA(kk) – (ArrayB(kk) / (T + ArrayC(kk))))

Next kk

If Abs(Sump – (1 / p)) < 0.005 Then

Worksheets(“Hoja1”).Cells(Fila, 1) = “El Punto de Rocío es T= ” & T & “oK” & ” ó ” & T – 273.2 & “oC”

Exit For

End If

Next T

For kk = 1 To Nc

Worksheets(“Hoja1”).Cells(Fila + kk, 2) = “x(” & kk & “)=” & (ArrayZ(kk) * p) / 10 ^ (ArrayA(kk) – (ArrayB(kk) / (T + ArrayC(kk))))

Summx = Summx + (ArrayZ(kk) * p) / 10 ^ (ArrayA(kk) – (ArrayB(kk) / (T + ArrayC(kk))))

Next kk

Worksheets(“Hoja1”).Range(“b36”).HorizontalAlignment = xlLeft

Worksheets(“Hoja1”).Cells(Fila + kk, 2) = Summx

End Sub

 

 

 

Acerca de

Professor of modeling and simulation, and process design at Escuela Politécnica Nacional, in Quito, Ecuador. . In the past I was a P4, P5, and D1 at the Organization for the Prohibition of Chemical Weapons, located in the Kingdom pf the Netherlands

Publicado en Sin categoría

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s

Enter your email address to follow this blog and receive notifications of new posts by email.

Únete a otros 1.772 seguidores

Categorías
Artículos y comentarios sobre modelado y simulación de plantas, y equipos de la industria química, con ejemplos
PREGUNTAS O INQUIETUDES
PUEDEN ENVIAR SUS PREGUNTAS/INQUIETUDES RESPECTO DE ARTÍCULOS DEL BLOG, O TEMAS AFINES A gasteaux@hotmail.com ASUNTO: BLOG DE INGENIERÍA QUÍMICA
A %d blogueros les gusta esto: