01.05.2013 Views

Modélisation, analyse mathématique et simulations numériques de ...

Modélisation, analyse mathématique et simulations numériques de ...

Modélisation, analyse mathématique et simulations numériques de ...

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.

tel-00656013, version 1 - 3 Jan 2012<br />

3 Partie III : un modèle d’écoulement sanguin 35<br />

L’idée est <strong>de</strong> se diriger vers un modèle d’écoulement plus réaliste, à savoir :<br />

• une partie permanente établie : profil <strong>de</strong> Poiseuille,<br />

• plus une perturbation périodique en temps due au pouls : profil <strong>de</strong> Wommersley.<br />

Précisément, nous considérons la géométrie plus simple <strong>de</strong> [40], illustrée à la Figure 1.7.<br />

Le problème <strong>de</strong> départ s’énonce comme suit.<br />

Trouver uε tel que : ⎧⎪<br />

∂tuε −∆uε +∇pε = 0 dans Ωε<br />

div uε = 0 dans Ωε<br />

⎨ uε = 0 sur Γ1 ∪Γε<br />

(3.2)<br />

⎪⎩<br />

pε = pin(t) sur Γin<br />

pε = pout = 0 sur Γout<br />

uε est x1-periodique<br />

On suppose que la pression est périodique en temps <strong>et</strong> ne dépend que <strong>de</strong> x1 en espace :<br />

on peut se restreindre à la résolution d’un problème scalaire sur la vitesse horizontale uε.<br />

Avant <strong>de</strong> suivre les étapes classiques pour l’obtention <strong>de</strong> lois <strong>de</strong> paroi (présentées brièvement<br />

à la section précé<strong>de</strong>nte), le point <strong>de</strong> départ est un passage en série <strong>de</strong> Fourier en temps,<br />

possible grâce aux hypothèses <strong>de</strong> périodicités faites. Nous sommes ainsi ramenés à résoudre<br />

le problème suivant, pour tout mo<strong>de</strong> k ∈ Z ∗ : trouver ûǫ,k tel que<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

(ik −∆)ûǫ,k = Ĉk dans Ωε,<br />

ûǫ,k = 0 sur Γε ∪Γ1,<br />

ûǫ,k x1 − périodique sur Γin ∪Γout.<br />

(3.3)<br />

où Ĉk est le k-ième coefficient <strong>de</strong> Fourier <strong>de</strong> la pression entrante. Les étapes suivantes, à<br />

savoir ansatz sur le développement <strong>de</strong> la solution, approximation d’ordre 0, û0,k, introduction<br />

<strong>de</strong> correcteurs, approximation d’ordre 1, Ûǫ,k (contenant la variable lente <strong>et</strong> la variable<br />

rapi<strong>de</strong>) dans le domaine rugueux, puis sa moyennisation par rapport à la variable rapi<strong>de</strong><br />

uε,k, <strong>et</strong> enfin la solution ˆ Vε,k <strong>de</strong> la loi <strong>de</strong> paroi implicite (<strong>de</strong> type Navier) dans le domaine<br />

lisse sont plus ou moins conservées. Il est toutefois intéressant <strong>de</strong> souligner <strong>de</strong>ux faits :<br />

• les mo<strong>de</strong>s <strong>de</strong> Fourier n’intéragissent pas avec la fréquence d’oscillation <strong>de</strong> la rugosité<br />

(heureusement!);<br />

• <strong>et</strong> dans ce cadre <strong>de</strong> modèle jou<strong>et</strong>, on obtient toutes les estimations, à chaque étape du<br />

processus d’approximation, <strong>de</strong> manière directe, sans invoquer <strong>de</strong> résultats théoriques

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

Saved successfully!

Ooh no, something went wrong!