S'abonner à un flux RSS
 

Jean-Michel Tanguy/SujetENTPE2022/CALLOCH - DECHAVASSINE - DIANOUX : Différence entre versions

De Wikibardig
(Solution analytique)
(Cas n°4 : Vague sphérique dans un domaine infini)
Ligne 70 : Ligne 70 :
  
 
== Cas n°4 : Vague sphérique dans un domaine infini ==
 
== Cas n°4 : Vague sphérique dans un domaine infini ==
 +
On étudie dans ce cas une surface dans un domaine de profondeur considérée comme infinie. Cette surface est agité par une onde originaire d'un cercle de rayon <math>r_0</math>. Le domaine est circulaire et l'onde sort librement en <math>r=R</math>. Cela se traduit par les conditions aux limites suivantes : <br />
 +
<math>\phi^{r_0} = 1</math> et <math>\phi_r^R = ik\phi</math>

Version du 24 mai 2022 à 10:02

Pour retourner à la page énoncé du sujet : Utilisateur:Jean-Michel Tanguy/SujetENTPE2022

Sommaire

Cadre de l'étude

On cherche à modéliser une houle dans quatre cas différents. Pour cela, on utilise le modèle de Berkhoff, dont l'équation est :$ \nabla.(CC_g\nabla \phi)+k^2CC_g\phi=0 $
$ \phi $ est le potentiel, $ k $ est le nombre d’onde, $ C $ est la célérité de l’onde et $ C_g $ est la célérité de groupe des vagues.
On se place dans le domaine des ondes longues, ce qui entraîne que $ C = C_g = \sqrt{gH} $, avec g l'accélération de la pesanteur et H la profondeur du canal. L'équation de Berkhoff devient ainsi :

$ \nabla^2 \phi + k^2\phi = 0 $.

L'évolution de la hauteur de la houle en fonction du temps est donnée par : $ h(x,t)=\Re \left (\phi e^{-i\omega t} \right ) $

Nous allons résoudre cette équation en utilisant la méthode d'homotopie, plus précisément la méthode HAM "Homotopie Perturbation Method". Le principe est de pouvoir passer, par une transformation continue, d'un milieu complexe dont il est difficile de connaitre le comportement à un milieu plus simple.
Elle consiste à résoudre l'équation $ (1-p)L(\phi-\psi)+pH(p)({d^2 \phi \over {dx^2}} + k^2 \phi)=0 $, avec L un opérateur linéaire (ici on prendra la dérivée seconde par rapport à x) et $ \psi $ la solution initiale. p est un paramètre qui varie entre 0 et 1, il représente la progression de la solution initiale ($ p = 0 $) à la solution finale ($ p = 1 $). La méthode HAM implique que $ H(p) = 1 $

Cas n°1 : Canal plat ouvert

On considère un canal à fond plat, ouvert aux deux extrémités, de longueur L et de profondeur H. On adopte la notation $ \phi_x $ pour la dérivée de $ \phi $ par rapport à x et $ \phi^x $ pour la valeur de $ \phi $ en x.
En aval, il entre une onde de fréquence unitaire $ \phi=1 $ (condition de Dirichlet) et sortie libre amont $ \phi_x=ik\phi $ (condition de Robin).
Diagramme canal cas 1.png

Résolution analytique

Soit à résoudre : $ \phi_{xx} + k^2\phi = 0 $
La solution est donc de la forme : $ \phi = Ae^{-ikx} + B e^{ikx} $
Les conditions aux limites de notre problème sont : $ \phi^0 = 1 $ et $ \phi_x^L = \phi^L $
Donc $ B = 1 $ et $ A = 0 $
Ainsi, $ \phi = e^{ikx} $. L'enjeu est maintenant de retrouver cette solution en utilisant la méthode d'homotopie.


Résolution homotopique

Avec une solution initiale nulle, la relation d'homotopie devient :

$ (1-p)\phi_xx + k^2\phi = 0 $
On décompose $ \phi $ en somme de fonctions : $ \phi = \sum_{i=0}^{+\infty} p^i\phi_i $ Cette relation est vrai quel que soit p, donc on peut établir que les coefficients de chaque puissance de p sont nuls aussi. On décompose par ordre, c'est-à-dire par puissance de p.

  • Ordre 0

$ \phi_{0,xx} = 0 $ d'où $ \phi_0 = Ax+B $
On sait que :
$ \phi_0(x = 0) = 1 $ donc $ B = 1 $
et $ \phi_{0,x}(x = L) = ik\phi_0(x =L) $ donc $ A = ik(AL + 1) $ donc $ A = {ik \over{1-ikL}} $
Ainsi, $ \phi_0(x) = {ik\over{1-ikL}} x + 1 $

  • Ordre 1

$ \phi_{1,xx} + k^2\phi_0 = 0 \iff \phi_1 = -k^2\iint \phi_0 dx^2 + Ax + B $
Nos conditions aux limites sont toujours : $ \phi_0^1=0 \ et \ \phi^L_{1,x}=ik\phi^L_1 $. Il vient : $ \phi_1= −k^2L{(k^2L^2+3ikL−3) \over {3(1−ikL)^2}}x−k^2({ik\over{6(1−ikL)}}x^3+ \dfrac {1} {2}x^2) $

  • Ordre n

On a : $ \forall n, \phi_{n,xx} + k^2\phi_{n-1} = 0 $
Ainsi, $ \phi_n = -k^2\iint \phi_{n-1} dx^2 +A_nx +B $. (B = 0 )
Comme $ \phi_{n,x}^L = ik\phi_{n-1}^L, A_n = {1\over {1-ikL}}(k^2\int_0^L\phi_{n-1}dx^2 - ik^3\iint_0^L\phi_{n-1}dx^2) $
Finalement, on obtient :
$ h(x,t) = \sum p^n\Re(\phi_ne^{-i\omega t}) $


On peut visualiser l'évolution de la solution homotopique selon la valeur de p grâce au logiciel Maxima:
Ezgif.com-gif-maker.gif
Sur le graphique ci-dessus, la solution analytique est représentée par les croix rouges. La solution homotopique est représentée par la courbe verte. On illustre ici la variation du paramètre p, de 0 (solution en bleue) à 1 (où elle est la plus proche de la solution analytique). La méthode par homotopie permet bien de retrouver la solution analytique.

Cas n°2 : Canal plat fermé

On étudie un canal similaire à celui du cas précédent, mais fermé en x=L. Cela se traduit par les conditions aux limites suivantes :
$ \phi_x=ik(2-\phi) $ en x=0 et $ \phi_x = 0 $ en x=L

Solution analytique

On cherche $ \phi $ sous la forme $ \phi = Ae^{-ikx}+Be^{ikx} $, car il y a réflexion totale.
Les conditions aux limites donnent respectivement B=1 et $ A = e^{2ikL} $
On a ainsi $ \phi = e^{ik(2L-x)}+e^{ikx} $

Solution par homotopie

Ordre 0

Cas n°3 : Canal en pente

Cas n°4 : Vague sphérique dans un domaine infini

On étudie dans ce cas une surface dans un domaine de profondeur considérée comme infinie. Cette surface est agité par une onde originaire d'un cercle de rayon $ r_0 $. Le domaine est circulaire et l'onde sort librement en $ r=R $. Cela se traduit par les conditions aux limites suivantes :
$ \phi^{r_0} = 1 $ et $ \phi_r^R = ik\phi $

Outils personnels