Pulse waves in the viscoelastic Kelvin-Voigt model: a revisited approach

Este artículo presenta una nueva solución integral para las ondas de pulso en un medio viscoelástico bajo el modelo de Kelvin-Voigt, la cual evita la transformada inversa de Laplace numérica y ofrece fórmulas más eficientes y expresiones asintóticas para excitaciones tipo delta y escalón.

Lo esencial

¡Claro que sí! Imagina que este artículo científico es como un manual de instrucciones mejorado para entender cómo se comportan ciertos materiales "extraños" cuando los golpeas o los estiras de repente.

Aquí tienes la explicación en español, usando analogías de la vida cotidiana:

1. ¿De qué trata el problema? (El "Sándwich" de Materiales)

Imagina que tienes un material que es una mezcla de dos cosas:

• Un resorte (Elástico): Como una goma elástica. Si lo estiras, vuelve a su sitio inmediatamente.
• Un chupete o miel (Viscoso): Como miel espesa. Si lo mueves, se resiste y se mueve lento.

La mayoría de los materiales reales (como el suelo, el asfalto o incluso el cuerpo humano) son una mezcla de ambos. A esto los científicos le llaman Kelvin-Voigt.

El problema que estudian los autores es: Si golpeas un extremo de una barra hecha de este material (como un terremoto o un golpe), ¿cómo viaja esa onda a través de ella?

2. El viejo método vs. La nueva solución

Antes, para calcular cómo se mueve esa onda, los científicos usaban una herramienta matemática muy potente pero difícil de usar llamada Transformada de Laplace.

• La analogía: Imagina que para saber dónde llegará una pelota lanzada al aire, tienes que resolver una ecuación en un "mundo paralelo" (el plano complejo) y luego traducir la respuesta de vuelta a nuestro mundo. Es como intentar adivinar el sabor de una sopa probando solo el vapor que sale de la olla, pero en un idioma que nadie entiende bien.
• El problema: Para obtener el resultado final, tenías que hacer cálculos numéricos muy complicados en ese "mundo paralelo". Era lento y propenso a errores de computadora.

¿Qué hacen estos autores?
Han encontrado una nueva fórmula (una integral) que es como tener una receta directa.

• En lugar de ir al "mundo paralelo" y traducir, ahora puedes calcular la respuesta del material directamente en nuestro mundo, paso a paso, usando una fórmula que es mucho más sencilla y rápida para las computadoras.

3. Los dos tipos de "Golpes" que probaron

Para verificar su nueva receta, probaron dos escenarios:

• El "Golpe de Martillo" (Pulso Delta): Imagina que das un golpe seco y rápido con un martillo.
  • Resultado: La onda viaja, se desvanece y se dispersa. La nueva fórmula calcula esto muy rápido y con mucha precisión.
• El "Empujón Constante" (Pulso Escalón): Imagina que empujas el material y lo mantienes empujando (como estirar una banda elástica y dejarla ahí).
  • Resultado: El material se deforma y luego se estabiliza. La nueva fórmula también predice esto perfectamente.

4. ¿Por qué es importante? (La analogía del Terremoto)

El papel menciona la sismología (el estudio de los terremotos).

• Cuando ocurre un terremoto, las ondas viajan a través de la corteza terrestre. La Tierra no es un bloque de vidrio perfecto; es un material viscoelástico (tiene partes rígidas y partes que se deforman como plastilina).
• Con la fórmula antigua, calcular cómo se siente el temblor en una ciudad lejana era lento y costoso.
• Con la nueva fórmula de los autores, los científicos pueden simular estos terremotos mucho más rápido y entender mejor cómo se disipa la energía. Es como pasar de calcular manualmente la trayectoria de un cohete a usar un software de simulación moderno.

5. El "Secreto" de los extremos (Asintóticos)

Además de la fórmula general, los autores descubrieron reglas simples para situaciones extremas:

• Si el tiempo es muy corto o la distancia muy larga: La onda se comporta de una manera muy específica (como si fuera una gota de tinta cayendo en agua).
• Si el tiempo es muy largo: El material se asienta y la onda desaparece.

Ellos encontraron fórmulas cortas y fáciles para predecir estos comportamientos sin tener que hacer los cálculos grandes.

En resumen

Este artículo es como si un chef famoso dijera: "Durante años, todos hemos cocinado este plato (la onda en materiales) usando un proceso de 50 pasos muy confuso. Nosotros hemos encontrado un atajo de 5 pasos que da exactamente el mismo sabor, pero es más rápido, más limpio y nos permite saber exactamente qué pasará si cambiamos un ingrediente."

Lo más valioso: Han creado una herramienta matemática nueva que es más eficiente para las computadoras y más fácil de entender para predecir cómo se mueven las ondas en materiales reales, desde el suelo bajo nuestros pies hasta materiales industriales.

Resumen técnico

Aquí tienes un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Resumen Técnico: Ondas de pulso en el modelo viscoelástico de Kelvin-Voigt: un enfoque revisado

Autores: Juan Luis González-Santander, Francesco Mainardi y Andrea Mentrelli.
Fuente: Mathematics (MDPI), Vol. 14 (2025).

1. Planteamiento del Problema

El artículo aborda el problema fundamental de determinar la respuesta mecánica $r(x, t)$ (que puede representar tensión, deformación, desplazamiento o velocidad de partícula) de un medio homogéneo semi-infinito inicialmente en reposo ($x \ge 0$) ante la aplicación de un pulso en el origen ($x=0$).

• Contexto: Este problema es crucial en la reología y la sismología para entender la propagación de ondas en medios lineales dispersivos con disipación.
• Modelo: Se centra específicamente en el modelo de Kelvin-Voigt, un cuerpo viscoelástico de tipo III que no exhibe elasticidad instantánea ($J_g = 0$) ni velocidad de frente de onda finita ($c \to \infty$), comportándose como un fenómeno de difusión.
• Limitación de la literatura previa: Las soluciones existentes suelen requerir la inversión de la transformada de Laplace mediante integración numérica en el plano complejo, un proceso computacionalmente costoso y complejo. Además, las formulaciones integrales anteriores (como las de Hanin o Morrison) presentan dificultades numéricas para ciertos límites o no permiten fácilmente el análisis asintótico.

2. Metodología

Los autores desarrollan una solución analítica en forma integral que evita la inversión directa de la transformada de Laplace. La metodología sigue estos pasos:

1. Formulación en el dominio de Laplace:

   • Se parte de la relación constitutiva de Kelvin-Voigt: $\sigma = G_e [\varepsilon + t_\varepsilon \partial_t \varepsilon]$, donde $t_\varepsilon$ es el tiempo de retardo y $G_e$ el módulo de equilibrio.
   • Se aplican las ecuaciones de movimiento y cinemática en el dominio de Laplace para obtener la solución general $\tilde{r}(x, s) = \tilde{r}_0(s) \exp\left(-\frac{s}{c'}\sqrt{1+t_\varepsilon s} \, x\right)$.

2. Adimensionalización:

   • Se introducen variables adimensionales para simplificar el sistema:
     • $\xi = \frac{x}{c' t_\varepsilon}$ (distancia adimensional).
     • $\tau = \frac{t}{t_\varepsilon}$ (tiempo adimensional).
     • $c' = \sqrt{G_e/\rho}$ (velocidad característica).

3. Inversión mediante Teorema de Convolución y Funciones Especiales:

   • En lugar de invertir directamente, se utiliza el teorema de convolución de Laplace.
   • Se define una función núcleo $f(\xi, \tau)$ relacionada con la transformada inversa de $\exp(-\xi\sqrt{s} + \xi/\sqrt{s})$.
   • Mediante expansiones en series y el uso de funciones de Hermite y la función hipergeométrica regularizada ${}_0F_1$, se deriva una representación integral cerrada para el núcleo $g(\xi, \tau) = \partial_\tau f(\xi, \tau)$.

4. Derivación de Soluciones Específicas:

   • Pulso Delta ($r_0(t) = \delta(t)$): La respuesta se expresa directamente en términos del núcleo $g(\xi, \tau)$.
   • Pulso Escalón ($r_0(t) = \theta(t)$): La respuesta se obtiene integrando el núcleo ponderado por la función de excitación, resultando en una expresión que combina integrales de ${}_0F_1$ y funciones de error complementarias ($\text{erfc}$).

5. Análisis Asintótico:

   • Se derivan fórmulas asintóticas explícitas para los límites de tiempo y distancia: $\tau \to 0, \infty$ y $\xi \to 0, \infty$.
   • Se utilizan expansiones de funciones hipergeométricas generalizadas (${}_0F_2$) y series asintóticas para aproximar el comportamiento en estos regímenes.

3. Contribuciones Clave

• Nueva Forma Integral: Se presenta una solución en forma integral para la respuesta mecánica en modelos de Kelvin-Voigt que es más simple y computacionalmente eficiente que las soluciones inversas de Laplace tradicionales o las formulaciones integrales anteriores (Hanin, Morrison, Dozio).
• Eficiencia Numérica: La nueva formulación evita la integración en el plano complejo. Las comparaciones numéricas realizadas con Mathematica muestran una discrepancia mínima ($< 10^{-9}$) entre la solución integral propuesta y la solución exacta por inversión de Laplace, validando su precisión.
• Análisis Asintótico Completo: Se obtienen fórmulas asintóticas simples y explícitas para ambos tipos de excitación (delta y escalón) en los límites de corto/largo tiempo y corta/larga distancia, algo que resultaba difícil con las formulaciones previas.
• Corrección de Literatura: Se señala que la solución integral propuesta por Dozio (1990) parece ser incorrecta según los experimentos numéricos de los autores, mientras que la de Morrison es difícil de calcular numéricamente para tiempos grandes ($\tau \gg 1$).

4. Resultados Principales

• Solución General: La respuesta $r(\xi, \tau)$ se expresa como una convolución de la excitación $r_0$ con un núcleo integral que involucra ${}_0F_1$.
• Comportamiento Asintótico:
  • Para $\tau \to 0$ o $\xi \to \infty$ (Pulso Delta): La respuesta decae exponencialmente con un factor de difusión gaussiana: $r \approx \frac{\xi}{2\sqrt{\pi}\tau^{3/2}} \exp\left(-\frac{\xi^2}{4\tau} - \tau\right)$.
  • Para $\tau \to \infty$ (Pulso Escalón): La respuesta tiende a 1, con correcciones que involucran funciones hipergeométricas ${}_0F_2$.
  • Para $\xi \to 0$ (Pulso Escalón): Se obtiene una expansión en serie que muestra cómo la respuesta se aproxima a la unidad con correcciones de orden $\xi$ y $\xi^2$.
• Validación: Los gráficos presentados en el artículo confirman la equivalencia numérica entre la nueva fórmula integral, la inversión de Laplace y las aproximaciones asintóticas en todo el dominio de interés.

5. Significado e Impacto

Este trabajo es significativo por varias razones:

1. Aplicabilidad en Sismología: Ofrece herramientas más eficientes para modelar la propagación de ondas sísmicas en medios viscoelásticos, donde la disipación es un factor crítico.
2. Herramienta Computacional: Proporciona una alternativa robusta para simulaciones numéricas, eliminando la necesidad de costosas integrales de contorno complejas.
3. Avance Teórico: Cierra la brecha entre la solución teórica exacta y su implementación práctica, facilitando el análisis de comportamiento asintótico que es esencial para entender la física del problema en regímenes extremos (tiempos muy cortos o muy largos).
4. Recursos Abiertos: Los autores han hecho disponible un cuaderno de Mathematica que permite a la comunidad científica verificar y utilizar sus resultados, fomentando la reproducibilidad.

En conclusión, el artículo revisa y mejora sustancialmente el tratamiento matemático de las ondas transitorias en medios de Kelvin-Voigt, ofreciendo soluciones integrales superiores tanto en eficiencia computacional como en utilidad para el análisis analítico.