A high-fidelity geometric multiscale hemodynamic model for predicting myocardial ischemia

BACKGROUND AND OBJECTIVES

Coronary computed tomography angiography (CCTA) derived fractional flow reserve (CT-FFR) requires a maximal hyperemic state to be modeled by assuming the total coronary resistance decreased to a constant 0.24 of that under the resting state. However, this assumption neglects the vasodilator capacity of individual patients. Herein, we proposed a high-fidelity geometric multiscale model (HFMM) to characterize coronary pressure and flow under the resting state, seeking to better predict myocardial ischemia by using CCTA-derived instantaneous wave-free ratio (CT-iFR).

METHODS

Fifty-seven patients (62 lesions) who had undergone CCTA and were then referred to invasive FFR were prospectively enrolled. The coronary microcirculation resistance hemodynamic model (RHM) under the resting condition was established on a patient-specific basis. Coupled with a closed-loop geometric multiscale model (CGM) of their individual coronary circulations, the HFMM model was established to non-invasively derive the CT-iFR from CCTA images.

RESULTS

With the invasive FFR being the reference standard, accuracy of the obtained CT-iFR in identifying myocardial ischemia was greater than those of the CCTA and non-invasively derived CT-FFR (90.32% vs. 79.03% vs. 84.3%). The overall computational time of CT-iFR was 61 ± 6 min, faster than that of the CT-FFR (8 h). The sensitivity, specificity, positive predictive value, and negative predictive value of the CT-iFR in discriminating an invasive FFR > 0.8 were 78% (95% CI: 40-97%), 92% (95% CI: 82-98%), 64% (95% CI: 39-83%), and 96% (95% CI:88-99%), respectively.

CONCLUSIONS

A high-fidelity geometric multiscale hemodynamic model was developed for rapid and accurate estimation of CT-iFR. Compared with CT-FFR, CT-iFR is of less computational cost and enables assessment of tandem lesions.

Accurate Calculation of FFR Based on a Physics-Driven Fluid‐Structure Interaction Model

Background: The conventional FFRct numerical calculation method uses a model with a multi-scale geometry based upon CFD, and rigid walls. Therefore, important interactions between the elastic vessel wall and blood flow are not routinely considered. Changes in the resistance of coronary microcirculation during hyperaemia are likewise not typically incorporated using a fluid–structure interaction (FSI) algorithm. It is likely that both have resulted in FFRct calculation errors.

Objective: In this study we incorporated both the influence of vascular elasticity and coronary microcirculatory structure on FFR, to improve the accuracy of FFRct calculation. Thus, in this study, a physics-driven 3D–0D coupled model including fluid–structure interaction was established to calculate accurate FFRct values.

Methods: Based upon a novel geometric multi-scale modeling technology, a FSI simulation approach was used. A lumped parameter model (0D) was used as the outlet boundary condition for the 3D FSI coronary artery model to incorporate physiological microcirculation, with bidirectional coupling between the two models.

Results: The accuracy, sensitivity, specificity, and both positive and negative predictive values of FFR_DC calculated based upon the coupled 3D–0D model were 86.7, 66.7, 84.6, 66.7, and 91.7%, respectively. Compared to the calculated value using the basic CFD model (MSE = 5.9%, accuracy rate = 80%), the FFR_CFD calculated based on the coupled 3D–0D model has a smaller MSE of 1.9%.

Conclusion: The physics-driven coupled 3D–0D model that incorporates fluid–structure interactions not only consider the influence of the elastic vessel wall on blood flow, but also provides reliable microvascular resistance boundary conditions for the 3D FSI model. This allows for a calculation that is based upon conditions that are closer to the physiological environment, and thus improves the accuracy of FFRct calculation. It is likely that more accurate information will provide an enhanced recommendation regarding percutaneous coronary intervention (PCI) in the clinic.

The Effects of the Mechanical Properties of Vascular Grafts and an Anisotropic Hyperelastic Aortic Model on Local Hemodynamics during Modified Blalock-Taussig Shunt Operation, Assessed Using FSI Simulation

Cardiovascular surgery requires the use of state-of-the-art artificial materials. For example, microporous polytetrafluoroethylene grafts manufactured by Gore-Tex^® are used for the treatment of cyanotic heart defects (i.e., modified Blalock–Taussig shunt). Significant mortality during this palliative operation has led surgeons to adopt mathematical models to eliminate complications by performing fluid–solid interaction (FSI) simulations. To proceed with FSI modeling, it is necessary to know either the mechanical properties of the aorta and graft or the rheological properties of blood. The properties of the aorta and blood can be found in the literature, but there are no data about the mechanical properties of Gore-Tex^® grafts. Experimental studies were carried out on the mechanical properties vascular grafts adopted for modified pediatric Blalock–Taussig shunts. Parameters of two models (the five-parameter Mooney–Rivlin model and the three-parameter Yeoh model) were determined by uniaxial experimental curve fitting. The obtained data were used for patient-specific FSI modeling of local blood flow in the “aorta-modified Blalock–Taussig shunt–pulmonary artery” system in three different shunt locations: central, right, and left. The anisotropic model of the aortic material showed higher stress values at the peak moment of systole, which may be a key factor determining the strength characteristics of the aorta and pulmonary artery. Additionally, this mechanical parameter is important when installing a central shunt, since it is in the area of the central anastomosis that an increase in stress on the aortic wall is observed. According to computations, the anisotropic model shows smaller values for the displacements of both the aorta and the shunt, which in turn may affect the success of preoperative predictions. Thus, it can be concluded that the anisotropic properties of the aorta play an important role in preoperative modeling.

Multiscale modeling of a modified Blalock-Taussig surgery in a patient-specific tetralogy of Fallot

Tetralogy of Fallot (TOF) is a congenital heart anomaly that causes a drastic reduction in the oxygen level. In this study, we coupled a lumped-parameter model with a patient-specific three-dimensional (3D) model which included a modified Blalock-Taussig (MBT) shunt. By forming a closed loop, we investigated the effects of certain parameters on the flow rates and the pressures at different locations of the developed network. A local sensitivity analysis on an initial zero-dimensional (0D) closed-loop model was conducted. The 0D lumped parameter (LP) model was then refined based on the results of the multiscale 0D-3D model and the local sensitivity analysis was repeated for the refined 0D model. It was shown that the maximum pressure of the pulmonary bed had the highest sensitivity of 94% to the diameter of MBT shunt. We observed that the existence of the flow in the shunt during the diastole caused an elevated wall shear stress (WSS) in the pulmonary artery. In this work, we calculated the flow velocity and pressure field in a 3D patient-specific aorta with an MBT shunt, and then we used the results to increase the accuracy of our LP model to simulate numerous 0D simulations in a significantly shorter time, which is potentially applicable for medical decision-making.

Numerical Simulation of the Effect of Pulmonary Vascular Resistance on the Hemodynamics of Reoperation After Failure of One and a Half Ventricle Repair

OBJECTIVE

The one and a half ventricle repair (1.5VR) is a common clinical choice for patients with right heart dysfunction. Considering the influence of blood circulation failure and reoperation in urgent need, this essay aims to explore the hemodynamic effects of different pulmonary vascular resistance (PVR) values on reoperation after 1.5VR failure.

METHODS

The lumped parameter model (LPM) was used to simulate the reoperation, including the return biventricular repair (2VR), ligation of azygos vein (1.5VR') and return single ventricular repair (1.0VR). Firstly, the debugging parameters were used to simulate the hemodynamics of 2VR. Secondly, the value of PVR was changed from one to four times while the other parameters remained unchanged. Finally, 15 cardiac cycles were simulated and the 15th result was obtained. In this work, the left and right ventricular stroke work and their sum (Plv, Prv, Ptotal), the left and right ventricular ejection fraction (LVEF, RVEF), the mean Cardiac Output (mCO) and the mean pressure and flow-rate ratio of superior and inferior vena cava (mPsvc\mPivc and mQsvc\mQivc), respectively, were used to describe the hemodynamics of reoperation.

RESULTS

With the change of PVR from one to four times, the values of Plv, Prv, Ptotal, LVEF, and RVEF gradually decreased. The change rate of Plv, Ptotal and LVEF of 1.0VR were the largest in the three kinds of reoperation. The change rate of Prv of 1.5VR' was larger than that of 2VR, but it was the opposite for their EF change rate. The mCO of 2VR, 1.5VR', and 1.0VR decreased by 18.53%, 37.58%, and 48.07%, respectively. The mPsvc\mPivc of 1.5VR' increased from 3.76 to 6.77 and the mQsvc\mQivc decreased from 0.55 to 0.36, while the mPsvc\mPivc and mQsvc\mQivc of 2VR and 1.0VR remained 1 and 0.67, respectively. The peak value of the tricuspid flow-rate (Qti) waveform of 2VR and 1.5VR' changed from "E peak" to "A peak."

CONCLUSION

The numerical results demonstrate the highly reoperation-dependent hemodynamic consequences and their responses to variations in PVR. Comprehensive analysis of EF, mCO and ventricular stroke work indicates that PVR has a greater impact on 1.5VR' and 1.0VR. Therefore, we suggest that the selection strategy of reoperation should focus on PVR.

HEMODYNAMICS-BASED LONG-TERM PATENCY OF DIFFERENT SEQUENTIAL GRAFTING: A PATIENT-SPECIFIC MULTI-SCALE STUDY

Background and aims: Sequential grafting is one of the common coronary artery bypass grafting (CABG) surgery. But the influence of the sequential grafting position on hemodynamics and the graft patency is still unclear. Materials and methods: The zero-dimensional/three-dimensional (0D/3D) coupling method was used to finalize the multi-scale simulation of two different sequential grafting models. First, a patient-specific 3D model was reconstructed based on coronary computed tomography angiography (CCTA) images. Two different sequential grafts were implemented on this patient-specific 3D model by using virtual surgery. Thus, two different postoperative 3D models were built. Then, a lumped parameter model (LPM; 0D) was built based on the patient physiological data to simulate the cardiovascular system. Finally, the 0D/3D coupling method was used to perform the numerical simulation by coupling a 0D LPM of the cardiovascular system and the patient-specific 3D models. Moreover, the long-term patency of these two different sequential grafts was discussed in this paper. Results: The coronary flow rate and the graft flow were calculated and illustrated. The instantaneous wave-free ratio (iFR) were calculated. Postoperative iFR values increase to over 0.90 for both sequential grafts. Some hemodynamics parameters were also illustrated, such as wall shear stress (WSS), oscillatory shear index (OSI). The area of low WSS in Model 1 was much less than that in Model 2. Two regions of high OSI exist in Model 2, while only one in Model 1. Conclusions: No significant differences exist on the short-term outcomes of two models. But the long-term patency of Model 2 was worse. The Model 1 may enhance long-term patency of grafting and should be priority when the sequential grafting need to be carried out.

SURGICAL DECISION OF CORONARY ARTERY BYPASS GRAFTING FOR NORMAL LEFT ANTERIOR DESCENDING ARTERY (LAD) AND LAD WITH STENOSIS: SEQUENTIAL GRAFT OR NOT

Sequential graft was used frequently in clinical studies. In this study, the hemodynamic effect of one kind of sequential graft in two different conditions of the lesion was discussed and some recommendations on the surgical procedures were made. A patient-specific three-dimensional (3D) model of left anterior descending artery (LAD) was reconstructed. A moderate stenosis exist in the trunk of LAD between the first and the second diagonal branch (D1 and D2). Another 3D model without stenosis was also reconstructed based on the patient-specific model. Sequential graft and single graft were applied on these two 3D model. Thus four 3D models were built so that the hemodynamic effect of sequential graft can be discussed. The zero-dimensional (0D)/3D coupling method was used to perform the numerical simulation by coupling the 3D artery model with a 0D lumped parameter model of the cardiovascular system. The flow rates in the branches of LAD and the graft flow were calculated and illustrated in this paper. The wall shear stress (WSS) and oscillatory shear index (OSI) were also calculated and depicted. If the native LAD is stenosis, sequential graft should be applied for the short-term outcomes. Moreover, the long-term patency of the sequential graft applied on the stenosis LAD is good. The long-term patency of the single graft was bad. But the short-term outcomes are almost the same when LAD is not stenosis. If no stenosis exist in the native LAD, a graft with smaller diameter should be applied to improve the long-term patency.

DIRECT CORONARY COUPLING APPROACH FOR COMPUTING FFRCT

With the advances in computational fluid dynamics (CFD) and image-based modeling techniques, fractional flow reserve (FFR) can be computed from coronary computed tomography angiography (CTA) scans (FFRCT). However, this non-invasive approach requires large-scale computational resources, which limits its application in routine clinical setting. A 3D–0D coupling approach is proposed to improve the coupling efficiency of FFRCT. Aortic–root is modeled by a lumped parameter model and connected with the models of left ventricle and systemic circulation. With this approach, the interested coronary regions can be directly coupled to the lumped parameter model, resulting in a significant reduction (up to 20 times reduction) in the volume of the CFD computing domain. The proposed approach is applied to a patient-specific model and compared with previous non-reduced approach. Results show that the computed coronary flow rates, pressure waveforms and FFRCT contours by the proposed approach are consistent well with that of the non-reduced approach. These results demonstrate that the proposed approach can reduce the CFD computing domain of FFRCT significantly while maintaining the similar accuracy as compared with the non-reduced approach, and it can be further employed to promote FFRCT in routine clinical setting.