27.06.2013 Views

Simulation numérique de l'essorage et du refroidissement d'un film ...

Simulation numérique de l'essorage et du refroidissement d'un film ...

Simulation numérique de l'essorage et du refroidissement d'un film ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

74 Chapitre 3 Description <strong>et</strong> validation <strong>de</strong>s métho<strong>de</strong>s <strong>numérique</strong>s<br />

où L0, t0, u0 <strong>et</strong> p0 sont respectivement la longueur <strong>de</strong> référence spatiale, le temps, la<br />

vitesse <strong>et</strong> la pression. La constante K assure que le terme <strong>de</strong> lagrangien augmenté est<br />

<strong>de</strong>ux ou trois ordres <strong>de</strong> magnitu<strong>de</strong> plus grand que les autres contributions physiques dans<br />

les équations <strong>du</strong> mouvement. La métho<strong>de</strong> <strong>numérique</strong> <strong>de</strong> minimisation pour résoudre le<br />

lagrangien est l’algorithme d’Uzawa.<br />

Remarque : Dans la plupart <strong>de</strong>s cas, 2 itérations <strong>du</strong> lagrangien suffisent à obtenir un<br />

rési<strong>du</strong> <strong>et</strong> une divergence corrects. Dans tous les cas il est intéressant d’imposer un nombre<br />

d’itérations <strong>du</strong> lagrangien plus important pour améliorer la robustesse <strong>du</strong> calcul. Cepen-<br />

dant, le calcul est d’autant plus long que le nombre d’itérations <strong>du</strong> lagrangien est impor-<br />

tant. Un bon compromis est <strong>de</strong> conserver 2 itérations <strong>du</strong> lagrangien mais d’augmenter le<br />

nombre d’itérations internes <strong>du</strong> solveur à 50 itérations pour obtenir une meilleure résolu-<br />

tion implicite <strong>de</strong>s équations.<br />

Schémas temporels<br />

Les équations sont discrétisées implicitement en temps par le schéma d’Euler d’ordre<br />

1 pour les écoulements stationnaires <strong>et</strong> par le schéma <strong>de</strong> Gear d’ordre 2 pour les écoule-<br />

ments instationnaires. La discrétisation est écrite <strong>de</strong> façon générale par l’utilisation <strong>de</strong>s<br />

coefficients g1, g2 <strong>et</strong> g3 :<br />

n+1 ∂φ<br />

=<br />

∂t<br />

g1φn+1 + g2φn + g3φn−1 ∆t<br />

(3.7)<br />

Un schéma d’Euler correspond à g1=1 <strong>et</strong> g2=-1, tandis qu’un schéma <strong>de</strong> Gear corres-<br />

pond à g1 = 3/2, g2 = −2, g3 = 1/2.<br />

Schémas spatiaux pour le terme d’inertie<br />

Concernant le traitement <strong>de</strong> l’inertie, il est connu que le schéma centré est précis au<br />

second ordre <strong>et</strong> préconisé pour les écoulements turbulents car seule la diffusion moléculaire<br />

intervient dans la dissipation <strong>de</strong> l’énergie. Cependant, il génère <strong>de</strong>s instabilités dès qu’il<br />

y a <strong>de</strong> forts gradients <strong>et</strong> pro<strong>du</strong>it <strong>de</strong>s on<strong>de</strong>s <strong>de</strong> dispersion. Sur les surfaces libres, <strong>de</strong> forts<br />

gradients <strong>de</strong> masse volumique in<strong>du</strong>isent <strong>de</strong> fortes disparités <strong>de</strong> flux dans l’équation <strong>de</strong><br />

Navier-Stokes.<br />

Nous utilisons dans un souci <strong>de</strong> stabilité un schéma hybri<strong>de</strong> Upwind/Centré. Le schéma<br />

hybri<strong>de</strong> est plus diffusif que le schéma centré mais plus stable pour <strong>de</strong>s calculs en dipha-<br />

sique [Patankar 80].

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!