S'abonner à un flux RSS
 

Jean-Michel Tanguy/SujetENTPE2022/BEAUDET - MIGEON - MOBISSON

De Wikibardig

Sae3.png

Sommaire

Introduction

Contexte

Aujourd'hui, le changement climatique et l'augmentation de la température constituent des enjeux majeurs pour la planète, et plus spécialement pour les populations vivant le long des côtes avec une importante hausse du niveau des océans Les deux principales causes de ce haussement sont la fonte des glaciers continentaux et la dilatation thermique des océans. Selon l'U.S Global Change Reserach Program for the Fourth National Climate Assessment, entre 1900 et 2018, le niveau moyen des océans s'est élevé d'environ 20 cm, et cela s'accélère à partie de 2020 allant jusqu'à plus de 3,5 mm par an. De nombreux scientifiques estiment que le niveau de la mer pourrait avoir augmenté de 1 à 3 mètres d'ici 2100. Une telle hausse du niveau de la mer pourrait entraîner la suppression de terres ainsi que le déplacement de populations côtières.

La houle correspond au mouvements de la mer de façon ondulatoire, provoqués des zones de vents lointaines. L'étude suivante a pour but d'étudier son impact sur les littoraux, entraînant de plus en plus de conséquences extrêmes. Ces études utilisent des modélisations numériques représentant la houle, tel que le modèle de Berkhoff constitué d'Equations aux Dérivées Partielles, en prenant en compte une représentation des phénomènes de shoaling, réfraction, diffraction et réflexion des différentes structures le long des côtes.

Phénomène de Houle

Lorsque la houle se rapproche des côtes la topographie sous-marine modifie la houle à travers différents phénomènes.

Réfraction : C'est le phénomène majeur qui a lieu. Plus la houle se rapproche de la côte et donc plus la profondeur diminue plus les lignes de crêtes vont être parallèles au rivage.

Diffraction : phénomène lorsque la houle rencontre un obstacle plus petit ou une ouverture.


Réflexion : Ce phénomène a lieu lorsque la houle rencontre une structure profonde et large par rapport à la longueur d'onde de cette dernière.

Shoaling : Effet des vagues de surfaces entrant dans des eaux moins profondes modifiant ainsi la profondeur des vagues.

Modèle de Berkhoff

Le modèle de Berkhoff permet de représenter les phénomènes de réfraction, diffraction et de réflexion que subit la houle. Cette équation a été introduite par Berkhoff en 1973 sous la forme :

$    \displaystyle   \nabla.(CC_g\nabla \phi)+k^2CC_g\phi=0  $

Avec : $ \displaystyle \phi $ : le potentiel

k : le nombre d’onde fonction de H (profondeur) et de ω (fréquence) avec $ \omega^2=gk \tanh(kH) $

C : la célérité de l'onde

Cg : célérité de groupe de vagues

Les hypothèses pour ce modèle :

Nous nous placerons pour la suite dans un modèle linéaire et donc nous approximerons :

- C = Cg = $ \sqrt{gH} $

- h(x,t) = $ \R\phi e^{-i\omega t} $ qui correspond à l'évolution dans le temps et la distance de la houle

Homotopie

On a suivante à résoudre :

$ L(u(x)=f(x)+N(u(x), x \in \Omega $

avec comme conditions limites: $ B(u,u_n=0), x \in \Gamma $

où L est un opérateur linéaire, N un opérateur non-linéaire et f les termes complémentaires de l’équation.

L'homotopie de Liao [9] est définie de la manière suivante :

$ (1-p)\left[ L(U(x,t);p)-L(u_0(x,t)) \right ] +H(p)\left[ L(U(x,t);p)-N(U(x,t),p)-f(x) \right] $

$ p \in \left[0,1\right] $ est un paramètre entre 0 et 1, u0 est une estimation initiale de la solution.

Lorsque p=0 on a la solution initiale U=u0 et lorsque p=1, on a la solution exacte.

La transformation de p de 0 à 1 qui à partir de U(x) de la solution estimée donne la solution exacte est issue de la transformation de U(x) en série de Taylor

$ U(x)=u_0(x)+\sum_{m=1}^\infty u_m(x)p^m $ avec $ u_m(x)=\dfrac{1}{m!}\dfrac{\partial^mU(x,t;p)}{\partial p^m}\Bigr]_{p=0} $

La mise en œuvre de l’homotopie nécessite de faire des hypothèses sur la fonction H(p) qui permet notamment un ajustement des non-linéarités. Pour la méthode HAM « Homotopie Perturbation Method », que nous mettrons en œuvre: $ H(p)=1 $.


Premier Cas

Nous nous plaçons dans le cas d'un canal unidimensionnel de longueur L avec 2 conditions limites.


Cylindre.png

La condition de Dirichlet : $ \phi $ = 1 à l'entrée du canal

La condition de Robin : $ \phi_{x} = ik\phi $ à la sortie du canal

Pour simplifier la résolution on suppose que ces conditions sont imposées en \x = 0

On peut alors simplifier l'écriture de l'équation de Berkhoff tel que :

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

Solution analytique

On a l'équation : $ \displaystyle\frac{\partial^2 ϕ}{\partial x^2} +k^2ϕ = 0 $

L'équation caractéristique associée s'écrit : $ x^2+k^2 = 0 $

Le déterminant de cette équation vaut : $ \Delta = -4k^2 < 0 $

La solution est donc complexe et s'écrit sous la forme : $ \phi(x) = A\mathrm{e}^{(-ikx)}+B\mathrm{e}^{(ikx)} $

On peut déterminer les constantes A et B à l'aide des conditions aux limites:

- $ \phi(x=0) = A+B = 1 $

- $ \dfrac{\partial \phi}{\partial x}(x=l) = ik\phi $

On en déduit: $ A = 0 ; B = 1 $

On a alors : $ \phi(x,t) = \mathrm{e}^{(i(kx-\omega t))} $

L'évolution de la hauteur de la houle s'écrit donc : $ h(x,t) = Re(\phi(x,t)) = cos(kx-\omega t) $

Solution par méthode d'homotopie

On obtient la relation par homotopie à l'aide :

- de la dérivée seconde qui correspond à l'opérateur linéaire

- d'une solution initiale nulle $ \displaystyle (1-p)\phi_{xx}+p(\phi_{xx}+k^2\phi)=0 $

On remplace ensuite dans la solution initiale nulle la décomposition en série entière (en rouge) et la dérivée de cette dernière (en vert) :

$ \displaystyle (1-p)({\color{Green}\phi_{0,xx}(x)+p\phi_{1,xx}(x)+p^2\phi_{2,xx}(x)+p^3\phi_{3,xx}(x)+...})+p[{\color{Green}\phi_{0,xx}(x)+p\phi_{1,xx}(x)+p^2\phi_{2,xx}(x)+p^3\phi_{3,xx}(x)}+...+k^2({\color{Red}\phi_0(x)+p\phi_0(x)+p^2\phi_2(x)+p^3\phi_3(x)+...})]=0 $

Il faut ensuite simplifier et écrire cette relation suivant les puissances de p croissantes. Cette relation étant valable quel que soit p, tous les coefficients devant les puissances de p sont donc nuls. ??

  • Ordre 0

On a donc p = 0 et $ \displaystyle \phi_{0,xx}=0 $

On obtient la relation suivante : $ \phi_0=Ax+B $ Avec les conditions initiales on a :

$ \phi_{0}(x=0)=1 $ et donc $ B = 1 $

$ \Leftrightarrow \dfrac{\partial \phi_{0}}{\partial x}(x=L)=ik \phi_{0}(x=L) $

$ \Leftrightarrow A=ik(AL+1) $

$ \Leftrightarrow A=\frac{ik}{1-ikL} $

donc $ \phi_0=\frac{ik}{1-ikL}x+1 $

  • Ordre 1

On a p = 1 et la relation : $ \displaystyle \phi_{1,xx}=0 $ $ \Leftrightarrow \phi_1= -k^2(\frac{ikx^3}{6(1-ikL)}+\frac{x^2}{2})+Cx+D $

A l'aide des conditions aux limites :

- $ \phi_1(x=0) = 0 $

- $ \dfrac{\partial \phi_{0}}{\partial x}(x=L)=ik \phi_1(x=L) $

On obtient: $ ik(-k^2(\frac{ikL^3}{6(1-ikL)}+\frac{L^2}{2})+CL) =-k^2(\frac{ikL^2}{2(1-ikL)}+L)+C $ $ \Leftrightarrow C = -\frac{k^2L(k^2L^2+3ikl-3)}{3(1-ikl)^2} $

On a alors : $ \phi_1 = -k^2(\frac{ikx^3}{6(1-ikL)}+\frac{x^2}{2}) -\frac{k^2L(k^2L^2+3ikl-3)}{3(1-ikl)^2}x $

Ordres 2 et supérieurs

En vue des calculs de plus en plus complexes à réaliser, on utilise le logiciel WSMAXIMA dans le but de trouver $ \phi $ à ces ordres. Afin de les déterminer, on s'aidera des relations :

- Longueur : $ L = 2\lambda $ (m)

- Profondeur : $ H = 40 $ (m)

- Longueur d'onde : $ \lambda = \frac{2\pi}{k} $ (m)

- Nombre d'onde : $ k = \frac{1}{100} $

- Vitesse de l'onde : $ c = \sqrt{gH} $

On superpose graphiquement la solution analytique (rouge) et la solution par méthode d'homotopie de l'ordre 0 à 5.

Cas1 Ordre1.png Cas1 Ordre1b.png Cas1 Ordre2.png

Cas1 Ordre3.png Cas1 Ordre4.png Cas1 Ordre5.png


On constate que la solution par méthode d'homotopie se rapproche de plus en plus de la solution analytique au fur et à mesure que l'ordre augmente.

Deuxième Cas

Ce deuxième correspond de même à un canal de longueur L. Néanmoins, nous avons désormais les conditions aux limites :

- flux en aval : $ \phi_x(x=0) = ik(2-\phi(x=0)) $

- réflexion totale en amont : $ \phi_x(x=L) = 0 $

Solution analytique

On a l'équation : $ \displaystyle\frac{\partial^2 ϕ}{\partial x^2} +k^2ϕ = 0 $

L'équation caractéristique associée s'écrit : $ x^2+k^2 = 0 $

Le déterminant de cette équation vaut : $ \Delta = -4k^2 < 0 $

La solution est donc complexe et s'écrit sous la forme : $ \phi(x) = A\mathrm{e}^{(-ikx)}+B\mathrm{e}^{(ikx)} $

On peut déterminer les constantes A et B à l'aide des conditions aux limites:

- $ \dfrac{\partial \phi}{\partial x}(x=0) = ik(2-\phi(x=0) \Leftrightarrow -ikA+ikB = ik(2-A-B) \Leftrightarrow B = 1 $

- $ \dfrac{\partial \phi}{\partial x}(x=l) = 0 \Leftrightarrow --ikA\mathrm{e}^{(-ikL)}+ikB\mathrm{e}^{(ikL)} = 0 \Leftrightarrow A = \mathrm{e}^{(2ikL)} $

On a donc : $ h(x,t) = Re(\phi(x,t)) = cos(k(2L-x) - \omega t) $


Solution par méthode d'homotopie

On procède de la même manière que pour la cas 1 en trouvant $ \phi $ aux différents ordres.

  • Ordre 0

On a pour ce cas p = 0 On a alors la relation : $ \phi_{0,xx} = \dfrac{\partial^2 \phi_1}{\partial x^2} = 0 $

$ \Leftrightarrow \phi = Ax+B $ A l'aide des conditions aux limites:

- $ \dfrac{\partial \phi_0}{\partial x}(x=0) = ik(2-\phi_0(x=0)) \Leftrightarrow ik(2-B) = 0 \Leftrightarrow B = 2 $

- $ \dfrac{\partial \phi_0}{\partial x}(x=L) = 0 \Leftrightarrow A = 0 $

On a donc $ \phi_0 = 2 $

  • Ordre 1

On a pour ce cas p = 1 On a alors la relation : $ \phi_{1,xx} = \dfrac{\partial^2 \phi_1}{\partial x^2}+k^2\phi_0 = 0 \Leftrightarrow \phi_1 = -k^2x^2+Cx+D $

A l'aide des conditions aux limites :

- $ \dfrac{\partial \phi_1}{\partial x}(x=0) = ik\phi_1(x=0) \Leftrightarrow ikD = 2k^2L \Leftrightarrow D = 2ikL $

- $ \dfrac{\partial \phi_1}{\partial x}(x=L) = 0 \Leftrightarrow C = 2k^2L $

On a donc $ \phi_1 = -k^2x^2+2k^2Lx+2ikL $

Ordres 2 et supérieurs

Dans le but de trouver $ \phi $ aux ordres 2 et supérieurs, on utilisera le logiciel WSMAXIMA. On pourra donc, par principe de superposition des solutions analytiques et de celle par méthode d'homotopie, vérifier les solutions obtenues.

On superpose la solution analytique (rouge) et la solution par méthode d'homotopie des ordres 0 à 5.

Cas2 Ordre0.png Cas2 Ordre1.png Cas2 Ordre2.png

Cas2 Ordre3.png Cas2 Ordre4.png Cas2 Ordre5.png

On constate que la solution par méthode d'homotopie se rapproche de plus en plus de la solution analytique au fur et à mesure que l'ordre augmente.

Troisième Cas

Pour ce troisième cas on se place dans un domaine de longueur L avec un fond dont la pente est constante. Nous avons deux conditions pour ce modèle qui sont $ \phi $ = 1 à l'entrée et $ \phi_{x} = ik\phi $ à la sortie.

On pose aussi les conditions suivantes :

- k non constant : $ \displaystyle k(x) = k_0\sqrt{\dfrac{H_0}{H(x)}} = k_0\sqrt{\dfrac{H_0}{H_0 - sx}} $

- $ CC_g = gH(x) $

On obtient alors pour l'équation de Berkhoff :

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

Solution analytique

A l'aide du changement de variable $ z(x) = H_0 - sx $ on va obtenir une équation de type bessel et par identification on pourra obtenir $ \alpha $.

On a alors :

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

On pose:

- $ J_0 $ la fonction de Bessel de 1ère espèce

- $ Y_0 $ la fonction de Bessel de 2ème espèce

- $ A $ et $ B $ constantes

La solution $ \phi $ de l'équation est de la forme : $ \phi(z) = AJ_0(2\alpha \sqrt{z}) + BY_0(2\alpha \sqrt{z}) $

A l'aide des condititions aux limites :

- $ \phi(0) = 1 $

- $ \phi_x{x=L} = ik\phi(x=L) $

On obtient alors:

- $ AJ_0(2\alpha \sqrt{H_1}) + BY_0(2\alpha \sqrt{H_0}) = 1 $

- $ AJ_1(2\alpha \sqrt{H_1})+ BY_1(2\alpha \sqrt{H_0}) = AJ_0(2\alpha \sqrt{H_1}) + BY_0(2 \alpha \sqrt{H_1}) $

Puis :

- $ A = \frac{Y_1^2 -Y_0^2}{J_0(Y_1^2 -Y_0^2)-Y_0(J_1^2 -iJ_0^2)} $

- $ \displaystyle B = \frac{iJ_0^2 -J_1^2}{J_0(Y_1^2 -Y_0^2)-Y_0(J_1^2 -iJ_0^2)} $

La solution s'écrit donc :

$ \displaystyle \phi(z) = \frac{Y_1^2- Y_0^2}{J_0(Y_1^2- Y_0^2)-Y_0(J_1^2 -iJ_0^2)} J_0(2\alpha \sqrt{z}) + \frac{iJ_0^2 -J_1^2}{J_0(Y_1^2 -Y_0^2)-Y_0(J_1^2 -iJ_0^2)}Y_0(2\alpha \sqrt{z}) $

Solution par méthode d'homotopie

On utilise l'équation: $ H(x)\phi_xx(x) + H'(x)\phi_x(x) + k^2H(x)\phi(x) = 0 $ avec $ H(x) = H_0 -sx $

Avec la relation d'homotopie, on revient à résoudre le système:

- $ \phi_{0,xx} = 0 $

- $ \phi_{1,xx} = -k^2(H_0-sx)\phi_0+s\phi_{0,x}-(H_0-sx-1)\phi_{0,xx} $

  • Ordre 0

On avait $ \phi_{0,xx} = 0 $ Donc $ \phi_0 = Ax+B $ avec $ A $ et $ B $ constantes à déterminer avec les conditions aux limites

A l'aide des conditions aux limites :

- $ \phi_0{x=0} = 1 \leftrightarrow B = 1 $

- $ \frac{\partial \phi_0}{\partial x}(x=L) = ik(AL+B) \leftrightarrow A = \frac{ik}{1-ikL} $

On obtient finalement :

$ \phi_0(x) = \frac{ik}{1-ikL}x+1 $

Quatrième Cas

Ce cas ci traite d'une vague sphérique générée par une source sinusoïdale périodique. La source est appliquée autour d'un cercle de rayon $ r_0 $ centré dans un domaine circulaire de rayon R, avec sortie libre de l'onde en r = R. Dans ce cas-là, l'équation de Berkhoff se simplifie en équation de Helmholtz, et peux s'exprimer en coordonnées polaires. on a ainsi les conditions suivantes :

- $ \phi_{rr} + \frac{1}{r}\phi_r + k^2\phi = 0 $ (symétrie de révolution, donc indépendant de $ \theta $ )

- $ \phi(r=r_0) = 1 $

- $ \phi_r(r=R) = ik \phi(r=R) $

Dans le cas des coordonnées polaires, l'équation de Berkhoff se simplifie en équation de Helmhotz: $ \Delta \phi +k^2\phi = 0 $

Ce qui donne donc l'équation : $ \dfrac{\partial^2 \phi}{\partial r^2} + \dfrac{\partial \phi}{\partial r} + k^2\phi = 0 $


Solution analytique

On pose :

- $ J_0 $ fonction de Bessel de première espèce

- $ Y_0 $ fonction de Bessel de 2ème espèce

- $ A $ et $ B $ constantes

On a alors $ \phi $ sous forme d'équation de Bessel : $ \phi(r) = AJ_0(r) + BY_0(r) $

A l'aide des conditions aux limites :

- $ AJ_0(r_0) + BY_0(r_0) = 1 $

- $ AJ'_0(R) + BY'_0(R) = ik(AJ_0(R) + BY_0(R)) $

- $ \dfrac{\partial r^n J_n(r)}{\partial r} = r^n J_{n-1}(r) $

- $ J_{-1}(n) = -J_1(n) $

On en déduit :

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

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

On a donc finalement :

$ \phi(r) = \frac{ikY_0(R)+Y_1(R)}{ikY_0(R)J_0(r_0)+Y_1(R)J_0(r_0)-ikJ_0(R)-J_1(R)Y_0(r_0)}J_0(r) + (1-J_0(r_0)\frac{ikY_0(R)+Y_1(R)}{ikY_0(R)J_0(r_0)+Y_1(R)J_0(r_0)-ikJ_0(R)-J_1(R)Y_0(r_0)}) \frac{1}{Y_0(r_0)}Y_0(r) $

Solution par méthode d'homotopie

On procède de la même manière que pour la cas 1 en trouvant $ \phi $ aux différents ordres.

  • Ordre 0

On a pour ce cas p = 0 On a la relation : $ \phi_0,rr = 0 \Leftrightarrow \phi_0 = Ar+B $

A l'aide des conditions aux limites :

- $ \phi_0(r=r_0) = 1 \Leftrightarrow Ar_0+B = 1 $

- $ \phi_r(r=R) = ik \phi(r=R) \Leftrightarrow A = ik(AR+B) $

On a donc : $ A = \frac{ik}{1+ik(r_0-R)} $ et $ B = \frac{1-ikR}{1+ik(r_0-R)} $

On en déduit alors : $ \phi_0 = \frac{ik}{1+ik(r_0-R)} r + \frac{1-ikR}{1+ik(r_0-R)} $

  • Ordre 1

On a pour ce cas la relation :

$ \phi_{1,rr} + \frac{1}{r}\phi_{0,r} + k^2\phi_0 = 0 \leftrightarrow \phi_1(r) + \frac{ik}{1+ik(r_0-R)}r(\ln(r)-1) + \frac{ik}{6(1+ik(r_0-R))}r^3 + \frac{1-ikR}{2(1+ik(r_0-R))}r^2 + Cr + D = 0 $

A l'aide des conditions aux limites suivantes :

- $ \phi_1(r=r_0) = 1 $

- $ \phi_r(r=R) = ik \phi(r=R) $

On en déduit donc :

- $ C = \frac{1}{1+ik(r_0-R)}(ik(AR(\ln(R)-1)\frac{ik}{1+ik(r_0-R)}+\frac{ik}{6(1+ik(r_0-R))}R^3 + \frac{1-ikR}{2(1+ik(r_0-R))}R^2 - A\ln(R) - \frac{ik}{2(1+ik(r_0-R))}R^3-\frac{1-ikR}{2(1+ik(r_0-R))}R)-A\ln(R) - \frac{ik}{2(1+ik(r_0-R))}R^3-\frac{1-ikR}{2(1+ik(r_0-R))}R) $

- $ D = \frac{1}{1+ik(r_0-R)}(A\ln(R) + \frac{ik}{2(1+ik(r_0-R))}R^3+\frac{1-ikR}{2(1+ik(r_0-R))}R - ik(AR(\ln(R)-1)\frac{ik}{1+ik(r_0-R)}+\frac{ik}{6(1+ik(r_0-R))}R^3 + \frac{1-ikR}{2(1+ik(r_0-R))}R^2 - Ar_0(\ln(r_0)-1)\frac{ik}{1+ik(r_0-R)}-\frac{ik}{6(1+ik(r_0-R))}r_0^3 - \frac{1-ikR}{2(1+ik(r_0-R))}r_0^2))r_0 -Ar_0(\ln(r_0)-1)\frac{ik}{1+ik(r_0-R)}-\frac{ik}{6(1+ik(r_0-R))}r_0^3 - \frac{1-ikR}{2(1+ik(r_0-R))} $


Ordre 2 et supérieurs

Dans le but de trouver $ \phi $ aux ordres 2 et supérieurs, on utilisera le logiciel WSMAXIMA. On pourra donc, par principe de superposition des solutions analytiques et de celle par méthode d'homotopie, vérifier les solutions obtenues.

Outils personnels