Analytical and Gegenbauer spectral solutions of damped wave equations with singular potentials in ℝn

Résumé Technique : Solutions Spectrales Analytiques et de Gegenbauer d'Équations d'Ondes Amorties avec Potentiels Singuliers dans $\mathbb{R}^n$

Formulation du Problème
Cette étude examine l'existence, l'unicité et l'approximation numérique de solutions d'une équation d'onde amortie multidimensionnelle définie sur un domaine rectangulaire $\Omega = [0, \pi]^n$. L'équation directrice est un opérateur mixte hyperbolo-parabolique incorporant un terme source singulier de la forme $|x|^{-\alpha}$ :
$$\frac{\partial^2 u}{\partial t^2} + \gamma \frac{\partial u}{\partial t} - \Delta u + \frac{u}{|x|^\alpha} = f(x, t)$$
Le problème est formulé sous conditions initiales de Cauchy ($u(x,0)=\phi(x)$, $u_t(x,0)=\psi(x)$) et conditions aux limites de Neumann non homogènes ($\partial u / \partial \nu = g_j$). L'analyse se concentre sur le régime où le paramètre de singularité satisfait $0 < \alpha < n$ et le coefficient d'amortissement $\gamma > 0$. La présence du potentiel singulier introduit des défis analytiques significatifs concernant l'intégrabilité et la régularité, tandis que la nature mixte de l'opérateur nécessite d'aborder simultanément la propagation des ondes et la dissipation diffusive.

Méthodologie
L'article emploie une double approche, à la fois analytique et numérique :

1. Cadre Analytique (Décomposition en fonctions propres) :

   • La solution est décomposée en une partie homogène transitoire ($u_h$) et une partie particulière en régime permanent ($u_p$).
   • La solution particulière est construite à l'aide d'une série de Fourier en temps, réduisant l'EDP dépendante du temps à un système d'équations elliptiques couplées pour chaque harmonique temporelle $k$.
   • En introduisant des fonctions spatiales à valeurs complexes $U_k(x) = A_k(x) - iB_k(x)$, le système est réduit à une équation de type Helmholtz complexe : $(L - \lambda_k)U_k = F_k$, où $L = -\Delta + |x|^{-\alpha}$ et $\lambda_k = k^2 - i\gamma k$.
   • La solution est explicitement représentée sous la forme d'une double série utilisant les fonctions propres $\Phi_m(x) = \prod \cos(m_j x_j)$ du Laplacien de Neumann. Les coefficients sont déterminés par projection sur la base orthogonale, tenant compte à la fois des sources internes et des flux aux limites.

2. Cadre Numérique (Collocation Spectrale de Gegenbauer) :

   • Une méthode spectrale d'ordre élevé est développée en utilisant des polynômes de Gegenbauer avec le paramètre $\lambda=1$ (ce qui réduit aux polynômes de Chebyshev pour une distribution optimale des nœuds).
   • La discrétisation spatiale utilise des nœuds de Chebyshev-Gauss-Lobatto (CGL) sur une grille de produit tensoriel pour gérer la géométrie du domaine.
   • Le potentiel singulier $|x|^{-\alpha}$ est régularisé près de l'origine à l'aide de $V_\epsilon(x) = (|x|^2 + \epsilon^2)^{-\alpha/2}$ afin de prévenir le dépassement numérique (overflow) tout en préservant le comportement asymptotique.
   • L'intégration temporelle est réalisée à l'aide d'un schéma Implicite-Explicite (IMEX) de second ordre, traitant le terme de diffusion rigide de manière implicite et les termes hyperboliques de manière explicite pour garantir la stabilité sans conditions CFL restrictives.

Contributions Clés et Résultats Théoriques
L'article établit les propriétés mathématiques rigoureuses de la solution :

• Convergence des Coefficients : Il est prouvé que les coefficients de l'expansion décroissent asymptotiquement comme $O(|m|^{-n-\alpha})$, garantissant que la série est bien définie malgré la singularité, à condition que $0 < \alpha < n$.
• Stabilité et Non-Résonance : Le coefficient d'amortissement $\gamma > 0$ est montré comme étant critique. Il assure que le dénominateur spectral complexe $|\mu_m - \lambda_k| \geq \gamma k > 0$, empêchant la division par zéro et éliminant les phénomènes de résonance pure.
• Régularité et Espaces de Sobolev :
  • La solution $u$ est montrée comme étant bornée et bien posée dans $L^2(\Omega)$.
  • Le gradient spatial $\nabla u$ reste intégrable au carré (c'est-à-dire $u \in H^1(\Omega)$) uniquement sous la condition plus stricte $0 < \alpha < n-2$. Pour $n-2 \leq \alpha < n$, le gradient présente une singularité à l'origine.
• Stabilité Globale : Une estimation de norme explicite est dérivée, montrant que la norme $L^2$ de la solution décroît exponentiellement selon $e^{-\gamma t/2}$, confirmant la stabilité globale et la dissipation d'énergie.
• Analyse d'Erreur : Des bornes d'erreur théoriques pour la méthode de collocation de Gegenbauer sont établies. L'erreur de discrétisation spatiale converge algébriquement à un taux de $O(N^{-(n-\alpha-1/2)})$, reflétant la perte de régularité induite par la singularité. Pour les singularités faibles, une convergence quasi-spectrale est observée.

Résultats Numériques
Des expériences numériques exhaustives ont été menées en une, deux et trois dimensions.

• Validation : Les simulations valident les taux de convergence théoriques. Pour les singularités faibles ($\alpha = 0.3$), la méthode atteint une convergence d'ordre élevé (taux $> 3.5$). Pour les singularités fortes proches du seuil critique ($\alpha \approx n - 0.5$), le taux de convergence se dégrade gracieusement vers la limite algébrique théorique.
• Comportement de l'Énergie : Les graphiques numériques confirment la décroissance exponentielle de l'énergie au cours du temps, conformément au paramètre d'amortissement $\gamma$.
• Mise à l'échelle Dimensionnelle : Les résultats démontrent que la méthode reste robuste et stable à travers les dimensions $n=1, 2, 3$, avec des coûts de calcul suivant une loi en $O(N^n)$.

Signification
L'article affirme combler une lacune dans la littérature en fournissant une représentation analytique explicite pour les opérateurs mixtes hyperbolo-paraboliques de dimension $n$ sous potentiels singuliers, un problème où les théorèmes d'existence standards échouent souvent à fournir des solutions constructives. Le travail démontre que le coefficient d'amortissement agit comme un stabilisateur fondamental contre la résonance en présence de singularités. De plus, l'approche intégrée analytique-numérique offre un modèle robuste pour caractériser la propagation des ondes et la diffusion dans des milieux non homogènes avec des singularités localisées. Les auteurs positionnent ces découvertes comme directement applicables à des domaines physiques tels que la thermoélasticité, l'interaction fluide-structure et la propagation des ondes électromagnétiques dans des milieux dissipatifs.