Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation

A novel multiscale mathematical and computational model of the pulmonary circulation is presented and used to analyse both arterial and venous pressure and flow. This work is a major advance over previous studies by Olufsen et al. (Ann Biomed Eng 28:1281–1299, 2012) which only considered the arterial circulation. For the first three generations of vessels within the pulmonary circulation, geometry is specified from patient-specific measurements obtained using magnetic resonance imaging (MRI). Blood flow and pressure in the larger arteries and veins are predicted using a nonlinear, cross-sectional-area-averaged system of equations for a Newtonian fluid in an elastic tube. Inflow into the main pulmonary artery is obtained from MRI measurements, while pressure entering the left atrium from the main pulmonary vein is kept constant at the normal mean value of 2 mmHg. Each terminal vessel in the network of ‘large’ arteries is connected to its corresponding terminal vein via a network of vessels representing the vascular bed of smaller arteries and veins. We develop and implement an algorithm to calculate the admittance of each vascular bed, using bifurcating structured trees and recursion. The structured-tree models take into account the geometry and material properties of the ‘smaller’ arteries and veins of radii ≥ 50 μ m. We study the effects on flow and pressure associated with three classes of pulmonary hypertension expressed via stiffening of larger and smaller vessels, and vascular rarefaction. The results of simulating these pathological conditions are in agreement with clinical observations, showing that the model has potential for assisting with diagnosis and treatment for circulatory diseases within the lung

Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

Lung Circulation Modeling: Status and Prospect

Mathematical modeling has been used to interpret anatomical and physiological data obtained from metabolic and hemodynamic studies aimed at investigating structure-function relationships in the vasculature of the lung, and how these relationships are affected by lung injury and disease. The indicator dilution method was used to study the activity of redox processes within the lung. A steady-state model of the data was constructed and used to show that pulmonary endothelial cells may play an important role in reducing redox active compounds and that those reduction rates can be altered with oxidative stress induced by exposure to high oxygen environments. In addition, a morphometric model of the pulmonary vasculature was described and used to detect, describe,and predict changes in vascular morphology that occur in response to chronic exposure to low-oxygen environments, a common model of pulmonary hypertension. Finally, the model was used to construct simulated circulatory networks designed to aid in evaluation of competing hypotheses regarding the relative contribution of various morphological and biomechanical changes observed with hypoxia. These examples illustrate the role of mathematical modeling in the integration of the emerging metabolic, hemodynamic, and morphometric databases

UV Exposed Optical Fibers with Frequency Domain Reflectometry for Device Tracking in Intra-Arterial Procedures

Shape tracking of medical devices using strain sensing properties in optical fibers has seen increased attention in recent years. In this paper, we propose a novel guidance system for intra-arterial procedures using a distributed strain sensing device based on optical frequency domain reflectometry (OFDR) to track the shape of a catheter. Tracking enhancement is provided by exposing a fiber triplet to a focused ultraviolet beam, producing high scattering properties. Contrary to typical quasi-distributed strain sensors, we propose a truly distributed strain sensing approach, which allows to reconstruct a fiber triplet in real-time. A 3D roadmap of the hepatic anatomy integrated with a 4D MR imaging sequence allows to navigate the catheter within the pre-interventional anatomy, and map the blood flow velocities in the arterial tree. We employed Riemannian anisotropic heat kernels to map the sensed data to the pre-interventional model. Experiments in synthetic phantoms and an in vivo model are presented. Results show that the tracking accuracy is suitable for interventional tracking applications, with a mean 3D shape reconstruction errors of 1.6 +/- 0.3 mm. This study demonstrates the promising potential of MR-compatible UV-exposed OFDR optical fibers for non-ionizing device guidance in intra-arterial procedures

Quantification of Pulmonary Arterial Wall Distensibility Using Parameters Extracted from Volumetric Micro-CT Images

Stiffening, or loss of distensibility, of arterial vessel walls is among the manifestations of a number of vascular diseases including pulmonary arterial hypertension. We are attempting to quantify the mechanical properties of vessel walls of the pulmonary arterial tree using parameters derived from high-resolution volumetric x-ray CT images of rat lungs. The pulmonary arterial trees of the excised lungs are filled with a contrast agent. The lungs are imaged with arterial pressures spanning the physiological range. Vessel segment diameters are measured from the inlet to the periphery, and distensibilities calculated from diameters as a function of pressure. The method shows promise as an adjunct to other morphometric techniques such as histology and corrosion casting. It possesses the advantages of being nondestructive, characterizing the vascular structures while the lungs are imaged rapidly and in a near-physiological state, and providing the ability to associate mechanical properties with vessel location in the intact tree hierarchy

Aging is associated with an earlier arrival of reflected waves without a distal shift in reflection sites

Background-Despite pronounced increases in central pulse wave velocity (PWV) with aging, reflected wave transit time (RWTT), traditionally defined as the timing of the inflection point (T-INF) in the central pressure waveform, does not appreciably decrease, leading to the controversial proposition of a "distal-shift" of reflection sites. T-INF, however, is exceptionally prone to measurement error and is also affected by ejection pattern and not only by wave reflection. We assessed whether RWTT, assessed by advanced pressure-flow analysis, demonstrates the expected decline with aging. Methods and Results-We studied a sample of unselected adults without cardiovascular disease (n=48; median age 48 years) and a clinical population of older adults with suspected/established cardiovascular disease (n=164; 61 years). We measured central pressure and flow with carotid tonometry and phase-contrast MRI, respectively. We assessed RWTT using wave-separation analysis (RWTTWSA) and partially distributed tube-load (TL) modeling (RWTTTL). Consistent with previous reports, T-INF did not appreciably decrease with age despite pronounced increases in PWV in both populations. However, aging was associated with pronounced decreases in RWTTWSA (general population -15.0 ms/decade, P<0.001; clinical population -9.07 ms/decade, P=0.003) and RWTTTL (general -15.8 ms/decade, P<0.001; clinical -11.8 ms/decade, P<0.001). There was no evidence of an increased effective reflecting distance by either method. TINF was shown to reliably represent RWTT only under highly unrealistic assumptions about input impedance. Conclusions-RWTT declines with age in parallel with increased PWV, with earlier effects of wave reflections and without a distal shift in reflecting sites. These findings have important implications for our understanding of the role of wave reflections with aging

Study of cardiovascular function using a coupled left ventricle and systemic circulation model

To gain insight into cardio-arterial interactions, a coupled left ventricle-systemic artery (LV–SA) model is developed that incorporates a three-dimensional finite-strain left ventricle (LV), and a physiologically-based one-dimensional model for the systemic arteries (SA). The coupling of the LV model and the SA model is achieved by matching the pressure and the flow rate at the aortic root, i.e. the SA model feeds back the pressure as a boundary condition to the LV model, and the aortic flow rate from the LV model is used as the input for the SA model. The governing equations of the coupled system are solved using a combined immersed-boundary finite-element (IB/FE) method and a Lax–Wendroff scheme. A baseline case using physiological measurements of healthy subjects, and four exemplar cases based on different physiological and pathological scenarios are studied using the LV–SA model. The results of the baseline case agree well with published experimental data. The four exemplar cases predict varied pathological responses of the cardiovascular system, which are also supported by clinical observations. The new model can be used to gain insight into cardio-arterial interactions across a range of clinical applications

Visualization and Correction of Automated Segmentation, Tracking and Lineaging from 5-D Stem Cell Image Sequences

Results: We present an application that enables the quantitative analysis of multichannel 5-D (x, y, z, t, channel) and large montage confocal fluorescence microscopy images. The image sequences show stem cells together with blood vessels, enabling quantification of the dynamic behaviors of stem cells in relation to their vascular niche, with applications in developmental and cancer biology. Our application automatically segments, tracks, and lineages the image sequence data and then allows the user to view and edit the results of automated algorithms in a stereoscopic 3-D window while simultaneously viewing the stem cell lineage tree in a 2-D window. Using the GPU to store and render the image sequence data enables a hybrid computational approach. An inference-based approach utilizing user-provided edits to automatically correct related mistakes executes interactively on the system CPU while the GPU handles 3-D visualization tasks. Conclusions: By exploiting commodity computer gaming hardware, we have developed an application that can be run in the laboratory to facilitate rapid iteration through biological experiments. There is a pressing need for visualization and analysis tools for 5-D live cell image data. We combine accurate unsupervised processes with an intuitive visualization of the results. Our validation interface allows for each data set to be corrected to 100% accuracy, ensuring that downstream data analysis is accurate and verifiable. Our tool is the first to combine all of these aspects, leveraging the synergies obtained by utilizing validation information from stereo visualization to improve the low level image processing tasks.Comment: BioVis 2014 conferenc

Automation Process for Morphometric Analysis of Volumetric CT Data from Pulmonary Vasculature in Rats

With advances in medical imaging scanners, it has become commonplace to generate large multidimensional datasets. These datasets require tools for a rapid, thorough analysis. To address this need, we have developed an automated algorithm for morphometric analysis incorporating A Visualization Workshop computational and image processing libraries for three-dimensional segmentation, vascular tree generation and structural hierarchical ordering with a two-stage numeric optimization procedure for estimating vessel diameters. We combine this new technique with our mathematical models of pulmonary vascular morphology to quantify structural and functional attributes of lung arterial trees. Our physiological studies require repeated measurements of vascular structure to determine differences in vessel biomechanical properties between animal models of pulmonary disease. Automation provides many advantages including significantly improved speed and minimized operator interaction and biasing. The results are validated by comparison with previously published rat pulmonary arterial micro-CT data analysis techniques, in which vessels were manually mapped and measured using intense operator intervention