A Stabilized Mixed Hybrid Finite Element Framework for Large-Deformation Poromechanics with Torsion: Application to Transmural Perfusion and Vascular Hemodynamics

Le cœur humain est une pompe qui fait bien plus que simplement presser le sang ; il se tord, s'étire et se remodèle à chaque battement. À l'intérieur de la paroi musculaire épaisse du ventricule gauche, un problème complexe de dynamique des fluides se joue. Le muscle cardiaque n'est pas un bloc de tissu solide, mais un matériau semblable à une éponge, rempli à environ quinze pour cent de fluide qui circule à travers de minuscules canaux. Lorsque le cœur se remplit de sang puis se tord pour l'éjecter, ce fluide interne est comprimé et poussé, créant une pression qui peut soit aider, soit entraver le flux sanguin à travers les propres petits vaisseaux du cœur. Si la pression à l'intérieur du muscle devient trop élevée, elle peut écraser ces vaisseaux délicats, privant le tissu cardiaque d'oxygène. Comprendre comment le mouvement de torsion du cœur interagit avec ce flux de fluide interne est crucial, pourtant les modèles informatiques traditionnels peinent souvent à le simuler avec précision, produisant fréquemment des erreurs confuses ou se "bloquant" lorsque le tissu est presque impossible à comprimer.

Un chercheur de l'Université de technologie d'Eindhoven a développé une nouvelle façon de simuler ce comportement, créant un cadre mathématique qui capture la torsion et l'étirement à grande échelle de la paroi cardiaque sans tomber dans ces pièges courants. En traitant le muscle cardiaque comme une éponge saturée de fluide et en utilisant une méthode spécialisée qui suit séparément le mouvement du tissu solide et l'écoulement du fluide, l'étude révèle comment le mouvement de torsion naturel du cœur agit comme un mécanisme protecteur. Les simulations montrent que lorsque le cœur est autorisé à se tordre librement pendant le remplissage, cela réduit considérablement le stress sur les fibres de la paroi interne et abaisse la pression interne qui menace de faire s'effondrer les vaisseaux sanguins. Inversement, lorsque la paroi du cœur devient anormalement épaisse en raison d'une tension artériuse élevée, le fluide se retrouve piégé et la pression augmente dangereusement, une condition que les modèles informatiques standards ne parviennent souvent pas à prédire correctement.

Le cœur de ce travail réside dans la résolution d'un problème de longue date concernant la modélisation des matériaux mous et humides comme le muscle cardiaque. Dans les simulations informatiques standard, lorsqu'un matériau est presque incompressible — signifiant qu'il ne peut pas être écrasé en un volume plus petit — les mathématiques échouent souvent, créant des pics de pression irréalistes ou des résultats oscillants et saccadés qui ne reflètent pas la réalité. C'est ce qu'on appelle le verrouillage volumétrique. Pour résoudre cela, le chercheur a formulé une approche « hybride mixte ». Au lieu d'essayer de tout calculer en même temps, cette méthode décompose le problème en morceaux plus petits, imposant des règles strictes qui garantissent que le fluide circule de manière continue à travers les limites de chaque morceau, tout comme il le doit dans le corps réel. Cette technique permet à l'ordinateur de résoudre séparément le mouvement du tissu et la pression aux bords des blocs de simulation, garantissant que le flux de fluide est conservé exactement et que les résultats restent fluides et physiquement réalistes, même sous une déformation extrême.

Pour tester ce nouveau cadre, le chercheur a d'abord dérivé une solution mathématique précise pour un cylindre à paroi épaisse qui est gonflé et tordu, servant de référence parfaite. Il a ensuite appliqué la méthode à un modèle de ventricule gauche, une section cylindrique de muscle cardiaque dont les dimensions correspondent à l'anatomie humaine. Les simulations ont exploré trois scénarios spécifiques. Premièrement, ils ont examiné comment différents niveaux de pression sanguine à l'intérieur des propres vaisseaux coronaires du cœur affectent la rigidité de la paroi. Ils ont découvert que lorsque les vaisseaux sont bien remplis de sang, la pression du fluide interne augmente, ce qui rigidifie la paroi cardiaque au repos. Ce phénomène, connu sous le nom d'effet « tuyau d'arrosage », signifie que le cœur nécessite plus de pression pour se remplir, limitant ainsi le volume de sang qu'il peut absorber. Le nouveau modèle a prédit avec précision ce changement, montrant qu'une pression de perfusion interne élevée réduit la contrainte solide dans la paroi musculaire de dix kilopascals à trois kilopascals, partageant ainsi la charge entre le fluide et le tissu.

Le deuxième scénario a examiné le rôle vital du mouvement de torsion du cœur. Dans un cœur sain, les fibres musculaires sont disposées en spirale, allant de la surface interne à la surface externe, changeant d'angle au fur et à mesure. Lorsque le cœur se remplit, cette architecture en spirale provoque la torsion de la paroi. Les simulations ont comparé un cœur autorisé à se tordre librement à un cœur maintenu rigide. Lorsque la torsion était empêchée, le stress sur les fibres musculaires à la paroi interne grimpait en flèche à 57,20 kilopascals, une concentration de force dangereuse. Cependant, lorsque le modèle était autorisé à se tordre naturellement, ce pic de contrainte chutait de 45,5 % pour atteindre 31,20 kilopascals. Cette redistribution de la force aidait également le fluide à s'écouler plus facilement à travers la paroi, abaissant la pression des pores internes d'environ 128 millimètres de mercure à 122 millimètres de mercure. L'étude démontre que la torsion du cœur n'est pas seulement un sous-produit de la contraction, mais un mécanisme nécessaire pour relâcher le stress et maintenir le flux sanguin vers les couches internes du muscle.

Le troisième et plus critique cas d'étude a traité le cœur soumis à la contrainte d'une hypertension sévère, où la paroi s'épaissit considérablement pour résister à des pressions de pompage élevées. Dans ce scénario, le rayon interne du cœur a été modélisé à 2,00 millimètres et le rayon externe à 5,00 millimètres, avec une pression de cavité de 180 millimètres de mercure. Dans ces conditions, la paroi épaissie emprisonne le fluide près de la surface interne car le fluide doit parcourir une plus longue distance pour s'échapper. Le nouveau modèle hybride mixte a montré que ce fluide piégé crée une zone de haute pression soutenue qui persiste même après la charge initiale. En revanche, lorsque le chercheur a exécuté le même scénario en utilisant des techniques de modélisation informatique standard, les résultats étaient chaotiques. La méthode standard produisait des oscillations de pression non physiques et spécieuses près de la paroi interne et ne parvenait pas à capturer la véritable distribution des contraintes, entraînant un taux d'erreur massif de 12,75 %. Le nouveau modèle, cependant, a produit une courbe de pression lisse et monotone avec une erreur de seulement 0,01 %, identifiant correctement que la paroi épaisse combinée au mouvement de torsion crée un obstacle mécanique qui ralentit le drainage et soutient une pression interne élevée.

Ces conclusions suggèrent que la capacité du cœur à se tordre et à se détordre est une partie fondamentale de sa santé hydraulique. Lorsque le cœur se tord librement, il homogénéise la contrainte à travers la paroi et facilite le drainage du fluide, protégeant les microvaisseaux de l'effondrement. Lorsque la paroi s'épaissit ou que la torsion est restreinte, le fluide se retrouve piégé, entraînant des pressions élevées qui peuvent compromettre le flux sanguin. L'étude confirme que le nouveau cadre computationnel est capable de capturer ces interactions subtiles mais vitales, offrant un outil fiable pour comprendre comment les forces mécaniques et le flux de fluide sont couplés dans le cœur. Bien que le travail actuel se concentre sur une géométrie cylindrique simplifiée pour prouver le concept, le succès de cette méthode à gérer de grandes déformations et des interactions complexes entre fluide et structure pointe vers un avenir où de tels modèles pourront être appliqués à des formes cardiaques tridimensionnelles plus complexes pour mieux comprendre l'insuffisance cardiaque et les effets de l'hypertension artérielle.