S'abonner à un flux RSS
 

Utilisateur:Jean-Michel Tanguy/SujetENTPE2020/trinome2: BUCHER-GIORGI-LOTHON : Différence entre versions

De Wikibardig
Ligne 1 : Ligne 1 :
Bon courage et n'hésitez pas à communiquer avec moi par mail: jm.tanguy@shf-hydro.org
+
 
 +
Dans le but de modéliser la houle, nous allons utiliser le  '''modèle de Berkhoff'''. C'est un modèle bi-dimensionnel qui a pour expression:
 +
 
 +
<math> \frac{\partial  }{\partial x_i}(CC_g \frac{\partial \phi}{\partial x_i})+ k^2CC_g\phi = 0 </math> avec ''i'' = {1,2}
 +
 
 +
 
 +
On note:  
 +
 
 +
<math>\phi</math>: le potentiel de vitesse,
 +
<math>C</math>: la célérité de l'onde,
 +
<math>C_g</math>: la célérité de groupe des vagues,
 +
<math>k</math>: le nombre d'onde,
 +
 
 +
 
 +
L'expression du nombre d'onde k est la suivante:
 +
 
 +
<math>\omega^2 = gk\tan(kH) </math>
 +
 
 +
On note:
 +
 
 +
H: la profondeur,
 +
<math>\omega</math>: la fréquence avec <math>\omega = \frac{2\pi}{T}</math>
 +
 
 +
 
 +
Nous nous plaçons dans le cadre des ondes longues, ce qui signifie que <math> C = C_g = \sqrt{gH}</math>.
 +
L'équation devient alors:
 +
 
 +
<math> \frac{\partial  }{\partial x_i}(H \frac{\partial \phi}{\partial x_i})+ k^2H\phi = 0 </math>
 +
 
 +
 
 +
== Étude cas ==
 +
 
 +
=== Cas n°1 : canal uniforme unidimensionnel plat de longueur L avec sortie libre amont===
 +
 
 +
Dans le cadre de l'hypothèse '''''canal uniforme plat''''', on a  <math> H = constante </math>. On a alors:
 +
 
 +
<math>  \frac{\partial  }{\partial x_i}(H \frac{\partial \phi}{\partial x_i})+ k^2H\phi = 0 \Leftrightarrow  H\frac{\partial  }{\partial x_i}( \frac{\partial \phi}{\partial x_i})+ k^2H\phi = 0 \Leftrightarrow \cancel{H}\frac{\partial  }{\partial x_i}( \frac{\partial \phi}{\partial x_i})+ k^2\cancel{H}\phi = 0  \Leftrightarrow \frac{\partial^2 \phi}{\partial x_i^2}+ k^2\phi = 0 </math>
 +
 
 +
On obtient alors une équation différentielle du second degré sans second membre à résoudre. L'équation d’onde est la suivante:
 +
<math>h(x,t)=a_0\cos(kx-\omega t)</math>
 +
 
 +
 
 +
==== Résolution ====
 +
<math> \phi(x)=A\cos(kx)+B\sin(kx) </math>
 +
 
 +
<math> k = \frac{\omega}{c} \Rightarrow k = \frac{2\pi}{T\sqrt{gH_0}}  </math>
 +
 
 +
<math> \phi(0)=1 </math>
 +
   
 +
<math>\phi(L)=ik\phi  </math>
 +
 
 +
<math> \phi(0)=A=1 \phi(L)=\cos(kL)+B\sin(kL)=ik\phi  </math>
 +
 
 +
<math> B=\frac{ik\phi-\cos(kL)}{\sin(kL)} </math>
 +
 
 +
<math> \phi(x)= \cos(kx)+\frac{ik\phi-\cos(kL)}{\sin(kL)}\sin(kx) </math>
 +
 
 +
<math>Re(\phi)= \cos(kx)-\frac{\cos(kL)}{\sin(kL)}\sin(kx)</math>
 +
 
 +
<math> |\phi|=\sqrt{(Re(\phi))^2+\left(\frac{k\phi\sin(kx)}{\sin(kL)}\right)^2}</math>

Version du 21 avril 2020 à 23:44

Dans le but de modéliser la houle, nous allons utiliser le modèle de Berkhoff. C'est un modèle bi-dimensionnel qui a pour expression:

$ \frac{\partial }{\partial x_i}(CC_g \frac{\partial \phi}{\partial x_i})+ k^2CC_g\phi = 0 $ avec i = {1,2}


On note:

$ \phi $: le potentiel de vitesse, $ C $: la célérité de l'onde, $ C_g $: la célérité de groupe des vagues, $ k $: le nombre d'onde,


L'expression du nombre d'onde k est la suivante:

$ \omega^2 = gk\tan(kH) $

On note:

H: la profondeur, $ \omega $: la fréquence avec $ \omega = \frac{2\pi}{T} $


Nous nous plaçons dans le cadre des ondes longues, ce qui signifie que $ C = C_g = \sqrt{gH} $. L'équation devient alors:

$ \frac{\partial }{\partial x_i}(H \frac{\partial \phi}{\partial x_i})+ k^2H\phi = 0 $


Étude cas

Cas n°1 : canal uniforme unidimensionnel plat de longueur L avec sortie libre amont

Dans le cadre de l'hypothèse canal uniforme plat, on a $ H = constante $. On a alors:

$   \frac{\partial  }{\partial x_i}(H \frac{\partial \phi}{\partial x_i})+ k^2H\phi = 0 \Leftrightarrow  H\frac{\partial  }{\partial x_i}( \frac{\partial \phi}{\partial x_i})+ k^2H\phi = 0 \Leftrightarrow \cancel{H}\frac{\partial  }{\partial x_i}( \frac{\partial \phi}{\partial x_i})+ k^2\cancel{H}\phi = 0  \Leftrightarrow \frac{\partial^2 \phi}{\partial x_i^2}+ k^2\phi = 0  $

On obtient alors une équation différentielle du second degré sans second membre à résoudre. L'équation d’onde est la suivante:

$ h(x,t)=a_0\cos(kx-\omega t) $


Résolution

$ \phi(x)=A\cos(kx)+B\sin(kx) $

$ k = \frac{\omega}{c} \Rightarrow k = \frac{2\pi}{T\sqrt{gH_0}} $

$ \phi(0)=1 $

$ \phi(L)=ik\phi $

$ \phi(0)=A=1 \phi(L)=\cos(kL)+B\sin(kL)=ik\phi $

$ B=\frac{ik\phi-\cos(kL)}{\sin(kL)} $

$ \phi(x)= \cos(kx)+\frac{ik\phi-\cos(kL)}{\sin(kL)}\sin(kx) $

$ Re(\phi)= \cos(kx)-\frac{\cos(kL)}{\sin(kL)}\sin(kx) $

$ |\phi|=\sqrt{(Re(\phi))^2+\left(\frac{k\phi\sin(kx)}{\sin(kL)}\right)^2} $

Outils personnels