UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE INGENIERÍA QUÍMICA
Métodos Numéricos Para Ingeniería Química (PI524B) – 2022-1
Fecha: 13-08-2022
EXAMEN SUSTITUTORIO
Tiempo: máximo 120 minutos
▪ El examen final consta de 4 preguntas. El puntaje total es 20.
▪ Las respuestas deben ser justificadas de manera detallada.
1. Ajuste de modelos a datos
Cuando una sustancia química ingresa a un cuerpo de agua puede que reaccione química o
bioquímicamente transformándose en otra sustancia quimica. En la mayoría de las reacciones
importantes en el agua el orden de la reacción es un entero. Las reacciones que con mayor
frecuencia ocurren en el agua son de la forma:
𝑑𝐶 𝑛
𝑑𝑡
=− 𝑘𝐶
donde 𝑛 es el orden de la reacción y 𝐶 es la concentración del reactante (sustancia química).
Para estudiar la fotodegradación de una muestra de agua con respecto a la salinidad, se disuelve una
pequeña cantidad de bromo en el agua, se coloca en un frasco de vidrio claro y limpio y se expone a
la luz solar. Se presenta los siguientes datos del estudio:
𝑡 (días) 10 20 30 40 50 60
𝐶 (mg/L) 3,52 2,48 1,75 1,23 0,87 0,61
a) A partir del método integral se obtienen los siguientes modelos para la reacción:
Modelo 1: 𝐶 = 𝐶0 − 𝑘𝑡
−𝑘𝑡
Modelo 2: 𝐶 = 𝐶0𝑒
1 1
Modelo 3: 𝐶
= 𝐶0
+ 𝑘𝑡
donde 𝐶0 es la concentración inicial y 𝑘 es la constante de la velocidad de reacción. Usando el
método de mínimos cuadrados determine 𝐶0 en mg/L y 𝑘 en las unidades correspondientes a
cada modelo. (3
puntos)
2
b) Calcule el coeficiente de determinación (𝑟 ) para cada modelo mediante:
2 𝑆𝑡−𝑆𝑟
𝑟 = 𝑆𝑡
𝑚 2 𝑚 2
(
donde 𝑆𝑡 = ∑ 𝐶 − 𝐶𝑖
𝑖=1
) ( *
)
y 𝑆𝑟 = ∑ 𝐶𝑖 − 𝐶𝑖 , 𝐶𝑖 es la concentración medida, 𝐶 es la media
𝑖=1
*
aritmética de los 𝑚 datos de concentración medida y 𝐶𝑖 es la concentración calculada con el
modelo. Indique que modelo se ajusta a los datos del estudio y justifique. (2
puntos)
2. Ecuación diferencial ordinaria
Las sustancias químicas 𝐴 y 𝐵 reaccionan para obtener una nueva sustancia química 𝐶, de acuedo a
la reacción química:
𝐴 + 𝐵→𝐶
La ecuación cinética de la reacción química es:
𝑑𝐶𝐴
𝑑𝑡
= 𝑘𝐶𝐴𝐶𝐵
Inicialmente las concentraciones de 𝐴, 𝐵 y 𝐶 son 100, 150 y 0 mmol/L, respectivamente. Si luego de
haber transcurrido 1 h se formó 15 mmol/L de 𝐶
a) Determine la constante de la velocidad de reacción 𝑘 en L/(mmol.h) (2
puntos)
b) Usando el método de Euler, determine la concentración de 𝐶 en mmol/L después de haber
transcurrido 2 h. Considere como tamaño de paso ∆𝑡 = 0, 1 h. (3 puntos)
3. Ecuaciones No Lineales
Para el flujo de una esfera a través de un fluido estancado Lapple y Shepherd (1940) presentaron la
relación entre el factor de fricción ( 𝑓 ) y el número de Reynolds ( 𝑅𝑒𝑝 )
24 18.5 5
𝑓 = { 𝑅𝑒 ; 𝑅𝑒𝑝 < 2 0.6 ; 2≤𝑅𝑒𝑝 < 500 0. 44 ; 500≤𝑅𝑒𝑝 < 2 . 10
𝑝 𝑅𝑒𝑝
En 1896 Turton y Levenspiel proponen la siguiente ecuación:
𝑓 =
24
𝑅𝑒𝑝 ( 1 + 0. 173 𝑅𝑒 ) +
0.657
𝑝
0.413
1+16300 𝑅𝑒𝑝
−1.09
Utilizando las ecuaciones mencionadas y empleando un método abierto, calcular la velocidad terminal ( 𝑉𝑡 ) de
una burbuja de agua de 5 mm de diámetro ( 𝐷𝑝 ) cayendo a través de aire estancado a 20°C, si se tiene los
siguientes datos adicionales (5
puntos)
𝐷𝑝 𝑉𝑡 ρ
𝑅𝑒𝑝 = µ
Propiedades a 20°C
𝐾𝑔
Para el agua: ρ𝑝 = 999 3
𝑚
𝐾𝑔 −6 𝐾𝑔
Para el aire: ρ = 1. 2047 3 ; μ = 18. 17 * 10 𝑚. 𝑠
𝑚
4. Integrales dobles
Para la siguiente integral doble:
𝑊 = ∬𝑅2 𝑥 − ( 2 2
3
𝑦+
𝑎
3 )𝑑𝐴
𝑅= { (𝑥 , 𝑦 ) ∈ 𝑅2 / 0≤𝑥≤6 ; − 2≤𝑦≤1 }
a) Determinar el valor de 𝑎 si se sabe que al utilizar integración numérica de Simpson 1/3 y
3/8 donde corresponda, con:
𝑥𝑖 = 0, 1, 2 , 4, 6 𝑦𝑖 = − 2 , − 1 . 0, 1
b) Con el valor de 𝑎 encontrado evalué la integral doble, empleando el método de
Gauss-Legendre con dos puntos
Comente sus resultados (5
puntos)
Los profesores