Adaptive homotopy continuation for robust dispersion curve computation in viscoelastic waveguides: guaranteed branch identity continuity

Este artículo introduce un marco de continuación de homotopía de material adaptativo que garantiza la continuidad de la identidad de las ramas y permite el cálculo robusto y automatizado de las curvas de dispersión en guías de onda viscoelásticas de sección transversal arbitraria, mapeando el problema no hermítico con pérdidas a uno auxiliar sin pérdidas mientras gestiona eficazmente los puntos excepcionales y los desafíos de seguimiento de modos.

Lo esencial

La Gran Imagen: Mapeando una Montaña Nebosa

Imagina que estás intentando dibujar un mapa de una cordillera. Esta montaña representa cómo viajan las ondas sonoras a través de un material (como un ala de fibra de carbono en un avión).

• La Montaña "Elástica" (Día Claro): En un material perfecto y sin pérdidas (como un resorte rígido), la montaña está clara. Puedes ver cada pico y valle perfectamente. Los caminos (ondas) son distintos y fáciles de seguir.
• La Montaña "Viscoelástica" (Día Neboso): Los materiales del mundo real (como la fibra de carbono con pegamento) absorben energía. Esto es como si una niebla espesa se deslizara. Los caminos se vuelven borrosos, se retuercen entre sí y, a veces, dos caminos parecen fusionarse en uno antes de separarse de nuevo. Esto se llama "desviación de modos".

El Problema:
Los métodos existentes para dibujar este mapa intentan navegar directamente por la montaña nebulosa. Comienzan en la niebla, adivinan dónde está un camino e intentan seguirlo. Pero como la niebla es tan espesa y los caminos se retuercen tan salvajemente, los cartógrafos a menudo se pierden. Podrían cambiar accidentalmente de seguir el Camino A al Camino B, o podrían perder un camino por completo. Esto resulta en un mapa roto e inexacto.

La Solución: El Ascensor de "Homotopía"

Los autores de este artículo proponen una nueva estrategia ingeniosa. En lugar de intentar navegar directamente por la niebla, construyen un ascensor que conecta el día claro con el día nebuloso.

1. Paso 1: Mapear el Día Claro Primero.
   Comienzan en la parte inferior del ascensor donde el aire está perfectamente claro (el estado "elástico"). Aquí, los caminos son rectos y distintos. Dibujan todo el mapa perfectamente, etiquetando cada camino individual (Modo 1, Modo 2, etc.) con un 100% de certeza.

2. Paso 2: El Viaje Lento hacia Arriba.
   Luego presionan lentamente el botón para subir en el ascensor. A medida que ascienden, la niebla (amortiguación/pérdida de material) se espesa gradualmente.

   • El Truco Mágico: Como se mueven lentamente y de forma continua, pueden observar los caminos que ya etiquetaron. Incluso a medida que la niebla se espesa, pueden ver que el "Camino A" sigue siendo el "Camino A", solo ligeramente distorsionado. No tienen que adivinar; simplemente siguen el rastro que ya hicieron.

3. Paso 3: Llegando a la Niebla.
   Para cuando llegan a la cima (el estado "viscoelástico"), tienen un mapa completo y preciso de la montaña nebulosa. Como nunca perdieron el rastro de los caminos durante el viaje, las etiquetas que dieron a los caminos en la parte inferior siguen siendo correctas en la cima.

Conceptos Clave Explicados Simplemente

1. La "Identidad de la Rama" (La Etiqueta de Nombre)
En el mundo nebuloso, los caminos pueden acercarse mucho y parecer que intercambian lugares.

• La Vieja Forma: Si miras la montaña nebulosa, podrías pensar: "Oh, ese camino parece estar cruzando al otro", y accidentalmente intercambias sus nombres.
• La Nueva Forma: Como los autores rastrearon los caminos desde el día claro, saben con certeza que el "Camino A" nunca realmente intercambió lugares con el "Camino B". Mantuvieron las etiquetas de nombre en los caminos correctos todo el tiempo.

2. Los "Puntos Excepcionales" (Los Remolinos Nebulosos)
A veces, la niebla se vuelve tan espesa que dos caminos se fusionan realmente en un único remolino antes de separarse de nuevo. Esto se llama un "Punto Excepcional".

• Tipo I (Zona Segura): En la mayoría de los materiales comunes, estos remolinos ocurren "al margen" en el mundo matemático. Los caminos en nuestro mapa simplemente se acercan, se retuercen y se pasan sin confundirse. El nuevo método maneja esto perfectamente.
• Tipo II (Zona de Peligro): Si el material es extremadamente disipativo (niebla muy espesa), el remolino podría moverse justo sobre el camino. En este caso raro, los caminos sí intercambian identidades. El artículo admite que si esto sucede, las etiquetas automáticas podrían mezclarse. Sin embargo, el método es lo suficientemente inteligente como para dar la alarma: "Oye, los caminos están haciendo algo extraño aquí; es posible que necesites intercambiar las etiquetas manualmente".

3. Por Qué Esto Es Mejor

• Métodos Antiguos: Como intentar caminar a través de un bosque denso en la oscuridad, tropezando con raíces y adivinando hacia dónde está el Norte. A menudo terminas perdido o caminando en círculos.
• Este Método: Como caminar a través del bosque a plena luz del día, memorizando la ruta y luego caminando de nuevo en la oscuridad sabiendo exactamente dónde está cada árbol.

Lo Que el Artículo Demostró Realmente

Los autores probaron este método de "ascensor" en varias formas diferentes de materiales (placas planas, apilados asimétricos y una barra en forma de L).

• El Resultado: Su método produjo mapas perfectos y continuos para casi todos los casos de prueba, incluso cuando el material era bastante "disipativo" (absorbía sonido).
• La Comparación: Lo compararon con el mejor software actual (llamado "Calculador de Dispersión"). El software antiguo a menudo se perdía, perdía caminos o dibujaba líneas irregulares y rotas en las áreas complicadas. El nuevo método dibujó líneas suaves y correctas cada vez.
• El Límite: El método funciona mejor cuando la "niebla" no es demasiado espesa. Si el material es extremadamente disipativo (un caso muy raro), las etiquetas automáticas podrían confundirse, pero las matemáticas detrás de escena siguen siendo precisas; solo necesitas corregir las etiquetas al final.

Resumen

Este artículo presenta una forma inteligente de calcular cómo se mueven las ondas sonoras a través de materiales complejos que absorben energía. En lugar de luchar para resolver el problema desordenado directamente, resuelve primero la versión limpia y luego transforma lentamente la solución en la versión desordenada. Esto garantiza que la "identidad" de las ondas nunca se pierda, resultando en un mapa mucho más confiable y preciso para los ingenieros.

Resumen técnico

Resumen Técnico: Continuación de Homotopía Adaptativa para el Cálculo Robusto de Curvas de Dispersión

Enunciado del Problema
El cálculo preciso de las curvas de dispersión en guías de onda viscoelásticas de sección transversal arbitraria presenta dos desafíos numéricos entrelazados: resolver problemas de autovalores no hermíticos y rastrear de manera fiable las ramas modales a través de la frecuencia. En materiales con pérdidas, como los polímeros reforzados con fibra de carbono (CFRP), la amortiguación material introduce aritmética compleja que destruye la estructura hermítica de las matrices del sistema. Esto conduce a un problema de autovalores polinómico en el dominio del número de onda complejo, donde los solucionadores iterativos estándar luchan con la selección del parámetro de desplazamiento, y los métodos de búsqueda de raíces se vuelven engorrosos en regiones de modos densamente poblados.

Crucialmente, la presencia de amortiguación material complica el rastreo de modos. En regiones de desviación modal (mode veering), los autovectores intercambian características rápidamente en intervalos de frecuencia estrechos. En sistemas no hermíticos, estas regiones pueden involucrar puntos excepcionales (EP) donde los autovalores y autovectores se fusionan, provocando que el jacobiano del sistema extendido se vuelva singular. Los métodos existentes, que típicamente fijan el estado material en la configuración viscoelástica objetivo e intentan una continuación directa en el dominio de la frecuencia, a menudo fallan en mantener la identidad correcta del modo. Esto resulta en "saltos de modo", cruces artificiales o la omisión completa de modos, particularmente en laminados asimétricos o escenarios de alta amortiguación.

Metodología
El artículo propone un marco de continuación de homotopía adaptativa (HC) que desacopla el difícil problema de rastreo de modos no hermíticos del proceso de solución de autovalores. La innovación central es una homotopía material parametrizada por un factor de atenuación escalar $s \in [0, 1]$, que mapea continuamente el problema viscoelástico objetivo (con pérdidas, no hermítico) en $s=1$ hacia un problema elástico auxiliar (sin pérdidas, hermítico) en $s=0$.

La metodología procede en dos etapas distintas:

1. Etapa 1: Solución Hermítica y Rastreo de Modos ($s=0$):

   • El marco resuelve primero el problema de dispersión para el límite elástico sin pérdidas. En este régimen, el sistema es hermítico, asegurando autovalores reales y autovectores mutuamente ortogonales.
   • El rastreo de modos se realiza utilizando el Criterio de Garantía Modal (MAC) combinado con una estrategia de refinamiento adaptativo del número de onda. Esta estrategia inserta puntos de muestreo en regiones de variación rápida de los autovectores (regiones de desviación) para asegurar una conectividad modal inequívoca.
   • Una estrategia de mapeo disperso selecciona un subconjunto de puntos de solución "clave" del conjunto de datos hermítico denso. Estos puntos se eligen para preservar la continuidad modal (a través de umbrales MAC) y la precisión geométrica (a través de umbrales de error de interpolación), reduciendo significativamente la carga computacional para la etapa subsiguiente.

2. Etapa 2: Rastreo Adaptativo de la Ruta de Homotopía ($s=0 \to s=1$):

   • Las soluciones clave identificadas se propagan desde el estado elástico hasta el estado viscoelástico objetivo utilizando un algoritmo de homotopía predictor-correktor con parametrización de longitud de arco.
   • La ruta de homotopía se define en el espacio de parámetros materiales ($s$) a frecuencias reales fijas, en lugar de en el dominio de la frecuencia. Esto desacopla el proceso de rastreo de las interacciones modales complejas (desviación) que ocurren en el dominio de la frecuencia en el estado objetivo.
   • Un control adaptativo del tamaño de paso ajusta el paso de homotopía $\Delta s$ basándose en la brecha de autovalores local y la curvatura de la ruta, asegurando robustez cerca de regiones de fuerte interacción modal.
   • El algoritmo emplea un paso de refinamiento hacia atrás para asegurar que la solución en $s=1$ se obtenga con precisión de Newton completa.

Fundamentos Teóricos
El marco se basa en la teoría de perturbaciones analítica, que garantiza que los pares de autovalores varían analíticamente con respecto al parámetro de homotopía $s$ siempre que la ruta no interseque un punto excepcional en el plano complejo $s$. Esto asegura una continuidad de identidad de rama—una correspondencia uno a uno entre las soluciones en $s=0$ y $s=1$.

El artículo distingue entre dos regímenes topológicos basados en la ubicación de los EP en el plano de frecuencia complejo con respecto al eje real:

• Topología Tipo I: Los dos EP que gobiernan una interacción modal se encuentran en lados opuestos del eje de frecuencia real. En este régimen, las etiquetas de modo físicas establecidas en $s=0$ se heredan automáticamente en $s=1$ sin post-procesamiento. Las curvas de dispersión resultantes exhiben un comportamiento característico "Tipo I": desviación de la parte real acompañada de cruce de la parte imaginaria.
• Topología Tipo II: Los EP se encuentran en el mismo lado del eje real. Aquí, la herencia automática de etiquetas falla, y las etiquetas de modo deben intercambiarse en el punto de cruce para mantener la continuidad física.

Contribuciones Clave

1. Síntesis Teórica: Los autores sintetizan la teoría de perturbaciones analítica, la física de los EP y la teoría de cruces no hermíticos para establecer que la identidad de la rama se preserva a lo largo de la ruta de homotopía material siempre que no se cruce ningún EP. Aclaran que para topologías Tipo I, las etiquetas de modo físicas se preservan automáticamente, eliminando la necesidad de post-procesamiento heurístico.
2. Marco Metodológico: Se desarrolla un marco computacional completo que realiza un rastreo de modos preciso en el régimen hermítico y propaga estas identidades al estado viscoelástico. El método es inherentemente paralelizable y se aplica tanto a modelos de amortiguación independientes de la frecuencia como dependientes de la frecuencia.
3. Herramientas de Diagnóstico: Para casos donde se viola la suposición Tipo I (régimen Tipo II), el artículo identifica dos firmas diagnósticas post-hoc motivadas físicamente:
   • Un cruce extremadamente agudo de la parte imaginaria.
   • Una discrepancia discernible entre la velocidad de grupo espectral ($v_g$) y la velocidad de flujo de energía ($v_e$).
     Estos indicadores alertan a los usuarios sobre posibles desajustes de etiquetas, guiando los intercambios de etiquetas manuales necesarios.

Resultados Numéricos y Validación
El marco se validó contra el ampliamente utilizado paquete de herramientas Dispersion Calculator (DC) en cinco ejemplos numéricos:

• Laminados Simétricos (Sym1, Sym2): El método resolvió correctamente los cruces protegidos por simetría y las regiones de desviación donde DC identificó erróneamente la desviación como cruces o produjo separaciones artificiales en el número de onda imaginario.
• Laminados Asimétricos (UnSym1, UnSym2): En estos casos, la desviación modal generalizada hizo que DC perdiera modos por completo o sufriera saltos de modo espurios. El marco HC capturó con éxito el espectro completo con la conectividad modal correcta, incluso en presencia de desviación densa.
• Sección Transversal Arbitraria (Barra en L): El método demostró aplicabilidad a geometrías 2D generales, capturando correctamente el comportamiento de desviación y el sutil desplazamiento entre $v_g$ y $v_e$ causado por la amortiguación no hermítica.
• Alta Amortiguación (UnSym3, $\eta=0.05$): En un régimen que se acerca a la topología Tipo II, el método HC continuó produciendo números de onda numéricamente precisos donde DC falló. Las firmas diagnósticas (cruce imaginario agudo y discrepancia $v_g$-$v_e$) señalaron con éxito la necesidad de intercambio de etiquetas.

Significado y Afirmaciones
El artículo afirma que el marco propuesto de continuación de homotopía adaptativa ofrece una alternativa robusta y teóricamente fundamentada a los métodos existentes de cálculo de dispersión. Al desacoplar el rastreo de modos del proceso de solución no hermítico, el método supera la fragilidad de la continuación directa en el dominio de la frecuencia en presencia de puntos excepcionales y desviación modal.

Los autores enfatizan que el marco garantiza la continuidad de identidad de rama a lo largo de la ruta de homotopía, independientemente de la topología EP del sistema objetivo. Si bien la herencia de etiquetas físicas es automática solo para topologías Tipo I, el método asegura que incluso en regímenes Tipo II, las soluciones numéricas permanezcan precisas, con herramientas de diagnóstico proporcionadas para guiar las correcciones de etiquetas necesarias. Este enfoque transforma al solucionador de dispersión de una "caja negra" propensa a errores irreversibles en una herramienta transparente capaz de manejar interacciones modales complejas en guías de onda viscoelásticas de sección transversal arbitraria. Se demuestra que el método es computacionalmente eficiente, particularmente para laminados simétricos, y escalable a grandes modelos de elementos finitos mediante procesamiento paralelo y mapeo disperso adaptativo.