Development and feasibility study of a two-dimensional ultrasonic-measurement-integrated blood flow analysis system for hemodynamics in carotid arteries

Towards multi-modal data fusion for super-resolution and denoising of 4D-Flow MRI

4D-Flow magnetic resonance imaging (MRI) has enabled in vivo time-resolved measurement of three-dimensional blood flow velocities in the human vascular system. However, its clinical use has been hampered by two main issues, namely, low spatio-temporal resolution and acquisition noise. While patient-specific computational fluid dynamics (CFD) simulations can address the resolution and noise issues, its fidelity is impacted by accuracy of estimation of boundary conditions, model parameters, vascular geometry, and flow model assumptions. In this paper a scheme to address limitations of both modalities through data-fusion is presented. The solutions of the patient-specific CFD simulation are characterized using proper orthogonal decomposition (POD). Next, a process of projecting the 4D-Flow MRI data onto the POD basis and projection coefficient mapping using generalized dynamic mode decomposition (DMD) enables simultaneous super-resolution and denoising of 4D-Flow MRI. The method has been tested using numerical phantoms derived from patient-specific aneurysmal geometries and applied to in vivo 4D-Flow MRI data.

Photoplethysmography and ultrasonic-measurement-integrated simulation to clarify the relation between two-dimensional unsteady blood flow field and forward and backward waves in a carotid artery

Understanding the spatiotemporal change in hemodynamics is essential for the basic research of atherosclerosis. The objective of this study was to establish a methodology to clarify the relation between a two-dimensional (2D) unsteady blood flow field and forward and backward propagating waves in a carotid artery. This study utilized photoplethysmography (PPG) for blood pressure measurement and two-dimensional ultrasonic-measurement-integrated (2D-UMI) simulation for flow field analysis. The validity of the methodology was confirmed in an experiment for a carotid artery of a healthy volunteer. Synchronization between the pressure measurement and flow field analysis was achieved with an error of <10 ms. A 2D unsteady blood flow field in the carotid artery was characterized in relation to forward and backward waves. 2D-UMI simulation reproduced the flow field in which the wall shear stress takes a maximum at the time of the backward wave superiority in the systolic phase, whereas 2D ordinary simulation failed to reproduce this feature because of poor reproducibility of velocity distribution. In conclusion, the proposed methodology using PPG and 2D-UMI simulation was shown to be a potential tool to clarify the relation between 2D unsteady blood flow field and the forward and backward waves in a carotid artery.

Computational fluid dynamics simulations of blood flow regularized by 3D phase contrast MRI

Background

Phase contrast magnetic resonance imaging (PC-MRI) is used clinically for quantitative assessment of cardiovascular flow and function, as it is capable of providing directly-measured 3D velocity maps. Alternatively, vascular flow can be estimated from model-based computation fluid dynamics (CFD) calculations. CFD provides arbitrarily high resolution, but its accuracy hinges on model assumptions, while velocity fields measured with PC-MRI generally do not satisfy the equations of fluid dynamics, provide limited resolution, and suffer from partial volume effects. The purpose of this study is to develop a proof-of-concept numerical procedure for constructing a simulated flow field that is influenced by both direct PC-MRI measurements and a fluid physics model, thereby taking advantage of both the accuracy of PC-MRI and the high spatial resolution of CFD. The use of the proposed approach in regularizing 3D flow fields is evaluated.

Methods

The proposed algorithm incorporates both a Newtonian fluid physics model and a linear PC-MRI signal model. The model equations are solved numerically using a modified CFD algorithm. The numerical solution corresponds to the optimal solution of a generalized Tikhonov regularization, which provides a flow field that satisfies the flow physics equations, while being close enough to the measured PC-MRI velocity profile. The feasibility of the proposed approach is demonstrated on data from the carotid bifurcation of one healthy volunteer, and also from a pulsatile carotid flow phantom.

Results

The proposed solver produces flow fields that are in better agreement with direct PC-MRI measurements than CFD alone, and converges faster, while closely satisfying the fluid dynamics equations. For the implementation that provided the best results, the signal-to-error ratio (with respect to the PC-MRI measurements) in the phantom experiment was 6.56 dB higher than that of conventional CFD; in the in vivo experiment, it was 2.15 dB higher.

Conclusions

The proposed approach allows partial or complete measurements to be incorporated into a modified CFD solver, for improving the accuracy of the resulting flow fields estimates. This can be used for reducing scan time, increasing the spatial resolution, and/or denoising the PC-MRI measurements.

Effects of inflow velocity profile on two-dimensional hemodynamic analysis by ordinary and ultrasonic-measurement-integrated simulations

Two-dimensional ultrasonic-measurement-integrated (2D-UMI) simulation correctly reproduces hemodynamics even with an inexact inflow velocity distribution. This study aimed to investigate which is superior, a two-dimensional ordinary (2D-O) simulation with an accurate inflow velocity distribution or a 2D-UMI simulation with an inaccurate one. 2D-O and 2D-UMI simulations were performed for blood flow in a carotid artery with four upstream velocity boundary conditions: a velocity profile with backprojected measured Doppler velocities (condition A), and velocity profiles with a measured Doppler velocity distribution, a parabolic one, and a uniform one, magnitude being obtained by inflow velocity estimation (conditions B, C, and D, respectively). The error of Doppler velocity against the measurement data was sensitive to the inflow velocity distribution in the 2D-O simulation, but not in the 2D-UMI simulation with the inflow velocity estimation. Among the results in conditions B, C, and D, the error in the worst 2D-UMI simulation with condition D was 31 % of that in the best 2D-O simulation with condition B, implying the superiority of the 2D-UMI simulation with an inaccurate inflow velocity distribution over the 2D-O simulation with an exact one. Condition A resulted in a larger error than the other conditions in both the 2D-O and 2D-UMI simulations.

Study of Estimation Method for Unsteady Inflow Velocity in Two-Dimensional Ultrasonic-Measurement-Integrated Blood Flow Simulation

Information on hemodynamics is essential for elucidation of mechanisms and development of novel diagnostic methods for circulatory diseases. Two-dimensional ultrasonic-measurement-integrated (2D-UMI) simulation can correctly reproduce an intravascular blood flow field and hemodynamics by feeding back an ultrasonic measurement to the numerical blood flow simulation. In this method, it is critically important to give the correct cross-sectional average inflow velocity (inflow velocity) as the boundary condition. However, systematic study has not been done on the relative validity and effectiveness of existing inflow velocity estimation methods for various target flow fields. The aim of this study was to examine the existing methods systematically and to establish a method to accurately estimate inflow velocities for various vessel geometries and flow conditions in 2D-UMI simulations. A numerical experiment was performed for 2D-UMI simulation of blood flow models in a straight vessel with inflow velocity profiles symmetric and asymmetric to the vessel axis using existing evaluation functions based on Doppler velocity error for the inflow velocity estimation. As a result, it was clarified that a significantly large estimation error occurs in the asymmetric flow due to a nonfeedback domain near the downstream end of the calculation domain. Hence, a new inflow velocity estimation method of 2D-UMI simulation is proposed in which the feedback and evaluation domains are extended to the downstream end. Further numerical experiments of 2D-UMI simulation for two realistic vessel geometries of a healthy blood vessel and a stenosed one confirmed the effectiveness of the proposed method.