MATEMÁTICAS ESPECIALES
Examen 3-2021-II
Departamento de Ingeniería Química
UNIVERSIDAD DE ANTIOQUIA
Docente: Luis Alberto Díaz M.
Integrantes: Angye Melissa Osorio Riveros, CC: 1006775397. Estefanía Restrepo Zapata, CC:
1038771888. Daniela Velásquez Cuartas, CC: 1152467646.
Este informe se realizó usando los parámetros ∝ = 5, 𝛽𝛽 = 5, 𝑦𝑦 = 7.
1. Se pide hallar la cantidad de oxígeno disuelto en agua desde una temperatura de 5°C hasta 30°C
aumentando en 5°C para una cantidad determinada de cloruro. Para ello, se aplicó el método
combinado de Simpson 3/8 y Simpson 1/3 (Tabla 1), tomando como límite superior 30°C, límite
inferior 5°C y un espaciamiento h de 5. El área total obtenida fue de 198.6042, que corresponde
a los miligramos de oxígeno disueltos por litro de solución.
T F(T) (O2 disuleto) Áreas Método
5 10.5
10 9.2 131.4375 Simpson 3/8
15 8.2
20 7.4
25 6.7 67.1667 Simpson 1/3
30 6.1
198.6042
Área total
Tabla 1: Aplicación del método de Simpson combinado para hallar cantidad de oxígeno en solución.
2. Para este problema, se desea hallar las dimensiones de un reactor de flujo pistón cilíndrico con
una relación del diámetro con respecto a la longitud de 3:1, que permita una conversión del
80% y un flujo molar de entrada determinado (FAO). Para ello se dan datos de velocidad de
reacción (-rA) como función de la conversión (X) a 300°F y 10 atm. La función para el volumen en
esta clase de reactor está dada por la ecuación 2.1.
(2.1)
Se aplicó el método del trapecio y el método de Simpson 1/3 (Tabla 2.1), para hallar la función
de X como el área debajo de una curva y resolver para el volumen.
Área (m2)
X -rA F(X) Trapecio Simpson 1/3
0 5.3 0.18868 0.01905
0.03860
0.1 5.2 0.19231 0.01962
0.2 5 0.20000 0.02111
0.04463
0.3 4.5 0.22222 0.02361
0.4 4 0.25000 0.02765
0.06207
0.5 3.3 0.30303 0.03515
0.6 2.5 0.40000 0.04778
0.11407
0.7 1.8 0.55556 0.06778
0.8 1.25 0.80000 0.26175 0.25937
Área total
Tabla 2.1: Aplicación del método de trapecio y Simpson 1/3 para hallar una integral como área bajo una curva.
Finalmente, se determinaron las dimensiones propuestas para el reactor a partir del volumen
total (Tabla 2.2).
3
Método Longitud [m] Diámetro [m] Radio [m] Volumen [m ]
Trapecio 1.3752 0.4584 0.2292 0.2269
Simpson 1/3 1.3710 0.4570 0.2285 0.2249
Tabla 2.2: Dimensiones para el reactor PFR cilíndrico.
3. Se debe hallar numéricamente la fuerza total ejercida sobre un mástil sobre un bote, descrita
por la ecuación 3.1.
(3.1)
Se usó el método de Simpson 1/3 para 32 subintervalos, con un espaciamiento h de 0.9375 para
hallar la función de z como el área bajo una curva. El área total hallada fue de 32075.67101, que
corresponde a la fuerza total ejercida sobre el mástil del bote.
4. En la primera parte de este problema, se pide hallar de forma numérica la distribución de la
densidad de masa (f(x)), utilizando diferencias finitas con orden de error 2 y para diámetros de
partícula de polvo (x) desde 0 hasta 150. Para ello se hizo uso de la ecuación de Rosin-Rammler-
Bennet (4.1), siendo su derivada con respecto a x, la distribución de la densidad de masa (Figura
4.1)
(4.1)
Figura 4.1: Distribución de densidad de masa con respecto al diámetro de partícula de polvo.
Posteriormente para hallar el área superficial por masa de polvo (Sm), descrita por la ecuación
4.2, se utilizó el método de Simpson 3/8 para 129 subintervalos.
(4.2)
Se encontró que el área superficial por masa de polvo para un diámetro mínimo (dmín) de 1 μm
hasta un diámetro de 130 μm y densidad de partícula (ρ) de 1 g/cm3, es de 0.4386 cm2/g.
5. Para dar solución a la cadena de desintegración radiactiva, en primer lugar, se construye la
matriz a desarrollar.
Figura 5.1: Matriz de ecuaciones para desintegración radiactiva.
La matriz puede interpretarse como una matriz que contiene diez ecuaciones diferenciales; la
solución al problema se hace mediante el método de Runge Kutta 2. La plantilla donde se ingresan
las ecuaciones diferenciales es de máximo 10 ecuaciones.
Figura 5.2. Solución del sistema matricial.
El número de subintervalos fue de 70 y el problema de valor inicial está condicionado por un
intervalo de tiempo de [0;5].
El comportamiento de desintegración radiactiva se evidencia en la siguiente figura.
Figura 5.3: Cadena de desintegración radiactiva en el tiempo.
6.
7. Se debe encontrar la solución numérica del problema dado por las ecuaciones (7.1), (7.2) y (7.3)
y sujeto a las condiciones (7.4), (7.5) y (7.6) desde t=0 hasta t=π/2 la solución analítica está sujeta
a las ecuaciones (7.7), (7.8) y (7.9) respectivamente.
(7.1)
(7.2)
(7.3)
(7.4)
(7.5)
(7.6)
(7.7)
(7.8)
(7.9)
Este ejercicio se desarrolló mediante el método de disparo con RK4 (Tradicional) para 24
subintervalos, e implementado una estrategia iterativa (Solver) para cambiar sistemáticamente
la condición inicial creada, hasta cumplir la condición final que nos daba el problema original,
como se muestra en la tabla 7.1.
Condición Valor
y1(0) [-] 4
y2(pi/2) [-] 163,835,924 →Inventada
y3(pi/2) [-] 328,003,173 →Inventada
Valor Solver
y2(pi/2) [-] 396,733,417
y3(pi/2) [-] 228,628,643
Tabla 7.1: Condiciones de frontera halladas con solver.
Los resultados para y encontrados de forma numérica por medio del método y los encontrados
de forma analítica se evidencian en las siguientes figuras.
Figura 7.1: Valores de Y1 numéricos y analíticos en el tiempo.
Figura 7.2: Valores de Y2 numéricos y analíticos en el tiempo.
Figura 7.3: Valores de Y3 numéricos y analíticos en el tiempo
8. En este ejercicio se propone una ecuación diferencial (8.1) la cual sirve para describir el
movimiento de un satélite con respecto a la tierra, bajo ciertas condiciones de frontera
determinadas. Para ello se utilizó el método de diferencias finitas con 20 subintervalos con los
parámetros y condiciones de frontera mostrados en la tabla 8.2.
(8.1)
Parámetros Valor
to 0
tf 3.1416
Ꜫ 0.1
#Sub.intervalos 20
h= Δx[-] 0.1571
Condiciones de fontera
x(0) [-] 0.9
y( π ) [-] 0
Tabla 8.1: Parámetros y condiciones de frontera para el movimiento satelital.
También se da la información necesaria para hallar la solución analítica del problema, los datos
obtenidos por medio de las dos formas de solución se observan en la figura 8.1
Figura 8.1: Movimiento satelital con respecto a la tierra en un tiempo determinado.
9. El modelo matemático de un catalizador cilíndrico que se encuentra en una reacción química
de orden n bajo condiciones isotérmicas puede escribirse por medio de las siguientes
ecuaciones.
(9.1)
(9.2)
(9.3)
Donde kCn es la velocidad de reacción por unidad de masa del catalizador, y De es la difusividad
efectiva, C0 es la concentración constante en el seno del fluido, y R es el radio de la partícula. Al
implementar el cambio de variable propuesto en las ecuaciones 9.4 y 9.5, la ecuación diferencial
se queda expresada como se muestra en 9.6.
(9.4)
(9.5)
(9.6)
para n=0. Encuentre el perfil de y vs. x Para Solucionar este ejercicio s e utilizan los siguientes
parámetros
Parámetro Valor
2
De(m /hora) 3.60E-06
3
P PK(mol/m *hora) 1
C0(mol/m3) 0.2
n 0
R 0.001
ϕ 0.3472
Este ejercicio se desarrolló mediante el método de disparo y RK2 (Raizton y Rabinowitz) con 20
subintervalos y haciendo la reducion de orden respectiva y el cambio de variable, se llega a la
solucion evidenciada en la Figura 9.1.
Figura 9.1: Valores de Y con respecto a X
10. Se pide encontrar la solución numérica del problema dado por las ecuaciones (10.1) a (10.3), y
sujeto a las condiciones (10.4), desde t = 0 min hasta t = 100 min.
(10.1)
(10.2)
(10.3)
(10.4)
Los valores de los parámetros requeridos para su solución se encuentran en la tabla (10.1). Las
unidades son consistentes con las ecuaciones.
Parámetros del problema
t0 0 W [kg/min] 100 T steam [°C] 250
T 1 (0) 20 Cp [kJ/kg°C] 2 M [Kg] 1000
T 2 (0) 20 To [°C] 20
T 3 (0) 20 UA [kJ/min°C] 10
tf 100
Tabla 10.1: Parámetros para problema de intercambio de calor en una serie de tanques
Para solucionar este ejercicio se ha pedido implementar el método de Runge-Kutta 3 con 105
subintervalos, seguido a esto se establece el proceso iterativo para el problema, partiendo de la
definición del método RK3 y con las condiciones iniciales dadas. El comportamiento de la
temperatura en el tiempo se evidencia en la figura 10.1
Tabla 10.1: Comportamiento de la temperatura en el tiempo.
11. En este ejercicio se debe encontrar la solución numérica del problema de la ecuación (11.1),
sujeta a las condiciones (11.2) y (11.3), cuya solución es (11.4), desde z = 0 hasta z = 0.001
(11.1)
(11.2)
(11.3)
(11.4)
Adicional a esto sus respectivos parámetros son los siguientes:
Parametros del problema
3
CAO [Kmol/m ] 0.2
-1
k [s ] 0.001
2
DAB [m /s] 1.2E-09
L [m] 0.001
Para la solución de este ejercicio se implementó el método de diferencias finitas para el Nodo 1
hacia adelante, del Nodo 2-14 diferencias finitas centradas, y el Nodo 15 diferencias finitas hacia
atrás mediante las siguientes ecuaciones.
𝑘 2
𝐶𝐴2 − 2 + ℎ 𝐶𝐴1 = −𝐶𝐴0 NODO 1
𝐷𝐴𝐵
𝑘 2
𝐶𝐴𝑖 +1 − 2 + ℎ 𝐶𝐴𝑖 + 𝐶𝐴𝑖 −1 = 0 NODO 2-14
𝐷𝐴𝐵
3𝐶𝐴15 − 4𝐶𝐴14 + 𝐶𝐴13 = 0 NODO 15
Mediante las cuales se logró obtener las diferentes soluciones numéricas y analítica para cada
sistema de ecuaciones que logramos evidenciar en la siguiente figura.
Figura 11.1: CAO con respecto a Z
12. Para un proceso de destilación binaria por lotes de dos compuestos designados como 1 y 2, se
puede expresar las moles de líquido remanente L como función de la fracción molar del
compuesto 2, x2, mediante las ecuaciones 12.1, 12.2 y mediante el modelo matemático 12.3.
(12.1)
(12.2)
(12.3)
Este ejercicio se desarrolló inicialmente mediante el método de Runge-Kutta 2 (punto medio)
luego se implementó la función Solver para obtener las moles necesarias hasta la pauta final
que nos daba el enunciado, dando como resultado, para una fracción de compuesto 2 de 0.8, la
cantidad de moles son 14.0924.
13. Encontrar la solución numérica del problema dado por las ecuaciones (13.1) y (13.2), y sujeto a
las condiciones (13.3) y (13.4), desde t = 0 hasta t = 20.
(13.1)
(13.2)
(13.2)
(13.3)
(13.4)
(13.5)
(13.6)
Para solucionar este ejercicio se realiza, partiendo de la definición del método de Runge-Kutta
2 (Heun) con 70 subintervalos y la solución se evidencia en la siguiente figura.
Figura 13.1: Comportamiento de Sustrato y Biomasa en el tiempo.
Luego mediante la función Solver, se encontró un t = 4,692864472 En este punto se puede
asumir que el sustrato se consumió complétame por lo tanto en la figura 13.1 solo se presentan
algunos resultados puesto que cuando se consume el sustrato ya no existiría un fenómeno que
se pueda representar físicamente.
14. Se tiene un sistema formado por una masa m (kg) que está unida a un resorte (k =2 N/m) que
parte del reposo desde 0.5 (m), con una constante de amortiguamiento q = 1,2 e impulsado por
una fuerza periódica dada p = (5cos(4t)). Determinar la posición y la velocidad de la masa en
t=1, t=3, t=5 y t=50.
Para hallar la solución de este ejercicio se utiliza el método de Runge-Kutta 3. Se realiza un
cambio de variable en la función inicial, esto con el fin de realizar una reducción de orden a la
ecuación diferencial. Una vez realizado se aborda el problema, el cual está condicionado a los
siguientes parámetros.
Parámetros Valor
t 0 [s] 0
t f [s] 50
u (0) 0.5
u[1] 0
# subintervalos 50
K (N/m) 2
q 1.2
p 5cos(4t)
Figura 14.1: Posición y velocidad en el tiempo.
15. Se analiza un tanque que contiene 100 galones de agua fresca, al cual se habilita un flujo
constante de entrada de dos galones por litro (Q), inicialmente de agua salada con una
concentración (C1) de 0.5 lb/gal y después de 10 minutos, se cambia el flujo por uno de agua
fresca (C1=0). El tanque también tiene un flujo de salida (Q2) de tal forma que el volumen se
mantenga constante. Para hallar la cantidad numérica de sal que habrá en el tanque luego de
20 minutos, se trabaja el problema en dos partes, haciendo uso del método de Runge-Kutta 3,
con la ecuación diferencial 15.1, y estableciendo los parámetros para cuando el flujo es de agua
salada (desde t0, hasta t1) y cuando el flujo es de agua fresca (desde t1, hasta t2), como se muestra
en la tabla 15.1. La concentración final en la primera parte del problema (X1), será la
concentración inicial en la segunda parte del problema.
(15.1)
Parámetros Valor
# subintervalos 80
h=Δt 0.25
Q 1 [gal/min] 2
Q 2 [gal/min] 2
C 1 para t0-t1 [lb/gal] 0.5
C 1 para t1-t2 [lb/gal] 0
v 0 [gal] 100
t 0 [min] 0
t 1 [min] 10
t 2 [min] 20
X0 0
X1 9.06346239
Tabla 315.1: Parámetros para comportamiento del agua en el tanque.
La concentración obtenida en la segunda parte del problema será la concentración total (X2total)
trascurridos los 20 minutos, el valor de dicha concentración fue de 7.4205 lb/gal.
El comportamiento de la concentración de sal en el tanque se evidencia en la siguiente figura.
Figura 15.1: Concentración de sal en el tanque transcurridos 20 minutos.
16. Se desea determinar el valor de x dentro de un intervalo dado, de modo que y = 50.
Para hallar la solución a la siguiente ecuación se utiliza el método de Runge Kutta 3, mediante
la ecuación 16.1 con un numero de subintervalos igual a veinte y condicionado inicialmente en
x a un intervalo desde x = -0.5 hasta x =-0.1
(16.1)
Para hallar el valor de x para el cual y = 50 se realiza por medio de interpolación lineal,
obteniendo para x un valor de -0.1563
Figura 16.1: Solución gráfica para los valores de y en función de x.
17. a. Se busca para la ecuación diferencial, hallar y (2).
Para hallar la solución de este ejercicio se utiliza el método de Runge Kutta 3 y utilizando la
ecuación mostrada a continuación.
(17.1)
La ecuación diferencial está condicionada por los siguientes parámetros.
Tabla 17.1: Condiciones para la ecuación diferencial.
Figura 17.1: Solución gráfica para los valores de Y(2).
b. Se busca para la ecuación diferencial, hallar y (0,2) - y’ (0,1)
Para hallar la solución de este ejercicio se utiliza el método de Runge Kutta 3, haciendo uso de
la ecuación 17.2. Se realiza un cambio de variable en la función inicial, con el fin de realizar una
reducción de orden a la ecuación diferencial.
(17.2)
El problema se encuentra condicionado a los siguientes parámetros:
Tabla 17.2: Condiciones iniciales para la ecuación diferencial.
Figura 17.3: Solución gráfica para los valores de Y (0,2) y Y’ (0,1) 17b.
c. Graficar y =f(t) y hallar y(2 π)
Para hallar la solución de este ejercicio se utiliza el método de Runge Kutta 3 y las siguientes
ecuaciones.
(17.3)
(17.4)
La ecuación diferencial está condicionada por los siguientes parámetros.
Tabla 17.3 Condiciones iniciales para la ecuación diferencial.
Figura 17.4: Solución grafica para la ecuación diferencial.
Referencias.
Chapra, S. C. C., & Raymond, P. (2007). Métodos numéricos para ingenieros.
México. DF McGraw-Hill.
Barrera, R & Diaz, L. (2020). Métodos numéricos en Excel y Matlab con
aplicaciones en ingeniería. Universidad de Antioquia.
Fogler, H. S (2010). Essentials of Chemical Reaction Engeneering. Pearson
Education.
Green, D. W. (1997). Perry´s Chemical Engineerps Handbook.
Reklaitis, G, V., & Schneider, D, R. (1983). Introduction to material and energy
balances (pp. 55-56). New York: Wiley.
