A parametric study of the effect of 3D plaque shape on local hemodynamics and implications for plaque instability

The vast majority of heart attacks occur when vulnerable plaques rupture, releasing their lipid content into the blood stream leading to thrombus formation and blockage of a coronary artery. Detection of these unstable plaques before they rupture remains a challenge. Hemodynamic features including wall shear stress (WSS) and wall shear stress gradient (WSSG) near the vulnerable plaque and local inflammation are known to affect plaque instability. In this work, a computational workflow has been developed to enable a comprehensive parametric study detailing the effects of 3D plaque shape on local hemodynamics and their implications for plaque instability. Parameterized geometric 3D plaque models are created within a patient-specific coronary artery tree using a NURBS (non-uniform rational B-splines)-based vascular modeling pipeline. Realistic blood flow features are simulated by using a Navier-Stokes solver within an isogeometric finite-element analysis framework. Near wall hemodynamic quantities such as WSS and WSSG are quantified, and vascular distribution of an inflammatory marker (VCAM-1) is estimated. Results show that proximally skewed eccentric plaques have the most vulnerable combination of high WSS and high positive spatial WSSG, and the presence of multiple lesions increases risk of rupture. The computational tool developed in this work, in conjunction with clinical data, -could help identify surrogate markers of plaque instability, potentially leading to a noninvasive clinical procedure for the detection of vulnerable plaques before rupture.

Exergy destruction in atrial fibrillation and a new 'Exergy Age Index'

The concept of exergy in living organisms has been widely used to explore correlations between exergy and different physiological conditions. Atrial fibrillation (AF) is an abnormal physiological condition that takes place inside the heart and is recognised as a common supraventricular arrhythmia. AF can significantly undermine heart function and subsequently circulatory system. Thus, exergy analysis of cardiac flow during AF is a procedure to quantify the long-term impact of persistent AF. The present study adopts the lumped modelling approach for considering cardiovascular circulation and thermoregulation of the body to evaluate the exergy consumption and destruction of the heart in AF. In order to assess the impact of AF, four common AF-associated characteristics including lack of atrial kick, left atrial remodelling, left ventricular systolic dysfunction, and high-frequency fibrillation are examined. The results show that among AF deficiencies, high-frequency fibrillation is the main cause of exergy destruction of the heart during AF. Moreover, a novel 'exergy age index' is proposed which has shown that high fibrillatory conditions in AF can significantly accelerate the heart ageing process, which in turn substantiates AF's adverse impact on the heart.

Multi-Objective Optimisation of a Novel Bypass Graft with a Spiral Ridge

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction of a helical component in the flow field, either by means of out-of-plane helicity graft geometry or a spiral ridge. While the latter has been found to lack in performance when compared to the out-of-plane helicity designs, recent findings support the idea that the existing spiral ridge grafts can be further improved in performance through optimising relevant design parameters. In the current study, robust multi-objective optimisation techniques are implemented, covering a wide range of possible designs coupled with proven and well validated computational fluid dynamics (CFD) algorithms. It is shown that the final set of suggested design parameters could significantly improve haemodynamic performance and therefore could be used to enhance the design of spiral ridge bypass grafts.

Challenges in Modeling Hemodynamics in Cerebral Aneurysms Related to Arteriovenous Malformations

PURPOSE

The significantly higher incidence of aneurysms in patients with arteriovenous malformations (AVMs) suggests a strong hemodynamic relationship between these lesions. The presence of an AVM alters hemodynamics in proximal vessels by drastically changing the distal resistance, thus affecting intra-aneurysmal flow. This study discusses the challenges associated with patient-specific modeling of aneurysms in the presence of AVMs.

METHODS

We explore how the presence of a generic distal AVM affects upstream aneurysms by examining the relationship between distal resistance and aneurysmal wall shear stress using physiologically realistic estimates for the influence of the AVM on hemodynamics. Using image-based computational models of aneurysms and surrounding vasculature, aneurysmal wall-shear stress is calculated for a range of distal resistances corresponding to the presence of AVMs of various sizes and compared with a control case representing the absence of an AVM.

RESULTS

In the patient cases considered, the alteration in aneurysmal wall shear stress due to the presence of an AVM is considerable, as much as 19 times the base case wall shear stress. Furthermore, the relationship between aneurysmal wall shear stress and distal resistance is shown to be highly geometry-dependent and nonlinear. In most cases, the range of physiologically realistic possibilities for AVM-related distal resistance are so large that patient-specific flow measurements are necessary for meaningful predictions of wall shear stress.

CONCLUSIONS

The presented work offers insight on the impact of distal AVMs on aneurysmal wall shear stress using physiologically realistic computational models. Patient-specific modeling of hemodynamics in aneurysms and associated AVMs has great potential for understanding lesion pathogenesis, surgical planning, and assessing the effect of treatment of one lesion relative to another. However, we show that modeling approaches cannot usually meaningfully quantify the impact of AVMs if based solely on imaging data from CT and X-ray angiography, currently used in clinical practice. Based on recent studies, it appears that 4D flow MRI is one promising approach to obtaining meaningful patient-specific flow boundary conditions that improve modeling fidelity.

Image-based patient-specific flow simulations are consistent with stroke in pediatric cerebrovascular disease

Moyamoya disease (MMD) is characterized by narrowing of the distal internal carotid artery and the circle of Willis (CoW) and leads to recurring ischemic and hemorrhagic stroke. A retrospective review of data from 50 pediatric MMD patients revealed that among the 24 who had a unilateral stroke and were surgically treated, 11 (45.8%) had a subsequent, contralateral stroke. There is no reliable way to predict these events. After a pilot study in Acta-/- mice that have features of MMD, we hypothesized that local hemodynamics are predictive of contralateral strokes and sought to develop a patient-specific analysis framework to noninvasively assess this stroke risk. A pediatric MMD patient with an occlusion in the right middle cerebral artery and a right-sided stroke, who was surgically treated and then had a contralateral stroke, was selected for analysis. By using an unsteady Navier-Stokes solver within an isogeometric analysis framework, blood flow was simulated in the CoW model reconstructed from the patient's postoperative imaging data, and the results were compared with those from an age- and sex-matched control subject. A wall shear rate (WSR) > 60,000 s-1 (about 12 × higher than the coagulation threshold of 5000 s-1 and 9 × higher than control) was measured in the terminal left supraclinoid artery; its location coincided with that of the subsequent postsurgical left-sided stroke. A parametric study of disease progression revealed a strong correlation between the degree of vascular morphology altered by MMD and local hemodynamic environment. The results suggest that an occlusion in the CoW could lead to excessive contralateral WSRs, resulting in thromboembolic ischemic events, and that WSR could be a predictor of future stroke.

Hemodynamic Abnormalities in the Aorta of Turner Syndrome Girls

Congenital abnormalities in girls and women with Turner syndrome (TS), alongside an underlying predisposition to obesity and hypertension, contribute to an increased risk of cardiovascular disease and ultimately reduced life expectancy. We observe that children with TS present a greater variance in aortic arch morphology than their healthy counterparts, and hypothesize that their hemodynamics is also different. In this study, computational fluid dynamic (CFD) simulations were performed for four TS girls, and three age-matched healthy girls, using patient-specific inlet boundary conditions, obtained from phase-contrast MRI data. The visualization of multidirectional blood flow revealed an increase in vortical flow in the arch, supra-aortic vessels, and descending aorta, and a correlation between the presence of aortic abnormalities and disturbed flow. Compared to the relatively homogeneous pattern of time-averaged wall shear stress (TAWSS) on the healthy aortae, a highly heterogeneous distribution with elevated TAWSS values was observed in the TS geometries. Visualization of further shear stress parameters, such as oscillatory shear index (OSI), normalized relative residence time (RRTn), and transverse WSS (transWSS), revealed dissimilar heterogeneity in the oscillatory and multidirectional nature of the aortic flow. Taking into account the young age of our TS cohort (average age 13 ± 2 years) and their obesity level (75% were obese or overweight), which is believed to accelerate the initiation and progression of endothelial dysfunction, these findings may be an indication of atherosclerotic disease manifesting earlier in life in TS patients. Age, obesity and aortic morphology may, therefore, play a key role in assessing cardiovascular risk in TS children.

Effects of Ageing on Aortic Circulation During Atrial Fibrillation; a Numerical Study on Different Aortic Morphologies

Atrial fibrillation (AF) can alter intra-cardiac flow and cardiac output that subsequently affects aortic flow circulation. These changes may become more significant where they occur concomitantly with ageing. Aortic ageing is accompanied with morphological changes such as dilation, lengthening, and arch unfolding. While the recognition of AF mechanism has been the subject of numerous studies, less focus has been devoted to the aortic circulation during the AF and there is a lack of such investigation at different ages. The current work aims to address the present gap. First, we analyse aortic flow distribution in three configurations, which attribute to young, middle and old people, using geometries constructed via clinical data. We then introduce two transient inlet flow conditions representative of key AF-associated defects. Results demonstrate that both AF and ageing negatively affect flow circulation. The main consequence of concomitant occurrence is enhancement of endothelial cell activation potential (ECAP) throughout the vascular domain, mainly at aortic arch and descending thoracic aorta, which is consistent with some clinical observations. The outcome of the current study suggests that AF exacerbates the vascular defects occurred due to the ageing, which increases the possibility of cardiovascular diseases per se.

A computational simulation platform for designing real-time monitoring systems with application to COVID-19

With the aim of contributing to the fight against the coronavirus disease 2019 (COVID-19), numerous strategies have been proposed. While developing an effective vaccine can take months up to years, detection of infected patients seems like one of the best ideas for controlling the situation. The role of biosensors in containing highly pathogenic viruses, saving lives and economy is evident. A new competitive numerical platform specifically for designing microfluidic-integrated biosensors is developed and presented in this work. Properties of the biosensor, sample, buffer fluid and even the microfluidic channel can be modified in this model. This feature provides the scientific community with the ability to design a specific biosensor for requested point-of-care (POC) applications. First, the validation of the presented numerical platform against experimental data and then results and discussion, highlighting the important role of the design parameters on the performance of the biosensor is presented. For the latter, the baseline case has been set on the previous studies on the biosensors suitable for SARS-CoV, which has the highest similarity to the 2019 nCoV. Subsequently, the effects of concentration of the targeted molecules in the sample, installation position and properties of the biosensor on its performance were investigated in 11 case studies. The presented numerical framework provides an insight into understanding of the virus reaction in the design process of the biosensor and enhances our preparation for any future outbreaks. Furthermore, the integration of biosensors with different devices for accelerating the process of defeating the pandemic is proposed.

Numerical framework for simulating bio-species transport in microfluidic channels with application to antibody biosensors

Diagnosis is a fundamental stage in health care and medical treatment. Microfluidic biosensors and lab-on-a-chip devices are amongst the few practical tools for achieving this goal. A new computational code, specifically for designing microfluidic-integrated biosensors is developed, the details of which is presented in this work. This new approach is developed using control-volume based finite-element (CVFEM) method and solves bio-recognition chemical reactions and full Navier–Stokes equations. The results of the proposed platform are validated against the experimental data for a microfluidic based biosensor, where excellent agreement is achieved. The properties of the biosensor, sample, buffer fluid and even the microfluidic channel can easily be modified in this platform. This feature provides the scientific community with the ability to design a specific biosensor for requested point-of-care applications.•

A new approach is developed using control-volume based finite-element (CVFEM) method for investigating flow inside a microfluidic-integrated biosensor. It is also used to study the influence of surface functionalization on binding cycle.

•

The proposed model solves bio-recognition chemical reactions as well as full Navier–Stokes and energy equations. Experimental-based or personalized equations of the chemical reactions and flow behaviour are adoptable to this code.

•

The developed model is Fortran-based and has the potential to be used in both industry and academia for biosensing technology.

Blood flow simulations in the pulmonary bifurcation in relation to adult patients with repaired tetralogy of Fallot

Understanding the haemodynamic environment of the pulmonary bifurcation is important in adults with repaired conotruncal congenital heart disease. In these patients, dysfunction of the pulmonary valve and narrowing of the branch pulmonary arteries are common and can have serious clinical consequences. The aim of this study was to numerically investigate the underlying blood flow characteristics in the pulmonary trunk under a range of simplified conditions. For that, an in-depth analysis was conducted in idealised two-dimensional geometries that facilitate parametric investigation of healthy and abnormal conditions. Subtle variations in morphology influenced the haemodynamic environment and wall shear stress distribution. The pressure in the left pulmonary artery was generally higher than that in the right and main arteries, but was markedly reduced in the presence of a local stenosis. Different downstream pressure conditions altered the branch flow ratio, from 50:50% to more realistic 60:40% ratios in the right and left pulmonary artery, respectively. Despite some simplifications, this study highlights some previously undocumented aspects of the flow in bifurcating geometries, by clarifying the role of the stagnation point location on wall shear stress and differential branch pressures. In addition, measurements of the mean pressure ratios in the pulmonary bifurcation are discussed in the context of a new haemodynamic index which could potentially contribute to the assessment of left pulmonary artery stenosis in tetralogy of Fallot patients. Further studies are required to confirm the results in patient-specific models with personalised physiological flow conditions.

Numerical Study of Atrial Fibrillation Effects on Flow Distribution in Aortic Circulation

Atrial fibrillation (AF) is the most common type of arrhythmia, which undermines cardiac function. Atrial fibrillation is a multi-facet malady and it may occur as a result of other diseases or it may trigger other problems. One of the main complications of AF is stroke due to the possibility of clot formation inside the atrium. However, the possibility of stroke occurrence due to the AF and the location from which an embolus dispatches are subject of debate. Another hypothesis about the embolus formation during AF is thrombus formation in aorta and carotid arteries, embolus detachment and its movement. To investigate the possibility of the latter postulation, the current work suggests a parametric study to quantify the sensitivity of aortic flow to four common AF traits including lack of atrial kick, atrial remodelling, left ventricle systolic dysfunction, and high frequency fibrillation. The simulation was carried out by coupling several in-house codes and ANSYS-CFX module. The results reveal that AF traits lower flow rate at left ventricular outflow tract, which in general lowers blood perfusion to systemic, cerebral and coronary circulations. Consequently, it leads to endothelial cell activation potential (ECAP) increase and variation of flow structure that both suggest predisposed areas to atherogenesis and thrombus formation in different regions in ascending aorta, aortic arch and descending thoracic aorta.

Optimisation of a Novel Spiral-Inducing Bypass Graft Using Computational Fluid Dynamics

Graft failure is currently a major concern for medical practitioners in treating Peripheral Vascular Disease (PVD) and Coronary Artery Disease (CAD). It is now widely accepted that unfavourable haemodynamic conditions play an essential role in the formation and development of intimal hyperplasia, which is the main cause of graft failure. This paper uses Computational Fluid Dynamics (CFD) to conduct a parametric study to enhance the design and performance of a novel prosthetic graft, which utilises internal ridge(s) to induce spiral flow. This design is primarily based on the identification of the blood flow as spiral in the whole arterial system and is believed to improve the graft longevity and patency rates at distal graft anastomoses. Four different design parameters were assessed in this work and the trailing edge orientation of the ridge was identified as the most important parameter to induce physiological swirling flow, while the height of the ridge also significantly contributed to the enhanced performance of this type of graft. Building on these conclusions, an enhanced configuration of spiral graft is proposed and compared against conventional and spiral grafts to reaffirm its potential benefits.