13.07.2015 Views

PDF (Intro, Chapitre 1, 2) - Les thèses en ligne de l'INP - Institut ...

PDF (Intro, Chapitre 1, 2) - Les thèses en ligne de l'INP - Institut ...

PDF (Intro, Chapitre 1, 2) - Les thèses en ligne de l'INP - Institut ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Chapitre</strong> 1. Médium Plan– Pour les faibles porosités, l’influ<strong>en</strong>ce <strong>de</strong>s fibres voisines n’est pas représ<strong>en</strong>tée uniquem<strong>en</strong>tpar la porosité. Dodd et al. [33] ont pris <strong>en</strong> compte cet aspect (cité par Li et Park [78]).– a=1 correspond au modèle <strong>de</strong> Spielman et Gor<strong>en</strong>.La transition <strong>en</strong>tre les <strong>de</strong>ux types d’écoulem<strong>en</strong>t est réalisée grâce à la continuité <strong>de</strong>s vitesseset <strong>de</strong>s contraintes tang<strong>en</strong>tielles et normales.Fibres parallèles à l’écoulem<strong>en</strong>t Il est nécessaire <strong>de</strong> résoudre l’équation <strong>de</strong> Brinkman pour¯r > a, avec ¯r = r R, pour un écoulem<strong>en</strong>t parallèle à la fibre. Soit, après adim<strong>en</strong>sionnem<strong>en</strong>t <strong>de</strong>svariables :(1 d¯r dū )= d ¯P¯r d¯r d¯r d¯z + α2 ū (1.11)où les variables avec une barre sont adim<strong>en</strong>sionnées : ū = u u fet ¯P = Rµu fPLa résolution <strong>de</strong> l’équation différ<strong>en</strong>tielle (1.11) permet <strong>de</strong> déterminer la vitesse :ū = CI 0 (α 2 ‖¯r) + DK 0(α 2 ‖¯r) − 1 α 2 ‖I 0 et K 0 sont les fonctions <strong>de</strong> Bessels modifiées à l’ordre 0 <strong>de</strong> première et secon<strong>de</strong> espècerespectivem<strong>en</strong>t et α ‖ = 1k ‖.Pour 1 < ¯r < a, nous avons :ū = ¯r24d ¯Pd¯z+ A ln ¯r + BA, B, C et D sont <strong>de</strong>s constantes d’intégration. Celles-ci sont déterminées grâces aux conditionsaux limites <strong>en</strong> ¯r = 1, ¯r = a et pour la limite <strong>en</strong> ¯r <strong>en</strong> l’infini. Voir <strong>de</strong> Li et Park [78] pourplus <strong>de</strong> détails.Il est possible alors <strong>de</strong> déduire la force <strong>de</strong> trainée F d , par unité <strong>de</strong> longueur, <strong>de</strong> l’écoulem<strong>en</strong>tsur le cylindre :¯F d = 2π dū∣d¯rLa chute <strong>de</strong> pression est alors reliée à F d comme pour les milieux cellulaires, voir équation(1.8). Ainsi, l’expression <strong>de</strong> la perméabilité d’un médium avec les fibres parallèles à l’écoulem<strong>en</strong>ts’écrit :∣r=1d ¯Pd¯z(11 − 1 ) ()1=2 1 − ε α‖2 − 1 4 + a24 − a22 ln a α ‖ aK 1 (α ‖ a)K 0 (α ‖ a) + α ‖ aK 1 (α ‖ a) ln a + a22(1.12)Nous obt<strong>en</strong>ons ainsi une fonction implicite <strong>de</strong> la perméabilité, avec K 1 est la fonction <strong>de</strong>Bessel modifiée à l’ordre 1 <strong>de</strong> secon<strong>de</strong> espèce. Le cas <strong>de</strong> Spielman et Gor<strong>en</strong> correpond au casparticulier a = 1. Dans ce cas, l’equation (1.12) se simplifie pour donner :11 − ε = 2K 1(α ‖ )α ‖ K 0 (α ‖ )(1.13)36

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

Saved successfully!

Ooh no, something went wrong!