Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 42 vistas

    Practica N4

    Cargado por

    ERIKA FERNANDA DUEÑAS CHANA
    El documento describe cuatro prácticas de simulación de procesos: 1) Determinación del equilibrio químico en un reactor mediante la aplicación de la ley de acción de masas y el método de Newton-Raphson. 2) Cálculo de concentraciones en tanques en serie usando el método de Runge-Kutta. 3) Simulación de un evaporador simple para concentrar jugo de naranja. 4) Cálculo de parámetros de un evaporador múltiple para disoluciones de NaOH.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Practica N4

    Cargado por

    ERIKA FERNANDA DUEÑAS CHANA
    42 vistas5 páginas
    El documento describe cuatro prácticas de simulación de procesos: 1) Determinación del equilibrio químico en un reactor mediante la aplicación de la ley de acción de masas y el método de Newton-Raphson. 2) Cálculo de concentraciones en tanques en serie usando el método de Runge-Kutta. 3) Simulación de un evaporador simple para concentrar jugo de naranja. 4) Cálculo de parámetros de un evaporador múltiple para disoluciones de NaOH.

    Título original

    PRACTICA N4

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe cuatro prácticas de simulación de procesos: 1) Determinación del equilibrio químico en un reactor mediante la aplicación de la ley de acción de masas y el método de Newton-Raphson. 2) Cálculo de concentraciones en tanques en serie usando el método de Runge-Kutta. 3) Simulación de un evaporador simple para concentrar jugo de naranja. 4) Cálculo de parámetros de un evaporador múltiple para disoluciones de NaOH.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    42 vistas5 páginas

    Practica N4

    Cargado por

    ERIKA FERNANDA DUEÑAS CHANA
    El documento describe cuatro prácticas de simulación de procesos: 1) Determinación del equilibrio químico en un reactor mediante la aplicación de la ley de acción de masas y el método de Newton-Raphson. 2) Cálculo de concentraciones en tanques en serie usando el método de Runge-Kutta. 3) Simulación de un evaporador simple para concentrar jugo de naranja. 4) Cálculo de parámetros de un evaporador múltiple para disoluciones de NaOH.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    Está en la página 1de 5

    PRACTICAS DE SIMULACION DE PROCESOS

    PRACTICA N° 4

    SIMULACIÓN DE FENOMENOS DE TRANSPORTE DE MASA

    1. DETERMINACIÓN DE EQUILIBRIO EN UN REACTOR


    En un reactor se efectúan las siguientes reacciones en fase gaseosa:
    A+B ⇄ C+D (1)
    A + C ⇄ 2E (2)
    A la temperatura de la reacción, las constantes de equilibrio son kp1=2.6 y kp2 = 3.1.
    Las composiciones iniciales son 2 mol/l de A y 1 mol/l de B.
    Calcule la composición a la salida del reactor, asumiendo que se alcanza el
    equilibrio. Con la aplicación de la ley de acción de masas se obtiene
    Para la reacción (1)
    (𝑥 − 𝑦)(𝑥)
    2.6 =
    (2 − 𝑥 − 𝑦)(1 − 𝑥)
    Para la reacción (2)
    (2𝑦)2
    3.1 =
    (2 − x − y)(x − y)
    que es un sistema de dos ecuaciones no lineales con dos incógnitas, cuya solución
    por el método de Newton-Raphson, por ejemplo, exige:
    ➢ Un vector inicial cercano a la solución, obtenible a partir de consideraciones
    físicas del problema.
    ➢ La matriz jacobiana, ampliada con el vector de funciones, que es
    relativamente fácil, puesto que las derivadas parciales son directas.
    Vector inicial. En virtud de las funciones y la existencia inicial de 2 moles de A y 1
    mol de B, se propone x = 0.8 y y = 0.4.
    Las derivadas parciales para la matriz Jacobiana se dan a continuacion.
    (𝑥 − 𝑦)(𝑥)
    𝑓1 (𝑥, 𝑦) = − 2.6
    (2 − 𝑥 − 𝑦)(1 − 𝑥)
    (2𝑦)2
    𝑓2 (𝑥, 𝑦) = − 3.1
    (2 − x − y)(x − y)
    𝜕𝑓1 (𝑥, 𝑦) (2 − 𝑥 − 𝑦)(1 − 𝑥)(2𝑥 − 𝑦) − (𝑥 − 𝑦)(𝑥)(−3 + 2𝑥 + 𝑦)
    =
    𝜕𝑥 ((2 − 𝑥 − 𝑦)(1 − 𝑥))2
    𝜕𝑓1 (𝑥, 𝑦) (2 − 𝑥 − 𝑦)(1 − 𝑥)(−𝑥) − (𝑥 − 𝑦)(𝑥)(𝑥 − 1)
    =
    𝜕𝑦 ((2 − 𝑥 − 𝑦)(1 − 𝑥))2
    𝜕𝑓2 (𝑥, 𝑦) −(2𝑦)2 (2 − 2𝑥)
    = 2
    𝜕𝑥 ((2 − x − y)(x − y))
    𝜕𝑓2 (𝑥, 𝑦) 8(2 − 𝑥 − 𝑦)(𝑥 − 𝑦)𝑦 − 8𝑦 2 (−1 + 𝑦)
    = 2
    𝜕𝑦 ((2 − x − y)(x − y))

    1
    DR. ROLANDO S. BASURCO CARPIO

    2. DETERMINACION DE CONCENTRACION EN TANQUES EN SERIE


    En el diagrama siguiente se toma una corriente de recirculación de 150 l/min a la
    salida del tanque 3 y se lleva al tanque 2, en tanto el volumen se mantiene constante
    en cada tanque e igual a 1000 litros, determine la concentración en cada tanque 10
    minutos después de iniciado el proceso.
    150 l/min
    300 l/min 300 l/min
    450 l/min
    50 g/l

    V1 V2 V3
    C1 C2 C3 300 l/min

    C1(0) =30 g/l C2(0) =30 g/l C3(0) =30 g/l


    V1 =1000 l V 2 =1000 l V3 =1000 l

    Hacemos un Balance de Masa en cada Tanque


    Acumulación = Entrada - Salida
    𝑑(𝑉𝐶1 )
    = 𝐸𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑆𝑎𝑙𝑖𝑑𝑎 = (300)(50) − (300)(𝐶1 )
    𝑑𝑡
    𝑑𝐶1
    = 15 − (0.3)(𝐶1 )
    𝑑𝑡
    𝑑(𝑉𝐶2 )
    = 𝐸𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑆𝑎𝑙𝑖𝑑𝑎 = (300)(𝐶1 ) + (150)(𝐶3 ) − (450)(𝐶2 )
    𝑑𝑡
    𝑑𝐶2
    = (0.3)(𝐶1 ) − (0.45)(𝐶2 ) + (0.15)(𝐶3 )
    𝑑𝑡
    𝑑(𝑉𝐶3 )
    = 𝐸𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑆𝑎𝑙𝑖𝑑𝑎 = (450)(𝐶2 ) − (450)(𝐶3 )
    𝑑𝑡
    𝑑𝐶3
    = (0.45)(𝐶2 ) − (0.45)(𝐶3 )
    𝑑𝑡
    Considerando h = 1 minuto, aplicamos el método de Runge Kutta, donde:
    𝑓1 = 15 − (0.3)(𝐶1 )
    𝑓2 = (0.3)(𝐶1 ) − (0.45)(𝐶2 ) + (0.15)(𝐶3 )
    𝑓3 = (0.45)(𝐶2 ) − (0.45)(𝐶3 )
    Considerando los puntos iniciales como C1(0)=30 g/l, C2(0)=30 g/l, C3(0)=30 g/l y el
    Intervalo a definir t= [0, 10], calcularemos las constantes de Runge Kutta para cada
    función:
    Es probable que uno de los procedimientos más difundidos y a la vez más exactos
    para obtener soluciones aproximadas al problema de valor inicial y’ = f(x, y), y(xo)
    = yo sea el método de Runge-Kutta de cuarto orden.
    El algoritmo recursivo es:
    𝑥𝑛+1 = 𝑥𝑛 + ℎ
    yn+1 = yn + (1/6) (k1 + 2k2 + 2k3 + k4) , n = 0, 1, 2, . . .

    2
    PRACTICAS DE SIMULACION DE PROCESOS

    k1 = hf(xn, yn)
    k2 = hf(xn + ½ h, yn + ½ k1)
    k3 = hf(xn + ½ h, yn + ½ k2)
    k4 = hf(xn + h, yn + k3)
    K11= f1(ti, C1i, C2i, C3i) = 15 – 0.3(30) = 6
    K12= f2(ti, C1i, C2i, C3i) = 0.3(30) – 0.45(30) +0.15(30) = 0
    K13= f2(ti, C1i, C2i, C3i) = 0.45(30) – 0.45(30) = 0
    K21= f1(ti, C1i+(h/2)(K11), C2i+(h/2)(K12), C3i+(h/2)(K13)) = 15 – 0.3(30 + 0.5*6) =
    5.1
    K22= f2(ti, C1i+(h/2)(K11), C2i+(h/2)(K12), C3i+(h/2)(K13)) = 0.3(30+0.5*6) – 0.45(30
    + 0.5*0) + 0.15(30 + 0.5*0) = 0.9
    K23= f2(ti, C1i+(h/2)(K11), C2i+(h/2)(K12), C3i+(h/2)(K13)) = 0.45(30+0.5*0) – 0.45(30
    + 0.5*0) = 0
    K31= f1(ti, C1i+(h/2)(K21), C2i+(h/2)(K22), C3i+(h/2)(K23)) = 15 – 0.3(30 + 0.5*30.6)
    = 5.2350
    K32= f2(ti, C1i+(h/2)(K21), C2i+(h/2)(K22), C3i+(h/2)(K23)) = 0.3(30+0.5*5.1) –
    0.45(30 + 0.5*0.9) + 0.15(30 + 0.5*0) = 0.5625
    K33= f3(ti, C1i+(h/2)(K21), C2i+(h/2)(K22), C3i+(h/2)(K23))= 0.45(30+0.5*0.9) –
    0.45(30 + 0.5*0) = 0.2025
    K41= f1(ti, C1i+(h)(K31), C2i+(h)(K32), C3i+(h/2)(K33)) = 15 – 0.3(30 + 1*5.2350) =
    4.4295
    K42= f2(ti, C1i+(h)(K31), C2i+(h(K32), C3i+(h)(K33)) = 0.3(30+1*5.2350) – 0.45(30 +
    1*0.5625) + 0.15(30 + 1*0.2025) = 1.3478
    K43= f3(ti, C1i+(h)(K31), C2i+(h)(K32), C3i+(h/2)(K33))= 0.45(30+1*0.5625) – 0.45(30
    + 1*0.2025) = 0.1620

    t1 = t0 + h = 0 + 1 =1
    C1 = C1(0) + (h/6) (K11 +2k21 +2 k31 + k41) = 30 + 1* (6 + 2*5.1 + 2*5.2350 +
    4.4295)/6 = 35.1883
    C2 = C2(0) + (h/6) (K12 +2k22 +2 k32 + k42) = 30 + 1* (0 + 2*0.9 + 2*0.5625 –
    1.3478) / 6 = 30.7121
    C3 = C3(0) + (h/6) (K13 +2k23 +2 k33 + k43) = 30 + 1* (0 + 2*0 + 2*0.2025 + 0.1620)
    / 6 = 30.0945
    3. SIMULACION DE UN EVAPORADOR
    Se debe diseñar un evaporador de efecto simple para concentrar una solución de
    jugo de naranja al 20% (por peso) hasta una solución al 50%. La solución diluida
    (el alimento) a 200ºF debe de alimentarse al evaporador a un flujo de 40000 lb/h
    para fines de calentamiento, se utiliza vapor de agua saturado a 350ºF. Se dispone
    de suficiente área en el condensador para mantener una presión de 0.9492 lb/pulg2
    (absoluta)en el espacio de vapor del condensador. En base a un coeficiente total
    de transferencia de calor de 300 Btu/h ft2 ºF calcular:
    3
    DR. ROLANDO S. BASURCO CARPIO

    La Entalpia del alimento es 262.8 BTU/lb y


    del licor concentrado 99.45 BTU/lb.
    a) El área de calentamiento requerida
    b) El consumo y la economía del vapor de
    agua.
    BALANCE DE MATERIA
    Balance del soluto: FXf= LxL
    Balance del solvente: F(1 – Xf ) = V + L(1 –
    xL)
    Balance de materia total: F = V + L
    BALANCE ENTÁLPICO ALREDEDOR DEL
    PROCESO
    F hf + Q – V HV – L hL = 0
    Pero como V = F - L
    F ( hf - hL) + Q – ( F – L ) (HV – hL ) = 0
    BALANCE DE LA ENTALPÍA EN EL VAPOR
    DE AGUA
    VoHo – Q – Voho = 0
    Q = Vo(Ho – ho) = Vo hfg
    VELOCIDAD DE TRANSFERENCIA DE CALOR Q
    Q = UA(To – T)
    4. TAREA DE LA PRÁCTICA
    1.- Problema 2 ejecutarlo en Matlab.
    2.- Evaporador múltiple para disoluciones de NaOH. - Una disolución de hidróxido
    sódico a 38 ºC y en cantidad de 5 Tm/hr de NaOH, se ha de concentrar desde el
    10% hasta el 50%, en peso, en un evaporador de triple efecto y circulación forzada;
    la alimentación es en cabeza, pasando sucesivamente por 1,2, 3 (Figura 1). El
    vapor fresco es de 1.6 atm abs. El condensador opera a un vacío de 710 mm
    referido a un barómetro de 760 mm. Los caudales serán tales que los coeficientes
    valdrán 5370, 2930 y 1950 Kcal/(hr)(m2) (ºC) en los efectos 1, 2 y 3,
    respectivamente.
    Calcular las temperaturas y las concentraciones en cada efecto. Si todos los efectos
    son iguales, calcular la superficie de calefacción de cada uno. ¿Cuál es la cantidad
    evaporada por kg de vapor de calefacción?

    4
    PRACTICAS DE SIMULACION DE PROCESOS

    DATOS TERMODINAMICOS
    42g en 100ml a 0ºC, 347g en 100ml a 100ºC
    F NaOH = 5000.kg.hr-1 z F = 0.1 x L3 = 0.5 y=0
    Número de efectos NE = 3 i = 1 a NE
    Presión en el condensador (= Presión de vapor en último efecto): p V3 = 710 torr
    Presión del vapor de entrada al primer efecto: pS = 1.6.atm
    Coeficientes globales de transmisión de calor
    Ui = 5370, 2930, 1950 kcal.m-2 hr-1 K-1

    También podría gustarte