Patient-Specific Multi-Scale Model Analysis of Hemodynamics Following the Hybrid Norwood Procedure for Hypoplastic Left Heart Syndrome: Effects of Reverse Blalock-Taussig Shunt Diameter

Impact of TAVR on coronary artery hemodynamics using clinical measurements and image-based patient-specific in silico modeling

In recent years, transcatheter aortic valve replacement (TAVR) has become the leading method for treating aortic stenosis. While the procedure has improved dramatically in the past decade, there are still uncertainties about the impact of TAVR on coronary blood flow. Recent research has indicated that negative coronary events after TAVR may be partially driven by impaired coronary blood flow dynamics. Furthermore, the current technologies to rapidly obtain non-invasive coronary blood flow data are relatively limited. Herein, we present a lumped parameter computational model to simulate coronary blood flow in the main arteries as well as a series of cardiovascular hemodynamic metrics. The model was designed to only use a few inputs parameters from echocardiography, computed tomography and a sphygmomanometer. The novel computational model was then validated and applied to 19 patients undergoing TAVR to examine the impact of the procedure on coronary blood flow in the left anterior descending (LAD) artery, left circumflex (LCX) artery and right coronary artery (RCA) and various global hemodynamics metrics. Based on our findings, the changes in coronary blood flow after TAVR varied and were subject specific (37% had increased flow in all three coronary arteries, 32% had decreased flow in all coronary arteries, and 31% had both increased and decreased flow in different coronary arteries). Additionally, valvular pressure gradient, left ventricle (LV) workload and maximum LV pressure decreased by 61.5%, 4.5% and 13.0% respectively, while mean arterial pressure and cardiac output increased by 6.9% and 9.9% after TAVR. By applying this proof-of-concept computational model, a series of hemodynamic metrics were generated non-invasively which can help to better understand the individual relationships between TAVR and mean and peak coronary flow rates. In the future, tools such as these may play a vital role by providing clinicians with rapid insight into various cardiac and coronary metrics, rendering the planning for TAVR and other cardiovascular procedures more personalized.

New Computational Approach to Shunt Design in Congenital Heart Palliation

Shunts are commonly used to redirect blood to pulmonary arteries in procedures that palliate congenital cardiovascular defects. Previous clinical studies and hemodynamic simulations reveal a critical role of shunt diameter in balancing flow to pulmonary versus systemic vessels, but the biomechanical process of creating the requisite anastomosis between the shunt and host vessel has received little attention. Here, we report a new Lagrange multiplier-based finite element approach that represents the shunt and host vessels as individual structures and predicts the anastomosis geometry and attachment force that result when the shunt is sutured at an incision in the host, followed by pressurization. Simulations suggest that anastomosis orifice opening increases markedly with increasing length of the host incision and moderately with increasing blood pressure. The host artery is further predicted to conform to common stiff synthetic shunts, whereas more compliant umbilical vessel shunts should conform to the host, with orifice area transitioning between these two extremes via a Hill-type function of shunt stiffness. Moreover, a direct relationship is expected between attachment forces and shunt stiffness. This new computational approach promises to aid in surgical planning for diverse vascular shunts by predicting in vivo pressurized geometries.

In-Silico and In-Vitro Analysis of the Novel Hybrid Comprehensive Stage II Operation for Single Ventricle Circulation

Single ventricle (SV) anomalies account for one-fourth of all congenital heart disease cases. The existing palliative treatment for this anomaly achieves a survival rate of only 50%. To reduce the trauma associated with surgical management, the hybrid comprehensive stage II (HCSII) operation was designed as an alternative for a select subset of SV patients with the adequate antegrade aortic flow. This study aims to provide better insight into the hemodynamics of HCSII patients utilizing a multiscale Computational Fluid Dynamics (CFD) model and a mock flow loop (MFL). Both 3D-0D loosely coupled CFD and MFL models have been tuned to match baseline hemodynamic parameters obtained from patient-specific catheterization data. The hemodynamic findings from clinical data closely match the in-vitro and in-silico measurements and show a strong correlation (r = 0.9). The geometrical modification applied to the models had little effect on the oxygen delivery. Similarly, the particle residence time study reveals that particles injected in the main pulmonary artery (MPA) have successfully ejected within one cardiac cycle, and no pathological flows were observed.

A Computational Study of Blood Flow Dynamics in the Pulmonary Arteries

In this work we study the blood dynamics in the pulmonary arteries by means of a 3D-0D geometric multiscale approach, where a detailed 3D model for the pulmonary arteries is coupled with a lumped parameters (0D) model of the cardiovascular system. We propose to investigate three strategies for the numerical solution of the 3D-0D coupled problem: the Splitting-Explicit and Implicit algorithms, where information are exchanged between 3D and 0D models at each time step at the interfaces, and the One-Way algorithm, where the 0D is solved first off-line. In our numerical experiments performed in a realistic patient-specific 3D domain with a physiologically calibrated 0D model, we discuss first the issue on instabilities that may arise when not suitable connections are considered between 3D and 0D models; second we compare the performance and accuracy of the three proposed numerical strategies. Finally, we report a comparison between a healthy and a hypertensive case, providing a preliminary result highlighting how our method could be used in future for clinical purposes.

Effects of Patent Ductus Arteriosus on the Hemodynamics of Modified Blalock-Taussig Shunt Based on Patient-Specific Simulation

The question of preserving the patent ductus arteriosus (PDA) during the modified Blalock–Taussig shunt (MBTS) procedure remains controversial. The goal of this study was to investigate the effects of the PDA on the flow features of the MBTS to help with preoperative surgery design and postoperative prediction. In this study, a patient with pulmonary atresia and PDA was included. A patient-specific three-dimensional model was reconstructed, and virtual surgeries of shunt insertion and ductus ligation were performed using computer-aided design. Computational fluid dynamics was utilized to analyze the hemodynamic parameters of varied models based on the patient-specific anatomy and physiological data. The preservation of the PDA competitively reduced the shunt flow but increased total pulmonary perfusion. The shunt flow and ductal flow collided, causing significant and complicated turbulence in the pulmonary artery where low wall shear stress, high oscillatory shear index, and high relative residence time were distributed. The highest energy loss was found when the PDA was preserved. The preservation of PDA is not recommended during MBTS procedures because it negatively influences hemodynamics. This may lead to pulmonary overperfusion, inadequate systemic perfusion, and a heavier cardiac burden, thus increasing the risk of heart failure. Also, it seems to bring no benefit in terms of reducing the risk for thrombosis.

Computational fluid dynamics investigation of the novel hybrid comprehensive stage II operation

Background

The hybrid comprehensive stage 2 (HCS2) procedure is a novel palliative operation applicable to a select subset of single ventricle patients with adequate native antegrade aortic flow to the upper body. Flow to the descending aorta, through the pulmonary outlet and ductal arch, is influenced by a stented intrapulmonary baffle connecting the branch pulmonary arteries. We used computational fluid dynamics (CFD) to elucidate the hemodynamic characteristics of this reconstruction.

Methods

We used multiscale CFD analysis of a synthetic, patient-derived HCS2 anatomic configuration with unsteady laminar flow conditions and a non-Newtonian blood model to quantify the resultant hemodynamics. The 3-dimensional CFD model was coupled to a 0-dimensional lumped parameter model of the peripheral circulation to determine the required boundary conditions.

Results

For the specific anatomy studied, the intrapulmonary baffle did not obstruct flow from the pulmonary trunk to ductal arch as long as the distance between the anterior pulmonary artery wall and baffle wall exceeded ∼7 mm. Vortex shedding off of the baffle wall did not develop, because of the short distance to the ductal arch. The stented baffle experienced significantly uneven "inward" loading from the systemic side. Pulmonary outlet flow separation distal to the baffle produced a low-speed recirculation region.

Conclusions

Hemodynamic patterns in this complex anatomy are generally favorable. Low flow recirculation could be mitigated by preoperative shape optimization. Calculated inward stresses on the pulmonary baffle can be used in the future to study baffle stent deformation, which is expected to be small.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.

A summary of second systemic pulmonary shunt for congenital heart disease with pulmonary hypoxemia

BACKGROUND

There has been an increasing number of children with congenital heart disease that undergo primary or second systemic-pulmonary shunt, while there are few reports on the second systemic-pulmonary shunt. Therefore, this study summarizes the experience of second systemic-pulmonary shunt for congenital heart disease in our hospital.

METHODS AND RESULTS

Sixty-five children with congenital heart disease who underwent systemic-pulmonary shunt for the second time in our hospital were analyzed. At the early stage after the operation, cyanosis improved and SpO2 significantly increased. One patient died in hospital (1.54%) and the causes of death were aggravated atrioventricular regurgitation, low cardiac output syndrome, and liver failure. Early complications occurred in 18 patients (27.7%). All the children were rechecked in our hospital every 3-6 months and the McGoon index significantly increased.

CONCLUSION

Systemic-pulmonary artery shunt can promote pulmonary vascular development, improve cyanosis symptoms, and increase the chance of radical treatment in children with pulmonary vascular dysplasia.

3D Simulation Analysis of Central Shunt in Patient-Specific Hemodynamics: Effects of Varying Degree of Pulmonary Artery Stenosis and Shunt Diameters

The objective of this study was to compare the effects of different shunt diameters and pulmonary artery (PA) stenosis grades on the hemodynamics of central shunts to determine an optimal surgical plan and improve the long-term outcomes of the operation. A 3D anatomical model was reconstructed based on the patient's clinical CT data. 3D computational fluid dynamics models were built with varying degrees of stenosis (the stenosis ratio α was represented by the ratio of blood flow through the main pulmonary artery to cardiac output, ranging from 0 to 30%; the smaller the value of α, the more severe the pulmonary artery stenosis) and varying shunt diameters (3, 3.5, 4, 4.5, and 5 mm). Our results show that the asymmetry of pulmonary artery flow increased with increasing shunt diameter and α, which will be more conducive to the development of the left pulmonary artery. Additionally, the pulmonary-to-systemic flow ratio (Q_P/Q_S) increases with the shunt diameter and α, and all the values exceed 1. When the shunt diameter is 3 mm and α = 0%, Q_P/Q_S reaches the minimum value of 1.01, and the oxygen delivery reaches the maximum value of 205.19 ml/min. However, increasing shunt diameter and α is beneficial to reduced power loss and smoother PA flow. In short, for patients with severe PA stenosis (α is small), a larger-diameter shunt may be preferred. Conversely, when the degree of PA stenosis is moderate, a smaller shunt diameter can be considered.