S'abonner à un flux RSS
 

Jean-Michel Tanguy/SujetENTPE2022/AMIR - COTTET-GAYDON - WUILMOT : Différence entre versions

De Wikibardig
(Ordres supérieurs)
(3èmecas)
Ligne 195 : Ligne 195 :
 
=='''3<sup>ème</sup>cas'''==
 
=='''3<sup>ème</sup>cas'''==
  
Dans ce troisième cas, nous nous plaçons dans l'hypothèse d'un canal monodimensionnel avec une pente de fond constante <math> \displaystyle s = 1/200 </math>, avec entrée par l'aval d'une onde de fréquence unitaire et sortie libre en amont.
+
Dans ce troisième cas, nous nous plaçons dans l'hypothèse d'un canal monodimensionnel avec une pente de fond constante <math> \displaystyle s = \frac{1}{200} </math>, avec entrée par l'aval d'une onde de fréquence unitaire et sortie libre en amont.
  
 
''Conditions aux limites :'' <br>
 
''Conditions aux limites :'' <br>

Version du 8 juin 2022 à 15:44

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


Sommaire

Introduction

De nos jours le changement climatique est de plus en plus important et risque, si nous ne faisons rien, de bousculer nos vies. Un des premiers secteurs touchés par ce réchauffement climatique est le secteur côtier. Notamment, dans son rapport de février 2022, le GIEC alertait qu'une montée des eaux comme nous ne l'avions pas connu depuis 3000ans allait se produire à cause de ce dernier. En effet, d’après l’OCDE, le niveau des mers pourrait augmenter de plus d’un mètre dans les 80 prochaines années. Une telle augmentation engendrerait des déplacements massifs de population ainsi que de nombreux dommages sur nos infrastructures.

Erosion cotes.jpg

A travers ce projet, nous étudierons l'impact du changement climatique sur les côtes et dans les estuaires.

Modèle de Berkhoff

La propagation des vagues est un problème tridimensionnel compliqué. Pour simplifier ce problème, Berkhoff (1972) vient réduire ce problème a 2 dimensions en combinant en un seul modèle les effets de la réfraction et de la diffraction.
Ainsi, pour modéliser la houle nous allons utiliser le modèle de Berkhoff, ou l'équation de pente douce.

En bi-dimensionnel : $ \boxed{\displaystyle\color{red}\nabla.(CC_g\nabla \phi)+k^2CC_g\phi=0} $
avec :

  • $ \phi $ : le potentiel,
  • $ k $ : 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 $ : la célérité de l’onde,
  • $ C_g $ : la célérité de groupe des vagues.

Pour simplifier le problème, nous nous plaçons dans le domaine des houles longues (dont la longueur d'onde est supérieur a 300m). Ainsi : $ C=C_g=\sqrt{gH} $

L'équation de Berkhoff devient : $ \boxed{\displaystyle\color{red}\Delta \phi+k^2\phi=0} $

L'équation donnant l'évolution dans le temps de la hauteur de houle est : $ h(x,t)=\mathrm{Re}(\phi(x,t)) $ avec $ \phi(x,t)=\phi(x)e^{-iwt} $

Remarque : Dans le cas de la houle, l'entrée se situe en aval et la sortie en amont

Principe d'homotopie

L'homotopie est une méthode permettant de passer d'une solution initiale à une solution exacte grâce au paramètre p qui varie entre 0 à 1. L'idée est d'assurer une déformation continue entre la solution initiale, lorsque p=0, et la solution exacte du système d'équations à résoudre, lorsque p=1. L'intérêt de la méthode est qu'elle permet de partir d'une solution connue relativement simple et de converger vers une solution complète et ainsi atteindre la solution finale avec peu de termes. De plus elle utilise la résolution formelle et peut donc facilement être programmée à l'aide d'outils de calcul formel tel que WXMAXIMA.

L'objectif du projet est alors de résoudre l'EDP de Berkhoff par la méthode semi-analytique d'homotopie grâce à des conditions limites et des conditions initiales ad-hoc fonction du domaine étudié et du régime de houle du projet.

1ercas - Canal monodimensionnel plat avec sortie libre en amont

Dans un premier temps, nous étudions le cas d'un canal plat de longueur L.

Conditions aux limites :

  • Condition de Dirichlet : entrée par l'aval d'une onde de fréquence unitaire $ \phi(x=0) = 1 $
  • Condition de Robin : sortie libre par l'amont $ \phi_{x}(x=L) =ik\phi(x=L) $

Pour simplifier le problème, nous imposons ces 2 conditions en x=0 (car ceci est équivalent au cas où les 2 conditions sont imposées aux 2 extrémités du domaines (c'est à dire en x=0 et x=L))

Simplification de l'équation de Berkhoff :
Dans ce premier cas, l'équation de Berkhoff devient $ \boxed{\color{red}\displaystyle\frac{\partial^2 ϕ}{\partial x^2} +k^2ϕ = 0} $

Solution analytique

L'équation caractéristique associée à $ \frac{\partial^{2}\phi}{\partial x^{2}} + k^{2}\phi = 0 $ est : $ x^{2} + k^{2} = 0 $.
Ainsi $ \Delta = -4k^{2} < 0 $
d'où $ \phi(x)=Ae^{-ikx}+Be^{ikx} $

Grâce aux CL on détermine les constantes A et B :

  • $ \phi(x=0)=A+B=1 $
  • $ \dfrac{\partial ϕ}{\partial x}(x=L)=ik \phi \quad \Longrightarrow A=0 $ et $ B=1 $

De cette manière : $ \phi(x)=e^{ikx} \Longrightarrow \phi(x,t)=e^{i(kx-wt)} $

L'évolution de la hauteur de la houle devient alors : $ \boxed{\color{red}h(x,t)=\mathrm{Re}(\phi(x,t))=\cos(kx-wt)} $

Solution par homotopie

La relation d'homotopie s'écrit en choisissant la dérivée seconde comme fonction auxiliaire linéaire et en partant d'une solution initiale nulle:

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

En injectant la décomposition en série entière $ \phi(x,p)=\phi_0(x)+p\phi_1(x)+p^2\phi_2(x)+p^3\phi_3(x)+... $ et sa seconde dérivée:

$ \phi_{0,xx}(x)+p\phi_{1,xx}(x)+p^2\phi_{2,xx}(x)+p^3\phi_{3,xx}(x)+... $

Nous obtenons:

$ (1-p)(\phi_{0,xx}(x)+p\phi_{1,xx}(x)+p^2\phi_{2,xx}(x)+p^3\phi_{3,xx}(x)+...)+p[\phi_{0,xx}(x)+p\phi_{1,xx}(x)+p^2\phi_{2,xx}(x)+p^3\phi_{3,xx}(x)+...+k^2(\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

En p=0, la relation d'homotopie devient : $ \phi_{0,xx} = \frac{\partial^{2}\phi}{\partial x^{2}} = 0 \Longleftrightarrow \phi_{0} = Ax+B $

Grâce au CL on détermine les constantes A et B :

  • $ \phi_{0}(x=0)=1 \Longleftrightarrow B=1 $
  • $ \dfrac{\partial \phi_{0}}{\partial x}(x=L)=ik \phi_{0}(x=L) \Longleftrightarrow A=ik(AL+B) \Longleftrightarrow A=\frac{ik}{1-ikL} $

D'où $ \boxed{\color{red}\phi_{0} =\frac{ik}{1-ikL}x+1} $

Ordre 1

En p=1, la relation d'homotopie devient : $ \phi_{1,xx} = \frac{\partial^{2}\phi_{1}}{\partial x^{2}} + k^2\phi_{0}= 0 \Longleftrightarrow \phi_{1} = -k^2\iint\phi_{0}\mathrm{d}x\mathrm{d}x + A'x + B' \Longleftrightarrow \phi_{1} = -k^2(\frac{ikx^3}{6(1-ikL)}+\frac{x^2}{2}) + A'x + B' $

Grâce au CL on détermine les constantes A' et B' :

  • $ \phi_{1}(x=0)=0 \Longleftrightarrow B'=0 $
  • $ \dfrac{\partial \phi_{1}}{\partial x}(x=L)=ik \phi_{1}(x=L) \Longleftrightarrow -k^2(\frac{ikL^2}{2(1-ikL)}+L) + A' = ik(-k^2(\frac{ikL^3}{6(1-ikL)}+\frac{L^2}{2}) + A'L) \Longleftrightarrow A' = -\frac{k^2L(k^2L^2+3ikL-3)}{3(1-ikL)^2} $

D'où $ \boxed{\color{red}\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 supérieurs

On utilise le logiciel WXMAXIMA pour trouver $ \phi $ à tous les ordres.

avec $ k=\frac{1}{100} $ (nombre d'onde en m-1), $ H=40 $ (profondeur en m), $ c=\sqrt{gH} $ (célérité de l'onde en m/s), $ \lambda=\frac{2\pi}{k} $ (longueur d'onde en m), $ L=2\lambda $ (longueur du domaine en m).

Puis on superpose la solution analytique et la solution par homotopie

Cas 1.gif

On constate que la solution par homotopie converge vers la solution analytique

Étude de la sensibilité de la solution en fonction du nombre d'onde

Dans ce premier cas, on constate que plus le produit $ kL $ est petit, plus la solution par homotopie converge rapidement (exemple avec $ kL=0,25 $, $ kL=1 $ et $ kL=2 $)

Cas 1 kL=0,25.png Cas 1 kL=1.png Cas 1 kL=2.png

2èmecas - Canal monodimensionnel plat avec réflexion totale en amont

Dans ce deuxième cas, nous nous plaçons toujours dans l'hypothèse d'un canal plat de longueur L mais avec des conditions aux limites différentes.

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 $

Simplification de l'équation de Berkhoff :
L'équation de Berkhoff se simplifie de la même façon que dans le 1ercas $ \boxed{\color{red}\displaystyle\frac{\partial^2 ϕ}{\partial x^2} +k^2ϕ = 0} $

Solution analytique

L'équation caractéristique associée à $ \frac{\partial^{2}\phi}{\partial x^{2}} + k^{2}\phi = 0 $ est : $ x^{2} + k^{2} = 0 $.
Ainsi $ \Delta = -4k^{2} < 0 $
d'où $ \phi(x)=Ae^{-ikx}+Be^{ikx} $

Grâce aux CL on détermine les constantes A et B :

  • $ \dfrac{\partial ϕ}{\partial x}(x=0)=ik(2- \phi(x=0)) \Longleftrightarrow -ikA + ikB = ik(2-(A+B)) \Longleftrightarrow B=1 $
  • $ \dfrac{\partial ϕ}{\partial x}(x=L)=0 \Longleftrightarrow -ikAe^{-ikL}+ikBe^{ikL} = 0 \Longleftrightarrow A=e^{2ikL} $

De cette manière : $ \phi(x)=e^{ik(2L-x)}+e^{ikx} \Longrightarrow \phi(x,t)=e^{i(k(2L-x)-wt)}+e^{i(kx-wt)} $

L'évolution de la hauteur de la houle devient alors : $ \boxed{\color{red}h(x,t)=\mathrm{Re}(\phi(x,t))=cos(k(2L-x)-wt)+cos(kx-wt)} $

Solution par homotopie

Ordre 0

En p=0, la relation d'homotopie devient : $ \phi_{0,xx} = \frac{\partial^{2}\phi}{\partial x^{2}} = 0 \Longleftrightarrow \phi = Ax+B $

Grâce au CL on détermine les constantes A et B :

  • $ \dfrac{\partial \phi_{0}}{\partial x}(x=L)=0 \Longleftrightarrow A=0 $
  • $ \dfrac{\partial \phi_{0}}{\partial x}(x=0)=ik(2- \phi_{0}(x=0)) \Longleftrightarrow 0=ik(2-B) \Longleftrightarrow B=2 $

D'où $ \boxed{\color{red}\phi_{0} = 2} $

Ordre 1

En p=1, la relation d'homotopie devient : $ \phi_{1,xx} = \frac{\partial^{2}\phi_{1}}{\partial x^{2}} + k^2\phi_{0}= 0 \Longleftrightarrow \phi_{1} = -k^2x^2 + A'x + B' $

Grâce au CL on détermine les constantes A' et B' :

  • $ \dfrac{\partial \phi_{1}}{\partial x}(x=L)=0 \Longleftrightarrow A'=2k^2L $
  • $ \dfrac{\partial \phi_{1}}{\partial x}(x=0)=ik\phi_{1}(x=0) \Longleftrightarrow 2k^2L=ikB' \Longleftrightarrow B'=2ikL $

D'où $ \boxed{\color{red}\phi_{1} = -k^2x^2 + 2k^2Lx + 2ikL} $

Ordres supérieurs

On utilise le logiciel WXMAXIMA pour trouver $ \phi $ à tous les ordres

Puis on superpose la solution analytique et la solution par homotopie

Cas 2 kL=0,001.gif

On constate que la solution par homotopie converge vers la solution analytique

Étude de la sensibilité de la solution en fonction du nombre d'onde

On remarque cependant que selon la valeur du produit $ kL $ (notamment $ kL \geqslant 1 $) la solution ne converge pas

Exemple pour $ kL = 1 $ :

Cas 2 kL=1.gif

Analyse des résultats, des limites et de l'interêt de la méthode d'homotopie dans les cas 1 et 2

Grâce à la méthode d'homotopie nous pouvons déterminer les solutions du modèle à partir de solutions simples et grâce au logiciel WXMAXIMA nous pouvons tracer l'évolution de la houle en fonction du temps.

Ce modèle nous permet d'étudier un cas simple de modélisation de la houle. Tout d'abord, nous faisons l'hypothèse d'une onde monodimensionnelle alors qu'en réalité l'onde est tridimensionnelle. Ensuite, nous n'avons pas pris en compte plusieurs paramètres tel que l'irrégularité des fonds océaniques et des littoraux.

3èmecas

Dans ce troisième cas, nous nous plaçons dans l'hypothèse d'un canal monodimensionnel avec une pente de fond constante $ \displaystyle s = \frac{1}{200} $, avec entrée par l'aval d'une onde de fréquence unitaire et sortie libre en amont.

Conditions aux limites :

  • En aval : $ \displaystyle \phi = 1 $ (Condition de Dirichlet)
  • En amont : $ \displaystyle \phi_{x} = ik\phi $ (Condition de Robin)

La résolution est plus complexe car $ \displaystyle H(x) $ n'est plus égal à $ \displaystyle Ho $ mais à $ \displaystyle H(x) = H_{0} - sx $

En reprenant le modèle de Berkhoff : $ \displaystyle ∇.(CCg∇ϕ)+k^2CCgϕ=0 $ avec $ \displaystyle CCg = gH(x) $, on obtient : $ \displaystyle \boxed{\color{red}\displaystyle H(x)\frac{\partial^2 \phi}{\partial x^2} + (-s)\frac{\partial \phi}{\partial x} + k_{0}^2H_{0} \phi = 0} $

Solution analytique

L'équation de Berkhoff se simplifie en prenant pour expression de la profondeur du fond : $ \displaystyle H(x) = H_{0} - sx $ et en prenant k (nombre d'onde) non constant et égal à : $ \displaystyle k(x) = k_{0}\sqrt{\frac{H_{0}}{H_{0}-sx}} $

On cherche donc à résoudre l'équation : $ \displaystyle (H_{0} - sx)\frac{\partial^2 \phi}{\partial x^2} + (-s)\frac{\partial \phi}{\partial x} + k_{0}^2H_{0} \phi = 0 $

Pour se faire, on cherche à mettre cette équation sous la forme d'une équation de type Bessel : $ \displaystyle z\frac{\partial^2 \phi}{\partial z^2} + \frac{\partial \phi}{\partial z} + \alpha^2 \phi = 0 $

En effectuant le changement de variable $ \displaystyle z(x) = H_{0} - sx $, on trouve : $ \displaystyle \boxed{\color{red}\displaystyle z(x)\frac{\partial^2 \phi}{\partial z^2} + \frac{\partial \phi}{\partial z} + \frac{k_{0}^2H_{0}}{s^2} \phi = 0} $

Donc : $ \displaystyle \alpha^2 = \frac{k_{0}^2H_{0}}{s^2} $

La solution générale de cette équation est : $ \displaystyle \phi(z) = C_{1} * J_{0}(2\alpha\sqrt{z}) + C_{2} * Y_{0}(2\alpha\sqrt{z}) $

Avec :

  • $ \displaystyle J_{0} $ : fonction de Bessel de première espèce
  • $ \displaystyle Y_{0} $ : fonction de Bessel de deuxième espèce
  • $ \displaystyle C_{1} $ et $ \displaystyle C_{2} $ : constantes à déterminer


On détermine les constantes $ \displaystyle C_{1} $ et $ \displaystyle C_{2} $ à l’aide des conditions aux limites :

  • $ \displaystyle \phi(z = H_{0}) = 1 $
  • $ \displaystyle -s \phi(z = H_{0} - sL) = ik\phi(z = H_{0} – sL) $

Et en utilisant le fait que :

  • $ \displaystyle \frac{\partial z^nJ_n(z)}{\partial z} = z^nJ_{n-1}(z) $
  • $ \displaystyle J_{-1}(z) = -J_1(z) $


On résout le système et on trouve :

  • $ \displaystyle C_{1} = \frac{Y_{1}^L – i Y_{0}^L}{J_{0}^0(Y_{1}^L – i Y_{0}^L) – Y_{0}^0(J_{1}^L – i J_{0}^L)} $
  • $ \displaystyle C_{2} = - \frac{ J_{1}^L – i J_{0}^L}{ J_{0}^0(Y_{1}^L – i Y_{0}^L) – Y_{0}^0(J_{1}^L – i J_{0}^L)} $

Avec :

  • $ \displaystyle J_{0}^0 = J_{0}(2\alpha\sqrt{H_{0}}) $
  • $ \displaystyle Y_{0}^0 = Y_{0}(2\alpha\sqrt{H_{0}}) $
  • $ \displaystyle J_{0}^L = J_{0}(2\alpha\sqrt{H_{0} - sL}) $
  • $ \displaystyle Y_{0}^L = Y_{0}(2\alpha\sqrt{H_{0} - sL}) $
  • $ \displaystyle J_{1}^L = J_{1}(2\alpha\sqrt{H_{0} - sL}) $
  • $ \displaystyle Y_{1}^L = Y_{1}(2\alpha\sqrt{H_{0} - sL}) $


D’où : $ \boxed{\color{red}\displaystyle \phi(z) = \frac{Y_{1}^L – i Y_{0}^L}{J_{0}^0(Y_{1}^L – i Y_{0}^L) – Y_{0}^0(J_{1}^L – i J_{0}^L)} * J_{0}(2\alpha\sqrt{z}) - \frac{ J_{1}^L – i J_{0}^L}{ J_{0}^0(Y_{1}^L – i Y_{0}^L) – Y_{0}^0(J_{1}^L – i J_{0}^L)} * Y_{0}(2\alpha\sqrt{z})} $

Solution par homotopie

Pour ce cas, la relation d’homotopie est la suivante :

$ \displaystyle (1 – p) \phi_{xx} + p ((1 - \epsilon x) \phi_{xx} - \epsilon \phi_{x} + ko^2 \phi) = 0 $

A l’aide de la décomposition en série entière de $ \displaystyle \phi $ et de ses dérivées, on obtient :

$ \displaystyle (1 – p) * (\phi_{0,xx}(x) + p\phi_{1,xx}(x) + p^2\phi_{2,xx}(x) + p^3\phi_{3,xx}(x) + …) + p * ((1 - \epsilon x) * (\phi_{0,xx} + p\phi_{1,xx}(x) + p^2\phi_{2,xx}(x) + p^3\phi_{3,xx}(x) + …) - \epsilon * (\phi_{0,x} + p\phi_{1,x}(x) + p^2\phi_{2,x}(x) + p^3\phi_{3,x}(x) + …) + ko^2 * (\phi_{0} + p\phi_{1}(x) + p^2\phi_{2}(x) + p^3\phi_{3}(x) + …)= 0 $


Ordre 0

A l’ordre 0 (correspondant à une puissance de p nulle), on a la relation :

$ \displaystyle \phi_{0,xx} = 0 \Longleftrightarrow \phi_{0} = A * x + B $


A l’aide des conditions aux limites, on détermine les constantes A et B :

  • $ \displaystyle \phi_{0}(x = 0) = 1 \Longleftrightarrow B = 1 $
  • $ \displaystyle \phi_{0,x}(x = L) = ik(x) \phi_{0}(x = L) \Longleftrightarrow A = \frac{ik(x)}{1 – ik(x)L} $

D’où : $ \boxed{\color{red}\displaystyle \phi_{0}(x) = 1 + \frac{ik(x)}{1 – ik(x)L}*x} $

Ordre 1

A l’ordre 1, on a la relation :

$ \displaystyle \phi_{1,xx} - \phi_{0,xx} + (1 - \epsilon x)\phi_{0,xx} - \epsilon\phi_{0,x} + k_{0}^2 \phi_{0} = 0 $

Cette expression se simplifie : $ \displaystyle \phi_{1,xx} - \epsilon\phi_{0,x} + k_{0}^2 \phi_{0} = 0 $

Comme on connait $ \displaystyle \phi_{0} $, en intégrant deux fois cette expression, on trouve $ \displaystyle \phi_{1} $ : $ \displaystyle \phi_1(x) = \int \int (\epsilon \phi_{0,x}-k_0^2\phi_0) \,dx\,dx +A_1x+B_1 $

A l’aide des conditions aux limites, $ \displaystyle \phi_{1}(x = 0) = 0 $ et $ \displaystyle \phi_{1,x}(x = L) = ik(x) \phi_{1}(x = L) $, on trouve les constantes grâce à la programmation sur le logiciel wxmaxima.

Ordres supérieurs

Pour n $ \geq 1 $, on a la relation :

$ \displaystyle \phi_{n,xx} - \phi_{n-1,xx} + (1 - \epsilon x)\phi_{n-1,xx} - \epsilon \phi_{n-1,x} + k_0^2 \phi_{n-1} = 0 $

Donc, $ \displaystyle \phi_n(x) = \int \int (\epsilon x \phi_{n-1,xx} + \epsilon\phi_{n-1,x} - k_0^2 \phi_0) \,dx\,dx + A_n x + B_n $

On peut résoudre cette équation par récurrence avec les conditions initiales :

  • $ \displaystyle \phi_{n}(x = 0) = 0 $
  • $ \displaystyle \phi_{n,x}(x = L) = ik(x) \phi_{n}(x = L) $

On utilise le logiciel WXMAXIMA pour trouver $ \displaystyle \phi $ à tous les ordres puis on superpose la solution analytique et la solution par homotopie

Webp.net-gifmaker.gif

Étude de la sensibilité de la solution en fonction du nombre d'onde

Analyse des résultats, des limites et de l'interêt de la méthode d'homotopie dans le cas 3

4èmecas

Dans ce 4ème cas, nous allons nous intéresser à une vague sphérique générée par une source périodique sinusoïdale.

Nous traitons ici de l'évolution de la surface libre 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 $.

Conditions aux limites:

  • $ \phi(r=r_0)=1 $
  • $ \frac{\partial \phi}{\partial r}(r=R)=ik\phi(r=R) $

Simplification de l'équation de Berkhoff :

L'équation de Berkhoff se simplifie alors en équation de Helmholtz et s'exprime en coordonnées polaires avec les conditions suivantes $ \displaystyle\boxed{\color{red} \Delta \phi + k^2\phi=0 } $

En coordonnées polaires, la relation ci-dessus s'écrit de manière simplifiée étant donné que le problème est caractérisé par une symétrie de révolution, donc est indépendant de $ \theta $.

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

Solution analytique

L'équation $ \frac{\partial^2 \phi}{\partial r^2} +\frac{\partial \phi}{\partial r} + k^2\phi=0 $ est une équation de Bessel. La solution générale de cette équation est donc de la forme :

$ \phi(r)=AJ_0(r)+BY_0(r) $

  • $ J_0 $ : fonction de Bessel de 1ère espèce
  • $ Y_0 $ : fonction de Bessel de 2ème espèce
  • A et B : constante

Grâce aux conditions limites, on détermine les constantes A et B.

  • $ AJ_0(r_0)+BY_0(r_0)=1 $
  • $ AJ_0'(R)+BY_0'(R)=ik\ (AJ_0(R)+BY_0(R)) $


En utilisant le fait que :

  • $ \displaystyle\frac{\partial r^nJ_n(r)}{\partial r} = r^nJ_{n-1}(r) $
  • $ J_{-1}(r) = -J_1(r) $

On trouve :

  • $ A = \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)} $
  • $ B = \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 par homotopie

Pour ce cas, la relation d'homotopie est la suivante : $ (1-p)\phi_{rr}+p(\phi_{rr}+\frac{1}{r}.\phi_{r}+k^2\phi)=0 $

Comme précédemment, on utilise la décomposition en série entière $ \phi(r,p)=\phi_0(r)+p\phi_1(r)+p^2\phi_2(r)+p^3\phi_3(r)+... $ et sa seconde dérivée: $ \phi_{0,rr}(r)+p\phi_{1,rr}(r)+p^2\phi_{2,rr}(r)+p^3\phi_{3,rr}(r)+... $

Nous obtenons:

$ (1-p)(\phi_{0,rr}(r)+p\phi_{1,rr}(r)+p^2\phi_{2,rr}(r)+p^3\phi_{3,rr}(r)+...)+p[\phi_{0,rr}(r)+p\phi_{1,rr}(r)+p^2\phi_{2,rr}(r)+p^3\phi_{3,rr}(r)+...+k^2(\phi_0(r)+p\phi_0(r)+p^2\phi_2(r)+p^3\phi_3(r)+...)]=0 $

Ordre 0

A l'ordre 0, on obtient la relation $ \phi_{0,rr}(r)=0 \Longleftrightarrow \phi_0(r)=Ar + B $

Grâce aux CL on détermine les constantes A et B :

  • $ Ar_0 + B = 1 \Longleftrightarrow A = \displaystyle \frac{ik}{1 + ik(r_0 - R)} $
  • $ A = ik(AR + B) \Longleftrightarrow B = \displaystyle \frac{1 - ikR}{1 + ik(r_0 - R)} $


D'où $ \boxed{\color{red}\displaystyle \phi_0(r) = \frac{ik}{1 + ik(r_0 - R)}r + \frac{1 - ikR}{1 + ik(r_0 - R)} } $

Ordre 1

A l'ordre 1, on trouve la relation $ \displaystyle \phi_{1,rr}(r) + \frac{1}{r}\phi_{0,r}(r) + k^2\phi_0(r) =0 \Longleftrightarrow \displaystyle\phi_1(r) + A r(ln(r)-1) + \frac{Ak^2}{6}r^3 + \frac{Bk^2}{2}r^2 +Cr + D = 0 $

avec :

  • $ A = \displaystyle \frac{ik}{1 + ik(r_0 - R)} $
  • $ B = \displaystyle \frac{1 - ikR}{1 + ik(r_0 - R)} $

En utilisant les conditions limites, on trouve que : $ \boxed{\displaystyle\color{red}\phi_1(r) + A r(ln(r)-1) + \frac{Ak^2}{6}r^3 + \frac{Bk^2}{2}r^2 +\frac{ik(\delta - \alpha) - \beta}{1 + ik(r_0 - R)}r + \frac{\beta - ik(\delta - \alpha)}{1 + ik(r_0 - R)}r_0 - \alpha = 0 } $

Sachant que :

  • $ \alpha = \displaystyle r_0(ln(r_0)-1)A + \frac{Ak^2}{6}r_0^3 + \frac{Bk^2}{2}r_0^2 $
  • $ \beta = \displaystyle Aln(R) + \frac{Ak^2}{2}R^2 + Bk^2R $
  • $ \delta = \displaystyle R(ln(R)-1)A + \frac{Ak^2}{6}R^3 + \frac{Bk^2}{2}R^2 $

Ordres supérieurs

Pour n $ \geq 1 $, on trouve la relation de récurrence suivante : $ \boxed{\displaystyle\color{red}\phi_{n,rr}(r)+\frac{1}{r}\phi_{n-1,r}(r)+k^2\phi_{n-1}(r) =0 } $.

On utilise ensuite le logiciel MAXIMA pour résoudre et trouver$ \phi $à un ordre quelconque.


Application temporelle

Dans cette partie, nous allons réaliser une application temporelle avec une fréquence de $ \omega=\sqrt{gk} $ car nous sommes en profondeur infinie.

Outils personnels