Kartik S. Sundareswaran

The total cavopulmonary connection resistance: a significant impact on single ventricle hemodynamics at rest and exercise

Little is known about the impact of the total cavopulmonary connection (TCPC) on resting and exercise hemodynamics in a single ventricle (SV) circulation. The aim of this study was to elucidate this mechanism using a lumped parameter model of the SV circulation. Pulmonary vascular resistance (1.96+/-0.80 WU) and systemic vascular resistances (18.4+/-7.2 WU) were obtained from catheterization data on 40 patients with a TCPC. TCPC resistances (0.39+/-0.26 WU) were established using computational fluid dynamic simulations conducted on anatomically accurate three-dimensional models reconstructed from MRI (n=16). These parameters were used in a lumped parameter model of the SV circulation to investigate the impact of TCPC resistance on SV hemodynamics under resting and exercise conditions. A biventricular model was used for comparison. For a biventricular circulation, the cardiac output (CO) dependence on TCPC resistance was negligible (sensitivity=-0.064 l.min(-1).WU(-1)) but not for the SV circulation (sensitivity=-0.88 l.min(-1).WU(-1)). The capacity to increase CO with heart rate was also severely reduced for the SV. At a simulated heart rate of 150 beats/min, the SV patient with the highest resistance (1.08 WU) had a significantly lower increase in CO (20.5%) compared with the SV patient with the lowest resistance (50%) and normal circulation (119%). This was due to the increased afterload (+35%) and decreased preload (-12%) associated with the SV circulation. In conclusion, TCPC resistance has a significant impact on resting hemodynamics and the exercise capacity of patients with a SV physiology.

Patient-specific surgical planning and hemodynamic computational fluid dynamics optimization through free-form haptic anatomy editing tool (SURGEM)

The first version of an anatomy editing/surgical planning tool (SURGEM) targeting anatomical complexity and patient-specific computational fluid dynamics (CFD) analysis is presented. Novel three-dimensional (3D) shape editing concepts and human–shape interaction technologies have been integrated to facilitate interactive surgical morphology alterations, grid generation and CFD analysis. In order to implement “manual hemodynamic optimization” at the surgery planning phase for patients with congenital heart defects, these tools are applied to design and evaluate possible modifications of patient-specific anatomies. In this context, anatomies involve complex geometric topologies and tortuous 3D blood flow pathways with multiple inlets and outlets. These tools make it possible to freely deform the lumen surface and to bend and position baffles through real-time, direct manipulation of the 3D models with both hands, thus eliminating the tedious and time-consuming phase of entering the desired geometry using traditional computer-aided design (CAD) systems. The 3D models of the modified anatomies are seamlessly exported and meshed for patient-specific CFD analysis. Free-formed anatomical modifications are quantified using an in-house skeletization based cross-sectional geometry analysis tool. Hemodynamic performance of the systematically modified anatomies is compared with the original anatomy using CFD. CFD results showed the relative importance of the various surgically created features such as pouch size, vena cave to pulmonary artery (PA) flare and PA stenosis. An interactive surgical-patch size estimator is also introduced. The combined design/analysis cycle time is used for comparing and optimizing surgical plans and improvements are tabulated. The reduced cost of patient-specific shape design and analysis process, made it possible to envision large clinical studies to assess the validity of predictive patient-specific CFD simulations. In this paper, model anatomical design studies are performed on a total of eight different complex patient specific anatomies. Using SURGEM, more than 30 new anatomical designs (or candidate configurations) are created, and the corresponding user times presented. CFD performances for eight of these candidate configurations are also presented.

A New Method for Registration-Based Medical Image Interpolation

A new technique is presented for interpolating between grey-scale images in a medical data set. Registration between neighboring slices is achieved with a modified control grid interpolation algorithm that selectively accepts displacement field updates in a manner optimized for performance. A cubic interpolator is then applied to pixel intensities correlated by the displacement fields. Special considerations are made for efficiency, interpolation quality, and compression in the implementation of the algorithm. Experimental results show that the new method achieves good quality, while offering dramatic improvement in efficiency relative to the best competing method.

Fontan hemodynamics: Importance of pulmonary artery diameter

OBJECTIVE We quantify the geometric and hemodynamic characteristics of extracardiac and lateral tunnel Fontan surgical options and correlate certain anatomic characteristics with their hemodynamic efficiency and patient cardiac index. METHODS AND RESULTS The study was conducted retrospectively on 22 patients undergoing Fontan operations (11 extracardiac and 11 lateral tunnel operations). Total cavopulmonary connection geometric parameters such as vessel areas, curvature, and offsets were quantified using a skeletonization method. Energy loss at the total cavopulmonary connection junction was available from previous in vitro experiments and computational fluid dynamic simulations for 5 and 9 patients, respectively. Cardiac index data were available for all patients. There was no significant difference in the mean and minimum cross-sectional vessel areas of the pulmonary artery between the extracardiac and lateral tunnel groups. The indexed energy dissipation within the total cavopulmonary connection was strongly correlated to minimum cross-sectional area of the pulmonary arteries (R(2) value of 0.90 and P < .0002), whereas all other geometric features, including shape characteristics, had no significant correlation. Finally, cardiac index significantly correlated with the minimum pulmonary artery area (P = .006), suggesting that total cavopulmonary connection energy losses significantly affect resting cardiac output. CONCLUSIONS The minimum outlet size of the total cavopulmonary connection (ie, minimum cross section of pulmonary artery) governs the energy loss characteristics of the total cavopulmonary connection more strongly than variations in the shapes corresponding to extracardiac and lateral tunnel configurations. Differences in pulmonary artery sizes must be accounted for when comparing energy losses between extracardiac and lateral tunnel geometries.

Correction of Pulmonary Arteriovenous Malformation Using Image-Based Surgical Planning

The objectives of this study were to develop an image-based surgical planning framework that 1) allows for in-depth analysis of pre-operative hemodynamics by the use of cardiac magnetic resonance and 2) enables surgeons to determine the optimum surgical scenarios before the operation. This framework is tailored for applications in which post-operative hemodynamics are important. In particular, it is exemplified here for a Fontan patient with severe left pulmonary arteriovenous malformations due to abnormal hepatic flow distribution to the lungs. Patients first undergo cardiac magnetic resonance for 3-dimensional anatomy and flow reconstruction. After analysis of the pre-operative flow fields, the 3-dimensional anatomy is imported into an interactive surgical planning interface for the surgeon to virtually perform multiple surgical scenarios. Associated hemodynamics are predicted by the use of a fully validated computational fluid dynamic solver. Finally, efficiency metrics (e.g., pressure decrease and hepatic flow distribution) are weighted against surgical feasibility to determine the optimal surgical option.

Neonatal Aortic Arch Hemodynamics and Perfusion During Cardiopulmonary Bypass

The objective of this study is to quantify the detailed three-dimensional (3D) pulsatile hemodynamics, mechanical loading, and perfusion characteristics of a patient-specific neonatal aortic arch during cardiopulmonary bypass (CPB). The 3D cardiac magnetic resonance imaging (MRI) reconstruction of a pediatric patient with a normal aortic arch is modified based on clinical literature to represent the neonatal morphology and flow conditions. The anatomical dimensions are verified from several literature sources. The CPB is created virtually in the computer by clamping the ascending aorta and inserting the computer-aided design model of the 10 Fr tapered generic cannula. Pulsatile (130 bpm) 3D blood flow velocities and pressures are computed using the commercial computational fluid dynamics (CFD) software. Second order accurate CFD settings are validated against particle image velocimetry experiments in an earlier study with a complex cardiovascular unsteady benchmark. CFD results in this manuscript are further compared with the in vivo physiological CPB pressure waveforms and demonstrated excellent agreement. Cannula inlet flow waveforms are measured from in vivo PC-MRI and 3 kg piglet neonatal animal model physiological experiments, distributed equally between the head-neck vessels and the descending aorta. Neonatal 3D aortic hemodynamics is also compared with that of the pediatric and fetal aortic stages. Detailed 3D flow fields, blood damage, wall shear stress (WSS), pressure drop, perfusion, and hemodynamic parameters describing the pulsatile energetics are calculated for both the physiological neonatal aorta and for the CPB aorta assembly. The primary flow structure is the high-speed canulla jet flow (approximately 3.0 m/s at peak flow), which eventually stagnates at the anterior aortic arch wall and low velocity flow in the cross-clamp pouch. These structures contributed to the reduced flow pulsatility (85%), increased WSS (50%), power loss (28%), and blood damage (288%), compared with normal neonatal aortic physiology. These drastic hemodynamic differences and associated intense biophysical loading of the pathological CPB configuration necessitate urgent bioengineering improvements--in hardware design, perfusion flow waveform, and configuration. This study serves to document the baseline condition, while the methodology presented can be utilized in preliminary CPB cannula design and in optimization studies reducing animal experiments. Coupled to a lumped-parameter model the 3D hemodynamic characteristics will aid the surgical decision making process of the perfusion strategies in complex congenital heart surgeries.

Blood flow distribution in a large series of patients having the Fontan operation: A cardiac magnetic resonance velocity mapping study

OBJECTIVES Our goal was to determine flow distribution in the cavopulmonary connections of patients with and without bilateral superior venae cavae who had the Fontan procedure. No large series exists that establishes the flow distributions in Fontan patients, which would be an important resource for everyday clinical use and may affect future surgical reconstruction. METHODS We studied 105 Fontan patients (aged 2-24 years) with through-plane phase contrast velocity mapping to determine flow rates in the inferior and superior venae cavae and left and right pulmonary arteries. Superior caval anastomosis type included 40 bidirectional Glenn shunts (of which 15 were bilateral) and 53 hemi-Fontan anastomoses; Fontan type included 69 intra-atrial baffles, 28 extracardiac conduits, and 4 atriopulmonary connections. RESULTS Total caval flow was 2.9 +/- 1.0 L x min(-1) x m(-2), with an inferior vena cava contribution of 59% +/- 15%. Total pulmonary flow was 2.5 +/- 0.8 L x min(-1) x m(-2), statistically less than caval flow and not explained by fenestration presence. The right pulmonary artery contribution (55% +/- 13%) was statistically greater than the left. In patients with bilateral superior cavae, the right cava accounted for 52% +/- 14% of the flow, with no difference in pulmonary flow splits (50% +/- 16% to the right). Age and body surface area correlated with percent inferior caval contribution (r = 0.60 and 0.74, respectively). Superior vena cava anastomosis and Fontan type did not significantly affect pulmonary flow splits. CONCLUSIONS Total Fontan cardiac index was 2.9 L x min(-1) x m(-2), with normal pulmonary flow splits (55% to the right lung). Inferior vena caval contribution to total flow increases with body surface area and age, consistent with data from healthy children.

Pulmonary hepatic flow distribution in total cavopulmonary connections: extracardiac versus intracardiac.

OBJECTIVE Pulmonary arteriovenous malformations can occur after the Fontan procedure and are believed to be associated with disproportionate pulmonary distribution of hepatic venous effluent. We studied the effect of total cavopulmonary connection geometry and the effect of increased cardiac output on distribution of inferior vena caval return to the lungs. METHODS Ten patients undergoing the Fontan procedure, 5 with extracardiac and 5 with intracardiac configurations of the total cavopulmonary connection, previously analyzed for power loss were processed for calculating the distribution of inferior vena caval return to the lungs (second-order accuracy). One idealized total cavopulmonary connection was similarly analyzed under parametric variation of inferior vena caval offset and cardiac output flow split. RESULTS Streaming of the inferior vena caval return in the idealized total cavopulmonary connection model was dependent on both inferior vena caval offset magnitude and cardiac output flow-split ratio. For patient-specific total cavopulmonary connections, preferential streaming of the inferior vena caval return was directly proportional to the cardiac output flow-split ratio in the intracardiac total cavopulmonary connections (P < .0001). Preferential streaming in extracardiac total cavopulmonary connections correlated to the inferior vena caval offset (P < .05) and did not correlate to cardiac output flow split. Enhanced mixing in intracardiac total cavopulmonary connections is speculated to explain the contrasting results. Exercising tends to reduce streaming toward the left pulmonary artery in intracardiac total cavopulmonary connections, whereas for extracardiac total cavopulmonary connections, exercising tends to equalize the streaming. CONCLUSIONS Extracardiac and intracardiac total cavopulmonary connections have inherently different streaming characteristics because of contrasting mixing characteristics caused by their geometric differences. Pulmonary artery diameters and inferior vena caval offsets might together determine hepatic flow streaming.

Individualized computer-based surgical planning to address pulmonary arteriovenous malformations in patients with a single ventricle with an interrupted inferior vena cava and azygous continuation

OBJECTIVE Pulmonary arteriovenous malformations caused by abnormal hepatic flow distribution can develop in patients with a single ventricle with an interrupted inferior vena cava. However, preoperatively determining the hepatic baffle design that optimizes hepatic flow distribution is far from trivial. The current study combines virtual surgery and numeric simulations to identify potential surgical strategies for patients with an interrupted inferior vena cava. METHODS Five patients with an interrupted inferior vena cava and severe pulmonary arteriovenous malformations were enrolled. Their in vivo anatomies were reconstructed from magnetic resonance imaging (n = 4) and computed tomography (n = 1), and alternate virtual surgery options (intracardiac/extracardiac, Y-grafts, hepato-to-azygous shunts, and azygous-to-hepatic shunts) were generated for each. Hepatic flow distribution was assessed for all options using a fully validated computational flow solver. RESULTS For patients with a single superior vena cava (n = 3), intracardiac/extracardiac connections proved dangerous, because even a small left or right offset led to a highly preferential hepatic flow distribution to the associated lung. The best results were obtained with either a Y-graft spanning the Kawashima to split the flow or hepato-to-azygous shunts to promote mixing. For patients with bilateral superior vena cavae (n = 2), results depended on the balance between the left and right superior inflows. When those were equal, connecting the hepatic baffle between the superior vena cavae performed well, but other options should be pursued otherwise. CONCLUSIONS This study demonstrates how virtual surgery environments can benefit the clinical community, especially for patients with a single ventricle with an interrupted inferior vena cava. Furthermore, the sensitivity of the optimal baffle design to the superior inflows underscores the need to characterize both preoperative anatomy and flows to identify the best option.

Fontan hemodynamics from 100 patient-specific cardiac magnetic resonance studies: a computational fluid dynamics analysis.

OBJECTIVES This study sought to quantify average hemodynamic metrics of the Fontan connection as reference for future investigations, compare connection types (intra-atrial vs extracardiac), and identify functional correlates using computational fluid dynamics in a large patient-specific cohort. Fontan hemodynamics, particularly power losses, are hypothesized to vary considerably among patients with a single ventricle and adversely affect systemic hemodynamics and ventricular function if suboptimal. METHODS Fontan connection models were created from cardiac magnetic resonance scans for 100 patients. Phase velocity cardiac magnetic resonance in the aorta, vena cavae, and pulmonary arteries was used to prescribe patient-specific time-averaged flow boundary conditions for computational fluid dynamics with a customized, validated solver. Comparison with 4-dimensional cardiac magnetic resonance velocity data from selected patients was used to provide additional verification of simulations. Indexed Fontan power loss, connection resistance, and hepatic flow distribution were quantified and correlated with systemic patient characteristics. RESULTS Indexed power loss varied by 2 orders of magnitude, whereas, on average, Fontan resistance was 15% to 20% of published values of pulmonary vascular resistance in single ventricles. A significant inverse relationship was observed between indexed power loss and both systemic venous flow and cardiac index. Comparison by connection type showed no differences between intra-atrial and extracardiac connections. Instead, the least efficient connections revealed adverse consequences from localized Fontan pathway stenosis. CONCLUSIONS Fontan power loss varies from patient to patient, and elevated levels are correlated with lower systemic flow and cardiac index. Fontan connection type does not influence hemodynamic efficiency, but an undersized or stenosed Fontan pathway or pulmonary arteries can be highly dissipative.