S'abonner à un flux RSS
 

Jean-Michel Tanguy/SujetENTPE2022/BENSAC - ECORSE - FLANDRIN : Différence entre versions

De Wikibardig
(Solution semi-analytique par homotopie)
(Solution semi-analytique par homotopie)
Ligne 342 : Ligne 342 :
 
On part de l'équation : <math> H(x) \phi_{xx}(x) + H'(x) \phi_{x}(x) + k^2 H(x) \phi(x) = 0 </math>. <br>
 
On part de l'équation : <math> H(x) \phi_{xx}(x) + H'(x) \phi_{x}(x) + k^2 H(x) \phi(x) = 0 </math>. <br>
  
Sachant que <math> H(x) = H_0 - sx </math>, la relation d'homotopie est la suivante : <br>
+
Sachant que <math> H(x) = H_0 - sx </math> et <math> k = k_0 \sqrt{\frac{H_0}{H(x)}} </math>, la relation d'homotopie est la suivante : <br>
<center> <math> (1-p)(\phi_{xx} - u_{0,xx}) + p[(H_0 - sx)\phi_{xx} - s\phi_{x} + k^2 (H_0-sx) \phi] =0 </math> </center> <br>
+
<center> <math> (1-p)(\phi_{xx} - u_{0,xx}) + p[(H_0 - sx)\phi_{xx} - s\phi_{x} + k_0^2 H_0 \phi] =0 </math> </center> <br>
  
 
En prenant <math> u_0 = 0 </math> pour simplifier les calculs, on obtient : <br>
 
En prenant <math> u_0 = 0 </math> pour simplifier les calculs, on obtient : <br>
<center> <math> (1-p)\phi_{xx} + p((H_0 - sx) \phi_{xx} - s\phi_{x} + k^2 (H_0-sx) \phi) =0 </math> </center> <br>
+
<center> <math> (1-p)\phi_{xx} + p((H_0 - sx) \phi_{xx} - s\phi_{x} + k_0^2 H_0 \phi) =0 </math> </center> <br>
  
  
 
La solution d'homotopie peut alors être transformée comme suit : <br>
 
La solution d'homotopie peut alors être transformée comme suit : <br>
<center> <math> \phi_{0,xx} + p\phi_{1,xx} + \cdots = - p[(H_0 -sx -1)\phi_{0,xx} - s\phi_{0,x} + k^2(H_0 -sx) \phi_0] - p^2[(H_0 -sx -1)\phi_{1,xx} - s\phi_{1,x} + k^2(H_0 -sx) \phi_1] - \cdots </math> </center> <br>
+
<center> <math> \phi_{0,xx} + p\phi_{1,xx} + \cdots = - p[(H_0 -sx -1)\phi_{0,xx} - s\phi_{0,x} + k_0^2 H_0 \phi_0] - p^2[(H_0 -sx -1)\phi_{1,xx} - s\phi_{1,x} + k_0^2 H_0 \phi_1] - \cdots </math> </center> <br>
  
 
La résolution de cette équation est alors celle du système linéaire suivant, par liberté de la famille <math> (1,p,p^2,p^3, \cdots) </math> : <br>  
 
La résolution de cette équation est alors celle du système linéaire suivant, par liberté de la famille <math> (1,p,p^2,p^3, \cdots) </math> : <br>  
Ligne 356 : Ligne 356 :
 
   \begin{eqnarray}
 
   \begin{eqnarray}
 
     \phi_{0,xx}(x) & = & 0 &&&(3.1)\\
 
     \phi_{0,xx}(x) & = & 0 &&&(3.1)\\
     \phi_{1,xx}(x) & = & s\phi_{0,x} - (H_0 -sx - 1) \phi_{0,xx}-k^2 (H_0 -sx) \phi_0(x) &&&(3.2)\\
+
     \phi_{1,xx}(x) & = & s\phi_{0,x} - (H_0 -sx - 1) \phi_{0,xx}-k_0^2 H_0 \phi_0(x) &&&(3.2)\\
     \phi_{2,xx}(x) & = & s\phi_{1,x} - (H_0 -sx - 1) \phi_{1,xx}-k^2 (H_0 -sx) \phi_1(x) &&&(3.3)\\
+
     \phi_{2,xx}(x) & = & s\phi_{1,x} - (H_0 -sx - 1) \phi_{1,xx}-k_0^2 H_0 \phi_1(x) &&&(3.3)\\
 
     &.& \\
 
     &.& \\
 
     &.& \\
 
     &.& \\
Ligne 393 : Ligne 393 :
 
! 1
 
! 1
 
| style="vertical-align:center" |On utilise l'équation <math> (3.2) </math> : <br>
 
| style="vertical-align:center" |On utilise l'équation <math> (3.2) </math> : <br>
<center> <math> \displaystyle \frac{\partial^2 \phi_1}{\partial x^2} = s\frac{\partial \phi_0}{\partial x} - (H_0 -sx - 1) \frac{\partial^2 \phi_0}{\partial x^2}-k^2 (H_0 -sx) \phi_0 </math> </center> <br>
+
<center> <math> \displaystyle \begin{align}
Donc d'après les calculs précédents : <br>
+
& \frac{\partial^2 \phi_1}{\partial x^2} = s\frac{\partial \phi_0}{\partial x} - (H_0 -sx - 1) \frac{\partial^2 \phi_0}{\partial x^2}-k_0^2 H_0 \phi_0 \\
<center> <math> \begin{align}
+
\iff & \frac{\partial \phi_1}{\partial x} = - (H_0 -sx - 1) \frac{\partial \phi_0}{\partial x}-k_0^2 H_0 \int \phi_0 .dx + \alpha_3' \\
& \frac{\partial^2 \phi_1}{\partial x^2} = s\phi_{0,x} - k^2 (H_0 -sx) \phi_0 \\
+
\iff & \phi_1 = - s \int \phi_0.dx - (H_0 -sx - 1) \phi_0 -k_0^2 H_0 \iint \phi_0 .dxdx + \alpha_3'x + \beta_3' \\
\iff & \frac{\partial \phi_{1}}{\partial x} = s\frac{ik}{1-ikL} x - k^2 (H_0 -sx) \left( \frac{ik}{1 - ikL}x + 1 \right) \\  
+
\end{align} </math> </center>
+
 
+
\iff & \phi_{1} = - (H_0 -sx - 1) \phi_0 - s \int \phi_0.dx -k^2 \left[(H_0 -sx) \iint \phi_0.dxdx + 2s \iiint \phi_{0}.dxdxdx \right] + \alpha_3'x + \beta_3' \\
+
 
\end{align} </math> </center> <br>
 
\end{align} </math> </center> <br>
Or : <math> \qquad \iint \phi_{0}.dxdx = \frac{x^2}{2} + \frac{ikx^3}{6(1 - ikL)} </math> <br>
+
Or : <math> \left\{ \begin{align}
Donc : <br>  
+
& \phi_{0} = \frac{ik}{1 - ikL}x +1 \\
<center> <math> \phi_{1}(x) = -k^2 \left( \frac{x^2}{2} + \frac{ikx^3}{6(1 - ikL)} \right) + \alpha_1' x + \beta_1' </math> </center>
+
& \int \phi_0 .dx = \frac{ik}{2(1-ikL)}x^2 + x \\
| <center> <math> \phi_{1}(x = 0) = \beta_1' = 0 </math> </center>
+
& \iint \phi_0.dxdx = \frac{ik}{6(1-ikL)}x^3 + \frac{x^2}{2} \\
 +
\end{align} \right. </math> <br>
 +
<br>
 +
D'où : <br>
 +
<center> <math> \begin{align}  
 +
\phi_1(x) = & - s \left( \frac{ikx^2}{2(1-ikL)} + x \right) - (H_0 -sx - 1) \left( \frac{ikx}{1 - ikL} + 1 \right) \\
 +
& - k_0^2 H_0 \left( \frac{ikx^3}{6(1-ikL)} + \frac{x^2}{2} \right) + \alpha_3'x + \beta_3' \end{align} </math> </center>
 +
| <center> <math> \begin{align}
 +
& \phi_{1}(0) = 0 \\
 +
\iff & - (H_0 - 1) + \beta_3' = 0 \\
 +
\iff & \beta_3' = H_0 -1 \\ \end{align} </math> </center>
 
| On a : <math> \qquad \qquad \qquad \qquad \frac{\partial \phi_{1}(L)}{\partial x} = ik\phi(L) </math> <br>  
 
| On a : <math> \qquad \qquad \qquad \qquad \frac{\partial \phi_{1}(L)}{\partial x} = ik\phi(L) </math> <br>  
 
C'est-à-dire : <br>
 
C'est-à-dire : <br>

Version du 8 juin 2022 à 14:49

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


Sommaire

Introduction

Avant.jpg
Après.jpg

Contextualisation

Dans la deuxième partie de son rapport publié le 28 février 2022, le GIEC alertait que le réchauffement climatique allait entrainer la plus rapide montée des eaux depuis 3000 ans. Cette montée menace de modifier profondément le trait de côte de nombreux pays, pouvant entrainer le déplacement de nombreuses populations, ainsi que la destruction d’ouvrages et d’habitations. La France est d'ailleurs l’un des pays européens les plus menacés par les inondations côtières et par cette augmentation du niveau des océans.

Ainsi, nous allons étudier et modéliser les effets de la houle sur le littoral avec la montée des eaux et afin de répondre aux besoins de prévoir dans les décennies qui viennent les modifications du trait de côte et les destructions possibles.

Le modèle qui va nous permettre de réaliser la modélisation de la houle est le modèle de Berkhoff qui est constitué d’une équation aux dérivées partielles. Nous étudierons différents cas, que nous résoudrons par une méthode analytique puis par homotopie.



Le modèle de Berkhoff

Le modèle de Berkhoff a pour but de combiner les effets de la réfraction et de la diffraction pour les vagues se propageant sur le fond marin dans un seul modèle, aussi appelé équation de pente douce. Son équation fonctionne sous l’hypothèse de houles de faibles amplitudes et sous l’approximation mild-slope (ce qui signifie qu’on considère que les fonds marins sont en pente douce). Berkhoff est le premier auteur à établir, dans la théorie potentielle linéaire, une équation de propagation de la houle couplant les effets de la diffraction et de la réfraction induits par les variations topographiques.

En bi-dimensionnel, il a pour équation : $ \boxed{\displaystyle\nabla.(CC_g\nabla \phi)+k^2CC_g\phi=0} $
où :
- $ \phi $ est le potentiel,
- $ k $ est le nombre d’onde fonction de la profondeur $ H $ et de la fréquence $ \omega $, par la relation implicite $ \omega^2=gk \tanh(kH) $ ,
- $ C $ est la célérité de l’onde,
- $ C_g $ est la célérité de groupe des vagues.

L'équation de pente douce est souvent utilisée dans l'ingénierie côtière pour calculer les changements de champ de vagues près des ports et des côtes. Certaines études portent sur sa réécriture pour simplifier sa résolution numérique. Nous choisirons la méthode d’homotopie pour résoudre cette équation.

Le principe de l'homotopie

L' homotopie est un « chemin » qui permet de passer d'une solution connue simple à une solution complète, au moyen d'un paramètre $ p $ variant entre 0 et 1. Le dit « chemin » est alors une déformation continue entre la première estimation de la solution simple (lorsque $ p = 0 $) et la solution complète (lorsque $ p=1 $) du système d'équations à résoudre.

Hypothèses de calcul

Pour simplifier le problème, nous nous intéresserons à la houle longue (dont la longueur d'onde est supérieure à 300 m). Sous cette hypothèse, on peut faire l'approximation que : $ C=C_g=\sqrt{gH} $.

Lorsque $ H = constante $, l'équation de Berkhoff est simplifiée comme suit :

$ \boxed{\displaystyle\Delta \phi+k^2\phi=0} $

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

Remarque générale : on traite le cas de la houle, donc l'entrée est située en aval, et la sortie en amont.

1er cas Canal monodimensionnel plat avec sortie libre en amont

La première géométrie étudiée est celle d'un canal monodimensionnel plat. On lui affecte une longueur L.

On modélise à son entrée, par l'aval, une onde de fréquence unitaire $ \phi = 1 $ (condition aux limites de Dirichlet), et une sortie libre en amont $ \phi_{x} =ik\phi $ (condition aux limites de Robin).

Par soucis de facilité, on impose ces deux conditions uniquement en $ x = 0 $. Il est possible de démontrer que cette situation est équivalente au cas où ces deux conditions sont appliquées aux deux extrémités du domaine, à savoir en $ x = 0 $ et en $ x = L $.


Simplification de l'équation de Berkhoff :

Dans cette situation, l'équation de Berkhoff simplifée peut être écrite, en domaine unidimensionnel :

$ \boxed{\displaystyle\frac{\partial^2 ϕ}{\partial x^2} +k^2ϕ = 0} $


Solution analytique

On peut réécrire l’équation $ \displaystyle \frac{\partial^{2}\phi}{\partial x^{2}} + k^{2}\phi = 0 $ sous sa forme caractéristique : $ X^{2} + k^{2} = 0 $.

En résolvant cette équation du second degré, on trouve que son discriminant $ \Delta = -4k^{2} $ est négatif. L'équation caractéristique admet donc deux racines complexes, qui sont $ X_{1} = - ik $ et $ X_{2} = ik $.

La solution générale est de la forme : $ \phi(x) = \alpha \, \exp^{X_{1}x} + \beta \, \exp^{X_{2}x} $.

On détermine les inconnues $ \alpha $ et $ \beta $ à l’aide des conditions aux limites énoncées précédemment : $ \left\{ \begin{align} \bullet &\text{En aval :} &\phi(x=0) & = 1 \quad &\Longrightarrow& \quad \alpha+\beta = 1 \\ \bullet &\text{En amont :} &\dfrac{\partial ϕ}{\partial x}(x=L) & = ik \, \phi(x = L) \quad &\Longrightarrow &\quad \alpha = 0 \quad \text{et} \quad \beta=1 \\ \end{align} \right. $

Solution semi-analytique par homotopie

Hypothèses et conditions limites

Soit $ p\in{[0,1]} $ un paramètre qui varie de 0 à 1.

En choisissant la dérivée seconde comme opérateur linéaire, la relation d'homotopie (avec $ p \in [0 \, ; 1] $) s'écrit :

$ (1-p)(\phi_{xx}-u_{0,xx}) + p(\phi_{xx} + k^2\phi) = 0 $


On rappelle que : $ \left\{ \begin{align} & \phi(x,p) = \sum_{n=0}^\infty p^n\phi_{n} \qquad \qquad \qquad \qquad \qquad &(\Leftrightarrow \phi(x,p) = \phi_0(x) + p\phi_1(x) + p^2\phi_2(x) + p^3\phi_3(x) + \cdots) \\ & \dfrac{\partial^2 ϕ}{\partial x^2}(x,p)= \phi_{xx}(x,p) = \sum_{n=0}^\infty p^n\phi_{n,xx} \qquad \qquad \qquad \qquad \qquad &(\Leftrightarrow \phi_{xx}(x,p) = \phi_{0,xx}(x) + p\phi_{1,xx}(x) + p^2\phi_{2,xx}(x) + p^3\phi_{3,xx}(x) + \cdots) \\ \end{align} \right. $


On part d'une solution initiale ($ p = 0 $) nulle : $ u_{0} = 0 $. La solution complète est obtenue lorsque $ p = 1 $.

La relation d'homotopie se simplifie alors comme suit :

$ \phi_{xx} + pk^2\phi = 0 $

Résolution aux différents ordres

La relation d’homotopie simplifiée ci-dessus peut-être développée de la manière suivante :

$ \phi_{0,xx}(x) + p\phi_{1,xx}(x) + p^2\phi_{2,xx}(x) + p^3\phi_{3,xx}(x) + \cdots = -k^2(p\phi_0(x) + p^2\phi_1(x) + p^3\phi_2(x) + p^4\phi_3(x) + \cdots) $


La résolution de cette équation est alors celle du système linéaire suivant, par liberté de la famille $ (1,p,p^2,p^3,\cdots) $ :

$ \left\{ \begin{eqnarray} \phi_{0,xx}(x) & = & 0 &&&(1.1)\\ \phi_{1,xx}(x) & = & -k^2 \phi_0(x) &&&(1.2)\\ \phi_{2,xx}(x) & = & -k^2 \phi_1(x) &&&(1.3)\\ \phi_{3,xx}(x) & = & -k^2 \phi_2(x) &&&(1.4)\\ &.& \\ &.& \\ &.& \\ \end{eqnarray} \right. $

On s’intéresse maintenant aux différents ordres de $ p $ :

Ordres Equation Conditions aux limites
En $ x = 0 $ En $ x = L $
0 On utilise l'équation $ (1.1) $ :
$ \begin{align} \phi_{0,xx} = \frac{\partial^{2}\phi_0}{\partial x^{2}} = 0 && \iff && \frac{\partial\phi_0}{\partial x} = \alpha_1 \\ && \iff && \phi_{0}(x) = \alpha_1 x + \beta_1 \end{align} $
$ \begin{align} \phi_0(0) = 1 \qquad \quad & \Longrightarrow & \beta_1 = 1 \quad \\ \left. \begin{eqnarray} \frac{\partial\phi_0}{\partial x} = ik\phi_0 \\ \frac{\partial\phi_0}{\partial x} = \alpha_1 \end{eqnarray} \right\} \quad & \Longrightarrow & \alpha_1 = ik\phi_0 \\ \end{align} $ $ \begin{align} & \text{On a :} & \alpha_1 = ik\phi_0 \quad \text{et} \quad \beta_1 = 1 \\ & \text{Or :} & \phi(x = L) = \alpha_1 L + \beta_1 \quad & \Longrightarrow & \; \alpha_1 = \frac{ik}{1-ikL} \qquad \\ \end{align} $

Ainsi :

$ \boxed{\displaystyle \phi_{0}(x) = 1 - \frac{kx}{2} + \frac{ikx}{2}} $
1 On utilise l'équation $ (1.2) $ :
$ \begin{align} \phi_{1,xx} = \frac{\partial^2\phi_{1}}{\partial x^2} = -k^2\phi_0 & \iff& \frac{\partial\phi_{1}}{\partial x} = -k^2 \int \phi_{0}.dx + \alpha_1' \\ & \iff& \phi_{1} = - k^2\iint \phi_{0}.dxdx + \alpha_1' x + \beta_1' \\ \end{align} $

Or : $ \qquad \iint \phi_{0}.dxdx = \frac{x^2}{2} + \frac{ikx^3}{6(1 - ikL)} $
Donc :

$ \phi_{1}(x) = -k^2 \left( \frac{x^2}{2} + \frac{ikx^3}{6(1 - ikL)} \right) + \alpha_1' x + \beta_1' $
$ \phi_{1}(x = 0) = \beta_1' = 0 $
On a : $ \qquad \qquad \qquad \qquad \frac{\partial \phi_{1}(L)}{\partial x} = ik\phi(L) $

C'est-à-dire :

$ -k^2 \left(\frac{ikL^2}{2(1-ikL)}+L\right) + \alpha_1' = ik \left(-k^2 \left(\frac{ikL^3}{6(1-ikL)} + \frac{L^2}{2} \right) + \alpha_1'L \right) $

$ \Longrightarrow \qquad \alpha_1' = \frac{k}{2} - \frac{ik}{3} $

Ainsi :

$ \boxed{\displaystyle \phi_{1}(x) = -k^2 \left( \frac{x^2}{2} + \frac{ikx^3}{6(1 - ikL)} \right) + \left( \frac{k}{2} - \frac{ik}{3} \right) x} $

On continue jusqu’à l’ordre n.
On peut représenter les courbes de la solution par homotopie des $ \phi $ aux différents ordres de $ p $ par rapport à la solution analytique.

Voici la représentation pour les premiers ordres :

Cas1 bis.gif

2e cas Canal monodimensionnel plat avec réflexion totale en amont

On étudie maintenant le cas d'un canal monodimensionnel plat de longueur L, cette fois avec réflexion totale en amont.

On modélise, en aval, un flux $ \phi_{x} = ik(2-\phi) $, et la réflexion totale en amont $ \phi_{x} = 0 $ .


Simplification de l'équation de Berkhoff :
Dans cette situation, l'équation de Berkhoff prend la même expression que dans le 1er cas :

$ \boxed{\displaystyle\frac{\partial^2 ϕ}{\partial x^2} +k^2ϕ = 0} $

Solution analytique

On reprend l'exact même raisonnement qu'au 1er cas, pour trouver la forme de la solution générale : $ \phi(x) = \alpha \, \exp^{-ikx} + \beta \, \exp^{ikx} $. On détermine les inconnues $ \alpha $ et $ \beta $ à l’aide des nouvelles conditions aux limites : $ \left\{ \begin{align} \bullet &\text{En aval :} &\dfrac{\partial ϕ}{\partial x}(x=0) = ik(2- \phi(x = 0)) \quad &\iff& \quad -ik \, \alpha + ik \, \beta = ik(2- (\alpha+\beta)) \quad \quad &\Longrightarrow& \quad \beta = 1 \\ \bullet &\text{Réflexion totale en amont :} &\dfrac{\partial ϕ}{\partial x}(x=L)=0 \quad &\Longleftrightarrow &\quad -ik\, \alpha \, e^{-ikL}+ik \, \beta \, e^{ikL} = 0 \quad &\Longrightarrow & \quad \alpha = e^{2ikL} \\ \end{align} \right. $

Solution semi-analytique par homotopie

On utilise le même raisonnement qu'au 1er cas.

Il convient d'étudier plus en détail la relation (vraie pour $ x = 0 $) : $ \dfrac{\partial ϕ}{\partial x} = ik(2- \phi) $.

On rappelle que : $ \left\{ \begin{align} & \phi(x,p) = \sum_{n=0}^\infty p^n\phi_{n} \qquad \qquad \qquad \qquad \qquad &(\Leftrightarrow \phi(x,p) = \phi_0(x) + p\phi_1(x) + p^2\phi_2(x) + p^3\phi_3(x) + \cdots) \\ & \dfrac{\partial ϕ}{\partial x}(x,p)= \phi_{x}(x,p) = \sum_{n=0}^\infty p^n\phi_{n,x} \qquad \qquad \qquad \qquad \qquad &(\Leftrightarrow \phi_{x}(x,p) = \phi_{0,x}(x) + p\phi_{1,x}(x) + p^2\phi_{2,x}(x) + p^3\phi_{3,x}(x) + \cdots) \\ \end{align} \right. $

Ainsi (lorsque $ x = 0 $) : $ \quad \phi_{0,x}(x) + p\phi_{1,x}(x) + p^2\phi_{2,x}(x) + \cdots = 2ik - ik(\phi_0(x) + p\phi_1(x) + p^2\phi_2(x) + \cdots) \quad \iff \quad \left\{ \begin{eqnarray} \phi_{0,x}(x) & = & ik(2-\phi_0(x)) &&&(2.1)\\ \phi_{1,x}(x) & = & -ik\phi_1(x) &&&(2.2)\\ \phi_{2,x}(x) & = & -ik\phi_2(x) &&&(2.3)\\ &.& \\ &.& \\ &.& \\ \end{eqnarray} \right. $

On s’intéresse maintenant aux différents ordres de $ p $ :

Ordres Equation Conditions aux limites
En $ x = L $ En $ x = 0 $
0 On utilise l'équation $ (1.1) $ :
$ \begin{align} \phi_{0,xx} = \frac{\partial^{2}\phi}{\partial x^{2}} = 0 && \iff && \frac{\partial\phi}{\partial x} = \alpha_2 \\ && \iff && \phi_{0}(x) = \alpha_2 x + \beta_2 \end{align} $
$ \begin{align} \frac{\partial\phi_0(L)}{\partial x} = 0 \quad \Longrightarrow \quad \alpha_2 = 0 \end{align} $
$ \begin{align} \frac{\partial\phi_0(0)}{\partial x} = ik(2- \phi_0(0)) \quad & \iff & \quad &0 = ik(2-\beta_2) \\ & \iff & &\beta_2 = 2 \\ \end{align} $

D'où :
$ \boxed{\displaystyle \phi_0(x) = 2} $
1 On utilise l'équation $ (1.2) $ :
$ \begin{align} \phi_{1,xx} = \frac{\partial^2\phi_{1}}{\partial x^2} = -k^2\phi_0 && \iff && \frac{\partial\phi_{1}}{\partial x} = -k^2 \int \phi_{0}.dx + \alpha_2' \\ && \iff && \phi_{1} = - k^2\iint \phi_{0}.dxdx + \alpha_2' x + \beta_2' \end{align} $

Or : $ \qquad \qquad \qquad \iint \phi_{0}.dxdx = x^2 $
Donc :

$ \phi_{1}(x) = -k^2 x^2 + \alpha_2' x + \beta_2' $
$ \begin{align} \frac{\partial\phi_1(L)}{\partial x} = 0 \quad \Longrightarrow \quad \alpha_2' = 2k^2L \\ \end{align} $
$ \begin{align} \text{On utilise l'equation } (2.2) : \qquad \quad & \frac{\partial\phi_1(0)}{\partial x} = ik\phi_1(0) \\ \iff & \quad 2k^2 L = ik\beta_2' \\ \iff & \quad \beta_2' = -2ikL \\ \end{align} $
D'où :
$ \boxed{\displaystyle \phi_1(x) = -k^2 x^2 +2k^2 Lx -2ikL} $

Voici la représentation sur Maxima pour les premiers ordres :

Cas2 bis.gif

3e cas Canal monodimensionnel de pente constante avec sortie libre en amont

On étudie maintenant le cas d'un canal monodimensionnel de longueur L, dont la pente de fond est constante égale à $ s $.

On modélise à son entrée, par l'aval, une onde de fréquence unitaire $ \phi = 1 $ (condition aux limites de Dirichlet), et une sortie libre en amont $ \phi_{x} =ik\phi $ (condition aux limites de Robin).


Simplification de l'équation de Berkhoff :

Dans cette situation, la profondeur du canal $ H $ n'est pas une constante : $ H = H(x) $. On conserve l'hypothèse $ C = C_g = \sqrt{gH} $. Pour simplifier l'équation de Berkhoff, il est nécessaire de repartir de son expression initiale :

$ \begin{align} & \nabla.(CC_g\nabla \phi)+k^2CC_g\phi=0 \\ \iff & \nabla.(gH(x)\, \nabla \phi) + k^2 gH(x)\phi = 0 \end{align} $

Finalement, l'équation de Berkhoff dans ce cas s'écrit :

$ \boxed{\displaystyle H \frac{\partial^2 ϕ}{\partial x^2} + \frac{\partial H}{\partial x} \frac{\partialϕ}{\partial x} + k^2 Hϕ = 0} $


Hypothèses :

On prendra $ H(x) = H_0 - sx $, où $ H_0 $ est la hauteur à l'origine, et $ s $ représente la pente.
On considère également que $ k = k_0 \sqrt{\frac{H_0}{H(x)}} $

Solution analytique

Les conditions aux limites :
$ \left\{ \begin{align} & \bullet \; \phi(0) = 1 \\ & \bullet \; \phi_{x}(x=L) = ik\phi(x=L) \\ \end{align} \right. $

Hypothèses et simplifications dans ce cas précis :

En effectuant le changement de variable $ z(x) = H_0 - sx $, l’équation prend la forme :

$ \displaystyle z(x) \frac{\partial^2 \phi}{\partial x^2} – s \frac{\partial \phi}{\partial x}+ H_0 k_0^2 \phi = 0 $


Afin de résoudre cette équation, on la simplifie en suivant les étapes ci-dessous :

$ \begin{align} & z(x) \frac{\partial^2 \phi}{\partial x^2} – s \frac{\partial \phi}{\partial x}+ H_0 k_0^2 \phi = 0 \\ \iff & z(x) \frac{\partial}{\partial x}(\frac{\partial \phi}{\partial z} \frac{\partial z}{\partial x}) -s \frac{\partial \phi}{\partial z} \frac{\partial z}{\partial x}+ H_0 k_0^2 \phi = 0 \\ \iff & -s z(x) \frac{\partial}{\partial x}(\frac{\partial \phi}{\partial z}) + s^2 \frac{\partial \phi}{\partial z}+ H_0 k_0^2 \phi = 0 \\ \end{align} $


$ \displaystyle \begin{align} \text{En posant } F = \frac{\partial \phi}{\partial z}\text{, on trouve :} \qquad & s^2 \frac{\partial \phi}{\partial z} -s z(x) \frac{\partial F}{\partial x} \frac{\partial z}{\partial x} + H_0 k_0^2 \phi = 0 \\ \iff & z(x) \frac{\partial^2 \phi}{\partial z^2} + \frac{\partial \phi}{\partial z} + \frac{H_0 k_0^2}{s^2} \phi= 0 \\ \end{align} $

En posant $ \alpha^2 = \frac{k_0^2 H_0}{s^2} $, on retrouve une équation de Bessel : $ \displaystyle z(x) \frac{\partial^2 \phi}{\partial z^2} + \frac{\partial \phi}{\partial z} + \alpha^2 \phi = 0 $


La solution de cette équation est alors $ \displaystyle \phi(z) = B J_0 (2 \alpha \sqrt{z}) + C Y_0 (2 \alpha \sqrt{z}) $, où : $ \left\{ \begin{align} & J_0 \text{ est une fonction de Bessel de 1ère espèce à l’ordre 0} \\ & Y_0 \text{ est une fonction de Bessel de 2ème espèce à l'ordre 0} \\ & B \text{ et } C \text{ sont des constantes à déterminer à l'aide des conditions limites} \\ \end{align} \right. $.


Avec les conditions limites, on obtient :

$ BJ_{0}(2\alpha \sqrt H_{0}) + CY_{0}(2\alpha \sqrt H_{0}) = 1 $

et :

$ BJ_{1}(2\alpha \sqrt H_{1}) + CY_{1}(2\alpha \sqrt H_{1}) = i(BJ_{0}(2\alpha \sqrt H_{1}) + CY_{0}(2\alpha \sqrt H_{1})) $


Ce qui revient à dire que :

$ B = \frac{(Y_{1}^{2}-Y_{0}^{2})}{J_{0}(Y_{1}^{2}-Y_{0}^{2}) - Y_{0}^{0}(J_{1}^{2}-iJ_{0}^{2})} $

et :

$ C = \frac{- (J_{1}^{2} - iJ_{0}^{2})}{J_{0}(Y_{1}^{2}-Y_{0}^{2}) - Y_{0}^{0}(J_{1}^{2}-iJ_{0}^{2})} $

Solution semi-analytique par homotopie

Dans ce cas précis pour un souci de réalisme, il faut veiller à ce que la hauteur de la houle ne soit pas supérieure à celle de la profondeur en eau.
La programmation de cette solution dans ce cas se fera selon le modèle de Bessel.
Nous prenons de manière arbitraire la valeur $ s =\frac{1}{200} $ de la pente du fond pour la superposition des solutions analytique et par homotopie dans l’animation temporelle.


On part de l'équation : $ H(x) \phi_{xx}(x) + H'(x) \phi_{x}(x) + k^2 H(x) \phi(x) = 0 $.

Sachant que $ H(x) = H_0 - sx $ et $ k = k_0 \sqrt{\frac{H_0}{H(x)}} $, la relation d'homotopie est la suivante :

$ (1-p)(\phi_{xx} - u_{0,xx}) + p[(H_0 - sx)\phi_{xx} - s\phi_{x} + k_0^2 H_0 \phi] =0 $

En prenant $ u_0 = 0 $ pour simplifier les calculs, on obtient :

$ (1-p)\phi_{xx} + p((H_0 - sx) \phi_{xx} - s\phi_{x} + k_0^2 H_0 \phi) =0 $


La solution d'homotopie peut alors être transformée comme suit :

$ \phi_{0,xx} + p\phi_{1,xx} + \cdots = - p[(H_0 -sx -1)\phi_{0,xx} - s\phi_{0,x} + k_0^2 H_0 \phi_0] - p^2[(H_0 -sx -1)\phi_{1,xx} - s\phi_{1,x} + k_0^2 H_0 \phi_1] - \cdots $

La résolution de cette équation est alors celle du système linéaire suivant, par liberté de la famille $ (1,p,p^2,p^3, \cdots) $ :

$ \left\{ \begin{eqnarray} \phi_{0,xx}(x) & = & 0 &&&(3.1)\\ \phi_{1,xx}(x) & = & s\phi_{0,x} - (H_0 -sx - 1) \phi_{0,xx}-k_0^2 H_0 \phi_0(x) &&&(3.2)\\ \phi_{2,xx}(x) & = & s\phi_{1,x} - (H_0 -sx - 1) \phi_{1,xx}-k_0^2 H_0 \phi_1(x) &&&(3.3)\\ &.& \\ &.& \\ &.& \\ \end{eqnarray} \right. $


Ordres Equation Conditions aux limites
En $ x = 0 $ En $ x = L $
0 On utilise l'équation $ (3.1) $ :
$ \begin{align} \phi_{0,xx} = \frac{\partial^{2}\phi_0}{\partial x^{2}} = 0 && \iff && \frac{\partial\phi_0}{\partial x} = \alpha_3 \\ && \iff && \phi_{0}(x) = \alpha_{3} x + \beta_3 \end{align} $
$ \begin{align} \phi_0(0) = 1 \qquad \quad & \Longrightarrow & \beta_3 = 1 \quad \\ \left. \begin{eqnarray} \frac{\partial\phi_0}{\partial x} = ik\phi_0 \\ \frac{\partial\phi_0}{\partial x} = \alpha_3 \end{eqnarray} \right\} \quad & \Longrightarrow & \alpha_3 = ik\phi_0 \\ \end{align} $ $ \begin{align} & \text{On a :} & \alpha_3 = ik\phi_0 \quad \text{et} \quad \beta_3 = 1 \\ & \text{Or :} & \frac{\partial \phi_0}{\partial x}(x = L) = ik(\alpha_3 L + \beta_3) \quad & \Longrightarrow & \; \alpha_3 = \frac{ik}{1-ikL} \qquad \\ \end{align} $

Ainsi :

$ \boxed{\displaystyle \phi_{0}(x) = \frac{ik}{1 - ikL}x + 1} $
1 On utilise l'équation $ (3.2) $ :
$ \displaystyle \begin{align} & \frac{\partial^2 \phi_1}{\partial x^2} = s\frac{\partial \phi_0}{\partial x} - (H_0 -sx - 1) \frac{\partial^2 \phi_0}{\partial x^2}-k_0^2 H_0 \phi_0 \\ \iff & \frac{\partial \phi_1}{\partial x} = - (H_0 -sx - 1) \frac{\partial \phi_0}{\partial x}-k_0^2 H_0 \int \phi_0 .dx + \alpha_3' \\ \iff & \phi_1 = - s \int \phi_0.dx - (H_0 -sx - 1) \phi_0 -k_0^2 H_0 \iint \phi_0 .dxdx + \alpha_3'x + \beta_3' \\ \end{align} $

Or : $ \left\{ \begin{align} & \phi_{0} = \frac{ik}{1 - ikL}x +1 \\ & \int \phi_0 .dx = \frac{ik}{2(1-ikL)}x^2 + x \\ & \iint \phi_0.dxdx = \frac{ik}{6(1-ikL)}x^3 + \frac{x^2}{2} \\ \end{align} \right. $

D'où :

$ \begin{align} \phi_1(x) = & - s \left( \frac{ikx^2}{2(1-ikL)} + x \right) - (H_0 -sx - 1) \left( \frac{ikx}{1 - ikL} + 1 \right) \\ & - k_0^2 H_0 \left( \frac{ikx^3}{6(1-ikL)} + \frac{x^2}{2} \right) + \alpha_3'x + \beta_3' \end{align} $
$ \begin{align} & \phi_{1}(0) = 0 \\ \iff & - (H_0 - 1) + \beta_3' = 0 \\ \iff & \beta_3' = H_0 -1 \\ \end{align} $
On a : $ \qquad \qquad \qquad \qquad \frac{\partial \phi_{1}(L)}{\partial x} = ik\phi(L) $

C'est-à-dire :

$ -k^2 \left(\frac{ikL^2}{2(1-ikL)}+L\right) + \alpha_1' = ik \left(-k^2 \left(\frac{ikL^3}{6(1-ikL)} + \frac{L^2}{2} \right) + \alpha_1'L \right) $

4e cas Vague sphérique générée par une source périodique sinusoïdale

On étudie désormais l'évolution de la surface libre d’une vague sphérique générée par une source périodique sinusoïdale, dans un domaine infini en grande profondeur. La source ponctuelle est appliquée autour d'un cercle de rayon $ r_{0} $ centré sur un domaine circulaire de rayon $ R $ qui laisse sortir librement cette onde en $ r = R $.

On modélise, en aval, une onde $ \phi = 1 $, et en amont, une onde $ \phi_x = ik\phi $.


Simplification de l'équation de Berkhoff :

Dans ce cas, l'équation de Berkhoff se transforme en équation de Helmholtz : $ \displaystyle \Delta \phi(r) + k^2\phi(r)=0 $.

Etant donné que le problème est caractérisé par une symétrie de révolution (il est indépendant de $ \theta $), la relation ci-dessus s'écrit, en coordonnées polaires :

$ \boxed{\displaystyle \frac{\partial^2 \phi}{\partial r^2} + \frac{1}{r} \frac{\partial \phi}{\partial r} + k^2 \phi = 0} $

Solution analytique

L'équation précédente étant une équation de Bessel, on utilise à nouveau la méthode de résolution du cas précédent, et on peut écrire :

$ \phi(r)= \alpha J_{0}(r)+\beta Y_{0}(r) $

avec : $ \left\{ \begin{align} & J_0 \text{ est une fonction de Bessel de 1ère espèce à l’ordre 0} \\ & Y_0 \text{ est une fonction de Bessel de 2ème espèce à l'ordre 0} \\ & \alpha \text{ et } \beta \text{ sont des constantes à déterminer à l'aide des conditions limites} \\ \end{align} \right. $.


On doit alors résoudre :

$ \phi_{rr}+\dfrac{1}{r}\phi_r + k^2\phi=0 $

Avec les conditions initiales suivantes :
$ \left\{ \begin{align} & \bullet \; \phi^{ r = r_{0}} = 1 \\ & \bullet \; \phi_r^{r=R} = ik\phi^{r=R} \\ \end{align} \right. $

On peut alors déterminer les constantes $ \; \alpha \; et \; \beta $ telles que : $ \left\{ \begin{align} & \alpha J_{0}(r_{0})+\beta Y_{0}(r_{0}) = 1 \\ & \alpha J_0'(R)+\beta Y_0'(R)=ik\ (\alpha J_{0}(R)+\beta Y_{0}(R)) \\ \end{align} \right. $


Or on sait que : $ \left\{ \begin{align} & \displaystyle\frac{\partial r^nJ_n(r)}{\partial r} = r^nJ_{n-1}(r) \\ & J_{-1}(r) = -J_1(r) \\ \end{align} \right. $

On trouve alors :

$ \alpha = \displaystyle\frac{Y_1(R) + ikY_0(R) }{Y_1(R)J_0(r_0) + ikY_0(R)J_0(r_0) - J_1(R)Y_0(r_0) - ikJ_0(R)} $

et :

$ \beta = \displaystyle\frac{1 }{Y_0(r_0)} \left(1-J_0(r_0)\frac{Y_1(R) + ikY_0(R) }{Y_1(R)J_0(r_0) + ikY_0(R)J_0(r_0) - J_1(R)Y_0(r_0) - ikJ_0(R)}\right) $

D'où : $ \displaystyle\phi(r)=\frac{Y_1(R) + ikY_0(R) }{Y_1(R)J_0(r_0) + ikY_0(R)J_0(r_0) - J_1(R)Y_0(r_0) - ikJ_0(R)}J_0(r)+\frac{1}{Y_0(r_0)}\left(1-J_0(r_0)\frac{Y_1(R) + ikY_0(R) }{Y_1(R)J_0(r_0) + ikY_0(R)J_0(r_0) - J_1(R)Y_0(r_0) - ikJ_0(R)}\right)Y_0(r) $

Solution semi-analytique par homotopie

Les hypothèses appliquées : On choisira les valeurs numériques suivantes: $ \bullet r_{0} = 1 m \bullet R = 100 m \bullet k = 0.1 m^{-1} $ Les résultats sont affichés suivant un rayon de $ r_{0} à R $ :

  1. le premier graphique superpose la solution analytique avec la solution par homotopie
  2. le second consiste en une application temporelle avec une fréquence de $ \omega=\sqrt{gk} $ car nous sommes en profondeur infiniée.
  3. on produira également une animation 2D si possible

La relation d'homotopie s'écrit de la manière suivante :

$ (1-p)(\phi_{rr} - u_{0,rr}) + p \left( \phi_{rr} + \frac{1}{r} \phi_{r} + k^2 \phi \right) =0 $

En prenant $ u_0 = 0 $ pour simplifier les calculs, on obtient :

$ \phi_{rr} + p \left( \frac{1}{r} \phi_{r} + k^2 \phi \right) =0 $


La solution d'homotopie peut alors être transformée comme suit :

$ \phi_{0,rr} + p\phi_{1,rr} + p^2 \phi_{2,rr} + \cdots = - p(\frac{1}{r} \phi_{0,r} + k^2 \phi_0) - p^2 (\frac{1}{r} \phi_{1,r} + k^2 \phi_1) - p^3 (\frac{1}{r} \phi_{2,r} + k^2 \phi_2) - \cdots $

La résolution de cette équation est alors celle du système linéaire suivant, par liberté de la famille $ (1,p,p^2,p^3, \cdots) $ :

$ \left\{ \begin{eqnarray} \phi_{0,rr}(r) & = & 0 &&&(4.1)\\ \phi_{1,rr}(r) & = & -\frac{1}{r} \phi_{0,r}(r) - k^2 \phi_0(r) &&&(4.2)\\ \phi_{2,rr}(r) & = & -\frac{1}{r} \phi_{1,r}(r) - k^2 \phi_1(r) &&&(4.3)\\ &.& \\ &.& \\ &.& \\ \end{eqnarray} \right. $


Ordres Equation Conditions aux limites
En $ r = r_0 $ En $ r = R $
0 On utilise l'équation $ (4.1) $ :
$ \begin{align} \phi_{0,rr} = \frac{\partial^{2}\phi_0}{\partial r^{2}} = 0 && \iff && \frac{\partial\phi_0}{\partial r} = \alpha_4 \quad & (*) \\ && \iff && \phi_{0}(r) = \alpha_{4} r + \beta_4 & (**) \\ \end{align} $
$ \begin{align} \phi_0(r_0) = 1 \qquad \quad & \iff & \alpha_{4} r_0 + \beta_4 = 1 \quad \\ & \iff & \beta_4 = 1 - \alpha_4 r_0 \\ \end{align} $
On a : $ \left\{ \begin{eqnarray} \frac{\partial\phi_0}{\partial r}(R) = ik\phi_0(R) \\ \frac{\partial\phi_0}{\partial r} = \alpha_4 \end{eqnarray} \right. \iff \alpha_4 = ik(\alpha_4 R + \beta_4) $

Or : $ \quad \beta_4 = 1 - \alpha_4 r_0 \quad $ d'où :

$ \left\{ \begin{eqnarray} \alpha_4 = \frac{ik}{1+ ik(r_0 - R)} \\ \beta_4 = \frac{1 - ikR}{1 + ik(r_0 - R)} \end{eqnarray} \right. $

Ainsi :

$ \boxed{\displaystyle \phi_{0}(r) = \frac{ik}{1+ ik(r_0 - R)}r + \frac{1 - ikR}{1 + ik(r_0 - R)}} $
1 On utilise l'équation $ (4.2) $ :
$ \displaystyle \frac{\partial^2 \phi_1}{\partial r^2} = -\frac{1}{r} \phi_{0,r}(r) - k^2 \phi_0(r) $

Or, d'après $ (*) $ et $ (**) $ : $ \displaystyle \phi_{1,rr} = -\frac{1}{r} \alpha_4 - k^2 (\alpha_4r + \beta_4) $
D'où :

$ \begin{align} & \phi_{1,r} = -ln(r) \alpha_4 - \alpha_4 \frac{k^2 r^2}{2} - k^2 \beta_4 r + \alpha_4' \\ \iff & \phi_1 = -(rln(r) + r)\alpha_4 - \alpha_4 \frac{k^2 r^3}{6} - \beta_4 \frac{k^2 r^2}{2} + \alpha_4' r + \beta_4' \\ \end{align} $

Avec : $ \displaystyle \alpha_4 = \frac{ik}{1+ ik(r_0 - R)} $ et $ \displaystyle \beta_4 = \frac{1 - ikR}{1 + ik(r_0 - R)} $

$ \displaystyle \begin{align} & \phi_{1}(r_0) = 0 \\ \iff & -(r_0 ln(r_0) + r_0)\alpha_4 - \alpha_4 \frac{k^2 r_0^3}{6} - \beta_4 \frac{k^2 r_0^2}{2} + \alpha_4' r_0 + \beta_4' = 0 \\ \iff & \beta_4' = (r_0 ln(r_0) + r_0)\alpha_4 + \alpha_4 \frac{k^2 r_0^3}{6} + \beta_4 \frac{k^2 r_0^2}{2} - \alpha_4' r_0 \\ \end{align} $

Pour la suite, on pose : $ \left\{ \begin{align} & F(r) = -\alpha_4 ln(r) - \frac{\alpha_4 k^2 r^2}{2} - k^2 \beta_4 r \\ & G(r) = -\alpha_4 r (ln(r) -1) - \frac{\alpha_4 k^2 r^3}{6} - \frac{\beta_4 k^2 r^2}{2} \\ \end{align} \right. $

On a :

$ \displaystyle \begin{align} & \frac{\partial \phi_{1}(R)}{\partial r} = ik\phi(R) \\ \iff & F(R) + \alpha_4' = ik(G(R) + \alpha_4' R + \beta_4') \\ \end{align} $

D'où :

$ \left\{ \begin{align} & \alpha_4' = \frac{ik(G(R) - G(r_0)) - F(R)}{1 - ik(r_0 - R)} \\ & \beta_4' = \frac{ik(G(R) - G(r_0)) - F(R)}{1 - ik(r_0 - R)} r_0 + G(r_0) \\ \end{align} \right. $

Outils personnels