84
3 Propri´ et´ es g´ en´ erales des ´
equations
repr´ esente en effet une bonne approximation de la solution `
a divergence non nulle.
⎧
⎪ ⎪ ⎨
⎪ ⎪ ⎩
ρ
� V ∗ − V n
Δt
+ V
n · ∇V
∗
�
− ∇ ·
� μ
� ∇V
∗
+ ∇
t
V
∗
��
− r ∇ (∇ · V
∗
) = −∇p
n
+ f
p
n+1
= p
n − r ∇ · V
∗
Suivant la raideur du probl` eme physique pos´ e, l’augmentation du param` etre r
permettra de converger vers la solution avec des pas de temps raisonnables. Dans
notre m´ ethode, `
a la diff´ erence des m´ ethodes de projection classiques, seule la correction de la vitesse V � sera prise en compte pour calculer le champ `
a divergence
nulle, V n+1 sera calcul´ ee directement :
V
n+1
= V
∗
+ V
�
L’introduction de ce changement de variable dans le syst` eme d’´ equations de
Navier-Stokes conduit, apr` es diff´ erence entre les deux syst` emes d’´ equations portant
sur V n+1 et sur V ∗ , et en prenant r → ∞, `
a la forme finale de l’´ equation correspondant `
a l’´ etape de correction :
∇
� ∇ · V
n+1
� = −∇ (∇ · V
∗
)
La divergence du champ de vitesse V � qui est la solution de cette ´
equation n’est
d´ efinie qu’` a une constante additive pr` es. Cette constante est maintenue `
a z´ ero en
assurant ∇ · V ∗ = 0 en un point de Γ et ∇ · V � = 0 en ce mˆ eme point.
Cette formulation est ´
equivalente `
a V = P ⊥
H V ∗ o` u H = H(div, Ω ) = {V ∈ L 2
Ω /∇ ·
V = 0} est le sous-espace de L 2 `
a divergence nulle et P ⊥
H la projection orthogonale
de L 2 sur H.
Les champs de vitesses V n+1 et V ∗ satisfont tous deux les conditions aux limites
physiques du probl` eme. Nous pouvons en d´ eduire les conditions aux limites sur V �
qui sont des conditions aux limites homog` enes.
3.7 D´ eg´ en´ erescences des ´
equations de Navier-Stokes
Les ´
equations compl` etes de Navier-Stokes quelle que soit la formulation adopt´ ee
sont non lin´ eaires et difficiles `
a r´ esoudre. Dans certaines situations les contraintes
de l’´ ecoulement, niveaux de vitesses, ´
etirement dans une direction privil´ egi´ ee, permettent d’adopter des approximations qui conduisent `
a des formes plus simples des
´
equations. Celles-ci ne donnent plus des solutions exactes du probl` eme initial pos´ e
mais uniquement des approximations de cette solution. Ce sont des ”mod` eles” appel´ ees aussi d´ eg´ en´ erescences des ´
equations de Navier-Stokes et dont l’utilisation
doit ˆ
etre effectu´ ee `
a bon escient.
3 Propri´ et´ es g´ en´ erales des ´
equations
repr´ esente en effet une bonne approximation de la solution `
a divergence non nulle.
⎧
⎪ ⎪ ⎨
⎪ ⎪ ⎩
ρ
� V ∗ − V n
Δt
+ V
n · ∇V
∗
�
− ∇ ·
� μ
� ∇V
∗
+ ∇
t
V
∗
��
− r ∇ (∇ · V
∗
) = −∇p
n
+ f
p
n+1
= p
n − r ∇ · V
∗
Suivant la raideur du probl` eme physique pos´ e, l’augmentation du param` etre r
permettra de converger vers la solution avec des pas de temps raisonnables. Dans
notre m´ ethode, `
a la diff´ erence des m´ ethodes de projection classiques, seule la correction de la vitesse V � sera prise en compte pour calculer le champ `
a divergence
nulle, V n+1 sera calcul´ ee directement :
V
n+1
= V
∗
+ V
�
L’introduction de ce changement de variable dans le syst` eme d’´ equations de
Navier-Stokes conduit, apr` es diff´ erence entre les deux syst` emes d’´ equations portant
sur V n+1 et sur V ∗ , et en prenant r → ∞, `
a la forme finale de l’´ equation correspondant `
a l’´ etape de correction :
∇
� ∇ · V
n+1
� = −∇ (∇ · V
∗
)
La divergence du champ de vitesse V � qui est la solution de cette ´
equation n’est
d´ efinie qu’` a une constante additive pr` es. Cette constante est maintenue `
a z´ ero en
assurant ∇ · V ∗ = 0 en un point de Γ et ∇ · V � = 0 en ce mˆ eme point.
Cette formulation est ´
equivalente `
a V = P ⊥
H V ∗ o` u H = H(div, Ω ) = {V ∈ L 2
Ω /∇ ·
V = 0} est le sous-espace de L 2 `
a divergence nulle et P ⊥
H la projection orthogonale
de L 2 sur H.
Les champs de vitesses V n+1 et V ∗ satisfont tous deux les conditions aux limites
physiques du probl` eme. Nous pouvons en d´ eduire les conditions aux limites sur V �
qui sont des conditions aux limites homog` enes.
3.7 D´ eg´ en´ erescences des ´
equations de Navier-Stokes
Les ´
equations compl` etes de Navier-Stokes quelle que soit la formulation adopt´ ee
sont non lin´ eaires et difficiles `
a r´ esoudre. Dans certaines situations les contraintes
de l’´ ecoulement, niveaux de vitesses, ´
etirement dans une direction privil´ egi´ ee, permettent d’adopter des approximations qui conduisent `
a des formes plus simples des
´
equations. Celles-ci ne donnent plus des solutions exactes du probl` eme initial pos´ e
mais uniquement des approximations de cette solution. Ce sont des ”mod` eles” appel´ ees aussi d´ eg´ en´ erescences des ´
equations de Navier-Stokes et dont l’utilisation
doit ˆ
etre effectu´ ee `
a bon escient.