Many cardiovascular diseases lead to local increases in relative pressure, reflecting the higher costs of driving blood flow. The utility of this biomarker for stratifying the severity of disease has thus driven the development of methods to measure these relative pressures. While intravascular catheterisation remains the most direct measure, its invasiveness limits clinical application in many instances. Non-invasive Doppler ultrasound estimates have partially addressed this gap; however only provide relative pressure estimates for a range of constricted cardiovascular conditions. Here we introduce a non-invasive method that enables arbitrary interrogation of relative pressures throughout an imaged vascular structure, leveraging modern phase contrast magnetic resonance imaging, the virtual work-energy equations, and a virtual field to provide robust and accurate estimates. The versatility and accuracy of the method is verified in a set of complex patient-specific cardiovascular models, where relative pressures into previously inaccessible flow regions are assessed. The method is further validated within a cohort of congenital heart disease patients, providing a novel tool for probing relative pressures in-vivo.

No abstract

Fluid-structure interaction (FSI) problems are pervasive in the computational engineering community. The need to address challenging FSI problems has led to the development of a broad range of numerical methods addressing a variety of application-specific demands. While a range of numerical and experimental benchmarks are present in the literature, few solutions are available that enable both verification and spatiotemporal convergence analysis. In this paper, we introduce a class of analytic solutions to FSI problems involving shear in channels and pipes. Comprised of 16 separate analytic solutions, our approach is permuted to enable progressive verification and analysis of FSI methods and implementations, in two and three dimensions, for static and transient scenarios as well as for linear and hyperelastic solid materials. Results are shown for a range of analytic models exhibiting progressively complex behavior. The utility of these solutions for analysis of convergence behavior is further demonstrated using a previously published monolithic FSI technique. The resulting class of analytic solutions addresses a core challenge in the development of novel FSI algorithms and implementations, providing a progressive testbed for verification and detailed convergence analysis.

Background Atrial fibrillation (AF) is responsible for almost one third of all strokes, with the left atrial appendage (LAA) being the primary thromboembolic source due to localised stimulation of prothrombotic mechanisms; blood stasis, hypercoagulability and endothelial damage, known as Virchow's triad. Aim We propose an in-silico modelling pipeline that leverages clinical imaging data to mechanistically assess patient thrombogenicity for all aspects of Virchow's triad to improve the prediction and prevention of AF-related stroke. Methods Two AF patients undergoing Cine magnetic resonance imaging (sinus rhythm (SR) N=1 or AF N=1 during imaging) were selected for 3D left atrial (LA) modelling with patient-specific myocardial deformation prescribed from image-derived wall motion. Blood stasis was quantified by computational fluid dynamics (CFD) simulations of 5 cardiac cycles [1]. Generation of three key coagulation proteins; thrombin, fibrinogen and fibrin, were modelled to represent thrombus growth and hypercoagulability [2]. Regions prone to thrombogenesis by endothelial damage were identified by the oscillatory shear index (OSI), time averaged wall shear stress (TAWSS) and endothelial cell activation potential (ECAP) metrics in the LAA [3]. Results Patient-specific LA simulations enabled the assessment of differences between SR and AF conditions, quantified as numerical characteristics of each aspect of Virchow's triad. In SR, blood flow velocities were in the range 0–2.6 m/s with mean of 0.85 m/s in the LA cavity, while AF had a range between 0–1.6 m/s with mean of 0.55 m/s. The peak and mean LAA velocities in SR were 0.85 m/s and 0.14 m/s, while AF had a peak LAA velocity of 0.32 m/s and mean of 0.09 m/s, showing a 38% decrease during AF. The thrombin concentration reached its steady state at 1.26 mmol/m3 in the AF case after 4.7 seconds, while thrombin was washed away from the initial injury site in SR. After 5 cardiac cycles of thrombus growth dynamics, the peak fibrin concentration in the LAA was 1.3 mmol/m3 in SR and 3.8 mmol/m3 in AF, with the thrombus area in AF being 40% larger. Fibrinogen concentration decreased at a rate equal to fibrin generation in both SR and AF solely in the area of thrombus formation. ECAP in the LAA had peak values of 2.9 in SR and 3.7 in AF, with the location at highest risk of thrombogenesis above the LAA entrance. LAA OSI had an average value of 0.45 in AF versus 0.36 in SR, showing a 26% increase. Similarly, the TAWSS was 3.5x10–3 Pa on average over the LAA in AF compared to 1.4x10–3 Pa in SR. Conclusions Patient-specific LA models combining these three quantitative characteristics can be used to predict the higher thrombogenic risk in AF. After further validation, this novel approach for quantitative assessment of AF patient thrombogenicity based on modelling all factors in Virchow's triad can personalise and improve management of AF patients with a risk of stroke. Funding Acknowledgement Type of funding sources: Public grant(s) – National budget only. Main funding source(s): UK Engineering and Physical Sciences Research Council