Interface problems of mixed spatial order

Résumé technique : Problèmes d'interface à ordres mixtes via la méthode de la transformée unifiée

Énoncé du problème
Ce travail traite du problème de la valeur initiale aux limites (IBVP) pour les équations aux dérivées partielles (EDP) d'évolution linéaire définies sur un domaine spatial divisé en deux demi-droites, $(-\infty, 0)$ et $(0, \infty)$, séparées par une interface en $x=0$. Contrairement aux études précédentes où les équations directrices de chaque côté de l'interface étaient identiques ou ne différaient que par des coefficients scalaires, cet article étudie des problèmes d'interface à « ordres mixtes ». Plus précisément, il considère l'interaction entre des équations d'évolution linéaire de ordres spatiaux et de comportements physiques fondamentalement différents (dispersif vs dissipatif).

Les exemples prototypiques analysés impliquent :

1. Demi-droite gauche ($x < 0$) : l'équation de Schrödinger linéaire (lS), l'équation de la chaleur (h), ou l'équation de Schrödinger biharmonique (biS).
2. Demi-droite droite ($x > 0$) : l'équation d'Airy ou l'équation de Korteweg-de Vries linéarisée (lKdV).

Ces systèmes sont couplés via des conditions d'interface en $x=0$. L'article considère deux ensembles de conditions :

• Condition (IfC.1) : continuité de la solution et de sa première dérivée spatiale ($C^1$).
• Condition (IfC.2) : continuité de la solution et de ses deux premières dérivées spatiales ($C^2$), requise pour l'équation de Schrödinger biharmonique d'ordre quatre.

Les valeurs d'interface (données aux limites à $x=0$) ne sont pas prescrites a priori ; ce sont des fonctions inconnues déterminées par le couplage des deux EDP distinctes.

Méthodologie
Les auteurs emploient la Méthode de la Transformée Unifiée (UTM), également connue sous le nom de méthode de Fokas, qui constitue une généralisation de la transformée de Fourier pour les IBVP sur des domaines avec des frontières. La procédure de résolution comprend les étapes suivantes :

1. Relations globales : Pour chaque problème de demi-droite, les auteurs dérivent des « relations globales » (identités spectrales) reliant les données initiales, la solution et les valeurs de bord inconnues dans le plan spectral complexe.
2. Changement de variables : Pour traiter les ordres mixtes, les auteurs introduisent des fonctions spécifiques de type racine ($\sigma, \rho, \tau, \gamma$) et un biholomorphisme $\nu$ (pour le cas lKdV) qui projettent le paramètre spectral $\lambda$ vers le plan complexe $k$ de manière à aligner les relations de dispersion des différentes équations (par exemple, en reliant $k^2, k^3$ et $k^4$).
3. Construction du système linéaire : En appliquant ces changements de variables aux relations globales, les auteurs construisent un système linéaire (typiquement $3\times3$ ou $4\times4$) pour les transformées des valeurs d'interface inconnues.
4. Représentation de la solution : La résolution de ce système linéaire permet d'éliminer les valeurs d'interface inconnues. Les auteurs réinjectent ensuite ces expressions dans les représentations intégrales des solutions pour chaque demi-droite.
5. Déformation de contour et asymptotique : Une contribution technique critique est la justification rigoureuse de la déformation des contours d'intégration (par exemple, de l'axe réel vers des contours comme $\tilde{R}$ dans le demi-plan inférieur) afin de contourner les singularités (spécifiquement à $\lambda=0$) et d'éliminer les intégrales impliquant les termes dépendants du temps ($u(x,t), v(x,t)$). Cela est réalisé grâce à des arguments sur mesure utilisant le comportement asymptotique spécifique des fonctions racines et des noyaux exponentiels, plutôt que par le lemme de Jordan standard, en raison du décalage des termes de l'ordre principal.

Contributions clés et résultats
L'article fournit des formules de solutions explicites pour quatre problèmes d'interface distincts, présentés dans les Propositions 1 à 4 :

• Proposition 1 : Résout l'interface entre l'équation de Schrödinger linéaire et l'équation d'Airy sous des conditions d'interface $C^1$.
• Proposition 2 : Résout l'interface entre l'équation de Schrödinger linéaire et l'équation de KdV linéarisée (incluant un terme d'ordre inférieur $av_x$) sous des conditions $C^1$. Cela nécessite l'introduction du biholomorphisme $\nu$ pour gérer la relation de dispersion non homogène.
• Proposition 3 : Résout l'interface entre l'équation de la chaleur (dissipative) et l'équation d'Airy (dispersive) sous des conditions $C^1$.
• Proposition 4 : Résout l'interface entre l'équation de Schrödinger biharmonique (dispersive d'ordre quatre) et l'équation d'Airy sous des conditions $C^2$. Ce cas implique un système linéaire $4\times4$ en raison de l'ordre supérieur de l'EDP.

Les formules explicites sont exprimées sous forme d'intégrales de contour impliquant les transformées des données initiales et des fonctions de noyau spécifiques dérivées des applications de type racine.

Évaluation numérique
Pour démontrer l'utilité des formules dérivées, les auteurs effectuent des évaluations numériques en utilisant la quadrature de Filon. Cette méthode est choisie pour gérer la nature hautement oscillatoire des intégrandes. L'article présente des graphiques des solutions pour diverses configurations de données initiales :

• Schrödinger-Airy : Montre la transition des ondes dispersives du régime unidirectionnel d'Airy vers le régime bidirectionnel de Schrödinger, et vice versa, mettant en évidence les motifs de type « roulette » et les effets de traînée causés par l'interface.
• Biharmonique Schrödinger-Airy : Illustre des interactions d'ondes similaires mais avec une dynamique régie par une relation de dispersion quartique, menant à l'émergence plus rapide de motifs complexes.
• Chaleur-Airy : Démontre l'interaction entre les ondes de diffusion et les ondes de dispersion, montrant une décroissance exponentielle du côté de la chaleur et une dispersion du côté d'Airy.

Signification
L'article affirme sa pertinence en étendant la Méthode de la Transformée Unifiée à une nouvelle classe de problèmes où les équations directrices de chaque côté d'une interface sont de types fondamentalement différents (ordres spatiaux et natures physiques distincts). Les applications précédentes de l'UTM aux problèmes d'interface étaient largement restreintes aux cas où les équations étaient identiques ou ne différaient que par des coefficients scalaires.

En résolvant avec succès ces problèmes d'ordres mixtes, les auteurs démontrent la robustesse de l'UTM pour traiter des structures spectrales complexes issues de relations de dispersion disparates. Ce travail constitue le premier cadre de solution analytique pour de telles interactions, qui sont pertinentes pour la modélisation de phénomènes physiques impliquant des milieux hétérogènes, tels que la croissance tumorale dans le tissu cérébral (comme cité dans l'introduction concernant les gliomes) ou les matériaux composites aux propriétés physiques variables. Les formules explicites dérivées servent de fondation pour des analyses futures, incluant de potentielles extensions aux problèmes non linéaires via des techniques de contraction.