Computational study of the hemodynamics in healthy and diseased pulmonary arteries

The purpose of this work is to provide a computational study of the hemodynamics in the pulmonary artery of two patients with a repaired congenital heart disease (Tetralogy of Fallot). The two patients, after living for years without a pulmonary valve, underwent a pulmonary valve replacement (PVR) operation. We adopt a 3D-0D geometric multiscale coupled fluid dynamics model of the heart circulation based on the Navier-Stokes equations for the 3D domain of the pulmonary artery, reconstructed thanks to the computed tomographies of the patients provided by the Niguarda Hospital of Milan, and on electrical circuit analogies for the rest of the circulation. We propose a novel comparison among a healthy, a pre-operative and a post-operative scenario in order to assess the efficacy of the valve replacement operation in an image-based patient-specific context. Moreover, we perform a study on the flow repartition at the main pulmonary bifurcation based on three different approaches for the imposition of boundary conditions. The method has been implemented in the Finite Element C++ library LifeX.

Lo scopo di questo lavoro è di fornire uno studio computazionale dell’emodinamica nell’arteria polmonare di due pazienti con una cardiopatia congenita corretta (la Tetralogia di Fallot). I due pazienti, dopo aver vissuto per anni senza la valvola polmonare, sono stati sottoposti ad un intervento di sostituzione della stessa. Adottiamo un modello 3D-0D geometrico multiscala e accoppiato della circolazione cardiaca basato sulle equazioni di Navier-Stokes per il dominio 3D dell’arteria polmonare, ricostruita grazie a delle TAC fornite dall’ospedale di Niguarda di Milano, e su delle analogie elettriche per il resto della circolazione. Proponiamo un innovativo paragone tra uno scenario sano, uno preoperatorio e uno postoperatorio per analizzare l’efficacia dell’intervento di sostituzione della valvola in un contesto image-based specifico sul paziente. Effettuiamo inoltre uno studio sulla ripartizione del flusso alla biforcazione principale dell’arteria polmonare basandoci su tre diversi approcci per l’imposizione delle condizioni al bordo del problema. Il metodo è stato implementato nella libreria C++ a elementi finiti LifeX.