Modélisation mathématique de l’athérosclérose

L'athérosclérose est une maladie inflammatoire qui commence quand les lipoprotéines de faible densité (LDL) entrent dans l'intima du vaisseau sanguin où elles sont oxydées (ox-LDL). Le ox-LDL est considéré comme un agent dangereux par le système immunitaire provoquant ainsi une réponse immunitaire. Cette réponse immunitaire déclenche le recrutement des monocytes dans l'intima où elles se transforment en macrophages et ensuite en cellules spumeuses. Ce dernier amplifie la production des cytokines inflammatoires et davantage de recrutement des monocytes. Ce processus auto-amplifié est compensé par la sécrétion de cytokines anti-inflammatoires (anti-inflammation biochimique) et la migration des cellules musculaires lisses pour former une chape fibreuse qui couvre le noyau lipidique. Cette chape fibreuse avec le noyau lipidique s'appellent la plaque d'athérosclérose. Celle-ci change la géométrie du vaisseau sanguin en le rétrécissant et interagit avec du flux sanguin. Cette interaction peut avoir des conséquences dangereuses liées à la rupture de plaque ou à la formation du caillot de sang. La thèse est consacrée à la modélisation mathématique de ces phénomènes. Elle est composée de deux parties : Nous développons des modèles mathématiques basés sur des équations de réaction diffusion afin de décrire le processus inflammatoire. Le premier modèle est unidimensionnel. Il nous permet d'expliquer comment le développement de l'athérosclérose dépend de la concentration en cholestérol (ox-LDL). Si cette concentration dans l'intima est basse, alors la maladie ne se développera pas. Les concentrations intermédiaires de ox-LDL peuvent mener au développement de la maladie dans certaines conditions. Nous montrons que l'inflammation se propage en front d'ondes de réaction-diffusion. Les concentrations élevées de ox-LDL engendre le développement de la maladie. Même une petite perturbation du cas non inflammatoire mène à une propagation d'ondes qui correspond à l'inflammation. Ensuite nous étudions un modèle bidimensionnel qui représente un système d'équations type réaction-diffusion sur une bande. La deuxième dimension correspond à la section transversale de l'intima et une condition aux limites non-linéaire décrit le recrutement des monocytes. Cette condition aux limites est une fonction des concentrations des cytokines. Nous démontrons l'existence des fronts de propagation d'onde et confirmons les résultats précédents qui montrent que l'athérosclérose se développe en tant qu'onde de réaction-diffusion. Les résultats théoriques des deux modèles sont confirmés par des simulations numériques qui montrent que le cas bidimensionnel converge vers le cas unidimensionnel quand l'épaisseur de l'intima tend vers zéro. Une fois la plaque se forme, elle interagit avec le flux sanguin engendrant de différentes conséquences mécaniques et biochimiques. Nous développons un modèle d'interaction fluide-structure. La plaque d'athérome composée d'un dépôt lipidique couvert par une chape fibreuse, les deux étant modélisés en tant que matériaux hyper-élastiques. Le sang est considéré comme un fluide non-Newtonien avec une viscosité variable modélisée selon la loi de Carreau. Les paramètres utilisés dans nos simulations sont tirés de données expérimentales mentionnées dans la littérature. Nous étudions les effets non-Newtoniens sur les recirculations du sang en aval de la plaque d'athérome et aussi sur les contraintes sur celle-ci. Les simulations montrent que le modèle Newtonien surestime les recirculations de manière significative par rapport au modèle non-Newtonien. Elles montrent aussi que le modèle Newtonien sous-estime légèrement les contraintes sur la plaque pour des taux de cisaillement usuels, mais cette sous-estimation devient importante pour des taux de cisaillement bas.