Impact of 4D-Flow CMR Parameters on Functional Evaluation of Fontan Circulation

We sought to evaluate the potential clinical role of 4D-flow cardiac magnetic resonance (CMR)-derived energetics and flow parameters in a cohort of patients' post-Fontan palliation. In patients with Fontan circulation who underwent 4D-Flow CMR, streamlines distribution was evaluated, as well a 4D-flow CMR-derived energetics parameters as kinetic energy (KE) and energy loss (EL) normalized by volume. EL/KE index as a marker of flow efficiency was also calculated. Cardiopulmonary exercise test (CPET) was also performed in a subgroup of patients. The population study included 55 patients (mean age 22 ± 11 years). The analysis of the streamlines revealed a preferential distribution of the right superior vena cava flow for the right pulmonary artery (62.5 ± 35.4%) and a mild preferential flow for the left pulmonary artery (52.3 ± 40.6%) of the inferior vena cave-pulmonary arteries (IVC-PA) conduit. Patients with heart failure (HF) presented lower IVC/PA-conduit flow (0.75 ± 0.5 vs 1.3 ± 0.5 l/min/m2, p = 0.004) and a higher mean flow-jet angle of the IVC-PA conduit (39.2 ± 22.8 vs 15.2 ± 8.9, p < 0.001) than the remaining patients. EL/KE index correlates inversely with VO2/kg/min: R: - 0.45, p = 0.01 peak, minute ventilation (VE) R: - 0.466, p < 0.01, maximal voluntary ventilation: R:0.44, p = 0.001 and positively with the physiological dead space to the tidal volume ratio (VD/VT) peak: R: 0.58, p < 0.01. From our data, lower blood flow in IVC/PA conduit and eccentric flow was associated with HF whereas higher EL/KE index was associated with reduced functional capacity and impaired lung function. Larger studies are needed to confirm our results and to further improve the prognostic role of the 4D-Flow CMR in this challenging population.

Noninvasive Advanced Cardiovascular Magnetic Resonance-Derived Fontan Hemodynamics Are Associated With Reduced Kidney Function But Not Albuminuria

BACKGROUND

Kidney disease is the most important predictor of death in patients with a Fontan circulation, yet its clinical and hemodynamic correlates have not been well established.

METHODS AND RESULTS

A total of 53 ambulatory patients with a Fontan circulation (median age, 16.2 years, 52.8% male patients) underwent advanced cardiovascular magnetic resonance assessment, including 4-dimensional flow imaging and computational fluid dynamics. Estimated glomerular filtration rate (eGFR) <90 mL/min per 1.73 m2 was observed in 20.8% and albumin-to-creatinine ratio >3 mg/mmol in 39.6%. The average eGFR decline rate was -1.83 mL/min per 1.73 m2 per year (95% CI, -2.67 to -0.99; P<0.001). Lower eGFR was associated with older age, larger body surface area at examination, longer time since Fontan procedure, and lower systemic ventricular ejection fraction. Higher albumin-to-creatinine ratio was associated with absence of fenestration at the Fontan operation, and older age and lower systemic ventricular ejection fraction at the assessment. Lower cross-sectional area of the Fontan conduit indexed to flow (r=0.32, P=0.038), higher inferior vena cava-conduit velocity mismatch factor (r=-0.35, P=0.022), higher kinetic energy indexed to flow in the total cavopulmonary connection (r=-0.59, P=0.005), and higher total cavopulmonary connection resistance (r=-0.42, P=0.005 at rest; r=-0.43, P=0.004 during exercise) were all associated with lower eGFR but not with albuminuria.

CONCLUSIONS

Kidney dysfunction and albuminuria are common among clinically well adolescents and young adults with a Fontan circulation. Advanced cardiovascular magnetic resonance-derived metrics indicative of declining Fontan hemodynamics are associated with eGFR and might serve as targets to improve kidney health. Albuminuria might be driven by other factors that need further investigation.

Hemi-Fontan and bidirectional Glenn operations result in flow-mediated viscous energy loss at the time of stage II palliation

Background

Superior cavopulmonary connection (SCPC) for stage II palliation of hypoplastic left heart syndrome (HLHS) is achieved most frequently by either a bidirectional Glenn (BDG) or hemi-Fontan (HF) operation. The comparison of flow hemodynamic efficiency at the region of surgical reconstruction and in proximal pulmonary arteries has been evaluated primarily using computational modeling techniques with conflicting reports. The purpose of this descriptive study was to compare flow hemodynamics following stage II (BDG vs HF) using 4-dimensional flow magnetic resonance imaging (4D-Flow MRI) with particular focus on flow-mediated viscous energy loss (EL') under matched hemodynamic conditions.

Methods

Patients with hypoplastic left heart syndrome (HLHS) who underwent either HF or BDG as part of stage II palliation underwent pre-Fontan 4D-Flow MRI. Patients were matched by the pulmonary vascular resistance index, net superior vena cava (SVC) flow, right pulmonary artery (RPA) and left pulmonary artery (LPA) size, and age. Maximum EL' throughout the cardiac cycle was calculated along the SVC-RPA and SVC-LPA tracts.

Results

Eight patients who underwent HF as part of their stage II single ventricle palliation were matched with 8 patients who underwent BDG. There were no differences between the 2 groups in median volumetric indices, including end-diastolic volume (P = .278) and end-systolic volume (P = .213). Moreover, no differences were observed in ejection fraction (P = .091) and cardiac index (P = .324). There also were no differences in peak EL' measured along the SVC-RPA tract (median, 0.05 mW for HF vs 0.04 mW for BDG; P = .365) or along the SVC-LPA tract (median, 0.05 mW vs 0.04 mW; P = .741).

Conclusions

The second stage of surgical palliation of HLHS using either HF or BDG results in similar flow-mediated viscous energy loss throughout the SCPC junction. 4D-Flow MRI and computational methods should be applied together to investigate flow hemodynamic patterns throughout the Fontan palliation and overall efficiency of the Fontan circuit.

Virtual surgery to predict optimized conduit size for adult Fontan patients with 16-mm conduits

OBJECTIVES

Recent evidence suggests that conduits implanted in Fontan patients at the age of 2-4 years become undersized for adulthood. The objective of this study is to use computational fluid dynamic models to evaluate the effect of virtual expansion of the Fontan conduit on haemodynamics and energetics of the total cavopulmonary connection (TCPC) under resting conditions and increased flow conditions.

METHODS

Patient-specific, magnetic resonance imaging-based simulation models of the TCPC were performed during resting and increased flow conditions. The original 16-mm conduits were virtually enlarged to 3 new sizes. The proposed conduit sizes were defined based on magnetic resonance imaging-derived conduit flow in each patient. Flow efficiency was evaluated based on power loss, pressure drop and resistance and thrombosis risk was based on flow stagnation volume and relative residence time (RRT).

RESULTS

Models of 5 adult patients with a 16-mm extracardiac Fontan connection were simulated and subsequently virtually expanded to 24-32 mm depending on patient-specific conduit flow. Virtual expansion led to a 40-65% decrease in pressure gradient across the TCPC depending on virtual conduit size. Despite improved energetics of the entire TCPC, the pulmonary arteries remained a significant contributor to energy loss (60-73% of total loss) even after virtual surgery. Flow stagnation volume inside the virtual conduit and surface area in case of elevated RRT (>20/Pa) increased after conduit enlargement but remained negligible (flow stagnation <2% of conduit volume in rest, <0.5% with exercise and elevated RRT <3% in rest, <1% with exercise).

CONCLUSIONS

Virtual expansion of 16-mm conduits to 24-32 mm, depending on patient-specific conduit flow, in Fontan patients significantly improves TCPC efficiency while thrombosis risk presumably remains low.

Impact of pulmonary artery flow distribution on Fontan hemodynamics and flow energetics

BACKGROUND

With improved life expectancy following Fontan palliation, there is an increasing population of patients with a total cavopulmonary connection. However, there is a poor understanding of which patients will experience Fontan failure and when. 4D flow MRI has identified several metrics of clinical interest, but longitudinal studies investigating hemodynamics in Fontan patients are lacking.

OBJECTIVE

We aimed to investigate the relationship between flow distribution to the pulmonary arteries and regional hemodynamic metrics in a unique cohort with follow-up 4D flow MRI.

MATERIALS AND METHODS

Patients with > 6 months of 4D flow MRI follow-up were included. Flow distribution from the caval veins to pulmonary arteries was measured in addition to regional measures of peak velocity, viscous energy loss (EL_mean and EL_tot), and kinetic energy.

RESULTS

Ten patients with total cavopulmonary connection (17.7 ± 8.8 years at baseline, follow-up: 4.4 ± 2.6 years) were included. Five subjects had unequal flow distribution from the IVC to the pulmonary arteries at baseline. Over time, these subjects tended to have larger increases in peak velocity (39.2% vs 6.6%), EL_mean (11.6% vs -38.3%), EL_tot (9.5% vs -36.2%), and kinetic energy (96.1% vs 36.3%) in the IVC. However, these differences were statistically insignificant. We found that changes in EL_mean and EL_tot were significantly associated with changes in peak velocity in the caval veins (R² > 0.5, P < 0.001).

CONCLUSION

Unequal flow distribution from the IVC may drive increasing peak velocities and viscous energy losses, which have been associated with worse clinical outcomes. Changes in peak velocity may serve as a surrogate measure for changes in viscous energy loss.

Particle Tracing Based on 4D Flow Magnetic Resonance Imaging: A Systematic Review into Methods, Applications, and Current Developments

BACKGROUND

Particle tracing based on 4D Flow MRI has been applied as a quantitative and qualitative postprocessing technique to study temporally evolving blood flow patterns.

PURPOSE

To systematically review the various methods to perform 4D Flow MRI-based particle tracing, as well as the clinical value, clinical applications, and current developments of the technique.

STUDY TYPE

The study type is systematic review.

SUBJECTS

Patients with cardiovascular disease (such as Marfan, Fontan, Tetralogy of Fallot), healthy controls, and cardiovascular phantoms that received 4D Flow MRI with particle tracing.

FIELD STRENGTH/SEQUENCE

Three-dimensional three-directional cine phase-contrast MRI, at 1.5 T and 3 T.

ASSESSMENT

Two systematic searches were performed on the PubMed database using Boolean operators and the relevant key terms covering 4D Flow MRI and particle tracing. One systematic search was focused on particle tracing methods, whereas the other on applications. Additional articles from other sources were sought out and included after a similar inspection. Particle tracing methods, clinical applications, clinical value, and current developments were extracted.

STATISTICAL TESTS

The main results of the included studies are summarized, without additional statistical analysis.

RESULTS

Of 127 unique articles retrieved from the initial search, 56 were included (28 for methods and 54 for applications). Most articles that described particle tracing methods used an adaptive timestep, a fourth order Runge-Kutta integration method, and linear interpolation in the time dimension. Particle tracing was applied in heart chambers, aorta, venae cavae, Fontan circulation, pulmonary arteries, abdominal vasculature, peripheral arteries, carotid arteries, and cerebral vasculature. Applications were grouped as intravascular, intracardiac, flow stasis, and research.

DATA CONCLUSIONS

Particle tracing based on 4D Flow MRI gives unique insight into blood flow in several cardiovascular diseases, but the quality depends heavily on the MRI data quality. Further studies are required to evaluate the clinical value of the technique for different cardiovascular diseases.

EVIDENCE LEVEL

5.

TECHNICAL EFFICACY

Stage 1.

Performance of Cardiac MRI in Pediatric and Adult Patients with Fontan Circulation

Cardiac MRI has become a widely accepted standard for anatomic and functional assessment of complex Fontan physiology, because it is noninvasive and suitable for comprehensive follow-up evaluation after Fontan completion. The use of cardiac MRI in pediatric and adult patients after completion of the Fontan procedure are described, and a practical and experience-based cardiac MRI protocol for evaluating these patients is provided. The current approach and study protocol in use at the authors' institution are presented, which address technical considerations concerning sequences, planning, and optimal image acquisition in patients with Fontan circulation. Additionally, for each sequence, the information that can be obtained and guidance on how to integrate it into clinical decision-making is discussed. Keywords: Pediatrics, MRI, MRI Functional Imaging, Heart, Congenital © RSNA, 2022.

4D flow cardiovascular magnetic resonance derived energetics in the Fontan circulation correlate with exercise capacity and CMR-derived liver fibrosis/congestion

AIM

This study explores the relationship between in vivo 4D flow cardiovascular magnetic resonance (CMR) derived blood flow energetics in the total cavopulmonary connection (TCPC), exercise capacity and CMR-derived liver fibrosis/congestion.

BACKGROUND

The Fontan circulation, in which both caval veins are directly connected with the pulmonary arteries (i.e. the TCPC) is the palliative approach for single ventricle patients. Blood flow efficiency in the TCPC has been associated with exercise capacity and liver fibrosis using computational fluid dynamic modelling. 4D flow CMR allows for assessment of in vivo blood flow energetics, including kinetic energy (KE) and viscous energy loss rate (EL).

METHODS

Fontan patients were prospectively evaluated between 2018 and 2021 using a comprehensive cardiovascular and liver CMR protocol, including 4D flow imaging of the TCPC. Peak oxygen consumption (VO2) was determined using cardiopulmonary exercise testing (CPET). Iron-corrected whole liver T1 (cT1) mapping was performed as a marker of liver fibrosis/congestion. KE and EL in the TCPC were computed from 4D flow CMR and normalized for inflow. Furthermore, blood flow energetics were compared between standardized segments of the TCPC.

RESULTS

Sixty-two Fontan patients were included (53% male, 17.3 ± 5.1 years). Maximal effort CPET was obtained in 50 patients (peak VO2 27.1 ± 6.2 ml/kg/min, 56 ± 12% of predicted). Both KE and EL in the entire TCPC (n = 28) were significantly correlated with cT1 (r = 0.50, p = 0.006 and r = 0.39, p = 0.04, respectively), peak VO2 (r = - 0.61, p = 0.003 and r = - 0.54, p = 0.009, respectively) and % predicted peak VO2 (r = - 0.44, p = 0.04 and r = - 0.46, p = 0.03, respectively). Segmental analysis indicated that the most adverse flow energetics were found in the Fontan tunnel and left pulmonary artery.

CONCLUSIONS

Adverse 4D flow CMR derived KE and EL in the TCPC correlate with decreased exercise capacity and increased levels of liver fibrosis/congestion. 4D flow CMR is promising as a non-invasive screening tool for identification of patients with adverse TCPC flow efficiency.

The Influence of Respiration on Blood Flow in the Fontan Circulation: Insights for Imaging-Based Clinical Evaluation of the Total Cavopulmonary Connection

Congenital heart disease is the most common birth defect and functionally univentricular heart defects represent the most severe end of this spectrum. The Fontan circulation provides an unique solution for single ventricle patients, by connecting both caval veins directly to the pulmonary arteries. As a result, the pulmonary circulation in Fontan palliated patients is characterized by a passive, low-energy circulation that depends on increased systemic venous pressure to drive blood toward the lungs. The absence of a subpulmonary ventricle led to the widely believed concept that respiration, by sucking blood to the pulmonary circulation during inspiration, is of great importance as a driving force for antegrade blood flow in Fontan patients. However, recent studies show that respiration influences pulsatility, but has a limited effect on net forward flow in the Fontan circulation. Importantly, since MRI examination is recommended every 2 years in Fontan patients, clinicians should be aware that most conventional MRI flow sequences do not capture the pulsatility of the blood flow as a result of the respiration. In this review, the unique flow dynamics influenced by the cardiac and respiratory cycle at multiple locations within the Fontan circulation is discussed. The impact of (not) incorporating respiration in different MRI flow sequences on the interpretation of clinical flow parameters will be covered. Finally, the influence of incorporating respiration in advanced computational fluid dynamic modeling will be outlined.

Segmental assessment of blood flow efficiency in the total cavopulmonary connection using four-dimensional flow magnetic resonance imaging: vortical flow is associated with increased viscous energy loss rate

Aims

To study flow-related energetics in multiple anatomical segments of the total cavopulmonary connection (TCPC) in Fontan patients from four-dimensional (4D) flow magnetic resonance imaging (MRI), and to study the relationship between adverse flow patterns and segment-specific energetics.

Methods and results

Twenty-six extracardiac Fontan patients underwent 4D flow MRI of the TCPC. A segmentation of the TCPC was automatically divided into five anatomical segments [conduit, superior vena cava (SVC), right/left pulmonary artery (LPA), and the Fontan confluence]. The presence of vortical flow in the pulmonary arteries or Fontan confluence was qualitatively scored. Kinetic energy (KE), viscous energy loss rate, and vorticity were calculated from the 4D flow MRI velocity field and normalized for segment length and/or inflow. Energetics were compared between segments and the relationship between vortical flow and segment cross-sectional area (CSA) with segment-specific energetics was determined. Vortical flow in the LPA (n = 6) and Fontan confluence (n = 12) were associated with significantly higher vorticity (P = 0.001 and P = 0.015, respectively) and viscous energy loss rate (P = 0.046 and P = 0.04, respectively) compared to patients without vortical flow. The LPA and conduit segments showed the highest KE and viscous energy loss rate, while most favourable energetics were observed in the SVC. Conduit CSA inversely correlated with KE (r = −0.614, P = 0.019) and viscous energy loss rate (r = −0.652, P = 0.011).

Conclusions

Vortical flow in the Fontan confluence and LPA associated with significantly increased viscous energy loss rate. Four-dimensional flow MRI-derived energetics may be used as a screening tool for direct, MRI-based assessment of flow efficiency in the TCPC.

Non-uniform mixing of hepatic venous flow and inferior vena cava flow in the Fontan conduit

Fontan patients require a balanced hepatic blood flow distribution (HFD) to prevent pulmonary arteriovenous malformations. Currently, HFD is quantified by tracking Fontan conduit flow, assuming hepatic venous (HV) flow to be uniformly distributed within the Fontan conduit. However, this assumption may be unvalid leading to inaccuracies in HFD quantification with potential clinical impact. The aim of this study was to (i) assess the mixing of HV flow and inferior vena caval (IVC) flow within the Fontan conduit and (ii) quantify HFD by directly tracking HV flow and quantitatively comparing results with the conventional approach. Patient-specific, time-resolved computational fluid dynamic models of 15 total cavopulmonary connections were generated, including the HV and subhepatic IVC. Mixing of HV and IVC flow, on a scale between 0 (no mixing) and 1 (perfect mixing), was assessed at the caudal and cranial Fontan conduit. HFD was quantified by tracking particles from the caudal (HFD_{caudal conduit}) and cranial (HFD_{cranial conduit}) conduit and from the hepatic veins (HFD_HV). HV flow was non-uniformly distributed at both the caudal (mean mixing 0.66 ± 0.13) and cranial (mean 0.79 ± 0.11) level within the Fontan conduit. On a cohort level, differences in HFD between methods were significant but small; HFD_HV (51.0 ± 20.6%) versus HFD_{caudal conduit} (48.2 ± 21.9%, p = 0.033) or HFD_{cranial conduit} (48.0 ± 21.9%, p = 0.044). However, individual absolute differences of 8.2–14.9% in HFD were observed in 4/15 patients. HV flow is non-uniformly distributed within the Fontan conduit. Substantial individual inaccuracies in HFD quantification were observed in a subset of patients with potential clinical impact.

Reduced scan time and superior image quality with 3D flow MRI compared to 4D flow MRI for hemodynamic evaluation of the Fontan pathway

Long scan times prohibit a widespread clinical applicability of 4D flow MRI in Fontan patients. As pulsatility in the Fontan pathway is minimal during the cardiac cycle, acquiring non-ECG gated 3D flow MRI may result in a reduction of scan time while accurately obtaining time-averaged clinical parameters in comparison with 2D and 4D flow MRI. Thirty-two Fontan patients prospectively underwent 2D (reference), 3D and 4D flow MRI of the Fontan pathway. Multiple clinical parameters were assessed from time-averaged flow rates, including the right-to-left pulmonary flow distribution (main endpoint) and systemic-to-pulmonary collateral flow (SPCF). A ten-fold reduction in scan time was achieved [4D flow 15.9 min (SD 2.7 min) and 3D flow 1.6 min (SD 7.8 s), p < 0.001] with a superior signal-to-noise ratio [mean ratio of SNRs 1.7 (0.8), p < 0.001] and vessel sharpness [mean ratio 1.2 (0.4), p = 0.01] with 3D flow. Compared to 2D flow, good–excellent agreement was shown for mean flow rates (ICC 0.82–0.96) and right-to-left pulmonary flow distribution (ICC 0.97). SPCF derived from 3D flow showed good agreement with that from 4D flow (ICC 0.86). 3D flow MRI allows for obtaining time-averaged flow rates and derived clinical parameters in the Fontan pathway with good–excellent agreement with 2D and 4D flow, but with a tenfold reduction in scan time and significantly improved image quality compared to 4D flow.

Hemodynamic interplay of vorticity, viscous energy loss, and kinetic energy from 4D Flow MRI and link to cardiac function in healthy subjects and Fontan patients

The purpose of this study was to directly assess (patho)physiology of intraventricular hemodynamic interplay between four-dimensional flow cardiovascular magnetic resonance imaging (4D Flow MRI)-derived vorticity with kinetic energy (KE) and viscous energy loss (EL) over the cardiac cycle and their association to ejection fraction (EF) and stroke volume (SV). Fifteen healthy subjects and thirty Fontan patients underwent whole heart 4D Flow MRI. Ventricular vorticity, KE, and EL were computed over systole (vorticity_vol_{avg systole}, KE_{avg systole}, and EL_{avg systole}) and diastole (vorticity_vol_{avg diastole}, KE_{avg diastole}, and EL_{avg diastole}). The association between vorticity_vol and KE and EL was tested by Spearman correlation. Fontan patients were grouped to normal and impaired EF groups. A significant correlation was found between SV and vorticity in healthy subjects (systolic: ρ = 0.84, P < 0.001; diastolic: ρ = 0.81, P < 0.001) and in Fontan patients (systolic: ρ = 0.61, P < 0.001; diastolic: ρ = 0.54, P = 0.002). Healthy subjects showed positive correlation between vorticity_vol versus KE (systole: ρ = 0.96, P < 0.001; diastole: ρ = 0.90, P < 0.001) and EL (systole: ρ = 0.85, P < 0.001; diastole: ρ = 0.84, P < 0.001). Fontan patients showed significantly elevated vorticity_vol compared with healthy subjects (vorticity_vol_{avg systole}: 3.1 [2.3-3.9] vs. 1.7 [1.3-2.4] L/s, P < 0.001; vorticity_vol_{avg diastole}: 3.1 [2.0-3.7] vs. 2.1 [1.6-2.8] L/s, P = 0.002). This elevated vorticity in Fontan patients showed strong association with KE (systole: ρ = 0.91, P < 0.001; diastole: ρ = 0.85, P < 0.001) and EL (systole: ρ = 0.82, P < 0.001; diastole: ρ = 0.89, P < 0.001). Fontan patients with normal EF showed significantly higher vorticity_vol_{avg systole} and EL_{avg systole}, but significantly decreased KE _{avg diastole}, in the presence of normal SV, compared with healthy subjects. Healthy subjects show strong physiological hemodynamic interplay between vorticity with KE and EL. Fontan patients demonstrate a pathophysiological hemodynamic interplay characterized by correlation of elevated vorticity with KE and EL in the presence of maintained normal stroke volume. Altered vorticity and energetic hemodynamics are found in the presence of normal EF in Fontan patients.NEW & NOTEWORTHY Physiologic intraventricular hemodynamic interplay/coupling is present in the healthy left ventricle between vorticity versus viscous energy loss and kinetic energy from four-dimensional flow cardiovascular magnetic resonance imaging (4D Flow MRI). Conversely, Fontan patients present compensatory pathophysiologic hemodynamic coupling by an increase in intraventricular vorticity that positively correlates to viscous energy loss and kinetic energy levels in the presence of maintained normal stroke volume. Altered vorticity and energetics are found in the presence of normal ejection fraction in Fontan patients.

4D flow MRI applications in congenital heart disease

Advances in the diagnosis and management of congenital heart disease (CHD) have resulted in a growing population of patients surviving well into adulthood and requiring lifelong follow-up. Flow quantification is a central component in the assessment of patients with CHD. 4D flow magnetic resonance imaging (MRI) has emerged as a tool that enables comprehensive study of flow. It involves the acquisition of a three-dimensional time-resolved volume with velocity encoding in all three spatial directions along the cardiac cycle. This allows flow quantification and visualization of blood flow patterns as well as the study of advanced hemodynamic parameters as kinetic energy and wall shear stress. 4D flow MRI-based study of flow has given insight into the altered hemodynamics in CHD particularly in bicuspid aortic valve disease and Fontan circulation. The aim of this review is to discuss the expanding clinical and research applications of 4D flow MRI in CHD as well its limitations.Key Points• Three-dimensional velocity encoding allows not only flow quantification but also the visualization of multidirectional flow patterns and the study of advanced hemodynamic parameters.• 4D flow MRI has added insight into the abnormal hemodynamics involved in congenital heart disease in particular in bicuspid aortic valve and Fontan circulation.• The main limitation of 4D flow MRI in congenital heart disease is the relatively long scan duration required for the complete coverage of the heart and great vessels with adequate spatiotemporal resolution.