Efficient Krylov solvers for inverse source problem in 2D space-time fractional diffusion equation

Résumé Technique : Solveurs de Krylov efficaces pour les problèmes de sources inverses dans les équations de diffusion fractionnaires espace-temps en 2D

Énoncé du Problème
Ce travail traite du problème de source inverse associé à une équation de diffusion fractionnaire (FDE) espace-temps en deux dimensions présentant des coefficients variables. L'équation directrice implique une dérivée fractionnaire temporel tempérée ($T D^{\xi,\rho}_t$) et des dérivées de Caputo spatiales dans deux dimensions ($CD^{\eta_1}_x, CD^{\eta_2}_y$). L'objectif est de reconstruire le terme de source inconnu $f(x, y)$ étant donné les conditions initiales et des données finales bruitées.

Le problème est intrinsèquement mal posé au sens de Hadamard, où de petites perturbations dans les données d'entrée (niveau de bruit $\epsilon$) peuvent entraîner des erreurs significatives dans la solution. Pour atténuer cela, les auteurs emploient la méthode de la quasi-valeur limite, une technique de régularisation non locale qui transforme le problème mal posé en un problème régularisé bien posé dépendant d'un paramètre.

Méthodologie
L'approche proposée combine la discrétisation numérique, l'analyse spectrale via la théorie des matrices de Toeplitz localement généralisées (GLT) et la conception de préconditionneurs efficaces pour les solveurs itératifs.

1. Discrétisation :

   • Le problème régularisé est discrétisé à l'aide d'un schéma de différences finies. La dérivée fractionnaire temporelle tempérée est approximée par la formule L1, tandis que les dérivées de Caputo spatiales sont également discrétisées via des approximations L1.
   • Ce processus génère un grand système linéaire de type « tout-en-une » (all-at-once) de la forme $A_{\tilde{N}} v_\lambda = z_\lambda$, où $\tilde{N} = (m+1)\nu(n)$ représente les degrés de liberté totaux (combinant les pas de temps $m$ et les points de grille spatiaux $\nu(n) = n_1 n_2$).
   • La matrice de coefficients résultante $A_{\tilde{N}}$ présente une structure complexe à blocs multi-niveaux. En raison du coefficient variable $a(x, y)$, la matrice n'est pas strictement de type Toeplitz mais possède une structure de type bloc de Toeplitz multi-niveaux perturbée par des matrices d'échantillonnage diagonales.

2. Analyse Spectrale (Théorie GLT) :

   • Les auteurs effectuent une analyse spectrale asymptotique des séquences de matrices en utilisant la théorie des matrices de Toeplitz localement généralisées (GLT). Ce cadre permet de caractériser les distributions de valeurs singulières et de valeurs propres des séquences de matrices issues des opérateurs différentiels discrétisés.
   • L'analyse identifie le symbole GLT de la séquence de matrices de coefficients. Le symbole est dérivé comme la somme du symbole associé au bloc de pas de temps (lié à la dérivée tempérée) et du bloc de l'opérateur spatial (impliquant le coefficient variable et les dérivées fractionnaires spatiales).
   • L'analyse confirme que la séquence de matrices est distribuée selon ce symbole, fournissant la base théorique pour le préconditionnement.

3. Stratégie de Préconditionnement :

   • Sur la base de l'analyse GLT, les auteurs proposent un préconditionneur triangulaire par blocs ($P_{\tilde{N}}$). Ce préconditionneur approxime la matrice originale $A_{\tilde{N}}$ en conservant la structure de bloc de sous-diagonale mais en simplifiant les blocs diagonaux.
   • Les résultats théoriques (Théorème 5.2) démontrent que la séquence de matrices préconditionnée $\{P_{\tilde{N}}^{-1} A_{\tilde{N}}\}$ est distribuée comme la fonction identité (symbole = 1) au sens de la GLT. Par conséquent, les valeurs propres du système préconditionné sont faiblement regroupées autour de 1, ce qui est une propriété souhaitable pour la convergence des méthodes de sous-espace de Krylov comme GMRES.
   • La résolution des systèmes linéaires impliquant le préconditionneur est gérée via des techniques standards pour les systèmes de type Toeplitz inférieurs triangulaires (par exemple, en utilisant des approximations circulantes ou des moyennes intégrales du coefficient variable).

Contributions Clés

• Formulation : L'article formule le problème de source inverse pour une FDE espace-temps 2D à coefficients variables et applique la méthode de la quasi-valeur limite pour assurer la stabilité face au bruit.
• Analyse Structurelle : Il fournit une analyse spectrale rigoureuse des matrices de discrétisation « tout-en-une » résultantes en utilisant la théorie GLT, en dérivant explicitement le symbole pour la structure de bloc multi-niveaux.
• Conception de Préconditionneur : Un nouveau préconditionneur triangulaire par blocs est construit, préservant la structure multi-niveaux de la discrétisation. Les auteurs prouvent que ce préconditionneur conduit à un regroupement des valeurs propres autour de l'unité, garantissant théoriquement une convergence efficace.
• Efficacité Computationnelle : La méthode proposée utilise la structure de la matrice pour réduire le coût des produits matrice-vecteur et des applications de préconditionneur à $O(m^2 \nu(n))$ opérations (avec des facteurs logarithmiques pour les résolutions de Toeplitz), ce qui la rend évolutive pour de grandes grilles.

Résultats Numériques
Des expériences numériques ont été menées sur le domaine $\Omega = (0, 1)^2$ avec divers ordres fractionnaires ($\xi, \eta_1, \eta_2$) et tailles de grille.

• Convergence : La méthode GMRES préconditionnée a nettement surpassé la version non préconditionnée. Par exemple, pour des tailles de grille $n_1=n_2=m=25$, le nombre d'itérations est passé de 186 (non préconditionné) à 39 (préconditionné) pour certains paramètres fractionnaires, représentant une réduction d'environ 80 %.
• Robustesse : Le nombre d'itérations pour le système préconditionné est resté presque constant à mesure que la discrétisation temporelle ($m$) augmentait et n'a augmenté que modérément lorsque la discrétisation spatiale ($n_1, n_2$) augmentait. En revanche, le système non préconditionné a montré une croissance rapide du nombre d'itérations avec l'affinement du maillage.
• Coût : Bien que l'application du préconditionneur ajoute une surcharge, le temps CPU total a été considérablement réduit grâce à la diminution drastique du nombre d'itérations GMRES nécessaires pour la convergence.

Signification et Revendications
L'article affirme que le cadre proposé offre une approche efficace et évolutive pour résoudre les problèmes de sources inverses dans les modèles de diffusion fractionnaire. En intégrant la régularisation, la discrétisation préservant la structure et l'analyse spectrale basée sur la GLT, la méthode accélère avec succès le solveur GMRES.

Les auteurs soulignent que le préconditionneur atténue efficacement la détérioration des propriétés spectrales typiquement associée à l'affinement du maillage dans les problèmes fractionnaires. Ils notent que bien que le regroupement soit « faible » (comme c'est souvent le cas pour les problèmes multi-niveaux), le comportement de convergence résultant est robuste à travers différents ordres fractionnaires et niveaux de discrétisation. Ce travail étend les études unidimensionnelles précédentes au cas plus complexe en deux dimensions, démontant que l'exploitation de structures matricielles avancées peut surmonter les défis computationnels posés par les opérateurs fractionnaires non locaux et les problèmes inverses mal posés.

L'article suggère modestement que des travaux futurs pourraient étendre ces méthodes à des géométries plus complexes, des éléments finis d'ordre supérieur, l'analyse isogéométrique et des dimensions spatiales plus élevées, notant que le cadre théorique présenté est adaptable à ces scénarios.