Métodos Numéricos para la Ingeniería Química- Grupo 5
UNIVERSIDAD NACIONAL DE TRUJILLO
"Año del Fortalecimiento de la Soberanía Nacional"
FACULTAD DE INGENIERÍA QUÍMICA
ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA QUÍMICA
TEMA:
Práctica de Clase
CURSO:
Métodos Numéricos para la Ingeniería Química
DOCENTE:
Guillermo David Evangelista Benítez
GRUPO: 5
AUTORES:
Díaz Linares, Ítalo (coordinador)
Cueva Gutiérrez, Astrid Elizabeth
Vásquez Delgado, Consuelo Sheyla
Colonia Chinchey, Yoselyn Judith
Horna Altamirano, Paola Lucía
Esquivel Hermenegildo, Melva Nerida
Sedano Ibáñez, Renzo Martin
Vásquez Cruz, Piero Cesar
Roncal Marocho, Gianella
Loyaga Aguirre Joel Alexis
Métodos Numéricos para la Ingeniería Química- Grupo 5
MÉTODOS NUMÉRICOS PARA LA INGENIERÍA QUÍMICA
EJERCICIO
SOLUCIÓN DEL PROBLEMA
Entrada -Salida - Lo que reacciona =Acumulación
𝑛 𝑑𝐶𝐴1
𝐹𝐶𝐴0 + 𝐹𝑅𝐶𝐴3 −(𝐹 + 𝐹𝑅 )𝐶𝐴1 −𝑘1 𝑉1 𝐶𝐴1 = 𝑑𝑡
𝑛 𝑑𝐶𝐴2
(𝐹 + 𝐹𝑅 )𝐶𝐴1 −(𝐹 + 𝐹𝑅 )𝐶𝐴2 −𝑘1 𝑉2 𝐶𝐴2 =
𝑑𝑡
𝑛 𝑑𝐶𝐴3
(𝐹 + 𝐹𝑅 )𝐶𝐴2 −(𝐹 + 𝐹𝑅 )𝐶𝐴3 −𝑘1 𝑉3 𝐶𝐴3 = 𝑑𝑡
Datos:
𝑭 = 10 𝑭𝑹𝟏 = 1 𝑭𝑹𝟐 = 2 𝑭𝑹𝟑 = 2 𝑪𝑨𝟎 = 1
𝑽𝟏 = 100 𝑽𝟐 = 50 𝑽𝟑 = 50 𝒌𝟏 = 0.1
Para 𝑪𝑨𝟏 :
𝑑𝐶𝐴1
𝐹. 𝐶𝐴0 + 𝐹𝑅1 . 𝐶𝐴3 − (𝐹 + 𝐹𝑅1 )𝐶𝐴1 − 𝑘1 𝑉1 𝐶𝐴1 =
𝑑𝑡
𝑑𝐶𝐴1
10(1) + 1. 𝐶𝐴3 − 11. 𝐶𝐴1 − 10. 𝐶𝐴1 =
𝑑𝑡
𝑑𝐶𝐴1
−21𝐶𝐴1 + 𝐶𝐴3 + 10 =
𝑑𝑡
Métodos Numéricos para la Ingeniería Química- Grupo 5
𝑑𝐶𝐴1
∫(−21 𝐶𝐴1 + 𝐶𝐴3 + 10) 𝑑𝑡 = ∫ 𝑑𝑡
𝑑𝑡
−21𝐶𝐴1 𝑡 + 𝐶𝐴3 𝑡 + 10𝑡 = 𝐶𝐴1
Si t = 1
−21𝐶𝐴1 + 𝐶𝐴3 + 10 = 𝐶𝐴1
−22𝐶𝐴1 + 𝐶𝐴3 = -10 ………………….. (1)
Para 𝑪𝑨𝟐 :
𝑑𝐶𝐴2
(𝐹 + 𝐹𝑅2 ). 𝐶𝐴1 − (𝐹 + 𝐹𝑅2 ). 𝐶𝐴2 − 𝑘1 𝑉2 𝐶𝐴2 = 𝑑𝑡
𝑑𝐶𝐴2
12𝐶𝐴1 − 12𝐶𝐴2 − 5𝐶𝐴2 =
𝑑𝑡
𝑑𝐶𝐴2
12𝐶𝐴1 − 17𝐶𝐴2 =
𝑑𝑡
𝑑𝐶𝐴2
∫(12 𝐶𝐴1 − 17𝐶𝐴2 ) 𝑑𝑡 = ∫ 𝑑𝑡
𝑑𝑡
12𝐶𝐴1 𝑡 − 17𝐶𝐴2 𝑡 = 𝐶𝐴2
Si t = 1
12𝐶𝐴1 − 17𝐶𝐴2 = 𝐶𝐴2
12𝐶𝐴2 − 18𝐶𝐴2 = 0 ……………………. (2)
Para 𝑪𝑨𝟑 :
𝑑𝐶𝐴3
(𝐹 + 𝐹𝑅3 )𝐶𝐴2 − (𝐹 + 𝐹𝑅3 )𝐶𝐴3 − 𝑘1 𝑉3 𝐶𝐴3 =
𝑑𝑡
𝑑𝐶𝐴3
12𝐶𝐴2 − 12𝐶𝐴3 − 5𝐶𝐴3 =
𝑑𝑡
𝑑𝐶𝐴3
12𝐶𝐴2 − 17𝐶𝐴3 =
𝑑𝑡
Métodos Numéricos para la Ingeniería Química- Grupo 5
𝑑𝐶𝐴3
∫(12 𝐶𝐴2 − 17𝐶𝐴3 ) 𝑑𝑡 = ∫ 𝑑𝑡
𝑑𝑡
12𝐶𝐴2 𝑡 − 17𝐶𝐴3 𝑡 = 𝐶𝐴3
Si t = 1
12𝐶𝐴2 − 17𝐶𝐴3 = 𝐶𝐴3
12𝐶𝐴2 − 18𝐶𝐴3 = 0 ………………………. (3)
Juntando las ecuaciones (1), (2) y (3) se obtiene la siguiente matriz:
−22𝐶𝐴1 + 0𝐶𝐴2 + 𝐶𝐴3 = −10
{ 12𝐶𝐴1 − 18𝐶𝐴2 + 0𝐶𝐴3 = 0
0𝐶𝐴1 + 12𝐶𝐴2 − 18𝐶𝐴3 = 0
−𝟐𝟐 𝟎 𝟏 −𝟏𝟎
𝟏𝟐 −𝟏𝟖 𝟎 𝟎
𝟎 𝟏𝟐 −𝟏𝟖 𝟎
Métodos Numéricos para la Ingeniería Química- Grupo 5
MATLAB
Código:
clear all, clc
disp('Curso: Métodos Numéricos para Ingeniería Química')
disp('Prof.: Dr. Guillermo Evangelista Benites')
disp('Clase: Solución de ecuaciones lineales')
disp('Tema : Método de Gauss')
disp('Fecha: Jueves, 9 de Junio del 2022')
disp(' ')
disp(' TAREA N°3')
disp('Ejercicio 3.32 - Pág.281 - Antonio Nieves y Federico Domínguez')
disp('----------------')
disp('Encontrar las concentraciones de A a régimen permanente')
disp('en cada reactor mostrado en la figura 3.24')
disp('----------------')
disp(' ')
disp('Solución: ')
disp('-------- ')
disp(' ')
disp('Sistema de ecuaciones de los reactores')
disp('-22*CA1 + 0*CA2 + 1*CA3 = -10')
disp(' 12*CA1 - 18*CA2 + 0*CA3 = 0')
disp(' 0*CA1 + 12*CA2 -18*CA3 = 0')
disp(' ')
disp('Matriz de coeficientes: ')
A=[-22 0 1; 12 -18 0; 0 12 -18]
disp('Vector de términos independientes: ')
B=[-10; 0; 0]
[n n] = size(A);
A = [A';B']';
x = zeros(n,1);
disp(' ')
disp('Resolviendo el sistema por Eliminación de Gauss:')
disp(' ')
for p = 1:n
for k = [1:p-1,p+1:n];
m = -A(k,p)/A(p,p);
A(k,:) = A(k,:)+m*A(p,:)
end
end
disp(' ')
disp('Vector solución: ')
C = A(:,n+1)./diag(A)
l=length(C);
for i=1:l
fprintf('Concentración de A en el reactor %d: %7.4f mol/L\n',i,C(i)')
end
Métodos Numéricos para la Ingeniería Química- Grupo 5
RESULTADO:
Métodos Numéricos para la Ingeniería Química- Grupo 5
Métodos Numéricos para la Ingeniería Química- Grupo 5
POLYMATH
Métodos Numéricos para la Ingeniería Química- Grupo 5
