Inverse problem for a multi-term time-fractional diffusion equation with the Caputo derivatives

1. Problème Étudié

L'article s'intéresse à un problème inverse de source pour une équation de diffusion temporelle fractionnaire à plusieurs termes, définie sur un domaine $\Omega \subset \mathbb{R}^N$. L'équation gouvernante est :

$$\sum_{j=1}^{M} q_j \partial_t^{\rho_j} u(x, t) + A(x, D)u(x, t) = f(x)g(t), \quad x \in \Omega, \, 0 < t \le T$$

où :

• $\partial_t^{\rho_j}$ désigne la dérivée fractionnaire de Caputo d'ordre $0 < \rho_M < \dots < \rho_1 < 1$.
• $A(x, D)$ est un opérateur différentiel elliptique formel auto-adjoint, positif, d'ordre pair $m \ge 2$.
• Le terme source est séparable : $f(x)g(t)$, où $g(t)$ est une fonction connue et $f(x)$ est la composante spatiale inconnue à reconstruire.
• Les conditions aux limites sont homogènes ($B_j u = 0$ sur $\partial\Omega$) et la condition initiale est $u(x, 0) = \phi(x)$.

Condition de surdétermination :
Le problème est clos par une condition supplémentaire donnée à un instant intérieur fixe $t_0 \in (0, T]$ :
$$u(x, t_0) = \Psi(x), \quad x \in \Omega$$

L'objectif est de déterminer l'existence, l'unicité et la régularité de la paire de solutions $(u(x, t), f(x))$.

2. Méthodologie

Les auteurs utilisent une approche spectrale basée sur la méthode de Fourier, combinée à une analyse asymptotique fine des fonctions spéciales.

1. Décomposition Spectrale :
   Le problème est réduit à une suite de problèmes aux valeurs propres pour l'opérateur $A(x, D)$ :
   $$A(x, D)v_k(x) = \lambda_k v_k(x)$$
   La solution $u(x, t)$ et la source $f(x)$ sont développées sur la base orthonormée des fonctions propres $\{v_k(x)\}$ associées aux valeurs propres $\lambda_k > 0$.

2. Réduction à des Équations Différentielles Ordinaires (EDO) :
   En projetant l'équation sur chaque mode spectral, le problème inverse devient une suite de problèmes de Cauchy pour des EDO fractionnaires à plusieurs termes, où les coefficients inconnus sont les coefficients de Fourier de $f(x)$, notés $f_k$.

3. Fonctions de Mittag-Leffler Multinomiales :
   Les solutions des EDO fractionnaires s'expriment naturellement à l'aide de la fonction de Mittag-Leffler multinomiale $E_{(\rho_1, \dots, \rho_M), \beta}(\cdot)$. La solution formelle fait intervenir un dénominateur spectral critique :
   $$b_{k, \rho_1}(t_0) = \int_0^{t_0} g(t_0 - \xi) \xi^{\rho_1-1} E_{\rho', \rho_1}(-\lambda_k \xi^{\rho_1}, \ast) \, d\xi$$
   La reconstruction de $f_k$ dépend de l'inversion de ce terme : $f_k = \frac{\Psi_k - \text{terme initial}}{b_{k, \rho_1}(t_0)}$.

4. Analyse Asymptotique :
   La clé de la démonstration réside dans l'étude du comportement asymptotique de la fonction de Mittag-Leffler multinomiale pour de grandes valeurs de $\lambda_k$. Les auteurs dérivent de nouvelles expansions asymptotiques précises pour établir des bornes inférieures uniformes pour le dénominateur $b_{k, \rho_1}(t_0)$.

3. Contributions Clés et Résultats

A. Analyse Asymptotique Nouvelle (Section 3)

Les auteurs établissent pour la première fois une expansion asymptotique précise de la fonction de Mittag-Leffler multinomiale $E_{\rho', \beta}(z_1, \dots, z_M)$ lorsque $|z_1| \to \infty$ et que les autres variables sont négatives.

• Résultat : Ils montrent que pour $\beta > 2\rho_1$, la fonction se comporte comme $-1/z_1 \Gamma(\beta-\rho_1)^{-1} + O(|z_1|^{-2})$.
• Importance : Cette estimation permet de contrôler le comportement du dénominateur spectral $b_{k, \rho_1}(t_0)$ pour les grands modes $k$.

B. Estimations du Dénominateur Spectral (Section 4)

En utilisant les asymptotiques ci-dessus et les propriétés de la fonction $g(t)$, les auteurs prouvent que le dénominateur $b_{k, \rho_1}(t_0)$ ne s'annule pas (ou s'annule seulement pour un nombre fini d'indices) et qu'il est borné inférieurement par $C/\lambda_k$.

• Cas 1 : Si $g(t)$ ne change pas de signe, alors $b_{k, \rho_1}(t_0) \neq 0$ pour tout $k$.
• Cas 2 : Si $g(t)$ change de signe mais $g(0) \neq 0$, alors $b_{k, \rho_1}(t_0) \neq 0$ pour tout $k$ suffisamment grand. L'ensemble des indices où il s'annule ($K_{0, \rho_1}$) est fini.

C. Existence et Unicité (Sections 5 et 6)

• Existence : Sous des hypothèses de régularité suffisante sur les données ($\phi, \Psi \in D(\hat{A}^{\tau+1})$ avec $\tau > N/2m$) et sur $g(t)$, l'article prouve l'existence d'une solution classique. La convergence absolue et uniforme des séries de Fourier (et de leurs dérivées) est rigoureusement justifiée grâce aux bornes inférieures du dénominateur.
• Unicité :
  • Si l'ensemble $K_{0, \rho_1}$ (indices où le dénominateur s'annule) est vide, la solution est unique.
  • Si $K_{0, \rho_1}$ n'est pas vide, l'unicité échoue. Cependant, les auteurs dérivent des conditions de compatibilité nécessaires et suffisantes ($\Psi_k = \dots$) pour l'existence d'une solution. Dans ce cas dégénéré, la solution n'est pas unique (les coefficients $f_k$ correspondants peuvent être choisis arbitrairement).

D. Cas Particulier du Tore (Section 7)

Les résultats sont appliqués au cas où l'opérateur $A$ est un opérateur elliptique homogène à coefficients constants sur le tore $N$-dimensionnel. Les auteurs identifient explicitement le domaine de l'opérateur fractionnaire avec les espaces de Sobolev $L^2_{\tau m}(\mathbb{T}^N)$, renforçant la clarté structurelle de l'approche spectrale.

4. Signification et Impact

Ce travail apporte plusieurs avancées significatives dans la théorie des problèmes inverses pour les équations fractionnaires :

1. Généralité de l'Opérateur : Contrairement à la plupart des études antérieures limitées à l'opérateur Laplacien ou ses puissances fractionnaires, ce papier traite un opérateur elliptique général auto-adjoint d'ordre $m \ge 2$.
2. Rigueur Mathématique : La dérivation de nouvelles estimations asymptotiques pour les fonctions de Mittag-Leffler multinomiales comble un vide théorique et fournit les outils nécessaires pour prouver la bien-poséité (existence, unicité, stabilité) dans des espaces de régularité élevée.
3. Traitement de l'Ill-Poséness : Bien que les problèmes inverses soient intrinsèquement mal posés, les auteurs montrent que sous des conditions structurelles appropriées (régularité des données et comportement de $g(t)$), des solutions classiques stables peuvent être obtenues.
4. Cadre Analytique Complet : Le papier fournit un cadre analytique complet, allant de la formulation spectrale à la justification de la convergence des séries, en passant par l'analyse des cas dégénérés (non-unicité).

En résumé, ce papier étend considérablement la théorie de bien-posé des problèmes inverses de source pour les équations de diffusion fractionnaires à plusieurs termes, en fournissant des preuves rigoureuses basées sur une analyse asymptotique fine des fonctions spéciales sous-jacentes.