Résolution des équations pour la dose.

Données physiques utilisées:

Les données sont inspirées du brouillon d'article que tu m'as passé (le deuxième article). Pour la masse m, je voulais avoir une fonction constante pour la condition aux bords, ce qui explique la définition choisie.

rayon de la lumière R: 15.0
épaisseur de la paroi E: 4.0
distance centre à centre des struts L: 7.0
longueur des struts : 1.5
diffusion: Dp = Dl = Id
nombre de struts (N) et masse (m) N = 4, m = 2e-6*(RN )
vitesse de l'écoulement u: u = (0,0,k(1-(x2+y2)/R2)), k=2

Mise en oeuvre

On formule le probleme en axi-symétrique. On programme dans MEF++ ce qui permet:

J'ai pris en aval et en amont de la section contenant le stent une section de même longueur (donc la longueur total du domaine est de trois fois la "longueur" du stent). Je n'ai pas fait de traitement particulier dans ces zones (maillage adapté etc). J'ai utilisé 6 maillages (structurés/non structurés):

Maillage nbre d'éléments nbre de sommets
Maillage non structuré artert 2244 1199
artert2 (raffinement de artert) 8976 4645
artert4 (raffinement de artert2) 35904 18269
Maillage structuré arterq 3078 1648
arterq2 (raffinement de arterq) 12312 6377
arterq4 (raffinement de arterq2) 49248 25069

Ici les temps de calculs ne sont pas à considèrer comme un facteur important. J'ai obtenu les solutions sur mon PC (linux 233 MHz) en moins de 3 minutes pour le plus gros cas. J'ai choisi une interpolation quadratique. Le problème des valeurs négatives semblent avoir disparu avec l'utilisation de ces maillages. Cliquer sur l'image pour l'agrandir.


Solution de
artert

Solution de
artert2

Solution de
artert4

Isovaleurs de
artert4
Isovaleurs de
artert4

Solution de
arterq

Solution de
arterq2

Solution de
arterq4

Isovaleurs de
arterq4
Isovaleurs de
arterq4