The branching pattern of epicardial coronary arteries is clearly three-dimensional, with correspondingly complex flow patterns. The objective of the present study was to perform a detailed hemodynamic analysis using a three-dimensional finite element method in a left anterior descending (LAD) epicardial arterial tree, including main trunk and primary branches, based on computed tomography scans. The inlet LAD flow velocity was measured in an anesthetized pig, and the outlet pressure boundary condition was estimated based on scaling laws. The spatial and temporal wall shear stress (WSS), gradient of WSS (WSSG), and oscillatory shear index (OSI) were calculated and used to identify regions of flow disturbances in the vicinity of primary bifurcations. We found that low WSS and high OSI coincide with disturbed flows (stagnated, secondary, and reversed flows) opposite to the flow divider and lateral to the junction orifice of the main trunk and primary branches. High time-averaged WSSG occurs in regions of bifurcations, with the flow divider having maximum values. Low WSS and high OSI were found to be related through a power law relationship. Furthermore, zones of low time-averaged WSS and high OSI amplified for larger diameter ratio and high inlet flow rate. Hence, different focal atherosclerotic-prone regions may be explained by different physical mechanism associated with certain critical levels of low WSS, high OSI, and high WSSG, which are strongly affected by the diameter ratio. The implications of the flow patterns for atherogenesis are enumerated.

The analysis and visualization of flows is a central problem in visualization. Topology based methods have gained increasing interest in recent years. This article describes a method for the detection of closed streamlines in flows. It is based on a special treatment of cases where a streamline reenters a cell to prevent infinite cycling during streamline calculation. The algorithm checks for possible exits of a loop of crossed edges and detects structurally stable closed streamlines. These global features are not detected by conventional topology and feature detection algorithms.

The paper presents a topology-based visualization method for time-dependent two-dimensional vector elds. A time interpolation enables the accurate tracking of critical points and closed orbits as well as the detection and identication of structural changes. This completely characterizes the topology of the unsteady ow. Bifurcation theory provides the theoretical framework. The results are conveyed by surfaces that separate subvolumes of uniform ow behavior in a three-dimensional space-time domain.

An accurate analysis of the spatial distribution of blood flow in any organ must be based on detailed morphometry (diameters, lengths, vessel numbers, and branching pattern) of the organ vasculature. Despite the significance of detailed morphometric data, there is relative scarcity of data on 3D vascular anatomy. One of the major reasons is that the process of morphometric data collection is labor intensive. The objective of this study is to validate a novel segmentation algorithm for semi-automation of morphometric data extraction. The utility of the method is demonstrated in porcine coronary arteries imaged by computerized tomography (CT). The coronary arteries of five porcine hearts were injected with a contrast-enhancing polymer. The coronary arterial tree proximal to 1 mm was extracted from the 3D CT images. By determining the centerlines of the extracted vessels, the vessel radii and lengths were identified for various vessel segments. The extraction algorithm described in this paper is based on a topological analysis of a vector field generated by normal vectors of the extracted vessel wall. With this approach, special focus is placed on achieving the highest accuracy of the measured values. To validate the algorithm, the results were compared to optical measurements of the main trunk of the coronary arteries with microscopy. The agreement was found to be excellent with a root mean square deviation between computed vessel diameters and optical measurements of 0.16 mm (<10% of the mean value) and an average deviation of 0.08 mm. The utility and future applications of the proposed method to speed up morphometric measurements of vascular trees are discussed.

The blood flow in the myocardium has significant spatial heterogeneity. The objective of this study was to develop a biophysical model based on detailed anatomical data to determine the heterogeneity of regional myocardial flow during diastole. The model predictions were compared with experimental measurements in a diastolic porcine heart in the absence of vessel tone using nonradioactive fluorescent microsphere measurements. The results from the model and experimental measurements showed good agreement. The relative flow dispersion in the arrested, vasodilated heart was found to be 44% and 48% numerically and experimentally, respectively. Furthermore, the flow dispersion was found to have fractal characteristics with fractal dimensions (D) of 1.25 and 1.27 predicted by the model and validated by the experiments, respectively. This validated three-dimensional model of normal diastolic heart will play an important role in elucidating the spatial heterogeneity of coronary blood flow, and serve as a foundation for understanding the interplay between cardiac mechanics and coronary hemodynamics.