Gilwoo Choi

Circumferential and longitudinal cyclic strain of the human thoracic aorta: Age-related changes

OBJECTIVE We developed a novel method using anatomic markers along the thoracic aorta to accurately quantify longitudinal and circumferential cyclic strain in nondiseased thoracic aortas during the cardiac cycle and to compute age-related changes of the human thoracic aorta. METHODS Changes in thoracic aorta cyclic strains were quantified using cardiac-gated computed tomography image data of 14 patients (aged 35 to 80 years) with no visible aortic pathology (aneurysms or dissection). We measured the diameter and circumferential cyclic strain in the arch and descending thoracic aorta (DTA), the longitudinal cyclic strain along the DTA, and changes in arch length and motion of the ascending aorta relative to the DTA. Diameters were computed distal to the left coronary artery, proximal and distal to the brachiocephalic trunk, and distal to the left common carotid, left subclavian, and the first and seventh intercostal arteries. Cyclic strains were computed using the Green-Lagrange strain tensor. Arch length was defined along the vessel centerline from the left coronary artery to the first intercostal artery. The length of the DTA was defined along the vessel centerline from the first to seventh intercostal artery. Longitudinal cyclic strain was quantified as the difference between the systolic and diastolic DTA lengths divided by the diastolic DTA length. Comparisons were made between seven younger (age, 41 +/- 7 years; 5 men) and seven older (age, 68 +/- 6 years; 5 men) patients. RESULTS The average increase of diameters of the thoracic aorta was 14% with age from the younger to the older (mean age, 41 vs 68 years) group. The average circumferential cyclic strain of the thoracic aorta decreased by 55% with age from the younger to the older group. The longitudinal cyclic strain decreased with age by 50% from the younger to older group (2.0% +/- 0.4% vs 1.0% +/- 1%, P = .03). The arch length increased by 14% with age from the younger to the older group (134 +/- 17 mm vs 152 +/- 10 mm, P = .03). CONCLUSIONS The thoracic aorta enlarges circumferentially and axially and deforms significantly less in the circumferential and longitudinal directions with increasing age. To our knowledge, this is the first quantitative description of in vivo longitudinal cyclic strain and length changes for the human thoracic aorta, creating a foundation for standards in reporting data related to in vivo deformation and may have significant implications in endoaortic device design, testing, and stability.

The Effect of Aging on Deformations of the Superficial Femoral Artery Resulting from Hip and Knee Flexion: Potential Clinical Implications

PURPOSE Vessel deformations have been implicated in endoluminal device fractures, and therefore better understanding of these deformations could be valuable for device regulation, evaluation, and design. The purpose of this study is to describe geometric changes of the superficial femoral artery (SFA) resulting from hip and knee flexion in older subjects. MATERIALS AND METHODS The SFAs of seven healthy subjects aged 50-70 years were imaged with magnetic resonance angiography with the legs straight and with hip and knee flexion. From geometric models constructed from these images, axial, twisting, and bending deformations were quantified. RESULTS There was greater shortening in the bottom third of the SFA than in the top two thirds (top, 5.9% +/- 3.0%; middle, 6.7% +/- 2.1%; bottom, 8.1% +/- 2.0% [mean +/- SD]; P < .05), significant twist in all sections (top, 1.3 degrees /cm +/- 0.8; middle, 1.8 degrees /cm +/- 1.1; bottom, 2.1 degrees /cm +/- 1.3), and greater curvature increase in the bottom third than in the top two thirds (top, 0.15 cm(-1) +/- 0.06; middle, 0.09 cm(-1) +/- 0.07; bottom, 0.41 cm(-1) +/- 0.22; P < .001). CONCLUSIONS The SFA tends to deform more in the bottom third than in the other sections, likely because of less musculoskeletal constraint distal to the adductor canal and vicinity of knee flexion. The SFAs of these older subjects curve off axis with normal joint flexion, probably resulting from known loss of arterial elasticity with age. This slackening of the vessel enables a method for noninvasive quantification of in vivo SFA strain, which may be valuable for treatment planning and device design. In addition, the spatially resolved arterial deformations quantified in this study may be useful for commercial and regulatory device evaluation.

In Vivo Quantification of Murine Aortic Cyclic Strain, Motion, and Curvature: Implications for Abdominal Aortic Aneurysm Growth

To develop methods to quantify cyclic strain, motion, and curvature of the murine abdominal aorta in vivo.

Coronary Artery Axial Plaque Stress and its Relationship With Lesion Geometry : Application of Computational Fluid Dynamics to Coronary CT Angiography

OBJECTIVES The purpose of this study was to characterize the hemodynamic force acting on plaque and to investigate its relationship with lesion geometry. BACKGROUND Coronary plaque rupture occurs when plaque stress exceeds plaque strength. METHODS Computational fluid dynamics was applied to 114 lesions (81 patients) from coronary computed tomography angiography. The axial plaque stress (APS) was computed by extracting the axial component of hemodynamic stress acting on stenotic lesions, and the axial lesion asymmetry was assessed by the luminal radius change over length (radius gradient [RG]). Lesions were divided into upstream-dominant (upstream RG > downstream RG) and downstream-dominant lesions (upstream RG < downstream RG) according to the RG. RESULTS Thirty-three lesions (28.9%) showed net retrograde axial plaque force. Upstream APS linearly increased as lesion severity increased, whereas downstream APS exhibited a concave function for lesion severity. There was a negative correlation (r = -0.274, p = 0.003) between APS and lesion length. The pressure gradient, computed tomography-derived fractional flow reserve (FFRCT), and wall shear stress were consistently higher in upstream segments, regardless of the lesion asymmetry. However, APS was higher in the upstream segment of upstream-dominant lesions (11,371.96 ± 5,575.14 dyne/cm(2) vs. 6,878.14 ± 4,319.51 dyne/cm(2), p < 0.001), and in the downstream segment of downstream-dominant lesions (7,681.12 ± 4,556.99 dyne/cm(2) vs. 11,990.55 ± 5,556.64 dyne/cm(2), p < 0.001). Although there were no differences in FFRCT, % diameter stenosis, and wall shear stress pattern, the distribution of APS was different between upstream- and downstream-dominant lesions. CONCLUSIONS APS uniquely characterizes the stenotic segment and has a strong relationship with lesion geometry. Clinical application of these hemodynamic and geometric indices may be helpful to assess the future risk of plaque rupture and to determine treatment strategy for patients with coronary artery disease. (Evaluation of FFR, WSS, and TPF Using CCTA; NCT01857687).

In vivo deformation of the human abdominal aorta and common iliac arteries with hip and knee flexion: implications for the design of stent-grafts.

Purpose: To quantify in vivo deformations of the abdominal aorta and common iliac arteries (CIAs) caused by musculoskeletal motion. Methods: Seven healthy subjects (age 34±11 years, range 24–50) were imaged in the supine and fetal positions (hip flexion angle 134.0°±9.7°) using contrast-enhanced magnetic resonance angiography. Longitudinal strain, twisting, and curvature change of the infrarenal aorta and CIAs were computed. The angle between the left and right CIAs and translation of the arteries were also computed. Results: Maximal hip flexion induced shortening (5.2%±4.6%), twisting (0.45±0.27 °/mm), and curvature changes (0.015±0.007 mm−1) of the CIAs. The angle between the CIAs increased by 17.6°±8.6°. The iliac arteries moved predominantly in the superior direction relative to the aortic bifurcation, which would induce compression and bending, thus increasing curvature and angle between the CIAs. The abdominal aorta also exhibited shortening (2.9%±2.1%) and twisting (0.07±0.05 °/mm) deformation associated with the hip flexion. Conclusion: Although this study was limited to a few healthy young adults, musculoskeletal motion, specifically hip flexion, caused significant in vivo morphological changes (shortening, twisting, and bending) of the arteries. Predominant superior translation of the CIAs was observed, which suggests that preclinical testing of cyclic superior-inferior translational motion may aid in predicting stent-graft fractures. In turn, stent-graft design could be improved, decreasing overall stent-graft–related complications.

Uncertainty quantification in coronary blood flow simulations: Impact of geometry, boundary conditions and blood viscosity.

Computational fluid dynamic methods are currently being used clinically to simulate blood flow and pressure and predict the functional significance of atherosclerotic lesions in patient-specific models of the coronary arteries extracted from noninvasive coronary computed tomography angiography (cCTA) data. One such technology, FFRCT, or noninvasive fractional flow reserve derived from CT data, has demonstrated high diagnostic accuracy as compared to invasively measured fractional flow reserve (FFR) obtained with a pressure wire inserted in the coronary arteries during diagnostic cardiac catheterization. However, uncertainties in modeling as well as measurement results in differences between these predicted and measured hemodynamic indices. Uncertainty in modeling can manifest in two forms - anatomic uncertainty resulting in error of the reconstructed 3D model and physiologic uncertainty resulting in errors in boundary conditions or blood viscosity. We present a data-driven framework for modeling these uncertainties and study their impact on blood flow simulations. The incompressible Navier-Stokes equations are used to model blood flow and an adaptive stochastic collocation method is used to model uncertainty propagation in the Navier-Stokes equations. We perform uncertainty quantification in two geometries, an idealized stenosis model and a patient specific model. We show that uncertainty in minimum lumen diameter (MLD) has the largest impact on hemodynamic simulations, followed by boundary resistance, viscosity and lesion length. We show that near the diagnostic cutoff (FFRCT=0.8), the uncertainty due to the latter three variables are lower than measurement uncertainty, while the uncertainty due to MLD is only slightly higher than measurement uncertainty. We also show that uncertainties are not additive but only slightly higher than the highest single parameter uncertainty. The method presented here can be used to output interval estimates of hemodynamic indices and visualize patient-specific maps of sensitivities.

Structural Mechanics Predictions Relating to Clinical Coronary Stent Fracture in a 5 Year Period in FDA MAUDE Database.

Endovascular stents are the mainstay of interventional cardiovascular medicine. Technological advances have reduced biological and clinical complications but not mechanical failure. Stent strut fracture is increasingly recognized as of paramount clinical importance. Though consensus reigns that fractures can result from material fatigue, how fracture is induced and the mechanisms underlying its clinical sequelae remain ill-defined. In this study, strut fractures were identified in the prospectively maintained Food and Drug Administration’s (FDA) Manufacturer and User Facility Device Experience Database (MAUDE), covering years 2006–2011, and differentiated based on specific coronary artery implantation site and device configuration. These data, and knowledge of the extent of dynamic arterial deformations obtained from patient CT images and published data, were used to define boundary conditions for 3D finite element models incorporating multimodal, multi-cycle deformation. The structural response for a range of stent designs and configurations was predicted by computational models and included estimation of maximum principal, minimum principal and equivalent plastic strains. Fatigue assessment was performed with Goodman diagrams and safe/unsafe regions defined for different stent designs. Von Mises stress and maximum principal strain increased with multimodal, fully reversed deformation. Spatial maps of unsafe locations corresponded to the identified locations of fracture in different coronary arteries in the clinical database. These findings, for the first time, provide insight into a potential link between patient adverse events and computational modeling of stent deformation. Understanding of the mechanical forces imposed under different implantation conditions may assist in rational design and optimal placement of these devices.

Respiration-induced Deformations of the Superior Mesenteric and Renal Arteries in Patients with Abdominal Aortic Aneurysms

PURPOSE To quantify respiration-induced deformations of the superior mesenteric artery (SMA), left renal artery (LRA), and right renal artery (RRA) in patients with small abdominal aortic aneurysms (AAAs). MATERIALS AND METHODS Sixteen men with AAAs (age 73 y ± 7) were imaged with contrast-enhanced magnetic resonance angiography during inspiratory and expiratory breath-holds. Centerline paths of the aorta and visceral arteries were acquired by geometric modeling and segmentation techniques. Vessel translations and changes in branching angle and curvature resulting from respiration were computed from centerline paths. RESULTS With expiration, the SMA, LRA, and RRA bifurcation points translated superiorly by 12.4 mm ± 9.5, 14.5 mm ± 8.8, and 12.7 mm ± 6.4 (P < .001), and posteriorly by 2.2 mm ± 2.7, 4.9 mm ± 4.2, and 5.6 mm ± 3.9 (P < .05), respectively, and the SMA translated rightward by 3.9 mm ± 4.9 (P < .01). With expiration, the SMA, LRA, and RRA angled upward by 9.7° ± 6.4, 7.5° ± 7.8, and 4.9° ± 5.3, respectively (P < .005). With expiration, mean curvature increased by 0.02 mm(-1) ± 0.01, 0.01 mm(-1) ± 0.01, and 0.01 mm(-1) ± 0.01 in the SMA, LRA, and RRA, respectively (P < .05). For inspiration and expiration, RRA curvature was greater than in other vessels (P < .025). CONCLUSIONS With expiration, the SMA, LRA, and RRA translated superiorly and posteriorly as a result of diaphragmatic motion, inducing upward angling of vessel branches and increased curvature. In addition, the SMA exhibited rightward translation with expiration. The RRA was significantly more tortuous, but deformed less than the other vessels during respiration.

Computational fluid dynamic measures of wall shear stress are related to coronary lesion characteristics

Objective To assess the distribution of pressure and shear-related forces acting on atherosclerotic plaques and their association with lesion characteristics using coronary CT angiography (cCTA)-based computational fluid dynamics (CFD) model of epicardial coronary arteries. Methods Patient-specific models of epicardial coronary arteries were reconstructed from cCTA in 80 patients (12 women, 63.8±9.0 years). The pressure and wall shear stress (WSS) in left anterior descending coronary arteries were assessed using CFD. High-risk plaques were defined as the presence of at least one of the following adverse plaque characteristics: low-density plaque, positive remodelling, napkin-ring sign and spotty calcification. Results At resting condition, 39.5% of stenotic segments (% diameter stenosis 52.3±14.4%) were exposed to high WSS (>40 dyne/cm2). When the stenotic lesion was subdivided into three segments, the distribution of WSS was different from that of pressure change and its magnitude was highest at minimal lumen area (p<0.001). High pressure gradient, proximal location, small lumen and short length were independent determinants of WSS (all p<0.05). The plaques exposed to the highest WSS tertile had a significantly greater proportion of high-risk plaques. The addition of WSS to % diameter stenosis significantly improved the measures of discrimination and reclassification of high-risk plaques (area under the curves from 0.540 to 0.718, p=0.031; net reclassification index 0.827, p<0.001). Conclusions The cCTA-based CFD method can improve the identification of high-risk plaques and the risk stratification for coronary artery disease patients by providing non-invasive measurements of WSS affecting coronary plaques.

Virtual interventions for image-based blood flow computation

Image-based blood flow computation provides great promise for evaluation of vascular devices and assessment of surgical procedures. However, many previous studies employ idealized arterial and device models or only patient-specific models from the image data after device deployment, since the tools for model construction are unavailable or limited and tedious to use. Moreover, in contrast to retrospective studies from existing data, there is a pressing need for prospective analysis with the goal of surgical planning. Therefore, it is necessary to construct models with deployed devices in a fast, virtual and interactive fashion. The goal of this paper is to develop new geometric methods to deploy stents or stent grafts virtually to patient-specific geometric models constructed from a 3D segmentation of medical images. A triangular surface representing the vessel lumen boundary is extracted from the segmentation. The diseased portion is either clipped and replaced by the surface of a deployed device or rerouted in the case of a bypass graft. For diseased arteries close to bifurcations, bifurcated device models are generated. A method to map a 2D strut pattern on the surface of a device is also presented. We demonstrate three applications of our methods in personalized surgical planning for aortic aneurysms, aortic coarctation, and coronary artery stenosis using blood flow computation. Our approach enables prospective model construction and may help to expand the throughput required by routine clinical uses in the future.