Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
Télécharger au format pdf ou txt
Vous êtes sur la page 1/ 22
Introduction aux volumes finis
Lyes DOUADJI Dpartement de Gnie lectrique & dlectronique Facult Sciences lIngnieur Universit de Boumerds
Introduction La mthode des Volumes Finis consiste intgrer, sur des volumes lmentaires, les quations crites sous forme intgrale. C'est une mthode particulirement bien adapte la discrtisation spatiale des lois de conservation, contrairement aux lments Finis, et est ainsi trs utilise en mcanique des fluides. Sa mise en oeuvre est simple si les volumes lmentaires ou "volumes de contrle" sont des rectangles en 2D ou des paralllpipdes en 3D. Cependant, la mthode des Volumes Finis permet d'utiliser des volumes de forme quelconque et donc de traiter des gomtries complexes, contrairement aux Diffrences Finies.
Volumes Finis pour une loi de conservation
Considrons une loi de conservation d'une grandeur physique w dans une maille de volume faisant intervenir un flux F(w) et un terme source S(w). Son expression sous forme intgrale est :
Appelons la surface de la maille, de normale extrieure n. Le thorme d'Ostrogradski conduit :
L'intgrale
reprsente la somme des flux travers chaque face de la maille. Le flux est suppos constant sur chaque face, l'intgrale se ramne une somme discrte sur chaque face de la maille. Il vient :
Volumes Finis pour une loi de conservation
La quantit Fface = F(wface) est une approximation du flux F sur une face de la maille, c'est le flux numrique sur la face considre.
La discrtisation spatiale revient calculer le bilan des flux sur une maille lmentaire. Ce bilan comprend la somme des contributions values sur chaque face de la maille. La manire dont on approche les flux numriques en fonction de l'inconnue discrte dtermine le schma numrique. L'criture du schma numrique peut galement utiliser des inconnues auxiliaires, par exemple le gradient de l'inconnue par maille.
Volumes Finis pour une loi de conservation
Explicitons maintenant le terme de drive temporelle. Un lment fondamental de la discrtisation en Volumes Finis est de supposer que la grandeur w est constante dans chaque maille et gale une valeur approche de sa moyenne sur la maille ou bien sa valeur au centre de la maille. D'autre part, le terme de drivation en temps est valu au moyen d'une mthode numrique d'intgration d'quation diffrentielle (Runge-Kutta, Euler explicite ou implicite...) et fait intervenir un pas de temps d'intgration t. Ce dernier peut tre constant ou variable. Pour fixer les ides, on crira la formulation avec une mthode d'Euler explicite. Notons w l'incrment de la grandeur w entre deux itrations temporelles successives. On peut ainsi crire :
Volumes Finis pour une loi de conservation
Finalement la loi de conservation discrtise avec la mthode des Volumes Finis peut s'crire :
Volumes Finis pour une loi de conservation
La mthodes des Volumes Finis consiste donc : Dcomposer la gomtrie en mailles lmentaires (laborer un maillage). Initialiser la grandeur w sur le domaine de calcul. Lancer le processus d'intgration temporelle jusqu' convergence avec : Calcul du bilan de flux par maille par un schma numrique. Calcul du terme source. Calcul de l'incrment temporel par une mthode numrique d'intgration. Application des conditions aux limites.
Cas monodimensionnel Considrons une loi de conservation 1D : O u est une grandeur physique fonction de la variable d'espace x et du temps t et f(u) est une fonction de u Le domaine de calcul est divis en N mailles de centre xi. Chaque maille a une taille hi = xi+1/2 - xi-1/2. Les indices demi-entier dsignent les interfaces de la maille avec les mailles voisines
Figure: Maillage 1D Le temps est discrtis en intervalles de pas constant t. La fonction u est suppose constante dans chaque maille et gale une valeur approche de la moyenne. Notons cette valeur moyenne dans la i-me maille de centre xi, l'instant t = n t. Ainsi :
Cas monodimensionnel Souvent, cette valeur approche de la moyenne est la valeur de la fonction u au centre xi de la maille, on parle alors de Volumes Finis Cell-Centered (et dans ce cas, Le discrtisation spatiale par les Volumes Finis consiste intgrer maille par maille la loi de conservation :
Soit pour la i-me maille de centre xi, au temps t = n t :
Ce qui s'intgre comme suit :
Cas monodimensionnel
dsigne une approximation du flux f(u) l'interface xi+1/2 et au temps n t. C'est le flux numrique au point xi+1/2. Ce flux numrique s'value en fonction des valeurs moyennes de u dans les mailles voisines, ce qui dtermine le schma numrique. Une mthode d'Euler explicite est utilise pour valuer la drive en temps (d'autres schmas peuvent tre utiliss, par exemple le schma de RungeKutta). La formulation discrtise en Volumes Finis de la loi de conservation est ainsi :
Cas bidimensionnel Considrons une loi de conservation d'une grandeur physique u(x; y; t) o x et y sont les deux directions d'espace. Le domaine gomtrique est divis en mailles lmentaires, par exemple en mailles rectangulaires comme reprsent sur la figure . La grandeur u est suppose constante dans chaque maille et gale une valeur approche de la moyenne sur la maille (ou encore la valeur au centre de la maille).
: reprsente la circulation sur un contour d'une maille lmentaire. Plaons-nous sur la maille de contour ABCD comme indiqu sur la figure. Le flux F est suppos constant sur chaque arte de la maille AB, BC, CD et AD.
Cas bidimensionnel L'intgrale se ramne une somme discrte sur chaque arte :
Ceci revient valuer le bilan des flux travers chaque facette de la maille.
Exemple simple 1D avec conditions de Dirichlet
Considrons l'quation diffrentielle suivante :
La fonction u(x) est suppos constante dans chaque maille et gale une valeur approche de la moyenne sur la maille considre La discrtisation spatiale par les Volumes Finis consiste intgrer maille par maille l'quation diffrentielle du problme, soit pour la i-me maille :
Ce qui donne aprs intgration :
pour i variant de 1 N
Exemple simple 1D avec conditions de Dirichlet
o dsigne la valeur moyenne de f sur la i-me maille:
Il reste maintenant exprimer u(xi-1/2) en fonction des inconnues ui. L'approximation la plus naturelle est de prendre la valeur moyenne de u(x) sur le segment [xi-1; xi], soit :
avec
Exemple simple 1D avec conditions de Dirichlet
Cette dernire expression n'est pas valable au bord gauche, pour i = 1, en x1/2 = 0, car elle fait intervenir le point x0 qui n'est pas dfini. Il se pose alors le problme du traitement des bords qui exige une formulation particulire. Une possibilit est de dfinir une maille fictive gauche de l'intervalle [0,1], et d'affecter une valeur moyenne de la fonction u dans cette maille. Une autre possibilit est de considrer la valeur moyenne de u(x1/2) non plus sur le segment [x0; x1] qui n'est pas dfini mais sur le segment [x1/2; x1], c'est ce que nous choisissons dans cet exemple. Ainsi on crit :
Exemple simple 1D avec conditions de Dirichlet
Et de mme pour le terme u(xi+1/2), on crit que:
La discrtisation en Volumes Finis est donc finalement :
Exemple simple 1D avec conditions de Dirichlet
Dans le cas particulier d'un maillage rgulier de pas h. La discrtisation en Volumes Finis devient :
Forme matricielle
Comparaison avec un schma aux Diffrences Finies
Nous introduisons un schma aux Diffrences Finies afin de comparer les deux mthodes de discrtisation. On se place sur le mme maillage construit pour les Volumes Finis. Les points xi seront considrs comme les noeuds du maillage pour les Diffrences Finies. Et la quantit ui dsignera alors la valeur de la fonction u au noeud x i. ATTENTION aux notations propres aux deux mthodes, dans un cas ui dsigne une valeur moyenne de u sur la i-me maille, et dans l'autre cas, ui dsigne la valeur de u en xi. L'quation rsoudre s'crit, sous forme discrte en chaque noeud xi :
Comparaison avec un schma aux Diffrences Finies
Approximons la drive seconde de u au moyen d'un schma l'ordre 2. On crit les dveloppements de Taylor de ui+1 et ui-1 au voisinage de xi :
Ce qui peut s'exprimer sous la forme :
Comparaison avec un schma aux Diffrences Finies
La discrtisation par des Diffrences Finies est la suivante :
Dans la cas particulier d'un maillage rgulier de pas h, l'quation discrtise s'crit :
Logiciels De nombreux codes de simulation numrique en mcanique des fluides reposent sur la mthode des diffrences finies : Fluent, StarCD, CFX, FineTurbo, elsA,
Rfrences C. Hirsch. Numerical computation of internal and external flows, volume I : fundamentals of numerical discretization. John Wiley & Sons, Chichester, New York, 1986.
