The Kadomtsev-Petviashvili equation in conformal variables for waves over topography

Cet article étend l'approche par application conforme aux écoulements potentiels avec une dépendance transversale faible pour dériver une équation de type Kadomtsev-Petviashvili en variables conformes, permettant de modéliser des ondes de surface sur un fond irrégulier sans nécessiter que la topographie soit une fonction lisse.

L'essentiel

🌊 L'histoire des vagues et du fond marin rugueux

Imaginez que vous êtes un océanographe. Votre travail consiste à prédire comment les vagues voyagent à la surface de l'océan. C'est facile si le fond de la mer est lisse comme une patinoire, mais la réalité est souvent différente : le fond peut être rempli de rochers, de fosses, de monts sous-marins, ou même de formes géométriques très bizarres (comme des marches d'escalier).

En physique, décrire le mouvement de l'eau sur un fond rugueux et irrégulier est un cauchemar mathématique. Les équations habituelles deviennent si complexes qu'elles sont impossibles à résoudre, surtout si le fond n'est pas une courbe douce mais une forme "cassée" ou discontinue.

🪄 La baguette magique : la "Carte Conformale"

C'est ici que les auteurs de ce papier (David Andrade et Marcelo Flamarion) apportent une solution ingénieuse. Ils utilisent une technique mathématique appelée mapping conforme (ou transformation conforme).

Pour faire simple, imaginez que vous avez une carte du monde dessinée sur un morceau de caoutchouc élastique.

• Le problème : Le fond de l'océan (le "caoutchouc") est tout tordu, avec des bosses et des creux.
• La solution : Vous étirez et déformez ce caoutchouc pour le rendre parfaitement plat et lisse, sans jamais le déchirer.

Dans cette nouvelle "carte plate" (appelée variables conformes), les équations de l'eau deviennent beaucoup plus simples à gérer. Le secret ? La rugosité du fond n'est plus dessinée par des lignes brisées, mais par une fonction de lissage (appelée $M$ ou "profondeur effective"). C'est comme si la magie mathématique prenait un fond de mer en béton armé et le transformait en un fond de sable lisse pour faire les calculs, tout en gardant la même physique.

📐 Deux nouvelles recettes pour les vagues

Les chercheurs ont utilisé cette astuce pour créer deux nouvelles équations (des "recettes" pour prédire les vagues) qui ressemblent à l'équation célèbre de Kadomtsev-Petviashvili (KP).

1. La recette pour les fonds qui changent doucement :
   Imaginez une pente de ski qui descend lentement. Même si le terrain est accidenté, la transformation mathématique rend la pente "lisse" dans notre carte magique. Cette équation permet de prédire comment une vague ralentit ou s'accélère en traversant ces zones.

2. La recette pour les petits obstacles :
   Imaginez des petits cailloux dispersés au fond. Si ils sont assez petits, on peut utiliser une autre version simplifiée de l'équation pour voir comment les vagues les contournent.

🎬 Ce que les simulations nous montrent

Les auteurs ont fait tourner ces équations sur ordinateur pour voir ce qui se passe avec une vague qui rencontre un obstacle (comme une série de rectangles sous l'eau, comme des marches géantes).

• Sans le fond (l'océan plat) : La vague avance tout droit, un peu comme un train sur des rails.
• Avec le fond rugueux :
  • La vague ralentit quand elle passe sur l'obstacle (comme une voiture qui monte une côte).
  • Elle se déforme : une grande vague principale est suivie d'une traînée de petites vagues oscillantes (comme le sillage d'un bateau).
  • Le plus important : les deux nouvelles équations donnent des résultats très précis, même si le fond est "cassé" ou irrégulier, là où les anciennes méthodes auraient échoué ou donné des résultats faux.

💡 Pourquoi c'est génial ?

Le message principal de ce papier est une révolution pour la modélisation :

Vous n'avez plus besoin d'un fond de mer parfait pour utiliser vos modèles de vagues.

Grâce à cette méthode, on peut prendre n'importe quel profil de fond (même très bizarre, même fait de marches carrées) et le transformer en une "profondeur effective" lisse. Cela permet d'utiliser des équations simples et rapides pour simuler des situations réalistes et complexes, comme les tsunamis traversant des archipels ou les vagues sur des récifs coralliens irréguliers.

En résumé : C'est comme si les chercheurs avaient inventé un traducteur universel qui transforme un paysage marin chaotique en un terrain de jeu lisse, permettant de prédire le comportement des vagues avec une précision incroyable, même là où la géographie est la plus sauvage.

Résumé technique

1. Problématique

L'article s'attaque au défi de modéliser la propagation d'ondes de surface faiblement non linéaires et faiblement dispersives sur des fonds marins (topographies) complexes, en particulier lorsque ces fonds ne sont pas des fonctions lisses ou régulières.

Les méthodes classiques de cartographie conforme sont puissantes pour les écoulements potentiels bidimensionnels (une dimension horizontale), mais leur extension à des dépendances transversales faibles (problèmes 3D) reste difficile. De plus, les modèles réduits existants (comme les équations de Korteweg-de Vries - KdV - ou Kadomtsev-Petviashvili - KP) supposent souvent que la topographie est une fonction lisse et lentement variable, ce qui limite leur applicabilité à des fonds rugueux ou discontinus (ex: fosses sous-marines, monts sous-marins).

L'objectif est de dériver une équation de type KP formulée en variables conformes qui puisse gérer des topographies « rugueuses » tout en restant un modèle réduit efficace.

2. Méthodologie

A. Cadre théorique et Cartographie Conforme
Les auteurs partent des équations d'Euler pour un écoulement potentiel incompressible et irrotationnel en 3D. Ils utilisent l'approche de cartographie conforme étendue au domaine 3D (développée précédemment par Andrade et Nachbin, 2018).

• Une transformation conforme $F$ est appliquée pour « aplatir » le domaine fluide physique (borné par une surface libre $\eta$ et un fond rugueux $H(x)$) vers un domaine de référence uniforme (une bande infinie).
• Les coordonnées physiques $(x, z)$ sont transformées en coordonnées harmoniques $(\xi, \zeta)$.
• Le Jacobien de la transformation, noté $J(\xi, \zeta)$, et sa valeur au fond $M(\xi) = \sqrt{J(\xi, 0)}$, jouent un rôle central. $M(\xi)$ encode l'information topographique.

B. Développement Asymptotique
Les auteurs appliquent une expansion asymptotique basée sur deux paramètres sans dimension :

• $\epsilon = a/h$ (paramètre de non-linéarité).
• $\mu = h/l$ (paramètre de dispersion).
• Une dépendance transversale faible est introduite via une échelle $\gamma = O(\mu)$.

Ils dérivent deux versions distinctes de l'équation KP selon l'hypothèse faite sur la topographie :

1. Cas de la profondeur lentement variable :

   • Hypothèse : Le coefficient $M(\xi)$ est une fonction lente de la variable conforme ($M(\xi) = M(\mu^2\xi)$).
   • Méthode : Utilisation des invariants de Riemann et d'une expansion perturbative pour les ondes unidirectionnelles.
   • Résultat : Une équation KP avec coefficients variables dépendant de $M(\xi)$ et de ses dérivées.

2. Cas de l'amplitude de topographie faible :

   • Hypothèse : La topographie est une petite perturbation ($M(\xi) = 1 + \epsilon m(\xi)$).
   • Méthode : Utilisation d'une expansion proposée récemment par Ludu et al. (2025).
   • Résultat : Une équation KP où les effets de la topographie apparaissent comme des termes de perturbation linéaires.

C. Concept de « Profondeur Effective »
Un point clé de la méthodologie est la définition de la profondeur effective $d(x) = M(\xi(x))$. Grâce à la régularité inhérente de la fonction $M(\xi)$ (qui est une version lissée de la topographie physique via une convolution implicite dans la transformation conforme), les auteurs montrent que même si la topographie physique $H(x)$ est discontinue (ex: rectangles), la profondeur effective $d(x)$ utilisée dans le modèle réduit reste lisse et bien définie.

3. Résultats Principaux

A. Équations Dérivées
Les auteurs obtiennent deux équations KP principales :

• Équation (41) : Pour une profondeur lentement variable. Elle inclut des termes de dérivées de $M$ et permet de retrouver l'équation KdV classique de Grimshaw/Johnson pour une profondeur variable, mais avec la profondeur effective à la place de la profondeur physique.
• Équation (53) : Pour une topographie de faible amplitude. Elle se rapproche des modèles de Ludu et al. (2025) mais dans un cadre conforme.

B. Validation Numérique
Des simulations numériques ont été réalisées en utilisant une méthode pseudo-spectrale (FFT) avec des conditions aux limites périodiques.

• Configuration : Propagation d'une impulsion gaussienne sur une topographie composée de rectangles (discontinus), illustrée à la Figure 1.
• Comparaison : Les auteurs comparent l'évolution temporelle de l'onde selon :
  1. L'équation KP pour profondeur variable (Eq. 41).
  2. L'équation KP pour faible amplitude (Eq. 53).
  3. L'équation KP classique (fond plat).
• Observations :
  • La topographie ralentit l'onde incidente.
  • L'interaction entre la topographie et la dispersion génère une trainée oscillatoire (wake) derrière le pulse principal.
  • La forme de cette trainée diffère significativement entre les deux modèles KP dérivés, montrant que le choix du modèle dépend de la nature de la variation de profondeur.

4. Contributions Clés et Signification

• Gestion des topographies « rugueuses » : La contribution majeure est la démonstration que les modèles réduits (KdV/KP) peuvent être appliqués à des topographies non lisses, voire discontinues, à condition d'utiliser la profondeur effective dérivée de la cartographie conforme. Cela lève une restriction majeure des modèles classiques qui exigeaient des fonctions $H(x)$ lisses.
• Extension des modèles existants : L'article fournit une extension cohérente des modèles d'ondes dispersives faiblement non linéaires, unifiant la théorie des écoulements potentiels 2D avec la dynamique 3D transversale.
• Nouveauté asymptotique : L'analyse asymptotique pour les profondeurs lentement variables dans le cadre des variables conformes est présentée comme une nouvelle approche, distincte des expansions classiques en coordonnées physiques.
• Implication pratique : Les résultats suggèrent que les ingénieurs et océanographes peuvent utiliser des profils de profondeur « quelconques » dans les modèles réduits existants, à condition de les remplacer par leur profondeur effective correspondante, élargissant ainsi le champ d'application de ces modèles pour la modélisation des tsunamis, des vagues internes ou de la génération d'ondes sur des fonds complexes.

En résumé, ce travail établit un pont rigoureux entre la géométrie complexe des fonds marins et les modèles d'ondes simplifiés, en utilisant la puissance de la cartographie conforme pour lisser mathématiquement les irrégularités physiques tout en conservant leur effet dynamique sur la propagation des ondes.