ABSTRACT
Advancement of subject-specific in silico medicine requires new imaging protocols tailored to specific anatomical features, paired with new constitutive model development based on structure/function relationships. In this study, we develop a new dual-velocity encoding coefficient (VENC) 4D flow MRI protocol that provides unprecedented spatial and temporal resolution of in vivo aortic deformation. All previous dual-VENC 4D flow MRI studies in the literature focus on an isolated segment of the aorta, which fail to capture the full spectrum of aortic heterogeneity that exists along the vessel length. The imaging protocol developed provides high sensitivity to all blood flow velocities throughout the entire cardiac cycle, overcoming the challenge of accurately measuring the highly unsteady nonuniform flow field in the aorta. Cross-sectional area change, volumetric flow rate, and compliance are observed to decrease with distance from the heart, while pulse wave velocity (PWV) is observed to increase. A nonlinear aortic lumen pressure-area relationship is observed throughout the aorta such that a high vessel compliance occurs during diastole, and a low vessel compliance occurs during systole. This suggests that a single value of compliance may not accurately represent vessel behavior during a cardiac cycle in vivo. This high-resolution MRI data provide key information on the spatial variation in nonlinear aortic compliance, which can significantly advance the state-of-the-art of in-silico diagnostic techniques for the human aorta.
Subject(s)
Magnetic Resonance Imaging , Pulse Wave Analysis , Aorta , Blood Flow Velocity , Humans , Imaging, Three-Dimensional , Phantoms, ImagingABSTRACT
OBJECTIVE: Retrograde aortic type A dissection (RTAD) is a known complication in patients with aortic type B dissection. The purpose of this computational fluid dynamics (CFD) study was to identify haemodynamic risk factors for the occurrence of RTAD. METHODS: Computed tomographic angiography (CTA) images of 10 patients with type B dissections, who subsequently developed a RTAD, were retrospectively analysed together with patients constituting a control group (n = 10) where no further vascular events after the initial type B dissection occurred. CFD simulations were conducted based on 3D surface models of the aortic lumen derived from CTA datasets. For both groups, pressures, velocity magnitudes and wall shear stress (WSS) were compared at the site of the future RTAD entry tear and the surrounding aortic wall. RESULTS: WSS at the site of the future entry tear was significantly elevated compared with the surrounding wall (15.10 Pa vs. 5.15 Pa, p < .001) and was significantly higher in the RTAD group than in the control group (6.05 Pa, p < .002). Pressures and velocity magnitudes were not significantly elevated at the entry tear (3825.8 Pa, 0.63 m/s) compared with the aortic arch (3549.8 Pa, 0.50 m/s) or control group (3501.7 Pa, 0.62 m/s). CONCLUSIONS: Increased WSS accompanies the occurrence of RTAD. The results merit the design for a prospective study to confirm whether WSS is a risk factor for the occurrence of RTAD.
Subject(s)
Aorta/physiopathology , Aortic Aneurysm/physiopathology , Aortic Dissection/physiopathology , Hemodynamics , Models, Cardiovascular , Patient-Specific Modeling , Aortic Dissection/diagnostic imaging , Aorta/diagnostic imaging , Aortic Aneurysm/diagnostic imaging , Aortography/methods , Blood Flow Velocity , Computed Tomography Angiography , Female , Humans , Hydrodynamics , Imaging, Three-Dimensional , Male , Middle Aged , Pilot Projects , Radiographic Image Interpretation, Computer-Assisted , Regional Blood Flow , Retrospective Studies , Risk Factors , Stress, MechanicalABSTRACT
Purpose: To correlate clinical and urodynamics parameters in Multiple Sclerosis patients (MS) presenting Lower Urinary Tract Symptoms (LUTS) with both Expanded Disability Status Scale (EDSS) and changes in white matter integrity as seen on Diffusion Tensor Images (DTI). LUTS worsen throughout MS, as does lesion burden. We investigated which symptoms correlated best with structural changes in white matter structure. Materials and Methods: Ten adult women >18 years were recruited with stable MS for ≥3 months and voiding dysfunction defined as %PVR/BV > 20%. Patients participated in a clinical Urodynamic Study (UDS) and completed several questionnaires (i.e., HAM, AUASS, NBS-QoL). DTI images were acquired using a 7-Tesla Siemens MAGNETOM Terra MRI scanner. DTI maps were constructed, and individual patients were co-registered with the ICBM-DTI-81 white matter atlas to extract fractional anisotropy (FA) and mean diffusivity (MD). Pearson's correlation test was performed between each WMT and clinical parameters and between clinical parameters and the EDSS score as well. P-values < 0.05 were considered significant. Results: Of the clinical parameters, %PVR/BV obtained from the average of multiple un-instrumented uroflow assessments had significant correlations to the greatest number of WMTs. Furthermore, we observed that in all recorded clinical parameters, %PVR/BV was the only significant parameter correlated to the EDSS score. Conclusion: This study demonstrates that %PVR/BV can be used as an objective parameter to gauge WMT changes and disease progression in MS patients. Future studies are needed to refine this model.
ABSTRACT
BACKGROUND: The aim of the study was to investigate the impact of fetal-type posterior cerebral artery (fPCA) variant on morphological parameters of posterior communicating artery (PComA) aneurysms for rupture risk assessment. MATERIALS AND METHODS: A total of 98 PComA aneurysms (62 ruptured and 36 unruptured) in 98 consecutive patients were reviewed. Morphological parameters were calculated including aneurysm size, aspect ratio (AR), size ratio (SR), dome-to-neck ratio, bottleneck factor and inflow angle. Performances of morphological parameters to discriminate rupture status were compared between aneurysms with or without fPCA. RESULTS: Fetal-type posterior cerebral artery variant was determined in 39 (39.8%, 25 ruptured and 14 unruptured) lesions. The ruptured group revealed a significantly larger size (p = 0.004), AR (p = 0.003), SR (p = 0.001), and inflow angle (p < 0.001). For the aneurysms without fPCA, all morphological parameters were significantly different between ruptured and unruptured aneurysms (p < 0.05); for the aneurysms with fPCA, only inflow angle (p = 0.001) was significantly related with the rupture status. Multivariate analysis showed that SR (p = 0.035 and p = 0.011) and inflow angle (p = 0.001 and p = 0.028) were independent rupture risk factors for the total cohort and the aneurysms without fPCA; while only inflow angle (p = 0.004) revealed to be independently related with rupture status of aneurysms without fPCA. CONCLUSIONS: The performances of morphological parameters to discriminate rupture status were different between PComA aneurysms with and without fPCA variants. Inflow angle might be a reliable predictor for rupture risk of PComA aneurysms.
Subject(s)
Aneurysm, Ruptured , Intracranial Aneurysm , Humans , Intracranial Aneurysm/complications , Aneurysm, Ruptured/etiology , Aneurysm, Ruptured/pathology , Retrospective Studies , Cerebral Angiography/adverse effectsABSTRACT
OBJECTIVES: Outcome prediction in DeBakey Type III aortic dissections (ADs) remains challenging. Large variations in AD morphology, physiology and treatment exist. Here, we investigate if computational fluid dynamics (CFD) can provide an initial understanding of pressure changes in an AD computational model when covering entry and exit tears and removing the intra-arterial septum (IS). DESIGN: A computational mesh was constructed from magnetic resonance images from one patient (one entrance and one exit tear) and CFD simulations performed (scenario #1). Additional meshes were derived by virtually (1) covering the exit tear (false lumen (FL) thrombus progression) (scenario #2), (2) covering the entrance tear (thoracic endovascular treatment, TEVAR) (scenario #3) and (3) completely removing the IS (fenestration) (scenario #4). Changes in flow patterns and pressures were quantified relative to the initial mesh. RESULTS: Systolic pressures increased for #2 (300 Pa increase) with largest inter-luminal differences distally (2500 Pa). In #3, false lumen pressure decreased essentially to zero. In #4, systolic pressure in combined lumen reduced from 2400 to 800 Pa. CONCLUSIONS: CFD results from computational models of a DeBakey type III AD representing separate coverage of entrance and exit tears correlated with clinical experience. The reported results present a preliminary look at a complex clinical problem.
Subject(s)
Aortic Aneurysm/surgery , Aortic Dissection/surgery , Computer Simulation , Endovascular Procedures , Hemodynamics , Hydrodynamics , Models, Cardiovascular , Aortic Dissection/physiopathology , Aortic Aneurysm/diagnosis , Aortic Aneurysm/physiopathology , Blood Pressure , Endovascular Procedures/adverse effects , Endovascular Procedures/instrumentation , Humans , Image Interpretation, Computer-Assisted , Imaging, Three-Dimensional , Magnetic Resonance Imaging, Cine , Male , Middle Aged , Regional Blood Flow , Surgical Mesh , Treatment OutcomeABSTRACT
PURPOSE: Computational fluid dynamics (CFD) simulations are increasingly used to model cerebral aneurysm hemodynamics. We investigated the capability of phase contrast magnetic resonance imaging (pcMRI), guided by specialized software for optimal slice definition (NOVA, Vassol Inc.) as a non-invasive method to measure intra-aneurysmal blood flow patterns in-vivo. In a novel approach, these blood flow patterns measured with pcMRI were qualitatively compared to the ones calculated with CFD. MATERIALS AND METHODS: The volumetric inflow rates into three unruptured cerebral aneurysms and the temporal variations of the intra-aneurysmal blood flow patterns were recorded with pcMRI. Transient CFD simulations were performed on geometric models of these aneurysms derived from 3D digital subtraction angiograms. Calculated intra-aneurysmal blood flow patterns were compared at the times of maximum and minimum arterial inflow to the ones measured with pcMRI and the temporal variations of these patterns during the cardiac cycle were investigated. RESULTS: In all three aneurysms, the main features of intra-aneurysmal flow patterns obtained with pcMRI consisted of areas with positive velocities components and areas with negative velocities components. The measured velocities ranged from approx. +/- 60 to +/- 100 cm/sec. Comparison with calculated CFD simulations showed good correlation with regard to the spatial distribution of these areas, while differences in calculated magnitudes of velocities were found. CONCLUSION: CFD simulations using inflow boundary conditions measured with pcMRI yield main features of intra-aneurysmal velocity patterns corresponding to intra-aneurysmal measurements performed with pcMRI. Thus, pcMRI may become a valuable complementary technique to CFD simulations to obtain in-vivo reference data for the study of aneurysmal hemodynamics. More data is needed to compare and fully explore the capabilities of both methods.
Subject(s)
Computer Simulation , Hemodynamics/physiology , Intracranial Aneurysm/diagnosis , Intracranial Aneurysm/physiopathology , Magnetic Resonance Angiography/methods , Blood Flow Velocity/physiology , Carotid Artery Diseases/diagnosis , Carotid Artery Diseases/physiopathology , Carotid Artery, Internal/physiopathology , Humans , Retrospective Studies , SoftwareABSTRACT
BACKGROUND AND PURPOSE: The multiple stent placement technique has largely improved the long-term outcomes of intracranial fusiform aneurysms, but the hemodynamic mechanisms remain unclear. In this study, we analyzed the hemodynamic changes caused by different stent-placement strategies in patient-specific models using the computational fluid dynamics technique, aiming to provide evidence for clinical decision-making. MATERIALS AND METHODS: Ten vertebral artery fusiform aneurysms were included, and their patient-specific computational fluid dynamics models were reconstructed. A fast virtual stent placement technique was used to simulate sequential multiple stent placements (from a single stent to triple stents) in the vertebral artery fusiform aneurysm models. Hemodynamic parameters, including wall shear stress, pressure, oscillatory shear index, relative residence time, and flow pattern, were calculated and compared among groups with different numbers of stents. RESULTS: Virtual stents were deployed in all 10 cases successfully, consistent with the real stent configuration. Wall shear stress decreased progressively by 7.2%, 20.6%, and 25.8% as the number of stents increased. Meanwhile, relative residence time and pressure increased on average by 11.3%, 15.4%, and 45.0% and by 15.7%, 21.5%, and 28.2%. The oscillatory shear index showed no stable variation trend. Flow patterns improved by weakening the intensity of the vortices and displacing the vortex center from the aneurysmal wall. CONCLUSIONS: Stent placement modifies hemodynamic patterns in vertebral artery fusiform aneurysms, which might favor thrombosis formation in the aneurysmal sac. This effect is amplified with the number of stents deployed. However, a potential risk of rupture or recanalization exists and should be considered when planning to use the multiple stent placement technique in vertebral artery fusiform aneurysms.
Subject(s)
Hemodynamics/physiology , Hydrodynamics , Intracranial Aneurysm/surgery , Models, Biological , Vertebral Artery/surgery , Aged , Aged, 80 and over , Female , Humans , Male , Middle Aged , Stents , Stress, MechanicalABSTRACT
Using data obtained from 3D digital subtraction angiography acquisitions, computational fluid dynamics techniques were used first to assess hemodynamic factors in geometrically correct models of 3 paraclinoid aneurysms and then again for assessment after virtual removal of the aneurysms and reconstruction of the parent artery. Simulations revealed an area of relatively low and rotating wall shear stresses at the location at which each aneurysm had developed. This phenomenon, to our knowledge, has not been previously described.
Subject(s)
Angiography, Digital Subtraction , Cerebral Arteries/diagnostic imaging , Cerebral Arteries/physiology , Hemodynamics , Imaging, Three-Dimensional , Intracranial Aneurysm/diagnostic imaging , Intracranial Aneurysm/physiopathology , Cerebral Arteries/physiopathology , Humans , Time FactorsABSTRACT
Type B aortic dissections are serious diseases with a 60 to 80 % 5-year survival rate. Although typically managed with a medical treatment, surgery may be necessary in the acute/subacute or the chronic phase if significant complications are encountered. For these patients, CT angiography is the first-line imaging modality, used for indicating and preparing the surgical procedure as well as for follow-up. Physicians in charge of these patients should be familiar with the key reading points. Visceral malperfusion is the most common acute complication, while aneurysmal dilatation of the false lumen is the most common chronic complication, with surgical management generally indicated when the axial diameter of the aorta exceeds 55mm. Endovascular treatment tends to replace open surgery: it requires precise measurements and identification of the entry tear (contribution of 4D-MRA).
Subject(s)
Aortic Aneurysm, Thoracic/diagnostic imaging , Aortic Dissection/diagnostic imaging , Preoperative Period , Aged , Aged, 80 and over , Aortic Dissection/surgery , Aorta, Thoracic/surgery , Aortic Aneurysm, Thoracic/surgery , Aortography , Computed Tomography Angiography , Endovascular Procedures , Female , Humans , Magnetic Resonance Angiography , Magnetic Resonance Imaging , Male , Middle Aged , Postoperative PeriodABSTRACT
BACKGROUND AND PURPOSE: To quantify the relationship and to demonstrate redundancies between hemodynamic and structural parameters before and after virtual treatment with a flow diverter device (FDD) in cerebral aneurysms. METHODS: Steady computational fluid dynamics (CFD) simulations were performed for 10 cerebral aneurysms where FDD treatment with the SILK device was simulated by virtually reducing the porosity at the aneurysm ostium. Velocity and pressure values proximal and distal to and at the aneurysm ostium as well as inside the aneurysm were quantified. In addition, dome-to-neck ratios and size ratios were determined. Multiple correlation analysis (MCA) and hierarchical cluster analysis (HCA) were conducted to demonstrate dependencies between both structural and hemodynamic parameters. RESULTS: Velocities in the aneurysm were reduced by 0.14m/s on average and correlated significantly (p<0.05) with velocity values in the parent artery (average correlation coefficient: 0.70). Pressure changes in the aneurysm correlated significantly with pressure values in the parent artery and aneurysm (average correlation coefficient: 0.87). MCA found statistically significant correlations between velocity values and between pressure values, respectively. HCA sorted velocity parameters, pressure parameters and structural parameters into different hierarchical clusters. HCA of aneurysms based on the parameter values yielded similar results by either including all (n=22) or only non-redundant parameters (n=2, 3 and 4). CONCLUSION: Hemodynamic and structural parameters before and after virtual FDD treatment show strong inter-correlations. Redundancy of parameters was demonstrated with hierarchical cluster analysis.
Subject(s)
Computer Simulation , Hemodynamics , Intracranial Aneurysm/therapy , Humans , Hydrodynamics , Intracranial Aneurysm/pathology , Intracranial Aneurysm/physiopathology , Pressure , StentsABSTRACT
OBJECTIVE: To quantitatively compare hemodynamics simulated with steady-state and transient computational fluid dynamics (CFD) simulations in cerebral aneurysms with single inflow, with focus at the aneurysm ostium. METHODS: Transient and steady-state CFD simulations were performed in 10 cerebral aneurysms. Distributions and average values for pressure, helicity, vorticity, and velocity were qualitatively compared at proximal and distal parent artery locations, at the ostium plane, and in the aneurysm, and scaling factors between the two kinds of simulations were determined. Relative inflow and outflow areas at the ostium were compared, as were average inflow and outflow velocities. In addition, values for the pressure-loss coefficient (PLC), a recently introduced parameter to assess aneurysm rupture risk, were compared for both kinds of simulation. RESULTS: Distributions of hemodynamic parameters had a similar shape but were lower for transient than for steady-state simulations. Averaged scaling factors over cases and anatomical locations showed differences for hemodynamic parameters (0.485 ± 0.01 for pressure, 0.33 ± 0.02 for helicity, 0.58 ± 0.06 for vorticity and 0.56 ± 0.04 for velocity). Good agreement between ratios of inflow and outflow areas at the aneurysm ostium was obtained (Pearson correlation coefficient >0.97, p<0.001) and for the PLC (linear regression slope 0.73 ± 0.14, R(2)=0.75). CONCLUSIONS: Steady-state simulations are a quick alternative to transient simulation for visualizing and quantifying inflow and outflow areas at the aneurysm ostium, potentially of value when planning flow diverter treatment and for quantifying the PLC, a potential indicator of aneurysm rupture.
Subject(s)
Computer Simulation , Hemodynamics/physiology , Hydrodynamics , Intracranial Aneurysm/physiopathology , Models, Cardiovascular , HumansABSTRACT
PURPOSE: The quantification of wall motion in cerebral aneurysms is of interest for the assessment of aneurysmal rupture risk, for providing boundary conditions for computational simulations and as a validation tool for theoretical models. MATERIALS AND METHODS: 2D cine phase contrast magnetic resonance imaging (2D pcMRI) in combination with quantitative magnetic resonance angiography (QMRA) was evaluated for measuring wall motion in 7 intracranial aneurysms. In each aneurysm, 2 (in one case 3) cross sections, oriented approximately perpendicular to each other, were measured. RESULTS: The maximum aneurysmal wall distention ranged from 0.16 mm to 1.6 mm (mean 0.67 mm), the maximum aneurysmal wall contraction was -1.91 mm to -0.34 mm (mean 0.94 mm), and the average wall displacement ranged from 0.04 mm to 0.31 mm (mean 0.15 mm). Statistically significant correlations between average wall displacement and the shape of inflow curves (p-value < 0.05) were found in 7 of 15 cross sections; statistically significant correlations between the displacement of the luminal boundary center point and the shape of inflow curves (p-value < 0.05) were found in 6 of 15 cross sections. CONCLUSION: 2D pcMRI in combination with QMRA is capable of visualizing and quantifying wall motion in cerebral aneurysms. However, application of this technique is currently restricted by its limited spatial resolution.
Subject(s)
Cerebral Angiography/methods , Image Enhancement/methods , Image Processing, Computer-Assisted/methods , Imaging, Three-Dimensional/methods , Intracranial Aneurysm/diagnosis , Intracranial Aneurysm/physiopathology , Magnetic Resonance Angiography/methods , Magnetic Resonance Imaging, Cine/methods , Muscle, Smooth, Vascular/physiology , Vasoconstriction/physiology , Vasodilation/physiology , Biomechanical Phenomena , Blood Flow Velocity/physiology , Cerebral Arteries/physiopathology , Humans , Mathematical Computing , Muscle Contraction/physiology , Pulsatile Flow/physiologyABSTRACT
Pulsed high intensity focused ultrasound (pHIFU) is a method for delivering ultrasound to tissue while avoiding high temperatures. The technique has been suggested for non-destructively enhancing local uptake of drugs. Side effects include thermal necrosis; therefore, real-time monitoring of tissue temperature is advantageous. This paper outlines a method for improving the treatment efficiency of pHIFU using the MR image-guided InSightec ExAblate® 2000 system, an ultrasound system integrated into a whole body human MRI scanner with the ability to measure temperature at the treatment location in near real time. Thermal measurements obtained during treatment of a tissue phantom were used to determine appropriate heating parameters, and compared to in vivo treatment of rabbit muscle. Optimization of the treatment procedure and ultrasound transducer steering patterns was then conducted with the goal of minimizing treatment time while avoiding overheating. The optimization was performed on the basis of approximate solutions to the standard bioheat equation. The commercial system software of the Exablate® system was modified to assist in this optimization. Depending on the size of the treatment volume, the presented results demonstrate that it is possible to use the technique described to cut treatment times significantly, up to one-third of that required by the current standard treatment cycle.