Taller final-Métodos

computacionales.
22/09/2012
Universidad Del Norte

Presentado por:
Andrea S. Arevalo ,
Maria Jose Mercado,
Lissy Lizarazo,
Julio Torres Trejos
 1. Una compañía aeroespacial se encuentra desarrollando un aditivo nuevo para un
combustible de aeronaves comerciales. El aditivo está compuesto de tres ingredientes: X,
Y y Z. Para un mayor rendimiento, la cantidad total de aditivo debe ser al menos de 6 mL/L
de combustible. Por razones de seguridad, la suma de los ingredientes X y Y altamente
flamables, no debe exceder los 2.5 mL/L. Además, la cantidad del ingrediente X siempre
debe ser mayor o igual a la de Y, y la de Z debe ser mayor que la mitad de la de Y. Si el
costo por mL para los ingredientes X, Y y Z es de $0.05, $0.025 y $0.15, respectivamente,
determine la mezcla de costo mínimo para un litro de combustible.

Con la ayuda de la herramienta de Excel, se desarrollo por medio de solver, estableciendo

como objetivo minimizar la mezcla de costo mínimo para un litro de combustible,
cambiando las celdas de cantidad de cada ingrediente para el aditivo ( como se muestra
claramente indicado en las tablas de Excel), y colocando como restricciones las mostradas
en la tabla.

Cantidad de cada ingrediente del aditivo (cambiando celda).

INGREDIENTES X Y Z RESTRICCIONES
1,25 1,25 3,5
RENDIMIENTO 1 1 1 6 >= 6
SEGURIDAD 1 1 0 2,5 <= 2,5
CANTIDAD 1 -1 0 0 >= 0
-0,5 1 2,875 >= 0

COSTO 0,05 0,025 0,15

MINIMIZAR – Mezcla de costo
mínimo para 1 Lt de 0,61875 Celda objetivo.
combustible
Ver archivo de Excel anexo.

Análisis de resultados: Al utilizar la herramienta de Excel para resolver el ejercicio se puede

observar que, efectivamente la cantidad del ingrediente Z es mucho mayor que los demás
ingredientes que son altamente inflamables, y se observa que al minimizar la suma producto
de los valores y la cantidad de cada ingrediente obtenemos un costo mínimo para 1 Lt de
combustible de $0,62/L lo que al comprarlo con los valores de cada ingrediente llega a ser
razonable.
2. A partir de un estudio de ingeniería del transporte para determinar el diseño apropiado de
pistas para bicicletas, se obtuvieron datos del ancho de las pistas y la distancia promedio
entre las bicicletas y los autos en circulación. Los datos de 9 calles son:

a) Grafique los datos.

b) Ajuste una línea recta a los datos con una regresión lineal. Agregue esta línea a la gráfica.
Esta es la regresión más adecuada para el conjunto de datos. Presente aproximaciones con
otras regresiones y brinde sus comentarios.
c) Si se considera que la distancia promedio mínima de seguridad entre bicicletas y los autos
en circulación es de 2 m, determine el ancho de pista mínimo correspondiente.

3,5

2,5 y = 0,7335x + 0,7167

R² = 0,8374
2
DISTANCIA

1,5 Series1
Lineal (Series1)
1

0,5

0
0 1 2 3 4
ANCHO DE PISTA

c. distancia promedio mínima= 2

La ecuación de la curva que se obtiene del ajuste de los datos es

Remplazando la distancia mínima promedio

( )
 3,5 y = 1,4411ln(x) + 1,2617
R² = 0,8479
3

2,5

DISTANCIA 2

1,5 Series1
Logarítmica (Series1)
1

0,5

0
0 1 2 3 4
ANCHO DE PISTA

c. distancia promedio mínima= 2

La ecuación de la curva que se obtiene del ajuste de los datos es

( )

Remplazando la distancia mínima promedio

( )

3,5

2,5 y = 1,0691e0,3427x
R² = 0,8322
DISTANCIA

1,5 Series1
Exponencial (Series1)
1

0,5

0
0 1 2 3 4
ANCHO DE PISTA
c. distancia promedio mínima= 2

La ecuación de la curva que se obtiene del ajuste de los datos es

Remplazando la distancia mínima promedio

( )

Análisis de resultados: Al graficar los datos Ancho de pista Vs. Distancia, la nube de datos se
ajusta razonablemente a una recta con pendiente positiva y utilizando regresión lineal y
comparándolo con la logarítmica y la exponencial se obtuvieron resultados muy similares,
aunque la regresión lineal sigue siendo la mejor para este conjunto de datos al observar que
su R² es el mas ajustado.

Con respecto al ancho mínimo correspondiente obtenido con cada una de las ecuaciones se
observa que tanto para la logarítmica como para la exponencial los resultados son los mismos,
se opta por elegir el obtenido de la regresión lineal que son 2, 18 m como ancho mínimo
correspondiente por considerarse adecuado.

3. Un flujo desarrollado por completo que pasa a través de un tubo de 40 cm de diámetro

tiene el perfil de velocidad que se presenta en la tabla siguiente. Encuentre la tasa de flujo
volumétrico, Q, con la relación ∫ 0, donde r es el eje radial del tubo, R es el
radio del tubo, y v es la velocidad. Resuelva el problema con dos enfoques diferentes.

r [cm] v [m/s]
0 0,914
2,5 0,89
5 0,847
7,5 0,795
10 0,719
12,5 0,543
15 0,427
17,5 0,204
20 0
a. Ajuste una curva polinomial a los datos de velocidad e intégrela de forma analítica.

Los datos que describen el perfil de velocidad, Tabla 1, de un flujo fueron graficados y
ajustados a una curva polinomial como se muestra en la Gráfica 1.

r [cm] r [m] v [m/s] 1

y = -24,36x2 + 0,3347x + 0,9049
0 0 0,914
0,8 R² = 0,9976
2,5 0,025 0,89
5 0,05 0,847 0,6

Velocidad
7,5 0,075 0,795
10 0,1 0,719 0,4
12,5 0,125 0,543 0,2
15 0,15 0,427
17,5 0,175 0,204 0
20 0,2 0 0 0,05 0,1 0,15 0,2 0,25
-0,2
Tabla 1 Datos de perfil de velocidad Radio

Gráfica 1 Perfil de velocidad y curva de ajuste

La ecuación de la curva que se obtiene del ajuste de los datos es

, donde y que es la variable dependiente representa la velocidad y
donde x es la variable independiente que es el radio del tubo. Explicado lo anterior, la
ecuación queda expresada como:

Ecuación 1

Por otro lado, se tiene la relación ∫ , donde r es el eje radial del tubo, R es el
radio del tubo, y v es la velocidad. La Ecuación 1 que resultó del ajuste se reemplaza en
esta relación y se procede a integrar analíticamente para hallar la tasa de flujo
volumétrico:

∫ ( )

∫ ( )

[ ∫ ∫ ∫ ]
 [ ]( )

[( ) ( )

( ) ]

[( ) ( ) ( )]
( )
⁄

b. Para la integración utilice una aplicación múltiple de la regla de Simpson 1/3.

La regla de simpson 1/3 múltiple se resumen en la siguiente ecuación:

( ) ∑ ( ) ∑ ( ) ( )
( )

Ecuación 2

Donde
es el límite superior
es el límite inferior
es el número de segmentos iguales en los que se dividirá el intervalo

Límites de Función de la Coeficiente regla de Simpson

segmentos Integral 4 si es impar; 2 si es par

1
0,000 0,000 0,000
2 0,050 0,270 1,081
3 0,100 0,436 0,872
4
0,150 0,383 1,532
0,200 -0,004 -0,004
SUMA=m 3,481
Tabla 2 Cálculo de la integral por el método de Simpson 1/3 múltiple

( )
 ( )

c. Encuentre el error porcentual con el uso de la integral del ajuste polinomial como el valor
más correcto.

Análisis de resultados. Desarrollando cada uno de los ítems de este ejercicio se puede concluir
que el método de Simpson como método numérico para la solución de integrales resulta ser muy
preciso al ser comparado con los resultados obtenidos con el desarrollo de la integración analítica.

Como el método de Simpson ajusta cada segmento a un polinomio de segundo grado y el ajuste
de la curva que se ingresa a la integral de Q fue una curva poligonal de orden 2 el grado de
precisión es exacto.

En cuanto a los resultados propios del ejercicio se obtuvo una tasa de flujo volumétrico de
⁄ para la tubería de 40 cm de diámetro.

4. En la investigación de un homicidio o muerte accidental, es importante estimar el tiempo

que ha transcurrido desde la muerte. De observaciones experimentales, se sabe que la
temperatura superficial de un objeto cambia con una tasa proporcional a la diferencia
entre la temperatura del objeto y la del ambiente circundante, o temperatura ambiente
(Ley de Newton del enfriamiento). Así, si T(t) es la temperatura del objeto al tiempo t, y Ta
es la temperatura ambiente constante:

( )

Donde K >0 es una constante de proporcionalidad. Suponga que en el momento t=0 se

descubre un cuerpo y se mide su temperatura, To. Se supone que en el momento de la
muerte, la temperatura del cuerpo, Td, era el valor normal de 37°C. Suponga que la
temperatura del cuerpo al ser descubierto era de 29.5°C, y que dos horas después era de
23.5 °C. La temperatura ambiente es de 20°C.
 a. Determine el valor de K y el tiempo de la muerte.

El primer paso para determinar los valores de k y el t correspondiente al tiempo de la

muerte es resolver la ecuación diferencial del enfriamiento como se muestra en la Tabla 3.

( )

∫ ∫
Solución de la Ecuación Diferencial ( )

( ) ( )

( )

Tabla 3 Solución de la ecuación diferencial de la Ley del Enfriamiento

Teniendo la ecuación de la Temperatura de un cuerpo expresada en función de la

temperatura inicial, la temperatura ambiente, el tiempo trascurrido después de la muerte
y la constante k, con los valores que ofrece el enunciado se procede a calcular las dos
variables que se piden. Este procedimiento se muestra en la Tabla 4:

( )

( ) ( ) ( ) ( )

Tabla 4 Cálculo de k y de t a partir de la solución de la ecuación diferencial

Los valores de k y t son 0.4993 y 1.166 horas respectivamente

 b. Resuelva la EDO en forma numérica y grafique los resultados.

Para resolver la ecuación diferencial ordinaria se empleo el método de Euler para la

solución de ecuaciones diferenciales ordinarias

( )

( )

Donde
( ) es la ecuación diferencial evaluada en y
es el tamaño de paso

La ley del enfriamiento definida por la ecuación diferencial puede ser resuelta por este
método. Haciendo uso del valor de k = 0.4993 calculado en el ítem anterior se obtiene
finalmente una función del tipo ( )

( )  ( ) 

Utilizando como valores iniciales la temperatura de 29.5º con la que fue hallado el cuerpo,
suponiendo este evento con un t=0 y un tamaño de paso de 0.1 en la escala del tiempo se
aplica el método de Euler y se obtienen los resultados que se muestran en la Tabla 5

T(t) t T
T(-1,2)= -1,2 37,05
T(-1,1)= -1,1 36,24
T(-1,0)= -1 35,46
T(-0,9)= -0,9 34,73
T(-0,8)= -0,8 34,03
T(-0,7)= -0,7 33,36
T(-0,6)= -0,6 32,73
T(-0,5)= -0,5 32,12
T(-0,4)= -0,4 31,54
T(-0,3)= -0,3 31,00
T(-0,2)= -0,2 30,47
T(-0,1)= -0,1 29,97
T(0) 0 29,50
T(0,1)= 0,1 29,03
T(0,2)= 0,2 28,58 38
T(0,3)= 0,3 28,15 36
T(0,4)= 0,4 27,74 34

Temperatura (ºC)
T(0,5)= 0,5 27,35
32
T(0,6)= 0,6 26,99
30
T(0,7)= 0,7 26,64
28
T(0,8)= 0,8 26,31
26
T(0,9)= 0,9 25,99
24
T(1)= 1 25,69
Antes del hallazgo 22 Después del
T(1,1)= 1,1 25,41 del cuerpo hallazgo del cuerpo
T(1,2)= 1,2 25,14 20
-1,5 -1 -0,5 0 0,5 1 1,5
Tiempo (horas)

Tabla 5 Solución numérica de la ecuación diferencial de la Ley del Enfriamiento. Condiciones iniciales
T=29.5º en t=0

Ahora, utilizando como valores iniciales la temperatura de 37ªC con la que muere el
cuerpo, suponiendo este evento con un t=0 y un tamaño de paso de 0.1 en la escala del
tiempo, se aplica el método de Euler y se obtienen los resultados que se muestran en la
Tabla 6.

T(t) t T
T(0) 0 37,00
T(0,1)= 0,1 36,15
T(0,2)= 0,2 35,34
T(0,3)= 0,3 34,58
T(0,4)= 0,4 33,85
T(0,5)= 0,5 33,16
T(0,6)= 0,6 32,50
T(0,7)= 0,7 31,88
T(0,8)= 0,8 31,28
T(0,9)= 0,9 30,72
T(1)= 1 30,19
T(1,1)= 1,1 29,68
T(1,2)= 1,2 29,19
T(1,3)= 1,3 28,74
T(1,4)= 1,4 28,30
T(1,5)= 1,5 27,88
T(1,6)= 1,6 27,49
 T(1,7)= 1,7 27,12 38,00
T(1,8)= 1,8 26,76 36,00
T(1,9)= 1,9 26,42
34,00
T(2,0)= 2 26,10

Temperatura (ºC)
T(2,1)= 2,1 25,80 32,00
T(2,2)= 2,2 25,51 30,00
T(2,3)= 2,3 25,23 28,00
T(2,4)= 2,4 24,97 26,00
24,00 Enfriamiento del cuerpo
desde su muerte
22,00
20,00
0 0,5 1 1,5 2 2,5
Tiempo (horas)

Tabla 6 Solución numérica de la ecuación diferencial de la Ley del Enfriamiento. Condiciones iniciales T=37º
en t=0

Análisis de resultados. Con el método de Euler se pudo desarrollar la ecuación diferencial

ordinaria de la ley del enfriamiento de una manera muy simple, uno de los cuidados que se
tuvieron presentes en la aplicación del método fue escoger una longitud de tamaño de paso
bastante pequeña para que los errores fuesen razonablemente pequeños. Otro paso que marca el
proceso con este método es definir las condiciones iniciales; como el problema nos brinda la
temperatura con la que muere el cuerpo, la temperatura con la que es hallado, la temperatura
ambiente y por otro lado se calculo el tiempo de muerte, se pueden tomar como condición inicial
el momento exacto de la muerte o el momento en que es hallado el cuerpo. El método fue
aplicado para ambos casos y los resultados son muy parecidos

5. El flujo a través de medios porosos queda descrito con la ecuación de Laplace:

donde h es la carga. Use métodos numéricos para determinar la distribución de la carga
para el sistema que se muestra en la figura.
Para poder determinar la distribucion de la carga en la placa, fue necesario dibidirlo en
forma de cuadricula como se ve en la siguente figura:

Se piensa determinar la carga de cada uno de estos nodos para y con estos valores se
puede tener una apreciación de como esta distribuida la carga en toda la placa.

Para hallar la carga en los nodos se utilizará la siguiente ecuación en forma general:

En puntos donde no haya valores de Ti-1 o de Tj-1 estos valores serán remplazados
utilizando la siguiente ecuación:

Que se puede entender mejor utilizando la siguiente gráfica:

Ya sabiendo esto, se halló una ecuación para cada uno de los 9 puntos mostrados en la
gráfica inicial, de los que no se conoce su carga, y así tener un sistema de nueve
ecuaciones con nueve incógnitas que fueron resueltas utilizando la herramienta solver
del programa Excel.

Al resolver estas ecuaciones se obtuvieron las siguientes cargas para cada uno de los
puntos de la placa (Las ecuaciones están planteadas y resueltas en el archivo de Excel de
este punto adjunto al trabajo):

 (1,1): 16,228
 (1,2): 16,297
 (1,3): 16,337
 (2,1): 17,159
 (2,2): 17,311
 (2,3): 17,377
 (3,1): 17,785
 (3,2): 18,412
 (3,3): 18,550

Y como ya se conoce la carga de los puntos (4,2) y (4,3) es de 20.

Análisis de resultados: De los resultados anteriores podemos apreciar como la

distribución de la carga en la placa va disminuyendo de derecha a izquierda y de arriba
hacia abajo, siendo el punto (1,1) el que tiene una carga menor siendo de 16,228
mientras que el punto (3,3) es el de mayor carga (sin tener en cuenta los puntos 4,2 y
4,3 cuya carga era conocida desde un principio) siendo esta de 18,550.