2º Parcial de Métodos Numéricos
Programa de Ingeniería Química
Prof. Juan Carlos Burgos, Ph.D.
Fecha de entrega: Antes de febrero 26 de 2019 a las 11:59 pm
Taller para realizar de manera individual. Entregar un reporte escrito donde se describa paso a
paso cualquier manipulación o derivación de ecuaciones antes de aplicar los métodos numéricos.
Adjuntar los archivos con los códigos usados.
1. (1.0) Cierto material sigue un comportamiento exponencial en tensión uniaxial como función de
la deformación, lo cual se representa mediante la siguiente expresión.
𝐸$ ()
𝜎= (𝑒 − 1)
𝑎
Donde 𝜎 = esfuerzo, 𝜀 = deformación, y 𝐸$ y 𝑎 son constantes del material que se determinan
experimentalmente. La tabla siguiente contiene los datos esfuerzo-deformación del material.
Datos 1 2 3 4 5 6 7 8 9 10 11 12 13 14
𝜎 (kN/m2) 0 87.8 96.6 176 263 351 571 834 1229 1624 2107 2678 3380 4258
𝜀 - 10-3 (m/m) 0 153 204 255 306 357 408 459 510 561 612 663 714 765
Evalúe las dos constantes (𝐸$ y 𝑎) del material usando valores de d𝜎/d𝜀 contra 𝜎 obtenidos a
partir de los datos 5 al 13 en la tabla. Aplique los métodos numéricos que considere
necesarios.
2. Un investigador informa sobre los datos tabulados a continuación de un experimento para
determinar la tasa de crecimiento de bacterias k (por día), como función de la concentración de
oxígeno c(mg/L). Se sabe que dichos datos pueden modelarse por medio de la siguiente
ecuación:
𝑘/(0 𝑐 3
𝑘=
𝑐4 + 𝑐 3
Donde cs y kmax son parámetros.
c 0.5 0.8 1.5 2.5 4
k 1.1 2.4 5.3 7.6 8.9
a) (1.0) Estime cs y kmax y pronostique la tasa de crecimiento para 2 mg/L utilizando el modelo
con los parámetros hallados.
b) (1.0) Utilice un método de interpolación de Lagrange de tercer grado para hallar la tasa de
crecimiento para 1, 2 y 3 mg/L.
3. (2.0) Un reactor esta estratificado térmicamente en la tabla siguiente:
Profundidad, m 0 0.5 1.0 1.5 2.0 2.5 3.0
Temperatura, ºC 70 68 55 22 13 11 10
Al graficar los valores de la tabla, se nota que el tanque puede idealizarse como dos zonas
separadas por un gradiente fuerte de temperatura, o termoclina.
La profundidad de este gradiente se define como el punto de inflexión de la curva
temperatura–profundidad. A esta profundidad, el flujo de calor de la superficie a la capa del
fondo se calcula con la ley de Fourier:
𝑑𝑇
𝐽 = −𝑘
𝑑𝑧
Use un ajuste de estos datos para determinar la profundidad de la termoclina. Si k=0.02
cal/(s*cm*ºC), calcule el flujo a través de esta interfaz