On the Dirichlet-Neumann operator for nearly spherical domains

1. Problématique et Contexte

L'article étudie l'opérateur de Dirichlet-Neumann (ODN), noté $G(h)$, associé à des domaines presque sphériques dans $\mathbb{R}^n$. Cet opérateur non local apparaît fréquemment dans les problèmes à frontière libre, notamment en dynamique des fluides, en optimisation de forme et en problèmes inverses.

Le domaine $\Omega$ est défini comme une perturbation de la sphère unité $S^{n-1}$ via une fonction de hauteur (ou d'élévation) $h : S^{n-1} \to \mathbb{R}$ :
$$\partial\Omega = \{ x(1 + h(x)) : x \in S^{n-1} \}, \quad \text{avec } 1+h > 0.$$
L'ODN est défini par la dérivée normale du potentiel harmonique $\Phi$ (solution du problème de Laplace avec condition de Dirichlet $\psi$) sur la frontière $\partial\Omega$. Le but est d'établir des estimations précises de la dépendance analytique de $G(h)$ par rapport à la forme du domaine (c'est-à-dire par rapport à $h$) et d'obtenir des estimations « tames » (linéaires par rapport aux normes les plus élevées) dans les espaces de Sobolev $H^s$.

2. Méthodologie et Innovations Techniques

La preuve repose sur une transformation du problème initial, défini sur un domaine variable $\Omega$, vers un problème défini sur une géométrie fixe (la boule unité $B_1$).

A. Changement de variables et extension harmonique
Les auteurs introduisent un difféomorphisme $\tilde{\gamma} : B_1 \to \Omega$ construit via l'extension harmonique de $h$ (opérateur de Poisson). Cette approche permet de transformer l'équation de Laplace en une équation elliptique à coefficients variables sur la boule fixe :
$$\text{div}(P \nabla u) = 0 \quad \text{dans } B_1, \quad u = \psi \quad \text{sur } S^{n-1},$$
où $u = \Phi \circ \tilde{\gamma}$ et $P$ est une matrice dépendant de $h$. La solution $u$ est ensuite développée en série de puissances de $h$.

B. Espaces de Sobolev non isotropes
L'innovation majeure réside dans l'utilisation d'une classe spécifique d'espaces de Sobolev non isotropes pour les fonctions définies sur la boule $B_1$.

• Régularité tangentielle : On considère une régularité élevée (éventuellement fractionnaire) dans les directions tangentes à la sphère.
• Régularité normale : On ne considère qu'une régularité entière et limitée (degré 1 ou 2) dans la direction normale.
  Cela permet d'éviter l'estimation de dérivées d'ordre élevé dans la direction normale, ce qui est la source principale de la dépendance des constantes par rapport à l'indice de régularité $s$ dans les travaux antérieurs.

C. Estimations de produits asymétriques
Un outil clé est l'établissement d'estimations de produits pour les normes de Sobolev de la forme :
$$\|fg\|_s \le C_{s_0} \|f\|_s \|g\|_{s_0} + C_s \|f\|_{s_0} \|g\|_s,$$
où la constante $C_{s_0}$ est indépendante de $s$ (le paramètre de haute régularité), tandis que $C_s$ dépend de $s$. Cette propriété est cruciale pour prouver que le rayon de convergence de la série de Taylor de l'opérateur ne dépend pas de la régularité $s$ choisie.

3. Résultats Principaux

Les auteurs démontrent deux théorèmes fondamentaux :

Théorème 1.1 (Estimations Tames et Analyticité)
Pour $s \ge 1/2$ et $s_0 > (n-1)/2$, il existe un rayon $\delta > 0$ indépendant de $s$ tel que, pour tout $h$ dans une boule de rayon $\delta$ dans $H^{s_0+1}$, l'opérateur $G(h)$ satisfait l'estimation tame suivante :
$$\|G(h)\psi\|_s \le C_0 \|\psi\|_{s+1} + (\chi_{s>s_0}) C_s (1 + \|h\|_{s+1}) \|\psi\|_{s_0+1}.$$
De plus, l'application $h \mapsto G(h)$ est analytique de l'espace des perturbations vers l'espace des opérateurs linéaires bornés.

• Perte de régularité : La perte de régularité pour la fonction de hauteur $h$ est de 1 (c'est-à-dire que si $h \in H^{s+1}$, le terme correctif dépend de $\|h\|_{s+1}$), ce qui est une amélioration par rapport aux résultats antérieurs (comme dans [5]) où la perte était de $3/2$.
• Indépendance du rayon : Le rayon de convergence $\delta$ et la constante $C_0$ ne dépendent pas de $s$.

Théorème 1.2 (Dérivées de l'opérateur)
Des estimations multilinéaires similaires sont établies pour les dérivées première et seconde de $G$ par rapport à $h$. Ces estimations confirment la régularité analytique de l'opérateur et fournissent des bornes précises pour les termes non linéaires.

4. Contributions Clés et Améliorations Techniques

Par rapport aux travaux précédents (notamment [5] pour la sphère en dimension 3), cet article apporte les améliorations suivantes :

1. Réduction de la perte de régularité : La perte de régularité sur $h$ passe de $3/2$ à $1$. Cela est obtenu grâce à une construction astucieuse du difféomorphisme $\tilde{\gamma}$ qui agit comme un inverse à droite de l'opérateur de trace près de la frontière.
2. Indépendance du rayon d'analyticité : Le rayon $\delta$ de la boule où les estimations tames et l'analyticité sont valables est indépendant de l'indice de régularité $s$. Dans les travaux antérieurs, ce rayon dépendait de $s$, ce qui rendait difficile l'application aux équations d'évolution où l'on fait tendre $s$ vers l'infini.
3. Généralité dimensionnelle : Les résultats sont valables pour toute dimension $n \ge 2$, alors que les résultats antérieurs étaient souvent restreints à $n=3$.
4. Paramètres réels : La régularité $s$ est traitée comme un nombre réel, et non seulement comme un entier.
5. Constantes optimales : La constante $C_0$ devant la norme la plus élevée de $\psi$ est indépendante de $s$.

5. Signification et Impact

Ce travail est significatif pour plusieurs raisons :

• Outils pour les équations d'évolution : Les estimations tames avec un rayon d'analyticité indépendant de $s$ sont essentielles pour prouver l'existence locale de solutions pour des équations d'évolution non linéaires (comme les équations d'eau peu profonde ou les équations d'Euler) utilisant des méthodes de type Nash-Moser ou des schémas itératifs.
• Optimisation de forme : La compréhension fine de la dépendance analytique de l'ODN par rapport à la géométrie du domaine est cruciale pour les problèmes d'optimisation de forme et les problèmes inverses.
• Rigueur technique : La démonstration de l'estimation de produit asymétrique (1.11) dans les espaces de Sobolev sur la boule, en utilisant des semi-normes non isotropes et des cartes locales, résout une difficulté technique majeure liée aux opérateurs d'extension standards dont les constantes dépendent de l'indice de régularité.

En résumé, cet article fournit un cadre analytique robuste et optimisé pour l'étude des opérateurs de Dirichlet-Neumann sur des domaines déformés, levant plusieurs obstacles techniques qui limitaient les applications précédentes en analyse non linéaire.