New imaging tools in cardiovascular medicine: computational fluid dynamics and 4D flow MRI

Establishing the distribution of cerebrovascular resistance using computational fluid dynamics and 4D flow MRI

Cerebrovascular resistance (CVR) regulates blood flow in the brain, but little is known about the vascular resistances of the individual cerebral territories. We present a method to calculate these resistances and investigate how CVR varies in the hemodynamically disturbed brain. We included 48 patients with stroke/TIA (29 with symptomatic carotid stenosis). By combining flow rate (4D flow MRI) and structural computed tomography angiography (CTA) data with computational fluid dynamics (CFD) we computed the perfusion pressures out from the circle of Willis, with which CVR of the MCA, ACA, and PCA territories was estimated. 56 controls were included for comparison of total CVR (tCVR). CVR were 33.8 ± 10.5, 59.0 ± 30.6, and 77.8 ± 21.3 mmHg s/ml for the MCA, ACA, and PCA territories. We found no differences in tCVR between patients, 9.3 ± 1.9 mmHg s/ml, and controls, 9.3 ± 2.0 mmHg s/ml (p = 0.88), nor in territorial CVR in the carotid stenosis patients between ipsilateral and contralateral hemispheres. Territorial resistance associated inversely to territorial brain volume (p < 0.001). These resistances may work as reference values when modelling blood flow in the circle of Willis, and the method can be used when there is need for subject-specific analysis.

4D flow MRI-derived energy loss and RV workload in adults with tetralogy of Fallot

PURPOSE

To assess flow energy loss (EL) pattern inside the pulmonary circulation in adult patients with repaired tetralogy of Fallot (TOF), particularly in TOF with pulmonary stenosis (PS) and pulmonary regurgitation (PR), as a cardiac workload parameter and its relationship to symptoms and major adverse cardiovascular events (MACE).

METHODS

Prospectively, 51 consecutive TOF adults after intracardiac repair, who underwent four-dimensional flow magnetic resonance imaging, were enrolled. All of them had significant PR (PR regurgitant fraction >25 %). TOF patients who had already reached the conventional criteria were excluded. We defined MACE as the following: 1) fatal arrhythmias, 2) sudden cardiac death, 3) surgical pulmonary valvular repair (PVR), 4) right heart failure (HF) needing diuretics and/or hospitalization within 2 years.

RESULTS

A total of 15 patients had MACE; 1) 10 patients underwent PVR within 2 years, 2) 2 patients had ventricular tachycardia, and 3) 6 patients developed right HF (overlapped). Right ventricular (RV) end diastolic volume index (EDVI), RV end systolic volume index (ESVI), average EL/cardiac output (CO), and diastolic EL/CO in patients with MACE were greater than ones without MACE. On a multivariate logistic analysis, the diastolic EL/CO ratio and RVEDVI had the highest odds with MACE in all TOF (odds ratio, 40.7 and 1.15. 95%CI, 1.83-905 and 1.02-13.0; p-value, 0.02 and 0.03). In sub-analysis within 29 patients with moderate PS with PR, and 10 patients with MACE showed higher diastolic EL/CO. Average and diastolic EL/CO negatively correlated with RV ejection fraction (EF) in this sub-analysis.

CONCLUSIONS

High EL, particularly, high diastolic EL/CO, were the important factors for MACE in adult TOF. Higher diastolic EL/CO was also related to lower RV EF and deteriorated RV function in adult TOF with PS and PR. Right-sided EL can be a sensitive marker of excessive cardiac workload which integrates both afterload and preload in adult patients with TOF, beyond the RV size.

Blood flow analysis with computational fluid dynamics in the left atrium after left atrial plication: a prospective study

OBJECTIVE

This study aimed to evaluate blood flow stagnation in an enlarged left atrium (LA) and prove that left atrial plication (LAP) could alleviate the stagnation.

METHODS

Five patients with chronic atrial fibrillation who underwent mitral valve surgery followed by LAP for an enlarged LA with a ≥ 60-mm diameter were included. We performed computational fluid dynamics (CFD) analysis using preoperative and postoperative computed tomography and four-dimensional flow magnetic resonance imaging. Additionally, computer graphics were used to create virtual left atrial appendage resection (LAAR) images. We performed CFD analysis to assess blood flow stagnation in the LA for three groups: preoperative, LAAR, and LAP.

RESULTS

When the average and constant stagnation volumes were both set to 100 preoperatively, the average stagnation volumes of the LAAR and LAP groups were 67.42 ± 18.64 and 35.88 ± 8.20, respectively. The constant stagnation volumes of these groups reduced to 45.01 ± 7.43 and 21.14 ± 7.70, respectively. The LAP group also had significantly lower average and constant stagnation volumes than those in the LAAR group (p = 0.006 and p = 0.033, respectively).

CONCLUSIONS

Blood flow stagnation was noted in the LAA and enlarged LA. CFD analysis revealed that LAP for the enlarged LA improved blood flow stagnation more than the virtual LAAR alone.

CLINICAL TRIAL REGISTRY NUMBER

UMIN000049923.

Recent Advances in Cardiac Magnetic Resonance for Imaging of Acute Myocardial Infarction

Acute myocardial infarction (AMI) is one of the primary causes of death worldwide, with a high incidence and mortality rate. Assessment of the infarcted and surviving myocardium, along with microvascular obstruction, is crucial for risk stratification, treatment, and prognosis in patients with AMI. Nonionizing radiation, excellent soft tissue contrast resolution, a large field of view, and multiplane imaging make cardiac magnetic resonance (CMR) a "one-stop" method for assessing cardiac structure, function, perfusion, and metabolism. Hence, this imaging technology is considered the "gold standard" for evaluating myocardial function and viability in AMI. This review critically compares the advantages and disadvantages of CMR with other cardiac imaging technologies, and relates the imaging findings to the underlying pathophysiological processes in AMI. A more thorough understanding of CMR technology will clarify their advanced clinical diagnosis and prognostic assessment applications, and assess the future approaches and challenges of CMR in the setting of AMI.

Advances in nanobased platforms for cardiovascular diseases: Early diagnosis, imaging, treatment, and tissue engineering

Cardiovascular diseases (CVDs) present a significant threat to health, with traditional therapeutics based treatment being hindered by inefficiencies, limited biological effects, and resistance to conventional drug. Addressing these challenges requires advanced approaches for early disease diagnosis and therapy. Nanotechnology and nanomedicine have emerged as promising avenues for personalized CVD diagnosis and treatment through theranostic agents. Nanoparticles serve as nanodevices or nanocarriers, efficiently transporting drugs to injury sites. These nanocarriers offer the potential for precise drug and gene delivery, overcoming issues like bioavailability and solubility. By attaching specific target molecules to nanoparticle surfaces, controlled drug release to targeted areas becomes feasible. In the field of cardiology, nanoplatforms have gained popularity due to their attributes, such as passive or active targeting of cardiac tissues, enhanced sensitivity and specificity, and easy penetration into heart and artery tissues due to their small size. However, concerns persist about the immunogenicity and cytotoxicity of nanomaterials, necessitating careful consideration. Nanoparticles also hold promise for CVD diagnosis and imaging, enabling straightforward diagnostic procedures and real-time tracking during therapy. Nanotechnology has revolutionized cardiovascular imaging, yielding multimodal and multifunctional vehicles that outperform traditional methods. The paper provides an overview of nanomaterial delivery routes, targeting techniques, and recent advances in treating, diagnosing, and engineering tissues for CVDs. It also discusses the future potential of nanomaterials in CVDs, including theranostics, aiming to enhance cardiovascular treatment in clinical practice. Ultimately, refining nanocarriers and delivery methods has the potential to enhance treatment effectiveness, minimize side effects, and improve patients' well-being and outcomes.

Differences in blood flow dynamics between balloon- and self-expandable valves in patients with aortic stenosis undergoing transcatheter aortic valve replacement

BACKGROUND

The differences in pre- and early post-procedural blood flow dynamics between the two major types of bioprosthetic valves, the balloon-expandable valve (BEV) and self-expandable valve (SEV), in patients with aortic stenosis (AS) undergoing transcatheter aortic valve replacement (TAVR), have not been investigated. We aimed to investigate the differences in blood flow dynamics between the BEV and SEV using four-dimensional flow cardiovascular magnetic resonance (4D flow CMR).

METHODS

We prospectively examined 98 consecutive patients with severe AS who underwent TAVR between May 2018 and November 2021 (58 BEV and 40 SEV) after excluding those without CMR because of a contraindication, inadequate imaging from the analyses, or patients' refusal. CMR was performed in all participants before (median interval, 22 [interquartile range (IQR) 4-39] days) and after (median interval, 6 [IQR 3-6] days) TAVR. We compared the changes in blood flow patterns, wall shear stress (WSS), and energy loss (EL) in the ascending aorta (AAo) between the BEV and SEV using 4D flow CMR.

RESULTS

The absolute reductions in helical flow and flow eccentricity were significantly higher in the SEV group compared in the BEV group after TAVR (BEV: - 0.22 ± 0.86 vs. SEV: - 0.85 ± 0.80, P < 0.001 and BEV: - 0.11 ± 0.79 vs. SEV: - 0.50 ± 0.88, P = 0.037, respectively); there were no significant differences in vortical flow between the groups. The absolute reduction of average WSS was significantly higher in the SEV group compared to the BEV group after TAVR (BEV: - 0.6 [- 2.1 to 0.5] Pa vs. SEV: - 1.8 [- 3.5 to - 0.8] Pa, P = 0.006). The systolic EL in the AAo significantly decreased after TAVR in both the groups, while the absolute reduction was comparable between the groups.

CONCLUSIONS

Helical flow, flow eccentricity, and average WSS in the AAo were significantly decreased after SEV implantation compared to BEV implantation, providing functional insights for valve selection in patients with AS undergoing TAVR. Our findings offer valuable insights into blood flow dynamics, aiding in the selection of valves for patients with AS undergoing TAVR. Further larger-scale studies are warranted to confirm the prognostic significance of hemodynamic changes in these patients.

Mechanism Analysis of Vascular Calcification Based on Fluid Dynamics

Vascular calcification is the abnormal deposition of calcium phosphate complexes in blood vessels, which is regarded as the pathological basis of multiple cardiovascular diseases. The flowing blood exerts a frictional force called shear stress on the vascular wall. Blood vessels have different hydrodynamic properties due to discrepancies in geometric and mechanical properties. The disturbance of the blood flow in the bending area and the branch point of the arterial tree produces a shear stress lower than the physiological magnitude of the laminar shear stress, which can induce the occurrence of vascular calcification. Endothelial cells sense the fluid dynamics of blood and transmit electrical and chemical signals to the full-thickness of blood vessels. Through crosstalk with endothelial cells, smooth muscle cells trigger osteogenic transformation, involved in mediating vascular intima and media calcification. In addition, based on the detection of fluid dynamics parameters, emerging imaging technologies such as 4D Flow MRI and computational fluid dynamics have greatly improved the early diagnosis ability of cardiovascular diseases, showing extremely high clinical application prospects.

Intracardiac Flow Analysis of the Right Ventricle in Pediatric Patients With Repaired Tetralogy of Fallot Using a Novel Color Doppler Velocity Reconstruction

BACKGROUND

Repaired tetralogy of Fallot (RTOF) patients will develop right ventricular (RV) dysfunction from chronic pulmonary regurgitation (PR). Cardiac magnetic resonance sequences such as four-dimensional flow can demonstrate altered vorticity and flow energy loss (FEL); however, they are not as available as conventional echocardiography (echo). The study determined whether a novel, vendor-independent Doppler velocity reconstruction (DoVeR) could measure RV intracardiac flow in conventional echo of RTOF patients. The primary hypothesis was that DoVeR could detect increased vorticity and diastolic FEL in RTOF patients.

METHODS

Repaired tetralogy of Fallot patients with echo were retrospectively paired with age-/size-matched controls. Doppler velocity reconstruction employed the stream function-vorticity equation to approximate intracardiac flow fields from color Doppler. A velocity field of the right ventricle was reconstructed from the apical 4-chamber view. Vortex strength (VS, area integral of vorticity) and FEL were derived from DoVeR. Cardiac magnetic resonance and exercise stress parameters (performed within 1 year of echo) were collected for analysis.

RESULTS

Twenty RTOF patients and age-matched controls were included in the study. Mean regurgitant fraction was 40.5% ± 7.6%, and indexed RV end-diastolic volume was 158 ± 36 mL/m2. Repaired tetralogy of Fallot patients had higher total, mean diastolic, and peak diastolic VS (P = .0013, P = .0012, P = .0032, respectively) and higher total, mean diastolic, and peak diastolic body surface area-indexed FEL (P = .0016, P = .0022, P < .001, respectively). Peak diastolic indexed FEL and peak diastolic VS had weak-to-moderate negative correlation with RV ejection fraction (r = -0.52 [P = .019] and r = -0.49 [P = .030], respectively) and left ventricular ejection fraction (r = -0.47 [P = .034] and r = -0.64 [P = .002], respectively). Mean diastolic indexed FEL and VS had moderate-to-strong negative correlation with percent predicted maximal oxygen consumption (r = -0.69 [P = .012] and r = -0.75 [P = .006], respectively).

CONCLUSIONS

DoVeR can detect alterations to intracardiac flow in RTOF patients from conventional color Doppler imaging. Echo-based measures of diastolic VS and FEL correlated with ventricular function. DoVeR has the potential to provide serial evaluation of abnormal flow dynamics in RTOF patients.

Super-resolution 4D flow MRI to quantify aortic regurgitation using computational fluid dynamics and deep learning

Changes in cardiovascular hemodynamics are closely related to the development of aortic regurgitation (AR), a type of valvular heart disease. Metrics derived from blood flows are used to indicate AR onset and evaluate its severity. These metrics can be non-invasively obtained using four-dimensional (4D) flow magnetic resonance imaging (MRI), where accuracy is primarily dependent on spatial resolution. However, insufficient resolution often results from limitations in 4D flow MRI and complex aortic regurgitation hemodynamics. To address this, computational fluid dynamics simulations were transformed into synthetic 4D flow MRI data and used to train a variety of neural networks. These networks generated super-resolution, full-field phase images with an upsample factor of 4. Results showed decreased velocity error, high structural similarity scores, and improved learning capabilities from previous work. Further validation was performed on two sets of in vivo 4D flow MRI data and demonstrated success in de-noising flow images. This approach presents an opportunity to comprehensively analyse AR hemodynamics in a non-invasive manner.

Pulmonary artery blood flow dynamics in chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension is caused by incomplete resolution and organization of thrombi. Blood flow dynamics are involved in thrombus formation; however, only a few studies have reported on pulmonary artery blood flow dynamics in patients with chronic thromboembolic pulmonary hypertension. Furthermore, the effects of treatment interventions on pulmonary artery blood flow dynamics are not fully understood. The aim of the study was to evaluate pulmonary artery blood flow dynamics in patients with chronic thromboembolic pulmonary hypertension before and after pulmonary endarterectomy and balloon pulmonary angioplasty, using computational fluid dynamics. We analyzed patient-specific pulmonary artery models of 10 patients with chronic thromboembolic pulmonary hypertension and three controls using computational fluid dynamics. In patients with chronic thromboembolic pulmonary hypertension, flow velocity and wall shear stress in the pulmonary arteries were significantly decreased, and the oscillatory shear index and blood stagnation volume were significantly increased than in controls. Pulmonary endarterectomy induced redistribution of pulmonary blood flow and improved blood flow dynamics in the pulmonary artery. Balloon pulmonary angioplasty improved pulmonary blood flow disturbance, decreased blood flow stagnation, and increased wall shear stress, leading to vasodilatation of the distal portion of the pulmonary artery following balloon pulmonary angioplasty treatment.

Intracardiac energy inefficiency during atrial fibrillation rhythm with heart failure: Assessment by echo-vector flow mapping

AIMS

Intracardiac dynamics during atrial fibrillation(AF) complicated by heart failure(HF) are not fully understood. The aim of this study was to evaluate the impact of intracardiac dynamics assessed by echo-vector flow mapping on AF complicated by HF.

METHODS AND RESULTS

We enrolled 76 AF patients receiving sinus restoration therapy and assessed energy loss(EL) by echo-vector flow mapping during both AF rhythm and sinus rhythm. Patients were divided into 2 groups according to serum NT-proBNP level: high NT-proBNP group(≥1800 pg/ml during AF rhythm: N = 19), and low NT-proBNP group(N = 57). Outcome measures were average ELs/stroke volume(SV) in left ventricle(LV) and left atrium(LA). Average EL/SVs during AF rhythm in the LV and LA were significantly larger in the high NT-proBNP group than the low NT-proBNP group(54.2 mE/m·L versus 41.2 mE/m·L, P = 0.02; 3.2 mE/m·L versus 1.9 mE/m·L, P = 0.01). The significantly larger EL/SV in the high NT-proBNP group was recorded for maximum EL/SVs. In patients with high NT-proBNP, large vortex formation with extreme EL was detected in the LV and LA during the diastolic phase. After sinus restoration, the average decrease of EL/SV in the LV and LA was larger in the high NT-proBNP group(-21.4 mE/m·L versus +2.6 mE/m·L, P = 0.04; -1.6 mE/m·L versus -0.3 mE/m·L, P = 0.02). Average EL/SV during sinus rhythm did not significantly differ between the high and low NT-proBNP groups in the LV and LA.

CONCLUSIONS

High EL during AF rhythm as intracardiac energy inefficiency was associated with high serum NT-proBNP levels and improved after sinus restoration.

Correlation between aortic valve protein levels and vector flow mapping of wall shear stress and oscillatory shear index in patients supported with continuous-flow left ventricular assist devices

BACKGROUND

Continuous-flow left ventricular assist devices commonly lead to aortic regurgitation, which results in decreased pump efficiency and worsening heart failure. We hypothesized that non-physiological wall shear stress and oscillatory shear index alter the abundance of structural proteins in aortic valves of left ventricular assist device (LVAD) patients.

METHODS

Doppler images of aortic valves of patients undergoing heart transplants were obtained. Eight patients had been supported with LVADs, whereas 10 were not. Aortic valve tissue was collected and protein levels were analyzed using mass spectrometry. Echocardiographic images were analyzed and wall shear stress and oscillatory shear index were calculated. The relationship between normalized levels of individual proteins and in vivo echocardiographic measurements was evaluated.

RESULTS

Of the 57 proteins of interest, there was a strong negative correlation between levels of 15 proteins and the wall shear stress (R < -0.500, p ≤ 0.05), and a moderate negative correlation between 16 proteins and wall shear stress (R -0.500 to -0.300, p ≤ 0.05). Gene ontology analysis demonstrated clusters of proteins involved in cellular structure. Proteins negatively correlated with WSS included those with cytoskeletal, actin/myosin, cell-cell junction and extracellular functions. C: In aortic valve tissue, 31 proteins were identified involved in cellular structure and extracellular junctions with a negative correlation between their levels and wall shear stress. These findings suggest an association between the forces acting on the aortic valve (AV) and leaflet protein abundance, and may form a mechanical basis for the increased risk of aortic leaflet degeneration in LVAD patients.

Changes in Internal Thoracic Artery Blood Flow According to the Degree of Stenosis of the Anterior Descending Branch of the Left Coronary Artery

PURPOSE

Computational fluid dynamics has enabled the evaluation of coronary flow reserve. The purpose of this study was to clarify the hemodynamic variation and reserve potential of the left internal thoracic artery (LITA).

METHODS

Four patients were selected on the basis of various native coronary stenosis patterns and graft design. The wall shear stress and oscillatory shear index were measured, and one patient was selected. Next, we created three hypothetical lesions with 75%, 90%, and 99% stenosis in front of the graft anastomosis, and compared the changes in LITA blood flow and coronary flow distribution.

RESULTS

In the 75% to 90% stenosis model, blood flow was significantly higher in the native coronary flow proximal to the coronary artery bypass anastomosis regardless of time phase. In the 99% stenosis model, blood flow from the LITA was significantly dominant compared to native coronary flow at the proximal site of anastomosis. The range of LITA flow variability was the largest at 99% stenosis, with a difference of 70 ml/min.

CONCLUSION

The 99% stenosis model showed the highest LITA flow. The range of LITA flow variability is large, suggesting that it may vary according to the rate of native coronary stenosis.

Subsecond Three-Dimensional Nitrogen-15 Magnetic Resonance Imaging Facilitated by Parahydrogen-Based Hyperpolarization

Magnetic resonance imaging (MRI) provides unique information about the internal structure and function of living organisms in a non-invasive way. The use of conventional proton MRI for the observation of real-time metabolism is hampered by the dominant signals of water and fat, which are abundant in living organisms. Heteronuclear MRI in conjunction with the hyperpolarization methods does not encounter this issue. In this work, we polarized 15N nuclei of [15N1]fampridine (a drug used for the treatment of multiple sclerosis) to the level of 4% in nuclear magnetic resonance (NMR) experiments and 0.7% in MRI studies using spin-lock-induced crossing combined with signal amplification by reversible exchange. Consequently, three-dimensional 15N MRI of the hyperpolarized 15N-labeled drug was acquired in 0.1 s with a signal-to-noise ratio of 70. In addition, the NMR signal enhancements for 15N-enriched fampridine and fampridine with a natural abundance of 15N nuclei were compared and an explanation for their difference was proposed.

An Oscillatory Shear Index-Based Model to Describe Progressive Carotid Artery Stenosis

Background and aims: This study describes and demonstrates the applicability of a novel in silico method for modeling progressive carotid artery stenosis using the oscillatory shear index (OSI) as the basis of stenosis. Methods: Three-dimensional reconstructions of 11 carotid arteries were generated using patient-derived magnetic resonance angiography and duplex ultrasound data. Computational fluid dynamic simulations were sequentially generated following computational stenosis assessment, and corresponding changes in OSI were observed and used as measure of morphological stabilization. Results: 6 carotid models showed progressive stenosis with statistically significant increases in regions of high OSI (OSI >.2, P < .05) with eventual carotid occlusion in 1 of the cases. Three models remained free or nearly free of increased OSI, whereas 1 model showed an overall decrease in high OSI regions (P < .05) and another trended in that direction but did not achieve statistical significance (P = .145). Conclusions: To our knowledge, this is the first computational model describing progressive stenosis in any peripheral artery including the carotid. Taken together, this study provides a novel framework for computational hemodynamic investigations on progressive atherosclerosis in the carotid artery.

The Effect of Ascending Aortic Repair on Left Ventricular Remodeling

Left ventricular (LV) hypertrophy is common in patients with thoracic aortic diseases and is associated with increased long-term mortality. Thoracic aortic aneurysms are reported to increase LV afterload because of kinetic energy loss within the aneurysm sac, which may improve after surgical repair. However, LV afterload may also increase because of the stiffness of prosthetics used for aortic repair. We sought to investigate the long-term effect of surgical aortic repair with prostheses on postsurgical LV mass. We reviewed patients who underwent ascending aortic replacement with a prosthesis at our institution from January 2008 to December 2018. We calculated the LV mass index based on pre- and postoperative echocardiogram measurements. The primary outcome was the change in LV mass index 6 months after aortic repair. Patients aged <18 years and those who had concomitant cardiac operations, severe aortic valve disease, or who had no echocardiographic data were excluded. Of 1,008 patients who underwent ascending aortic replacement, 134 (51 with acute aortic dissections) were included. The median baseline and follow-up LV mass index were 107 (90 to 135) g/m² and 101 (83 to 123) g/m², respectively. Overall, there was a significant reduction of LV mass index over time (p = 0.03). LV mass index decreased in 77 patients (59%). Presentation due to acute aortic dissection (p = 0.03) and baseline LV mass index (p <0.001) were significant predictors of LV mass reduction. In conclusion, LV mass index may significantly decrease over time after the aortic repair, but the course is highly variable. The largest decrease occurred in patients who presented because of aortic dissections rather than for elective repair of aneurysms.

Evaluation of pulmonary blood flow in bilateral bidirectional Glenn shunts: value of 4‐D flow cardiac magnetic resonance in the evaluation of pulmonary artery confluence stenosis

Bilateral bidirectional Glenn shunts are associated with the risk of developing pulmonary artery bifurcation stenosis, resulting in variable pulmonary blood flow to either lung. This could negatively impact the subsequent stages of the single ventricle palliation pathway. This report highlights the value of 4D flow sequence from the cardiac magnetic resonance imaging in demonstrating the pulmonary blood flow characteristics following a bilateral bidirectional Glenn procedure. Mapping the blood flow pattern and its quantification to each lung provide objective insights into the possible predisposing factors for the development of pulmonary bifurcation stenosis.

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

Hemodynamic Parameters for Cardiovascular System in 4D Flow MRI: Mathematical Definition and Clinical Applications

Blood flow imaging becomes an emerging trend in cardiology with the recent progress in computer technology. It not only visualizes colorful flow velocity streamlines but also quantifies the mechanical stress on cardiovascular structures; thus, it can provide the detailed inspections of the pathophysiology of diseases and predict the prognosis of cardiovascular functions. Clinical applications include the comprehensive assessment of hemodynamics and cardiac functions in echocardiography vector flow mapping (VFM), 4D flow MRI, and surgical planning as a simulation medicine in computational fluid dynamics (CFD).For evaluation of the hemodynamics, novel mathematically derived parameters obtained using measured velocity distributions are essential. Among them, the traditional and typical parameters are wall shear stress (WSS) and its related parameters. These parameters indicate the mechanical damages to endothelial cells, resulting in degenerative intimal change in vascular diseases. Apart from WSS, there are abundant parameters that describe the strength of the vortical and/or helical flow patterns. For instance, vorticity, enstrophy, and circulation indicate the rotating flow strength or power of 2D vortical flows. In addition, helicity, which is defined as the cross-linking number of the vortex filaments, indicates the 3D helical flow strength and adequately describes the turbulent flow in the aortic root in cases with complicated anatomies. For the description of turbulence caused by the diseased flow, there exist two types of parameters based on completely different concepts, namely: energy loss (EL) and turbulent kinetic energy (TKE). EL is the dissipated energy with blood viscosity and evaluates the cardiac workload related to the prognosis of heart failure. TKE describes the fluctuation in kinetic energy during turbulence, which describes the severity of the diseases that cause jet flow. These parameters are based on intuitive and clear physiological concepts, and are suitable for in vivo flow measurements using inner velocity profiles.

Does the Pulsatile Non-uniform Flow Matter in MR Flowmetry?

3D cine phase-contrast (4D flow) MRI is a sequence with great potential for non-invasive time-resolved 3D flowmetry at arbitrary vessel sections in various blood vessels. However, it is not widely known that the flowmetry with 4D flow MRI is vulnerable to pulsatile and non-uniform flow. Due to the limited spatial and temporal resolutions, averaging within the 3D voxel is occurring during the flowmetry. A simple solution is to avoid setting the measurement plane in the area where non-uniform flow is dominant, which is possible with an aid of streamline depictions generated by computational fluid dynamics (CFD) or 4D flow MRI data. Unlike 4D flow MRI, flowmetry in CFD simulation can use higher spatial and temporal resolution depending on computer performance; therefore, it is robust to fluctuating non-uniform flow. However, the performance of CFD simulations might be limited due to inlet conditions with low temporal resolution. Difficulty applying complex blood flow such as reflection flow from periphery may also limit accurate simulation. Caution should be taken when comparing the result of CFD simulation to that of 4D flow measurement.

Computational Modeling of Right Ventricular Motion and Intracardiac Flow in Repaired Tetralogy of Fallot

PURPOSE

Patients with repaired Tetralogy of Fallot (rTOF) will develop dilation of the right ventricle (RV) from chronic pulmonary insufficiency and require pulmonary valve replacement (PVR). Cardiac MRI (cMRI) is used to guide therapy but has limitations in studying novel intracardiac flow parameters. This pilot study aimed to demonstrate feasibility of reconstructing RV motion and simulating intracardiac flow in rTOF patients, exclusively using conventional cMRI and an immersed-boundary method computational fluid dynamic (CFD) solver.

METHODS

Four rTOF patients and three normal controls underwent cMRI including 4D flow. 3D RV models were segmented from cMRI images. Feature-tracking software captured RV endocardial contours from cMRI long-axis and short-axis cine stacks. RV motion was reconstructed via diffeomorphic mapping (Deformetrica, deformetrica.org), serving as the domain boundary for CFD. Fully-resolved direct numerical simulations were performed over several cardiac cycles. Intracardiac vorticity, kinetic energy (KE) and turbulent kinetic energy (TKE) was measured. For validation, RV motion was compared to manual tracings, results of KE were compared between CFD and 4D flow.

RESULTS

Diastolic vorticity and TKE in rTOF patients were 4.12 ± 2.42 mJ/L and 115 ± 27/s, compared to 2.96 ± 2.16 mJ/L and 78 ± 45/s in controls. There was good agreement between RV motion and manual tracings. The difference in diastolic KE between CFD and 4D flow by Bland-Altman analysis was - 0.89910 to 2 mJ/mL (95% limits of agreement: - 1.351 × 10^-2 mJ/mL to 1.171 × 10^-2 mJ/mL).

CONCLUSION

This CFD framework can produce intracardiac flow in rTOF patients. CFD has the potential for predicting the effects of PVR in rTOF patients and improve the clinical indications guided by cMRI.

Initial experience with intensity distribution analysis of hemodynamic parameters in the thoracic aorta using four-dimensional magnetic resonance imaging: A comparison between groups with different ejection fractions

The purpose of this study was to investigate whether there were significant differences in the intensity distributions of thoracic aorta hemodynamic parameters between groups with different ejection fractions (EF) using four-dimensional flow magnetic resonance imaging and to investigate the relationships between each parameter.A total of 26 patients, 13 each with EF of >60% and <30%, underwent cardiac four-dimensional flow magnetic resonance imaging (EF >60%: mean age: 54 ± 11.6 years, EF <30%: mean age: 49.2 ± 17.2 years). The thoracic aorta was divided into the proximal and distal ascending aorta (AAo), aortic arch, and the proximal and distal descending aorta, and each section was further divided into the anterior wall, posterior wall, lesser curvature, and greater curvature. The intensity distributions of wall shear stress (WSS), energy loss (EL), and vorticity (Vort) (hemodynamic parameters) and the concordance rates between these distributions were analyzed.The concordance rate between the intensity distributions of EL and Vort was high. Only the intensity distributions of EL and Vort in the distal AAo differed significantly between the groups (P < .001). In the EF >60% group, these intensity distributions showed higher values in the greater curvature of the AAo, whereas in the EF <30% group higher values were seen in the lesser curvature of the AAo.Although there was no significant intergroup difference in the WSS intensity distribution, in the EF <30% group the WSS intensity distribution tended to exhibit higher values in the lesser curvature of the distal AAo, and the WSS intensity distribution values for the greater curvature tended to gradually increase from the arch to the proximal descending aorta.The only significant differences between the EF groups were found in the intensity distributions of EL and Vort in the distal AAo. This suggests that the distributions of atherosclerosis may be EF-dependent.

Disturbed flow's impact on cellular changes indicative of vascular aneurysm initiation, expansion, and rupture: A pathological and methodological review

Aneurysms are malformations within the arterial vasculature brought on by the structural breakdown of the microarchitecture of the vessel wall, with aneurysms posing serious health risks in the event of their rupture. Blood flow within vessels is generally laminar with high, unidirectional wall shear stressors that modulate vascular endothelial cell functionality and regulate vascular smooth muscle cells. However, altered vascular geometry induced by bifurcations, significant curvature, stenosis, or clinical interventions can alter the flow, generating low stressor disturbed flow patterns. Disturbed flow is associated with altered cellular morphology, upregulated expression of proteins modulating inflammation, decreased regulation of vascular permeability, degraded extracellular matrix, and heightened cellular apoptosis. The understanding of the effects disturbed flow has on the cellular cascades which initiate aneurysms and promote their subsequent growth can further elucidate the nature of this complex pathology. This review summarizes the current knowledge about the disturbed flow and its relation to aneurysm pathology, the methods used to investigate these relations, as well as how such knowledge has impacted clinical treatment methodologies. This information can contribute to the understanding of the development, growth, and rupture of aneurysms and help develop novel research and aneurysmal treatment techniques.

Prediction for future occurrence of type A aortic dissection using computational fluid dynamics

OBJECTIVES

The actual underlying mechanisms of acute type A aortic dissection (AAAD) are not well understood. The present study aimed to elucidate the mechanism of AAAD using computational fluid dynamics (CFD) analysis.

METHODS

We performed CFD analysis using patient-specific computed tomography imaging in 3 healthy control cases and 3 patients with AAAD. From computed tomography images, we made a healthy control model or pre-dissection model for CFD analysis. Pulsatile cardiac flow during one cardiac cycle was simulated, and a three-dimensional flow streamline was visualized to evaluate flow velocity, wall shear stress and oscillatory shear index (OSI).

RESULTS

In healthy controls, the transvalvular aortic flow was parallel to the ascending aorta. There was no spotty high OSI area at the ascending aorta. In pre-dissection patients, accelerated transvalvular aortic flow was towards the posterolateral ascending aorta. The vortex flow was observed on the side of the lesser curvature in mid-systole and expanded throughout the entire ascending aorta during diastole. Systolic wall shear stress was high due to the accelerated aortic blood flow on the side of the greater curvature of the ascending aorta. On the side of the lesser curvature, high OSI areas were observed around the vortex flow. In all pre-dissection cases, a spotty high OSI area was in close proximity to the actual primary entry site of the future AAAD.

CONCLUSIONS

The pre-onset high OSI area with vortex flow is closely associated with the future primary entry site. Therefore, we can elucidate the mechanism of AAAD with CFD analysis.

Computational evaluation of inferior vena cava filters through computational fluid dynamics methods

Numerical simulation is growing in its importance toward the design, testing and evaluation of medical devices. Computational fluid dynamics and finite element analysis allow improved calculation of stress, heat transfer, and flow to better understand the medical device environment. Current research focuses not only on improving medical devices, but also on improving the computational tools themselves. As methods and computer technology allow for faster simulation times, iterations and trials can be performed faster to collect more data. Given the adverse events associated with long-term inferior vena cava (IVC) filter placement, IVC filter design and device evaluation are of paramount importance. This work reviews computational methods used to develop, test, and improve IVC filters to ultimately serve the needs of the patient.

Estimating the Haemodynamic Streamline Vena Contracta as the Effective Orifice Area Measured from Reconstructed Multislice Phase-contrast MR Images for Patients with Moderately Accelerated Aortic Stenosis

PURPOSE

In aortic stenosis (AS), the discrepancy between moderately accelerated flow and effective orifice area (EOA) continues to pose a challenge. We developed a method of measuring the vena contracta area as hemodynamic EOA using cardiac MRI focusing on AS patients with a moderately accelerated flow to solve the problem that AS severity can currently be determined only by echocardiography.

METHODS

We investigated 40 patients with a peak transvalvular velocity > 3.0 m/s on transthoracic echocardiography (TTE). The patients were divided into highly accelerated and moderately accelerated AS groups according to whether or not the peak transvalvular velocity was ≥ 4.0 m/s. From the multislice 2D cine phase-contrast MRI data, the cross-sectional area of the vena contracta of the reconstructed streamline in the Valsalva sinus was defined as MRI-EOAs. Patient symptoms and echocardiography data, including EOA (defined as TTE-EOA), were derived from the continuity equation using TTE.

RESULTS

All participants in the highly accelerated AS group (n = 19) showed a peak velocity ≥ 4.0 m/s in MRI. Eleven patients in the moderately accelerated AS group (n = 21) had a TTE-EOA < 1.00 cm². In the moderately accelerated AS group, MRI-EOAs demonstrated a strong correlation with TTE-EOAs (r = 0.76, P < 0.01). Meanwhile, in the highly accelerated AS group, MRI-EOAs demonstrated positivity but a moderate correlation with TTE-EOAs (r = 0.63, P = 0.004). MRI-EOAs were overestimated compared to TTE-EOAs. In terms of the moderately accelerated AS group, the best cut-off value for MRI-EOAs was < 1.23 cm², compatible with TTE-EOAs < 1.00 cm², with an excellent prediction of the New York Heart Association classification ≥ III (sensitivity 87.5%, specificity 76.9%).

CONCLUSION

MRI-EOAs may be an alternative to conventional echocardiography for patients with moderately accelerated AS, especially those with discordant echocardiographic parameters.

Blood flow dynamics with four-dimensional flow cardiovascular magnetic resonance in patients with aortic stenosis before and after transcatheter aortic valve replacement

BACKGROUND

Pre- and post-procedural hemodynamic changes which could affect adverse outcomes in aortic stenosis (AS) patients who undergo transcatheter aortic valve replacement (TAVR) have not been well investigated. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) enables accurate analysis of blood flow dynamics such as flow velocity, flow pattern, wall shear stress (WSS), and energy loss (EL). We sought to examine the changes in blood flow dynamics of patients with severe AS who underwent TAVR.

METHODS

We examined 32 consecutive severe AS patients who underwent TAVR between May 2018 and June 2019 (17 men, 82 ± 5 years, median left ventricular ejection fraction 61%, 6 self-expanding valve), after excluding those without CMR because of a contraindication or inadequate imaging from the analyses. We analyzed blood flow patterns, WSS and EL in the ascending aorta (AAo), and those changes before and after TAVR using 4D flow CMR.

RESULTS

After TAVR, semi-quantified helical flow in the AAo was significantly decreased (1.4 ± 0.6 vs. 1.9 ± 0.8, P = 0.002), whereas vortical flow and eccentricity showed no significant changes. WSS along the ascending aortic circumference was significantly decreased in the left (P = 0.038) and left anterior (P = 0.033) wall at the basal level, right posterior (P = 0.011) and left (P = 0.010) wall at the middle level, and right (P = 0.012), left posterior (P = 0.019) and left anterior (P = 0.028) wall at the upper level. EL in the AAo was significantly decreased (15.6 [10.8-25.1 vs. 25.8 [18.6-36.2]] mW, P = 0.012). Furthermore, a significant negative correlation was observed between EL and effective orifice area index after TAVR (r = - 0.38, P = 0.034).

CONCLUSIONS

In severe AS patients undergoing TAVR, 4D flow CMR demonstrates that TAVR improves blood flow dynamics, especially when a larger effective orifice area index is obtained.

Blood flow energy loss: a predictor for the recovery of left ventricular function after bioprosthetic aortic valve replacement

OBJECTIVES

It is difficult to estimate the improvement in left ventricular (LV) function after aortic valve replacement (AVR). The present study aimed to evaluate whether energy loss (EL) can predict the postoperative LV function after AVR.

METHODS

Nine patients who underwent AVR with a bioprosthetic valve were enrolled in the present study. Porcine prostheses were used in 5 patients and bovine pericardial prostheses were used in 4 patients. The aortic flow pattern was visualized and EL and cardiac output (CO) were measured using 4-dimensional flow magnetic resonance imaging from the LV to the descending aorta; the EL/CO ratio in the extracted area was calculated as total EL/CO ratio.

RESULTS

With a porcine valve, a severe helical flow was observed in the ascending aorta during the holosystolic phase. In contrast, with a bovine pericardial valve, straight transvalvular aortic flow was observed in the early systolic phase and 2 large vortical flows occurred on both sides of the greater and lesser curvature of the ascending aorta after the mid-systolic period. The total EL/CO ratio was strongly correlated with LV ejection fraction improvement after AVR (r = 0.74, P = 0.02).

CONCLUSIONS

The aortic flow pattern is different between the porcine valve and bovine pericardial valve. The total EL/CO ratio is a valuable tool for evaluating the postoperative LV ejection fraction improvement after AVR. Optimization of total EL/CO ratio would have potential to improve haemodynamic performances after AVR.

Four-Dimensional Flow MRI of Abdominal Veins: A Systematic Review

The aim of this systematic review is to provide an overview of the use of Four-Dimensional Magnetic Resonance Imaging of vector blood flow (4D Flow MRI) in the abdominal veins. This study was composed according to the PRISMA guidelines 2009. The literature search was conducted in MEDLINE, Cochrane Library, EMBASE, and Web of Science. Quality assessment of the included studies was performed using the QUADAS-2 tool. The initial search yielded 781 studies and 21 studies were included. All studies successfully applied 4D Flow MRI in abdominal veins. Four-Dimensional Flow MRI was capable of discerning between healthy subjects and patients with cirrhosis and/or portal hypertension. The visual quality and inter-observer agreement of 4D Flow MRI were rated as excellent and good to excellent, respectively, and the studies utilized several different MRI data sampling strategies. By applying spiral sampling with compressed sensing to 4D Flow MRI, the blood flow of several abdominal veins could be imaged simultaneously in 18-25 s, without a significant loss of visual quality. Four-Dimensional Flow MRI might be a useful alternative to Doppler sonography for the diagnosis of cirrhosis and portal hypertension. Further clinical studies need to establish consensus regarding MRI sampling strategies in patients and healthy subjects.

Expanded Polytetrafluoroethylene Conduits With Bulging Sinuses and a Fan-Shaped Valve in Right Ventricular Outflow Tract Reconstruction

We developed a handmade expanded polytetrafluoroethylene (ePTFE) pulmonary valvular conduit (PVC) with bulging sinuses and a fan-shaped ePTFE valve for right ventricular outflow tract (RVOT) reconstruction. We aimed to investigate the results of this device in this multicenter study. From 2001 to 2020, 1776 patients underwent RVOT reconstruction using ePTFE PVCs at 65 institutions in Japan. The median age and body weight were 4.1 years (range, 3 days to 67.1 years) and 13.3 (range, 1.8-91.3) kg, respectively. The median PVC size was 18 (range, 8-24) mm. The median Z-value of the ePTFE PVC was 1.1 (range, -3.8 to 5.0). The ePTFE PVC conditions were investigated by cardiac echocardiography and catheterization. The median follow-up period was 3.3 years (range, 0 day to 16.2 years). There were only 9 cases (0.5%) with PVC-related unknown deaths. Reintervention was performed in 283 patients (15.9%), and 190 patients (10.7%) required explantation. Freedom from reintervention and explantation at 5/10 years were 86.7/61.5% and 93.0/69.1%, respectively. At the latest echocardiography, PVC regurgitation grade was better than mild in 88.4% patients. The average peak RVOT gradient was 15.7 ± 15.9 mm Hg at the latest cardiac catheterization. ePTFE PVC infection was detected in only 8 patients (0.5%). Relative stenosis due to somatic growth was the most common cause of PVC explantation. The performance of ePTFE in terms of durability, valvular performance, and the resistance against infection is considerable and may replace conventional prosthetic materials. Further improvement of the ePTFE membrane is essential to prevent valvular dysfunction.

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.

Influence of aneurysmal aortic root geometry on mechanical stress to the aortic valve leaflet

AIMS

While mechanical stress caused by blood flow, e.g. wall shear stress (WSS), and related parameters, e.g. oscillatory shear index (OSI), are increasingly being recognized as key moderators of various cardiovascular diseases, studies on valves have been limited because of a lack of appropriate imaging modalities. We investigated the influence of aortic root geometry on WSS and OSI on the aortic valve (AV) leaflet.

METHODS AND RESULTS

We applied our novel approach of intraoperative epi-aortic echocardiogram to measure the haemodynamic parameters of WSS and OSI on the AV leaflet. Thirty-six patients were included, which included those who underwent valve-sparing aortic root replacement (VSARR) with no significant aortic regurgitation (n = 17) and coronary artery bypass graft (CABG) with normal AV (n = 19). At baseline, those who underwent VSARR had a higher systolic WSS (0.52 ± 0.12 vs. 0.32 ± 0.08 Pa, respectively, P < 0.001) and a higher OSI (0.37 ± 0.06 vs. 0.29 ± 0.04, respectively, P < 0.001) on the aortic side of the AV leaflet than those who underwent CABG. Multivariate regression analysis revealed that the size of the sinus of Valsalva had a significant association with WSS and OSI. Following VSARR, WSS and OSI values decreased significantly compared with the baseline values (WSS: 0.29 ± 0.12 Pa, P < 0.001; OSI: 0.26 ± 0.09, P < 0.001), and became comparable to the values in those who underwent CABG (WSS, P = 0.42; OSI, P = 0.15).

CONCLUSIONS

Mechanical stress on the AV gets altered in correlation with the size of the aortic root. An aneurysmal aortic root may expose the leaflet to abnormal fluid dynamics. The VSARR procedure appeared to reduce these abnormalities.

Hemodynamic Modeling of Biological Aortic Valve Replacement Using Preoperative Data Only

Objectives: Prediction of aortic hemodynamics after aortic valve replacement (AVR) could help optimize treatment planning and improve outcomes. This study aims to demonstrate an approach to predict postoperative maximum velocity, maximum pressure gradient, secondary flow degree (SFD), and normalized flow displacement (NFD) in patients receiving biological AVR.

Methods: Virtual AVR was performed for 10 patients, who received actual AVR with a biological prosthesis. The virtual AVRs used only preoperative anatomical and 4D flow MRI data. Subsequently, computational fluid dynamics (CFD) simulations were performed and the abovementioned hemodynamic parameters compared between postoperative 4D flow MRI data and CFD results.

Results: For maximum velocities and pressure gradients, postoperative 4D flow MRI data and CFD results were strongly correlated (R² = 0.75 and R² = 0.81) with low root mean square error (0.21 m/s and 3.8 mmHg). SFD and NFD were moderately and weakly correlated at R² = 0.44 and R² = 0.20, respectively. Flow visualization through streamlines indicates good qualitative agreement between 4D flow MRI data and CFD results in most cases.

Conclusion: The approach presented here seems suitable to estimate postoperative maximum velocity and pressure gradient in patients receiving biological AVR, using only preoperative MRI data. The workflow can be performed in a reasonable time frame and offers a method to estimate postoperative valve prosthesis performance and to identify patients at risk of patient-prosthesis mismatch preoperatively. Novel parameters, such as SFD and NFD, appear to be more sensitive, and estimation seems harder. Further workflow optimization and validation of results seems warranted.

Cardiovascular Magnetic Resonance in Right Heart and Pulmonary Circulation Disorders

Right heart and pulmonary circulation disorders are generally caused by right ventricle (RV) pressure overload, volume overload, and cardiomyopathy, and they are associated with distinct clinical courses and therapeutic approaches, although they often may coexist. Cardiac magnetic resonance (CMR) provides a noninvasive accurate and reproducible multiplanar anatomic and functional assessment, tissue characterization, and blood flow evaluation of the right heart and pulmonary circulation. This article reviews the current status of the CMR, the most recent techniques, the new parameters and their clinical utility in diagnosis, prognosis, and therapeutic management in the right heart and pulmonary circulation disorders.

Long-term results of large-calibre expanded polytetrafluoroethylene-valved conduits with bulging sinuses

OBJECTIVES

In Japan, homograft and bovine jugular vein are available in very limited institutions for the reconstruction of the right ventricular outflow tract, and handmade expanded polytetrafluoroethylene (ePTFE)-valved conduits have been widely used instead. This study aimed to clarify the long-term outcomes and the durability of the ePTFE-valved conduits purely by narrowing down to those with large sizes to eliminate the influence of the body growth.

METHODS

Between January 2002 and December 2015, patients who underwent right ventricular outflow tract reconstruction in 34 Japanese institutions using ePTFE-valved conduits with a diameter of ≥18 mm were included. All the valved conduits were made in the authors' institution and delivered to each participating institution.

RESULTS

Overall, 502 patients were included. Early mortality was 1.4% and not related to conduit failure. The overall survival rate was 98.2% at 5 years and 96.6% at 10 years. Freedom from conduit explantation was 99.5% at 5 years and 89.0% at 10 years. Three patients (0.13 per 100 patient-years) developed infective endocarditis of the conduit, and only 1 patient required conduit removal. Pulmonary insufficiency was mild or less in 480 (96%) patients, and conduit stenosis was mild or less in 436 (88%) patients at the latest follow-up.

CONCLUSIONS

By narrowing the analyses down to only ePTFE conduits with a large size, satisfactory long-term outcomes of these conduits with a fan-shaped valve and bulging sinuses were shown. These conduits would be among the optimal choices for right ventricular outflow tract reconstruction.

From Early Morphometrics to Machine Learning-What Future for Cardiovascular Imaging of the Pulmonary Circulation?

Imaging plays a cardinal role in the diagnosis and management of diseases of the pulmonary circulation. Behind the picture itself, every digital image contains a wealth of quantitative data, which are hardly analysed in current routine clinical practice and this is now being transformed by radiomics. Mathematical analyses of these data using novel techniques, such as vascular morphometry (including vascular tortuosity and vascular volumes), blood flow imaging (including quantitative lung perfusion and computational flow dynamics), and artificial intelligence, are opening a window on the complex pathophysiology and structure-function relationships of pulmonary vascular diseases. They have the potential to make dramatic alterations to how clinicians investigate the pulmonary circulation, with the consequences of more rapid diagnosis and a reduction in the need for invasive procedures in the future. Applied to multimodality imaging, they can provide new information to improve disease characterization and increase diagnostic accuracy. These new technologies may be used as sophisticated biomarkers for risk prediction modelling of prognosis and for optimising the long-term management of pulmonary circulatory diseases. These innovative techniques will require evaluation in clinical trials and may in themselves serve as successful surrogate end points in trials in the years to come.

Computerized virtual surgery based on computational fluid dynamics simulation for planning coronary revascularization with aortic root replacement in adult congenital heart disease: a case report

A 38-year-old woman presented with exertional dyspnea and chest compression. She had undergone repair of congenital supravalvular aortic stenosis at 8 years of age. Contrast-enhanced computed tomography showed re-stenosis in the ascending aorta, bilateral coronary arterial aneurysm, and a highly thickened left ventricular wall. Release of stenosis was necessary to avoid left ventricular functional deterioration; however, it could cause demand–supply mismatch in coronary flow due to substantial left ventricular hypertrophy. Sufficient statistical evidence was not available in this situation; therefore, computerized virtual surgery based on computational fluid dynamics (CFD) was performed to predict the postoperative hemodynamics. Consequently, root replacement with in situ Carrel patch coronary reconstruction was considered a better option than coronary artery graft bypass in the left-side coronary flow supply. The patient underwent root replacement with in situ Carrel patch coronary reconstruction as planned based on CFD without any complication and was discharged 15 days postoperatively.

Computational fluid dynamics of internal mammary artery-left anterior descending artery anastomoses

OBJECTIVES

The aim of this study was to elucidate the remodelling of the internal mammary artery (IMA)-left anterior descending artery anastomosis and compare 2 different anastomosis techniques (end-to-side versus side-to-side) using computational fluid dynamics.

METHODS

This study included 9 patients. Computed tomography (CT) angiography was performed immediately after coronary artery bypass grafting (CABG) and at 3-6 months later. The computational fluid dynamics models were made using the CT data. The pulsatile 3-dimensional blood flow was achieved with the finite volume method to evaluate the postoperative morphological and haemodynamic changes at the anastomosis in each patient. Flow velocity distribution, wall shear stress (WSS) and its fluctuation oscillatory shear index were measured.

RESULTS

No early or mid-term graft occlusion was observed in the study series. In the side-to-side anastomosis, pouch formation at the distal end of IMA caused a vortex flow with low WSS immediately after CABG. However, at 3-6 months after surgery, this pouch disappeared. As a result, the laminar straight flow with uniform WSS distribution was achieved inside the anastomosis. In the end-to-side anastomosis, the anastomosis shape was remodelled, resulting in a laminar flow pattern with uniform WSS distribution. A patchy high oscillatory shear index was detected at the IMA wall on the top of anastomosis in either anastomosis techniques immediately after the surgery, but it disappeared at 3-6 months after surgery.

CONCLUSIONS

Regardless of the anastomosis technique used, a successful remodelling of the IMA-left anterior descending artery anastomosis shape was achieved a few months after surgery, resulting in a straightforward flow streamline, with uniform WSS distribution and minimal oscillatory shear index.

Computational fluid dynamics simulations of flow distribution and graft designs in apicoaortic bypass

OBJECTIVE

Apicoaortic bypass has double outlets and its graft design is similar to that of a left ventricular assist device (LVAD). The left ventricular apex to the descending aorta (LV-DsAo) bypass is widely used in apicoaortic bypass. In contrast, the left ventricular apex to the ascending aorta (LV-AsAo) bypass is standard in LVAD surgery. This study aimed to evaluate the graft designs of apicoaortic bypass and their effects on flow distribution and energy loss (EL).

METHODS

A simulation study using computational fluid dynamics was performed on the geometry and hemodynamics data obtained from a 30-year-old patient who underwent a LV-DsAo bypass. The ratio of the cardiac output (CO) through the ascending aorta (AsAo) and apicoaortic conduit was set at 50:50, 30:70, and 10:90. Regional blood flow (RBF) and EL were calculated for the different distribution ratios. As an alternative to the LV-DsAo bypass, a virtual LV-AsAo bypass surgery was performed, and each parameter was compared with that of the LV-DsAo bypass.

RESULTS

At a distribution ratio of 50:50, the RBF to the head and EL were 16.4% of the total CO and 62.0 mW in the LV-DsAo bypass, and 32.3% and 81.5 mW in the LV-AsAo bypass, respectively. The RBF to the head decreased with the CO through the AsAo in the LV-DsAo bypass, but it was constant in the LV-AsAo bypass. The EL increased inversely with the CO through the AsAo in both graft designs.

CONCLUSION

The regional blood flow distribution was different, but the trend of the EL which increased inversely with the CO through the AsAo was similar between the LV-DsAo and LV-AsAo bypasses.

Personalized Interventions: A Reality in the Next 20 Years or Pie in the Sky

There is no better representation of the need for personalization of care than the breadth and complexity of congenital heart disease. Advanced imaging modalities are now standard of care in the field, and the advancements being made to three-dimensional visualization technologies are growing as a means of pre-procedural preparation. Incorporating emerging modeling approaches, such as computational fluid dynamics, will push the limits of our ability to predict outcomes, and this information may be both obtained and utilized during a single procedure in the future. Artificial intelligence and customized devices may soon surface as realistic tools for the care of patients with congenital heart disease, as they are showing growing evidence of feasibility within other fields. This review illustrates the great strides that have been made and the persistent challenges that exist within the field of congenital interventional cardiology, a field which must continue to innovate and push the limits to achieve personalization of the interventions it provides.

A novel in vivo assessment of fluid dynamics on aortic valve leaflet using epi-aortic echocardiogram

BACKGROUND

Mechanical stress caused by blood flow, such as wall shear stress (WSS) and its related parameters, is key moderator of endothelial degeneration. However, an in vivo method to measure WSS on heart valves has not been developed.

METHODS

We developed a novel approach, based on vector flow mapping using intraoperative epi-aortic echocardiogram, to measure WSS and oscillatory shear index (OSI) on the aortic valve. We prospectively enrolled 15 patients with normal valves, who underwent coronary artery bypass graft.

RESULTS

Systolic WSS on the ventricularis (2.40 ± 0.44 Pa [1.45-3.00 Pa]) was higher than systolic WSS on the fibrosa (0.33 ± 0.08 Pa [0.14-0.47 Pa], P < .001) and diastolic WSS on the ventricularis (0.18 ± 0.07 Pa [0.04-0.28 Pa], P < .001). Oscillatory shear index on the fibrosa was higher than on the ventricularis (0.29 ± 0.04 [0.24-0.36] vs 0.05 ± 0.03 [0.01-0.12], P < .001). A pilot study involving two patients with severe aortic regurgitation showed significantly different values in fluid dynamics.

CONCLUSION

Vector flow mapping method using intraoperative epi-aortic echocardiogram is an effective way of measuring WSS and OSI on normal aortic leaflet in vivo, allowing for better understanding of the pathophysiology of aortic valve diseases.

Biomechanical Forces and Atherosclerosis: From Mechanism to Diagnosis and Treatment

The article provides an overview of current views on the role of biomechanical forces in the pathogenesis of atherosclerosis. The importance of biomechanical forces in maintaining vascular homeostasis is considered. We provide descriptions of mechanosensing and mechanotransduction. The roles of wall shear stress and circumferential wall stress in the initiation, progression and destabilization of atherosclerotic plaque are described. The data on the possibilities of assessing biomechanical factors in clinical practice and the clinical significance of this approach are presented. The article concludes with a discussion on current therapeutic approaches based on the modulation of biomechanical forces.

Fluid dynamics approach to airway obstruction

BACKGROUND

Fluid dynamics theory, which is a fundamental underlying concept applied to fluid management, has not been introduced to analyze the human respiratory system. We hypothesized that one of the potential mechanisms that promotes airflow limitation in patients with airway obstructive disease would be elucidated by using fluid dynamics theory.

METHODS

We calculated the values of pressure loss and static pressure change under virtual tracheal stenotic conditions using the fluid dynamics approach.

RESULTS

Under normal conditions, the absolute values of pressure loss and static pressure change are very low. However, once airway stenosis occurs, it is confirmed that they would be dramatically elevated.

CONCLUSIONS

The fluid dynamics approach to airway obstruction is very constructive. The treatment strategy for airway obstruction and the reasons for airflow limitation are well explained by using this approach.

Optimal 4DFlow MR sequence parameters for the assessment of internal carotid artery stenosis: a simulation study

PURPOSE

In patients with ICA stenosis, increased peak systolic velocity is a marker of stenosis at risk of ischemic stroke. 4DFlow MRI is a reproducible technique to evaluate velocities in ICA stenosis, although it seems to underestimate velocities as compared with Doppler ultrasonography. The purpose of our study was to confirm that velocities were underestimated on a new set of data acquired with a clinical 4DFlow sequence, and to devise optimal acquisition parameters for ICA stenosis exploration based on a numerical simulation.

METHODS

After review board approval, 15 healthy controls and 12 patients presenting ICA stenosis were explored with Doppler ultrasonography and 4DFlow MRI. We created a 2-dimensional simulation of ICA stenosis and its corresponding 4DFlow acquisition, and compared its mean peak systolic velocity underestimation to real MRI and Doppler. We then simulated the acquisition for voxel size ranging from 0.5 to 1.25 mm and number of phases per cardiac cycle ranging from 10 to 25.

RESULTS

On acquired data, 4DFlow MR underestimated peak systolic velocities (mean difference between Doppler and 4DFlow: - 35 cm/s), especially high velocities. With spatial and temporal resolutions equivalent to MR acquisition, our simulation yielded similar underestimation (mean difference: - 31 cm/s, P = 0.30). Simulations showed that 0.7-mm resolution and 20 phases per cardiac cycle would be necessary to record peak systolic velocities up to 250 cm/s.

CONCLUSION

Higher spatial resolution can provide accurate peak systolic velocities measurement with 4DFlow MRI, thus allowing better ICA stenosis assessment. Further studies are needed to validate the proposed parameters.

Innovative Modeling Techniques and 3D Printing in Patients with Left Ventricular Assist Devices: A Bridge from Bench to Clinical Practice

Left ventricular assist devices (LVAD) cause altered flow dynamics that may result in complications such as stroke, pump thrombosis, bleeding, or aortic regurgitation. Understanding altered flow dynamics is important in order to develop more efficient and durable pump configurations. In patients with LVAD, hemodynamic assessment is limited to imaging techniques such as echocardiography which precludes detailed assessment of fluid dynamics. In this review article, we present some innovative modeling techniques that are often used in device development or for research purposes, but have not been utilized clinically. Computational fluid dynamic (CFD) modeling is based on computer simulations and particle image velocimetry (PIV) employs ex vivo models that helps study fluid characteristics such as pressure, shear stress, and velocity. Both techniques may help elaborate our understanding of complications that occur with LVAD and could be potentially used in the future to troubleshoot LVAD-related alarms. These techniques coupled with 3D printing may also allow for patient-specific device implants, lowering the risk of complications increasing device durability.

Exploration of time sequential, patient specific 3D heart unlocks clinical understanding

Objectives

The purpose was to create a time sequential three-dimensional virtual reality model, also referred to as a four-dimensional model, to explore its possible benefit and clinical applications. We hypothesized that this novel solution allows for the visuospatial benefits of the 3D model and the dynamic benefits of other existing imaging modalities.

Background

We have seen how 3D models hold great value in medical decision making by eliminating the variable visuospatial skills of practitioners. They have proved especially invaluable concerning the correction of congenital heart defects and have altered the course of many surgeries. There are, however, limitations to three-dimensional models. The static models only show what the heart looks like in one snapshot of its cycle and do not allow for an understanding of the physiological and dynamic processes.

Methods

This solution segments a 3D heart derived from a 2D image stack, times the 18 phases of a cardiac cycle and creates a 4D model that can be manipulated in space and time through the use of virtual reality.

Results

We believe the 4D heart provides a unique understanding of in situ cardiac anatomy not possible with other imaging techniques. Our expanding case series of clinician experiences and their immediate recognition of the potency of this technique is highly encouraging and reveals the future of functional and dynamic 4D representations of anatomy.

Conclusions

The 4D heart improved our understanding around complex 3D relationships over time. We propose time and effort dedicated to 4D cardiac imaging analysis of dynamic cardiac pathologies such as hypertrophic obstructive cardiomyopathy or a pre-op Rastelli repair with a narrow outflow tract could offer tremendous insight into the medical decision-making process.

Electronic supplementary material

The online version of this article (10.1186/s41205-018-0034-7) contains supplementary material, which is available to authorized users.