Existence et approximation numérique d'un système de Boussinesq unidimensionnel à coefficients variables sur un intervalle fini

Informations de base

• ID de l'article : 2510.10829
• Titre : Existence and numerical approximation of a one-dimensional Boussinesq system with variable coefficients on a finite interval
• Auteurs : Deissy Marcela Pizo, Juan Carlos Muñoz Grajales (Universidad del Valle, Colombie)
• Classification : math.NA, cs.NA, math.AP
• Date de publication : 14 octobre 2025 (prépublication arXiv)
• Lien de l'article : https://arxiv.org/abs/2510.10829

Résumé

Cet article étudie le problème du caractère bien posé d'un modèle non linéaire dispersif à coefficients variables, qui décrit l'évolution de la propagation des ondes de surface dans un canal d'eau peu profonde de longueur finie unidimensionnelle présentant une topographie de fond irrégulière. Pour compléter l'analyse théorique, les auteurs utilisent un solveur numérique développé dans leurs travaux antérieurs (méthode des éléments finis) pour approximer les solutions du modèle sur un intervalle spatial borné, en considérant diverses valeurs de paramètres et formes fonctionnelles des coefficients variables du système de Boussinesq. De plus, l'article propose des expériences numériques préliminaires pour résoudre un problème inverse : reconstruire la hauteur d'onde initiale et la vitesse du fluide à partir des données mesurées au temps final. Ceci est réalisé en construisant un problème d'optimisation où les conditions initiales sont estimées comme solution minimisant une fonctionnelle objectif, qui quantifie la différence entre l'état final observé et la solution numérique évoluant à partir des conditions initiales d'essai.

Contexte et motivation de la recherche

1. Problèmes à résoudre :
   • Étudier le caractère bien posé mathématique (existence, unicité, dépendance continue) du système de Boussinesq unidimensionnel à coefficients variables
   • Développer des méthodes numériques efficaces pour résoudre ce système à coefficients variables
   • Explorer le problème inverse de reconstruction des conditions initiales à partir des données d'observation au temps final
2. Importance du problème :
   • Le système de Boussinesq est un modèle mathématique important pour décrire la propagation des ondes en eau peu profonde, avec des applications importantes en génie océanique et protection côtière
   • Les domaines aquatiques réels présentent souvent une topographie de fond non uniforme, rendant les modèles à coefficients variables plus proches de la réalité
   • Les problèmes inverses ont une valeur pratique importante dans la prévision des vagues et la surveillance océanique
3. Limitations des méthodes existantes :
   • La plupart des recherches se concentrent sur les systèmes de Boussinesq à coefficients constants
   • L'analyse théorique pour le cas à coefficients variables est relativement insuffisante
   • La recherche sur les problèmes inverses pour ces systèmes dispersifs reste largement inexploitée
4. Motivation de la recherche :
   • Combler les lacunes dans l'analyse théorique du système de Boussinesq à coefficients variables
   • Fournir des modèles mathématiques plus précis et des outils numériques pour les applications pratiques en génie océanique
   • Ouvrir de nouvelles directions de recherche sur les problèmes inverses pour les systèmes d'ondes dispersives

Contributions principales

1. Établissement de la théorie du caractère bien posé du système de Boussinesq unidimensionnel à coefficients variables : preuve de l'existence, l'unicité et la dépendance continue de la solution par rapport aux données initiales
2. Développement d'une méthode d'analyse théorique basée sur la technique de la fonction de Green : transformation du problème en problème de point fixe d'un opérateur non linéaire
3. Fourniture d'une loi de conservation de l'énergie : preuve que la fonctionnelle d'énergie du système reste conservée dans le temps
4. Implémentation d'un solveur numérique par éléments finis efficace : développement d'un schéma numérique stable combinant la bibliothèque FEniCS
5. Première étude du problème inverse pour ce type de système dispersif : proposition d'une méthode de reconstruction des conditions initiales basée sur l'optimisation variationnelle

Détails des méthodes

Définition du problème

Étude du système de Boussinesq à coefficients variables suivant : $\begin{cases} N_t + [(1 + \alpha c^2(\xi)N)V]_\xi - \frac{\beta}{6}N_{\xi\xi t} = 0 \\ V_t + (c(\xi)N)_\xi + \frac{1}{2}\alpha(c^2V^2)_\xi - \frac{\beta}{6}V_{\xi\xi t} = 0 \end{cases}$

où $(\xi,t) \in [0,L] \times [0,T]$, avec les conditions initiales : $N(\xi,0) = N_0(\xi), \quad V(\xi,0) = V_0(\xi)$

et les conditions aux limites de Dirichlet : $N(0,t) = N(L,t) = 0, \quad V(0,t) = V(L,t) = 0$

Méthodes d'analyse théorique

1. Technique de la fonction de Green Réécriture du système sous la forme : $\left(I - \frac{\beta}{6}\partial_\xi^2\right)N_t = -\partial_\xi[(1 + \alpha c^2(\xi)N)V]$$\left(I - \frac{\beta}{6}\partial_\xi^2\right)V_t = -\partial_\xi\left[c(\xi)N + \frac{1}{2}\alpha c^2V^2\right]$

Utilisation de la fonction de Green de l'opérateur $P = I - \frac{\beta}{6}\partial_\xi^2$ : $G(\xi,s) = \frac{1}{2\sqrt{\beta/6}} \frac{\cosh\left(\frac{L-|\xi-s|}{\sqrt{\beta/6}}\right) - \cosh\left(\frac{L-(\xi+s)}{\sqrt{\beta/6}}\right)}{\sinh\left(\frac{L}{\sqrt{\beta/6}}\right)}$

2. Théorème du point fixe Transformation du problème en système d'équations intégrales, définition d'un opérateur non linéaire $A$, et preuve de l'existence et l'unicité de la solution via le théorème du point fixe de Banach.

Méthodes numériques

Discrétisation par éléments finis Utilisation de la méthode de Galerkin, avec la formulation faible : $\partial_t\left(\langle N,v_1\rangle_{L^2} + \frac{\beta}{6}\langle N_\xi, v_{1\xi}\rangle_{L^2}\right) = \langle[(1+\alpha c^2N)V], v_{1\xi}\rangle_{L^2}$

Discrétisation temporelle Utilisation de la méthode $\theta$ ($\theta = 1/2$) pour l'intégration temporelle, assurant la stabilité numérique.

Méthode du problème inverse

Définition de la fonctionnelle objectif : $J(N,V) := \frac{1}{2}\int_0^L \left(|N(\xi,T;N_0,V_0) - N_T(\xi)|^2 + |V(\xi,T;N_0,V_0) - V_T(\xi)|^2\right)d\xi$

Minimisation de cette fonctionnelle à l'aide de l'algorithme L-BFGS-B, avec implémentation de la différentiation automatique via la bibliothèque Dolfin-Adjoint.

Configuration expérimentale

Paramètres numériques

• Domaine spatial : $[-20, 40]$
• Nombre de points de maillage : 3000 points spatiaux
• Résolution spatiale : $\Delta\xi \approx 0.02$
• Pas de temps : $\Delta t = 8/3000 \approx 2.7 \times 10^{-3}$
• Paramètres du modèle : $\alpha = \beta = 0.1$

Conditions initiales

Utilisation d'un profil gaussien localisé : $N(\xi,0) = V(\xi,0) = e^{-(\xi-18)^2}$

Fonctions à coefficients variables

Cas 1 : Oscillation lisse avec perturbation gaussienne $c(\xi) = 1 + 0.3\sin\left(\frac{\pi}{5}\xi\right) + 0.6e^{-(\xi-8)^2}$

Cas 2 : Fonction par morceaux constante (simulation d'un milieu stratifié) $c(\xi) = \begin{cases} 0.8, & \xi < 20 \\ 1.5, & 20 < \xi \leq 22 \\ 2.0, & 22 < \xi \leq 27 \\ \vdots \end{cases}$

Résultats expérimentaux

Résultats théoriques

Théorème 3.1 (Caractère bien posé local) : Il existe $T_0 > 0$ tel que le système admette une solution unique $(N,V)$ dans $L^2_{T_0} \times L^2_{T_0}$, satisfaisant l'estimation : $\|(N,V)\|_{L^2_{T_0} \times L^2_{T_0}} \leq \|(N_0,V_0)\|_{H^1_0 \times H^1_0} T_0^{1/2} \exp\left(\frac{L^{1/2}}{(\beta/6)^{1/2}}DT_0\right)$

Théorème 3.2 (Conservation de l'énergie) : $\|V\|_{L^2}^2 + \alpha\|c|N|^{1/2}V\|_{L^2}^2 + \||c|^{1/2}N\|_{L^2}^2 = \text{constante}$

Résultats des expériences numériques

1. Résolution du problème direct :
   • Simulation réussie de la propagation des ondes dans un milieu à coefficients variables
   • Observation des changements de vitesse d'onde, de la modulation d'amplitude et des effets dispersifs dus aux variations de topographie du fond
   • Démonstration d'une bonne stabilité et précision de la méthode numérique
2. Reconstruction du problème inverse :
   • L'algorithme L-BFGS-B améliore significativement la reconstruction des conditions initiales en 4 itérations
   • La hauteur d'onde et la vitesse du fluide reconstruites convergent progressivement vers le profil cible
   • Validation de l'efficacité de la stratégie d'assimilation variationnelle

Vérification numérique

• La loi de conservation de l'énergie est bien préservée dans les calculs numériques
• La cohérence des résultats à différentes résolutions de maillage valide la convergence de la méthode numérique

Travaux connexes

1. Théorie des systèmes de Boussinesq :
   • Bona & Chen (1998, 2002, 2004) : établissement de la théorie fondamentale pour le cas à coefficients constants
   • Cet article étend la théorie au cas à coefficients variables, comblant une lacune théorique
2. Méthodes numériques :
   • Application de la méthode des éléments finis aux équations d'ondes dispersives
   • Cet article implémente un solveur efficace combinant la bibliothèque FEniCS
3. Problèmes inverses :
   • Méthodes traditionnelles pour les problèmes inverses d'équations d'ondes
   • Cet article applique pour la première fois l'optimisation variationnelle au problème inverse du système de Boussinesq

Conclusions et discussion

Conclusions principales

1. Établissement réussi des fondations mathématiques du système de Boussinesq unidimensionnel à coefficients variables
2. Développement de méthodes numériques efficaces pour la résolution
3. Ouverture d'une nouvelle direction de recherche sur les problèmes inverses pour les systèmes dispersifs
4. Fourniture d'outils théoriques et computationnels pour les applications pratiques en génie océanique

Limitations

1. Aspect théorique :
   • Seul le caractère bien posé local a été prouvé ; l'existence globale reste à étudier
   • Absence d'analyse théorique sur l'unicité et la stabilité du problème inverse
2. Aspect numérique :
   • Les expériences du problème inverse ne considèrent que le cas idéal sans bruit
   • Absence d'analyse de la sensibilité au bruit d'observation
3. Aspect applicatif :
   • Limitation au cas unidimensionnel ; l'extension à deux dimensions est nécessaire pour les applications pratiques
   • Choix relativement simples des fonctions à coefficients variables

Directions futures

1. Développement théorique :
   • Théorie du caractère bien posé global
   • Analyse de l'unicité et de la stabilité du problème inverse
   • Théorie de la régularisation
2. Améliorations numériques :
   • Méthodes robustes de reconstruction avec données bruitées
   • Extension aux cas de dimension supérieure
   • Intégration des techniques d'apprentissage automatique
3. Applications pratiques :
   • Validation avec des données océaniques réelles
   • Intégration aux systèmes d'observation océanique

Évaluation approfondie

Points forts

1. Contribution théorique significative : première étude systématique du caractère bien posé du système de Boussinesq à coefficients variables
2. Innovation méthodologique : utilisation ingénieuse de la technique de la fonction de Green pour traiter les problèmes à coefficients variables
3. Implémentation numérique complète : fourniture d'un cadre complet de résolution numérique
4. Problématique novatrice : ouverture d'un nouveau domaine de recherche sur les problèmes inverses pour les systèmes dispersifs
5. Rédaction claire : dérivations mathématiques rigoureuses et exposition claire

Insuffisances

1. Profondeur théorique limitée : établissement uniquement de la théorie locale, absence d'analyse globale
2. Vérification expérimentale insuffisante : expériences du problème inverse relativement simples, manque de vérification dans des scénarios complexes
3. Applicabilité à valider : nécessité de davantage de vérifications avec des données réelles
4. Absence d'analyse de complexité computationnelle : pas d'analyse de la complexité algorithmique

Impact

1. Valeur académique : pose les fondations pour le développement théorique des systèmes dispersifs à coefficients variables
2. Perspectives d'application : potentiel d'application en génie océanique, prévision des vagues et domaines connexes
3. Contribution méthodologique : la combinaison de la méthode de la fonction de Green et de l'optimisation variationnelle offre de nouvelles perspectives pour les problèmes connexes

Domaines d'application

1. Modélisation de la propagation des ondes en eau peu profonde
2. Analyse des vagues en génie côtier
3. Problèmes d'inversion dans la surveillance de l'environnement océanique
4. Assimilation de données en océanographie numérique

Références

L'article cite 20 références importantes couvrant la théorie des systèmes de Boussinesq, les méthodes numériques, l'analyse fonctionnelle et d'autres domaines connexes, fournissant une base théorique solide pour la recherche.