Blow-up of solutions to the Euler-Poisson-Darbox equation with critical power nonlinearity

1. Problème Étudié

L'article se concentre sur l'analyse du comportement des solutions de l'équation d'Euler-Poisson-Darboux (EPD) semi-linéaire avec un terme d'amortissement singulier et une source non linéaire critique. Le problème de Cauchy singulier considéré dans $\mathbb{R}^n$ ($n \ge 1$) est donné par :

$$\begin{cases} u_{tt} - \Delta u + \frac{\mu}{t} u_t = t^\alpha |u|^p, & t > 0, \, x \in \mathbb{R}^n, \\ u(0, x) = u_0(x), \\ u_t(0, x) = 0, \end{cases}$$

où :

• $\mu > 0$ est un paramètre d'amortissement (le coefficient $\frac{\mu}{t}$ devient singulier à $t=0$).
• $\alpha \ge 0$ et $p > 1$ sont des exposants non linéaires.
• Les données initiales $u_0$ sont positives, appartiennent à $H^1(\mathbb{R}^n)$ et ont un support compact.

L'objectif principal est d'établir un résultat d'explosion en temps fini (blow-up) et une estimation supérieure de la durée de vie (lifespan) pour le cas critique, où l'exposant $p$ correspond à la racine positive $p_S(n+\mu, \alpha)$ de l'équation quadratique caractéristique :
$$\gamma(n, \mu, \alpha, p) := 2 + (n + \mu + 1 + 2\alpha)p - (n + \mu - 1)p^2 = 0.$$

Ce cas est considéré comme critique car il représente la frontière entre l'existence globale et l'explosion en temps fini. Le cas sous-critique ($p < p_S$) avait déjà été traité dans un travail précédent des mêmes auteurs.

2. Méthodologie

La démonstration repose sur une méthode variationnelle améliorée, utilisant des fonctions tests spécifiques et une argumentation itérative. Les étapes clés sont :

A. Construction d'une Fonction Test Améliorée

Pour surmonter la difficulté posée par la singularité de l'amortissement à $t=0$ et obtenir une croissance logarithmique supplémentaire nécessaire au cas critique, les auteurs construisent une fonction test $b_q(t, x)$ sophistiquée.

• Fonction de base : Ils utilisent la fonction de Bessel modifiée $K_{\frac{\mu-1}{2}}(t)$ pour définir une solution $h(t)$ de l'équation ODE associée à l'opérateur d'amortissement.
• Intégration sur l'échelle : La fonction test est définie comme une intégrale paramétrée sur une échelle $\lambda$ :
  $$b_q(t, x) = \int_0^1 h(\lambda t) \phi(\lambda x) \lambda^{q-1} d\lambda,$$
  où $\phi(x)$ est une fonction radiale liée aux harmoniques sphériques.
• Avantage : Cette construction permet d'obtenir un comportement asymptotique précis ($b_q \sim t^{-q}$) et d'éliminer les singularités à l'origine grâce aux propriétés de récurrence des fonctions de Bessel.

B. Fonctionnelle d'Explosion et Estimations

Les auteurs définissent une fonctionnelle $Y(M)$ intégrant le terme non linéaire sur l'espace-temps, pondérée par la fonction test $b_q$ et une fonction de coupure $\eta$.

• Majoration (Upper Bound) : En utilisant l'inégalité de Hölder et les propriétés asymptotiques de $b_q$ et de sa dérivée temporelle, ils montrent que $Y(t)$ est contrôlé par une intégrale de la non-linéarité multipliée par un facteur logarithmique $(\ln t)^{1/p'}$.
• Minoration (Lower Bound) : Grâce à la construction de la fonction test et à l'hypothèse de données initiales positives, ils établissent une minoration de l'intégrale du terme non linéaire. Dans le cas critique, cette minoration présente une croissance logarithmique supplémentaire : $Y(t) \gtrsim \varepsilon^p \ln t$.

C. Argumentation Différentielle

En combinant la majoration et la minoration, ils obtiennent une inégalité différentielle de la forme :
$$t Y'(t) \ge C (\ln t)^{1-p} Y(t)^p.$$
L'intégration de cette inégalité conduit à une contradiction si le temps d'existence est infini, prouvant ainsi l'explosion en temps fini.

3. Résultats Principaux

Le théorème principal (Théorème 1.2) établit les conditions suivantes pour l'explosion en temps fini :

• Conditions sur les paramètres :
  • Si $n \ge 3$, alors $\alpha \ge 0$.
  • Si $n = 2$, alors $\alpha > 0$.
  • $0 < \mu \le \mu^*(n, \alpha)$, où $\mu^*$ est une borne supérieure explicite dépendant de $n$ et $\alpha$.
  • $p = p_S(n+\mu, \alpha)$ (le cas critique).
• Estimation de la durée de vie (Lifespan) :
  Pour des données initiales petites de taille $\varepsilon$, la durée de vie maximale $T(\varepsilon)$ de la solution satisfait :
  $$T(\varepsilon) \le \exp\left( C \varepsilon^{-p(p-1)} \right),$$
  où $C$ est une constante positive indépendante de $\varepsilon$.

Ce résultat confirme que dans le régime critique, la durée de vie croît de manière exponentielle inversement proportionnelle à une puissance de la taille des données initiales, ce qui est typique des problèmes critiques (par opposition au cas sous-critique où la durée de vie est souvent de type puissance).

4. Contributions Clés

1. Résolution partielle d'un problème ouvert : L'article répond partiellement au "Problème ouvert 1" posé par D'Abbicco (2021) concernant l'absence de résultats complets d'explosion pour l'équation EPD semi-linéaire.
2. Traitement du cas critique : C'est la première fois qu'une estimation précise de la durée de vie est obtenue pour le cas critique de l'équation EPD avec amortissement singulier et source $t^\alpha |u|^p$.
3. Nouvelle fonction test : L'introduction d'une fonction test construite via une intégration sur un paramètre d'échelle de fonctions de Bessel modifiées permet de gérer la singularité à $t=0$ et d'extraire le terme logarithmique crucial pour le cas critique.
4. Gestion de la singularité : La méthode proposée contourne efficacement la difficulté technique liée à l'intégration par parties avec le coefficient $\frac{\mu}{t}$ à l'origine, un obstacle majeur dans les travaux antérieurs sur les problèmes réguliers.

5. Signification et Impact

Ce travail est significatif pour la théorie des équations aux dérivées partielles hyperboliques non linéaires :

• Il complète la classification des comportements de solutions (existence globale vs explosion) pour l'équation EPD, un modèle important en physique (théorie des surfaces, propagation du son, champs gravitationnels).
• Il démontre que la présence d'un amortissement singulier $\frac{\mu}{t}$ modifie la nature critique de l'équation (définie par $p_S(n+\mu, \alpha)$) par rapport à l'équation des ondes classique.
• La technique de construction de fonction test par intégration paramétrique pourrait être appliquée à d'autres modèles d'ondes avec coefficients singuliers ou non linéarités critiques.

En résumé, l'article fournit une preuve rigoureuse de l'explosion en temps fini pour le cas critique de l'équation EPD, en établissant une borne supérieure précise pour la durée de vie, comblant ainsi une lacune importante dans la littérature mathématique actuelle sur ce sujet.