Patient-specific assessment of cardiovascular function by combination of clinical data and computational model with applications to patients undergoing Fontan operation

The Influence of Aortic Valve Disease on Coronary Hemodynamics: A Computational Model-Based Study

Aortic valve disease (AVD) often coexists with coronary artery disease (CAD), but whether and how the two diseases are correlated remains poorly understood. In this study, a zero-three dimensional (0-3D) multi-scale modeling method was developed to integrate coronary artery hemodynamics, aortic valve dynamics, coronary flow autoregulation mechanism, and systemic hemodynamics into a unique model system, thereby yielding a mathematical tool for quantifying the influences of aortic valve stenosis (AS) and aortic valve regurgitation (AR) on hemodynamics in large coronary arteries. The model was applied to simulate blood flows in six patient-specific left anterior descending coronary arteries (LADs) under various aortic valve conditions (i.e., control (free of AVD), AS, and AR). Obtained results showed that the space-averaged oscillatory shear index (SA-OSI) was significantly higher under the AS condition but lower under the AR condition in comparison with the control condition. Relatively, the overall magnitude of wall shear stress was less affected by AVD. Further data analysis revealed that AS induced the increase in OSI in LADs mainly through its role in augmenting the low-frequency components of coronary flow waveform. These findings imply that AS might increase the risk or progression of CAD by deteriorating the hemodynamic environment in coronary arteries.

Model-based analysis of the sensitivities and diagnostic implications of FFR and CFR under various pathological conditions

Although fractional flow reserve (FFR) and coronary flow reserve (CFR) are both frequently used to assess the functional severity of coronary artery stenosis, discordant results of diagnosis between FFR and CFR in some patient cohorts have been reported. In the present study, a computational model was employed to quantify the impacts of various pathophysiological factors on FFR and CFR. In addition, a hyperemic myocardial ischemic index (HMIx) was proposed as a reference for comparing the diagnostic performances of FFR and CFR. Obtained results showed that CFR was more susceptible than FFR to the influence of many pathophysiological factors unrelated to coronary artery stenosis. In particular, the numerical study proved that increasing hyperemic coronary microvascular resistance significantly elevated FFR while reducing CFR despite fixed severity of coronary artery stenosis, whereas introducing aortic valve disease only caused a significant decrease in CFR with little influence on FFR. These results provided theoretical evidence for explaining some clinical observations, such as the increased risk of discordant diagnostic results between FFR and CFR in patients with increased hyperemic microvascular resistance, and significant increase in CFR after surgical relief of severe aortic valve disease. When evaluated with respect to the predictive value for hyperemic myocardial ischemia, the performance of FFR was found to be considerably compromised in the presence of severe coronary vasodilation dysfunction or aortic valve disease, whereas the relationship between CFR and HMIx remained relatively stable, suggesting that CFR may be a more reliable indicator of myocardial ischemia under complex pathophysiological conditions.

Impact of Atrial Fibrillation on Fontan Circulation: Fontan Computational Model

BACKGROUND

Patients with Fontan circulation may develop heart failure resulting in atrial fibrillation during the late phase. Inotropic effects to ameliorate hemodynamics on the Fontan circulation are not well understood, especially when in atrial fibrillation. This study was performed to determine whether dobutamine therapy in patients with Fontan circulation has limited effects on improving hemodynamics.

METHODS

Lumped computational models (sinus and atrial fibrillation) were employed including biventricular, atriopulmonary connection (APC), and extracardiac total cavopulmonary connection (TCPC) Fontan models. The condition of atrial fibrillation including lack of atrial beat, irregular ventricular contraction, and time-varying elastance for the ventricle was introduced. A different dose of dobutamine was given by varying the elastance of the ventricle, heart rate, and peripheral resistance.

RESULTS

In all models, the cardiac output decreased by 22.5% to 25.8% in atrial fibrillation. At 10 μg/kg/min of dobutamine in sinus rhythm, the cardiac output increased by 32.3% in the biventricular model but by only 9.2% (p<0.001) and 9.1% (p<0.001) in the APC and TCPC Fontan models, respectively. At 10 μg/kg/min of dobutamine in atrial fibrillation, the percent increase in the cardiac output in the Fontan circulation [11.8% increase in APC (p<0.001) and 11.9% increase in TCPC (p<0.001)] was significantly less than that in the biventricular circulation (32.3% increase).

CONCLUSIONS

In the Fontan circulation, atrial fibrillation itself reduced the cardiac output by approximately 25%, and dobutamine had a limited effect on increasing the cardiac output, especially when in atrial fibrillation. Maintaining sinus rhythm in patients with Fontan circulation is very important.

Fluid-Structure Interaction Simulation of an Intra-Atrial Fontan Connection

Simple Summary

A fluid-structure interaction (FSI) simulation of an intra-atrial Fontan connection was performed. Power loss and pressure drop results fluctuated less during the FSI simulation than during the simulation run with rigid walls, but there were no observable differences in time-averaged pressure drop, connection power loss or hepatic flow distribution. These results suggested that employing a rigid wall is a reasonable assumption when evaluating time-averaged hemodynamic quantities of the Fontan connection under resting breath-held flow conditions.

Abstract

Total cavopulmonary connection (TCPC) hemodynamics has been hypothesized to be associated with long-term complications in single ventricle heart defect patients. Rigid wall assumption has been commonly used when evaluating TCPC hemodynamics using computational fluid dynamics (CFD) simulation. Previous study has evaluated impact of wall compliance on extra-cardiac TCPC hemodynamics using fluid-structure interaction (FSI) simulation. However, the impact of ignoring wall compliance on the presumably more compliant intra-atrial TCPC hemodynamics is not fully understood. To narrow this knowledge gap, this study aims to investigate impact of wall compliance on an intra-atrial TCPC hemodynamics. A patient-specific model of an intra-atrial TCPC is simulated with an FSI model. Patient-specific 3D TCPC anatomies were reconstructed from transverse cardiovascular magnetic resonance images. Patient-specific vessel flow rate from phase-contrast magnetic resonance imaging (MRI) at the Fontan pathway and the superior vena cava under resting condition were prescribed at the inlets. From the FSI simulation, the degree of wall deformation was compared with in vivo wall deformation from phase-contrast MRI data as validation of the FSI model. Then, TCPC flow structure, power loss and hepatic flow distribution (HFD) were compared between rigid wall and FSI simulation. There were differences in instantaneous pressure drop, power loss and HFD between rigid wall and FSI simulations, but no difference in the time-averaged quantities. The findings of this study support the use of a rigid wall assumption on evaluation of time-averaged intra-atrial TCPC hemodynamic metric under resting breath-held condition.

Hemodynamic Effects of A Simplified Venturi Conduit for Fontan Circulation: A Pilot, In Silico Analysis

Objectives: To study the effects of a self-powered Fontan circulation in both idealized Fontan models and patient-specific models. Methods: In silico, a conduit with a nozzle was introduced from ascending aorta into the anastomosis of superior vena cava and pulmonary artery. Computational fluid dynamics (CFD) simulation was applied to calculate the fluid fields of models. Three 3-dimentional idealized models with different offsets were reconstructed by computer-aided design to evaluate the effects of the self-powered conduit. Furthermore, to validate the effects in patient-specific models, the conduit was introduced to three reconstructed models with different offsets. Results: The pressures at superior venae cavae and inferior venae cavae were decreased in both idealized models (0.4 mmHg) and patient-specific models (0.7 mmHg). In idealized models, the flows to left lungs were decreased (70%) by the jets from the conduits. However, in patient-specific models, the reductions of blood to the left lungs were relatively limited (30%) comparing to idealized models. Conclusions: CFD simulation was applied to analyze the effectiveness of the Fontan self-powered conduit. This self-powered conduit may help to decrease the venae cavae pressures and increase the flow to pulmonary arteries.

Developing transmission line equations of oxygen transport for predicting oxygen distribution in the arterial system

The oxygen content in the arterial system plays a significant role in determining the physiological status of a human body. Understanding the oxygen concentration distribution in the arterial system is beneficial for the prevention and intervention of vascular disease. However, the oxygen concentration in the arteries could not be noninvasively monitored in clinical research. Although the oxygen concentration distribution in a vessel could be obtained from a three-dimensional (3D) numerical simulation of blood flow coupled with oxygen transport, a 3D numerical simulation of the systemic arterial tree is complicated and requires considerable computational resources and time. However, the lumped parameter model of oxygen transport derived from transmission line equations of oxygen transport requires fewer computational resources and less time to numerically predict the oxygen concentration distribution in the systemic arterial tree. In this study, transmission line equations of oxygen transport are developed according to the theory of oxygen transport in the vessel, and fluid transmission line equations are used as the theoretical reference for the development. The transmission line equations of oxygen transport could also be regarded as the theoretical basis for developing lumped parameter models of other substances in blood.

Determinant Factors for Arterial Hemodynamics in Hypertension: Theoretical Insights From a Computational Model-Based Study

Hypertension is a well-documented predictive factor for cardiovascular events. Clinical studies have extensively demonstrated the differential hemodynamic consequences of various antihypertensive drugs, but failed to clearly elucidate the underlying mechanisms due to the difficulty in performing a quantitative deterministic analysis based on clinical data that carry confounding information stemming from interpatient differences and the nonlinearity of cardiovascular hemodynamics. In the present study, a multiscale model of the cardiovascular system was developed to quantitatively investigate the relationships between hemodynamic variables and cardiovascular properties under hypertensive conditions, aiming to establish a theoretical basis for assisting in the interpretation of clinical observations or optimization of therapy. Results demonstrated that heart period, central arterial stiffness, and arteriolar radius were the major determinant factors for blood pressures and flow pulsatility indices both in large arteries and in the microcirculation. These factors differed in the degree and the way in which they affect hemodynamic variables due to their differential effects on wave reflections in the vascular system. In particular, it was found that the hemodynamic effects of varying arteriolar radius were considerably influenced by the state of central arterial stiffness, and vice versa, which implied the potential of optimizing antihypertensive treatment by selecting proper drugs based on patient-specific cardiovascular conditions. When analyzed in relation to clinical observations, the simulated results provided mechanistic explanations for the beneficial pressure-lowering effects of vasodilators as compared to β-blockers, and highlighted the significance of monitoring and normalizing arterial stiffness in the treatment of hypertension.

Impact of the location of the fenestration on Fontan circulation haemodynamics: a three-dimensional, computational model study

OBJECTIVES

There is no consensus or theoretical explanation regarding the optimal location for the fenestration during the Fontan operation. We investigated the impact of the location of the fenestration on Fontan haemodynamics using a three-dimensional Fontan model in various physiological conditions.

METHODS

A three-dimensional Fontan model was constructed on the basis of CT images, and a 4-mm-diameter fenestration was located between the extracardiac Fontan conduit and the right atrium at three positions: superior, middle, and inferior part of the conduit. Haemodynamics in the Fontan route were analysed using a three-dimensional computational fluid dynamic model in realistic physiological conditions, which were predicted using a lumped parameter model of the cardiovascular system. The respiratory effect of the caval flow was taken into account. The flow rate through the fenestration, the effect of lowering the central venous pressure, and wall shear stress in the Fontan circuit were evaluated under central venous pressures of 10, 15, and 20 mmHg. The pulse power index and pulsatile energy loss index were calculated as energy loss indices.

RESULTS

Under all central venous pressures, the middle-part fenestration demonstrated the most significant effect on enhancing the flow rate through the fenestration while lowering the central venous pressure. The middle-part fenestration produced the highest time-averaged wall shear stress, pressure pulse index, and pulsatile energy loss index.

CONCLUSIONS

Despite slightly elevated energy loss, the middle-part fenestration most significantly increased cardiac output and lowered central venous pressure under respiration in the Fontan circulation.

A computational study of the Fontan circulation with fenestration or hepatic vein exclusion

Fontan patients may undergo additional surgical modifications to mitigate complications like protein-losing enteropathy, liver cirrhosis, and other issues in their splanchnic circulation. Recent case reports show promise for several types of modifications, but the subtle effects of these surgeries on the circulation are not well understood. In this paper, we employ mathematical modeling of blood flow to systematically quantify the impact of these surgical changes on extracardiac Fontan hemodynamics. We investigate two modifications: (1) the fenestrated Fontan and (2) the Fontan with hepatic vein exclusion. Closed-loop hemodynamic models are used, which consist of one-dimensional networks for the major vessels and zero-dimensional models for the heart and organ beds. Numerical results suggest the hepatic vein exclusion has the greatest overall impact on the hemodynamics, followed by the largest sized fenestration. In particular, the hepatic vein exclusion drastically lowers portal venous pressure while the fenestration decreases pulmonary artery pressure. Both modifications increase flow to the intestines, a finding consistent with their utility in clinical practice for combating complications in the splanchnic circulation.

Creating virtual ARDS patients

This paper presents the methodology used in patient-specific calibration of a novel highly integrated model of the cardiovascular and pulmonary pathophysiology associated with Acute Respiratory Distress Syndrome (ARDS). We focus on data from previously published clinical trials on the static and dynamic cardio-pulmonary responses of three ARDS patients to changes in ventilator settings. From this data, the parameters of the integrated model were identified using an optimization-based methodology in multiple stages. Computational simulations confirm that the resulting model outputs accurately reproduce the available clinical data. Our results open up the possibility of creating in silico a biobank of virtual ARDS patients that could be used to evaluate current, and investigate novel, therapeutic strategies.

Patient-specific fracture risk assessment of vertebrae: A multiscale approach coupling X-ray physics and continuum micromechanics

While in clinical settings, bone mineral density measured by computed tomography (CT) remains the key indicator for bone fracture risk, there is an ongoing quest for more engineering mechanics-based approaches for safety analyses of the skeleton. This calls for determination of suitable material properties from respective CT data, where the traditional approach consists of regression analyses between attenuation-related grey values and mechanical properties. We here present a physics-oriented approach, considering that elasticity and strength of bone tissue originate from the material microstructure and the mechanical properties of its elementary components. Firstly, we reconstruct the linear relation between the clinically accessible grey values making up a CT, and the X-ray attenuation coefficients quantifying the intensity losses from which the image is actually reconstructed. Therefore, we combine X-ray attenuation averaging at different length scales and over different tissues, with recently identified 'universal' composition characteristics of the latter. This gives access to both the normally non-disclosed X-ray energy employed in the CT-device and to in vivo patient-specific and location-specific bone composition variables, such as voxel-specific mass density, as well as collagen and mineral contents. The latter feed an experimentally validated multiscale elastoplastic model based on the hierarchical organization of bone. Corresponding elasticity maps across the organ enter a finite element simulation of a typical load case, and the resulting stress states are increased in a proportional fashion, so as to check the safety against ultimate material failure. In the young patient investigated, even normal physiological loading is probable to already imply plastic events associated with the hydrated mineral crystals in the bone ultrastructure, while the safety factor against failure is still as high as five. Copyright © 2016 John Wiley & Sons, Ltd.

Computational modeling of blood flow steal phenomena caused by subclavian stenoses

The study of steal mechanisms caused by vessel obstructions is of the utmost importance to gain understanding about their pathophysiology, as well as to improve diagnosis and management procedures. The goal of this work is to perform a computational study to gain insight into the hemodynamic forces that drive blood flow steal mechanisms caused by subclavian artery stenosis. Such condition triggers a flow disorder known as subclavian steal. When this occurs in patients with internal thoracic artery anastomosed to the coronary vessels, the phenomenon includes a coronary-subclavian steal. True steal can exist in cases of increased arm blood flow, potentially resulting in neurological complications and, in the case of coronary-subclavian steal, graft function failure. In this context, the anatomically detailed arterial network (ADAN) model is employed to simulate subclavian steal and coronary-subclavian steal phenomena. Model results are verified by comparison with published data. It is concluded that this kind of model allows us to effectively address complex hemomdynamic phenomena occurring in clinical practice. More specifically, in the studied conditions it is observed that a regional brain steal occurs, primarily affecting the posterior circulation, not fully compensated by the anterior circulation. In the case of patients with coronary revascularization, it is concluded that there is a large variability in graft hemodynamic environments, which physically explain both the success of the procedure in cases of severe occlusive disease, and the reason for graft dysfunction in mildly stenosed left anterior descending coronary artery, due to alternating graft flow waveform signatures.

Comparison of the Windkessel model and structured-tree model applied to prescribe outflow boundary conditions for a one-dimensional arterial tree model

One-dimensional (1D) modeling is a widely adopted approach for studying wave propagation phenomena in the arterial system. Despite the frequent use of the Windkessel (WK) model to prescribe outflow boundary conditions for 1D arterial tree models, it remains unclear to what extent the inherent limitation of the WK model in describing wave propagation in distal vasculatures affect hemodynamic variables simulated at the arterial level. In the present study, a 1D model of the arterial tree was coupled respectively with a WK boundary model and a structured-tree (ST) boundary model, yielding two types of arterial tree models. The effective resistances, compliances and inductances of the WK and ST boundary models were matched to facilitate quantitative comparisons. Obtained results showed that pressure/flow waves simulated by the two models were comparable in the aorta, whereas, their discrepancies increased towards the periphery. Wave analysis revealed that the differences in reflected waves generated by the boundary models were the major sources of pressure wave discrepancies observed in large arteries. Additional simulations performed under aging conditions demonstrated that arterial stiffening with age enlarged the discrepancies, but with the effects being partly counteracted by physiological aortic dilatation with age. These findings suggest that the method adopted for modeling the outflow boundary conditions has considerable influence on the performance of a 1D arterial tree model, with the extent of influence varying with the properties of the arterial system.

Total Cavopulmonary Connection is Superior to Atriopulmonary Connection Fontan in Preventing Thrombus Formation: Computer Simulation of Flow-Related Blood Coagulation

The classical Fontan route, namely the atriopulmonary connection (APC), continues to be associated with a risk of thrombus formation in the atrium. A conversion to a total cavopulmonary connection (TCPC) from the APC can ameliorate hemodynamics for the failed Fontan; however, the impact of these surgical operations on thrombus formation remains elusive. This study elucidates the underlying mechanism of thrombus formation in the Fontan route by using a two-dimensional computer hemodynamic simulation based on a simple blood coagulation rule. Hemodynamics in the Fontan route was simulated with Navier-Stokes equations. The blood coagulation and the hemodynamics were combined using a particle method. Three models were created: APC with a square atrium, APC with a round atrium, and TCPC. To examine the effects of the venous blood flow velocity, the velocity at rest and during exercise (0.5 and 1.0 W/kg) was measured. The total area of the thrombi increased over time. The APC square model showed the highest incidence for thrombus formation, followed by the APC round, whereas no thrombus was formed in the TCPC model. Slower blood flow at rest was associated with a higher incidence of thrombus formation. The TCPC was superior to the classical APC in terms of preventing thrombus formation, due to significant blood flow stagnation in the atrium of the APC. Thus, local hemodynamic behavior associated with the complex channel geometry plays a major role in thrombus formation in the Fontan route.

Patient-Specific Surgical Planning, Where Do We Stand? The Example of the Fontan Procedure

The Fontan surgery for single ventricle heart defects is a typical example of a clinical intervention in which patient-specific computational modeling can improve patient outcome: with the functional heterogeneity of the presenting patients, which precludes generic solutions, and the clear influence of the surgically-created Fontan connection on hemodynamics, it is acknowledged that individualized computational optimization of the post-operative hemodynamics can be of clinical value. A large body of literature has thus emerged seeking to provide clinically relevant answers and innovative solutions, with an increasing emphasis on patient-specific approaches. In this review we discuss the benefits and challenges of patient-specific simulations for the Fontan surgery, reviewing state of the art solutions and avenues for future development. We first discuss the clinical impact of patient-specific simulations, notably how they have contributed to our understanding of the link between Fontan hemodynamics and patient outcome. This is followed by a survey of methodologies for capturing patient-specific hemodynamics, with an emphasis on the challenges of defining patient-specific boundary conditions and their extension for prediction of post-operative outcome. We conclude with insights into potential future directions, noting that one of the most pressing issues might be the validation of the predictive capabilities of the developed framework.

An anatomically detailed arterial network model for one-dimensional computational hemodynamics

Simulation platforms are increasingly becoming complementary tools for cutting-edge cardiovascular research. The interplay among structural properties of the arterial wall, morphometry, anatomy, wave propagation phenomena, and ultimately, cardiovascular diseases continues to be poorly understood. Accurate models are powerful tools to shed light on these open problems. We developed an anatomically detailed computational model of the arterial vasculature to conduct 1-D blood flow simulations to serve as simulation infrastructure to aid cardiovascular research. An average arterial vasculature of a man was outlined in 3-D space to serve as geometrical substrate for the mathematical model. The architecture of this model comprises almost every arterial vessel acknowledged in the medical/anatomical literature, with a resolution down to the luminal area of perforator arteries. Over 2000 arterial vessels compose the model. Anatomical, physiological, and mechanical considerations were employed for the set up of model parameters and to determine criteria for blood flow distribution. Computational fluid dynamics was used to simulate blood flow and wave propagation phenomena in such arterial network. A sensitivity analysis was developed to unveil the contributions of model parameters to the conformation of the pressure waveforms. In addition, parameters were modified to target model to a patient-specific scenario. On the light of the knowledge domain, we conclude that the present model features excellent descriptive and predictive capabilities in both patient-generic and patient-specific cases, presenting a new step toward integrating an unprecedented anatomical description, morphometric, and simulations data to help in understanding complex arterial blood flow phenomena and related cardiovascular diseases.