Estoy leyendo un libro de física computacional [1] donde se estudia en profundidad el péndulo no lineal accionado. Esta es la ecuación derivada en el libro:
Los autores saben bajo qué condiciones el péndulo cambia a un régimen caótico, probablemente porque utilizaron literatura previa o porque experimentaron numéricamente con este sistema.
Mirando la ODE, ¿podemos predecir la existencia de un comportamiento caótico? En caso afirmativo, ¿es posible saber, cualitativamente, cómo se deben ajustar los parámetros para encontrar este comportamiento caótico?
[1] Física computacional , NJ Giordano y H. Nakanishi, segunda edición, Pearson Prentice Hall, 2006
En general, no. Es posible reconocer un (sistema de) ecuación ( es) diferencial(es) como no lineal (es) simplemente por inspección, pero hay muchos sistemas no lineales no caóticos. El caos es una condición más fuerte (y, desafortunadamente, no bien definida).
Además, muchos (¿la mayoría?) de los sistemas, incluido el famoso atractor de Lorenz , solo exhiben caos bajo ciertas condiciones. El atractor de Lorenz, por ejemplo, sufre una bifurcación en ; por el origen es un atractor estable . No hay forma de saber que la bifurcación está ahí, o en qué región el sistema es caótico, simplemente observando la forma de las ecuaciones, porque las mismas ecuaciones pueden describir tanto una solución no caótica como una caótica.
La respuesta es un claro no. En primer lugar, las EDO no lineales pueden no tener solución para algunas condiciones iniciales o, por el contrario, tener varias soluciones. Para la unicidad y existencia se necesita que la derivada sea continua lo que puede no ser obvio a simple vista para sistemas de varias EDO no lineales.
Una vez garantizada la unicidad y la existencia, en los sistemas caóticos siempre hay parámetros de control. Un sistema ODE no lineal dado tendrá soluciones simples o complejas pero no caóticas para un rango de parámetros de control y puede tener soluciones caóticas para otro rango de parámetros de control.
La ecuación logística aunque no sea una EDO (X(n+1) = µ.Xn.(1-Xn)) mostrará un comportamiento caótico solo para algunos valores de µ. No se conoce una regla general y simple que permita saber si para un sistema EDO no lineal dado existe un conjunto de valores de los parámetros de control para los cuales las soluciones son caóticas.
Sin embargo, es un poco más fácil desde un punto de vista físico. Las órbitas caóticas en física tienen 2 propiedades: divergencia exponencial (llamada sensibilidad a las condiciones iniciales) y disipación que asegura que las órbitas no exploten hasta el infinito. Esto se ve topológicamente como "estiramiento" y "plegamiento" en el espacio de fase. Entonces, si toma un sistema forzado que será disipativo Y sus ecuaciones de movimiento no serán lineales, entonces tendrá una buena oportunidad de encontrar regímenes caóticos.
El caso del caos hamiltoniano (problema del cuerpo gravitacional N) es diferente, ya que no implica estiramientos y plegamientos, sino inestabilidades orbitales en un toro.
Más allá de estos 2 casos en los que se pueden buscar soluciones caóticas, está por supuesto el dominio de las EDP que da lugar a soluciones caóticas espaciotemporales (en contraposición al mero caos temporal para el caso de las EDO).
Incluso si la pregunta no concierne al PDE, la respuesta es obviamente no para ellos también, incluso si el comportamiento caótico es más frecuente en el espacio temporal que en el dominio temporal.
Benjamín Hodgson
usuario11543
Benjamín Hodgson