Heat Equation driven by mixed local-nonlocal operators with non-regular space-dependent coefficients

---

1. Problématique et Contexte

L'article s'intéresse à la résolution du problème de Cauchy pour une équation de la chaleur régie par un opérateur de diffusion mixte local-nonlocal avec des coefficients spatiaux irréguliers (distributions).

L'équation étudiée est de la forme :
$$\begin{cases} \partial_t u(t, x) + L u(t, x) = 0, & (t, x) \in (0, T) \times \mathbb{R}^d, \\ u(0, x) = u_0(x), & x \in \mathbb{R}^d, \end{cases}$$
où l'opérateur $L$ est défini par :
$$L u = -\nabla \cdot (a(x) \nabla u) + (-\Delta)^{s/2} \big( b(x) (-\Delta)^{s/2} u \big) + c(x)u,$$
avec $0 < s < 1$.

Défis principaux :

• Opérateur Mixte : Le terme $-\nabla \cdot (a(x) \nabla u)$ modélise la diffusion locale classique (mouvement brownien), tandis que le terme fractionnaire $(-\Delta)^{s/2} (b(x) (-\Delta)^{s/2} u)$ modélise des interactions à longue portée et des sauts de type Lévy.
• Irrégularité des Coefficients : Les coefficients $a, b, c$ et la donnée initiale $u_0$ sont supposés appartenir à l'espace des distributions $\mathcal{D}'(\mathbb{R}^d)$. Cela inclut des cas où les coefficients peuvent être singuliers (par exemple, des mesures de Dirac).
• Obstruction de Schwartz : Lorsque les coefficients sont des distributions, le produit de distributions n'est pas bien défini en général. Par conséquent, l'équation ne peut pas être interprétée dans le sens classique ou même dans le sens faible standard, car les termes non linéaires (produits de coefficients et de solutions) deviennent illégaux.

2. Méthodologie

Pour surmonter l'obstruction liée au produit de distributions, les auteurs adoptent la théorie des solutions très faibles (very weak solutions), introduite par Garetto et Ruzhansky. La démarche se déroule en trois étapes :

1. Régularisation de Friedrichs :
   Les coefficients singuliers $a, b, c$ et la donnée initiale $u_0$ sont régularisés par convolution avec un noyau de lissage (mollifier) $\psi_{\omega(\varepsilon)}$ dépendant d'un paramètre $\varepsilon \in (0, 1]$. On obtient ainsi des familles de coefficients lisses $(a_\varepsilon, b_\varepsilon, c_\varepsilon)$ et de données initiales $(u_{0,\varepsilon})$.

2. Résolution du Problème Régularisé :
   Pour chaque $\varepsilon > 0$, le problème régularisé est bien posé dans le cadre classique des équations paraboliques. Les auteurs utilisent la théorie variationnelle (formulation faible) pour établir l'existence et l'unicité d'une solution $u_\varepsilon$ dans des espaces d'énergie appropriés (espaces de Sobolev $H^1$ et $H^s$).

3. Analyse des Réseaux de Solutions (Solution Nets) :
   Au lieu de chercher une limite ponctuelle unique (qui pourrait ne pas exister), on étudie le comportement du réseau de solutions $(u_\varepsilon)_\varepsilon$ lorsque $\varepsilon \to 0$.

   • Modération : On définit la notion de "réseau modéré" (moderate net), où la croissance de la norme de la solution (et de ses dérivées) par rapport à $\varepsilon$ est contrôlée par une puissance de $\omega(\varepsilon)^{-N}$.
   • Négligeabilité : On définit la notion de "réseau négligeable", où la norme tend vers zéro plus vite que toute puissance de $\omega(\varepsilon)$.
   • Solution Très Faible : Une solution très faible est définie comme un réseau modéré de solutions régularisées, unique à un terme négligeable près.

3. Résultats Clés

A. Bien-poséness Classique (Cas Régulier)

Dans la Section 3, les auteurs établissent d'abord les résultats de base pour des coefficients bornés et mesurables ($a, b, c \in L^\infty$) satisfaisant des conditions de positivité uniforme ($a \ge a_0 > 0$, etc.).

• Théorème : Existence et unicité d'une solution faible unique dans l'espace $L^\infty(0, T; H^1(\mathbb{R}^d)) \cap C([0, T]; L^2(\mathbb{R}^d))$.
• Estimation d'énergie : Une estimation a priori est dérivée, assurant la stabilité de la solution par rapport aux données initiales.

B. Existence de Solutions Très Faibles (Cas Irrégulier)

Dans la Section 4.2, le cadre est étendu aux coefficients distributionnels.

• Théorème 4.1 : Sous l'hypothèse que les coefficients régularisés sont strictement positifs de manière uniforme et modérés, il existe un réseau de solutions $(u_\varepsilon)_\varepsilon$ qui est modéré dans $L^\infty(0, T; H^1)$. Ce réseau constitue une solution très faible du problème original.

C. Unicité (à négligeabilité près)

Dans la Section 4.3, la question de l'unicité est traitée.

• Théorème 4.2 : Si deux réseaux de régularisation différents (mais équivalents au sens de la négligeabilité) sont utilisés, les réseaux de solutions correspondants diffèrent uniquement par un terme négligeable. Cela garantit que la solution très faible est unique dans l'algèbre des solutions modérées modulo les termes négligeables.

D. Cohérence avec la Théorie Classique

Dans la Section 4.4, les auteurs prouvent que cette nouvelle notion de solution est cohérente avec la théorie classique.

• Théorème 4.3 : Si les coefficients et la donnée initiale sont réguliers (et non distributionnels), la solution très faible obtenue par régularisation converge vers la solution faible classique unique lorsque $\varepsilon \to 0$. La convergence a lieu dans $C([0, T]; L^2) \cap L^2(0, T; H^1)$.

4. Contributions Majeures

1. Extension aux Coefficients Singuliers : C'est la première étude unifiée traitant des équations de la chaleur avec des opérateurs mixtes local-nonlocal et des coefficients spatiaux distributionnels.
2. Cadre Analytique Rigoureux : L'article fournit un cadre mathématique solide (basé sur la régularisation et l'analyse asymptotique des réseaux) pour manipuler des produits de distributions dans des équations paraboliques complexes.
3. Généralité de l'Opérateur : La méthode s'applique à un opérateur général combinant diffusion elliptique, diffusion fractionnaire et potentiel, avec des coefficients variables et irréguliers.
4. Preuve de Cohérence : La démonstration que la solution très faible se réduit à la solution classique lorsque les données deviennent régulières valide la pertinence physique et mathématique de l'approche.

5. Signification et Impact

Ce travail est significatif pour plusieurs domaines :

• Modélisation Biologique et Écologique : Il permet de modéliser des phénomènes de dispersion d'espèces où les coefficients (taux de diffusion, mortalité) peuvent être hétérogènes de manière extrême ou singulière (ex: zones de pollution, barrières physiques).
• Physique des Milieux Complexes : Il offre des outils pour étudier le transport anormal dans des milieux poreux ou turbulents où les propriétés matérielles varient de façon irrégulière.
• Analyse Mathématique : Il ouvre la voie à l'étude d'autres équations aux dérivées partielles (EDP) non linéaires ou hyperboliques avec des coefficients singuliers et des opérateurs non locaux, un domaine où les méthodes classiques échouent souvent.

En résumé, cet article établit une fondation robuste pour l'analyse théorique et numérique de modèles de diffusion hybrides dans des environnements hautement hétérogènes et singuliers.