A novel method for describing biomechanical properties of the aortic wall based on the three-dimensional fluid-structure interaction model

A Novel Computational Biomechanics Framework to Model Vascular Mechanopropagation in Deep Bone Marrow

The mechanical stimuli generated by body exercise can be transmitted from cortical bone into the deep bone marrow (mechanopropagation). Excitingly, a mechanosensitive perivascular stem cell niche is recently identified within the bone marrow for osteogenesis and lymphopoiesis. Although it is long known that they are maintained by exercise-induced mechanical stimulation, the mechanopropagation from compact bone to deep bone marrow vasculature remains elusive of this fundamental mechanobiology field. No experimental system is available yet to directly understand such exercise-induced mechanopropagation at the bone-vessel interface. To this end, taking advantage of the revolutionary in vivo 3D deep bone imaging, an integrated computational biomechanics framework to quantitatively evaluate the mechanopropagation capabilities for bone marrow arterioles, arteries, and sinusoids is devised. As a highlight, the 3D geometries of blood vessels are smoothly reconstructed in the presence of vessel wall thickness and intravascular pulse pressure. By implementing the 5-parameter Mooney-Rivlin model that simulates the hyperelastic vessel properties, finite element analysis to thoroughly investigate the mechanical effects of exercise-induced intravascular vibratory stretching on bone marrow vasculature is performed. In addition, the blood pressure and cortical bone bending effects on vascular mechanoproperties are examined. For the first time, movement-induced mechanopropagation from the hard cortical bone to the soft vasculature in the bone marrow is numerically simulated. It is concluded that arterioles and arteries are much more efficient in propagating mechanical force than sinusoids due to their stiffness. In the future, this in-silico approach can be combined with other clinical imaging modalities for subject/patient-specific vascular reconstruction and biomechanical analysis, providing large-scale phenotypic data for personalized mechanobiology discovery.

A novel computational fluid dynamic method and validation for assessing distal cerebrovascular microcirculatory resistance

BACKGROUND AND OBJECTIVE

The non-invasive assessment of microcirculatory resistance could improve the treatment of cerebrovascular stenosis. This study aimed to validate a novel computational strategy for determining the reference value of microcirculatory resistance in patients with cerebrovascular stenosis.

METHODS

We reconstructed a patient-specific 3-dimensional model of the extracranial-intracranial arteries. A computational strategy incorporating patient-specific pressure-wire measurements was developed to estimate the blood flow rate and microcirculatory resistance. Throughout the computational fluid dynamics (CFD) simulation, the boundary conditions were adjusted according to the developed algorithm. Pearson correlation and Bland-Altman analyses were used to quantify the correlation and agreement between CFD calculations and transcranial Doppler (TCD) assessment.

RESULTS

A strong correlation was found between the CFD-based and invasive distal pressure measurements (P<0.0001). Meanwhile, the CFD and TCD-based flow measurements were highly correlated (r = 0.853; P = 0.001). Furthermore, there was a correlation between the mean velocity measured by CFD and the mean velocity measured by TCD (r = 0.777; P<0.001). Good agreement was observed between the mass flow by CFD simulation and volumetric flow by TCD (P = 0.0266, mean difference: -0.7814 mmHg, limits of agreement, -4.0905 - 2.5276). However, the mean velocities from CFD simulation were in less agreement with those from the TCD assessment (P = 0.3992, mean difference, -0.0485; limits of agreement, -0.6141 - 0.5170). Results of the CFD simulation indicate that the flow resistance varies greatly between individuals.

CONCLUSIONS

The computational strategy of incorporating patient-specific pressure-wire measurements may serve as an effective approach to evaluate the actual reference values of microcirculatory resistance. In addition, an individualized assessment of non-invasive flow resistance is necessary for the accurate determination of non-invasive cerebrovascular pressure.

Numerical Models Can Assist Choice of an Aortic Phantom for In Vitro Testing

(1) Background: The realization of appropriate aortic replicas for in vitro experiments requires a suitable choice of both the material and geometry. The matching between the grade of details of the geometry and the mechanical response of the materials is an open issue that deserves attention. (2) Methods: To explore this issue, we performed a series of Fluid–Structure Interaction simulations, which compared the dynamics of three aortic models. Specifically, we reproduced a patient-specific geometry with a wall of biological tissue or silicone, and a parametric geometry based on in vivo data made in silicone. The biological tissue and the silicone were modeled with a fiber-oriented anisotropic and isotropic hyperelastic model, respectively. (3) Results: Clearly, both the aorta’s geometry and its constitutive material contribute to the determination of the aortic arch deformation; specifically, the parametric aorta exhibits a strain field similar to the patient-specific model with biological tissue. On the contrary, the local geometry affects the flow velocity distribution quite a lot, although it plays a minor role in the helicity along the arch. (4) Conclusions: The use of a patient-specific prototype in silicone does not a priori ensure a satisfactory reproducibility of the real aorta dynamics. Furthermore, the present simulations suggest that the realization of a simplified replica with the same compliance of the real aorta is able to mimic the overall behavior of the vessel.

Transient numerical simulation of the right coronary artery originating from the left sinus and the effect of its acute take-off angle on hemodynamics

Background

An anomalous origin of the right coronary artery from the left coronary artery sinus is usually characterized by an acute take-off angle. Most affected patients have no clinical symptoms; however, some patients have decreased blood flow into the right coronary artery during exercise, which can lead to symptoms such as myocardial ischemia. Most researchers who have studied an anomalous origin of the right coronary artery from the left coronary artery sinus have done so through clinical cases. In this study, we used numerical simulation to evaluate the hemodynamics of this condition and the effect of an acute take-off angle on hemodynamic parameters. We expect that the results of this study will help in further understanding the clinical symptoms of this anomaly and the hemodynamic impact of an acute take-off angle.

Methods

Three-dimensional models were reconstructed based on the computed tomography images from 16 patients with a normal right coronary artery and 26 patients with an anomalous origin of the right coronary artery from the left coronary artery sinus. A numerical simulation of a two-way fluid-structure interaction was executed with ANSYS Workbench software. The blood was assumed to be an incompressible Newtonian fluid, and the vessel was assumed to be an isotropic, linear elastic material. Hemodynamic parameters and the effect of an acute take-off angle were statistically analyzed.

Results

During the systolic period, the wall pressure in the right coronary artery was significantly reduced in patients with an anomalous origin of the right coronary artery (t =1.32 s, P=0.0001; t =1.34-1.46 s, P<0.0001). The wall shear stress in the abnormal group was higher at the beginning of the systolic period (t =1.24 s, P=0.0473; t =1.26 s, P=0.0193; t =1.28 s, P=0.0441). The acute take-off angle was smaller in patients with clinical symptoms (27.81°±4.406°) than in patients without clinical symptoms (31.86°±2.789°; P=0.017). In the symptomatic group, pressure was negatively correlated with the acute take-off angle (P=0.0185-0.0341, r=-0.459 to -0.4167).

Conclusions

This study shows that an anomalous origin of the right coronary artery from the left coronary artery sinus causes changes in hemodynamic parameters, and that an acute take-off angle in patients with this anomaly is associated with terminal ischemia of the right coronary artery.

Spatial Configuration of Abdominal Aortic Aneurysm Analysis as a Useful Tool for the Estimation of Stent-Graft Migration

The aim of this study was to prepare a self-made mathematical algorithm for the estimation of risk of stent-graft migration with the use of data on abdominal aortic aneurysm (AAA) size and geometry of blood flow through aneurysm sac before or after stent-graft implantation. AngioCT data from 20 patients aged 50–60 years, before and after stent-graft placement in the AAA was analyzed. In order to estimate the risk of stent-graft migration for each patient we prepared an opposite spatial configuration of virtually reconstructed stent-graft with long body or short body. Thus, three groups of 3D geometries were analyzed: 20 geometries representing 3D models of aneurysm, 20 geometries representing 3D models of long body stent-grafts, and 20 geometries representing 3D models of short body stent-graft. The proposed self-made algorithm demonstrated its efficiency and usefulness in estimating wall shear stress (WSS) values. Comparison of the long or short type of stent-graft with AAA geometries allowed to analyze the implants’ spatial configuration. Our study indicated that short stent-graft, after placement in the AAA sac, generated lower drug forces compare to the long stent-graft. Each time shape factor was higher for short stent-graft compare to long stent-graft.

Recent Advances in Biomechanical Characterization of Thoracic Aortic Aneurysms

Thoracic aortic aneurysm (TAA) is a focal enlargement of the thoracic aorta, but the etiology of this disease is not fully understood. Previous work suggests that various genetic syndromes, congenital defects such as bicuspid aortic valve, hypertension, and age are associated with TAA formation. Though occurrence of TAAs is rare, they can be life-threatening when dissection or rupture occurs. Prevention of these adverse events often requires surgical intervention through full aortic root replacement or implantation of endovascular stent grafts. Currently, aneurysm diameters and expansion rates are used to determine if intervention is warranted. Unfortunately, this approach oversimplifies the complex aortopathy. Improving treatment of TAAs will likely require an increased understanding of the biological and biomechanical factors contributing to the disease. Past studies have substantially contributed to our knowledge of TAAs using various ex vivo, in vivo, and computational methods to biomechanically characterize the thoracic aorta. However, any singular approach typically focuses on only material properties of the aortic wall, intra-aneurysmal hemodynamics, or in vivo vessel dynamics, neglecting combinatorial factors that influence aneurysm development and progression. In this review, we briefly summarize the current understanding of TAA causes, treatment, and progression, before discussing recent advances in biomechanical studies of TAAs and possible future directions. We identify the need for comprehensive approaches that combine multiple characterization methods to study the mechanisms contributing to focal weakening and rupture. We hope this summary and analysis will inspire future studies leading to improved prediction of thoracic aneurysm progression and rupture, improving patient diagnoses and outcomes.

Shape and Enhancement Analysis as a Useful Tool for the Presentation of Blood Hemodynamic Properties in the Area of Aortic Dissection

The aim of this study was to create a mathematical approach for blood hemodynamic description with the use of brightness analysis. Medical data was collected from three male patients aged from 45 to 65 years with acute type IIIb aortic dissection that started proximal to the left subclavian artery and involved the renal arteries. For the recognition of wall dissection areas Digital Imaging and Communications in Medicine (DICOM) data were applied. The distance from descending aorta to the diaphragm was analyzed. Each time Feret (D_F) and Hydraulic (D_Hy) diameter were calculated. Moreover, an average brightness (B_AV) was analyzed. Finally, to describe blood hemodynamic in the area of aortic wall dissection, mathematical function combining difference in brightness value and diameter for each computed tomography (CT) scan was calculated. The results indicated that D_F described common duct more accurately compare to D_Hy. While, D_Hy described more accurately true and false ducts. Each time when connection of true and false duct appeared, true duct had lower brightness compare to common duct and false duct. Moreover, false duct characterized with higher brightness compare to common duct. In summary, the proposed algorithm mimics changes in brightness value for patients with acute type IIIb aortic dissection.

A novel vision-based system for quantitative analysis of abdominal aortic aneurysm deformation

BACKGROUND

In clinical diagnostics, combination of different imaging techniques is applied to assess spatial configuration of the abdominal aortic aneurysm (AAA) and deformation of its wall. As deformation of aneurysm wall is crucial parameter in assessing wall rupture, we aimed to develop and validate a Non-Invasive Vision-Based System (NIVBS) for the analysis of 3D elastic artificial abdominal aortic models. 3D-printed elastic AAA models from four patients were applied for the reconstruction of real hemodynamic. During experiments, the inlet boundary conditions included the injection volume and frequency of pulsation averaged from electrocardiography traces. NIVBS system was equipped with nine cameras placed at a constant distance to record wall movement from 360^o angle and a dedicated set of artificial lights providing coherent illumination. Additionally, self-prepared algorithms for image acquisition, processing, segmentation, and contour detection were used to analyze wall deformation. Finally, the shape deformation factor was applied to evaluate aorta's deformation. Experimental results were confronted with medical data from AngioCT and 2D speckle-tracking echocardiography (2DSTE).

RESULTS

Image square analyses indicated that the optimal distance between the camera's lens and the investigated object was in the range of 0.30-0.35 m. There was approximately 1.44% difference observed in aneurysm diameters between NIVBS (86.57 ± 5.86 mm) and AngioCT (87.82 ± 6.04 mm) (p = 0.7764). The accuracy of developed algorithm for the reconstruction of the AAA deformation was equal to 98.56%. Bland-Altman analysis showed that the difference between clinical data (2DSTE) and predicted wall deformation (NIVBS) for all patients was 0.00 mm (confidence interval equal to 0.12 mm) for aneurysm size, 0.01 mm (confidence interval equal to 0.13 mm) and 0.00 mm (confidence interval equal to 0.09 mm) for the anterior and posterior side, as well as 0.01 mm (confidence interval equal to 0.18 mm) and 0.01 mm (confidence interval equal to 0.11 mm) for the left and right side. The optimal range of camera's lens did not affect acquired values.

CONCLUSIONS

The NIVBS with proposed algorithm that reconstructs the pressure from surrounding organs is appropriate to analyze the AAAs in water environment. Moreover, NIVBS allowed detailed quantitative analysis of aneurysm sac wall deformation.