Numerical study of a Transmission Problem in Elasticity with kind damping

1. Problématique et Contexte

L'article étudie un problème de transmission dans un système élastique unidimensionnel composé de deux matériaux distincts ($(-L, 0)$ et $(0, L)$). Contrairement aux études antérieures sur l'amortissement frictionnel localisé (qui conduit souvent à une décroissance exponentielle) ou l'amortissement viscoélastique de Kelvin-Voigt (qui peut être analytique ou instable selon la localisation), ce travail se concentre sur l'introduction d'un amortissement fractionnaire (dérivée de Caputo d'ordre $\alpha \in (0, 1)$) agissant dans tout le domaine.

Le système modélisé est le suivant :
$$\begin{cases} \rho_1 u_{tt} - k_1 u_{xx} + \partial_t^{\alpha, \eta} u = 0 & \text{dans } (-L, 0) \times (0, \infty) \\ \rho_2 v_{tt} - k_2 v_{xx} + \partial_t^{\alpha, \eta} v = 0 & \text{dans } (0, L) \times (0, \infty) \end{cases}$$
avec des conditions aux limites de Dirichlet aux extrémités et des conditions de transmission continues en $x=0$ (continuité du déplacement et de la contrainte). Le terme $\partial_t^{\alpha, \eta}$ représente une dérivée fractionnaire avec un poids exponentiel $\eta \ge 0$.

L'objectif principal est d'analyser le comportement asymptotique du semi-groupe associé $(S_A(t))_{t \ge 0}$, en particulier le taux de décroissance de l'énergie, et de valider ces résultats théoriques par une étude numérique.

2. Méthodologie

A. Modélisation Augmentée (Approche par champ auxiliaire)

Pour traiter la difficulté mathématique posée par la dérivée fractionnaire non locale en temps, les auteurs utilisent une méthode d'augmentation. En s'appuyant sur les travaux de Choi et MacCamy et sur une représentation spectrale, ils remplacent l'opérateur fractionnaire par un système couplé incluant une équation de diffusion dans un domaine auxiliaire (variable $\xi \in \mathbb{R}$).
Le terme de dissipation est réécrit comme :
$$\partial_t^{\alpha, \eta} u = C \int_{\mathbb{R}} \mu(\xi) \phi(x, t, \xi) d\xi$$
où $\phi$ satisfait une équation de type chaleur avec source dépendante de la vitesse. Cela permet de transformer le problème en un système d'équations aux dérivées partielles classiques (hyperbolique couplé à une équation parabolique), facilitant l'analyse fonctionnelle.

B. Analyse Théorique (Théorie des Semi-groupes)

Le système est formulé comme un problème d'évolution abstrait $\dot{U} = AU$ dans un espace de Hilbert approprié $\mathcal{H}$.

1. Bien-poséité : L'opérateur $A$ est montré comme étant le générateur d'un semi-groupe de contractions $C_0$, garantissant l'existence et l'unicité des solutions.
2. Stabilité Forte : En utilisant le théorème d'Arendt-Batty, les auteurs prouvent que le semi-groupe est fortement stable (convergence vers zéro pour toute donnée initiale dans l'espace d'énergie). Cela repose sur l'absence de valeurs propres imaginaires pures pour l'opérateur $A$.
3. Stabilité Polynomiale : L'analyse du taux de décroissance est effectuée en étudiant la norme de la résolvante $(i\lambda I - A)^{-1}$ sur l'axe imaginaire.
   • Le cas $\eta = 0$ (dérivée fractionnaire pure) conduit à une décroissance polynomiale de l'énergie de l'ordre de $O(t^{-1})$.
   • Le cas $\eta > 0$ (avec poids exponentiel) conduit à une décroissance polynomiale de l'ordre de $O(t^{-\frac{2}{1-\alpha}})$. Pour $\alpha=0.5$, le taux est $O(t^{-4})$.
   • La décroissance exponentielle n'est pas atteinte, ce qui est cohérent avec la nature de l'amortissement fractionnaire localisé.

C. Approximation Numérique

Pour valider les résultats théoriques, les auteurs développent un schéma numérique hybride :

• Discrétisation spatiale : Méthode des Volumes Finis (FVM) pour les variables d'espace $u$ et $v$, assurant la conservation des flux aux interfaces.
• Discrétisation temporelle : Un schéma Newmark (ordre 2) est utilisé pour l'équation d'onde, couplé à une méthode Crank-Nicolson pour l'équation de diffusion auxiliaire $\phi$. Ce couplage garantit la conservation (ou la décroissance non oscillatoire) de l'énergie discrète.
• Approximation de l'intégrale : L'intégrale sur $\xi$ est approchée par une quadrature numérique (somme de Riemann) sur un domaine borné $[-R, R]$.

3. Résultats Clés

1. Décroissance Polynomiale Optimale :

   • Pour $\eta = 0$, l'énergie $E(t)$ décroît comme $t^{-1}$.
   • Pour $\eta > 0$, l'énergie décroît comme $t^{-\frac{2}{1-\alpha}}$.
   • Les auteurs démontrent que ces taux sont optimaux (ils ne peuvent pas être améliorés en exponentiel).

2. Validation Numérique :

   • Les simulations numériques reproduisent fidèlement les phénomènes physiques : réflexion, réfraction des ondes à l'interface de transmission, et amortissement progressif.
   • Les graphiques en échelle log-log montrent une correspondance parfaite entre les courbes de décroissance de l'énergie simulée et les bornes théoriques ($y = C_0 t^{-1}$ pour $\eta=0$ et $y = C_4 t^{-4}$ pour $\eta > 0$ avec $\alpha=0.5$).
   • L'étude confirme la sensibilité du schéma numérique au paramètre $\eta$ : des valeurs trop élevées peuvent déstabiliser le schéma, nécessitant un choix prudent.

4. Contributions et Significativité

• Apport Théorique : L'article établit rigoureusement le taux de stabilité polynomial pour un problème de transmission élastique avec amortissement fractionnaire, comblant un vide entre les études sur l'amortissement frictionnel (exponentiel) et viscoélastique.
• Apport Méthodologique : La combinaison d'une formulation augmentée (pour éviter les dérivées fractionnaires directes) et d'un schéma numérique conservatif (Newmark-Crank-Nicolson) offre une approche robuste pour simuler ce type de systèmes complexes.
• Significativité : Ce travail est crucial pour la modélisation de matériaux réels présentant des comportements viscoélastiques complexes ou des mécanismes d'amortissement non locaux. Il fournit des outils de prédiction précis pour la conception de structures soumises à des vibrations, en montrant que l'amortissement fractionnaire, bien que moins rapide que l'exponentiel, assure une stabilité asymptotique garantie avec un taux de convergence calculable.

En résumé, l'article démontre que l'amortissement fractionnaire dans un système de transmission élastique conduit à une stabilité polynomiale optimale, et fournit une validation numérique convaincante de ces résultats théoriques.