03/02/2010

Modélisation des
réservoirs souterrains
A. Dassargues
[fichier ModRes4 2010]

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Méthode des éléments finis pour les

écoulements en milieu souterrain
ECOULEMENT PERMANENT
Principes généraux
Minimum d’une fonctionnelle
Evaluation de la solution
Implémentation pratique
Généralisation

ECOULEMENT TRANSITOIRE

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

1
 03/02/2010

Ecoulement permanent
Principes généraux

„ écoulement 2D en milieu poreux

„ flux d’infiltration
K'
e'
„ flux de drainance h h1
K'
L= (h1 − h)
e'
Equation de l’écoulement
en milieu poreux
∂h ∂ ⎛ ∂h ⎞ ∂ ⎛ ∂h ⎞
S. = I + ⎜ T . ⎟ + ⎜⎜ T . ⎟⎟ + L
∂t ∂x ⎝ ∂x ⎠ ∂y ⎝ ∂y ⎠
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Principes généraux (2)

„ éléments triangulaires ou rectangulaires

„ les inconnues (hauteurs piézométriques) ne sont
calculées qu’aux nœuds: la valeur en tout point est
calculée en utilisant les fonctions de base ou les
fonctions d’interpolation
„ la méthode convient particulièrement pour les
domaines très hétérogènes et anisotropes et aux
frontières irrégulières

Le fonctionnement de la méthode des éléments finis

est démontré ci-dessous en utilisant l’ “approximation
par minimisation d’une fonctionnelle”.

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

2
 03/02/2010

Principes généraux (3)

Les étapes suivantes sont adoptées :

„ l’hypothèse que la solution est unique

„ la recherche d’une fonctionnelle qui sera minimum
pour la solution
„ la preuve que la solution trouvée est bien la meilleure
possible
„ l’utilisation des fonctions d’interpolation pour trouver
la valeur locale de la fonctionnelle
„ l’expression aux différents nœuds du minimum de la
fonctionnelle par rapport à la variable principale
conduira automatiquement à un système d’équations
algébriques

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Principes généraux (4)

Pour l’écoulement permanent:

∂ ⎛ ∂h ⎞ ∂ ⎛ ∂h ⎞ (h − h1 )
⎜ T . ⎟ + ⎜⎜ T . ⎟⎟ + I − =0
∂x ⎝ ∂x ⎠ ∂y ⎝ ∂y ⎠ c
c = e' / K '
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

3
 03/02/2010

Principes généraux (5)

Cette équation doit être intégrée sur le domaine étudié en
tenant compte des conditions aux frontières:
„ hauteur piézométrique imposée sur une partie S1 de la
frontière h= f
„ flux nul imposé sur une partie S2 de la frontière
S2

S1
∂h ∂n = 0
h=f(x,y,t) waterdichte grenzen
frontières imperméables

S2 h=f(x,y,t)

S1
x
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Minimum d’une fonctionnelle

Pour notre problème d’écoulement 2D, nous pourrons
montrer que la fonctionnelle suivante peut être utilisée :

1 ⎡ ⎛ ∂h ⎞ h 2 − 2.h.h1 ⎤
2
⎛ ∂h ⎞
2

U = ∫∫ ⎢T .⎜ ⎟ + T .⎜⎜ ⎟⎟ − 2.I .h + ⎥ dxdy

2 R ⎢ ⎝ ∂x ⎠ ∂
⎝ ⎠y c ⎥⎦
⎣
jusqu’à présent, h(x,y) est une fonction quelconque excepté que :

h= f sur S1
∂h ∂n = 0 sur S2

Si U = U* pour h = h* U* doit être le minimum de U

h = h * +α . v( x, y )

avec v( x, y ) une fonction quelconque, mais v = 0 sur S1

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

4
 03/02/2010

Minimum d’une fonctionnelle (2)

… en entrant cette valeur dans la fonctionnelle :

⎡ ⎛ ∂h * ⎞ ⎛ ∂v ⎞ ⎛ ∂h * ⎞ ⎛ ∂v ⎞ h * .v − v.h1 ⎤
U = U * +α .∫∫ ⎢T .⎜ ⎟.⎜ ⎟ + T .⎜⎜ ⎟⎟.⎜⎜ ⎟⎟ − I .v + ⎥ dxdy
⎣ ⎝ ∂x ⎠ ⎝ ∂x ⎠ ⎝ ∂y ⎠ ⎝ ∂y ⎠ c ⎦
R

⎡ ⎛ ∂v ⎞
2 2
⎛ ∂v ⎞ v 2⎤
+ α 2 .∫∫ ⎢T .⎜ ⎟ + T .⎜⎜ ⎟⎟ + ⎥ dxdy
R
⎢⎣ ⎝ ∂x ⎠ ⎝ ∂y ⎠ c⎥
⎦
U = U * +α . A1 + α 2 . A2
… développement de A1 A2 >0
⎡∂ ⎛ ∂h * ⎞ ∂ ⎛ ∂h * ⎞⎤
A1 = ∫∫ ⎢ .⎜ T .v. ⎟ + .⎜⎜ T .v. ⎟⎥dxdy
∂y ⎟⎠⎦
R ∂x =0
⎣ ⎝ ∂x ⎠ ∂y ⎝
⎡ ∂ ⎛ ∂h * ⎞ ∂ ⎛ ∂h * ⎞ h * − h1 ⎤
− ∫∫ v.⎢ ⎜ T . ⎟ + ⎜⎜ T . ⎟⎟ + I − ⎥ dxdy
⎣ ∂x ⎝ ∂x ⎠ ∂y ⎝ ∂y ⎠ c ⎦
R
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Minimum d’une fonctionnelle (3)

⎛ ∂v x ∂v y ⎞
… par le théorème de Green : ∫∫ ⎜⎜⎝ ∂x +
∂y ⎟⎠
⎟dxdy = ∫ vn .dS
S

∂h * ∂h *
A1 = ∫ T .v. .dS + ∫ T .v. .dS
S1 ∂n S2 ∂n
∂h *
avec ∫ S1
T .v.
∂n
.dS = 0 car v=0 sur S1

∂h * ∂h ∂n = 0 sur S2
∫S2
T .v.
∂n
.dS = 0 car

A1 = 0 U > U*

… le minimum de U est trouvé pour h = h*(x,y)

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

5
 03/02/2010

Evaluation de la solution
… maintenant il faut prouver que h*(x,y) est bien la meilleure
des solutions (au vu du type d’approximation choisi dans
la méthode numérique)
… sur un exemple 2
∂ h
T. +I =0
I ∂x 2
y

h
l
h0
x

x = 0 → h = h0 x = l → ∂h ∂x = 0
I .l I 2
solution analytique h = h0 + .x − .x
T 2T Terug
Hydrogéologie et Géologie
naar eerste pagina
de l’Environnement

Evaluation de la solution (2)

Fonctions d’interpolation linéaire (entre les nœuds)
∂h
h = h0 + C.x =C
∂x

U=
1 l
2 ∫0
[
T .C 2
− 2 . I .h0 − 2 . I .C . x ]
dx =
1
2
T .C 2
.l − I .h0 .l − I .C .
l2
2
∂ 2U
… le minimum de U par rapport à C pour : ∂U = 0 >0
∂C ∂C 2
∂U I .l 2 I .l
= T .C.l − =0 C=
∂C 2 2T
∂ 2U … solution numérique trouvée :
= T .l > 0
∂C 2
I .l
h = h0 + .x
2T Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

6
 03/02/2010

Evaluation de la solution (3)

Fonctions d’interpolation linéaire (2) solution numérique :

I .l
h = h0 + .x
2T
solution analytique
analytische oplossing

rechte
droitelijn

l
h0
x

I .l I 2
solution analytique : h = h0 + .x − .x
T 2T
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Evaluation de la solution (4)

Fonctions d’interpolation parabolique
∂h
h = h0 + C.x + D.x 2 = C + 2.D.x
∂x
⎛ ∂h ⎞
2

⎜ ⎟ = C + 4.C.D.x + 4.D .x
2 2 2

⎝ ∂x ⎠

1⎡ ⎛ 2 4 ⎞ ⎛ 1 1 ⎞⎤
U= ⎢T .⎜ C .l + 2.C.D.l 2 + . D 2 . l 3 ⎟ − 2.I .⎜ h0 .l + .C.l 2 + . D. l 3 ⎟⎥
2⎣ ⎝ 3 ⎠ ⎝ 2 3 ⎠⎦
… le minimum de U par rapport à C et D :

∂U I .l 2 ∂ 2U
= T .C.l + T .D.l 2 − =0 = T .l > 0
∂C 2 ∂C 2
∂U 4 I .l 3 ∂ 2U 4 3
= T .C.l 2 + .T .D.l 3 − =0 = .T .l > 0
∂D 3 3 ∂D 2 3 Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

7
 03/02/2010

Evaluation de la solution (5)

Fonctions d’interpolation parabolique (2)
solution numérique :
I .l I 2
C=
I .l
D=−
I h = h0 + .x − .x
T 2.T T 2.T
solution analytique
analytische oplossing
= numerieke
= solution numérique
oplossing

l
h0
x

I .l I 2
solution analytique : h = h0 + .x − .x
T 2T Hydrogéologie et Géologie
Terug naar eerste pagina
de l’Environnement

Implémentation en pratique

… pour des éléments finis triangulaires

Pour notre problème (considéré sans flux L), la fonctionnelle

peut s’écrire (dans l’élément j ) :
1 ⎡ ⎛ ∂h ⎞ 2 ⎛ ∂h ⎞
2
⎤
2 ∫∫R j ⎢ ⎝ ∂x ⎠
Uj = ⎢T .⎜ ⎟ + T .⎜ ⎟
⎜ ∂y ⎟ − 2.I .h ⎥ dxdy
⎣ ⎝ ⎠ ⎥⎦
U j = U 1j + U 2j
1 ⎡ ⎛ ∂h ⎞ 2 ⎛ ∂h ⎞ ⎤
2

2 ∫∫R j ⎢ ⎝ ∂x ⎠
avec U = 1
⎢T .⎜ ⎟ + T .⎜⎜ ⎟⎟ ⎥ dxdy
⎝ ∂y ⎠ ⎥⎦
j
⎣

U 2j = − ∫∫ I .[Link]
Rj

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

8
 03/02/2010

Implémentation en pratique (2)

… dans chaque élément, h est supposé une fonction
linéaire des coordonnées (interpolation linéaire) :
h = p.x + q. y + r
… le champ de la piézométrie est décrit dans chaque
élément fini par un plan, le champ complet est approximé
par l’ensemble des plans des différents éléments
Chaque plan est défini par les 3 valeurs de h aux noeuds
vlak
plan
z
h3
h4
h1
h2

y
x
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Implémentation en pratique (3)

dans l’élément j , les deux termes de l’intégrale
peuvent être calculés séparément en entrant la
valeur de h :
1
[
U 1j = .∫∫ T .( p 2 + q 2 ) dxdy
2 Rj
]
p et q seront uniquement dépendants des coordonnées
et des valeurs de piézométrie aux 3 noeuds de l’élément

1 1
U 1j = .T j .( p 2 + q 2 ).∫∫ dxdy = .T j .( p 2 + q 2 ). A j
2 Rj 2
1 1
avec A j = . x1.( y2 − y3 ) + x2 .( y3 − y1 ) + x3 .( y1 − y2 ) = .∆
2 2

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

9
 03/02/2010

Implémentation en pratique (4)

U 2j = − I j .∫∫ [Link]
Rj

approximation du terme d’infiltration, en utilisant

une valeur moyenne de h sur l’élément

1 1 1 3
h= .∫∫ [Link] = .(h1 + h2 + h3 ) = .∑ hk
Aj R j 3 3 k =1
I j .Aj 3
U 2j = − .∑ hk
3 k =1

… donc, dans l’élément j :

1 I .A 3
U j = .T j .( p 2 + q 2 ). A j − j j .∑ hk
2 3 k =1
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Implémentation en pratique (5)

pour l’élément j , le plan est calculé sur base des 3
valeurs de h ( h1, h2, h3) aux trois noeuds :

⎧h1 = p.x1 + q. y1 + r
⎪
⎨h2 = p.x2 + q. y2 + r
⎪h = p.x + q. y + r ⎧b1 = y2 − y3
⎩ 3 3 3
⎪
⎨b2 = y3 − y1
⎧ h1.b1 + h2 .b2 + h3 .b3 1 3 ⎪b = y − y
⎪ p = = .∑ bk .hk ⎩ 3
⎪ ∆ ∆ k =1 1 2

⎨ ⎧c1 = x3 − x2
⎪q = h1.c1 + h2 .c2 + h3 .c3 = 1 . c .h
3
⎪
⎪⎩ ∆
∑
∆ k =1
k k ⎨c2 = x1 − x3
⎪c = x − x
⎩ 3 2 1

avec ∆ = x1.( y2 − y3 ) + x2 .( y3 − y1 ) + x3 .( y1 − y2 ) = x1.b1 + x2 .b2 + x3 .b3

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

10
 03/02/2010

Implémentation en pratique (6)

en entrant p et q dans Uj :

1 1 ⎡⎛ 3 ⎞⎛ 3 ⎞ ⎛ 3 ⎞⎛ 3 ⎞⎤ I j . A j 3
U j = .T j . 2 . A j ⎢⎜ ∑ bk .hk ⎟.⎜ ∑ bl .hl ⎟ + ⎜ ∑ ck .hk ⎟.⎜ ∑ cl .hl ⎟⎥ − .∑ hk
2 ∆ ⎣⎝ k =1 ⎠ ⎝ l =1 ⎠ ⎝ k =1 ⎠ ⎝ l =1 ⎠⎦ 3 k =1

1 1 ⎡⎛ 3 ⎞⎛ 3 ⎞ ⎛ 3 ⎞⎛ 3 ⎞⎤ I j . A j 3
U j = .T j . ⎢⎜ ∑ bk .hk ⎟.⎜ ∑ bl .hl ⎟ + ⎜ ∑ ck .hk ⎟.⎜ ∑ cl .hl ⎟⎥ − .∑ hk
4 ∆ ⎣⎝ k =1 ⎠ ⎝ l =1 ⎠ ⎝ k =1 ⎠ ⎝ l =1 ⎠⎦ 3 k =1

1 3 3 3
U j = .∑∑ Pkl .hk .hl − ∑ Qk .hk
2 k =1 l =1 k =1

Tj 1
avec Pkl = .(bk .bl + ck .cl ) et Qk = .I j . A j
2.∆ 3
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Implémentation en pratique (7)

les coefficients Pkl sont uniquement dépendants des
coordonnées des noeuds et de la transmissivité de
l’élément concerné
les coefficients Qk sont uniquement dépendants des
coordonnées, de la valeur d’infiltration et d’une valeur
moyenne de la hauteur piézométrique sur l’élément fini
concerné
Pour l’ensemble du domaine, à la place de k = 1 à 3 (ou l = 1 à 3 ),
on tient compte de la contribution des n noeuds des m éléments
pour obtenir la valeur totale de la fonctionnelle U:

1 n n n
U = .∑ ∑ Pkl .hk .hl − ∑ Q k .hk
2 k =1 l =1 k =1

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

11
 03/02/2010

Implémentation en pratique (8)

le minimum de cette fonctionnelle sera trouvé pour
∂U ∂hi = 0.
Cette dérivée est calculée terme par terme:
∂ ⎛1 n n ⎞ 1⎛ n n ⎞
⎜ .∑∑ Pkl .hk .hl ⎟ = .⎜⎜ ∑ Pji .h j + ∑ Pij .h j ⎟⎟
∂hi ⎝ 2 k =1 l =1 ⎠ 2 ⎝ j =1 j =1 ⎠
∂ ⎛ n ⎞ ∂h avec Pij = Pji
⎜ ∑ Qk .hk ⎟ = Qi . i = Qi
∂hi ⎝ k =1 ⎠ ∂hi
… le système d’équations algébriques:

∂U ⎛ n ⎞
= ⎜⎜ ∑ Pij .h j ⎟⎟ − Qi = 0
∂hi ⎝ j =1 ⎠

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Implémentation en pratique (9)

⎛ n ⎞
⎜ ∑ Pij .h j ⎟ = Qi
⎜ ⎟
⎝ j =1 ⎠
système de n équations algébriques avec n inconnues :
h1, h2, …, hn. ⎛ P11 P12 ... ... ... ... ... P1n ⎞ ⎛ h1 ⎞ ⎛ Q1 ⎞
⎜ ⎟⎜ ⎟ ⎜ ⎟
⎜ P21 P22 ... ... ... ... ... ... ⎟ ⎜ h2 ⎟ ⎜ Q2 ⎟
⎜ ... ... ⎟ ⎜ ... ⎟ ⎜ ... ⎟
⎜ ⎟⎜ ⎟ ⎜ ⎟
⎜ ... ⎟.⎜ ... ⎟ = ⎜ ... ⎟
⎜ ⎟⎜ ⎟ ⎜ ⎟
⎜ ... ⎟ ⎜ ... ⎟ ⎜ ... ⎟
⎜ ... ⎟ ⎜ ... ⎟ ⎜ ... ⎟
⎜⎜ ⎟⎟ ⎜ ⎟ ⎜ ⎟
⎝ Pn1 ... ... ... ... ... ... Pnn ⎠ ⎝ hn ⎠ ⎝ Qn ⎠
avec Pij matrice de ‘rigidité’

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

12
 03/02/2010

Implémentation en pratique: exemple

2 4 6
h1 = h2 = 10 m
1 3
h =h =0m
10 0 m T5 = T6 = T = T = 2 m2/jour
1 2 3 4
m 2 4 h3 = ? h4 = ? ⎧b1 = y2 − y3
⎪
1 3 5
Tj ⎨b2 = y3 − y1
Pkl = .(bk .bl + ck .cl ) ⎪⎩b3 = y1 − y2
Elément 1
2.∆ ⎧c1 = x3 − x2
⎪
2 (2) 4
(3) ∆ = x1b1 + x2b2 + x3b3 ⎨c2 = x1 − x3
⎪c = x − x
1 x1 = 0 y1 = 0 b1 = 0 c1 = 1 ⎩ 3 2 1
x2 = 0 y2 = 1 b2 = 1 c2 = -1
x3 = 1 y3 = 1 b3 = -1 c3 = 0 P ij
matrice 6 x 6
∆ = -1 A=½ ⎛ 1 −1 0 ⎞
1 (1) ⎜ ⎟
⎜ −1 −1 ⎟
⎛ 1 −1
2
0⎞ ⎜ ⎟
⎜ ⎟ ⎜ ⎟
P = ⎜ −1
1
2 −1 ⎟ ⎜ 0 −1 1 ⎟
⎜ 0 −1 1 ⎟ ⎜ ⎟
⎝ ⎠ ⎜
⎜ Hydrogéologie
Terug
⎟ Géologie
⎟
et
naar eerste pagina
⎝ ⎠
de l’Environnement

Implémentation en pratique: exemple

2 4 6
Tj
Pkl = .(bk .bl + ck .cl )
1 3 2.∆
⎧b1 = y2 − y3
10 0m ∆ = x1b1 + x2b2 + x3b3 ⎪
⎨b2 = y3 − y1
m 2 4 ⎪b = y − y
⎩ 3 1 2

1 3 5 ⎧c1 = x3 − x2
Elément 2 ⎪
4 ⎨c2 = x1 − x3
(3) x1 = 0 y1 = 0 b1 = -1 c1 = 0 ⎪c = x − x
⎩ 3 2 1
x2 = 1 y2 = 0 b2 = 1 c2 = -1
x3 = 1 y3 = 1 b3 = 0 c3 = 1
2 ∆=1 A=½ ⎛⎛ 21+1 − 1 −-11 0+0 0 ⎞ ⎞
⎜⎜ ⎟ ⎟
1(1) 3 (2) ⎛ 1 −1 0⎞ ⎜⎜−−11 2 − 1− 1 ⎟ ⎟
⎜ ⎟ ⎜⎜−-11 22 − -11 ⎟ ⎟
P = ⎜ −1
2
2 −1 ⎟ ⎜⎜ ⎟ ⎟
⎜ 0 −1 1 ⎟ ⎜⎜ 00 − 1 −-11 1+12
+0 ⎟ ⎟
⎝ ⎠ ⎜⎜ ⎟ ⎟
⎜⎜ ⎟ ⎟
⎜⎜ ⎟ ⎟
⎝⎝ ⎠ ⎠
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

13
 03/02/2010

Implémentation en pratique: exemple

2 4 6

1 3 ⎧b1 = y2 − y3
10 0m ⎪
⎨b2 = y3 − y1
m Tj ⎪b = y − y
2 4 Pkl = .(bk .bl + ck .cl ) ⎩ 3 1 2
2.∆ ⎧c1 = x3 − x2
1 3 5 ⎪
∆ = x1b1 + x2b2 + x3b3 ⎨c2 = x1 − x3
Elément 3 ⎪c = x − x
4 (2) 6 ⎩ 3 2 1
(3)

⎛ ⎛ 2 2 −−1 1 −−11 00
3
⎞⎞
⎜⎜ ⎟⎟
⎛ 1 −1 0⎞ ⎜ ⎜−−1 1 2 2 −−11 ⎟⎟
⎜ ⎟ ⎜ ⎜−−1 1 ⎟⎟
3(1) P = ⎜ −1
3
2 −1 ⎟ ⎜⎜
32+1 −−21-1 00
⎟⎟
⎜ 0 −1 1 ⎟ ⎜ ⎜ 0 0 −−1 1 −−21-1 42+2 −-11⎟ ⎟
⎝ ⎠ ⎜⎜ ⎟⎟
⎜⎜ ⎟⎟
⎜⎜ ⎟⎟
⎝⎝ 00 - 1-1 1⎠ ⎠
1

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Implémentation en pratique: exemple

2 4 6
Tj
Pkl = .(bk .bl + ck .cl )
1 3 2.∆ ⎧b1 = y2 − y3
⎪
10 0m ∆ = x1b1 + x2b2 + x3b3 ⎨b2 = y3 − y1
m ⎪b = y − y
2 4 ⎩ 3 1 2

⎧c1 = x3 − x2
1 3 5 ⎪
⎨c2 = x1 − x3
Elément 4 ⎪c = x − x
6 ⎩ 3 2 1
(3)

⎛ 2 −1 −1 0 ⎞⎞
⎜ ⎟⎟
⎜ −1 2 −1 ⎟⎟
5 ⎜ −1 34 +1
−2 −-11 0 ⎟0 ⎟
+0
3(1) (2) ⎜ ⎟⎟
⎜ 0 −1 − 2 4 − 1− ⎟1 ⎟
⎛ 1 −1 0⎞ ⎜ −-11 22 -1− ⎟⎟1 ⎟⎟
⎜ ⎟ ⎜
P = ⎜ −1
4
2 −1 ⎟ ⎜ 0+0 - 1 −-11 1+1⎟⎠2 ⎟⎠
⎝
⎜ 0 −1 1 ⎟
⎝ ⎠
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

14
 03/02/2010

Implémentation en pratique: exemple

4 6
2 ⎛ 2 −1 −1 0 ⎞
⎜ ⎟
1 3 ⎜ −1 2 −1 ⎟
10 0m ⎜ −1 4 −2 −1 0⎟
⎜ ⎟
m 2 4 ⎜ 0 −1 − 2 4 −1⎟
⎜ −1 2 − 1 ⎟⎟
1 3 5 ⎜
⎜ −1 2 ⎟⎠
⎝ 0 -1

⎧− h1 + 0.h2 + 4h3 − 2h4 − h5 + 0.h6 = 0 ⎧− 10 + 4h3 − 2h4 = 0

⎪0.h − h − 2h + 4h − 0.h − h = 0 ⎪− 10 − 2h + 4h = 0
⎪ 1 2 3 4 5 6 ⎪
⎨ ⎨
3 4
h
⎪ 1 = h = 10 h = h = 10
2
⎪ 1 2
⎪⎩h5 = h6 = 0 ⎪⎩h5 = h6 = 0

Méthode itérative (Gauss-Seidel) :

h3 = 0 3.75 4.69 4.92 4.98 4.995 4.999 …
h4 = 2.5 4.375 4.84 4.96 4.99 4.998 …

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Généralisation
Les variables dépendantes aux noeuds sont appelées les
degrés de liberté de l’élément:
les valeurs piézométriques aux 3 noeuds

… le système symétrique d’équations algébriques

pour un écoulement en régime permanent :

[S ].{h} = {Q}+ {G}+ {B}

avec [S ] matrice de rigidité
{Q} vecteur source/ disparition

{G} vecteur de la pesanteur

{B} vecteur des conditions aux frontières

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

15
 03/02/2010

Généralisation (2)
Pour un écoulement permanent, on peut
généraliser la méthodes des éléments :
n
h( x, y , z ) ≅ ∑ h j . N j ( x, y , z )
j =1

avec hj la hauteur piézométrique au noeud j

Nj la fonction d’interpolation, de base, ou de
forme relative au noeud j
… le plus souvent une forme particulière de la fonction de
forme est proposée. Les fonctions de forme Nj prennent
la valeur 1 au noeud j et une variation linéaire sur tous
les éléments qui touchent ce noeud j. En dehors: Nj = 0

la matrice de rigidité sera une matrice-bande

propriété très importante pour la résolution numérique du
système
Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Ecoulement transitoire
„ Schéma d’intégration temporel par différences finies
„ Méthode de Galerkin + schéma d’intégration de Galerkin
„ Des fonctions de base et de forme sont développées en
incluant l’intégration temporelle par différences finies:
n
h( x, y, z , t ) ≅ ∑ h j (t ) . N j ( x, y, z )
j =1

… en général :

[M ]. ⎧⎨ dh ⎫⎬ + [S ].{h} = {Q}+ {G}+ {B}

⎩ dt ⎭
avec [M ] matrice du terme d’emmagasinement
⎧ dh ⎫
⎨ ⎬ vecteur des n dh dt
⎩ dt ⎭ Hydrogéologie et Géologie
Terug naar eerste pagina
de l’Environnement

16
 03/02/2010

Discrétisation : FDM vs FEM

FEM : Maillages triangulaires (irréguliers)

FDM : Maillage orthogonaux (réguliers) Hydrogéologie et Géologie
Terug naar eerste pagina
de l’Environnement

Comparaison : FDM vs FEM

FDM : Concept simple, facile à programmer

Plus difficile de définir géométries irrégulières
Raffinement du maillage peu efficace

FEM : Solution continue dans le domaine

Flexibilité pour géométries irrégulières
Plus complexe à programmer (connectivité
des éléments)

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

17
 03/02/2010

Maillage
„ Orientation selon anisotropie (FDM)
‹ Edwards Aquifer (roc fracturé)

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Maillage
„ Orientation (FEM) Edwards Aquifer

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

18
 03/02/2010

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

Hydrogéologie
Terug et Géologie
naar eerste pagina
de l’Environnement

19