L'évaluation clinique de la performance

            valvulaire est basée sur la mesure de la vitesse maximale dans le jet transvalvulaire.

            Un modèle théorique basé sur la modélisation de profils minces et la théorie des

            écoulements potentiels est développé pour accéder aux vitesses dans ce jet. Les

            résultats sont validés par des simulations numériques. Ce modèle présente un excellent

            rapport  précision des résultats / temps CPU. On montrera comment son utilisation permet

            d'améliorer le diagnostic clinique

DEPLANO, Valérie

GUIVIER-CURIEN, Carine

PÉLISSIER, Robert

French

Colloque avec actes et comité de lecture. Internationale.International audienceL'évaluation clinique de la performance valvulaire est basée sur la mesure de la vitesse maximale dans le jet transvalvulaire. Un modèle théorique basé sur la modélisation de profils minces et la théorie des écoulements potentiels est développé pour accéder aux vitesses dans ce jet. Les résultats sont validés par des simulations numériques. Ce modèle présente un excellent rapport  précision des résultats / temps CPU. On montrera comment son utilisation permet d'améliorer le diagnostic clinique

Guivier-Curien, Carine

Deplano, Valérie

Pélissier, Robert

HAL AMU

Modèle d’écoulement autour d’une prothèse valvulaire pour le diagnostic clinique

I-Revues

19ème Congrès Français de Mécanique                                                                           Marseille, 24-28 août 2009 
  1 
Modèle d’écoulement autour d’une prothèse valvulaire pour le 
diagnostic clinique 
C. GUIVIER -CURIENa, V. DEPLANOa, R. PELISSIERa  
a. Institut de Recherche sur les Phénomènes Hors Equilibre (IRPHE UMR6594), Equipe de Biomécanique 
Cardiovasculaire, Marseille (France) 
Résumé :  
L’évaluation clinique de la performance valvulaire est basée sur la mesure de la vitesse maximale dans le jet 
transvalvulaire. Un modèle théorique basé sur la modélisation des profils minces et la théorie des 
écoulements potentiels est développé pour accéder aux profils de vitesses dans ce jet. Les résultats sont
validés par des simulations numériques. Ce modèle présente un excellent rapport entre la précision des
résultats et le temps de calcul CPU. On montrera comment son utilisation permet d’améliorer le diagnostic 
clinique. 
Abstract : 
The assessment of the valvular performance is mainly based on maximum velocity measurements recorded in 
the transvalvular jet. A theoretical model.based on the modelling of thin airfoils and the theory of potential 
flows is developed to assess the velocity profiles in this jet. The results are validated by numerical 
simulations. The advantage of this model is the good ratio that is obtained between the accuracy of the 
results and the CPU time. We will show how its use could improve the clinical diagnosis. 
Mots clefs : écoulement potentiel, modèle théorique, prothèse val ulaire biclapet, performance 
valvulaire 
1 Introduction 
Une prothèse commerciale biclapet, en position d’ouvert re maximale dans l’écoulement sanguin présente 
des similitudes avec deux profils minces placés dans u  écoulement potentiel. Sur la base de cette 
observation, un modèle théorique bidimensionnel de l’écoulement a été développé afin de rendre compte des 
profils de vitesse à l’entrée de la prothèse. L’évaluation de la performance valvulaire s’appuie aujourd’hui 
sur ces mesures de vitesses, et plus particulièrement sur la vitesse maximale enregistrée à l’entrée de la 
prothèse. Au-delà du modèle proposé, les résultats obtenus sont discutés en terme d’exploitation clinique et 
de diagnostic. 
2 Modèles et Méthodes 
2.1 Modèle théorique 
2.1.1 Hypothèses fluide et géométrie 
Le fluide est supposé incompressible et parfait. L’écoulement est 2-D et stationnaire. En amont de la 
prothèse, le profil de vitesse est de type plateau avec une amplitude U0=0.9m/s correspondant aux valeurs 
enregistrées dans la voie de chasse du ventricule ga che au pic systolique. 
La prothèse commerciale considérée est une prothèse de type St Jude. Les deux clapets sont respectivement 
nommés L1 et L2. Ils forment chacun un angle α de 5° avec l’axe longitudinal du repère principal, dans le 
plan d’intérêt XOZ (figure 1). Le bord d’attaque dechacun des clapets est l’origine d’un repère local lié à 
celui-ci. La distance entre les deux parois confinat l’écoulement est 2a. 
19ème Congrès Français de Mécanique                                                                           Marseille, 24-28 août 2009 
  2 
  
FIG. 1 – Géométries et notations. A gauche, prothèse bicapet St Jude. Au centre, schéma du modèle 
géométrique. Le repère principal est noté XOZ, les deux repères locaux attachés à chaque clapet sont notés 
X iOZi, avec i={1,2}. A droite, zoom sur le repère local lié à L1 et notations pour l’expression des vitesses. 
2.1.2 Expressions des composantes longitudinale Uz et transversale Ux de la 
vitesse en tout point du fluide 
Chaque clapet est modélisé par une distribution tourbillonnaire continue. A partir de l’expression de la 
fonction de courant ψi, il est alors possible [1] d’exprimer les composantes longitudinales Uzi,Li et 
transversales Uxi,Li de la vitesse induite par chacun des clapets en tout point P(xiP,ziP) du plan, dans le repère 
lié au clapet Li (i={1,2}). 
                   
 
Où γi est l’intensité de la circulation définie pour le clapet Li à partir d’un développement en série et de 
considérations géométriques. Cette intensité s’exprime par : 
 
Le confinement de l’écoulement dû à la présence des parois cardiaques et artérielles est pris en compte au 
travers la méthode des images [1]. Quatre générations de deux images chacune sont considérées dans 
l’expression finale des composantes de la vitesse. Finalement, le courant amont uniforme dirigé selon l’axe 
Oz s’ajoute aux vitesses induites pour exprimer les composantes de la vitesse Uz et Ux du fluide en tout point 
P(xP,zP) du repère principal : 
19ème Congrès Français de Mécanique                                                                           Marseille, 24-28 août 2009 
  3 
 
Où ILin
k désigne la kème image dans la nème génération image du clapet Li. ’implémentation du modèle 
théorique et les simulations sont réalisées à l’aide du logiciel Matlab. 
2.2 Modèle numérique 
Afin de valider les résultats obtenus, un modèle numérique a été développé. Pour s’assurer qu’une 
description des profils de vitesse par le modèle théorique 2-D est suffisante en vue d’une exploitation 
clinique, le modèle numérique géométrique choisi est 3-D. Les autres hypothèses géométriques et fluide ont 
été conservées. Les simulations ont été menées à l’aide du logiciel commercial Fluent pour résoudre les 
classiques équations de Navier-Stokes. 
3 Résultats 
Les valeurs des composantes de la vitesse longitudinale (Uz) et transversale (Ux) dans le plan XOZ obtenues 
par les deux modèles sont comparées.  
Sur la figure 2, sont représentées les valeurs de la composante Uz de la vitesse le long de l’axe longitudinal 
OZ. La différence moyenne et maximale entre les deux courbes est de 1.8% et 4.9% respectivement. D’un 
point de vue clinique, on s’intéresse plus particulièrement à la valeur de la vitesse maximale sur l’axe 
longitudinal. Elle est localisée à l’abscisse z=-1.36 mm et z=-1.91 mm respectivement dans le modèle 
théorique et numérique. La zone d’intérêt clinique s’étend entre les abscisses z=-2.4mm (bord d’attaque du 
clapet) et z=0mm (en aval du maximum de vitesse). 
 
FIG. 2 – Comparaison de la composante Uz de la vitesse obtenue au travers le modèle numérique (-) et 
théorique (x) le long de l’axe longitudinal Oz du repère principal. A droite, un zoom est réalisé au niveau des 
clapets. La zone délimitée par les pointillés correspond au jet central tranvalvulaire. 
Sont représentés sur la figure 3, les profils de vit sse longitudinaux (à gauche) et transversaux (à droite) 
obtenus sur différents diamètres compris dans la zone d’intérêt clinique. Une bonne concordance est 
observée entre les différentes courbes. La différence moyenne enregistrée sur la composante Uz de la vitesse 
est de 2.3%, 3.0% et 2.7% pour les abscisses z=-2.4mm, -1.36mm et 0mm respectivement. 
La composante Ux de la vitesse n’est pas de premier intérêt dans l’évaluation clinique de la performance 
valvulaire. Néanmoins, la comparaison proposée ici montre une nouvelle fois des profils similaires. L’entrée 
des clapets crée un convergent immédiatement suivi par un divergent. Les vecteurs vitesse sont donc dirigés 
vers l’axe central avant de longer les clapets. Ceci explique le changement de signe observé de part et d’autre 
de l’axe central, entre les profils Ux à z=-2.4mm et ceux plus en aval. 
19ème Congrès Français de Mécanique                                                                           Marseille, 24-28 août 2009 
  4 
Ces résultats valident ainsi le modèle théorique et montrent, comme déjà constaté expérimentalement dans la 
littérature [3], qu’un modèle 2D est suffisant pour rendre compte des écoulements trans- et post-valvulaires 
dans le cas d’une prothèse biclapet. 
 
FIG. 3 – Comparaison des profils de vitesses longitudinaux (à gauche) et transversaux (à droite) obtenus 
dans la zone d’intérêt clinique, par le modèle théorique (-) et numérique (x) 
4 Discussion 
L’évaluation de la performance valvulaire de la valve aortique native ou prothétique est réalisée par le calcul 
d’un index appelé EOA (Effective Orifice Area), obtenu lors d’un examen Doppler [2]. Ce calcul revient à 
estimer l’aire de la vena contracta observée en aval de tout rétrécissement de section. En routine clinique, 
l’EOA est calculé à partir de l’équation de continuité : 
 
Où ALV et VLV sont respectivement l’aire et la vitesse enregistrée dans la voie de chasse du ventricule 
gauche. Vjet correspond à la vitesse maximale enregistrée sur le faisceau Doppler balayant l’axe longitudinal 
19ème Congrès Français de Mécanique                                                                           Marseille, 24-28 août 2009 
  5 
central. Bien que cette définition soit tout à fait justifiable dans le cas d’une valve native, elle soulève des 
interrogations dans le cas d’une prothèse présentant plus qu’un orifice, comme c’est le cas pour la prothèse 
biclapet. Ce calcul n’a plus aucun sens physique puisqu’un jet central mais également deux jets latéraux sont 
caractéristiques des écoulements trans- et post-valvulaire d’une telle prothèse. D’un point de vue de la 
mécanique des fluides, l’équation de continuité ne peut donc plus être appliquée en l’état. Le calcul de 
l’EOA est faussé et peut donc être la cause de mauvaises interprétations cliniques. 
Cette observation est à relier aux différences d’EOA relevées dans la littérature et peuvent être expliquées 
par le modèle théorique proposé. Considérons par exemple les EOAs mesurées pour une prothèse 
commerciale de type St Jude (P1) et de type On-X (P2). La principale différence entre ces deux valves 
biclapets et l’angle maximal d’ouverture existant er  le clapet et l’axe central longitudinal : il est de 5° 
pour P1 et proche de zéro pour P2. La figure 4 présente le profil de vitesse longitudinal dont l’amplitude au 
centre est maximale pour deux prothèses : P1, avec un angle de 5° et P2, avec un angle de 0.5°. Si les autres 
données géométriques restent identiques, l’intensité de la circulation est plus faible dans le cas de P2 que de 
P1 entraînant par conséquent des vitesses induites plus faibles. L’EOA calculé pour P2 est donc plus grand 
que celui calculé pour P1. Or le fonctionnement, ie. la performance valvulaire, de ces deux prothèses e t 
optimal puisque les deux clapets de chacune de ces deux prothèses sont supposés avoir atteint leur angle 
d’ouverture maximale. On montre ainsi qu’une variation de 5 degrés sur l’angle d’inclinaison correspond à 
une variation relative de 25% sur la valeur de 1/Vjet et par conséquent sur la valeur de l’EOA. 
 
FIG. 4 – Profils de vitesses longitudinaux pour deux prothèses P1 (angle α= 5°) et P2 (angle α= 0.5°). A 
droite, variation de la valeur 1/Vjet en fonction de α. 
Par ailleurs, supposant que la position des pivots reste la même pour ces deux prothèses, l’espace entre l s 
bords d’attaque de chacun des deux clapets s’accroît dans le cas de P2, diminuant la vitesse induite dans le 
jet central. L’EOA et le design initial de la prothèse sont donc étroitement liés. Ainsi, comparer les valeurs 
optimales de l’EOA obtenus pour différentes prothèses commerciales revient à comparer les géométries de 
ces prothèses plutôt que d’établir un diagnostic de performance de la prothèse in situ.  
Le modèle théorique nous permet également de considérer des cas pathologiques comme le cas d’un 
dysfonctionnement de l’un des clapets. L’angle d’inclinaison de L1 est fixe et vaut α1=5°. Celui de L2 varie 
de α2=5° à 20°. La figure 5 présente les profils de vitesse longitudinaux obtenus. 
19ème Congrès Français de Mécanique                                                                           Marseille, 24-28 août 2009 
  6 
 
FIG. 5 – Profils de vitesses longitudinaux pour unemême prothèse mais dont les angles d’inclinaisons 
varient, simulant un dysfonctionnement de la prothèse 
Les profils de vitesse se dissymétrisent rapidement au fur et à mesure que la configuration géométrique elle-
même se dissymétrise. Le jet central est quasiment plateau dans le cas α1=α2 : une différence de 0.1m/s est 
observée entre la vitesse sur l’axe et la vitesse en proche parois de L1 et L2. Lorsque α2=20°, des valeurs 
élevées de la vitesse en proche paroi de L2 sont enregistrées, dissymétrisant le jet central. En supposant que 
l’on veuille continuer d’estimer le débit central, il devient alors difficile de faire l’hypothèse d’un profil 
plateau pour désigner le débit comme le produit de l’EOA et de la vitesse Vjet. Néanmoins, les différents 
profils obtenus soulignent l’intérêt d’une mesure Doppler des vitesses de part et d’autre des clapets pour 
rapidement détecter un dysfonctionnement de la prothèse et aider à l’interprétation de la valeur de l’EOA 
alors calculée. 
5 Conclusion 
Un modèle théorique basé sur la théorie des profils minces et des écoulements potentiels bidimensionnels a 
été proposé afin de rendre compte de l’écoulement à l’entrée d’une prothèse commerciale de type biclapet. 
Le très bon rapport entre le temps de calcul (inférieur à la minute) et la précision des résultats est le principal 
avantage de ce modèle en comparaison aux modèles numériq es de type interaction fluide-structure de la 
littérature. De par le paramétrage réalisé, ce modèle peut s’appliquer à de nombreux cas de figures cliniques. 
En mettant ainsi en évidence plusieurs problèmes liés à l’évaluation correcte de la performance valvulire 
prothétique, ce modèle est une base de travail en vu  de l’amélioration de ce diagnostic. 
References  
[1] Milne-Thomson L.M., Theoretical aerodynamics, 2nd edition, 1952. 
[2] Oh J.K., The echo manual, 2nd edition, 1999 
[3] Kaminsky R., PIV validation of blood-heart valve leaflet interaction modelling, International Journal of 
Artificial Organs, 30, 640-648, 2007. 
 


Modèle d’écoulement autour d’une prothèse valvulaire pour le diagnostic

            clinique

http://documents.irevues.inist.fr/bitstream/2042/36652/1/571.pdf

Modèle d’écoulement autour d’une prothèse valvulaire pour le diagnostic clinique

Abstract

Similar works

Full text

Available Versions

HAL AMU

I-Revues