Search

Introduction: The increased survival rate of patients with congenital heart disease (CHD) has made it likely that 70%-95% of infants with CHDs surviving into adulthood often require careful follow-up and (repeat) interventions. Patients with CHDs often have abnormal blood flow patterns, due to both primary cardiac defect and the consequent surgical or endovascular repair., Area Covered: Computational fluid dynamics (CFD) alone or coupled with advanced imaging tools can assess blood flow patterns of CHDs to both understand their pathophysiology and anticipate the results of surgical or interventional repair., Expert Opinion: CFD is a mathematical technique that quantifies and describes the characteristics of fluid flow using the laws of physics. Through dedicated software based on virtual reconstruction and simulation and patients' real data coming from computed tomography, magnetic resonance imaging, and 3/4 D-ultrasound, reconstruction of models of circulation of most CHD can be accomplished. CFD can provide insights about the pathophysiology of coronary artery anomalies, interatrial shunts, coarctation of the aorta and aortic bicuspid valve, tetralogy of Fallot and univentricular heart, with the capability in some cases of simulating different types of surgical or interventional repair and tailoring the treatment on the basis of these findings.

Superior cavopulmonary circulation (SCPC) can be achieved by either the Hemi-Fontan (hF) or Bidirectional Glenn (bG) connection. Debate remains as to which results in best hemodynamic results. Adopting patient-specific multiscale computational modeling, we examined both the local dynamics and global physiology to determine if surgical choice can lead to different hemodynamic outcomes. Six patients (age: 3-6 months) underwent cardiac magnetic resonance imaging and catheterization prior to SCPC surgery. For each patient: (1) a finite 3-dimensional (3D) volume model of the preoperative anatomy was constructed to include detailed definition of the distal branch pulmonary arteries, (2) virtual hF and bG operations were performed to create 2 SCPC 3D models, and (3) a specific lumped network representing each patient's entire cardiovascular circulation was developed from clinical data. Using a previously validated multiscale algorithm that couples the 3D models with lumped network, both local flow dynamics, that is, power loss, and global systemic physiology can be quantified. In 2 patients whose preoperative imaging demonstrated significant left pulmonary artery (LPA) stenosis, we performed virtual pulmonary arterioplasty to assess its effect. In one patient, the hF model showed higher power loss (107%) than the bG, while in 3, the power losses were higher in the bG models (18-35%). In the remaining 2 patients, the power loss differences were minor. Despite these variations, for all patients, there were no significant differences between the hF and bG models in hemodynamic or physiological outcomes, including cardiac output, superior vena cava pressure, right-left pulmonary flow distribution, and systemic oxygen delivery. In the 2 patients with LPA stenosis, arterioplasty led to better LPA flow (5-8%) while halving the power loss, but without important improvements in SVC pressure or cardiac output. Despite power loss differences, both hF and bG result in similar SCPC hemodynamics and physiology outcome. This suggests that for SCPC, the pre-existing patient-specific physiology and condition, such as pulmonary vascular resistance, are more deterministic in the hemodynamic performance than the type of surgical palliation. Multiscale modeling can be a decision-assist tool to assess whether an extensive LPA reconstruction is needed at the time of SCPC for LPA stenosis., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Single ventricle heart defects involve pathologies in which the heart has only one functional pumping chamber. In these conditions, treatment consists of three staged procedures. At stage 1 pulmonary flow is provided through an artificial shunt from the systemic circulation. Representative hemodynamics models able to explore different virtual surgical options can be built based on pre-operative imaging and patient data. In this context, the specification of boundary conditions is necessary to compute pressure and flow in the entire domain. However, these boundary conditions are rarely the measured variables. Moreover, to take into account the rest of the circulation outside of the three-dimensional modeled domain, a number of reduced order models exist. A simplified method is presented to iteratively, but automatically, tune reduced model parameters from hemodynamic data clinically measured before stage 2 surgery. Patient-specific local hemodynamics around the distal systemic-to-pulmonary shunt anastomosis and the connected pulmonary arteries are also analyzed. Multi-scale models of pre-stage 2 single ventricle patients are developed, including a 3D model of shunt-pulmonary connection and a number of pulmonary arteries. For each pulmonary outlet a total downstream resistance is identified, consistent with measured flow split and pressures. Target pressures such as minimum, maximum or average over one or both lungs are considered, depending on the clinical measurement. When possible, both steady and pulsatile identifications are performed. The methodology is demonstrated with six patient-specific models: the clinical target data are well-matched, except for one case where clinical data were subsequently found inconsistent. Inhomogeneous pressure, swirling blood flow patterns and very high wall shear stress 3D maps highlight similarities and differences among patients. Steady and pulsatile tuning results are similar. This work demonstrates (1) how to use routine clinical data to define boundary conditions for patient-specific 3D models in pre-stage 2 single ventricle circulations and (2) how simulations can help to check the coherence of clinical data, or provide insights to clinicians that are otherwise difficult to measure, such as in the presence of kinks. Finally, the choice of steady vs. pulsatile tuning, limitations and possible extensions of this work are discussed.

Newborns with single ventricle physiology are usually palliated with a multi-staged procedure. When cardiovascular complications e.g., collateral vessel formation occur during the inter-stage periods, further treatments are required. An 8-month-old patient, who underwent second stage (i.e., bi-directional Glenn, BDG) surgery at 4 months, was diagnosed with a major veno-venous collateral vessel (VVC) which was endovascularly occluded to improve blood oxygen saturations. Few clinical data were collected at 8 months, whereas at 4 months a more detailed data set was available. The aim of this study is threefold: (i) to show how to build a patient-specific model describing the hemodynamics in the presence of VVC, using patient-specific clinical data collected at different times; (ii) to use this model to perform virtual VVC occlusion for quantitative hemodynamics prediction; and (iii) to compare predicted hemodynamics with post-operative measurements. The three-dimensional BDG geometry, resulting from the virtual surgery on the first stage model, was coupled with a lumped parameter model (LPM) of the 8-month patient's circulation. The latter was developed by scaling the 4-month LPM to account for changes in vascular impedances due to growth and adaptation. After virtual VVC closure, the model confirmed the 2 mmHg BDG pressure increase, as clinically observed, suggesting the importance of modeling vascular adaptation following the BDG procedure.

Objective: Currently, no quantitative guidelines have been established for treatment of left pulmonary artery (LPA) stenosis. This study aims to quantify the effects of LPA stenosis on postoperative hemodynamics for single-ventricle patients undergoing stage II superior cavopulmonary connection (SCPC) surgery, using a multiscale computational approach., Methods: Image data from 6 patients were segmented to produce 3-dimensional models of the pulmonary arteries before stage II surgery. Pressure and flow measurements were used to tune a 0-dimensional model of the entire circulation. Postoperative geometries were generated through stage II virtual surgery; varying degrees of LPA stenosis were applied using mesh morphing and hemodynamics assessed through coupled 0-3-dimensional simulations. To relate metrics of stenosis to clinical classifications, pediatric cardiologists and surgeons ranked the degrees of stenosis in the models. The effects of LPA stenosis were assessed based on left-to-right pulmonary artery flow split ratios, mean pressure drop across the stenosis, cardiac pressure-volume loops, and other clinically relevant parameters., Results: Stenosis of >65% of the vessel diameter was required to produce a right pulmonary artery:LPA flow split <30%, and/or a mean pressure drop of >3.0 mm Hg, defined as clinically significant changes., Conclusions: The effects of <65% stenosis on SCPC hemodynamics and physiology were minor and may not justify the increased complexity of adding LPA arterioplasty to the SCPC operation. However, in the longer term, pulmonary augmentation may affect outcomes of the Fontan completion surgery, as pulmonary artery distortion is a risk factor that may influence stage III physiology., (Copyright © 2015 The American Association for Thoracic Surgery. Published by Elsevier Inc. All rights reserved.)

Background: Computational fluid dynamics has been increasingly used in congenital heart surgery to simulate pathophysiological blood flow, investigate surgical options, or design medical devices. Several commercial and research computational or numerical codes have been developed. They present different approaches to numerically solve the blood flow equations, raising the question whether these numerical codes are equally reliable to achieve accurate simulation results. Accordingly, we sought to examine the influence of numerical code selection in several complex congenital cardiac operations., Material and Methods: The main steps of blood flow simulations are detailed (geometrical mesh, boundary conditions, and solver numerical methods) for congenital cardiac operations of increasing complexity. The first case tests different numerical solutions against an analytical, or exact, solution. In the second case, the three-dimensional domain is a patient-specific superior cavopulmonary anastomosis. As an analytical solution does not exist in such a complex geometry, different numerical solutions are compared. Finally, a realistic case of a systemic-to-pulmonary shunt is presented with both geometrically and physiologically challenging conditions. For all, solutions from a commercially available code and an open-source research code are compared., Results: In the first case, as the mesh or solver numerical method is refined, the simulation results for both codes converged to the analytical solution. In the second example, velocity differences between the two codes are greater when the resolution of the mesh were lower and less refined. The third case with realistic anatomy reveals that the pulsatile complex flow is very similar for both codes., Conclusions: The precise setup of the numerical cases has more influence on the results than the choice of numerical codes. The need for detailed construction of the numerical model that requires high computational cost depends on the precision needed to answer the biomedical question at hand and should be assessed for each problem on a combination of clinically relevant patient-specific geometry and physiological conditions., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Single ventricle hearts are congenital cardiovascular defects in which the heart has only one functional pumping chamber. The treatment for these conditions typically requires a three-staged operative process where Stage 1 is typically achieved by a shunt between the systemic and pulmonary arteries, and Stage 2 by connecting the superior venous return to the pulmonary circulation. Surgically, the Stage 2 circulation can be achieved through a procedure called the Hemi-Fontan, which reconstructs the right atrium and pulmonary artery to allow for an enlarged confluence with the superior vena cava. Based on pre-operative data obtained from two patients prior to Stage 2 surgery, we developed two patient-specific multi-scale computational models, each including the 3D geometrical model of the surgical junction constructed from magnetic resonance imaging, and a closed-loop systemic lumped-parameter network derived from clinical measurements. "Virtual" Hemi-Fontan surgery was performed on the 3D model with guidance from clinical surgeons, and a corresponding multi-scale simulation predicts the patient's post-operative hemodynamic and physiologic conditions. For each patient, a post-operative active scenario with an increase in the heart rate (HR) and a decrease in the pulmonary and systemic vascular resistance (PVR and SVR) was also performed. Results between the baseline and this "active" state were compared to evaluate the hemodynamic and physiologic implications of changing conditions. Simulation results revealed a characteristic swirling vortex in the Hemi-Fontan in both patients, with flow hugging the wall along the SVC to Hemi-Fontan confluence. One patient model had higher levels of swirling, recirculation, and flow stagnation. However, in both models, the power loss within the surgical junction was less than 13% of the total power loss in the pulmonary circulation, and less than 2% of the total ventricular power. This implies little impact of the surgical junction geometry on the SVC pressure, cardiac output, and other systemic parameters. In contrast, varying HR, PVR, and SVR led to significant changes in theses clinically relevant global parameters. Adopting a work-flow of customized virtual planning of the Hemi-Fontan procedure with patient-specific data, this study demonstrates the ability of multi-scale modeling to reproduce patient specific flow conditions under differing physiological states. Results demonstrate that the same operation performed in two different patients can lead to different hemodynamic characteristics, and that modeling can be used to uncover physiologic changes associated with different clinical conditions., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

The objective of this work is to perform a virtual planning of surgical repairs in patients with congenital heart diseases--to test the predictive capability of a closed-loop multi-scale model. As a first step, we reproduced the pre-operative state of a specific patient with a univentricular circulation and a bidirectional cavopulmonary anastomosis (BCPA), starting from the patient's clinical data. Namely, by adopting a closed-loop multi-scale approach, the boundary conditions at the inlet and outlet sections of the three-dimensional model were automatically calculated by a lumped parameter network. Successively, we simulated three alternative surgical designs of the total cavopulmonary connection (TCPC). In particular, a T-junction of the venae cavae to the pulmonary arteries (T-TCPC), a design with an offset between the venae cavae (O-TCPC) and a Y-graft design (Y-TCPC) were compared. A multi-scale closed-loop model consisting of a lumped parameter network representing the whole circulation and a patient-specific three-dimensional finite volume model of the BCPA with detailed pulmonary anatomy was built. The three TCPC alternatives were investigated in terms of energetics and haemodynamics. Effects of exercise were also investigated. Results showed that the pre-operative caval flows should not be used as boundary conditions in post-operative simulations owing to changes in the flow waveforms post-operatively. The multi-scale approach is a possible solution to overcome this incongruence. Power losses of the Y-TCPC were lower than all other TCPC models both at rest and under exercise conditions and it distributed the inferior vena cava flow evenly to both lungs. Further work is needed to correlate results from these simulations with clinical outcomes.

Total cavopulmonary connection is a surgical procedure adopted to treat complex congenital malformations of the right heart. It consists basically in a connection of both venae cavae directly to the right pulmonary artery. In this paper a three-dimensional model of this connection is presented, which is based on in vivo measurements performed by means of magnetic resonance. The model was developed by means of computational fluid dynamics techniques, namely the finite element method. The aim of this study was to verify the capability of such a model to predict the distribution of the blood flow into the pulmonary arteries, by comparison with in vivo velocity measurements. Different simulations were performed on a single clinical case to test the sensitivity of the model to different boundary conditions, in terms of inlet velocity profiles as well as outlet pressure levels. Results showed that the flow distribution between the lungs is slightly affected by the shape of inlet velocity profiles, whereas it is influenced by different pressure levels to a greater extent.

The bidirectional cavopulmonary anastomosis is used as a staged procedure or a definitive palliation of univentricular hearts. It is often performed in the presence of an additional blood flow arising from the native pulmonary outflow tract. In this paper, the effects of the severity of the pulmonary outflow obstruction and the pulmonary arteriolar resistance are analysed with regard to the haemodynamics in the superior vena cava and the blood distribution into the lungs. A computer model has been developed, which can represent both the preoperative and the postoperative (systemic and pulmonary) circulations in a patient with a double-outlet univentricular heart. It is particularly detailed in the region of the large vessels and includes components that account for local three-dimensional effects due to the actual shape of the anastomosis. Results have indicated that the mean pressure in the superior vena cava increases from 8.2 to 19.2 mmHg with pulmonary arteriolar resistance ranging from 0.8 to 7.9 Woods units and pulmonary outflow obstruction ranging from 50 to 100%. The percentage flow distribution to the right lung has turned out to be heavily affected by the flow competition and has ranged from 43 to 50% of the total flow to the lungs in the systolic phase, and from 51 to 62% in the diastolic phase. The model allows routinely used clinical indices to be computed, as well as the evaluation of new indices, which is potentially helpful in the clinical assessment of postoperative haemodynamics (e.g. the right-to-left lung flow ratio and the superior vena cava-to-pulmonary flow ratio).