AUTOMATIZACION Y CONTROL EN UN BIORREACTOR

AEROBICO

CONTROL POR RETROALIMENTACION

El esquema de control que se muestra en la figura se conoce como control por retroalimentación, también
se le llama circuito de control por retroalimentación. Esta técnica la aplicó por primera vez James Watt
hace casi 200’ años, para controlar un proceso industrial; consistía en mantener constante la velocidad de
una máquina de vapor con carga variable; se trataba de una aplicación del control regulador. En ese
procedimiento se toma la variable controlada y se retroalimenta al controlador para que este pueda tomar
una decisión.

TRANSDUCTOR I/P

BIORREACTOR AEROBICO

Medidores de temperatura

La temperatura es una de las principales variables que afectan el curso de los

procesos químicos, por tal razón esta variable debe ser medida con la mayor
exactitud posible para poder controlarla adecuadamente. Dentro de los principales
instrumentos que se utilizan para la medición de temperatura se tiene:
Termopares. Se basan en el hecho de que una corriente fluye en un circuito
continuo de dos alambres metálicos diferentes, si las dos junturas, o uniones, se
encuentran a diferentes temperaturas.

Termómetros de resistencia. Se basan en el hecho de que los metales cambian

su resistencia eléctrica cuando se someten a un cambio de temperaturas.

Termómetros llenos. Los termómetros de sistema lleno se diseñan para

proporcionar una indicación de la temperatura a cierta distancia del punto de
medición. El elemento sensible o medición (bulbo o ampolla) tiene un gas o un
líquido que cambia de volumen, presión o presión de vapor con la temperatura. Este
cambio se comunica por medio de un tubo capilar al Tubo de Bourdon u otro
dispositivo sensible a la presión y el volumen.

Estos dispositivos debido a su simplicidad se utilizan con frecuencia en los procesos

industriales.

Termómetros bimetálicos. El bimetal termostático se define como un material

compuesto que consta de tiras de dos ó más metales unidos entre sí. Debido a los
diferentes índices de expansión de sus componentes, Esta composición tiende a
cambiar de curvatura cuando se somete a una variación de temperatura. Los
termostatos bimetálicos se destinan a utilizarse para temperaturas que oscilan entre
1000º F y –300º F e incluso a niveles inferiores.

Termómetros de líquido en capilares de vidrio. Las tres formas de termómetros

de líquido en capilares de vidrio son:

1. Los totalmente hechos de vidrio (de cuello grabado o de escala cerrada).

2. De tubo y escala.
3. Industriales.
Estos termómetros no se utilizan en sistemas de control automático pero si se
utilizan profundamente como dispositivo de medición para el control manual y en
laboratorios de control.
Pirómetros. “Pirometría de Radiación”, es la determinación de la temperatura de
un objeto por medio de la cantidad y la naturaleza de la energía que irradia. Estos
dispositivos se clasifican en:

1. Pirómetros ópticos; basados en la brillantez de un objeto caliente.

2. Pirómetros de radiación; que miden el índice de emisión de energía por unidad
de área

VARIABLES
1.2.3 Variables de Procesos

Para sistematizar el estudio dinámico de procesos, las variables se clasifican en las

siguientes categorías:

1.2.4 Variables de perturbación

Variable que ocasiona que la variable de control se desvíe del punto de control.

1.2.5 Variables manipulación

Es la variable que se utiliza para mantener a la variable controlada en el punto de

control (punto de fijación o de régimen).
1.2.6 Variables de control

Variable que se debe mantener o controlar dentro de algún valor deseado.

u : vector de agentes de control

d : vector de perturbaciones

x : vector de variables de estado

y : vector de variables de salida

VARIABLES A CONTROLAR: T, h, S,

VARIABLES DE PERTURBACION: Qf,

VARIABLES A MANIPULAR: V, Qe
 BALANCE DE MATERIA
CONTROL DE NIVEL

q  k h ......(1)
dh
A  qe  q .......(2)
dt
(1)  (2)
dh
A  qe  k h ........(3)
dt
Linealizando (3)
k
k h k h ( h  h)
2 h
q  qe .......( Estado  estacionario)

ECUACION
dh k
A  qe  ( qe  ( h  h))
dt 2 h

BALANCE DE MATERIA DE SUSTRATO

BALANCE GENERAL:

𝑑 (𝑉. 𝑆)
= 𝐹𝑜 𝑆𝑜 − 𝜌 𝑆𝑆 + 𝑉. 𝑟. 𝑆
𝑑𝑡

𝑑 𝑉𝑆
= 𝐹𝑜 𝑆𝑜 − 𝐹𝑆𝑆 + 𝑉. 𝑟
𝑑𝑡

𝑑 𝑉𝑆 𝑑𝑆 𝑑𝑉
=𝑉 +𝑆
𝑑𝑡 𝑑𝑡 𝑑𝑡

𝑑𝑉
= 𝐹𝑜 − 𝐹
𝑑𝑡
 𝑑𝑆
𝑉 + 𝑆(𝐹𝑜 − 𝐹) = 𝐹𝑜 𝑆𝑜 − 𝐹𝑆𝑆 − 𝑉. 𝑟
𝑑𝑡

Consideramos que: 𝑟 = −𝑆𝐾

𝑑𝑆
𝑉 + 𝑆(𝐹𝑜 − 𝐹) = 𝐹𝑜 𝑆𝑜 − 𝐹𝑆𝑆 + 𝑉. 𝑆𝐾
𝑑𝑡

𝑑𝑆
𝑉 = 𝐹𝑜 (𝑆𝑜 − 𝑆) + 𝑉. 𝑆𝐾
𝑑𝑡

𝑑𝑆 𝐹𝑜
= (𝑆 − 𝑆) + 𝑉. 𝑆𝐾 … 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 𝐺𝑒𝑛𝑒𝑟𝑎𝑙
𝑑𝑡 𝑉 𝑜

Linealización:

𝑑𝑆
𝑉 = 𝐹𝑜 (𝑆𝑜 − 𝑆) − 𝑉 2 . 𝑆𝐾
𝑑𝑡

𝑑𝑆
𝑉 = 𝐹𝑜 𝑆𝑜 − 𝐹𝑜 𝑆 − 𝑉 2 . 𝑆𝐾 … 𝑒𝑐𝑢𝑎𝑐𝑖ó𝑛 1
𝑑𝑡

𝑑𝑆
𝑉 = 𝐹𝑜 𝑆𝑜̅ − 𝐹̅𝑜 (𝑆𝑜 − 𝑆𝑜̅ ) + 𝑆𝑜̅ (𝐹𝑜 − 𝐹̅𝑜 ) − [𝐹𝑜 𝑆̅ + 𝐹̅𝑜 (𝑆 − 𝑆̅) + 𝑆̅(𝐹𝑜 + 𝐹̅𝑜 )] −
𝑑𝑡
̅ 𝑆̅ + 𝐾
[𝑉 2 𝐾 ̅ (𝑆 − 𝑆̅) + 𝑆̅(𝐾 − 𝐾
̅ )]

𝑑𝑆
𝑉 𝑑 𝑡 = 𝐹𝑜 𝑆𝑜̅ − 𝐹̅𝑜 𝑆𝑜 − 𝐹̅𝑜 𝑆𝑜̅ + 𝑆𝑜̅ 𝐹𝑜 − 𝑆𝑜̅ 𝐹̅𝑜 − 𝐹𝑜 𝑆̅ − 𝐹̅𝑜 𝑆 + 𝐹̅𝑜 𝑆̅ − 𝑆̅𝐹𝑜 − 𝑆̅𝐹̅𝑜 − 𝑉 2 𝐾
̅ 𝑆̅ −

𝐾 ̅ 𝑆̅ − 𝑆̅𝐾 + 𝑆̅𝐾
̅𝑆 + 𝐾 ̅

𝑑𝑆
𝑉 𝑑 𝑡 = − 𝐹̅𝑜 𝑆𝑜 − 𝐹𝑜 𝑆̅ − 𝐹̅𝑜 𝑆 − 𝑆̅𝐹𝑜 − [𝑉 2 𝐾
̅ 𝑆̅ + 𝐾
̅ 𝑆 − 𝑆̅𝐾
̅]

𝑑𝑆
𝑉 = − 𝐹̅𝑜 𝑆𝑜 − 𝐹𝑜 𝑆̅ − 𝐹̅𝑜 𝑆 − 𝑆̅𝐹𝑜 − 𝑉 2 𝐾
̅ 𝑆̅ − 𝐾
̅ 𝑆 + 𝑆̅𝐾
̅ … 𝑒𝑐𝑢𝑎𝑐𝑖ó𝑛 2
𝑑𝑡

𝑑𝑆
𝑉 = 𝐹̅𝑜 (𝑆𝑜 − 𝑆𝑜̅ ) + 𝐹𝑜 (𝑆𝑜̅ − 𝑆̅) + 𝐹𝑜 (𝑆𝑜 − 𝑆) − 𝐾
̅ 𝑆̅(𝑉 2 − 1) − 𝐾
̅𝑆
𝑑𝑡

Sumamos la ecuación 1 y la ecuación 2:

𝑑 𝑆′
𝑉 = − 𝐹𝑜′ 𝑆𝑜 − 𝐹𝑜 𝑆 ′ − 𝐹̅𝑜 𝑆𝑜̅ + 𝑆̅𝐹𝑜 − 𝐹𝑜 𝑆 − 𝑉 2 𝐾
̅ ′𝑆̅′ − 𝐾
̅ 𝑆 + 𝑆̅𝐾
̅
𝑑𝑡
 𝑑 𝑆′
𝑉 = − 𝐹𝑜′ 𝑆𝑜 − 𝐹𝑜 𝑆 ′ + 𝑆𝑜̅ (𝐹𝑜 − 𝐹̅𝑜 ) + 𝐹̅𝑜 𝑆𝑜̅ + 𝐹𝑜 (𝑆̅ − 𝑆) − 𝑉 2 𝐾
̅ ′𝑆̅′ − 𝐾
̅ (𝑆 + 𝑆̅)
𝑑𝑡
𝑑 𝑆′
{𝑉 +𝐹𝑜 𝑆 ′ =− 𝐹𝑜′ 𝑆𝑜 +𝑆𝑜̅ (𝐹𝑜 −𝐹̅𝑜 )+𝐹𝑜 (𝑆̅−𝑆)−𝑉 2 𝐾
̅ ′𝑆̅′−𝐾
̅ (𝑆+𝑆̅)}
𝑑𝑡
𝐹̅𝑜

𝑉 𝑑 𝑆′ 𝑆𝑜 ′ 𝑆𝑜̅ 𝐹𝑜 𝑉 2 ′𝑆̅′ 𝐾̅
+ 𝑆′ = 𝐹𝑜 + (𝐹𝑜 − 𝐹̅𝑜 ) + (𝑆̅ − 𝑆) − ̅ + (𝑆̅ − 𝑆)
𝐾
𝐹̅𝑜 𝑑 𝑡 𝐹̅𝑜 𝐹̅𝑜 𝐹̅𝑜 𝐹̅𝑜 𝐹̅𝑜

𝜏 𝐾𝑆 , 𝐹𝑜

𝑑 𝑆′ 𝑆 ̅ ̅)
(𝐹 +𝐾 𝑉 2
𝜏 + 𝑆 ′ = 𝐾𝑆 , 𝐹𝑜′ + 𝐹̅𝑜 (𝐹𝑜 − 𝐹̅𝑜 ) + 𝑜𝐹̅ (𝑆̅ − 𝑆) − 𝐹̅ 𝑆′𝐾′
𝑑𝑡 𝑜 𝑜 𝑜

𝐾𝑆 , 𝑆𝑜 𝐾𝑆 , 𝐾

𝑑 𝑆′ (𝐹𝑜 − 𝐾̅)
𝜏 + 𝑆 ′ = (𝐾𝑆 , 𝐹𝑜 ) + 𝑆𝑜̅ (𝐾𝑆 , 𝑆𝑜 ) + (𝑆̅ − 𝑆) − (𝐾𝑆 , 𝐾)𝐾′
𝑑𝑡 𝐹̅𝑜

Ecuación de transferencia:

𝐹𝑜 𝑆𝑜 V 𝐹, 𝑆

𝑑𝑆
𝑉 = 𝐹𝑜 (𝑆𝑜 − 𝑆) + 𝑉 2 . 𝑆𝐾
𝑑𝑡

𝐹𝑜
𝑆. 𝑆(𝑆) = (𝑆 − 𝑆(𝑆) ) + 𝑉 2 . 𝑆(𝑆) 𝐾
𝑉 𝑜 (𝑆)

𝐹𝑜 𝐹𝑜
𝑆. 𝑆(𝑆) − 𝑉 2 . 𝑆(𝑆) 𝐾 = 𝑆𝑜 (𝑆) − 𝑆(𝑆)
𝑉 𝑉

𝐹𝑜 𝐹𝑜
𝑆. 𝑆(𝑆) − 𝑉 2 . 𝑆(𝑆) 𝐾 + 𝑆(𝑆) = 𝑆
𝑉 𝑉 𝑜 (𝑆)

𝐹𝑜 𝐹𝑜
𝑆(𝑆) [𝑆 − [𝑉𝐾 + ]] = 𝑆
𝑉 𝑉 𝑜 (𝑆)

𝐹𝑜
𝑆(𝑆) 𝑉
=
𝑆𝑜 (𝑆) 𝐹
[𝑆 − [𝑉𝐾 + 𝑉𝑜 ]]
Igualación de la ecuación general:

𝑆(𝑆)
= 𝐺(𝑆)
𝑆𝑜 (𝑆)

𝐹𝑜
𝑉
𝐺(𝑆) = 𝐹
𝑆𝑉−𝑉 2 𝐾+ 𝑜 𝑉
𝑉

𝐹𝑜
𝐺(𝑆) =
𝑆𝑉 − 𝑉 2 𝐾 + 𝐹𝑜

𝐹𝑜
(𝐺(𝑆) = )
𝑆𝑉 − 𝑉 2 𝐾 + 𝐹𝑜
𝐹𝑜

1
𝐺(𝑆) =
𝑉
(𝑆
𝐹𝑜 − 𝑉𝐾) − 1

1
𝐺(𝑆) = … 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑎 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎
𝑉
(𝑉𝐾 − 𝑆) + 1
𝐹𝑜

BALANCE DE ENERGIA

Representación gráfica del BIORREACTOR

Sabemos:

 
Velocidad de 
   Flujo   Flujo  Velocidad  Velocidad
q, T, ρ


 acumulación    que    que    de    de 
     sale   generación   consumo 
 en el sistema   entra       
 
 

Entonces:
𝜕𝑇
𝐹𝜌𝐴 𝐶𝑃 − 𝑉(∆𝐻𝑟 )(𝑟̅𝐴 + 𝐶1 (𝑇 − 𝑇̅) + 𝐶2 (𝐶𝐴 − ̅̅̅
𝐶𝐴 )) − 𝑈𝐴[𝑇 − 𝑇𝑅0 ] − 𝐹𝜌𝐴 𝐶𝑃 𝑇 = 𝑉𝐶 𝜌𝐶 𝐶𝑉
𝜕𝑡
 𝜕𝑇
𝐹𝜌𝐴 𝐶𝑃 𝑇0 − 𝑉(∆𝐻𝑟 )(𝑟̅𝐴 + 𝐶1 (𝑇 − 𝑇̅) + 𝐶2 (𝐶𝐴 − ̅̅̅
𝐶𝐴 )) − 𝑈𝐴𝑇 + 𝑈𝐴𝑇𝑅0 − 𝐹𝜌𝐴 𝐶𝑃 𝑇 = 𝑉𝐶 𝜌𝐶 𝐶𝑉
𝜕𝑡
𝜕𝑇
𝐹𝜌𝐴 𝐶𝑃 𝑇0 ′ − 𝑉(∆𝐻𝑟 )(𝑟̅𝐴 + 𝐶1 (𝑇 − 𝑇̅) + 𝐶2 (𝐶𝐴 − ̅̅̅
𝐶𝐴 )) − 𝑈𝐴𝑇 + 𝑈𝐴𝑇𝑅0 − 𝐹𝜌𝐴 𝐶𝑃 𝑇 = 𝑉𝐶 𝜌𝐶 𝐶𝑉
𝜕𝑡
𝜕𝑇′
𝐹𝜌𝐴 𝐶𝑃 𝑇0′ − 𝑉𝐶1 (∆𝐻𝑟 ) − 𝑉𝐶2 𝐶𝐴′ (∆𝐻𝑟 ) − 𝑈𝐴𝑇 ′ − 𝐹𝜌𝐴 𝐶𝑃 𝑇′ = 𝑉𝐶 𝜌𝐶 𝐶𝑉
𝜕𝑡
𝜕𝑇′
𝐹𝜌𝐴 𝐶𝑃 𝑇0 ′ − (𝑉𝐶1 (∆𝐻𝑟 ) + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃 )𝑇′ − 𝑉𝐶2 𝐶𝐴′ (∆𝐻𝑟 ) = 𝑉𝐶 𝜌𝐶 𝐶𝑉
𝜕𝑡
𝐹𝜌𝐴 𝐶𝑃 𝑉𝐶2 (∆𝐻𝑟 )
𝑇𝐴0 ′ − 𝑇′ − 𝐶′
𝑉𝐶1 (∆𝐻𝑟 ) + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃 𝑉𝐶1 (∆𝐻𝑟 ) + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃 𝐴
𝑉𝐶 𝜌𝐶 𝐶𝑉 𝜕𝑇′
=
𝑉𝐶1 (∆𝐻 )
𝑟 + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃 𝜕𝑡

𝜕𝑇′
𝜏2 + 𝑇′ = 𝐾3 𝑇𝐴0 ′ − 𝐾4 𝐶𝐴′
𝜕𝑡
𝑉𝐶 𝜌𝐶 𝐶𝑉
𝜏2 =
𝑉𝐶1 (∆𝐻𝑟 ) + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃

𝐹𝜌𝐴 𝐶𝑃
𝐾3 =
𝑉𝐶1 (∆𝐻𝑟 ) + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃

𝑉𝐶2 (∆𝐻𝑟 )
𝐾4 =
𝑉𝐶1 (∆𝐻𝑟 ) + 𝑈𝐴 + 𝐹𝜌𝐴 𝐶𝑃

−7782,35 ̅̅̅
𝐶𝐴 −7781.35
𝑟𝐴 = 𝑟̅𝐴 + 1,350 × 1017 . 𝑒 𝑇 . 2 (𝑇 − 𝑇̅) + 1,3315 × 1013 . 𝑒 𝑇 (𝐶𝐴 − ̅̅̅
𝐶𝐴 )
𝑇̅

BALANCE DE ENERGIA: Se realizara en el fluido frio

Suposiciones

a) La densidad de flujo es constante para ambos fluidos.

b) El flujo del fluido frio es contante.
c) Solo se manipulara el flujo de fluido caliente (vapor).
d) Variable controlable es la temperatura de salida de la mezcla.
e) La capacidad calorífica del fluido frio es constante.
f) Suponemos un estado pseudo-estacionario.

Balance de energía
 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝐶𝑎𝑙𝑜𝑟
[𝑑𝑒 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖𝑜𝑛] = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎] − [ 𝑑𝑒 𝑠𝑎𝑎𝑙𝑖𝑑𝑎 ] + [ 𝑐𝑒𝑑𝑖𝑑𝑜 ]
𝑑𝑒 𝑒𝑛𝑒𝑟𝑔𝑖𝑎 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔𝑖𝑎 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔𝑖𝑎 𝑜 𝑔𝑎𝑛𝑎𝑑𝑜

Entonces aplicando el principio anterior a nuestro volumen de control

𝑑𝑇𝑓𝑠 𝑑𝑇𝑓𝑠
𝑚𝑓 𝑐𝑝𝑓 + 𝑚𝑖𝑐 𝑐𝑝𝑖𝑐 = 𝐹𝑐𝑝𝑓 𝑇𝑓𝑒 − 𝐹𝑐𝑝𝑓 𝑇𝑓𝑠 + 𝜆𝑤 (1)
𝑑𝑡 𝑑𝑡
Donde

𝑚𝑓 = 𝑚𝑎𝑠𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟𝑖𝑜 (𝑚𝑒𝑧𝑐𝑙𝑎 𝑒𝑡𝑎𝑛𝑜𝑙 − 𝑎𝑔𝑢𝑎)

𝑚𝑖𝑐 = 𝑚𝑎𝑠𝑎 𝑑𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜𝑟 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟

𝑐𝑝𝑓 = 𝑐𝑎𝑙𝑜𝑟 𝑒𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐𝑜 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟𝑖𝑜

𝑐𝑝𝑖𝑐 = 𝑐𝑎𝑙𝑜𝑟 𝑒𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐𝑜 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 𝑑𝑒𝑙 𝑚𝑎𝑡𝑒𝑟𝑖𝑎𝑙 𝑑𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜𝑟 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟

𝐹 = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝑚𝑎𝑠𝑖𝑐𝑜 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟𝑖𝑜

𝑇𝑓𝑒 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑒𝑛𝑡𝑟𝑒𝑑𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟𝑖𝑜 𝑑𝑒 𝑙𝑎 𝑚𝑒𝑧𝑐𝑙𝑎

𝑇𝑓𝑠 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟𝑖𝑜 𝑑𝑒 𝑙𝑎 𝑚𝑒𝑧𝑐𝑙𝑎

𝜆 = 𝑐𝑎𝑙𝑜𝑟 𝑙𝑎𝑡𝑒𝑛𝑡𝑒 𝑑𝑒 𝑣𝑎𝑝𝑜𝑟𝑖𝑧𝑎𝑐𝑖𝑜𝑛 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑒

𝑤 = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝑚𝑎𝑠𝑖𝑐𝑜 𝑑𝑒 𝑣𝑎𝑝𝑜𝑟

Agrupando la ecuación (1)

𝑑𝑇𝑓𝑠
(𝑚𝑓 𝑐𝑝𝑓 + 𝑚𝑖𝑐 𝑐𝑝𝑖𝑐 ) = 𝐹𝑐𝑝𝑓 (𝑇𝑓𝑒 − 𝑇𝑓𝑠 ) + 𝜆𝑤 (2)
𝑑𝑡
Dividiendo la ecuación (2) entre 𝐹𝐶𝑝𝑓 𝜆

(𝑚𝑓 𝑐𝑝𝑓 + 𝑚𝑖𝑐 𝑐𝑝𝑖𝑐 ) 𝑑𝑇𝑓𝑠 1 1

= (𝑇𝑓𝑒 − 𝑇𝑓𝑠 ) + 𝑤
𝐹𝑐𝑝𝑓 𝜆 𝑑𝑡 𝜆 𝐹𝑐𝑝𝑓

Hacemos cambios de variables para facilitar los cálculos

(𝑚𝑓 𝑐𝑝𝑓 + 𝑚𝑖𝑐 𝑐𝑝𝑖𝑐 )

𝑘1 =
𝐹𝑐𝑝𝑓 𝜆

1
𝑘2 =
𝜆
 1
𝑘3 =
𝐹𝑐𝑝𝑓

𝑇 = 𝑇𝑓𝑠

𝑇𝑖 = 𝑇𝑓𝑒

Obtenemos la ecuación de balance de energía:

𝑑𝑇
𝑘1 = 𝑘2 (𝑇𝑖 − 𝑇) + 𝑘3 𝑤 (3)
𝑑𝑡
𝑑𝑇
Despejando 𝑑𝑡

𝑑𝑇 𝑘2 𝑘3
= (𝑇𝑖 − 𝑇) + 𝑤
𝑑𝑡 𝑘1 𝑘1