calculo de perfiles de concentración con MATLAB

Anexo III
Cálculo de perfiles de concentraciones mediante Matlab
En un reactor PFR no hay ningún tipo de mezcla en la dirección del flujo, por lo
que el fluido pasa sin mezclarse como si semoviera a lo largo del reactor
impulsado por un pistón. Además en la dirección radial hay mezcla completa,
por lo tanto todas las propiedades del fluido en la sección non uniformes. Por lo
tanto en cadadV el reactor PFR, este se comporta como un CSTR Es por eso
que podemos hacer la consideración de que un reactor PFR se comporta como
muchos reactores CSTR en serie, como se puede ver en el Gráfico10, que
relaciona el volumen total del sistema con las concentraciones de un reactivo
A, donde a medida que aumenta el número de CSTR’s en serie, su
comportamiento se va asemejando al de un reactorPFR.
-----------Como simplificación se halló el perfil de concentraciones como si se utilizara un
reactor PFR ideal, ya que se disponía de mayor conocimiento de Matlab para la
resolución del tipode sistema de ecuaciones que se obtenía al realizar los
cálculos pertinentes. Como se aprecia en el Grafico 10, para una determinada
concentración de entrada, la concentración de salida coincidepara ambos dos
modelos de reacción, por lo que la simplificación se considera válida.

Si planteamos el balance de materia en un PFR para los triglicéridos
obtenemos lo siguiente:
𝐹 𝑇𝐺 − (𝐹 𝑇𝐺 + 𝑑𝐹𝑇𝐺 ) = (−𝑟 𝑇𝐺 ) ∗ 𝑑𝑉
−𝑑𝐹 𝑇𝐺 = (−𝑟 𝑇𝐺 ) ∗ 𝑑𝑉
−𝑑𝐶 𝑇𝐺 ∗ 𝑞 = (−𝑟 𝑇𝐺 ) ∗ 𝑑𝑉
−𝑑𝐶 𝑇𝐺 ∗ 𝑞
= (−𝑟 𝑇𝐺 )
𝑑𝑉
−

𝑑𝐶 𝑇𝐺
= (−𝑟 𝑇𝐺 )
𝑑𝜏

𝒅𝑪 𝑻𝑮
= −(−𝒓 𝑻𝑮 )
𝒅𝝉

Para los diglicéridos:
𝐹 𝐷𝐺 − (𝐹 𝐷𝐺 + 𝑑𝐹 𝐷𝐺) = (−𝑟 𝐷𝐺 ) ∗ 𝑑𝑉
−𝑑𝐹 𝐷𝐺 = (−𝑟 𝐷𝐺 ) ∗ 𝑑𝑉
−𝑑𝐶 𝐷𝐺 ∗ 𝑞 = (−𝑟 𝐷𝐺 ) ∗ 𝑑𝑉
−𝑑𝐶 𝐷𝐺 ∗ 𝑞
= (−𝑟 𝐷𝐺 )
𝑑𝑉
𝑑𝐶 𝐷𝐺
= (−𝑟 𝐷𝐺 )
𝑑𝜏

−

𝒅𝑪 𝑫𝑮
= −(−𝒓 𝑫𝑮 )
𝒅𝝉
Para los monoglicéridos:
𝐹 𝑀𝐺 − (𝐹 𝑀𝐺 + 𝑑𝐹 𝑀𝐺) = (−𝑟 𝑀𝐺 ) ∗ 𝑑𝑉
−𝑑𝐹 𝑀𝐺 = (−𝑟 𝑀𝐺 ) ∗ 𝑑𝑉
−𝑑𝐶 𝑀𝐺 ∗ 𝑞 = (−𝑟 𝑀𝐺 ) ∗ 𝑑𝑉
−𝑑𝐶 𝑀𝐺 ∗ 𝑞
= (−𝑟 𝑀𝐺 )
𝑑𝑉
−

𝑑𝐶 𝑀𝐺
= (−𝑟 𝑀𝐺 )
𝑑𝜏

−

𝒅𝑪 𝑴𝑮
= (−𝒓
𝒅𝝉

𝑴𝑮 )

Para la glicerina:
𝐹 𝐺 − (𝐹 𝐺 +...