🔬 mesoscale physics

HarmonicBalance.jl: A Julia suite for nonlinear dynamics using harmonic balance

HarmonicBalance.jl es un paquete gratuito de Julia que simplifica el análisis de sistemas periódicos no lineales al combinar el método de equilibrio armónico con álgebra simbólica y continuación de homotopía para determinar numéricamente todas las soluciones de estado estacionario y su dinámica de fluctuación, con una interfaz gráfica de usuario completa para explorar las topologías de las soluciones.

Autores originales: Jan Košata, Javier del Pino, Toni L. Heugel, Oded Zilberberg

Publicado 2026-09-17

📖 4 min de lectura☕ Lectura para el café

CC BY 4.0

Autores originales: Jan Košata, Javier del Pino, Toni L. Heugel, Oded Zilberberg

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Muchos sistemas físicos, desde el balanceo de un puente con el viento hasta el flujo de electricidad a través de un circuito, están regidos por reglas que cambian con el tiempo. Cuando estos sistemas son simples, su comportamiento es predecible y constante. Sin embargo, cuando diferentes partes de un sistema interactúan de formas no lineales y complejas, pueden producir resultados sorprendentes. Pueden saltar repentinamente a un nuevo estado, oscilar en patrones inesperados o establecerse en uno de varios estados estables posibles. Los científicos saben desde hace tiempo que estos sistemas no lineales están en todas partes, apareciendo en todo, desde la dinámica de fluidos hasta los procesos biológicos. Sin embargo, predecir exactamente cómo se comportará tal sistema es notoriamente difícil. Los métodos tradicionales suelen basarse en suponer un punto de partida y observar cómo evoluciona el sistema, pero este enfoque tiene un fallo importante: generalmente solo encuentra un resultado, perdiendo los otros estados posibles que podrían existir con la misma validez. Si un sistema puede establecerse en tres patrones estables diferentes, una simulación estándar podría revelar solo uno, dejando la imagen completa incompleta y potencialmente engañosa para ingenieros e investigadores que intentan comprender las verdaderas capacidades del sistema.

Para resolver este problema, un equipo de investigadores ha desarrollado una nueva herramienta de software llamada HarmonicBalance.jl, diseñada para mapear cada uno de los posibles estados estacionarios de estos sistemas complejos y variables en el tiempo. En lugar de observar cómo evoluciona un sistema desde un único punto de partida, la herramienta trata el problema como un rompecabezas para encontrar todas las soluciones posibles a la vez. Funciona asumiendo que el comportamiento del sistema puede describirse como una combinación de ondas repetitivas, o armónicos, que cambian muy lentamente a lo largo del tiempo. Al traducir las complejas ecuaciones móviles del sistema en un conjunto de ecuaciones algebraicas estáticas, el software puede entonces utilizar una poderosa técnica matemática llamada continuación de homotopía. Este método actúa como un puente continuo, partiendo de una versión simple y resoluble del problema y transformándola lentamente en la versión compleja del mundo real. A medida que se transforma, rastrea cada camino que conduce a una solución, asegurando que ningún estado estacionario posible se quede atrás. Esto permite a los investigadores ver todo el panorama de posibilidades, incluyendo estados que podrían ser inestables o difíciles de alcanzar en un experimento real.

El equipo probó este enfoque en varios ejemplos clásicos de la física no lineal, como un oscilador impulsado que responde a fuerzas externas con un giro de no linealidad. En estas pruebas, el software identificó con éxito múltiples estados estables donde los métodos tradicionales habrían encontrado solo uno. Por ejemplo, en un sistema de dos osciladores acoplados, la herramienta reveló un complejo panorama de soluciones, mostrando cómo el sistema podría saltar entre diferentes estados dependiendo de la dirección en la que se cambie un parámetro. Este fenómeno, conocido como histéresis, fue mapeado en detalle, mostrando exactamente dónde el sistema saltaría de un estado a otro. Los investigadores también demostraron que la herramienta podía manejar sistemas cada vez más complejos, como una cadena de cinco osciladores acoplados, encontrando cientos de soluciones distintas en cuestión de minutos en una computadora estándar. Si bien el número de soluciones teóricas puede crecer exponencialmente con el tamaño del sistema, el software filtra eficientemente las imposibles, centrándose en los estados físicos reales que realmente existen.

Más allá de simplemente encontrar estos estados, la herramienta también analiza su estabilidad, determinando qué soluciones son robustas y cuáles colapsarían ante una pequeña perturbación. Incluso puede simular cómo respondería un sistema a cambios lentos en su entorno, imitando la forma en que un experimentador podría girar gradualmente una perilla para recorrer diferentes condiciones. Esta capacidad es crucial para comprender fenómenos como la ruptura espontánea de la simetría, donde un sistema elige un estado sobre otro sin ningún impulso externo. El software está escrito en Julia, un lenguaje conocido por su velocidad y facilidad de uso, y está disponible gratuitamente para la comunidad científica. Al proporcionar una forma de ver la gama completa de comportamientos posibles en sistemas no lineales, esta herramienta ofrece un nuevo nivel de claridad para los investigadores en campos que van desde la ingeniería estructural hasta la óptica. Transforma el estudio de estos sistemas complejos de un juego de azar, donde uno espera tropezar con el estado correcto, en una exploración exhaustiva donde se tiene en cuenta cada posibilidad.

Resumen Técnico de "HarmonicBalance.jl: Una suite de Julia para la dinámica no lineal mediante el equilibrio armónico"

Planteamiento del Problema
Las ecuaciones diferenciales ordinarias (EDO) no lineales gobiernan la evolución temporal de diversos sistemas físicos, desde la dinámica de fluidos y la robótica hasta la optomecánica y las redes de osciladores. Mientras que los sistemas lineales suelen relajarse hacia un estado estacionario único, los sistemas no lineales a menudo exhiben múltiples soluciones de estado estacionario, lo que conduce a fenómenos como la histéresis, la ruptura espontánea de la simetría y el cambio inducido por ruido. Un desafío significativo en el análisis de estos sistemas es que los resolvedores numéricos de EDO estándar (problemas de valor inicial) típicamente encuentran solo un estado estacionario por ejecución, el cual depende de las condiciones iniciales. En consecuencia, una exploración completa del paisaje de soluciones requiere un muestreo infinito del espacio de condiciones iniciales, lo cual es computacionalmente inviable. Además, aunque transformar sistemas forzados en un marco rotatorio puede hacer que los forzamientos sean estacionarios, encontrar todas las raíces de las ecuaciones algebraicas no lineales acopladas resultantes sigue siendo difícil para sistemas multivariables, ya que los algoritmos de búsqueda de raíces estándar (por ejemplo, Newton-Raphson) solo localizan soluciones cercanas a una suposición inicial.

Metodología
Los autores presentan HarmonicBalance.jl, un paquete de Julia diseñado para encontrar todas las soluciones de estado estacionario para sistemas no lineales con forzamiento armónico. La metodología integra tres componentes principales:

Ansatz de Equilibrio Armónico: El paquete asume que, para tiempos suficientemente largos, la respuesta del sistema puede aproximarse mediante una suma finita de términos armónicos con amplitudes de variación lenta. Para un sistema de $N$ EDO de segundo orden con dependencia temporal armónica, la solución $x_i(t)$ se expande como:
$x_i(t) = \sum_{j=1}^{M} [u_{i,j}(T)\cos(\omega_{i,j}t) + v_{i,j}(T)\sin(\omega_{i,j}t)]$
donde $T$ es una escala de tiempo de grano grueso. Al sustituir este ansatz en las EDO originales y filtrar las oscilaciones rápidas (descartando las derivadas temporales de segundo orden de las amplitudes), el problema se reduce a un conjunto de EDO de primer orden autónomas para las variables armónicas $u$ y $v$ . Los estados estacionarios se encuentran igualando a cero las derivadas temporales de estas variables, lo que resulta en un sistema de ecuaciones algebraicas polinómicas acopladas.
Continuación de Homotopía: Para resolver los sistemas polinómicos resultantes, el paquete emplea el método de continuación de homotopía a través de la librería HomotopyContinuation.jl. A diferencia de los métodos iterativos de búsqueda de raíces, la continuación de homotopía deforma un sistema analíticamente soluble (con un número conocido de raíces) hacia el sistema objetivo. Este enfoque garantiza el descubrimiento de todas las raíces complejas. El paquete utiliza un algoritmo de dos pasos para gestionar el costo computacional: una fase de "calentamiento" (warm-up) que rastrea trayectorias desde un sistema genérico para filtrar trayectorias singulares, seguida de una fase de homotopía de parámetros para rastrear las trayectorias no singulares restantes hacia soluciones físicas.
Integración Simbólica y Numérica: El paquete utiliza Symbolics.jl para la manipulación simbólica con el fin de derivar automáticamente las ecuaciones armónicas y las matrices Jacobianas. Se conecta con DifferentialEquations.jl para simulaciones dependientes del tiempo, permitiendo a los usuarios verificar los resultados de estado estacionario y estudiar barridos adiabáticos de parámetros (histéresis).

Contribuciones Clave

Búsqueda Exhaustiva de Soluciones: La contribución principal es la capacidad de determinar numéricamente todas las soluciones de estado estacionario (reales y complejas) para sistemas con forzamiento armónico, en lugar de una única solución dependiente de las condiciones iniciales.
Marco Unificado: El paquete unifica el álgebra simbólica, el equilibrio armónico y la continuación de homotopía en un único entorno Julia de código abierto, ofreciendo alto rendimiento y flexibilidad.
Análisis de Estabilidad y Fluctuaciones: Proporciona herramientas para clasificar automáticamente las soluciones como reales o complejas, y estables o inestables, basándose en los autovalores de la matriz Jacobiana. También soporta el cálculo de espectros de respuesta lineal para puntos fijos estables.
Visualización: La suite incluye una interfaz gráfica de usuario y rutinas de trazado para visualizar diagramas de fase 1D y 2D, ramas de solución y estructuras de bifurcación, facilitando la exploración de topologías de solución complejas.

Resultos y Ejemplos
El artículo demuestra las capacidades del paquete a través de varios ejemplos:

Oscilador de Duffing Forzado: El paquete identifica con éxito las múltiples soluciones de estado estacionario de un único oscilador de Duffing, incluyendo la respuesta superarmónica ( $3\omega_d$ ) generada por la no linealidad. Reproduce curvas de resonancia estándar y diagramas de fase que muestran regiones con una, dos o tres soluciones.
Osciladores de Duffing Acoplados: Para un sistema de dos osciladores acoplados, el paquete mapea el paisaje de soluciones a través del espacio de parámetros, revelando topologías intrincadas. Contrasta el conjunto de soluciones "completo" hallado vía homotopía con el subconjunto limitado accesible mediante barridos temporales adiabáticos, demostrando la histéresis y los saltos entre ramas.
Oscilador con Forzamiento Paramétrico: El paquete analiza un oscilador de Duffing con forzamiento paramétrico (ecuación de Mathieu), identificando regiones de soluciones estables únicas y dobles, y caracterizando la bifurcación de tipo pitchfork donde el número de estados estables cambia.
Escalabilidad del Rendimiento: Las pruebas de rendimiento en una cadena de $N$ osciladores de Duffing acoplados muestran que, si bien el límite de Bézout (el máximo teórico de soluciones) escala exponencialmente ( $3^{2N}$ ), el número de soluciones reales físicas es significativamente menor. El paquete puede resolver sistemas con $N=5$ (que involucran cientos de trayectorias) en un solo núcleo de CPU en cuestión de minutos, destacando la eficiencia de la estrategia de seguimiento de trayectorias.

Significancia
Los autores posicionan a HarmonicBalance.jl como una herramienta para cerrar la brecha entre el análisis teórico y la observación experimental en la dinámica no lineal. Al permitir la exploración del paisaje completo de soluciones, el paquete permite a los investigadores identificar fenómenos observables experimentalmente, tales como la histéresis y la ruptura de simetría, que a menudo se pasan por alto en las simulaciones estándar de evolución temporal. El software está diseñado para ser accesible a la comunidad académica, ofreciendo un enfoque "ab-initio" donde se pueden definir ecuaciones de movimiento arbitrarias sin las restricciones de una GUI especializada. Los autores señalan que la herramienta es aplicable a una amplia gama de campos, incluyendo la dinámica estructural, la óptica no lineal, la optomecánica y las máquinas de Ising. Aunque la versión actual se centra en la dinámica determinista y los estados estacionarios, el diseño modular está destinado a soportar extensiones futuras, como la detección de ciclos límite y métodos de promediado de orden superior.

Más como este