Simulations reveal adverse hemodynamics in patients with multiple systemic to pulmonary shunts.

For newborns diagnosed with pulmonary atresia or severe pulmonary stenosis leading to insufficient pulmonary blood flow, cyanosis can be mitigated with placement of a modified Blalock-Taussig shunt (MBTS) between the innominate and pulmonary arteries. In some clinical scenarios, patients receive two systemic-to-pulmonary connections, either by leaving the patent ductus arteriosus (PDA) open or by adding an additional central shunt (CS) in conjunction with the MBTS. This practice has been motivated by the thinking that an additional source of pulmonary blood flow could beneficially increase pulmonary flow and provide the security of an alternate pathway in case of thrombosis. However, there have been clinical reports of premature shunt occlusion when more than one shunt is employed, leading to speculation that multiple shunts may in fact lead to unfavorable hemodynamics and increased mortality. In this study, we hypothesize that multiple shunts may lead to undesirable flow competition, resulting in increased residence time (RT) and elevated risk of thrombosis, as well as pulmonary overcirculation. Computational fluid dynamics-based multiscale simulations were performed to compare a range of shunt configurations and systematically quantify flow competition, pulmonary circulation, and other clinically relevant parameters. In total, 23 cases were evaluated by systematically changing the PDA/CS diameter, pulmonary vascular resistance (PVR), and MBTS position and compared by quantifying oxygen delivery (OD) to the systemic and coronary beds, wall shear stress (WSS), oscillatory shear index (OSI), WSS gradient (WSSG), and RT in the pulmonary artery (PA), and MBTS. Results showed that smaller PDA/CS diameters can lead to flow conditions consistent with increased thrombus formation due to flow competition in the PA, and larger PDA/CS diameters can lead to insufficient OD due to pulmonary hyperfusion. In the worst case scenario, it was found that multiple shunts can lead to a 160% increase in RT and a 10% decrease in OD. Based on the simulation results presented in this study, clinical outcomes for patients receiving multiple shunts should be critically investigated, as this practice appears to provide no benefit in terms of OD and may actually increase thrombotic risk.

Numerical blood flow simulation in surgical corrections: what do we need for an accurate analysis?

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.

Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge.

Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

The assisted bidirectional Glenn: a novel surgical approach for first-stage single-ventricle heart palliation.

BACKGROUND: Outcomes after a modified Blalock-Taussig shunt (mBTS) in neonates with single-ventricle physiology remain unsatisfactory. However, initial palliation with a superior cavopulmonary connection, such as a bidirectional Glenn (BDG), is discouraged, owing to potential for inadequate pulmonary blood flow (PBF). We tested the feasibility of a novel surgical approach, adopting the engineering concept of an ejector pump, whereby the flow in the BDG is "assisted" by injection of a high-energy flow stream from the systemic circulation. METHODS: Realistic 3-dimensional models of the neonatal mBTS and BDG circulations were created. The "assisted" bidirectional Glenn (ABG) consisted of a shunt between the right innominate artery and the superior vena cava (SVC), with a 1.5-mm clip near the SVC anastomosis to create a Venturi effect. The 3 models were coupled to a validated hydraulic circulation model, and 2 pulmonary vascular resistance (PVR) values (7 and 2.3 Wood units) were simulated. RESULTS: The ABG provided the highest systemic oxygen saturation and oxygen delivery at both PVR levels. In addition to achieving higher PBF than the BDG, the ABG produced a lower single-ventricular workload than mBTS. SVC pressure was highest in the ABG model (ABG: 15; Glenn: 11; mBTS: 3 mm Hg; PVR = 7 Wood units), but at low PVR, the SVC pressure was significantly lower (ABG: 8; Glenn: 6; mBTS: <3 mm Hg). CONCLUSIONS: Adopting the principle of an ejector pump, with additional flow directed into the SVC in a BDG, the ABG appears to increase PBF with a modest increase in SVC and pulmonary arterial pressure. Although multiscale modeling results demonstrate the conceptual feasibility of the ABG circulation, further technical refinement and investigations are necessary, especially in an appropriate animal model.

In Vitro Assessment of the Assisted Bidirectional Glenn Procedure for Stage One Single Ventricle Repair.

This in vitro study compares the hemodynamic performance of the Norwood and the Glenn circulations to assess the performance of a novel assisted bidirectional Glenn (ABG) procedure for stage one single ventricle surgery. In the ABG, the flow in a bidirectional Glenn procedure is assisted by injection of a high-energy flow stream from the systemic circulation using an aorta-caval shunt with nozzle. The aim is to explore experimentally the potential of the ABG as a surgical alternative to current surgical practice. The experiments are directly compared against previously published numerical simulations. A multiscale mock circulatory system was used to measure the hemodynamic performance of the three circulations. For each circulation, the system was tested using both low and high values of pulmonary vascular resistance. Resulting parameters measured were: pressure and flow rate at left/right pulmonary artery and superior vena cava (SVC). Systemic oxygen delivery (OD) was calculated. A parametric study of the ratio of ABG nozzle to shunt diameter was done. We report time-based comparisons with numerical simulations for the three surgical variants tested. The ABG circulation demonstrated an increase of 30-38% in pulmonary flow with a 2-3.7 mmHg increase in SVC pressure compared to the Glenn and a 4-14% higher systemic OD than either the Norwood or the Glenn. The nozzle/shunt diameter ratio affected the local hemodynamics. These experimental results agreed with those of the numerical model: mean flow values were not significantly different (p > 0.05) while mean pressures were comparable within 1.2 mmHg. The results verify the approaches providing two tools to study this complicated circulation. Using a realistic experimental model we demonstrate the performance of a novel surgical procedure with potential to improve patient hemodynamics in early palliation of the univentricular circulation.

Thrombotic risk stratification using computational modeling in patients with coronary artery aneurysms following Kawasaki disease.

Kawasaki disease (KD) is the leading cause of acquired heart disease in children and can result in life-threatening coronary artery aneurysms in up to 25 % of patients. These aneurysms put patients at risk of thrombus formation, myocardial infarction, and sudden death. Clinicians must therefore decide which patients should be treated with anticoagulant medication, and/or surgical or percutaneous intervention. Current recommendations regarding initiation of anticoagulant therapy are based on anatomy alone with historical data suggesting that patients with aneurysms [Formula: see text]8 mm are at greatest risk of thrombosis. Given the multitude of variables that influence thrombus formation, we postulated that hemodynamic data derived from patient-specific simulations would more accurately predict risk of thrombosis than maximum diameter alone. Patient-specific blood flow simulations were performed on five KD patients with aneurysms and one KD patient with normal coronary arteries. Key hemodynamic and geometric parameters, including wall shear stress, particle residence time, and shape indices, were extracted from the models and simulations and compared with clinical outcomes. Preliminary fluid structure interaction simulations with radial expansion were performed, revealing modest differences in wall shear stress compared to the rigid wall case. Simulations provide compelling evidence that hemodynamic parameters may be a more accurate predictor of thrombotic risk than aneurysm diameter alone and motivate the need for follow-up studies with a larger cohort. These results suggest that a clinical index incorporating hemodynamic information be used in the future to select patients for anticoagulant therapy.

A non-discrete method for computation of residence time in fluid mechanics simulations.

Cardiovascular simulations provide a promising means to predict risk of thrombosis in grafts, devices, and surgical anatomies in adult and pediatric patients. Although the pathways for platelet activation and clot formation are not yet fully understood, recent findings suggest that thrombosis risk is increased in regions of flow recirculation and high residence time (RT). Current approaches for calculating RT are typically based on releasing a finite number of Lagrangian particles into the flow field and calculating RT by tracking their positions. However, special care must be taken to achieve temporal and spatial convergence, often requiring repeated simulations. In this work, we introduce a non-discrete method in which RT is calculated in an Eulerian framework using the advection-diffusion equation. We first present the formulation for calculating residence time in a given region of interest using two alternate definitions. The physical significance and sensitivity of the two measures of RT are discussed and their mathematical relation is established. An extension to a point-wise value is also presented. The methods presented here are then applied in a 2D cavity and two representative clinical scenarios, involving shunt placement for single ventricle heart defects and Kawasaki disease. In the second case study, we explored the relationship between RT and wall shear stress, a parameter of particular importance in cardiovascular disease.

Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery.

We present a computational framework for multiscale modeling and simulation of blood flow in coronary artery bypass graft (CABG) patients. Using this framework, only CT and non-invasive clinical measurements are required without the need to assume pressure and/or flow waveforms in the coronaries and we can capture global circulatory dynamics. We demonstrate this methodology in a case study of a patient with multiple CABGs. A patient-specific model of the blood vessels is constructed from CT image data to include the aorta, aortic branch vessels (brachiocephalic artery and carotids), the coronary arteries and multiple bypass grafts. The rest of the circulatory system is modeled using a lumped parameter network (LPN) 0 dimensional (0D) system comprised of resistances, capacitors (compliance), inductors (inertance), elastance and diodes (valves) that are tuned to match patient-specific clinical data. A finite element solver is used to compute blood flow and pressure in the 3D (3 dimensional) model, and this solver is implicitly coupled to the 0D LPN code at all inlets and outlets. By systematically parameterizing the graft geometry, we evaluate the influence of graft shape on the local hemodynamics, and global circulatory dynamics. Virtual manipulation of graft geometry is automated using Bezier splines and control points along the pathlines. Using this framework, we quantify wall shear stress, wall shear stress gradients and oscillatory shear index for different surgical geometries. We also compare pressures, flow rates and ventricular pressure-volume loops pre- and post-bypass graft surgery. We observe that PV loops do not change significantly after CABG but that both coronary perfusion and local hemodynamic parameters near the anastomosis region change substantially. Implications for future patient-specific optimization of CABG are discussed.