Computational fluid dynamics (CFD) are the gold standard in studying blood flow dynamics. However, CFD results are dependent on the boundary conditions and the computation model. The purpose of this study was to validate CFD methods using comparison with actual measurements of the blood flow vector obtained with four-dimensional (4D) flow magnetic resonance imaging (MRI). 4D Flow MRI was performed on a healthy adult and a child with double-aortic arch. The aortic lumen was segmented to visualize the blood flow. The CFD analyses were performed for the same geometries based on three turbulent models: laminar, large eddy simulation (LES), and the renormalization group k–ε model (RNG k–ε). The flow-velocity vector components, namely the wall shear stress (WSS) and flow energy loss (EL), of the MRI and CFD results were compared. The flow rate of the MRI results was underestimated in small vessels, including the neck vessels. Spiral flow in the ascending aorta caused by the left ventricular twist was observed by MRI. Secondary flow distal to the aortic arch was well realized in both CFD and MRI. The average correlation coefficients of the velocity vector components of MRI and CFD for the child were the highest for the RNG k–ε model (0.530 in ascending aorta, 0.768 in the aortic arch, 0.584 in the descending aorta). The WSS and EL values of MRI were less than half of those of CFD, but the WSS distribution patterns were quite similar. The WSS and EL estimates were higher in RNG k–ε and LES than in the laminar model because of eddy viscosity. The CFD computation realized accurate flow distal to the aortic arch, and the WSS distribution was well simulated compared to actual measurement using 4D Flow MRI. However, the helical flow was not simulated in the ascending aorta. The accuracy was enhanced by using the turbulence model, and the RNG k–ε model showed the highest correlation with 4D Flow MRI.

Blood flow imaging is a novel technology in cardiovascular medicine and surgery. Today, two types of blood flow imaging tools are available: measurement-based flow visualization including 4D flow MRI (or 3D cine phase-contrast magnetic resonance imaging), or echocardiography flow visualization software, and computer flow simulation modeling based on computational fluid dynamics (CFD). MRI and echocardiography flow visualization provide measured blood flow but have limitations in temporal and spatial resolution, whereas CFD flow calculates the flow according to assumptions instead of flow measurement, and it has sufficiently fine resolution up to the computer memory limit, and it enables even virtual surgery when combined with computer graphics. Blood flow imaging provides profound insight into the pathophysiology of cardiovascular diseases, because it quantifies and visualizes mechanical stress on the vessel walls or heart ventricle. Wall shear stress (WSS) is a stress on the endothelial wall caused by the near wall blood flow, and it is thought to be a predictor of atherosclerosis progression in coronary or aortic diseases. Flow energy loss (EL) is the loss of blood flow energy caused by viscous friction of turbulent diseased flow, and it is expected to be a predictor of ventricular workload on various heart diseases including heart valve disease, cardiomyopathy, and congenital heart diseases. Blood flow imaging can provide useful information for developing predictive medicine in cardiovascular diseases, and may lead to breakthroughs in cardiovascular surgery, especially in the decision-making process.

The purpose of this study was to investigate the accuracy of a quantitative method for estimating arterial hepatic blood flow and portal hepatic blood flow separately using a dual-input single-compartment model compared with the maximum slope method using computer simulations and clinical data. In computer simulations, the rate constants for the transfer of contrast agent (CA) from the hepatic artery to the liver (K(1a)), from the portal vein to the liver (K(1p)) and from the liver to the blood (k(2)) were estimated from simulated time-density curves with various transit times of CA from the aorta to the liver (tau(a)) and from the portal vein to the liver (tau(p)) using the linear least-squares (LLSQ) method. In clinical studies, dynamic CT data were acquired from 27 patients, and parametric maps of K(1a), K(1p) and k(2) were generated by applying the LLSQ method pixel by pixel. In simulation studies, tau(a) and tau(p) were found to have a large and a small effect on the estimates of K(1a) and K(1p), respectively. In clinical studies, the K(1a) and K(1p) values estimated with the maximum slope method were underestimated by 60+/-29% and 37+/-12%, respectively, compared with those estimated by the LLSQ method. In conclusion, our results suggest that correction of tau(a) is necessary for accurately estimating K(1a) and K(1p). Our method is therefore promising for the evaluation of hepatic blood flow in various liver diseases because it allows us to evaluate arterial hepatic blood flow and portal hepatic blood flow separately and visually.

Mitral valve morphology after mitral valve surgery affects postoperative intraventricular flow patterns and long-term cardiac performance. We visualized ventricular flow by echocardiography vector flow mapping (VFM) to reveal the impact of different mitral valve procedures. Eleven cases of mechanical mitral valve replacement (nine in the anti-anatomical and two in the anatomical position), three bioprosthetic mitral valve replacements, and four mitral valve repairs were evaluated. The mean age at the procedure was 57.4 ± 17.8 year, and the echocardiography VFM in the apical long-axis view was performed 119.9 ± 126.7 months later. Flow energy loss (EL), kinetic pressure (KP), and the flow energy efficiency ratio (EL/KP) were measured. The cases with MVR in the anatomical position and with valve repair had normal vortex directionality ("Clockwise"; N = 6), whereas those with MVR in the anti-anatomical position and with a bioprosthetic mitral valve had the vortex in the opposite direction ("Counterclockwise"; N = 12). During diastole, vortex direction had no effect on EL ("Clockwise": 0.080 ± 0.025 W/m; "Counterclockwise": 0.083 ± 0.048 W/m; P = 0.31) or KP ("Clockwise": 0.117 ± 0.021 N; "Counterclockwise": 0.099 ± 0.057 N; P = 0.023). However, during systole, the EL/KP ratio was significantly higher in the "Counterclockwise" vortex than that in the "Clockwise" vortex (1.056 ± 0.463 vs. 0.617 ± 0.158; P = 0.009). MVP and MVR with a mechanical valve in the anatomical position preserve the physiological vortex, whereas MVR with a mechanical valve in the anti-anatomical position and a bioprosthetic mitral valve generate inefficient vortex flow patterns, resulting in a potential increase in excessive cardiac workload.

Quantitative measurement of hepatic perfusion has the potential to provide important information in the assessment and management of various liver diseases. The utility of hepatic perfusion characterization relies on the resolution of each component of its dual blood supply, i.e. the hepatic artery and portal vein. In this study, a linear equation was derived by integrating the differential equation describing the kinetic behaviour of contrast agent (CA) in a dual-input single-compartment model, from which the kinetic parameters can be easily obtained using the linear least-squares method. The usefulness of this method was investigated using computer simulations, in comparison with the non-linear least-squares (NLSQ) method. This method calculated the kinetic parameters faster than the NLSQ method by a factor of approximately 10, with almost the same accuracy as the NLSQ method. This method will be useful for analysing the kinetic behaviour of CA in the unique liver environment, especially by generating the functional images of kinetic parameters.