UNIVERSIDAD NACIONAL DE INGENIERIA

FACULTAD DE INGENIERIA QUIMICA

DEPARTAMENTO DE OPERACIONES UNITARIAS

INGENIERIA QUIMICA.

TRABAJO DE CURSO.

Elaborado por:
Alejandra Cáceres Jarquín.
Jassir Barahona Colindres.

Grupo:
5T2-IQ

Docente:
Ing. Javier Enrique Ramírez Meza

Asignatura:
Control de Procesos.

Fecha:
09 de julio 2020
 1. Introducción

1.1. Descripción del proceso

Ejercicio 9.4; Corripio

El evaporador que se ilustra en la figura 9-27 se utiliza para concentrar una solución de azúcar
desde una fracción de masa de alimentación XF, hasta una fracción de masa de producto XP.
Como medio de calefacción se dispone de vapor saturado a 20 psia; el agua de enfriamiento
entra al condensador barométrico con una temperatura TW. Se puede suponer que la presión
en el evaporador es la presión de vapor del agua a la temperatura de salida del condensador,
TC, que la mezcla del contenido del evaporador es perfecta y que la temperatura de los tubos
de calandria es igual a la del vapor, que se condensa, TS. Las pérdidas de calor se pueden
despreciar.

El área de transferencia del evaporador es de 309.5 pies2, y su capacidad es de 4.2 pies3 de

solución por pie de nivel. La capacidad calórica de los tubos de calandria es 154 Btu/oF, el
coeficiente de transferencia de calor es 150 Btu/hpie2- oF. El volumen de la solución abajo de
los tubos de transferencia de calor es de 4.68 pies3.

Las válvulas de control tienen características lineales, constantes de tiempo de 5 seg y se

dimensionan para una sobrecapacidad del 100%. La caída de presión a través de la válvula
de control de alimentación es constante a 16 psi, y a través de la válvula de producto es igual
a la presión en el evaporador, debido a que con un soporte barométrico se hace que la presión,
en el lado “corriente abajo” de la válvula de producto, sea cero absolutos. El sensor de nivel
es un flotador con rango de 0 a 2 pies por arriba del fondo de los tubos. El rango del transmisor
de concentración (AT20) es de 0.3 a 0.8 unidades de fracción de masa, y la constante de
tiempo, de 25 seg. El controlador de nivel (LC20) es proporcional, con un punto de control
de 25% del rango, y se ajusta para un control estricto de nivel (banda proporcional entre 1 y
10%). El controlador de concentración (ARC20) es proporcional-integra (PI).

Las condiciones de diseño son:

Fracción de masa de alimentación: 0.10

Temperatura de alimentación: 90°F

Fracción de masa del producto (punto de control): 0.70

Temperatura de salida del condensador: 105oF

Temperatura de entrada del agua de enfriamiento: 90°F

Las siguientes correlaciones se pueden utilizar para calcular las propiedades físicas de la
solución y del vapor:

Entalpia de la solución: H = (1 - 0.7x)(T - TR) Btu/lb

Densidad de la solución: p = 62.384 + 23.14x+ 11.27x2 lb/pies3

Elevación del punto de ebullición: BPR = 1.13 - 6.59x + 26.9x2 oF

donde: x es la fracción de masa de azúcar

T es la temperatura, oF

TR es la temperatura de referencia, oF

La capacidad calórica del agua líquida es 1.0 Btu/lb- oF, y la del vapor de agua es 0.46 Btu/lb-
o
F; estas capacidades se pueden suponer constantes y se utilizan para corregir el calor de
vaporización latente del agua por cambios en la temperatura. La ecuación de Antoine se
puede utilizar para calcular la presión de vapor del agua en función de la temperatura, y
viceversa

3816.44
Para el agua 𝑙𝑛𝑃 𝑜 = 18.3036 −
𝑇+227

Donde Po está en mm de mercurio T en C

El evaporador se encuentra inicialmente a las condiciones de diseño de estado estacionario.

Las posibles perturbaciones son la razón del vapor, la fracción de masa de alimentación, la
razón de agua para enfriamiento que entra al condensador, la temperatura del agua para
enfriamiento y el punto de control de la composición del producto.
a) Se deben escribir las ecuaciones del modelo y calcular las condiciones de diseño de estado
estacionario. Se deben programar las ecuaciones para resolverlas en una computadora y
verificar que las condiciones iniciales constituyen un estado estacionario.

b) El controlador de nivel se debe ajustar para una respuesta de razón de asentamiento de un

cuarto mediante ensayo y error. El controlador de composición se debe ajustar con el método
de ganancia última de Ziegler-Nichols.

c) Se deben obtener las respuestas a entradas escalón’ en el punto de control de la

composición y en cada una de las perturbaciones que se listan en el enunciado del problema.
Los resultados se deben comentar brevemente.

Datos:

Pvsat = 20 Psia  =  seg

Ttcalan = TS Válvula dimensionada para una sobrecapacidad
del 100%
Las pérdidas de calor se desprecian. Pvalv = 16 psi
Aevap = 309.5 pies2 Rflot = 0 a 2 pies
Vevap = 4.2 pies3 de solución por pie RAT20 = 0.3 a 0.8 unidades de fracción de masa
de nivel
Cptcalan = 154 Btu/oF  = 25 seg
U = 4.68 pies3
A través de la válvula de producto es igual a la presión en el evaporador, debido a que con
un soporte barométrico se hace que la presión, en el lado “corriente abajo” de la válvula de
producto, sea cero absolutos.

El controlador de nivel (LC20) es proporcional, con un punto de control de 25% del rango,
y se ajusta para un control estricto de nivel (banda proporcional entre 1 y 10%). El
controlador de concentración (ARC20) es proporcional-integra (PI).

Las condiciones de diseño son: XF = 0.10; TF = 90°F; XP = 0.70; TC = 105oF y TW = 90°F

Las siguientes correlaciones se pueden utilizar para calcular las propiedades físicas de la
solución y del vapor:
Propiedad Ecuación unidades
Entalpia de la solución: H = (1 - 0.7x) (T - TR) [Btu/lb]
Densidad de la solución: p = 62.384 + 23.14x+ 11.27x2 [lb/pies3]
Elevación del punto de ebullición: BPR = 1.13 - 6.59x + 26.9x2 [oF]

donde: x es la fracción de masa de azúcar (-), T es la temperatura, (oF), y TR es la temperatura

de referencia, (oF)
2. Implementación de balances para obtener el modelo
matemático

2.1. Identificación del objetivo de control, las variables de control,

variables perturbada, y posibles variables manipuladas.

Objetivo de Control

Mantener el nivel de la solución en el evaporador y la mantener la fracción másica en la

salida en el setpoint.

Variable de Control:

• El nivel (h) de la solución en el evaporador.

• La concentración (Xp) de la solución a la salida del evaporador.

Variable Manipulada:

• El flujo de la solución en la alimentación.

• El flujo de la solución en la salida.

Variable Perturbadora:

• La razón del vapor.

• La razón de agua para enfriamiento que entra al condensador.
• La temperatura del agua para enfriamiento.
2.2. Obtención del modelo matemático

2.2.1. Breve explicación del por qué se seleccionó el balance que implementó, esto debe
estar relacionado con el objetivo de control.

Se selecciono un Balance Masico en el evaporador debido que se controlara el nivel de la

solución en él, así mismo un balance de solutos ya que se debe controla la fracción másica a
la salida del evaporador.

2.2.2. Implementación de los balances de materia y/o energía acorde con la(s)
variable(s) de control identificadas.

Balance de Materia Total en el Evaporador

Estado estacionario

𝑑𝑚̇𝑠
= 𝐹𝑠 − 𝐿𝑃𝑠 − 𝐺2 𝑠 𝐸𝑐. 1
𝑑𝑡

El flujo masico es:

𝑚̇ = 𝜌 ∗ 𝑉 𝐸𝑐. 2

𝑑𝑚̇𝑠 𝑑(𝑝𝑠 ∗ 𝑉𝑠 )
= = 𝐹𝑠 − 𝐿𝑃𝑠 − 𝐺2 𝐸𝑐. 3
𝑑𝑡 𝑑𝑡

𝑉 = 𝐴∗ℎ 𝐸𝑐. 4

Derivada de un producto:

𝑑𝑚̇ 𝑑(𝜌 ∗ ℎ)
=𝐴 = 𝐴[(𝜌′ 𝑠 ∗ ℎ𝑠 ) + (𝜌𝑠 ∗ ℎ′ 𝑠 )] 𝐸𝑐. 5
𝑑𝑡 𝑑𝑡
 • Sustituyendo Ec.5 en Ec.3 se obtiene:

𝐴[(𝜌′ 𝑠 ∗ ℎ𝑠 ) + (𝜌𝑠 ∗ ℎ′ 𝑠 )] = 𝐹𝑠 − 𝐿𝑃𝑠 𝐸𝑐. 6

(𝐴 ∗ ℎ𝑠 )𝜌′ 𝑠 + (𝐴 ∗ 𝜌𝑠 )ℎ′ 𝑠 = 𝐹𝑠 − 𝐿𝑃𝑠 𝐸𝑐. 7

Estado no estacionario

• Aplicando la Ec.7 para el balance en estado no estacionario, se obtiene:

(𝐴 ∗ ℎ𝑠 )𝜌′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ′ (𝑡) = 𝐹(𝑡) − 𝐿𝑃 (𝑡) 𝐸𝑐. 8

El flujo de líquido a través de una válvula de globo está dado por la siguiente ecuación:

𝐹(𝑡) = 500𝐶𝑉 (𝑣𝑝(𝑡))√Δ𝑃(𝑡) ∗ 𝐺 𝐸𝑐. 9

𝐿𝑃 (𝑡) = 500𝐶𝑉 (𝑣𝑝(𝑡))√Δ𝑃(𝑡) ∗ 𝐺 𝐸𝑐. 10

• Sustituyendo en Ec.8 la Ec.9 y Ec.10 donde se obtiene:

(𝐴 ∗ ℎ𝑠 )𝜌′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ′ (𝑡) = [500𝐶𝑉 (𝑣𝑝(𝑡))√Δ𝑃(𝑡) ∗ 𝐺] [500𝐶𝑉 (𝑣𝑝(𝑡))√Δ𝑃(𝑡) ∗ 𝐺] 𝐸𝑐. 11

𝐴𝑙𝑖𝑚 𝑃𝑟𝑜𝑑

• Sustituyendo en Ec.10 la Ec.11 se obtiene:

(𝐴 ∗ ℎ𝑠 )𝜌′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ′ (𝑡) = [500𝐶𝑉 (𝑣𝑝(𝑡))√Δ𝑃(𝑡) ∗ 𝐺] − [500𝐶𝑉 (𝑣𝑝(𝑡))√Δ𝑃(𝑡) ∗ 𝐺] − 𝐺2 𝐸𝑐. 12

𝐹 𝑃
Linealizando

Para el flujo a través de la válvula de alimentación:

𝜑𝐹 𝜑𝐹
𝐹(𝑡) ≈ 𝐹̅ + | (𝑣𝑝(𝑡) − 𝑣𝑝
̅̅̅̅) + (ℎ(𝑡) − ℎ̅) 𝐸𝑐. 13
𝜑𝑣𝑝 𝜑ℎ
1
̅̅̅̅
500𝐶𝑉 Δ𝑃(𝑡)𝐺𝑣𝑝
𝐹(𝑡) ≈ 𝐹̅ + 500𝐶𝑉 √Δ𝑃(𝑡) ∗ 𝐺 ∗ (𝑣𝑝(𝑡) − 𝑣𝑝
̅̅̅̅) + [Δ𝑃(𝑡) ∗ 𝐺]−2 (ℎ(𝑡) − ℎ̅ ) 𝐸𝑐. 14
2

Para simplificar la notación:

𝐶1 = 500𝐶𝑉 √Δ𝑃(𝑡) ∗ 𝐺 𝐸𝑐. 15

500𝐶𝑉 Δ𝑃(𝑡)𝐺𝑣𝑝
̅̅̅̅ 1
𝐶2 = [Δ𝑃(𝑡) ∗ 𝐺]−2 𝐸𝑐. 16
2

• Sustituyendo en Ec.14 la Ec.15 y Ec.16 se obtiene:

̅̅̅̅) + 𝐶2 (ℎ(𝑡) − ℎ̅)

𝐹(𝑡) ≈ 𝐹̅ + 𝐶1 (𝑣𝑝(𝑡) − 𝑣𝑝 𝐸𝑐. 17

Para el flujo a través de la válvula de producto:

𝜑𝐿𝑃 𝜑𝐿𝑃
𝐿𝑃 (𝑡) ≈ 𝑃̅ + | (𝑣𝑝(𝑡) − 𝑣𝑝
̅̅̅̅) + (ℎ(𝑡) − ℎ̅) 𝐸𝑐. 18
𝜑𝑣𝑝 𝜑ℎ

500𝐶𝑉 Δ𝑃(𝑡)𝐺𝑣𝑝
̅̅̅̅ 1
𝐿𝑃 (𝑡) ≈ 𝐿̅𝑃 + 500𝐶𝑉 √Δ𝑃(𝑡) ∗ 𝐺 ∗ (𝑣𝑝(𝑡) − ̅𝑣𝑝
̅̅̅) + [Δ𝑃(𝑡) ∗ 𝐺]−2 (ℎ(𝑡) − ℎ̅ ) 𝐸𝑐. 19
2

Para simplificar la notación:

𝐶3 = 500𝐶𝑉 √Δ𝑃(𝑡) ∗ 𝐺 𝐸𝑐. 20

500𝐶𝑉 Δ𝑃(𝑡)𝐺𝑣𝑝
̅̅̅̅ 1
𝐶4 = [Δ𝑃(𝑡) ∗ 𝐺]−2 𝐸𝑐. 21
2
 • Sustituyendo en Ec.19 la Ec.20 y Ec.21 se obtiene:

̅̅̅̅) + 𝐶4 (ℎ(𝑡) − ℎ̅)

𝐿𝑃 (𝑡) ≈ 𝐿̅𝑃 + 𝐶3 (𝑣𝑝(𝑡) − 𝑣𝑝 𝐸𝑐. 22

• Sustituyendo en Ec.13 la Ec.18 y Ec.23 se obtiene:

̅̅̅) + 𝐶2 (ℎ(𝑡) − ℎ̅)] − [𝐿̅𝑃 + 𝐶3 (𝑣𝑝(𝑡) − ̅𝑣𝑝

(𝐴 ∗ ℎ𝑠 )𝜌 ′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ′ (𝑡) = [𝐹̅ + 𝐶1 (𝑣𝑝(𝑡) − ̅𝑣𝑝 ̅̅̅) + 𝐶4 (ℎ(𝑡) − ℎ̅)]−𝐺2 (𝑡) 𝐸𝑐. 23

• Al sustraer la ecuación de estado estacionario Ec.6 en la Ec.23 se obtiene:

(𝐴 ∗ ℎ𝑠 )𝜌̅′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ̅′ (𝑡) = [𝐹̅ + 𝐶1(𝑣𝑝(𝑡) − 𝑣𝑝

̅̅̅̅) + 𝐶2(ℎ(𝑡) − ℎ̅ )] − [𝐿̅𝑃 + 𝐶3(𝑣𝑝(𝑡) − 𝑣𝑝
̅̅̅̅) + 𝐶4(ℎ(𝑡) − ℎ̅)] − 𝐺2(𝑡) 𝐸𝑐.24

(𝐴 ∗ ℎ𝑠 )𝜌̅′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ̅′ (𝑡) = 𝐹̅ + 𝐶1(𝑣𝑝(𝑡) − 𝑣𝑝

̅̅̅̅) + 𝐶2(ℎ(𝑡) − ℎ̅) − 𝐿̅𝑃 − 𝐶3(𝑣𝑝(𝑡) − 𝑣𝑝
̅̅̅̅) − 𝐶4(ℎ(𝑡) − ℎ̅) − ̅̅̅
𝐺2(𝑡) 𝐸𝑐.25

• Definiendo las variables de desviación

̅̅̅̅(𝒕)𝑭 = 𝑣𝑝(𝑡) − 𝑣𝑝
𝒗𝒑 ̅̅̅̅

̅̅̅̅(𝒕)𝑷 = 𝑣𝑝(𝑡) − 𝑣𝑝
𝒗𝒑 ̅̅̅̅

̅ (𝒕) = ℎ(𝑡) − ℎ̅
𝒉

• Sustituyendo las variables de desviación en la Ec.25

(𝐴 ∗ ℎ𝑠 )𝜌̅′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ̅′ (𝑡) = 𝐹̅ + 𝐶1 𝒗𝒑 ̅ (𝒕) − 𝐿̅𝑃 − 𝐶3 𝒗𝒑

̅̅̅̅(𝒕)𝑭 + 𝐶2 𝒉 ̅ (𝒕) − ̅̅̅
̅̅̅̅(𝒕)𝑷 − 𝐶4 𝒉 𝐺2 (𝑡) 𝐸𝑐. 26

(𝐴 ∗ ℎ𝑠 )𝜌′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ′ (𝑡) = 𝐹̅ + 𝐶1 𝒗𝒑

̅̅̅̅(𝒕)𝑭 − 𝐿̅𝑃 − 𝐶3 𝒗𝒑 ̅ (𝒕) − 𝐺
̅̅̅̅(𝒕)𝑷 − (𝐶4 − 𝐶2 )𝒉 ̅̅̅2 (𝑡) 𝐸𝑐. 27

• Dividiendo la Ec.27 entre (𝐶2 − 𝐶4 ) se obtiene:

(𝐴 ∗ ℎ𝑠 ) ′ (𝐴 ∗ 𝜌𝑠 ) ′ 𝐹̅ ̅̅̅̅(𝒕)
𝒗𝒑 𝑭 𝐿̅𝑃 ̅̅̅̅(𝒕)𝑷
𝒗𝒑 ̅̅̅
𝐺2 (𝑡)
(𝐶2 − 𝐶4 )
𝜌 (𝑡) +
(𝐶2 − 𝐶4 )
ℎ (𝑡) = + − −
(𝐶2 − 𝐶4 ) (𝐶2 − 𝐶4 ) (𝐶2 − 𝐶4 ) (𝐶2 − 𝐶4 )
̅ (𝒕) −
−𝒉
(𝐶2 − 𝐶4 )
𝐸𝑐. 28
Donde:

(𝐴 ∗ ℎ𝑠 )
𝜏1 = 𝐸𝑐. 29
( 𝐶2 − 𝐶4 )

La densidad de la solución está dada por la relación 𝜌 = 62.384 + 23.14𝑥 + 11.27𝑥 2

𝐴 ∗ [62.384 + 23.14𝑥 + 11.27𝑥 2 ]

𝜏2 = 𝐸𝑐. 30
( 𝐶2 − 𝐶4 )

1
𝐾𝑃1 = 𝐸𝑐. 31
( 𝐶2 − 𝐶4 )

• Sustituyendo Ec.29, Ec.30 y Ec.31 en Ec.28 se obtiene:

𝜏1 𝜌′ (𝑡) + 𝜏2 ℎ′ (𝑡) = 𝐾𝑃1 [𝐹̅ + 𝒗𝒑 ̅̅̅̅(𝒕)𝑷 − ̅̅̅

̅̅̅̅(𝒕)𝑭 − 𝐿̅𝑃 − 𝒗𝒑 ̅( 𝒕 )
𝐺2 (𝑡) ] − 𝒉 𝐸𝑐. 32

• Aplicando Transformada de Laplace a la Ec.22 se obtiene:

𝜏1 [𝑠𝜌̅(𝑠) − 𝜌̅(0)] + 𝜏2 [𝑠ℎ̅(𝑠) − ℎ̅(0)] = 𝐾𝑃1 [𝐹̅ + 𝒗𝒑 ̅̅̅(𝒕)𝑷 − 𝐺̅2 (𝑠)] − 𝒉̅ (𝒔) 𝐸𝑐. 33
̅̅̅(𝒔)𝑭 − 𝐿̅𝑃 − 𝒗𝒑

̅ (𝒔) + 𝒉
̅(𝒔) + 𝜏2 𝑠𝒉
𝜏1 𝑠𝝆 ̅ (𝒔) = 𝐾𝑃1 [𝐹̅ + 𝒗𝒑
̅̅̅̅(𝒔)𝑭 − 𝐿̅𝑃 − 𝒗𝒑 ̅̅̅2 (𝑠)]
̅̅̅̅(𝒕)𝑷 − 𝐺 𝐸𝑐. 34

̅ (𝒔) = 𝐾𝑃1 [𝐹̅ + 𝒗𝒑

(𝜏2 𝑠 + 1 ) 𝒉 ̅̅̅̅(𝒔)𝑭 − 𝐿̅𝑃 − 𝒗𝒑
̅̅̅̅(𝒔)𝑷 − ̅̅̅ ̅(𝒔)
𝐺2 (𝑠)] − 𝜏1 𝑠𝝆 𝐸𝑐. 35

𝐾𝑃1 𝜏1 𝑠
̅ (𝒔) =
𝒉 [𝐹̅ + 𝒗𝒑 ̅̅̅̅(𝒔)𝑷 − ̅̅̅
̅̅̅̅(𝒔)𝑭 − 𝐿̅𝑃 − 𝒗𝒑 𝐺2 (𝑠)] − ̅ (𝒔)
𝝆 𝐸𝑐. 36
(𝜏2 𝑠 + 1 ) (𝜏2 𝑠 + 1 )
Balance de Materia de solutos

En Estacionario

𝑑(𝑚̇𝑠 ∗ 𝑥𝑃 𝑠 )
= 𝐹𝑠 ∗ 𝑥𝐹 𝑠 − 𝐿𝑃𝑠 ∗ 𝑥𝑃 𝑠 𝐸𝑐. 37
𝑑𝑡

𝑑(𝑚̇𝑠 ∗ 𝑥𝑝 𝑠 ) 𝑑[(𝐴ℎ𝑠 )𝑥𝑝 𝑠 ]

= = 𝐴[(𝜌′ 𝑠 ∗ ℎ𝑠 ) + (𝜌𝑠 ∗ ℎ′ 𝑠 )] + 𝑚 𝑇 𝑥𝑝′ 𝑠 𝐸𝑐. 38
𝑑𝑡 𝑑𝑡

• Sustituyendo Ec.38 en Ec.37 se obtiene:

𝐴[(𝜌′ 𝑠 ∗ ℎ𝑠 ) + (𝜌𝑠 ∗ ℎ′ 𝑠 )] + 𝑚 𝑇 𝑥𝑝′ 𝑠 = 𝐹 ∗ 𝑥𝐹 𝑠 − 𝐿𝑃 ∗ 𝑥𝑃 𝑠 𝐸𝑐. 39

Estado no estacionario

• Aplicando la Ec.40 para el balance en estado no estacionario, se obtiene:

𝐴[(𝜌′ (𝑡) ∗ ℎ𝑠 ) + (𝜌𝑠 ∗ ℎ′ (𝑡))] + 𝑚 𝑇 𝑥𝑝′ (𝑡) = 𝐹 ∗ 𝑥𝐹 (𝑡) − 𝐿𝑃 ∗ 𝑥𝑃 (𝑡) 𝐸𝑐. 40

• Sustituyendo Ec.17 y Ec.22 en Ec.40

(𝐴 ∗ ℎ𝑠 )𝜌′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ′ (𝑡) + 𝑚 𝑇 𝑥𝑝′ (𝑡) = 𝐹𝑥𝐹 (𝑡)−𝐿𝑃 𝑥𝑃 (𝑡) 𝐸𝑐. 41

• Al sustraer la ecuación de estado estacionario Ec.37 en la Ec.41 se obtiene:

(𝐴 ∗ ℎ𝑠 )𝜌̅ ′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ̅′ (𝑡) + 𝑚 𝑇 𝑥̅ 𝑝′ (𝑡) = 𝐹𝑥𝐹 (𝑡)−𝐿𝑃 𝑥𝑃 (𝑡) 𝐸𝑐. 42

• Sustituyendo las variables de desviación en la Ec.42

(𝐴 ∗ ℎ𝑠 )𝜌̅ ′ (𝑡) + (𝐴 ∗ 𝜌𝑠 )ℎ̅′ (𝑡) + 𝑚 𝑇 𝑥̅ 𝑝′ (𝑡) = 𝐹𝑥̅ 𝐹 (𝑡)−𝐿𝑃 𝑥̅𝑃 (𝑡) 𝐸𝑐. 43
 • Dividiendo la Ec.43 entre [𝐿𝑃 ] se obtiene:

(𝐴 ∗ ℎ𝑠 ) ′ (𝐴 ∗ 𝜌𝑠 ) ′ 𝑚𝑇 ′ 𝐹
𝜌̅ (𝑡) + ℎ̅ (𝑡) + 𝑥̅ 𝑝 (𝑡) = 𝑥̅𝐹 (𝑡) − 𝑥̅𝑃 (𝑡) 𝐸𝑐. 44
𝐿𝑃 𝐿𝑃 𝐿𝑃 𝐿𝑃

Donde:

(𝐴 ∗ ℎ𝑠 )
𝜏3 = 𝐸𝑐. 45
𝐿𝑃

(𝐴 ∗ 𝜌𝑠 )
𝜏4 = 𝐸𝑐. 46
𝐿𝑃

𝑚𝑇
𝜏5 = 𝐸𝑐. 47
𝐿𝑃

𝐹
𝐾𝑃2 = 𝐸𝑐. 48
𝐿𝑃

• Sustituyendo Ec.45, Ec.46, Ec.47 y Ec.48 en la Ec.44 se obtiene:

𝜏3 𝜌̅ ′ (𝑡) + 𝜏4 ℎ̅′ (𝑡) + 𝜏5 𝑥̅𝑝′ (𝑡) = 𝐾𝑃2 𝑥̅𝐹 (𝑡) − 𝑥̅𝑃 (𝑡) 𝐸𝑐. 49

• Aplicando Transformada de Laplace a la Ec.49 se obtiene:

𝜏3 [𝑠𝜌̅ (𝑠) − 𝜌̅ (0)] + 𝜏4 [𝑠ℎ̅(𝑠) − ℎ̅(0)] + 𝜏5 [𝑠𝑥̅𝑝 (𝑠) − 𝑥̅𝑝 (0)] = 𝐾𝑃2 𝑥̅𝐹 (𝑠) − 𝑥̅𝑃 (𝑠) 𝐸𝑐. 50

𝜏3 𝑠𝜌̅ (𝑠) + 𝜏4 𝑠ℎ̅(𝑠) + 𝜏5 𝑠𝑥̅𝑝 (𝑠) = 𝐾𝑃2 𝑥̅𝐹 (𝑠) − 𝑥̅𝑃 (𝑠) 𝐸𝑐. 51
 𝜏5 𝑠𝑥̅𝑝 (𝑠) + 𝑥̅𝑃 (𝑠) = 𝐾𝑃2 𝑥̅𝐹 (𝑠) − 𝜏3 𝑠𝜌̅ (𝑠) − 𝜏4 𝑠ℎ̅(𝑠) 𝐸𝑐. 52

̅𝑷 (𝒔) = 𝐾𝑃2 𝑥̅𝐹 (𝑠) − 𝜏3 𝑠𝜌̅ (𝑠) − 𝜏4 𝑠ℎ̅(𝑠)

(𝜏5 𝑠 + 1)𝒙 𝐸𝑐. 53

𝐾𝑃2 𝜏3 𝑠 𝜏4 𝑠
̅𝑷 (𝒔) =
𝒙 𝑥̅𝐹 (𝑠) − 𝜌̅ (𝑠) − ℎ̅(𝑠) 𝐸𝑐. 54
(𝜏5 𝑠 + 1) (𝜏5 𝑠 + 1) (𝜏5 𝑠 + 1)

Balance de Energía

Evaporador

Donde:

Entalpia de la solución: (1 - 0.7x) (T - TR) Btu/lb

La densidad de la solución: 62.384+23.14x+11.7x2

• Linealizando la Entalpia y La Entalpia por la Densidad, se obtiene:

𝜑𝐻
= (−0.7𝑥)(𝑇 − 𝑇𝑅 ) = 𝐶 𝐸𝑐. 55
𝜑𝑥
𝜑𝐻
= (1 − 0.7𝑥) = 𝐵 𝐸𝑐. 56
𝜑𝑇
𝐻𝑙𝑖𝑛𝑒𝑙𝑖𝑧𝑎𝑑𝑎 = (−0.7𝑥)(𝑇 − 𝑇𝑅 )(𝑥 − 𝑥𝑠 ) + (1 − 0.7𝑥𝑠 )(𝑇 − 𝑇𝑠 ) 𝐸𝑐. 57

𝐻𝑙𝑖𝑛𝑒𝑙𝑖𝑧𝑎𝑑𝑎 = 𝐶𝒙 ̅
̅ + 𝐵𝑻 𝐸𝑐. 58
(𝐻 ∗ 𝜌) = [(1 − 0.7𝑥)(𝑇 − 𝑇𝑅 )] ∗ [62.384 + 23.14𝑥 + 11.7𝑥 2 ] 𝐸𝑐. 59
(𝐻 ∗ 𝜌) = [62.38 − 20.52𝑥 − 4.49𝑥 2 − 8.19𝑥 3 ](𝑇 − 𝑇𝑅 ) 𝐸𝑐. 60
 𝜑𝐻𝜌
[ ] = (20.52 − 8.91𝑥𝑠 − 24.57𝑥𝑠 2 )(𝑇𝑠 − 𝑇𝑅 ) = 𝐷𝑠 𝐸𝑐. 61
𝜑𝑥 𝑠
𝜑𝐻𝜌
[ ] = 62.38 − 20.52𝑥𝑠 − 4.49𝑥𝑠 2 − 8.19𝑥𝑠 3 = 𝐸𝑠 𝐸𝑐. 62
𝜑𝑇 𝑠

𝐻𝜌𝑙𝑖𝑛𝑒𝑙𝑖𝑧𝑎𝑑𝑎 = 𝐷(𝑥 − 𝑥𝑠 ) + 𝐸𝑠 (𝑇 − 𝑇𝑠 ) 𝐸𝑐. 63

𝐻𝜌𝑙𝑖𝑛𝑒𝑙𝑖𝑧𝑎𝑑𝑎 = 𝐷𝒙 ̅=𝑿
̅ + 𝐸𝑠 𝑻 𝐸𝑐. 64

Estado estacionario:

𝑑(𝑋𝑠 𝑉𝑠 )
= 𝐹𝑠 ∗ 𝐻𝐹𝑠 − 𝐺2 𝑠 𝐻𝐺2 𝑠 − 𝐿𝑃 𝑠 ∗ 𝐻𝑃𝑠 + 𝑈𝐴[𝑇𝑉 𝑠 − 𝑇𝑠 ] 𝐸𝑐. 65
𝑑𝑡

𝑑(𝑋𝑠 𝑉𝑠 ) 𝑑(𝑋𝑠 ℎ𝑠 )
=𝐴 = 𝐹𝑠 𝐻𝐹𝑠 − 𝐺2 𝑠 𝐻𝐺2 𝑠 − 𝐿𝑃 𝑠 𝐻𝑃𝑠 + 𝑈𝐴[𝑇𝑉 𝑠 − 𝑇𝑠 ] 𝐸𝑐. 66
𝑑𝑡 𝑑𝑡

• Sustituyendo Ec.58 en Ec.66, se obtiene:

𝑑(𝑋𝑠 ℎ𝑠 )
𝐴 = 𝐹𝑠 (𝐶𝒙 ̅ 𝑠 ) − 𝐺2 𝑠 𝐻𝐺2 𝑠 − 𝐿𝑃 (𝐶𝒙
̅𝒔 + 𝐵𝑻 ̅𝑠 ) + 𝑈𝐴[𝑇𝑉 𝑠 − 𝑇𝑠 ]
̅𝒔 + 𝐵𝑻 𝐸𝑐. 67
𝑑𝑡 𝑠

• Derivada de un producto:

𝑑(𝑋𝑠 ℎ𝑠 )
𝐴 = 𝐴{[𝑋 ′ 𝑠 ∗ ℎ𝑠 ] + [𝑋𝑠 ∗ ℎ′ 𝑠 ]} 𝐸𝑐. 68
𝑑𝑡

• Sustituyendo Ec.68 en Ec.67, se obtiene:

𝐴{[𝑋 ′ 𝑠 ∗ ℎ𝑠 ] + [𝑋𝑠 ∗ ℎ′ 𝑠 ]} = 𝐹𝑠 ∗ 𝐻𝐹𝑠 − 𝐺2 𝑠 𝐻𝐺2 𝑠 − 𝐿𝑃 𝑠 ∗ 𝐻𝑃𝑠 + 𝑈𝐴[𝑇𝑉 𝑠 − 𝑇𝑠 ] 𝐸𝑐. 69

Estado no estacionario:

𝑑[𝑋(𝑡) ∗ 𝑉(𝑡)]
= 𝐹(𝑡)𝐻𝐹 (𝑡) − 𝐺2 (𝑡)𝐻𝐺2 (𝑡) − 𝐿𝑃 (𝑡)𝐻𝑃 (𝑡) + 𝑈𝐴[𝑇𝑉 (𝑡) − 𝑇(𝑡)] 𝐸𝑐. 70
𝑑𝑡

𝐴{[𝑋 ′ (𝑡) ∗ ℎ𝑠 ] + [𝑋𝑠 ∗ ℎ′ (𝑡)]} = 𝐹(𝑡)𝐻𝐹 (𝑡) − 𝐺2 (𝑡)𝐻𝐺2 (𝑡) − 𝐿𝑃 (𝑡) ∗ 𝐻𝑃 (𝑡) + 𝑈𝐴[𝑇𝑉 (𝑡) − 𝑇(𝑡)] 𝐸𝑐. 71

• Sustituyendo la Ec.58 en Ec.71, se obtiene:

𝐴{[𝑋 ′ (𝑡) ∗ ℎ𝑠 ] + [𝑋 ′ 𝑠 ∗ ℎ′ (𝑡)]} = {𝐹(𝑡)[𝐶𝒙 ̅ (𝑡)]} − 𝐺2 (𝑡)𝐻𝐺2 (𝑡) − {𝐿𝑝 (𝑡)[𝐶𝒙

̅(𝑡) + 𝐵𝑻 ̅ (𝑡)]} + 𝑈𝐴[𝑇𝑉 (𝑡) − 𝑇(𝑡)]
̅(𝑡) + 𝐵𝑻 𝐸𝑐. 72

• Al sustraer la ecuación de estado estacionario Ec.69 en la Ec.72 se obtiene:

𝐴{[𝑋̅ ′ (𝑡) ∗ ℎ𝑠 ] + [𝑋 ′ 𝑠 ∗ ℎ̅′ (𝑡)]} = {𝐹̅ (𝑡)[𝐶𝒙 ̅ (𝑡)]} − 𝐺̅2 (𝑡)𝐻

̅(𝑡) + 𝐵𝑻 ̅𝐺2 (𝑡) − {𝐿̅𝑝 (𝑡)[𝐶𝒙 ̅ (𝑡)]} + 𝑈𝐴[𝑇̅𝑉 (𝑡) − 𝑇̅(𝑡)]
̅(𝑡) + 𝐵𝑻 𝐸𝑐. 73

(𝐴ℎ𝑠 )𝑋̅ ′ (𝑡) + (𝐴𝑋 ′ 𝑠 )ℎ̅′ (𝑡) = {𝐹̅ (𝑡)[𝐶𝒙 ̅ (𝑡)]} − 𝐺̅2 (𝑡)𝐻
̅(𝑡) + 𝐵𝑻 ̅𝐺2 (𝑡) − {𝐿̅𝑝 (𝑡)[𝐶𝒙 ̅ (𝑡)]} + 𝑈𝐴𝑇̅𝑉 (𝑡) − 𝑈𝐴𝑇̅(𝑡)
̅(𝑡) + 𝐵𝑻 𝐸𝑐. 74

• Dividiendo la Ec.74 entre UA, se obtiene:

(𝐴ℎ𝑠 ) ′ (𝐴𝑋 ′ 𝑠 ) ′ {𝐹̅ (𝑡)[[𝐶𝒙 ̅ (𝑡)]]} 𝐺̅2 (𝑡)𝐻

̅(𝑡) + 𝐵𝑻 ̅̅̅𝑝 (𝑡)[[𝐶𝒙
̅𝐺2 (𝑡) {𝐿 ̅ (𝑡)]]}
̅(𝑡) + 𝐵𝑻
𝑈𝐴
𝑋̅ (𝑡) +
𝑈𝐴
ℎ̅ (𝑡) =
𝑈𝐴
−
𝑈𝐴
−
𝑈𝐴
+ ̅̅̅
𝑇𝑉 (𝑡) − 𝑇̅(𝑡) 𝐸𝑐. 75

Donde:

𝐴ℎ𝑠
𝜏6 = 𝐸𝑐. 76
𝑈𝐴

𝐴𝑋 ′ 𝑠
𝜏7 = 𝐸𝑐. 77
𝑈𝐴

1
𝐾𝑝3 = 𝐸𝑐. 78
𝑈𝐴
 • Sustituyendo Ec.76, Ec.77, Ec.78 en la Ec.75, se obtiene:

𝜏6 𝑋̅ ′ (𝑡) + 𝜏7 ℎ̅′ (𝑡) = 𝐾𝑝3 [{𝐹̅ (𝑡)[[𝐶𝒙 ̅ (𝑡)]]} − {𝐿

̅(𝑡) + 𝐵𝑻 ̅̅̅𝑝 (𝑡)[[𝐶𝒙 ̅ (𝑡)]]} − 𝐺̅2 (𝑡)𝐻
̅(𝑡) + 𝐵𝑻 ̅𝐺2 (𝑡)] + ̅̅̅
𝑇𝑉 (𝑡) − 𝑇̅(𝑡) 𝐸𝑐. 79

• Aplicando Transformada de Laplace a la Ec.66 se obtiene:

𝜏6[𝑠𝑋̅(𝑠) − 𝑋̅(0)] + 𝜏7[𝑠ℎ̅(𝑡) − ℎ̅(0)] = 𝐾𝑝3 [{𝐹̅ (𝑠)[[𝐶𝒙 ̅ (𝑠)]]} − {𝐿

̅(𝑠) + 𝐵𝑻 ̅̅̅𝑝 (𝑠)[[𝐶𝒙 ̅(𝑠)]]} − 𝐺̅2 (𝑠)𝐻
̅(𝑠) + 𝐵𝑻 ̅𝐺2(𝑡)] + 𝑇
̅̅̅ ̅
𝑉 (𝑠) − 𝑇 (𝑠) 𝐸𝑐. 80

𝜏6 𝑠𝑿 ̅ (𝒕) = 𝐾𝑝3 [{𝐹̅ (𝑠)[[𝐵𝒙

̅ (𝒔) + 𝜏7 𝑠𝒉 ̅ (𝑠)]]} − {𝐿
̅(𝑠) + 𝐴𝑻 ̅̅̅𝑝 (𝑠)[[𝐵𝒙 ̅ (𝑠)]]} − 𝐺̅2 (𝑠)𝐻
̅(𝑠) + 𝐴𝑻 ̅𝐺2 (𝑡)] + ̅̅̅
𝑇𝑉 (𝑠) − 𝑇̅(𝑠) 𝐸𝑐. 81

̅ (𝒔) = 𝐾𝑝3 [{𝐹̅ (𝑠)[[𝐶𝒙

𝜏6 𝑠𝑿 ̅ (𝑠)]]} − {𝐿
̅(𝑠) + 𝐵𝑻 ̅̅̅𝑝 (𝑠)[[𝐶𝒙 ̅ (𝑠)]]} − 𝐺̅2 (𝑠)𝐻
̅(𝑠) + 𝐵𝑻 ̅𝐺2 (𝑡)] + ̅̅̅ ̅ (𝒕)
𝑇𝑉 (𝑠) − 𝑇̅(𝑠) − 𝜏7 𝑠𝒉 𝐸𝑐. 82

̅ (𝒔) =
𝑿
𝐾𝑝3
𝜏6 𝑠
[{𝐹̅ (𝑠)[[𝐶𝒙 ̅ (𝑠)]]} − {𝐿
̅(𝑠) + 𝐵𝑻 ̅̅̅𝑝 (𝑠)[[𝐶𝒙 ̅𝐺2 (𝑡)] + 1 [𝑇
̅ (𝑠)]]} − 𝐺̅2 (𝑠)𝐻
̅(𝑠) + 𝐵𝑻 ̅̅̅(𝑠) − 𝑇̅(𝑠)] − 𝜏7 𝒉
𝜏 𝑠 𝑉 6 𝜏
̅ (𝒕)
6
𝐸𝑐. 83

Para los tubos calandria

Estado Estacionario

𝑑𝑇𝑉 𝑠
𝐶𝑚 = 𝐺1 𝑠 𝜆 − 𝑈𝐴[𝑇𝑉 𝑠 − 𝑇𝑠 ] 𝐸𝑐. 84
𝑑𝑡

Estado no Estacionario

𝑑𝑇𝑉 (𝑡)
𝐶𝑚 = 𝐺1 (𝑡)𝜆 − 𝑈𝐴[𝑇𝑉 (𝑡) − 𝑇(𝑡)] 𝐸𝑐. 85
𝑑𝑡

• Sustrayendo la ecuación de estado estacionario Ec.84 en la del estado no

estacionario Ec.85, se obtiene:
 𝑑𝑇̅𝑉 (𝑡)
𝐶𝑚 = 𝐺1̅ (𝑡)𝜆 − 𝑈𝐴[𝑇̅𝑉 (𝑡) − 𝑇̅(𝑡)] 𝐸𝑐. 86
𝑑𝑡

𝑑𝑇̅𝑉 (𝑡)
𝐶𝑚 = 𝐺1̅ (𝑡)𝜆 − 𝑈𝐴𝑇̅𝑉 (𝑡) + 𝑈𝐴𝑇̅(𝑡) 𝐸𝑐. 87
𝑑𝑡

• Dividiendo la Ec.87 entre UA, se obtiene:

𝐶𝑚 𝑑𝑇̅𝑉 (𝑡) 𝜆
= 𝐺1̅ (𝑡) − 𝑇̅𝑉 (𝑡) + 𝑇̅(𝑡) 𝐸𝑐. 88
𝑈𝐴 𝑑𝑡 𝑈𝐴

Donde:

𝐶𝑚
𝜏8 = 𝐸𝑐. 89
𝑈𝐴

𝜆
𝐾𝑃4 = 𝐸𝑐. 90
𝑈𝐴

• Sustituyendo Ec.89 y Ec.90 en la Ec.88, se obtiene:

𝑑𝑇̅𝑉 (𝑡)
𝜏8 = 𝐾𝑃4 𝐺1̅ (𝑡) − 𝑇̅𝑉 (𝑡) + 𝑇̅(𝑡) 𝐸𝑐. 91
𝑑𝑡

• Aplicando Transformada de Laplace a la Ec.91 se obtiene:

𝜏8 [𝑠𝑇̅𝑉 (𝑠) − 𝑇̅𝑉 (0)] = 𝐾𝑃4 𝐺1̅ (𝑡) − 𝑇̅𝑉 (𝑠) + 𝑇(𝑠) 𝐸𝑐. 92
 𝜏8 𝑠𝑇̅𝑉 (𝑠) + 𝑇̅𝑉 (𝑠) = 𝐾𝑃4 𝐺1̅ (𝑡) + 𝑇̅(𝑠) 𝐸𝑐. 93

𝐾𝑃4 1
̅̅̅
𝑻𝑽̅(𝒔) = 𝐺1̅ + 𝑇̅(𝑠) 𝐸𝑐. 94
𝜏8 𝑠 + 1 𝜏8 𝑠 + 1

Condensador

Para el vapor que entra al condensador:

Estado estacionario

𝑑𝑄𝑐𝑠
= 𝐺2 𝑠 𝐻𝐺2 𝑠 𝐸𝑐. 95
𝑑𝑡

𝑑𝑄𝑐𝑠
= 𝐺2 𝑠 [𝐶𝑃 (𝑇𝑉1 𝑠 − 𝑇𝑠 )] 𝐸𝑐. 96
𝑑𝑡

Con la ecuación de Antoine se puede obtener la relación para 𝑇𝑉1

3816.44
𝑇𝑉1 = − 227 𝐸𝑐. 97
18.3036 − 𝑙𝑛𝑃1

Estado no Estacionario

𝑑𝑄𝑐(𝑡)
= 𝐺2 (𝑡)𝐻𝐺2 (𝑡) 𝐸𝑐. 98
𝑑𝑡

𝑑𝑄𝑐(𝑡)
= 𝐺2 (𝑡)[𝐶𝑃 (𝑇𝑉1 (𝑡) − 𝑇(𝑡))] 𝐸𝑐. 99
𝑑𝑡
 • Sustrayendo la ecuación del estado estacionario Ec.96 en la del estado no
estacionario Ec.99 se obtiene:

̅̅̅̅(𝑡)
𝑑𝑄𝑐
= ̅̅̅
𝐺2 (𝑡)[𝐶𝑃 𝐺2 (𝑇̅𝑉1 (𝑡) − 𝑇̅(𝑡))] 𝐸𝑐. 100
𝑑𝑡

Agua de enfriamiento:

Estado Estacionario

𝑑𝑄𝑎𝑠
= 𝐿𝑊2 𝑠 (𝐻𝐶 𝑠 −𝐻𝑊2 𝑠 ) 𝐸𝑐. 101
𝑑𝑡

𝑑𝑄𝑎𝑠
= 𝐿𝑊2 𝑠 [𝐶𝑃𝑊2 (𝑇𝐶 𝑠 −𝑇𝑊2 𝑠 )] 𝐸𝑐. 102
𝑑𝑡

Estado no Estacionario

𝑑𝑄𝑎(𝑡)
= 𝐿𝑊2 (𝑡)[𝐶𝑃𝑊2 (𝑇𝐶 (𝑡)−𝑇𝑊2 (𝑡))] 𝐸𝑐. 103
𝑑𝑡

• Sustrayendo la ecuación del estado estacionario Ec.102 en la del estado no

estacionario Ec.103 se obtiene:

̅̅̅̅(𝑡)
𝑑𝑄𝑎
= 𝐿̅𝑊 (𝑡)[𝐶𝑃𝑊2 (𝑇̅𝐶 (𝑡)−𝑇̅𝑊2 (𝑡))] 𝐸𝑐. 104
𝑑𝑡

̅̅̅̅ (𝑡) 𝑑𝑄𝑎

𝑑𝑄𝑐 ̅̅̅̅(𝑡)
= 𝐸𝑐. 105
𝑑𝑡 𝑑𝑡

• Calor cedido = Calor absorbido, se obtiene:

 ̅̅̅
𝐺2 (𝑡)[𝐶𝑃 𝐺2 (𝑇̅𝑉1 (𝑡) − 𝑇̅(𝑡))] = 𝐿̅𝑊 (𝑡)[𝐶𝑃𝑊2 (𝑇̅𝐶 (𝑡)−𝑇̅𝑊2 (𝑡))] 𝐸𝑐. 106

̅̅̅
𝐺2 (𝑡)[𝐶𝑃 𝐺2 (𝑇̅𝑉1 (𝑡) − 𝑇̅(𝑡))]
𝐿̅𝑊 (𝑡) = 𝐸𝑐. 107
[𝐶𝑃𝑊 (𝑇̅𝐶 (𝑡)−𝑇̅𝑊 (𝑡))]

• Aplicando transformada de Laplace a la Ec.107, se obtiene:

̅̅̅
𝐺2 (𝑠)[𝐶𝑃 𝐺2 (𝑇̅𝑉1 (𝑠) − 𝑇̅(𝑠))]
𝐿̅𝑊 (𝑠) = 𝐸𝑐. 108
[𝐶𝑃𝑊2 (𝑇̅𝐶 (𝑠)−𝑇̅𝑊 (𝑠))]

• Sustituyendo Ec.97 en Ec.108, se obtiene:

̅̅̅ 3816.44
𝐺2 (𝑠) [𝐶𝑃 𝐺2 ([ − 227] − 𝑇̅(𝑠))]
18.3036 − 𝑙𝑛𝑃1
𝐿̅𝑊 (𝑠) = 𝐸𝑐. 109
[𝐶𝑃𝑊2 (𝑇̅𝐶 (𝑠)−𝑇̅𝑊 (𝑠))]

Debido a que es vapor saturado se puede obtener los valores de entalpia en las tablas.
3. Determinación del comportamiento dinámico del proceso en
lazo abierto

3.1. Función de transferencia

3.1.1. Modelo Matemático

Sensor-Trasmisor de Nivel

Tomando la Ec.36:

𝐾𝑃1 𝜏1 𝑠
̅ (𝒔) =
𝒉 [𝐹̅ + 𝒗𝒑 ̅̅̅̅(𝒔)𝑷 − ̅̅̅
̅̅̅̅(𝒔)𝑭 − 𝐿̅𝑃 − 𝒗𝒑 𝐺2 (𝑠)] − ̅ (𝒔)
𝝆
(𝜏2 𝑠 + 1 ) (𝜏2 𝑠 + 1 )

Sensor de composición

Tomando la Ec.54:

𝐾𝑃2 𝜏3 𝑠 𝜏4 𝑠
̅𝑷 (𝒔) =
𝒙 𝑥̅𝐹 (𝑠) − 𝜌̅ (𝑠) − ℎ̅(𝑠)
(𝜏5 𝑠 + 1) (𝜏5 𝑠 + 1) (𝜏5 𝑠 + 1)

3.2. Solución de las Funciones de Transferencia

Sensor-Trasmisor de Nivel

Tomando la Ec.36 para anti transformar:

𝐾𝑃1 −𝑡/𝜏 𝜏1 −𝜏𝑡

̅ (𝒕) =
𝒉 𝑒 ̅ + 𝒗𝒑
2 ∗ [𝐹 ̅̅̅̅(𝒕)𝑭 − 𝐿̅𝑃 − 𝒗𝒑
̅̅̅̅(𝒕)𝑷 − ̅̅̅
𝐺2 (𝑡)] − ̅(𝒕) 𝐸𝑐. 110
𝑒 2 ∗𝝆
𝜏2 𝜏2
Sensor de composición

Tomando la Ec.54 para anti transformar:

𝐾𝑃2 −𝑡/𝜏 𝜏3 −𝜏𝑡 𝜏4 −𝜏𝑡

̅𝑷 (𝒕) =
𝒙 𝑒 5 ∗ 𝑥̅ (𝑡) −
𝐹 𝑒 5 ∗ 𝜌̅ (𝑡) − 𝑒 5 ∗ ℎ̅(𝑡) 𝐸𝑐. 111
𝜏5 𝜏5 𝜏5

3.2.1. Diagramas de Lazo abierto

3.2.1.1. Diagrama de bloque de Lazo abierto de Nivel

Figura 1: Diagrama de bloque de lazo abierto para el nivel del Evaporador.

 Figura 2: Diagrama del bloque de salida

Figura 3: Diagrama del bloque de alimentación.

 Figura 4: Resultado de la simulación del diagrama de lazo abierto de Nivel.

3.2.1.2. Diagrama de Lazo abierto de Composición

Figura 5: Diagrama de bloque de lazo abierto para la composición a la salida del Evaporador.
 Figura 6: Resultado de la simulación del diagrama de lazo abierto de Composición.

3.2.2 Respuesta de la(s) variable(s) de control considerando un cambio del

10% en la(s) variable(s) perturbación, con respecto al valor de estado
estacionario inicial.

Se obtuvieron las siguientes respuestas con el cambio:

Lazo abierto de Nivel:

 Figura 7: Cambio con 10% de las variables perturbadoras.

Lazo abierto de Composición:

Figura 8: Cambio con 10% de las variables perturbadoras.

4. Determinación del comportamiento de la (s) variable (s) de
control para el sistema de control presentado a cada problema

Función Para el Gp:

1 1
𝐺𝑝 (𝑠) = = 𝐸𝑐. 112
(𝜏5 𝑠 + 1 ) 3.34𝑠 + 1

Controlador de nivel LC20:

➢ El sensor de nivel es un flotador con rango de 0 a 2 pies por arriba del fondo de los
tubos.

➢ El controlador de nivel (LC20) es proporcional, con un punto de control de 25% del

rango, y se ajusta para un control estricto de nivel (banda proporcional entre 1 y
10%).

𝐺𝑇1 (𝑠) = 𝐾𝑇 𝐸𝑐. 113

𝑉𝑎.𝑡𝑢𝑏𝑜𝑠 4.68 𝐹𝑡 3
ℎ𝑠𝑜𝑙.𝑡𝑢𝑏𝑜𝑠 = = = 0.015 𝐹𝑡 𝐸𝑐. 114
𝐴 309.5 𝐹𝑡 2

(0.015 + 0)
𝐾𝑇 = = 0.007 𝐸𝑐. 115
(0.015 + 2) − (0.015 + 0)

Controlador Proporcional (P) de Lc20

𝐺𝐶1 (𝑠) = 𝐾𝑐 𝐸𝑐. 116

100%
𝐾𝑐 = = 𝐸𝑐. 117
𝑃𝐵
Para la Válvula de alimentación:

➢ Las válvulas de control tienen características lineales, constantes de tiempo de 5

segundos y se dimensionan para una sobrecapacidad del 100%

𝐾𝑉1
𝐺𝑉1 (𝑠) = 𝐸𝑐. 118
𝜏𝑉 𝑠 + 1

𝜏𝑉 = 5 𝑠𝑒𝑔 𝐸𝑐. 119

𝐶𝑉𝐹𝑀𝑎𝑥 𝜌𝐻2𝑂
𝐾𝑉1 = ∗ √Δ𝑃(𝑡) ∗ 𝐸𝑐. 120
100%𝐶𝑂 ρF

𝐹
𝐶𝑉𝐹𝑀𝑎𝑥 = 200% ∗ 𝐸𝑐. 121
500 ∗ 𝑣𝑝 ∗ √Δ𝑃(𝑡) ∗ 𝐺

𝑙𝑏
509.41 𝑙𝑏
𝐶𝑉𝐹𝑀𝑎𝑥 = 200% ∗ ℎ = 0.99 𝐸𝑐. 122
500 ∗ 0.5 ∗ √16 psia ∗ 1.038 ℎ√𝑝𝑠𝑖𝑎

𝑙𝑏
𝜌𝐻2𝑂 62.4 3
= √16 𝑝𝑠𝑖𝑎 ∗ 𝐹𝑡
√Δ𝑃(𝑡) ∗ 𝑙𝑏
= 15.4√𝑝𝑠𝑖𝑎 𝐸𝑐. 123
ρF 64.81 3
𝐹𝑡

𝑙𝑏
0.99
ℎ√𝑝𝑠𝑖𝑎 𝑙𝑏
𝐾𝑉1 = ∗ 15.4√𝑝𝑠𝑖𝑎 = 0.15 𝐸𝑐. 124
100%𝐶𝑂 ℎ%𝐶𝑂

0.15
𝐺𝑉1 (𝑠) = 𝐸𝑐. 125
5𝑠 + 1
Sensor de Composición AT20:

➢ El rango del transmisor de concentración (AT20) es de 0.3 a 0.8 unidades de

fracción de masa, y la constante de tiempo, de 25 seg.

𝐾𝑇2
𝐺𝑇2 (𝑠) = 𝐸𝑐. 126
𝜏𝑣 𝑠 + 1

0.3
𝐾𝑇2 = = 0.6 𝐸𝑐. 127
0.8 − 0.3

𝐾𝑇2 0.6
𝐺𝑇2 (𝑠) = = 𝐸𝑐. 128
𝜏𝑣 𝑠 + 1 25𝑠 + 1

Controlador ARC20 PI:

1
𝐺𝐶1 (𝑠) = 𝐾𝑐 [1 + ] 𝐸𝑐. 129
𝜏𝐼 𝑠

Para la Válvula de salida:

➢ Las válvulas de control tienen características lineales, constantes de tiempo de 5 seg

y se dimensionan para una sobrecapacidad del 100%

𝐾𝑉2
𝐺𝑉2 (𝑠) = 𝐸𝑐. 130
𝜏𝑉 𝑠 + 1

𝜏𝑉 = 5 𝑠𝑒𝑔 𝐸𝑐. 131

𝐶𝑉𝐹𝑀𝑎𝑥 𝜌𝐻2𝑂
𝐾𝑉1 = ∗ √Δ𝑃(𝑡) ∗ 𝐸𝑐. 132
100%𝐶𝑂 ρF

𝐿𝑝
𝐶𝑉𝐹𝑀𝑎𝑥 = 200% ∗ 𝐸𝑐. 133
500 ∗ 𝑣𝑝 ∗ √Δ𝑃(𝑡) ∗ 𝐺
 𝑙𝑏
152.36 𝑙𝑏
𝐶𝑉𝐹𝑀𝑎𝑥 = 200% ∗ ℎ = 0.99 𝐸𝑐. 134
500 ∗ 0.5 ∗ √1.103 psia ∗ 1.347 ℎ√𝑝𝑠𝑖𝑎

𝑙𝑏
𝜌𝐻2𝑂 62.4 3
= √1.103 𝑝𝑠𝑖𝑎 ∗ 𝐹𝑡 = 0.904 𝑝𝑠𝑖𝑎
√Δ𝑃(𝑡) ∗ 𝑙𝑏
√ 𝐸𝑐. 135
ρF 84.10 3
𝐹𝑡

𝑙𝑏
0.99
ℎ√𝑝𝑠𝑖𝑎 𝑙𝑏
𝐾𝑉2 = ∗ 0.904√𝑝𝑠𝑖𝑎 = 0.0089 𝐸𝑐. 136
100%𝐶𝑂 ℎ%𝐶𝑂

0.0089
𝐺𝑉2 (𝑠) = 𝐸𝑐. 137
5𝑠 + 1
4.1. Diagrama de bloque de lazo cerrado

Lazo cerrado del Proceso

Figura 9: Diagrama de Lazo cerrado del Proceso.

 Figura 10: Respuesta Grafica del Lazo cerrado.

̅ (𝒔)
Función de transferencia del Lazo cerrado Nivel 𝒉

𝐾𝑝1 𝐾𝑝1 𝐾𝑝1 𝐾𝑝1 𝐺𝐶1 𝐺𝑉1

̅(𝒔) =
𝒉 𝐹̅ + ̅̅̅̅(𝒔)𝑭 −
𝒗𝒑 ̅𝑃 −
𝐿 ̅̅̅
𝐺2 (𝑠) − ̅
𝒉 (𝑠) 𝐸𝑐. 138
1 + 𝐺𝑜𝑙1 1 + 𝐺𝑜𝑙1 1 + 𝐺𝑜𝑙1 1 + 𝐺𝑜𝑙1 1 + 𝐺𝑜𝑙1 𝑆𝑒𝑡𝑝𝑜𝑖𝑛𝑡

̅ (𝒔)
𝒉 𝐺𝐶1 𝐺𝑉1
= 𝐸𝑐. 139
̅ 𝑆𝑒𝑡𝑝𝑜𝑖𝑛𝑡 (𝑠)
𝒉 1 + 𝐺𝑜𝑙1

𝐺𝑜𝑙1 = 𝐺𝐶1 𝐺𝑉1 𝐺𝑇1 𝐸𝑐. 140

1 + 𝐺𝑜𝑙1 = 1 + 𝐺𝐶1 𝐺𝑉1 𝐺𝑇1 𝐸𝑐. 141

̅𝒑 (𝒔)
Función de transferencia del Lazo cerrado Composición 𝒙

𝐺𝑝 𝐾𝑝2 𝜏3 𝐺𝑃 𝐺𝐶2 𝐺𝑉2 𝐺𝑃2

̅𝒑 (𝒔) =
𝒙 𝑥̅𝐹 (𝑠) − 𝜌̅ (𝑠) + 𝑥̅ (𝑠) 𝐸𝑐. 142
1 + 𝐺𝑜𝑙1 1 + 𝐺𝑜𝑙2 1 + 𝐺𝑜𝑙2 𝑝 𝑆𝑒𝑡𝑝𝑜𝑖𝑛𝑡
 ̅𝒑 (𝒔)
𝒙 𝐺𝐶2 𝐺𝑉2 𝐺𝑃2
= 𝐸𝑐. 143
𝑥̅𝑝 𝑆𝑒𝑡𝑝𝑜𝑖𝑛𝑡 (𝑠) 1 + 𝐺𝑜𝑙2

𝐺𝑜𝑙2 = 𝐺𝐶2 𝐺𝑉2 𝐺𝑇2 𝐺𝑃 𝐸𝑐. 144

1 + 𝐺𝑜𝑙2 = 1 + 𝐺𝐶2 𝐺𝑉2 𝐺𝑇2 𝐺𝑃 𝐸𝑐. 145

4.2. Rango de los Parámetros

Para Composición (Controlador P.I):

• Sustituyendo las Variables correspondientes de la Ec.140:

1 0.0089 0.6 1
1 + 𝐺𝑜𝑙2 = 1 + [𝐾𝑐 (1 + )] ∗ ∗ ∗ 𝐸𝑐. 146
𝜏𝐼 𝑠 5𝑠 + 1 25𝑠 + 1 3.34𝑠 + 1

0.0053𝐾𝑐 𝜏𝐼 𝑠 + 0.0053𝜏𝐼 𝑠
1 + 𝐺𝑜𝑙2 = 1 + 𝐸𝑐. 147
417.5𝜏𝐼 𝑠 4 + 225.2𝜏𝐼 𝑠 3 + 33.34𝜏𝐼 𝑠 2 + 𝜏𝐼 𝑠

• Para tI = 0.1
o Sustituyendo tI en la Ec.147, se obtiene:

0.00053𝐾𝑐 + 0.00053
1 + 𝐺𝑜𝑙2 = 1 + 𝐸𝑐. 148
41.75𝑠 4 + 22.52𝑠 3 + 3.33𝑠 2 + 0.1𝑠

0 = 41.75𝑠4 + 22.52𝑠3 + 3.33𝑠2 + 0.1𝑠 + 0.00053𝐾𝑐 + 1 + 0.00053 𝐸𝑐. 149

𝑎4 𝑠 4 = 41.75 ; 𝑎3 𝑠 3 = 22.52 ; 𝑎2 𝑠 2 = 3.33 ; 𝑎1 = 0.1 ; 𝑎0 = 0.00053𝐾𝑐 + 1.00053 𝐸𝑐. 150
• Utilizando el método Routh-Hurwitz

Fila Coeficientes del polinomio

1 41.75 3.33 0.00053𝐾𝑐 + 1.00053
2 22.52 0.1 0
3 𝑏1 𝑏2 0

𝑎3 ∗ 𝑎2 − 𝑎4 ∗ 𝑎1 (22.52 ∗ 3.33) − (41.75 ∗ 0.1)

𝑏1 = = 𝐸𝑐. 151
𝑎3 225.2

74.99 − 4.175 70.815

𝑏1 = = = 3.14 𝐸𝑐. 152
22.52 22.52

𝑎3 ∗ 𝑎0 − 𝑎4 (22.52 ∗ 0.00053𝐾𝑐 + 1.00053 ) − (41.75)

𝑏2 = = 𝐸𝑐. 153
𝑎3 22.52

22.53 + 0.011𝐾𝑐 − 41.75 −19.21 + 0.011𝐾𝑐

𝑏2 = = 𝐸𝑐. 154
22.52 22.52

𝑏1 ∗ 𝑎1 − 𝑎3 ∗ 𝑏2 (3.14 ∗ 0.1) − [22.52 ∗ (−19.21 + 0.011𝐾𝑐 )]

𝐶1 = = 𝐸𝑐. 155
𝑏1 ( 3.14 ∗ 22.52)

0.314 + 432.6 − 0.24𝐾𝑐 432.914 − 0.24𝐾𝑐

𝐶1 = = 𝐸𝑐. 156
70.71 70.71
• Igualando a 0 para encontrar a Kc

432.914 − 0.24𝐾𝑐
≤0 𝐸𝑐. 157
70.71
432.914 − 0.24𝐾𝑐 ≤ 0 𝐸𝑐. 158

𝐾𝑐 ≤ 1803.80 𝐸𝑐. 159

• Para tI = 50
o Sustituyendo tI en la Ec.142, se obtiene:

0.265𝐾𝑐 + 0.265
1 + 𝐺𝑜𝑙2 = 1 + 𝐸𝑐. 160
20875𝑠 4 + 11260𝑠 3 + 1667𝑠 2 + 50𝑠

0 = 20875𝑠4 + 11260𝑠3 + 1667𝑠2 + 50𝑠 + 0.265𝐾𝑐 + 0.265 + 1 𝐸𝑐. 161

𝑎4 𝑠 4 = 20875 ; 𝑎3 𝑠 3 = 11260 ; 𝑎2 𝑠 2 = 1667 ; 𝑎1 = 50 ; 𝑎0 = 0.265𝐾𝑐 + 1.265 𝐸𝑐. 162

Fila Coeficientes del polinomio

1 20875 1667 0.265𝐾𝑐 + 1.265
2 11260 50 0
3 𝑏1 𝑏2 0

𝑎3 ∗ 𝑎2 − 𝑎4 ∗ 𝑎1 (11260 ∗ 1667) − (20875 ∗ 50)

𝑏1 = = 𝐸𝑐. 163
𝑎3 11260

18770420 − 1043750 17726670

𝑏1 = = = 1574.30 𝐸𝑐. 164
11260 11260

𝑎3 ∗ 𝑎0 − 𝑎4 (11260 ∗ 0.265𝐾𝑐 + 1.265 ) − (20875)

𝑏2 = = 𝐸𝑐. 165
𝑎3 11260
 2983.9𝐾𝑐 + 14243.9 − 20875 2983.9𝐾𝑐 − 6631.1
𝑏2 = = 𝐸𝑐. 166
11260 11260

𝑏1 ∗ 𝑎1 − 𝑎3 ∗ 𝑏2 (1574.30 ∗ 50) − [11260 ∗ (2983.9𝐾𝑐 − 6631.1)]

𝐶1 = = 𝐸𝑐. 167
𝑏1 ( 50 ∗ 1574.30)

78715 + 74666186 − 33588580𝐾𝑐 74744901 − 33588580𝐾𝑐

𝐶1 = = 𝐸𝑐. 168
78715 78715

74744901 − 33588580𝐾𝑐
≤0 𝐸𝑐. 169
78715

74744901 − 33588580𝐾𝑐 ≤ 0 𝐸𝑐. 170

𝐾𝑐 ≤ 2.22 𝐸𝑐. 171

• Los rangos para Kc del controlador PI obtenidos son:

2.22 ≤ 𝐾𝑐 ≤ 1803.80
0.1 ≤ 𝜏𝐼 ≤ 50

• Utilizando los limites obtienido, para observar si el sistema estable, observando

las siguientes graficas obtenidas:
 Figura 11: Respuesta del controlador PI para un Kc=2.2 y I=50.

Figura 12: Respuesta del controlador PI para un Kc=1803.8 y I=50.

 Figura 13.a: Respuesta del controlador PI para un Kc=1803.8 y I=0.1.

Figura 13.b: Respuesta del controlador PI para un Kc=1803.8 y I=50 (acercamiento).

 Figura 14: Respuesta del controlador PI para un Kc=2.2 y I=0.1.

• Utilizando Valores dentro de los rangos de Kc y I se encontró la estabilidad en:

Figura 15: Respuesta del controlador PI para un Kc=111.5 y I=1.5.

Para Nivel (Controlador P):

• Utilizando la Ec.117:

o PB: 1 a 10%
100%
𝐾𝑐 = = 100 𝐸𝑐. 172
1%

100%
𝐾𝑐 = = 10 𝐸𝑐. 173
10%

• Los rangos para Kc del controlador P obtenidos son:

10 ≤ 𝐾𝑐 ≤ 100

• Utilizando los limites obtienido, para observar si el sistema estable, observando

las siguientes graficas obtenidas:

Figura 16: Respuesta del controlador P para un Kc=100.

 Figura 17: Respuesta del controlador P para un Kc=10.

• Utilizando Valores dentro de los rangos de Kc se encontró la estabilidad en:

Figura 18: Respuesta del controlador P para un Kc=30.

4.3. Método de sintonización del controlador Ziegler-Nichols

• PI

o Utilizando la Ec.145:

1 0.0089 0.6 1
1 + 𝐺𝑜𝑙2 = 1 + [𝐾𝑐 (1 + )] ∗ ∗ ∗ 𝐸𝑐. 174
𝜏𝐼 𝑠 5𝑠 + 1 25𝑠 + 1 3.34𝑠 + 1

• Primeramente, se utilizó el método de substitución directa para calcular la

ganancia y el período de oscilación últimos de un controlador proporcional; Donde
1
=0
𝜏𝐼 𝑠

0 = 417.5𝑠 3 + 225.2𝑠 2 + 33.34𝑠 + 1 + 0.0053𝐾𝑐 𝐸𝑐. 175

• A continuación, se substituyo 𝒔 = 𝒊𝒘𝒖 y 𝑲𝒄 = 𝑲𝒄𝒖

0 = 417.5𝑖𝑤𝑢 3 + 225.2𝑖𝑤𝑢 2 + 33.34𝑖𝑤𝑢 + 1 + 0.0053𝐾𝑐𝑢 𝐸𝑐. 176

• Para obtener la siguiente simplificación, que se muestra en siguiente sistema de

ecuaciones:

−225.2𝑤𝑢 2 + 1 + 0.0053𝐾𝑐𝑢 = 0 𝐸𝑐. 177

−417.5𝑤𝑢 3 + 33.34𝑤𝑢 = 0 𝐸𝑐. 178

• Para 𝑤𝑢 = 0

1
𝐾𝑐𝑢 = − = −188.67 %𝐶𝑂%𝑇𝑂 𝐸𝑐. 179
0.0053
 33.34 𝑟𝑎𝑑
𝑤𝑢 = √ = 0.282 𝐸𝑐. 180
417.5 𝑠

• Para 𝑤𝑢 = 0.282 𝑟𝑎𝑑/𝑚𝑖𝑛

16.69
𝐾𝑐𝑢 = = 3190.05 %𝐶𝑂%𝑇𝑂 𝐸𝑐. 181
0.0053

La primera solución corresponde a la inestabilidad monotónica causada de tener la

acción errónea en el controlador.

2𝜋
𝑇𝑢 = = 22.28 min 𝐸𝑐. 182
0.282

De acuerdo con la tabla 6-1 de (Smith & Corripio, 1991):

𝐾𝑐𝑢
𝐾𝑐 = = 1450.02 %𝐶𝑂%𝑇𝑂 𝐸𝑐. 183
2.2

𝑇𝑢
𝜏𝐼 = = 18.57 min 𝐸𝑐. 184
1.2
 Figura 19: Respuesta del controlador PI con parámetros sintonizados.

4.3.1. Respuesta de la variable de control para un cambio del 10% en la variable

perturbadora, usando los parámetros sintonizados.
 Figura 20: Respuesta del controlador PI con parámetros sintonizados y 10% de cambio en una
variable perturbadora.

4.3.2. Respuesta de la variable de control para un cambio del 5% en el set-point

Figura 19: Respuesta del controlador PI con parámetros sintonizado y 5% de cambio en el set-point.
 5. Resultados y Análisis de Resultados

Los resultados de los balances y otras ecuaciones se muestran en el siguiente acápite con su
análisis:

5.1. Acápite

Variable Valor Unidades Notas

F 509.41 lb/h Flujo de alimentación (válvula lineal
Cv=vp=50%).
Lp 152.36 lb/h Flujo de producción (válvula lineal
Cv=vp=505).
m 509.41 lb/h Flujo masico total en el evaporador = al flujo
máximo de alimentación en el evaporador.
G1 5980.1 lb/h
G2 357.05 lb/h
Lw 1.423 lb/h
Lz 5980.1 lb/h Flujo de condensado de los tubos calandria = al
flujo de vapor que entra G1=Lz
Lc 358.47 lb/h
GF 1.038 -
GP 1.347 -
xF 0.1 -
xP 0.7 -
HF 12.09 BTU/lb
HP 14.28 BTU/lb
H1 1156.41 BTU/lb
H2 1107.11 BTU/lb
𝛒𝐅 64.81 lb/Ft3
𝛒𝐩 84.10 lb/Ft3
 𝑪𝒑𝒘 1 BTU/lb oF
𝑪𝒑𝑮𝟐 0.46 BTU/lb oF
𝚫𝐏𝐯𝐚𝐥𝐅 16 psia
𝚫𝐏𝐯𝐚𝐥𝐏 1.103 psia
PV 20 psia
P1 57.0446 mmHg La presión del evaporador en mmHg
PEvap 1.103 psia La presión a la Tc de vapor de agua = la presión
del evaporador
T 105 °F Temperatura a la presión del vapor en el
evaporador
TR 77 °F
Tv 227.92 °F Temperatura de Vapor saturado a 20 psia
Tv1 105.13 °F
Tw 90 °F
Tc 105 °F
A 309.5 Ft2
Vevap 6.05 Ft3 Volumen de solución en el evaporador con un
flujo de alimentación 509.41 y densidad de
84.10
Va. tubos 4.68 Ft3
h 1.44 Ft Sí su capacidad es de 4.2 pies3 de solución por
pie de nivel, entonces para un volumen de 6.05
Ft3
CvF 50 % Debido a que son válvulas de características
lineales Cv=vp; El flujo variara con respecto al
CvP 50 % porcentaje de apertura; Se tomo una abertura en
el rango de control de válvulas de 15% al 95%

vpF 50 %
vpP 50 %
U 150 BTU/hFt2-°F
 Cm 154 BTU/°F
𝛌 960.127 BTU/lb Calor latente tomado a una presión de vapor
saturado de 20 psia.

valv 5 s

Rflot 0-2 Ft

flot 25 s

RAT20 0.3-0.8 unidades de

fracción de masa
KP1 0.038
KP2 3.34
 17.31
 780.49
 2.92
 131.65
 3.34

F: Flujo masico en la alimentación lb/h.

Lp: Flujo masico en la salida lb/h.

m: Flujo masico total del evaporador lb/h.

G1: Flujo masico de vapor en la alimentación del evaporador lb/h.

G1: Flujo masico de vapor en la salida del evaporador lb/h.

Lz: Flujo de consensado a la salida del evaporador lb/h.

LW: Flujo de agua de enfriamiento lb/h.

GF: Gravedad específica del líquido que fluye a través de la válvula en la alimentación.

GP: Gravedad específica del líquido que fluye a través de la válvula en la salida.

𝑥𝐹 : Fraccion másica de la solución en la alimentación del evaporador.

𝑥𝑃 : Fracción másica de la solución en el producto a la salida de evaporador.

HF: Entalpia del flujo de solución en la alimentación del evaporador Btu/lb- oF.

Hp: Entalpia del flujo de solución en la salida del evaporador Btu/lb- oF.

H1: Entalpia del flujo del vapor en la entrada del evaporador Btu/lb- oF.

H2: Entalpia del flujo del vapor en la salida del evaporador Btu/lb- oF.

HW1: Entalpia del flujo del condensado en la salida del evaporador Btu/lb- oF.

HW2: Entalpia del flujo del condensado en la salida del condensador Btu/lb- oF.

𝜌𝐹 : Densidad de la solución en la alimentación lb/Ft3.

𝜌𝑝 : Densidad de la solución en la salida lb/Ft3.

PvalvF: Caída de presión de la válvula en alimentación psi.

PvalvP: Caída de presión de la válvula en la salida psia.

PV: Presión de vapor en la alimentación del evaporador psia.

P1: Presión de vapor en la entrada del condensador mmHg.

PEvap: Presión del evaporador psia.

T: Temperatura del evaporador oF.

𝑇𝑉 : Temperatura del vapor que se alimenta en los tubos calandria oF.

𝑇𝑉1 : Temperatura del vapor que entra al condensador oF.

Tw: Temperatura del agua de enfriamiento que entra al condensador oF.

Tc: Temperatura del condensado en la salida del condensador oF.

Ts=Tz: Temperatura del condensado en la salida del evaporador oF.

A: Área del evaporador Ft2.

h: Nivel del tanque Ft.

𝐶𝑉 = coeficiente de la válvula.

𝑣𝑝(𝑡)𝐹 = posición de la válvula en la alimentación. Este término representa la fracción de

apertura de la válvula.

𝑣𝑝(𝑡)𝑃 = posición de la válvula en la salida.

U: Coeficiente global de transferencia de calor Btu/h*Ft2- oF.

𝐶𝑚 : Capacidad calórica de los tubos de calandria BTU/oF.

𝑄𝑎: Calor absorbido BTU.

𝑄𝑐: Calor cedido BTU.

6. Conclusión

En síntesis, los procesos industriales pueden ser afectados por varios factores, ya sean por las
características de los equipos, materiales, mecanismo y ambiente donde se desarrolla el
proceso. Por ello, Mantener el funcionamiento del proceso en sus valores óptimos es trabajo
del ingeniero; donde deberá desarrollar una fuente de control de estas variables, de manera
que afecten en lo mínimo el proceso, así mismo, mantener el objetivo de control en los puntos
deseados.

En cuanto al problema planteado, el controlar el nivel y concentración de la solución en el

evaporador (Objetivo de control y variables de control), los cuales son casos comunes de
control en la industria; Al analizarlo se encontró las variables que causan efectos en estas,
dichas variables se dividen en las manipulables y perturbadoras, las cuales fueron
identificadas en la resolución del planteamiento y desarrollo del problema.

Seguidamente, se obtuvieron las ecuaciones donde se puede observar la dependencia que

poseen las variables de control con respecto a las manipuladas y perturbadoras, de igual
manera, se realizó los diagramas de lazo abierto y cerrado para el proceso en estudio, con
ello, se logró obtener la ecuación característica que maneja y manipula el proceso.
Finalmente, se trabajo con los controladores que contenía el problema, uno proporcional para
el control de nivel y uno proporcional-integral para el control de composición, estos
controladores fueron estabilizados con el método de Ruth, así mismo se aplican la estabilidad
con el requerido en el problema, los resultados fueron mostrados en las figuras de manera
gráfica.

7. Bibliografía

Barders, A. V. (1986). Operaciones Unitarias para la Industria. Mexico, D.F: Editorial

Limusa, S.A de C.V.

Seborg, D. E., Edgar, T. F., & Mellichamp, D. A. (2011). Process Dynmics and Control.
USA: Jhon Wiley and Sons, Inc.

Smith, C. A., & Corripio, A. B. (1991). Control Automático de Procesos Teoría y Práctica.
Mexico, DF: Editorial Limusa, S.A de C.V.

8. Apéndice

Tabla utilizada para el método de sintonización de Zielgle-Nicoles