2025-11-16T20:16:12.924612

On linear waves with memory in a Bessel-like medium

Giusti, Colombaro, Mentrelli
We discuss the propagation of harmonic and transient waves for systems governed by a wave equation with memory whose integral kernel involves ratios of modified Bessel functions of the first kind in the Laplace domain. In particular, the investigation of transient waves is carried out by means of a fully numerical approach based on the Talbot method for the numerical inversion of Laplace transforms.
academic

Sur les ondes linéaires avec mémoire dans un milieu de type Bessel

Informations fondamentales

  • ID de l'article: 2510.11493
  • Titre: On linear waves with memory in a Bessel-like medium
  • Auteurs: A. Giusti, I. Colombaro, A. Mentrelli
  • Classification: math-ph math.MP
  • Date de publication: 14 octobre 2025
  • Lien de l'article: https://arxiv.org/abs/2510.11493

Résumé

Cet article traite de la propagation des ondes harmoniques et transitoires dans des systèmes gouvernés par une équation d'onde avec mémoire, dont le noyau intégral implique le rapport des fonctions de Bessel modifiées de première espèce dans le domaine de Laplace. En particulier, l'étude des ondes transitoires est menée par une méthode entièrement numérique basée sur la méthode de Talbot, utilisée pour l'inversion numérique de la transformée de Laplace.

Contexte et motivation de la recherche

Description du problème

Le problème fondamental abordé par cette recherche concerne le comportement de la propagation des ondes dans un milieu de type Bessel présentant des effets de mémoire. Plus précisément, l'objet d'étude est l'équation d'onde avec mémoire de la forme suivante :

∂2Y∂t2−c2[1−Φν(t)∗]∂2Y∂x2=0\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi_\nu(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0

où le noyau de mémoire est représenté dans le domaine de Laplace par : Φ~ν(s)=2(ν+1)sτIν+1(sτ)Iν(sτ)\tilde{\Phi}_\nu(s) = \frac{2(\nu+1)}{\sqrt{s\tau}}\frac{I_{\nu+1}(\sqrt{s\tau})}{I_\nu(\sqrt{s\tau})}

Importance et applications

Ces équations présentent des applications importantes dans plusieurs domaines physiques :

  1. Hémodynamique: description de la propagation des ondes dans les conduits élastiques remplis de liquide
  2. Mécanique viscoélastique: modélisation des effets de mémoire des matériaux
  3. Systèmes électriques: analyse des circuits présentant des caractéristiques de mémoire

Limitations des méthodes existantes

Les recherches antérieures ont principalement utilisé l'algorithme semi-analytique de Buchen-Mainardi, mais cette méthode n'est fiable que près du front d'onde et ne peut pas fournir une image complète de l'évolution spatio-temporelle.

Motivation de la recherche

Cet article vise à surmonter les limitations des méthodes existantes par des approches numériques, en fournissant une analyse complète du comportement de la propagation des ondes dans tout le domaine spatio-temporel.

Contributions principales

  1. Établissement d'un cadre complet d'analyse de dispersion: pour le cas ν=0, analyse détaillée de la relation de dispersion des ondes harmoniques avec atténuation spatiale et périodicité temporelle (SATP)
  2. Dérivation des lois de dispersion exactes et des facteurs d'atténuation: expressions explicites données par les fonctions de Kelvin
  3. Implémentation d'une analyse numérique complète des ondes transitoires: utilisation de la méthode de Talbot pour le calcul numérique de l'inversion de Laplace
  4. Révélation des caractéristiques de dispersion anormale: démonstration que le système présente un phénomène de dispersion anormale où la vitesse de phase est inférieure à la vitesse de groupe

Détail des méthodes

Définition de la tâche

Étude des caractéristiques de propagation des ondes dans un milieu de mémoire de type Bessel, incluant :

  • Analyse des caractéristiques de dispersion et d'atténuation des ondes harmoniques
  • Calcul de la réponse transitoire aux échelons
  • Détermination des vitesses de phase et de groupe

Architecture du modèle

1. Équation d'onde fondamentale

Considération du cas particulier ν=0 : ∂2Y∂t2−c2[1−Φ(t)∗]∂2Y∂x2=0\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0

où : Φ~(s)=2sτI1(sτ)I0(sτ)\tilde{\Phi}(s) = \frac{2}{\sqrt{s\tau}}\frac{I_1(\sqrt{s\tau})}{I_0(\sqrt{s\tau})}

2. Méthode d'analyse de dispersion

Utilisation de la forme standard de solution d'onde harmonique : Y(t,x)=Aei(ωt−kx)Y(t,x) = A e^{i(\omega t - kx)}

Obtention de la relation de dispersion : (iω)2+k2c2[1−Φ^(ω)]=0(i\omega)^2 + k^2 c^2[1 - \hat{\Phi}(\omega)] = 0

3. Analyse des ondes SATP

Pour les ondes avec atténuation spatiale et périodicité temporelle, définition :

  • ω∈R\omega \in \mathbb{R} (fréquence réelle)
  • k=κ−iδattk = \kappa - i\delta_{att} (nombre d'onde complexe)

Résolution du système d'équations non linéaires :

undefined