Resolución del Sistema de Lorenz con Método RK4
Contenido del curso
Recordemos ahora que el sistema de Lorenz está dado por el siguiente conjunto de ecuaciones:
Aquí las variables (X₁, X₂, X₃) son proporcionales a (ω, Ta Tb), respectivamente y cada celda de convección que se forma tiene una velocidad característica de rotación ω y una variación en temperatura horizontal y vertical Ta y Tb, respectivamente. Recordemos que el sistema de Lorenz es una simplificación de modelado de un fenómeno que llamamos convección, el cual es el responsable de muchas cosas en nuestra vida diaria como el clima y el agua hirviendo en el fogón de nuestra cocina y ahora en esta clase resolveremos estas ecuaciones usando RK4. El código que escribimos en la clase anterior para implementar un paso del algoritmo RK4 también sirve para resolver un sistema de EDOs, en ese caso tanto y como y₀ serán estructuras tipo array() donde cada elemento corresponde a una de las variables (X₁, X₂, X₃).
Ahora, vamos a implementar una clase que contiene métodos para construir el sistema de EDOs de Lorenz y el uso de RK4 para obtener la solución numérica del mismo. Esta clase contiene lo siguiente:
def f(t, y): return y
Una vez definida esta clase es posible ejecutar la solución numérica unas pocas líneas, pero antes de hacerlo es importante notar que aquí el valor de los parámetros es crucial para determinar el tipo de comportamiento que el sistema tendrá, colocaremos dos ejemplos particulares:
La única diferencia entre ambas soluciones es que el valor del parámetro σ cambia de 5 (a la izquierda) a 10 (a la derecha). En ambos casos graficamos la velocidad de rotación de la celda de convección, y a partir de eso vemos que:
Pero es importante resaltar que ambos casos consideran la misma condición inicial. Y la conclusión de esto es que el sistema puede mostrar un comportamiento regular para ciertas condiciones físicas particulares y para otras puede mostrar un comportamiento turbulento, que es lo que usualmente sucede en la vida real. Incluso, otro aspecto importante de este sistema dinámico es que presenta una propiedad interesante denominada Caos. ¿Qué es el Caos?, bueno para entender esto vamos a comparar dos simulaciones que tengan los mismos parámetros pero cuyas condiciones iniciales sean diferentes tan solo por una valor muy pequeño:
Del gráfico resultante vemos que una pequeña variación en la condición inicial conduce a comportamientos totalmente diferentes y a esta característica de un sistema dinámico la denominamos Caos. En general, el caos es un término que hace referencia a que una pequeña modificación de un sistema implica consecuencias muy grandes. A menudo encontramos referencias de esto en la cultura pop, como aquellas películas del Efecto Mariposa. Si no las has visto y aún no crees estar seguro de lo que significa el caos, recomiendo que las veas para que tengas más conciencia del significado de caos.
Así pues, cuando hay caos presente en un sistema dinámico es extremadamente difícil hacer predicciones de este, porque en la vida real los datos se obtienen a partir de instrumentos que tienen una precisión finita, así que medir con precisión absoluta una condición inicial para la formación de un huracán (por ejemplo) es imposible, y esto ya representa una limitación para poder hacer predicciones del clima en general.
En nuestras próximas clases entraremos a la sección final de nuestro curso donde veremos como podemos usar métodos numéricos para modelar de manera más acertada la propagación de epidemias como el SARS-COV-2, por ahora cerramos esta clase con el notebook definitivo de lo que acabamos de hacer en este link.
Edson Luis Pérez Castillo
Eduardo Monzón
Giancarlo Arturo Cespedes Mas
En la comparativa del sistema al variar las condiciones iniciales hay un pequeño error: la segunda gráfica corresponde a uno de los valores de temperatura (x2) en lugar del valor de la velocidad de rotación (x1). La gráfica correcta sería la siguiente:
Si bien son bastante similares, corresponden a variables diferentes del sistema, es por ello que comento el error de graficado.
Gracias por la corrección.
Recomiendo mucho repasar métodos numéricos para llegar a este punto sin tantas complicaciones. De momento soy el único comentario de esta clase.