On the impact of vessel wall stiffness on quantitative flow dynamics in a synthetic model of the thoracic aorta

Toward accurate and fast velocity quantification with 3D ultrashort TE phase-contrast imaging

PURPOSE

Traditional phase-contrast MRI is affected by displacement artifacts caused by non-synchronized spatial- and velocity-encoding time points. The resulting inaccurate velocity maps can affect the accuracy of derived hemodynamic parameters. This study proposes and characterizes a 3D radial phase-contrast UTE (PC-UTE) sequence to reduce displacement artifacts. Furthermore, it investigates the displacement of a standard Cartesian flow sequence by utilizing a displacement-free synchronized-single-point-imaging MR sequence (SYNC-SPI) that requires clinically prohibitively long acquisition times.

METHODS

3D flow data was acquired at 3T at three different constant flow rates and varying spatial resolutions in a stenotic aorta phantom using the proposed PC-UTE, a Cartesian flow sequence, and a SYNC-SPI sequence as reference. Expected displacement artifacts were calculated from gradient timing waveforms and compared to displacement values measured in the in vitro flow experiments.

RESULTS

The PC-UTE sequence reduces displacement and intravoxel dephasing, leading to decreased geometric distortions and signal cancellations in magnitude images, and more spatially accurate velocity quantification compared to the Cartesian flow acquisitions; errors increase with velocity and higher spatial resolution.

CONCLUSION

PC-UTE MRI can measure velocity vector fields with greater accuracy than Cartesian acquisitions (although pulsatile fields were not studied) and shorter scan times than SYNC-SPI. As such, this approach is superior to traditional Cartesian 3D and 4D flow MRI when spatial misrepresentations cannot be tolerated, for example, when computational fluid dynamics simulations are compared to or combined with in vitro or in vivo measurements, or regional parameters such as wall shear stress are of interest.

Hemodynamic effects of entry and exit tear size in aortic dissection evaluated with in vitro magnetic resonance imaging and fluid-structure interaction simulation

Understanding the complex interplay between morphologic and hemodynamic features in aortic dissection is critical for risk stratification and for the development of individualized therapy. This work evaluates the effects of entry and exit tear size on the hemodynamics in type B aortic dissection by comparing fluid-structure interaction (FSI) simulations with in vitro 4D-flow magnetic resonance imaging (MRI). A baseline patient-specific 3D-printed model and two variants with modified tear size (smaller entry tear, smaller exit tear) were embedded into a flow- and pressure-controlled setup to perform MRI as well as 12-point catheter-based pressure measurements. The same models defined the wall and fluid domains for FSI simulations, for which boundary conditions were matched with measured data. Results showed exceptionally well matched complex flow patterns between 4D-flow MRI and FSI simulations. Compared to the baseline model, false lumen flow volume decreased with either a smaller entry tear (- 17.8 and - 18.5%, for FSI simulation and 4D-flow MRI, respectively) or smaller exit tear (- 16.0 and - 17.3%). True to false lumen pressure difference (initially 11.0 and 7.9 mmHg, for FSI simulation and catheter-based pressure measurements, respectively) increased with a smaller entry tear (28.9 and 14.6 mmHg), and became negative with a smaller exit tear (- 20.6 and - 13.2 mmHg). This work establishes quantitative and qualitative effects of entry or exit tear size on hemodynamics in aortic dissection, with particularly notable impact observed on FL pressurization. FSI simulations demonstrate acceptable qualitative and quantitative agreement with flow imaging, supporting its deployment in clinical studies.

Mechanical testing and comparison of porcine tissue, silicones and 3D-printed materials for cardiovascular phantoms

Background: Cardiovascular phantoms for patient education, pre-operative planning, surgical training, haemodynamic simulation, and device testing may help improve patient care. However, currently used materials may have different mechanical properties compared to biological tissue. Methods/Aim: The aim of this study was to investigate the mechanical properties of 3D-printing and silicone materials in comparison to biological cardiovascular tissues. Uniaxial cyclic tension testing was performed using dumbbell samples from porcine tissue (aorta, pulmonary artery, right and left ventricle). Flexible testing materials included 15 silicone (mixtures) and three 3D-printing materials. The modulus of elasticity was calculated for different deformation ranges. Results: The modulus of elasticity (0%-60%) for the aorta ranged from 0.16 to 0.18 N/mm2, for the pulmonary artery from 0.07 to 0.09 N/mm2, and for the right ventricle as well as the left ventricle short-axis from 0.1 to 0.16 N/mm2. For silicones the range of modulus of elasticity was 0.02-1.16 N/mm2, and for the 3D-printed materials from 0.85 to 1.02 N/mm2. The stress-strain curves of all tissues showed a non-linear behaviour in the cyclic tensile testing, with a distinct toe region, followed by exponential strain hardening behaviour towards the peak elongation. The vessel samples showed a more linear behaviour comparted to myocardial samples. The silicones and 3D printing materials exhibited near-linearity at higher strain ranges, with a decrease in stiffness following the initial deformation. All samples showed a deviation between the loading and unloading curves (hysteresis), and a reduction in peak force over the first few cycles (adaptation effect) at constant deformation. Conclusion: The modulus of elasticity of silicone mixtures is more in agreement to porcine cardiovascular tissues than 3D-printed materials. All synthetic materials showed an almost linear behaviour in the mechanical testing compared to the non-linear behaviour of the biological tissues, probably due to fibre recruitment mechanism in the latter.

Validating MRI-Derived Myocardial Stiffness Estimates Using In Vitro Synthetic Heart Models

Impaired cardiac filling in response to increased passive myocardial stiffness contributes to the pathophysiology of heart failure. By leveraging cardiac MRI data and ventricular pressure measurements, we can estimate in vivo passive myocardial stiffness using personalized inverse finite element models. While it is well-known that this approach is subject to uncertainties, only few studies quantify the accuracy of these stiffness estimates. This lack of validation is, at least in part, due to the absence of ground truth in vivo passive myocardial stiffness values. Here, using 3D printing, we created soft, homogenous, isotropic, hyperelastic heart phantoms of varying geometry and stiffness and simulate diastolic filling by incorporating the phantoms into an MRI-compatible left ventricular inflation system. We estimate phantom stiffness from MRI and pressure data using inverse finite element analyses based on a Neo-Hookean model. We demonstrate that our identified softest and stiffest values of 215.7 and 512.3 kPa agree well with the ground truth of 226.2 and 526.4 kPa. Overall, our estimated stiffnesses revealed a good agreement with the ground truth ([Formula: see text] error) across all models. Our results suggest that MRI-driven computational constitutive modeling can accurately estimate synthetic heart material stiffnesses in the range of 200-500 kPa.

Validation of the Reduced Unified Continuum Formulation Against In Vitro 4D-Flow MRI

We previously introduced and verified the reduced unified continuum formulation for vascular fluid-structure interaction (FSI) against Womersley's deformable wall theory. Our present work seeks to investigate its performance in a patient-specific aortic setting in which assumptions of idealized geometries and velocity profiles are invalid. Specifically, we leveraged 2D magnetic resonance imaging (MRI) and 4D-flow MRI to extract high-resolution anatomical and hemodynamic information from an in vitro flow circuit embedding a compliant 3D-printed aortic phantom. To accurately reflect experimental conditions, we numerically implemented viscoelastic external tissue support, vascular tissue prestressing, and skew boundary conditions enabling in-plane vascular motion at each inlet and outlet. Validation of our formulation is achieved through close quantitative agreement in pressures, lumen area changes, pulse wave velocity, and early systolic velocities, as well as qualitative agreement in late systolic flow structures. Our validated suite of FSI techniques offers a computationally efficient approach for numerical simulation of vascular hemodynamics. This study is among the first to validate a cardiovascular FSI formulation against an in vitro flow circuit involving a compliant vascular phantom of complex patient-specific anatomy.

Particle Image Velocimetry to Evaluate Pulse Wave Velocity in Aorta Phantom with the lnD-U Method

PURPOSE

Pulse wave velocity (PWV) is an indicator of arterial stiffness used in the prediction of cardiovascular disease such as atherosclerosis. Non-invasive methods performed with ultrasound probes allow one to compute PWV and aortic stiffness through the measurement of the aortic diameter (D) and blood flow velocity (U) with the lnD-U method. This technique based on in vivo acquisitions lacks validation since the aortic elasticity modulus cannot be verified with mechanical strength tests.

METHOD

In the present study, an alternative validation is carried out on an aorta phantom hosted in an aortic flow simulator which mimics pulsatile inflow conditions. This in vitro setup included a Particle Image Velocimetry device to visualize flow in a 2D longitudinal section of the phantom, compute velocity fields (U), and track wall displacements in the aorta phantom to measure the apparent diameter (AD) variations throughout cycles.

RESULTS

The lnD-U method was then applied to evaluate PWV (5.79 ± 0.33 m/s) and calculate the Young's modulus of the aorta phantom (0.56 ± 0.12 MPa). This last value was compared to the elasticity modulus (0.53 ± 0.07 MPa) evaluated with tensile strength tests on samples cut from the silicone phantom.

CONCLUSION

The PIV technique PWV measurement showed good agreement with the direct tensile test method with a 5.6% difference in Young's modulus. Considering the uncertainties from the two methods, the measured elasticities are consistent and close to a 50-60 years old male aortic behavior. The choice of silicone for the phantom material is a relevant and promising option to mimic the human aorta on in vitro systems.