Ada 1 2022

MÉTODOS NUMÉRICOS AVANZADOS
Actividad de aprendizaje 1. Compendio de ejercicios de la unidad I. Solución numérica de

ecuaciones diferenciales ordinarias y diferenciales parciales.
Resultado de Resuelve problemas de ingeniería que involucren encontrar la solución

aprendizaje numérica de ecuaciones diferenciales ordinarias con valores iniciales o
valores en la frontera de ecuaciones diferenciales parciales, utilizando
herramientas computacionales.
Equipos La actividad se realizará en equipos de 5 estudiantes.
Metodología A partir de la información brindada por el profesor, recursos de la
plataforma TEAMS y consulta de la bibliografía pertinente.
Recursos A continuación, se presentan algunos recursos que se pueden consultar:
-Mathews JH, Kurtis DF. (2007). Métodos numéricos con MATLAB.
Madrid: Pearson Educación.
- Chapra SC, Canale RP. (2007). Métodos numéricos para ingenieros. 5ª
ed. México: McGraw-Hill.
-Nakamura, S.(1992). Métodos numéricos aplicados con software.
México. Prentice-Hall Hispanoamericana.
Formato El trabajo deberá contener los siguientes apartados:
-Presentación
-Índice
-Antecedentes (Referente a los métodos utilizados)
-Metodología (Código de Matlab con breve explicación de que es lo que
hace, detalle de funciones de Matlab que utilice)
-Resultados y comentarios de cada ejercicio (Tablas, gráficas, etc.)
-Conclusiones del trabajo (a partir de los resultados obtenidos)
-Referencias Bibliográficas
Formato Word o PDF.
Fecha de Entrega Lunes 26 de Septiembre de 2022 (11:59 pm)
Criterios de - Contenido, código Matlab/Octave y resultados y conclusiones
evaluación -Se evaluará en una escala de 0 100.
- Trabajos duplicados, serán cancelados para ambos equipos.
-Las entregas posteriores a la fecha de entrega tienen una penalización
de 5% por día de atraso.
TEMA I. Problemas de valor inicial
1. El movimiento de un sistema acoplado masa resorte (véase la figura 1) está descrito por la
ecuación diferencial ordinaria que sigue:
Figura 1. Sistema
masa-resorte
donde x = desplazamiento desde la posición de equilibrio (m), t = tiempo (s), m = 20 kg masa, y c =

coeficiente de amortiguamiento (N · s/m). El coeficiente de amortiguamiento c adopta tres valores,
5 (sub amortiguado), 40 (amortiguamiento crítico), y 200 (sobre amortiguado). La constante del
resorte es k = 20 N/m. La velocidad inicial es de cero y el desplazamiento inicial es x = 1 m. Resuelva
esta ecuación con el uso de un método numérico durante el periodo de tiempo 0 ≤ t ≤ 15 s. Grafique
el desplazamiento versus el tiempo para cada uno de los tres valores del coeficiente de
amortiguamiento sobre la misma curva.
2. A model for a batch reactor in which penicillin is produced by fermentation has been derived
as follows (Constantantinides et al. 1970) for cell production and penicillin synthesis,
respectively:
where
y1 = dimensionless concentration of cell mass
y2 = dimensionless concentration of penicillin
t = dimensionless time, 0 ≤ t ≤ 1
Experiments have determined that: b1=13.1, b2= 0.94, b3=1.71
In all cases, experiment with the time step to assure an accurate solution. Also, graph y1 and y2 versus time
with appropriate annotations.
a) Solve this ODE-IVP using the Euler, Heun and RK method and determine the % error of the results
obtain for both methods
b) In all cases, experiment with the time step to assure an accurate solution. Also, graph y1 and y2
versus time with appropriate annotations.
.
3. Suponga que la siguiente reacción química se lleva a cabo en un reactor agitado continuo
(CSTR):
Donde las constantes de velocidad son las siguientes:
k1 = 1 min−1, k2 = 0 min−1, k3 = 2 min−1, k4 = 3 min−1
Las condiciones iniciales en el reactor son (mol/L):
CA0 = 1 CB0 = 0 CC0 = 0
Un balance de masa en estado inestable en cada componente conduce al siguiente conjunto de

EDO:
a) Resuelva el siguiente sistema de ecuaciones por el método de Euler para un valor de

h=0.01
b) Por el método de Runge-kutta 4 orden
4. A radioactive material (A) decomposes according to the series reaction:
where k1, and k2, are the rate constants and B and C are the intermediate and final products,
respectively. The rate equations are
where CA, CB, and CC, are the concentration of materials A, B, and C, respectively. The values of the
rate constants are: k1= 3 s-1 and k2= 1 s-1
a) At time t =1 s and t = 10 s, what are the concentrations of A, B, and C. Compare the result apply
Euler and Runge Kutta method.
b) Sketch the concentration profiles for A, B, and C.
5. The “Monod” model for bioreaction kinetics can be expressed as:
Where:
s = Growth limiting substrate concentration (ML−3)
x = Biomass concentration (ML−3)
k = Maximum specific uptake rate of the substrate (T−1) = 5
ks = Half saturation constant for growth (ML−3) = 20
y = Yield coefficient (MM−1) = 0.05
b = Decay coefficient (T−1) = 0.01
Initial conditions are s0 = 1000 and xo = 100.
a) Solve this ODE-IVP using the Euler method and RK method and determine the % error of
the results obtain for both methods
Note. An appropriate maximum time over which to integrate these equations must be determined
by experimentation. Also, experiment with the time step to guarantee good results. Graph s and x
versus time with appropriate annotations (do this with data for only an appropriate time step).
6. Método de newton para sistemas no lineales. Considere que las siguientes reacciones:
A + B →C
A+B→D
que se llevan a cabo en un reactor CSTR, con las siguientes condiciones(Tabla 2.4). Las leyes de
velocidad para las reacciones son:
Las ecuaciones derivadas de los balances por componentes, para el reactor CSTR, en el estado
estacionario son:
TEMA II. PROBLEMAS CON VALOR EN LA FRONTERA
7. La temperatura de distribución en una placa es descrita por la siguiente ecuación diferencial

ordinaria:
Las condiciones fronteras son T(x=0)=300 K y T(x=10)=400 K. La constante h’=0.05m2 es la constante

de transferencia de calor que representa las perdidas por conducción/convección y σ’=2.7x109 K-3
es la constante que considera las perdidas por radiación/convección. La temperatura ambiente(Ta)
es de 200 K. Resuelva por el método del disparo(ecuación no lineal) y por el método de diferencias
finitas(Δx=0.5) para un sistema de ecuaciones no lineales(usando el método de Newton). Puede
emplear recursos de Matlab como ode45, fzero, etc.
8. Utilice diferencias finitas para resolver la ecuación diferencial ordinaria con valores en la frontera
con condiciones de frontera u(0) = 10 y u(2) = 1. Grafique los resultados de u versus x. Utilice Δx =
0.1.
9. Utilice el método del disparo para solucionar
Obtenga una solución para las condiciones de frontera: T(0) = 200 y T(0.5) = 100. Y el método de
diferencias finitas con dx=1
10. A steady-state heat balance for a rod can be represented as
Obtain a solution for a 10-m rod with T (0) = 240 and T (10) = 150 (a) analytically, (b) with the
shooting method, and (c) using the finite-difference approach with Δ x = 1.
11. An insulated heated rod with a uniform heat source can be modeled with the Poisson equation:
Given a heat source f (x) = 25 °C/m2 and the boundary conditions T (x = 0) = 40 °C and T (x = 10) =
200 °C, solve for the temperature distribution with (a) the shooting method and (b) the finite-
difference method (Δ x = 2).
TEMA III. ECUACIONES DIFERENCIAS PARCILAES
12. Con el método de Crank-Nicholson calcule la distribución de temperatura en una barra larga y
delgada que tiene una longitud de 10 cm y los siguientes valores: k′ = 0.49 cal/(s * cm * °C), Δx
= 2 cm y Δt = 0.1 s. En t = 0, la temperatura de la barra es cero, y las condiciones de frontera se
fijan para todos los tiempos en T(0)= 100°C y T(10) = 50°C. Considere que la barra es de aluminio
con C = 0.2174 cal/(g * °C) y r = 2.7 g/cm3. Por lo tanto, k = 0.49/(2.7 * 0.2174) = 0.835 cm2/s y
λ = 0.835(0.1)/(2)2 = 0.020875
Resuelva en un rango de t= [0,12] y grafique los puntos específicos para los tiempos de t=[3, 6, 9,
12]
13. Use el método de Liebmann para resolver cuál sería la temperatura de la placa cuadrada
calentada que se ilustra en la figura 1, pero con la condición de frontera superior incrementada a
120º y la frontera izquierda disminuida a 60ºC. Utilice un factor de relajamiento de 1.2 para iterar a
es = 1%.
b. Repita nuevamente el ejercicio, pero emplee 49 nodos interiores (es decir, dx = dy = 5
cm).
c. Replique lo realizado en el b) pero considerando que el extremo inferior está aislado.

Ada 1 2022

Cargado por

Copyright:

Formatos disponibles

Ada 1 2022

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ada 1 2022

Cargado por

Copyright:

Formatos disponibles

MÉTODOS NUMÉRICOS AVANZADOS

Actividad de aprendizaje 1. Compendio de ejercicios de la unidad I. Solución numérica de

Resultado de Resuelve problemas de ingeniería que involucren encontrar la solución

donde x = desplazamiento desde la posición de equilibrio (m), t = tiempo (s), m = 20 kg masa, y c =

Donde las constantes de velocidad son las siguientes:

k1 = 1 min−1, k2 = 0 min−1, k3 = 2 min−1, k4 = 3 min−1

Las condiciones iniciales en el reactor son (mol/L):

CA0 = 1 CB0 = 0 CC0 = 0

Un balance de masa en estado inestable en cada componente conduce al siguiente conjunto de

a) Resuelva el siguiente sistema de ecuaciones por el método de Euler para un valor de

4. A radioactive material (A) decomposes according to the series reaction:

5. The “Monod” model for bioreaction kinetics can be expressed as:

Initial conditions are s0 = 1000 and xo = 100.

7. La temperatura de distribución en una placa es descrita por la siguiente ecuación diferencial

Las condiciones fronteras son T(x=0)=300 K y T(x=10)=400 K. La constante h’=0.05m2 es la constante

9. Utilice el método del disparo para solucionar

10. A steady-state heat balance for a rod can be represented as

TEMA III. ECUACIONES DIFERENCIAS PARCILAES

También podría gustarte