Practica 2

COMPUTACIN Y CLCULO NUMRICO (F-314)
Ejercitacin N 2: Ecuaciones diferenciales ordinarias

Mtodos numricos para ingenieros (5 edicin)
S.C. Chapra y R.P. Canale
Ed.: McGraw-Hill Interamericana
Adems de los ejercicios listados a continuacin, se recomienda analizar los ejercicios
resueltos a modo de ejemplo en el libro, resolver al menos uno de los incluidos en la seccin
Estudio de Casos y aquellos planteados durante las clases tericas.
Problema 1: Resuelva en forma analtica el problema de valores iniciales siguiente, en el intervalo
de x=0 a 2:
donde y(0) = 1. Grafique la solucin.

Problema 2: Utilice el mtodo de Euler con h = 0.5 y 0.25, para resolver el problema 1. Grafique
los resultados en la misma grfica para comparar en forma visual la exactitud de los dos tamaos
de paso.
Problema 3: Calcule la velocidad de un paracaidista en cada libre con el empleo del mtodo de
Euler para el caso en que m=80 kg y c=10 kg/s. Lleve a cabo el clculo desde t=0 hasta t=20 s con
un tamao de paso de 1s. Use una condicin inicial en que el paracaidista tiene una velocidad hacia
arriba de 20 m/s en t=0. Suponga que el paracadas se abre instantneamente en t=10 s, de modo
que el coeficiente de arrastre sube a 50 kg/s.
Problema 4: La ley de Newton del enfriamiento establece que la temperatura de un cuerpo cambia
con una tasa que es proporcional a la diferencia de su temperatura y la del medio que lo rodea
(temperatura ambiente).
donde T=temperatura del cuerpo (C), t=tiempo (min), k=constante de proporcionalidad (por
minuto), y Ta=temperatura del ambiente (C). Suponga que una tasa de caf tiene originalmente una
temperatura de 68C. Emplee el mtodo de Euler para calcular la temperatura desde t=0 hasta 10
min, usando un tamao de paso de 1 min, si Ta=21C y k=0.017min-1.
Problema 5: Emplee el mtodo de Heun con h=0.5 para resolver el problema 1. Itere el corrector
hasta obtener s= 1%.
Problema 6: Emplee el mtodo del punto medio con h=0.5 y 0.25, para resolver el problema 1.
Problema 7: Resuelva la ecuacin que se presenta a continuacin, de t=0 a 3, con h=0.1, con los
mtodos de a) Heun (sin corrector), y b) RK y Ralston de segundo orden:
Problema 8: Use los mtodos de a) Euler, y b) RK de cuarto orden, para resolver:
en el rango de x=0 a 1, con un tamao de paso de 0.2, con y(0) = 2, y z(0) = 4.

Problema 9: Implementar en cdigos Fortran, la solucin del problema 4, usando los mtodos de
integracin de Euler, Heun y Runge Kutta de cuarto orden. Para varios valores del paso de
integracin con cada mtodo, realizar slo una iteracin utilizando varios pasos h, y graficar el
logaritmo del error en funcin del logaritmo de h. Interpretar los resultados obteinidos.
Problema 10: El movimiento de un sistema acoplado
masa resorte (ver la figura) est descrito por la ecuacin
diferencial ordinaria que sigue:
donde x=desplazamiento desde la posicin de equilibrio (m), t=tiempo (s), m=20 kg masa, y c=
coeficiente de amortiguamiento (Ns/m). El coeficiente de amortiguamiento c adopta tres valores, 5
(subamortiguado), 40 (amortiguamiento crtico), y 200 (sobreamortiguado). 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
ecuacin con el uso de un mtodo numrico durante el periodo de tiempo 0 t 15 s. Grafique (en
una misma grfica) el desplazamiento versus el tiempo para cada uno de los tres valores del
coeficiente de amortiguamiento.
Problema 11: Dada
a) Estime el tamao de paso requerido para mantener la estabilidad con el uso del mtodo de Euler
explcito.
b) Si y(0)=0, utilice el mtodo de Euler implcito para obtener la solucin desde x=0 hasta 2, con un
tamao de paso de 0.1.
Problema 12: Dado que
Si y(0)=1, emplee el mtodo de Euler implcito para obtener una solucin de t=0 a 4, con un tamao
de paso de 0.4.
Problema 13: El balance de calor de estado estacionario de una barra se representa como:
Obtenga una solucin analtica para una barra de 10 m con T(0) =240 y T(10) = 150.
Problema 14: Use el mtodo del disparo para resolver el problema 13.
Problema 14: Emplee el mtodo del disparo para resolver:
con las condiciones de frontera y(0) = 5 y y(20) = 8.

Problema 15: Resuelva para la EDO no dimensionada, por medio del mtodo de diferencias finitas,
que describa la distribucin de la temperatura en una barra circular con fuente interna de calor S.
en el rango 0 r 1, con las condiciones de frontera
para S=1, 10 y 20 K/m2. Grafique la temperatura en funcin del radio.

Problema 16: Una barra calentada con una fuente de calor uniforme se puede modelar con la
ecuacin de Poisson:
Dada una fuente de calor f(x)=25 y las condiciones en la frontera T(x=0)=40 y T(x=10)=200,
resuelva para la distribucin de temperatura con a) el mtodo del disparo, y b) el mtodo de
diferencias finitas (x=2).
Conectarse al cluster (skynet):
- Si se encuentra fuera de las dependencias de la FCEIA (por ejemplo desde su casa, o algn bar),
primero debe conectarse a la red de FCEIA, con el siguiente comando, desde una teminal:
ssh [email protected]
password; F314ecen
- Una vez conectado a la red de FCEIA, o en caso de estar en la red de FCEIA (por ejemplo desde
Felix, o el edificio del reactor), ejecutar:
ssh [email protected]
password: ccnfis-ecen

Practica 2

Cargado por

Copyright:

Formatos disponibles

Practica 2

Cargado por

Información del documento

Descripción 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

Practica 2

Cargado por

Copyright:

Formatos disponibles

COMPUTACIN Y CLCULO NUMRICO (F-314)

Ejercitacin N 2: Ecuaciones diferenciales ordinarias

donde y(0) = 1. Grafique la solucin.

Problema 8: Use los mtodos de a) Euler, y b) RK de cuarto orden, para resolver:

en el rango de x=0 a 1, con un tamao de paso de 0.2, con y(0) = 2, y z(0) = 4.

con las condiciones de frontera y(0) = 5 y y(20) = 8.

en el rango 0 r 1, con las condiciones de frontera

para S=1, 10 y 20 K/m2. Grafique la temperatura en funcin del radio.

También podría gustarte