S'abonner à un flux RSS
 

ABET / CASTELLAN / KERGOAT : Différence entre versions

De Wikibardig
(Page créée avec « Bonjour, Vous voici rendu dans votre espace. Nous vous recommandons de prendre en mains l'outil WIKHYDRO assez rapidement de manière à pouvoir commencer à écrire en... »)
 
Ligne 1 : Ligne 1 :
Bonjour,
+
== IMPACT DU CHANGEMENT CLIMATIQUE SUR LES COTES ET DANS LES ESTUAIRES ==
  
Vous voici rendu dans votre espace.
+
=== Contexte ===
  
Nous vous recommandons de prendre en mains l'outil WIKHYDRO assez rapidement de manière à pouvoir commencer à écrire en LATEX, entrer des images et des vidéos.
 
  
Cette page fait partie intégrante du site du ministère de l'écologie. Elle est donc visible par tout internaute.Prenez-donc soin d'elle et faites en sorte qu'elle soit agréable à lire.
+
”Environ 1 milliard de personnes pourraient vivre d’ici 2050 dans des zones
 +
côtières menacées par la montée des eaux et les épisodes de submersions marines
 +
lors des tempêtes” alertent les experts de l’ONU.
  
Vous trouverez un [[Tutoriel]] et un [https://fr.wikipedia.org/wiki/Aide:Formules_TeX didacticiel LATEX]
+
C’est dans ce climat d’incertitude et d’insécurité qu’évoluent les populations
 +
cotières, c’est à dire plus de 60 pourcents de la population mondiale. C’est
 +
pourquoi des outils mathématiques sont nécessaires pour prévoir l’impact des
 +
houles sur le littoral.
  
Ce n'est pas une Mission Impossible, mais vous pouvez effacer cette introduction après lecture
 
  
Jean-Michel Tanguy
+
=== Etude des différents cas ===
 +
 
 +
 
 +
==== Cas n°1 ====
 +
 
 +
 
 +
On se place dans le cas d’un canal monodimensionnel plat de longueur L avec
 +
une fréquence d’entrée au niveau de l’aval φ = 1 et une sortie libre au niveau
 +
de l’amont avec φx = ikφ<\math> qui correspond à la condition de Robin
 +
 
 +
 
 +
===== ▷ Résolution Analytique ◁ =====
 +
 
 +
 
 +
On souhaite résoudre l’équation du modèle de Berkhoff ∇.(CCg ∇φ) + k2φ = 0
 +
 
 +
On pose une solution du potentiel φ = Ae^ikx + Be−ikx et on résoud l’équation
 +
à l’aide des conditions aux limites
 +
 
 +
 
 +
- Condition 1 :
 +
 
 +
φ(x = 0) = 1
 +
 
 +
φ(x = 0) = A + B = 1 ⇒  A = 1 − B
 +
 
 +
 
 +
- Condition 2 :
 +
 
 +
φx(L) = ikφ(L)
 +
 
 +
φx(L) = ikAeikL − ikBe−ikL = ik(AeikL + Be−ikL)
 +
 
 +
 
 +
Il vient:
 +
2ikBe−ikL ⇒ B = 0 ⇒  A = 1
 +
 
 +
 
 +
Finalement
 +
φ(x) = eikx = cos(kx) + isin(kx)
 +
 
 +
Ce qui nous donne:
 +
h(x, t) = Re(φ(x)×e−iωt = cos(kx)cos(ωt)−sin(kx)cos(ωt) =cos(kx − ωt)
 +
 
 +
 
 +
 
 +
===== ▷ Résolution par homotopie ◁ =====
 +
 
 +
 
 +
 
 +
On reprend les mêmes conditions initiales mais cette fois ci on part de la relation
 +
d’homotopie (1 − p)φxx + p(φxx+k2φ) = 0
 +
 
 +
 
 +
On introduit la décomposition de φ telle que:
 +
φ(x, p) = φ0(x) + pφ1(x) + p2φ3(x) + p3φ3(x) + ...
 +
 
 +
φxx(x) = φ0,xx(x) + pφ1,xx(x) + p2φ2,xx(x) + p3φ3,xx(x) + ...
 +
 
 +
 
 +
 
 +
L’équation générale d’homotopie avec la décomposition s’écrit donc
 +
 
 +
(1-p)(φ0,xx(x)+pφ1,xx(x)+p2φ2,xx(x)+p3φ3,xx(x)+...)+p(φ0,xx(x)+pφ1,xx(x)+p2φ2,xx(x) + p3φ3,xx(x) + ... + k2 + φ0(x) + pφ1(x) + p2φ3(x) + p3φ3(x) + ...) = 0
 +
 
 +
 
 +
* Ordre 0
 +
 
 +
 
 +
On reprend l’équation à l’ordre 0 telle que:
 +
(1 − p)φ0,xx + pφ0,xx = 0 ⇒ φ0,xx = 0 → φ0 = Ax + B
 +
 
 +
 
 +
Conditions aux limites
 +
 
 +
  - φ0:
 +
0 = 1 = B
 +
 
 +
  - φL:
 +
0,x = ikφL
 +
0 = ik(AL + 1) = A ⇒ A = \frac{ik}{1−ikL}
 +
On a donc
 +
φ0(x) =  \frac{ik}{1−ikL} x + 1
 +
 
 +
 
 +
 
 +
* Ordre 1
 +
 
 +
 
 +
On reprend l’équation à l’ordre 1 telle que:
 +
 
 +
(1 − p)(φ0,xx + pφ1,xx) + p(φ0,xx + pφ1,xx + k2φ0) = 0
 +
 
 +
On a vu à l’ordre 0 que φ0,xx = 0
 +
donc
 +
φ1,xx = −k2φ0 ⇒ φ1 = −k2 ∫ ∫ φ0dxdx + Ax + B
 +
 
 +
 
 +
 
 +
Conditions aux limites
 +
 
 +
- φ0
 +
2 = B
 +
 
 +
- φL
 +
2,x = ikφL
 +
2
 +
φL
 +
2,x = ik[−k2 ∫ ∫ (αx−k2(βx3 + 1
 +
2 x2))dxdx+AL] = ik[−k2(α 1
 +
6 L3 −k2(β 1
 +
20 L5 +
 +
1
 +
24 L4 + AL)]
 +
φL
 +
2,x = −k2[α 1
 +
2 L2 − k2(β 1
 +
4 L4 + 1
 +
6 L3)] + A
 +
A = ik
 +
1−ikL [−k2( 1
 +
6 αL3 − k2( 1
 +
20 βL5 + 1
 +
24 L4) − 1
 +
2 αL2 + k2(( 1
 +
4 βL4 + 1
 +
6 L3))]
 +
2
 +
φ2 = −k2[ 1
 +
6 αx3 − k2( 1
 +
20 βx5 + 1
 +
24 x4)] + Ax
 +
 
 +
 
 +
==== Cas n°2 ====
 +
 
 +
 
 +
 
 +
On se place dans le cas d’un domaine monodimensionnel plat de longueur L
 +
avec entrée par l’aval d’une onde de fréquence unitaire et une condition de flux
 +
aval φx = ik(2 − φ) et réflexion totale amont φx = 0.
 +
 
 +
 
 +
===== ▷ Résolution Analytique ◁ =====
 +
 
 +
 
 +
Comme dans le cas n°1:
 +
On obtient
 +
φ = Aeikx + Be−ikx
 +
et on résoud l’équation à l’aide des conditions aux limites:
 +
 
 +
 
 +
- En aval:
 +
φ(x = 0) = ik(2 − φ)
 +
ce qui nous donne
 +
φx(0) = ikA − ikB = ik(2 − A − B) ⇒ A = 1
 +
 
 +
 
 +
- En amont:
 +
φx(x = L) = 0
 +
ce qui nous donne
 +
φx(L) = ikAeikL − ikBe−ikL = 0 ⇒ B = e2ikL
 +
 
 +
 
 +
Finalement
 +
φ(x) = eikL − eik(2L−x)
 +
Et
 +
h(x, t) = Re(φe−iωt) = −2sin(kL − ωt)sin(kL + kx)
 +
 
 +
 
 +
===== ▷ Résolution par homotopie ◁ =====
 +
 
 +
 
 +
* Ordre 0
 +
 
 +
 
 +
On a
 +
φ0(x) = Ax + B et φ0,x = A
 +
 
 +
Ce qui nous donne φ0
 +
0,x = A = ik(2 − B) et φL
 +
0,x = A = 0
 +
 
 +
Donc B = 2 ⇒ φ0 = 2

Version du 27 mars 2023 à 16:24

Sommaire

IMPACT DU CHANGEMENT CLIMATIQUE SUR LES COTES ET DANS LES ESTUAIRES

Contexte

”Environ 1 milliard de personnes pourraient vivre d’ici 2050 dans des zones côtières menacées par la montée des eaux et les épisodes de submersions marines lors des tempêtes” alertent les experts de l’ONU.

C’est dans ce climat d’incertitude et d’insécurité qu’évoluent les populations cotières, c’est à dire plus de 60 pourcents de la population mondiale. C’est pourquoi des outils mathématiques sont nécessaires pour prévoir l’impact des houles sur le littoral.


Etude des différents cas

Cas n°1

On se place dans le cas d’un canal monodimensionnel plat de longueur L avec une fréquence d’entrée au niveau de l’aval φ = 1 et une sortie libre au niveau de l’amont avec φx = ikφ<\math> qui correspond à la condition de Robin


▷ Résolution Analytique ◁

On souhaite résoudre l’équation du modèle de Berkhoff ∇.(CCg ∇φ) + k2φ = 0

On pose une solution du potentiel φ = Ae^ikx + Be−ikx et on résoud l’équation à l’aide des conditions aux limites


- Condition 1 :

φ(x = 0) = 1
φ(x = 0) = A + B = 1 ⇒  A = 1 − B


- Condition 2 :

φx(L) = ikφ(L)
φx(L) = ikAeikL − ikBe−ikL = ik(AeikL + Be−ikL)


Il vient:

2ikBe−ikL ⇒ B = 0 ⇒  A = 1


Finalement

φ(x) = eikx = cos(kx) + isin(kx)

Ce qui nous donne:

h(x, t) = Re(φ(x)×e−iωt = cos(kx)cos(ωt)−sin(kx)cos(ωt) =cos(kx − ωt)


▷ Résolution par homotopie ◁

On reprend les mêmes conditions initiales mais cette fois ci on part de la relation d’homotopie (1 − p)φxx + p(φxx+k2φ) = 0


On introduit la décomposition de φ telle que:

φ(x, p) = φ0(x) + pφ1(x) + p2φ3(x) + p3φ3(x) + ...
φxx(x) = φ0,xx(x) + pφ1,xx(x) + p2φ2,xx(x) + p3φ3,xx(x) + ...


L’équation générale d’homotopie avec la décomposition s’écrit donc

(1-p)(φ0,xx(x)+pφ1,xx(x)+p2φ2,xx(x)+p3φ3,xx(x)+...)+p(φ0,xx(x)+pφ1,xx(x)+p2φ2,xx(x) + p3φ3,xx(x) + ... + k2 + φ0(x) + pφ1(x) + p2φ3(x) + p3φ3(x) + ...) = 0


  • Ordre 0


On reprend l’équation à l’ordre 0 telle que:

(1 − p)φ0,xx + pφ0,xx = 0 ⇒ φ0,xx = 0 → φ0 = Ax + B


Conditions aux limites

 - φ0:
0 = 1 = B
 - φL:
0,x = ikφL
0 = ik(AL + 1) = A ⇒ A = \frac{ik}{1−ikL}

On a donc

φ0(x) =  \frac{ik}{1−ikL} x + 1


  • Ordre 1


On reprend l’équation à l’ordre 1 telle que:

(1 − p)(φ0,xx + pφ1,xx) + p(φ0,xx + pφ1,xx + k2φ0) = 0

On a vu à l’ordre 0 que φ0,xx = 0 donc

φ1,xx = −k2φ0 ⇒ φ1 = −k2 ∫ ∫ φ0dxdx + Ax + B


Conditions aux limites

- φ0
2 = B
- φL

2,x = ikφL 2 φL 2,x = ik[−k2 ∫ ∫ (αx−k2(βx3 + 1 2 x2))dxdx+AL] = ik[−k2(α 1 6 L3 −k2(β 1 20 L5 + 1 24 L4 + AL)] φL 2,x = −k2[α 1 2 L2 − k2(β 1 4 L4 + 1 6 L3)] + A A = ik 1−ikL [−k2( 1 6 αL3 − k2( 1 20 βL5 + 1 24 L4) − 1 2 αL2 + k2(( 1 4 βL4 + 1 6 L3))] 2 φ2 = −k2[ 1 6 αx3 − k2( 1 20 βx5 + 1 24 x4)] + Ax


Cas n°2

On se place dans le cas d’un domaine monodimensionnel plat de longueur L avec entrée par l’aval d’une onde de fréquence unitaire et une condition de flux aval φx = ik(2 − φ) et réflexion totale amont φx = 0.


▷ Résolution Analytique ◁

Comme dans le cas n°1: On obtient

φ = Aeikx + Be−ikx

et on résoud l’équation à l’aide des conditions aux limites:


- En aval: 
φ(x = 0) = ik(2 − φ) 

ce qui nous donne

φx(0) = ikA − ikB = ik(2 − A − B) ⇒ A = 1


- En amont: 
φx(x = L) = 0

ce qui nous donne

φx(L) = ikAeikL − ikBe−ikL = 0 ⇒ B = e2ikL


Finalement 
φ(x) = eikL − eik(2L−x) 

Et

h(x, t) = Re(φe−iωt) = −2sin(kL − ωt)sin(kL + kx)


▷ Résolution par homotopie ◁
  • Ordre 0


On a

φ0(x) = Ax + B et φ0,x = A

Ce qui nous donne φ0

0,x = A = ik(2 − B) et φL
0,x = A = 0
Donc B = 2 ⇒ φ0 = 2
Outils personnels