Search

Oxygenator thrombosis during extracorporeal membrane oxygenation (ECMO), is a complication that necessitates component replacement. ECMO centers monitor clot burden by intermittent measurement of pressure drop across the oxygenator. An increase in pressure drop at a preset flow rate suggests an increase in resistance/clot formation within the oxygenator. This monitoring method comes with inherent disadvantages such as monitoring gaps, and increased risk of air embolism and infection. We explored utilizing flow measurement, which avoids such risks, as an indicator of ECMO circuit obstructions. The hypothesis that flow rate through a shunt tube in the circuit will increase as distal resistances in the circuit increases was tested. We experimentally simulated controlled levels of oxygenator obstructions using glass microspheres in an ex vivo veno-venous ECMO circuit and measured the change in shunt flow rate using over the tube ultra-sound flow probes. A mathematical model was also used to study the effect of distal resistances in the ECMO circuit on shunt flow. Results of both the mathematical model and the experiments showed a clear and measurable increase in shunt flow with increasing levels of oxygenator obstruction. Therefore, flow monitoring appears to be an effective non-contact and continuous method to monitor for obstruction during ECMO.

The lobes of the lungs slide relative to each other during breathing. Quantifying lobar sliding can aid in better understanding lung function, better modeling of lung dynamics, and a better understanding of the limits of image registration near fissures. We propose a novel measure to characterize lobe sliding in the lung based on the displacement field obtained from image registration of CT scans. When two sliding lobes are modeled as a continuum, the discontinuity in the displacement field at the fissure will manifest as elevated maximum shear – the proposed measure – which is capable of capturing both the level and orientation of sliding. Six human lungs were analyzed using scans spanning functional residual capacity (FRC) to total lung capacity (TLC). The lung lobes were segmented and registered on a lobe-by-lobe basis to obtain the displacement field from which the proposed sliding measure was calculated. The sliding measure was found to be insignificant in the parenchyma, as relatively little tissue shear occurs here. On the other hand, it was elevated along the fissures. Thus a map of the proposed sliding measure of the entire lung clearly delineates and quantifies sliding between lung lobes. Sliding is a key aspect of lung deformation during breathing. The proposed measure may help resolve artifacts introduced by sliding in deformation analysis techniques used for radiotherapy.

Predicting pressure induced wall stress in intracranial aneurysms continues to be of interest for aneurysm safety assessment. In quasi-static analysis, there are two distinct approaches that one may take, the forward approach and the inverse approach. The inverse approach starts from a deformed configuration and thus is naturally suited to image-based, patient-specific analysis. Early studies by the authors’ team suggested that the inverse approach, in the context of estimating the wall stress in cerebral aneurysms, depends weakly on the material description. In this article, we present a population study to further demonstrate the inverse method, in particular, the remarkable feature of insensitivity to material properties. Twenty six aneurysm models derived from patient-specific images were employed in the study. Wall stresses were predicted in both the inverse and forward approaches using three material models. Results showed that, while forward computation yielded up to ~100% stress difference between some materials, the inverse solutions stayed close across materials. The inverse method, in addition to being methodologically accurate in dealing with pre-deformations, has the added convenience of insensitivity to uncertainties in wall tissue properties. New insight into the stress-geometry relation was also discussed.

Stress analyses of patient-specific vascular structures commonly assume that the reconstructed in vivo configuration is stress free although it is in a pre-deformed state. We submit that this assumption can be obviated using an inverse approach, thus increasing accuracy of stress estimates. In this paper, we introduce an inverse approach of stress analysis for cerebral aneurysms modeled as nonlinear thin shell structures, and demonstrate the method using a patient-specific aneurysm. A lesion surface derived from medical images, which corresponds to the deformed configuration under the arterial pressure, is taken as the input. The wall stress in the given deformed configuration, together with the unstressed initial configuration, are predicted by solving the equilibrium equations as opposed to traditional approach where the deformed geometry is assumed stress free. This inverse approach also possesses a unique advantage, that is, for some lesions it enables us to predict the wall stress without accurate knowledge of the wall elastic property. In this study, we also investigate the sensitivity of the wall stress to material parameters. It is found that the in-plane component of the wall stress is indeed insensitive to the material model.

Arterial aneurysms are in a pre-deformed state in vivo under non-zero pressure. The ability to determine their zero pressure geometry may help in improving accuracy of determination of stress distribution and reverse estimation of material properties from dynamic imaging data. An approximate method to recover the zero pressure geometry of the AAA is proposed. This method is motivated by the observation that the patterns in displacement field for a given AAA are strikingly consistent in an AAA under all physiological pressures. The basic principle is to leverage this observation to iteratively identify the geometry that when subjected to the in vivo pressure, will recover the geometry reconstructed from in vivo imaging. The methodology is demonstrated and validated using patient-specific AAA models.

The risk of rupture of cerebral aneurysms has been correlated with the size of the aneurysm sac. It is conceivable that geometrical shape, not just size may also be related to aneurysm rupture potential. Further, aneurysm shape may also be a factor in deciding on treatment modalities, i.e., to clip or coil. However, our ability to make use of available information on aneurysm shape remains poor. In this study, methods were developed to quantify the seemingly arbitrary three-dimensional geometry of the aneurysm sac, using differential and computational geometry techniques. From computed tomography angiography (CTA) data, the three-dimensional geometry of five unruptured human cerebral aneurysms was reconstructed. Various indices (maximum diameter, neck diameter, height, aspect ratio, bottleneck factor, bulge location, volume, surface area, Gaussian and mean curvatures, isoperimetric ratio, and convexity ratio) were utilized to characterize the geometry of these aneurysm surfaces and four size-matched hypothetical control aneurysms. The physical meanings of various indices and their possible role as prognosticators for rupture risk and presurgical planning were discussed.

Abdominal aortic aneurysm (AAA) is a local, progressive dilation of the distal aorta that risks rupture until treated. Using the law of Laplace, in vivo assessment of AAA surface geometry could identify regions of high wall tensions as well as provide critical dimensional and shape data for customized endoluminal stent grafts. In this study, six patients with AAA underwent spiral computed tomography imaging and the inner wall of each AAA was identified, digitized, and reconstructed. A biquadric surface patch technique was used to compute the local principal curvatures, which required no assumptions regarding axisymmetry or other shape characteristics of the AAA surface. The spatial distribution of AAA principal curvatures demonstrated substantial axial asymmetry, and included adjacent elliptical and hyperbolic regions. To determine how much the curvature spatial distributions were dependent on tortuosity versus bulging, the effects of AAA tortuosity were removed from the three-dimensional (3D) reconstructions by aligning the centroids of each digitized contour to the z axis. The spatial distribution of principal curvatures of the modified 3D reconstructions were found to be largely axisymmetric, suggesting that much of the surface geometric asymmetry is due to AAA bending. On average, AAA surface area increased by 56% and abdominal aortic length increased by 27% over those for the normal aorta. Our results indicate that AAA surface geometry is highly complex and cannot be simulated by simple axisymmetric models, and suggests an equally complex wall stress distribution. © 1999 Biomedical Engineering Society.

Knowledge of the biomechanical behavior of abdominal aortic aneurysm (AAA) as compared to nonaneurysmal aorta may provide information on the natural history of this disease. We have performed uniaxial tensile testing of excised human aneurysmal and nonaneurysmal abdominal aortic specimens. A new mathematical model that conforms to the fibrous structure of the vascular tissue was used to quantify the measured elastic response. We determined for each specimen the yield (sigma y) and ultimate (sigma u) strengths, the separate contribution to total tissue stiffness by elastin (EE) and collagen (EC) fibers, and a collagen recruitment parameter (A), which is a measure of the tortuosity of the collagen fibers. There was no significant difference in any of these mechanical properties between longitudinal and circumferential AAA specimens, nor in EE and EC between longitudinally oriented aneurysmal and normal specimens. A, sigma y, and sigma u were all significantly higher for the normal than for the aneurysmal group: A = 0.223 +/- 0.046 versus A = 0.091 +/- 0.009 (mean +/- SEM; p0.0005), sigma y = 121.0 +/- 32.8 N/cm2 versus sigma y = 65.2 +/- 9.5 N/cm2 (p0.05), and sigma u = 201.4 +/- 39.4 N/cm2 versus sigma u = 86.4 +/- 10.2 N/cm2 (p0.0005), respectively. Our findings suggest that the AAA tissue is isotropic with respect to these mechanical properties. The observed difference in A between aneurysmal and normal aorta may be due to the complete recruitment and loading of collagen fibers at lower extensions in the former. Our data indicate that AAA rupture may be related to a reduction in tensile strength and that the biomechanical properties of AAA should be considered in assessing the severity of an individual aneurysm.

Little attention has been given to the stresses within the wall of bioresorbable vascular prostheses and how they might affect the resorption process. We modeled the graft "complex" (inner tissue capsule, residual graft, and outer tissue capsule) as a three-layered compound tube under internal pressure. Using this biomechanical model, we studied the effects of alterations in the geometry (i.e., radius and thickness) and mechanical properties of each stratum on the overall transmural stress distribution. Hypothetical simulations were performed to investigate the possible sequence of and alterations in the radial and circumferential stresses during the resorption process. Our results suggest that early in the resorption phase, the inner tissue capsule is subjected to compressive hoop stresses and concentrated, large-magnitude compressive radial stresses. This distribution gives way to the more typical distribution for a thick-walled tube when equilibration (i.e., complete resorption) is approached. The prediction of the compressive stresses in the pseudo-intima during early resorption parallels findings of an elevated mitotic index in that region at that time. This leads to a new hypothesis, namely, that compressive stresses, both in-plane and out-of-plane with respect to the regenerated vascular cells, participate in the resorption process of bioresorbable vascular grafts by modulating elevated cellular proliferative activity and may play an important role in other aspects of vascular cell biology. Results of recent experimentation support this hypothesis.

Calculation of residual stress in arteries, using the analytical approach has been quite valuable in our understanding of its critical role in vascular mechanics. Stresses are calculated at the central section of an infinitely long tube by imposing a constant axial stretch while deforming the artery from the stress-free state to its unloaded state. However, segments used to perform opening-angle measurements have finite lengths. Further, the stress-free artery configuration is assumed to be circular. Experiments show that they are slightly noncircular. The numerical approach to residual stress calculation can allow us to study both these issues. Using 3D cylindrical geometries and an isotropic material model, we investigated how segment length can affect residual stress calculations and identified the appropriate segment length for experiments. Further, we recorded and used the true noncircular stress-free state of an artery segment, computed the residual stress distribution, and compared it to that from a similar, but circular segment. Our findings suggest that segment length must be ten times the wall thickness for it to be “long” enough. We also found that the circularity assumption may be a reasonable approximation for typical arteries.