Universidad de Sucre & Departamento de Matemáticas

Curso: Métodos Numéricos, septiembre 2020

Guía de clase 1. Modelación matemática y métodos numéricos

Modelación matemática: es un proceso dinámico en el que se identifican los factores más relevantes de un fenómeno natural o artificial, los cuales
se expresan en una relación funcional que debe ser contrastada con la realidad, para realizar ajustes que permitan una mejor predicción del
fenómeno.

El proceso de modelación matemática no termina con la formulación funcional de la relación entre factores, 𝑦 = 𝑓(𝑥1 , 𝑥2 , … , 𝑥𝑛 ), es cíclico por
cuanto siempre surgen nuevos datos y/o contextos que permiten generalizar el modelo matemático.

Los modelos matemáticos incorporan conocimientos científicos como las leyes generales de la naturaleza y también conocimientos tecnológicos,
por ejemplo, mediciones de campo, propiedades de los suelos o la composición química de un lugar, ete.
Un modelo matemático de un fenómeno es esquema simplificado e idealizado de aquel, constituido por símbolos y operaciones matemáticas. Es
una traducción de la realidad física en términos matemáticos, es decir, una forma de representar cada uno de los factores que intervienen en un
cierto proceso físico mediante objetos matemáticos y sus interrelaciones.

Ejemplo 1 Modelo matemático para el periodo de oscilación de un péndulo.

Ejemplo 2 Modelo matemático para la distribución de la temperatura sobre una placa conductora de calor, en régimen estacionario.

Para modelar el proceso de distribución de temperatura se tiene que resolver una ecuación diferencial parcial; esto puede hacerse de forma
analítica ó numérica, mediante un sistema de ecuaciones, 𝐴𝑈 = 𝑏, que puede tener muchas incógnitas. En la secundaria a veces resolvíamos
sistemas de ecuaciones de tres incógnitas y hacer las operaciones para resolverlas era tedioso, como aquí son muchísimas incógnitas, desarrollamos
un método numérico, que implementado en una computadora, permite obtener una buena aproximación a la matriz 𝑈.
Para que la computadora haga este trabajo, se desarrollan algoritmos y software. Por eso se habla de modelación matemática y computación
científica.

Modelación matemática y computacional: la respuesta al anhelo ancestral de predecir a la naturaleza.

Solución del problema 1.3: en vez de la relación lineal de la ecuación (1.7), elija modelar la fuerza hacia arriba sobre el paracaidista como una
relación de segundo orden, 𝐹𝑈 = −𝑐 ′ 𝑣 2 , donde 𝑐 ′ =coeficiente de arrastre de segundo orden (kg/m), a. Con el empleo del cálculo, obtenga la
solución de forma cerrada para el caso en que al inicio el saltador se encuentra en reposo, (𝑣 = 0 en 𝑡 = 0), b. Repita el cálculo numérico del
ejemplo 1.2 con los mismos valores de condición inicial y de parámetros. Utilice un valor de 0.225 kg/m para 𝑐 ′ .

Solución analítica:

𝑑𝑣 𝑚𝑔−𝑐 ′ 𝑣 2 𝑐′ 𝑑𝑣
𝑎(𝑡) = = =𝑔− 𝑣2 ↔ 𝑐′
= 𝑑𝑡
𝑑𝑡 𝑚 𝑚 𝑔− 𝑣 2
𝑚

𝑚
𝑐 ′ 𝑑𝑣 = 𝑐 ∫
𝑣 𝑡 𝑣 ′ 𝑣
𝑑𝑣 𝑑𝑣 𝑚 𝑣 𝑑𝑣 𝑚 𝑣
∫ ′
= ∫ 𝑑𝑡 = 𝑡 ↔ ∫ 𝑚𝑔 𝑚𝑔 = ′
∫ = arctanh( )
0 𝑐 2 0 0 − 𝑣2 𝑚 0 ′ − 𝑣2 𝑐 0 𝑚𝑔 2 𝑚𝑔 𝑚𝑔
𝑔 −𝑚𝑣 𝑐 ′ 𝑐 (√ ′ ) − 𝑣 2 ′
2𝑐 √ ′ √ ′
𝑐 𝑐 𝑐

𝑚𝑔 𝑚𝑔 𝑚𝑔
2𝑐 ′ √ 𝑐 ′ 𝑣 2𝑐 ′ √ 𝑐 ′ 𝑣 𝑚𝑔 2𝑐 ′ √ 𝑐 ′
𝑡 = arctanh ↔ tanh 𝑡 = → 𝑣(𝑡) = √ ′ tanh 𝑡
𝑚 𝑚𝑔 𝑚 𝑚𝑔 𝑐 𝑚
√ √ ′
( 𝑐′ ) ( ) 𝑐 ( )

Para los valores, 𝑚 = 68.1, 𝑐 ′ = 0.225, 𝑔 = 9.81, 𝑣(𝑡) = 54.490tanh(0.360 ∗ 𝑡), modelo matemático para el paracaidista.
 𝑑𝑣 ∆𝑣 𝑣(𝑡𝑖+1 )−𝑣(𝑡𝑖 ) 𝑚𝑔−𝑐 ′ 𝑣 2 (𝑡𝑖 ) 𝑐′
Solución numérica: a(t)= ≈ = = =𝑔− (𝑣 (𝑡𝑖 ))2
𝑑𝑡 ∆𝑡 𝑡𝑖+1 −𝑡𝑖 𝑚 𝑚

𝑐′ 2
Esta ecuación se reordena para obtener 𝑣 (𝑡𝑖+1 ) = 𝑣 (𝑡𝑖 ) + [𝑔 − (𝑣(𝑡𝑖 )) ] (𝑡𝑖+1 − 𝑡𝑖 ), 𝑖 = 0,1,2, . . , 𝑁
𝑚

Ahora, al empezar los cálculos con 𝑡(0) = 0, la velocidad del paracaidista es igual a cero, 𝑣 (0) = 0 y (𝑡𝑖+1 − 𝑡𝑖 ) = 2, se obtienen los
siguientes valores de velocidad
0.225
𝑡1 = 2, 𝑣 (2) = 𝑣(0) + [9.81 − 𝑣(0)2 ] (2) = 0.000 + 19.620 = 19.620
68.1
0.225
𝑡2 = 4, 𝑣(4) = 𝑣(2) + [9.81 − 𝑣(2)2 ] (2) = 19.620 + 17.076 = 36.696
68.1
0.225
𝑡3 = 6, 𝑣(6) = 𝑣(4) + [9.81 − 𝑣(4)2 ] (2) = 36.696 + 10.722 = 47.418
68.1
Los valores restantes se obtienen de modo similar.

Estos resultados pueden obtenerse de modo más fácil y en menor tiempo, mediante un programa informático, codificado en un software para
computo científico. A continuación, se muestra el programa VelocPar, codificado en el software Scilab, que permite obtener los valores de tiempo
y velocidad de caída del paracaidista.
Cuestión: El modelo anterior solo considera la velocidad de caída si abrir el paracaídas, ¿Cómo puede estimarse el efecto del paracaídas sobre la
velocidad de caída?
Actividades para el estudiante

1. En el problema 1.8, página 17 del libro guía, determine: a. el modelo matemático asociado al proceso de transferencia de calor, b. calcule la
temperatura final del salón al cabo de 15 minutos.

𝑄 𝑃𝑉𝑀𝑤𝑡
Solución: a. Modelo matemático, ∆𝑇 = 𝑇𝑓 − 𝑇𝑖 = , dado que 𝑄 = 𝑁𝑒 𝑞𝑒 𝑡 y 𝑚 = , entonces,
𝑚𝐶𝑣 𝑅𝑇

𝑁𝑒 𝑞𝑒 𝑡 𝑅𝑇 𝑁𝑒 𝑞𝑒 𝑡 𝑅𝑇
∆𝑇 = 𝑇𝑓 − 𝑇𝑖 = ( ) 𝑃𝑉𝑀 = (0.001 60) 𝑃(𝑉 −𝑁 , ∆𝑇, en °K ó °C
𝐶𝑣 𝑤𝑡 𝐶𝑣 𝑠 𝑒 𝑉𝑒 )𝑀𝑤𝑡

Los parámetros son:

𝑁𝑒 = 35, 𝑞𝑒 = 80, 𝐶𝑣 = 0.718, 𝑡 = 20 ∗ 60 = 1200, 𝑅 = 8.314, 𝑇 = 20 + 273.15 = 293.15, 𝑃 = 101.325, 𝑉𝑠 = 11 ∗ 8 ∗ 3 = 264, 𝑉𝑒 = 0.075, 𝑀𝑤𝑡 = 28.97

𝑇𝑓 = 𝑇𝑖 + ∆𝑇 = 20 + 14.866 = 34.866 °𝐶
 1. Obtenga una solución analítica y numérica del problema 1.4, propuesto en la página 16 del libro guía.

𝑑𝑣 ∆𝑣 𝑣(𝑡𝑖+1 )−𝑣(𝑡𝑖 ) 𝑚𝑔−𝑐𝑣(𝑡𝑖 ) 𝑐

Solución numérica: ≈ = = =𝑔− 𝑣(𝑡𝑖 )
𝑑𝑡 ∆𝑡 𝑡𝑖+1 −𝑡𝑖 𝑚 𝑚
𝑐
Esta ecuación se reordena para obtener 𝑣 (𝑡𝑖+1 ) = 𝑣 (𝑡𝑖 ) + [𝑔 − 𝑣(𝑡𝑖)] (𝑡𝑖+1 − 𝑡𝑖 ), 𝑖 = 0,1,2, . . , 𝑁
𝑚

Paracaidista 1:
Empezar los cálculos con 𝑡0 = 0, la velocidad del paracaidista es igual a cero, 𝑣 (0) = 0 y (𝑡𝑖+1 − 𝑡𝑖 ) = 1, por lo que se obtiene la
siguiente función de velocidad
12
𝑣(𝑡𝑖+1 ) = 𝑣 (𝑡𝑖 ) + [𝑔 − 𝑣(𝑡𝑖 )] (𝑡𝑖+1 − 𝑡𝑖 ), 𝑖 = 0,1,2, . . , 𝑁
70
12
𝑡1 = 1, 𝑣 (2) = 𝑣(0) + [9.81 − 𝑣(0)] (1) = 9.810
70
12
𝑡2 = 2, 𝑣 (4) = 𝑣(2) + [9.81 − 𝑣(1)] (1) = 17.938
70
12
𝑡3 = 3, 𝑣 (6) = 𝑣(4) + [9.81 − 𝑣(2)] (1) = 24.673
70
12
𝑡4 = 2, 𝑣(4) = 𝑣(2) + [9.81 − 𝑣(3)] (1) = 30.253
70

Los valores restantes se obtienen de modo similar.

Paracaidista 2:

Empezar los cálculos con 𝑡0 = 0, la velocidad del paracaidista es igual a cero, 𝑣 (0) = 0 y (𝑡𝑖+1 − 𝑡𝑖 ) = 1, por lo que se obtiene la
siguiente función de velocidad
15
𝑣(𝑡𝑖+1 ) = 𝑣 (𝑡𝑖 ) + [𝑔 − 𝑣(𝑡𝑖 )] (𝑡𝑖+1 − 𝑡𝑖 ), 𝑖 = 0,1,2, . . , 𝑁
80
15
𝑡1 = 1, 𝑣 (1) = 𝑣(0) + [9.81 − 𝑣(0)] (1) = 9.810
80
15
𝑡2 = 2, 𝑣 (2) = 𝑣(1) + [9.81 − 𝑣(1)] (1) = 17.781
80
15
𝑡3 = 3, 𝑣 (3) = 𝑣(2) + [9.81 − 𝑣(2)] (1) = 24.257
80
15
𝑡4 = 4, 𝑣 (4) = 𝑣(3) + [9.81 − 𝑣(3)] (1) = 29.519
80

Los valores restantes se obtienen de modo similar.

Estos resultados pueden obtenerse de modo más fácil y en menor tiempo, mediante el programa VelocPar, codificado en el software Scilab, que
permite obtener los valores de tiempo y velocidad de caída del paracaidista, para diferentes valores de 𝑚 y 𝑐.
A los 9 segundos, el paracidista 1 ha adquirido, una velocidad, 𝑉 = 46.69 m/seg.
Entonces el paracaidista 2, requiere aproximadamente, 𝑡 = 10.8 seg., para adquirir una velocidad 𝑉 ≈ 46.69 m/seg.
P 1.10 Para el mismo tanque de almacenamiento que se describe en el problema 1.9, supónga qu el flujo de salida no es constante sino que la tasa
depende de la profundidad. Para este caso, la ecuación diferencial para la profundidad puede escribirse como

Use el método de Euler para resolver cuál sería la profundidad 𝑦, desde 𝑡 = 0 hasta 10 d, con un tamaño de paso de 0.5 d. Los valores de los
parámetros son 𝐴 = 1250 m2, 𝑄 = 450 m3/d, y 𝛼 = 150. Suponga que la condición inicial es 𝑦 = 0.

Solución:

𝑑𝑦 𝑄 𝛼 (1 + 𝑦)1.5 1
= 3 𝑠𝑒𝑛2 (𝑡) − → 𝑑𝑦 = (3𝑄𝑠𝑒𝑛2 (𝑡) − 150(1 + 𝑦)1.5 )𝑑𝑡
𝑑𝑡 𝐴 𝐴 𝐴
1
Para valores de 𝑡𝑖+1 y 𝑡𝑖 cercanos, 𝑑𝑦 ≈ ∆𝑦: = 𝑦𝑖+1 − 𝑦𝑖 = (3𝑄𝑠𝑒𝑛2 (𝑡𝑖 ) − 150(1 + 𝑦𝑖 )1.5 )(∆𝑡 = 𝑡𝑖+1 − 𝑡𝑖 ),
𝐴

1
𝑦𝑖+1 = 𝑦𝑖 + (3𝑄𝑠𝑒𝑛2 (𝑡) − 150(1 + 𝑦𝑖 )1.5 )(∆𝑡), 𝑖 = 0,1,2, … , 𝑁
𝐴

Seguidamente utilizamos el siguiente código para obtener valores de 𝑦(𝑡):

Actividad 2, para el estudiante: resuelva de modo analítico y numérico, los ejercicios 1.13 y 1.14, propuestos en la página 17 del libro guía.

Bibliografía

1. Chapra Steven, R. Canales. Métodos numéricos para ingenieros, edición 7 ma. Editorial Mc. Graw-Hill.
2. María Lucía Brito-Vallina, Isidro Alemán-Romero, Elena Fraga-Guerra, José Luís Para-García, Ruth Irene Arias-de Tapia. Papel de la modelación
matemática en la formación de los ingenieros. Ingeniería Mecánica vol.14 no.2 La Habana Mayo-ago. 2011
3. http://ciencia.unam.mx/leer/115/Modelacion_matematica_y_computacional_la_respuesta_al_anhelo_ancestral_de_predecir_a_la_nat
uraleza
4. C. Henry Edward, David E. Penney. Ecuaciones diferenciales y problemas con valores en la frontera, cómputo y modelado. Edición 4 ta.,
editorial Prentice-Hall.