Background
Properly understanding the origin and progression of the thoracic aortic aneurysm (TAA) can help prevent its growth and rupture. For a better understanding of this pathogenesis, the aortic blood flow has to be studied and interpreted in great detail. We can obtain detailed aortic blood flow information using magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) with a prescribed motion of the aortic wall.

Methods
We performed two different types of simulations—static (rigid wall) and dynamic (moving wall) for healthy control and a patient with a TAA. For the latter, we have developed a novel morphing approach based on the radial basis function (RBF) interpolation of the segmented 4D-flow MRI geometries at different time instants. Additionally, we have applied reconstructed 4D-flow MRI velocity profiles at the inlet with an automatic registration protocol.

Results
The simulated RBF-based movement of the aorta matched well with the original 4D-flow MRI geometries. The wall movement was most dominant in the ascending aorta, accompanied by the highest variation of the blood flow patterns. The resulting data indicated significant differences between the dynamic and static simulations, with a relative difference for the patient of 7.47±14.18% in time-averaged wall shear stress and 15.97±43.32% in the oscillatory shear index (for the whole domain).

Conclusions
In conclusion, the RBF-based morphing approach proved to be numerically accurate and computationally efficient in capturing complex kinematics of the aorta, as validated by 4D-flow MRI. We recommend this approach for future use in MRI-based CFD simulations in broad population studies. Performing these would bring a better understanding of the onset and growth of TAA.@en

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.

@en

Thoracic aortic aneurysm (TAA) is a local dilation of the thoracic aorta. Although universally used, aneurysm diameter alone is a poor predictor of major complications such as rupture. There is a need for better biomarkers for risk assessment that also reflect the aberrant flow patterns found in TAAs. Furthermore, hypertension is often present in TAA patients and may play a role in progression of aneurysm. The exact relation between TAAs and hypertension is poorly understood. This study aims to create a numerical model of hypertension in the aorta by using computational fluid dynamics. First, a normotensive state was simulated in which flow and resistance were kept unaltered. Second, a hypertensive state was modeled in which blood inflow was increased by 30%. Third, a hypertensive state was modeled in which the proximal and peripheral resistances and capacitance parameters from the three-element Windkessel boundary condition were adjusted to mimic an increase in resistance of the rest of the cardiovascular system. One patient with degenerative TAA and one healthy control were successfully simulated at hypertensive states and were extensively analyzed. Furthermore, three additional TAA patients and controls were simulated to validate our method. Hemodynamic variables such as wall shear stress, oscillatory shear index, endothelial cell activation potential (ECAP), vorticity and helicity were studied to gain more insight on the effects of hypertension on flow patterns in TAAs. By comparing a TAA patient and a control at normotensive state at peak-systole, helicity and vorticity were found to be lower in the TAA patient throughout the entire domain. No major changes in flow and flow derived quantities were observed for the TAA patient and control when resistance was increased. When flow rate was increased, regions with high ECAP values were found to reduce in TAA patients in the aneurysm region which could reduce the risk of thrombogenesis. Thus, it may be important to assess cardiac output in patients with TAA.

@en

OBJECTIVES: To date, it is not known if 16–20-mm extracardiac conduits are outgrown during somatic growth from childhood to adolescence. This study aims to determine total cavopulmonary connection (TCPC) haemodynamics in adolescent Fontan patients at rest and during simulated exercise and to assess the relationship between conduit size and haemodynamics. METHODS: Patient-specific, magnetic resonance imaging-based computational fluid dynamic models of the TCPC were performed in 51 extracardiac Fontan patients with 16–20-mm conduits. Power loss, pressure gradient and normalized resistance were quantified in rest and during simulated exercise. The cross-sectional area (CSA) (mean and minimum) of the vessels of the TCPC was determined and normalized for flow rate (mm²/l/min). Peak (predicted) oxygen uptake was assessed. RESULTS: The median age was 16.2 years (Q1–Q3 14.0–18.2). The normalized mean conduit CSA was 35–73% smaller compared to the inferior and superior vena cava, hepatic veins and left/right pulmonary artery (all P < 0.001). The median TCPC pressure gradient was 0.7 mmHg (Q1–Q3 0.5–0.8) and 2.0 (Q1–Q3 1.4–2.6) during rest and simulated exercise, respectively. A moderate–strong inverse nonlinear relationship was present between normalized mean conduit CSA and TCPC haemodynamics in rest and exercise. TCPC pressure gradients of >_1.0 at rest and >_3.0 mmHg during simulated exercise were observed in patients with a conduit CSA <_ 45 mm²/l/min and favourable haemodynamics (<1 mmHg during both rest and exercise) in conduits >_125 mm²/l/min. Normalized TCPC resistance correlated with (predicted) peak oxygen uptake. CONCLUSIONS: Extracardiac conduits of 16–20 mm have become relatively undersized in most adolescent Fontan patients leading to suboptimal haemodynamics.

@en

Background: Electroanatomical voltage mapping (EAVM) has been compared with late gadolinium enhancement cardiovascular magnetic resonance (LGE-CMR), which cannot delineate diffuse fibrosis. T1-mapping CMR overcomes the limitations of LGE-CMR, but it has not been directly compared against EAVM. Objectives: This study aims to assess the relationship between left ventricular (LV) endocardial voltage obtained by EAVM and extracellular volume (ECV) obtained by T1 mapping. Methods: The study investigated patients who underwent endocardial EAVM for ventricular arrhythmias (CARTO 3, Biosense Webster) together with preprocedural contrast-enhanced T1 mapping (Ingenia 3T, Philips Healthcare). After image integration, EAVM datapoints were projected onto LGE-CMR and ECV-encoded images. Average values of unipolar voltage (UV), bipolar voltage (BV), LGE transmurality, and ECV were merged from corresponding cardiac segments (6 per slice) and pooled for analysis. Results: The analysis included data from 628 segments from 18 patients (57 ± 13 years of age, 17% females, LV ejection fraction 48% ± 14%, nonischemic/ischemic cardiomyopathy/controls: 8/6/4 patients). Based on the 95th and 5th percentile values obtained from the controls, ECV >33%, BV <2.9 mV, and UV <6.7 mV were considered abnormal. There was a significant inverse association between voltage and ECV, but only in segments with abnormal ECV. Increased ECV could predict abnormal BV and UV with acceptable accuracy (area under the curve of 0.78 [95% CI: 0.74-0.83] and 0.84 [95% CI: 0.79-0.88]). Conclusions: This study found a significant inverse relationship between LV endocardial voltage and ECV. Real-time integration of T1 mapping may guide catheter mapping and may allow identification of areas of diffuse fibrosis potentially related to ventricular arrhythmias.

@en

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.

@en

There is a pressing need to establish novel biomarkers to predict the progression of thoracic aortic aneurysm (TAA) dilatation. Aside from hemodynamics, the roles of oxygen (O₂) and nitric oxide (NO) in TAA pathogenesis are potentially significant. As such, it is imperative to comprehend the relationship between aneurysm presence and species distribution in both the lumen and aortic wall. Given the limitations of existing imaging methods, we propose the use of patient-specific computational fluid dynamics (CFD) to explore this relationship. We have performed CFD simulations of O₂ and NO mass transfer in the lumen and aortic wall for two cases: a healthy control (HC) and a patient with TAA, both acquired using 4D-flow magnetic resonance imaging (MRI). The mass transfer of O₂ was based on active transport by hemoglobin, while the local variations of the wall shear stress (WSS) drove NO production. Comparing hemodynamic properties, the time-averaged WSS was considerably lower for TAA, while the oscillatory shear index and endothelial cell activation potential were notably elevated. O₂ and NO showed a non-uniform distribution within the lumen and an inverse correlation between the two species. We identified several locations of hypoxic regions for both cases due to lumen-side mass transfer limitations. In the wall, NO varied spatially, with a clear distinction between TAA and HC. In conclusion, the hemodynamics and mass transfer of NO in the aorta exhibit the potential to serve as a diagnostic biomarker for TAA. Furthermore, hypoxia may provide additional insights into the onset of other aortic pathologies.

@en

Hybrid e-health support was given to 11 insulin-dependent type 2 diabetes mellitus (DM2) patients, with electronic support plus a multi-disciplinary health support team. Challenges were low ICT and health literacy. After 50 weeks, attractiveness and feasibility of the intervention were perceived as high: recommendation 9.5 out of 10 and satisfaction 9.6 out of 10. Technology acceptance model (TAM) surveys showed high usefulness and feasibility. Acceptance and health behaviours were reinforced by the prolonged health results: aerobic and strength capacity levels were improved at 50 weeks, plus health related quality of life (plus biometric benefits and medication reductions, reported elsewhere). Regarding e-health theory, we conclude that iterative skill growth cycles are beneficial for long-term adoption and e-relationships. Next, the design analysis shows opportunities for additional affective and social support, on top of the strong benefits already apparent from the direct progress feedback loops used within the health coach processes.

@en

OBJECTIVES: Adequacy of 16-20mm extracardiac conduits for adolescent Fontan patients remains unknown. This study aims to evaluate conduit adequacy using the inferior vena cava (IVC)-conduit velocity mismatch factor along the respiratory cycle. METHODS: Real-time 2D flow MRI was prospectively acquired in 50 extracardiac (16-20mm conduits) Fontan patients (mean age 16.9 ± 4.5 years) at the subhepatic IVC, conduit and superior vena cava. Hepatic venous flow was determined by subtracting IVC flow from conduit flow. The cross-sectional area (CSA) was reported for each vessel. Mean flow and velocity was calculated during the average respiratory cycle, inspiration and expiration. The IVC-conduit velocity mismatch factor was determined as follows: Vconduit/VIVC, where V is the mean velocity. RESULTS: Median conduit CSA and IVC CSA were 221 mm2 (Q1-Q3 201-255) and 244 mm2 (Q1-Q3 203-265), respectively. From the IVC towards the conduit, flow rates increased significantly due to the entry of hepatic venous flow (IVC 1.9, Q1-Q3 1.5-2.2) versus conduit (3.3, Q1-Q3 2.5-4.0 l/min, P < 0.001). Consequently, mean velocity significantly increased (IVC 12 (Q1-Q3 11-14 cm/s) versus conduit 25 (Q1-Q3 17-31 cm/s), P < 0.001), resulting in a median IVC-conduit velocity mismatch of 1.8 (Q1-Q3 1.5-2.4), further augmenting during inspiration (median 2.3, Q1-Q3 1.8-3.0). IVC-conduit mismatch was inversely related to measured conduit size and positively correlated with conduit flow. The normalized IVC-conduit velocity mismatch factor during expiration and the entire respiratory cycle correlated with peak VO2 (r = -0.37, P = 0.014 and r = -0.31, P = 0.04, respectively). CONCLUSIONS: Important blood flow accelerations are observed from the IVC towards the conduit in adolescent Fontan patients, which is related to peak VO2. This study, therefore, raises concerns that implanted 16-20mm conduits have become undersized for older Fontan patients and future studies should clarify its effect on long-term outcome.

@en

In this study, we analyzed turbulent flows through a phantom (a 180^∘ bend with narrowing) at peak systole and a patient-specific coarctation of the aorta (CoA), with a pulsating flow, using magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). For MRI, a 4D-flow MRI is performed using a 3T scanner. For CFD, the standard k- ϵ, shear stress transport k- ω, and Reynolds stress (RSM) models are applied. A good agreement between measured and simulated velocity is obtained for the phantom, especially for CFD with RSM. The wall shear stress (WSS) shows significant differences between CFD and MRI in absolute values, due to the limited near-wall resolution of MRI. However, normalized WSS shows qualitatively very similar distributions of the local values between MRI and CFD. Finally, a direct comparison between in vivo 4D-flow MRI and CFD with the RSM turbulence model is performed in the CoA. MRI can properly identify regions with locally elevated or suppressed WSS. If the exact values of the WSS are necessary, CFD is the preferred method. For future applications, we recommend the use of the combined MRI/CFD method for analysis and evaluation of the local flow patterns and WSS in the aorta.

@en

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 (HFDcaudal conduit) and cranial (HFDcranial conduit) conduit and from the hepatic veins (HFDHV). 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; HFDHV (51.0 ± 20.6%) versus HFDcaudal conduit (48.2 ± 21.9%, p = 0.033) or HFDcranial 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.

@en

Aortic aneurysm is associated with aberrant blood flow and wall shear stress (WSS). This can be studied by coupling magnetic resonance imaging (MRI) with computational fluid dynamics (CFD). For patient-specific simulations, extra attention should be given to the variation in segmentation of the MRI data-set and its effect on WSS. We performed CFD simulations of blood flow in the aorta for ten different volunteers and provided corresponding WSS distributions. The aorta of each volunteer was segmented four times. The same inlet and outlet boundary conditions were applied for all segmentation variations of each volunteer. Steady-state CFD simulations were performed with inlet flow based on phase-contrast MRI during peak systole. We show that the commonly used comparison of mean and maximal values of WSS, based on CFD in the different segments of the thoracic aorta, yields good to excellent correlation (0.78–0.95) for rescan and moderate to excellent correlation (0.64–1.00) for intra- and interobserver reproducibility. However, the effect of geometrical variations is higher for the voxel-to-voxel comparison of WSS. With this analysis method, the correlation for different segments of the whole aorta is poor to moderate (0.43–0.66) for rescan and poor to good (0.48–0.73) for intra- and interobserver reproducibility. Therefore, we advise being critical about the CFD results based on the MRI segmentations to avoid possible misinterpretation. While the global values of WSS are similar for different modalities, the variation of results is high when considering the local distributions.

@en

Objectives Degenerative thoracic aortic aneurysm (TAA) patients are known to be at risk of life-threatening acute aortic events. Guidelines recommend preemptive surgery at diameters of greater than 55 mm, although many patients with small aneurysms show only mild growth rates and more than half of complications occur in aneurysms below this threshold. Thus, assessment of hemodynamics using 4-dimensional flow magnetic resonance has been of interest to obtain more insights in aneurysm development. Nonetheless, the role of aberrant flow patterns in TAA patients is not yet fully understood. Materials and Methods A total of 25 TAA patients and 22 controls underwent time-resolved 3-dimensional phase contrast magnetic resonance imaging with 3-directional velocity encoding (ie, 4-dimensional flow magnetic resonance imaging). Hemodynamic parameters such as vorticity, helicity, and wall shear stress (WSS) were calculated from velocity data in 3 anatomical segments of the ascending aorta (root, proximal, and distal). Regional WSS distribution was assessed for the full cardiac cycle. Results Flow vorticity and helicity were significantly lower for TAA patients in all segments. The proximal ascending aorta showed a significant increase in peak WSS in the outer curvature in TAA patients, whereas WSS values at the inner curvature were significantly lower as compared with controls. Furthermore, positive WSS gradients from sinotubular junction to midascending aorta were most prominent in the outer curvature, whereas from midascending aorta to brachiocephalic trunk, the outer curvature showed negative WSS gradients in the TAA group. Controls solely showed a positive gradient at the inner curvature for both segments. Conclusions Degenerative TAA patients show a decrease in flow vorticity and helicity, which is likely to cause perturbations in physiological flow patterns. The subsequent differing distribution of WSS might be a contributor to vessel wall remodeling and aneurysm formation.

@en