//img.uscri.be/pth/f8c4a1c06220a7b57693e74ead703325a3c9bbc9
Cet ouvrage fait partie de la bibliothèque YouScribe
Obtenez un accès à la bibliothèque pour le lire en ligne
En savoir plus

Eléments finis stabilisés pour des écoulements diphasiques compressible-incompressible

De
193 pages
Sous la direction de Boniface Nkonga
Thèse soutenue le 27 novembre 2009: Bordeaux 1
Dans cette thèse, nous nous intéressons à la simulation numérique d'écoulements instationnaires de deux fluides visqueux non miscibles, séparés par une interface mobile. Plus particulièrement des écoulements sans choc constitués d'une phase gazeuse et d'une phase liquide sont considérés. Pour modéliser de tels écoulements, une approche dans laquelle le gaz est décrit par les équations de Navier-Stokes compressible et le liquide par les équations de Navier-Stokes incompressible est proposée. C'est le couplage de ces deux modèles qui constitue l'originalité et l'enjeu principal de de cette thèse. Pour traiter cette difficulté majeure, une méthode globale (i.e. la même dans chaque phase) et simple à mettre en oeuvre est élaborée. L'utilisation des équations de Navier-Stokes formulées de façon unifiée pour les inconnues primitives (pression, vitesse et température) constitue le point de départ pour la construction de notre méthode qui repose sur les composants suivants: une méthode d'éléments finis stabilisés pour la discrétisation spatiale des équations de Navier-Stokes; une approche Level Set pour représenter précisément l'interface dont l'équation de transport a été résolue par une méthode de type Galerkin Discontinu; et des grandeurs moyennes pour traiter les discontinuités à l'interface. Le bon comportement de notre approche est illustré sur différents tests mono et bi-dimensionnels.
-Ecoulement diphasique
-Interface
-Navier-Stokes Incompressible
-Navier-Stokes compressible
-Eléments finis stabilisés
-Méthode Level Set
In this work, we are interested in the numerical simulation of instationnary viscous flows of two immiscible fluids, separated by a mobile interface. In particular, flows without shock composed of a gas phase and a liquid phase are considered. In order to modelize such flows, an approach in which the gaz is described by compressible Navier-Stokes equations and the liquid by incompressible Navier-Stokes équations is proposed. The coupling between these two models is the originality and the stake of this thesis. To treat this important difficulty, a global (i.e. the same for each phase) and simple method is elaborated. In our procedure we propose, using the Navier-Stokes equations formulated in set of primitives unknowns (pressure, velocity and temperature), to elaborate a strategy that relies on the follow components: the stabilized finite element method to discretize spatially the Navier-Stokes equations; the Level Set method for tracking the interface precisely with a discontinuous Galerkin method to solve the associated transport equation; and some averaged quantities to treat the discontinuities at the interface. The good behaviour of this approach is performed on both one and two spatial dimensions.
Source: http://www.theses.fr/2009BOR13872/document
Voir plus Voir moins

[
g
^







O



H

W



p



d

H

{

h

g





x

K

b



K

m

J

^





m

w

K



q

d

K

=

l

d

h

|



N
!

"
w

{
#
K

H
$


h
%
q
&
E



+
'
b
(
f
)
K


*

(
h
'
q

b
(
g




)
K

O
*
m

[
+
ƒ

e


,
h

H


+
l

o

]
-
g
.
H
/
Ž

ˆ

r



\
-
K
0
H
(


Š
+


t

v
-
_
.
y
/
h

i
1

2
/
3

4
E
2

3

5
l
6
o
7
o
8
D
9

9
j
:
G
;

<
\

`
(
_
/
c

h
(
i
=





-
j

F
)
J

l
>
o

o

t
*
v

_
*
y

h
<
i
?

@
/
A

B
E
;

C


l

o

o

_
'

-
s

c
)
{
-




J

O


+

D

E
J
F


Š

q
G
[
D
]
E
\
H
e
I
-
J


J
K

D

L
h
L

M
p
H
+
I

J


{

h

i
K


/
E

N
E
D

M
h
N
q
J
r
M

K

N
]
H
u

a
O
t
E
x
H
f
P
z
I
|
J






K

J
Q
h

O

W



H
R




D

E

K

k

m
$
n
!
p

r

*

W

O

M

}
S

"
i




#




[

]
T
_
"
a

]

c

d
Q
e

g




#
K

J
S
h

N
$
M



M
U
E

N

O

I
T

"
k
Q
m
S
n

p

r

+

s

]
!
u

a

t
V
x

f

z
T
|
"

Q

S


K

J

h

O

W
!


H





D
J
E
E
K

k
O
m
J
n

p
M
r

+

~




W


c
G
y
X
\

f

|

h
Y

Y
i
Z
K


L


/

H


K
K
D
E
W

H

K


D



h
+
ˆ
[

\

]
p
^
r
_
+
`
\
a
_
b

]
Œ
_
c
c
d
d
s
d
h
c
O
e

f
O
g

h


'



'
i



K

K
l

Š
J
l

q
h
o

[
N
a

b
M


_



|
M
c
j

E



N
K
F
J
O
h
H
O
I
W
J


H
h

k

l
H
m

m

n
o
o

p
q
q
o
r

o

+

[
sK
C
l

K

K

F

8

;

G



B

=

K

K

b

8

^

:





K

K

K

K

K

?

;

H

K

K

K

K

@

V

?

<

K







9

K

]

K

;

K

I

;

K

B

K

B

K

9

O

r

;

K

b

K

K

;

K

:

K

9

K

9

?

8

K

K

K

K

h









(

B

K

K
!
K

?

K

K
"
H
#
:

I

K
$
]

;

K
%
O

B
&
?

K

H
'
K

K
(
O
)
:

>

K

O

H

q

K

;
*
F
+
K
,
O
,
G

K
,
K
-
F
+
T
,
K
.
;

?
)
K

K

=

Q

K
-
b

K
,
K
/
:
)
?
0
=

=
,
R

9
1
>
2
:

?

>
3
:
,
K
-
S

=

K
4
K
5
K
6
K
7
K
8

9

:

;

<

=

>

?

@
k
8
(
A
l
:
,
?
n
B
N
C
?
D
K
9
K
;
K
E
K
?
o
=
:
F
6
;
K
G
K
H
K
B
K
?
K
=
]
=
B
:
C
I
?
J
K
?
K
K
K
K
H
K
G
K
J
K
?
K
K
K
K
K
K
K
P
K
F
K
:
K
;
K
I
K
K
K
K
K
M
K
F
K
=
K
K
K
K
K
K
K
K
4
9
4
G
4
=
K
7
4
9
L
K
;
K
B
K
G
N
?
>
F
B
;
E
<
K
=
K
?
K
B
N
A
>
8
B
9
7
:
K
A
K
?
K
K
N
K
=
K
J
K
K
K
K
K
K
K
f
K
]
K
=
K
]
K
J
K
K
K
K
K
a
K
J
K
;
K
K
K
K
K
K
K
K
K
a
K
9
K
8
K
<
K
K
K
K
K
K
K
d
K
?
K
:
4
K
4
K
4
K
K
K
M
>
L
A
;
C
B
E
G
<
?
?
=
K
6
K
7
K
8
<
=
7
:
;
C
J
J
9
:
?
<
F
N
F
?
9
8
=
:
:
B
8
?
<
=
K
K
K
K
V
K
J
K
;
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
O
K
g
K

K
i
K

K

K
"
K

K

K

K

K

K
i
K

K

4

O

4

K
[
O
'
P

7

9

B
,
?

=
+
Q
4
;
I
B
G
G
6
?
K
=
K
6
K
7
K
8
K
=
K
:
K
C
K
J
K
:
4
<
=
N
H
?
<
8
]
:
?
B
=
?
K
=
K
K
K
K
K
K
K
K
K
K
K
K
M
K
P
K
;
K
:
K
G
K
=
K
J
K
<
K
H
K
=
K
?
K
K
K
K
K
K
K
K
K
M
K
P
K
;
K
F
4
B
R
:
M
C
5
H
6
=
7
?
8
K
9
K
:
K
;
K
<
K
=
K
>
H
?
]
@
G
8
=
A
J
:
<
?
H
B
=
C
?
D
K
9
K
;
K
E
K
?
a
=
a
:
B
<
?
F
H
;
?
G
K
H
K
B
K
?
K
=
K
=
K
:
K
I
a
J
N
?
;
K
<
K
?
K
7
K
;
K
?
K
A
K
:
K
B
K
?
K
K
K
K
K
K
K
K
K
a
K
4
K
9
K
?
K
9
K
q
K
<
K
8
4
<
S
K
M
K
K
K
4
K
L
K
;
a
B
M
G
9
7
?
J
9
8
q
9
<
:
H
;
:
<
K
=
K
9
K
8
K
<
K
>
a
8
O
B
9
>
?
K
<
K
:
K
r
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
?
K
>
K
;
K
?
K
?
K
>
K
A
K
?
K
K
K
K
K
K
K
K
K
K
K
a
K
G
K
9
K
r
4
8
S
K
M
K
K
K
M
K
T
K
:
K
G
K
:
K
9
R
?
4
:
>
<
;
F
?
;
9
G
c
H
?
B
F
?
K
=
K
=
K
:
K
I
K
J
K
?
K
>
>
?
?
=
t
N
?
?
K
6
K
7
K
8
K
9
K
:
K
;
K
<
R
=
:
>
?
?
8
@
:
8
B
A
D
:
;
?
?
B
7
C
:
D
N
9
?
;
K
E
K
?
K
=
O
F
M
;
M
G
;
H
=
B
B
?
F
=
:
=
<
:
?
I
8
J
8
?
B
K
F
K
>
K
=
4
;
U
e
O
:
V
<
;
=
<
?
F
9
J
I
7
J
=
=
:
9
;
;
<
K
K
K
K
K
K
K
K
O
K
O
K
;
K
F
K
7
K
:
K
<
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
f
K

K

K

K

K

K

K

K

M

4

W

X

+
i
,



Y





,



,

-

+

,

(

+

,

'

,


j

"
-

*




%
-
&
-



)
)
Z



(



-
Z
+
m
,
'
-

W

[
(
4
,
\

N
_
?
;
<
J
N
?
?
8
B
:
8
7
J
K
:
K
9
K
N
K
?
K
=
K
=
K
7
K
B
K
J
K
8
K
G
K
N
K
?
K
9
K
]
K
;
K
>
K
?
a
>
M
?
?
=
F
N
:
?
9
J
;
N
>
?
F
G
8
?
7
<
H
9
8
=
?
^
K
<
K
:
K
=
K
=
K
9
K
8
K
I
K
:
K
J
K
:
K
=
K
N
K
?
K
=
a
K
M
K
4
K
H
K
B
K
F
K
?
K
<
K
;
K
H
K
?
K
=
K
I
K
?
K
:
M
F
O
G
4
B
K
=
4
:
_
J
B
K
;
K
I
K
J
K
`
K
?
K
G
K
?
K
=
K
G
a
;
M
>
M
`
H
?
B
J
F
?
?
=
;
K
H
K
?
K
=
K
I
K
?
K
F
K
G
K
B
K
=
K
:
K
J
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
a
K
M
K
O
K
H
K
B
K
F
K
?
K
;
K
H
M
?
a
=
4
I
K
?
M
:
b
F
N
G
?
B
9
=
]
:
;
J
>
K
?
K
>
K
?
K
\
K
8
K
J
K
?
K
B
K
E
K
:
O
<
K
K
P
K
H
K
;
K
]
K
H
K
;
K
;
K
N
K
?
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
a
K
b
K
?
K
]
K
>
M
=
S
7
4
N
K
B
O
6
@
?
N
H
?
7
F
J
?
=
=
:
=
:
:
c
9
<
N
?
?
Q
>
F
?
K
=
K
9
K
8
K
I
K
:
K
J
K
:
K
=
K
?
K
B
R
K
K
K
b
K
?
K
]
K
>
K
=
K
?
K
p
K
?
K
Q
K
;
K
9
K
B
K
F
K
:
K
r
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
R
M
K
S
b
4
?
K
]
a
>
b
=
N
?
?
p
9
?
]
Q
;
;
>
9
?
8
=
9
>
B
c
<
N
q
?
K
J
K
N
K
?
K
G
K
?
K
<
K
9
K
=
K
^
K
<
S
:
K
=
b
=
?
9
]
8
>
I
=
:
8
J
=
:
8
=
J
N
8
?
?
=
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
K
a
K
O
K
a
K
N
K
9
K
;
K
?
K
]
K
:
O
:
d
s
M
N
5
9
J
;
N
?
?
?
G
?
?
D
<
9
9
K
=
K
^
K
<
K
:
K
=
K
=
K
9
K
8
K
I
K
:
f
J
L
:
B
=
7
N
8
?
:
=
<
H
J
;
I
7
J
B
K
J
K
?
K
=
K
N
K
?
K
6
K
7
K
8
K
9
K
:
K
;
K
<
K
=
K
>
K
?
K
@
K
8
d
A
K
:
V
?
<
B
:
D
:
9
<
;
>
E
=
?
7
=
`
K
J
K
:
K
9
K
B
K
8
K
?
K
K
K
K
K
K
K
K
K
K
O
K
R
K
M
K
K
K
4
K
L
K
;
R
B
a
G
M
7
;
J
`
8
J
9
I
:
C
;
7
<
>
\
K
T
K
D
K
H
K
;
K
7
K
B
K
J
K
?
K
=
K
N
K
?
K
6
K
7
K
8
K
9
K
:
4
;
:
<

=2

?





7