Computational modeling based on confocal imaging predicts changes in osteocyte and dendrite shear stress due to canalicular loss with aging

Exercise and physical activity exert mechanical loading on the bones which induces bone formation. However, the relationship between the osteocyte lacunar-canalicular morphology and mechanical stress experienced locally by osteocytes transducing signals for bone formation is not fully understood. In this study, we used computational modeling to predict the effect of canalicular density, the number of fluid inlets, and load direction on fluid flow shear stress (FFSS) and bone strains and how these might change following the microstructural deterioration of the lacunar-canalicular network that occurs with aging. Four distinct computational models were initially generated of osteocytes with either ten or eighteen dendrites using a fluid-structure interaction method with idealized geometries. Next, a young and a simulated aged osteocyte were developed from confocal images after FITC staining of the femur of a 4-month-old C57BL/6 mouse to estimate FFSS using a computational fluid dynamics approach. The models predicted higher fluid velocities in the canaliculi versus the lacunae. Comparison of idealized models with five versus one fluid inlet indicated that with four more inlets, one-half of the dendrites experienced FFSS greater than 0.8 Pa, which has been associated with osteogenic responses. Confocal image-based models of real osteocytes indicated a six times higher ratio of canalicular to lacunar surface area in the young osteocyte model than the simulated aged model and the average FFSS in the young model (FFSS = 0.46 Pa) was three times greater than the aged model (FFSS = 0.15 Pa). Interestingly, the surface area with FFSS values above 0.8 Pa was 23 times greater in the young versus the simulated aged model. These findings may explain the impaired mechano-responsiveness of osteocytes with aging.

Hemodynamic study of blood flow in the aorta during the interventional robot treatment using fluid-structure interaction

An interventional robot is a means for vascular diagnosis and treatment, and it can perform dredging, releasing drug and operating. Normal hemodynamic indicators are a prerequisite for the application of interventional robots. The current hemodynamic research is limited to the absence of interventional devices or interventional devices in fixed positions. Considering the coupling effect of blood, vessels and robots, based on the bi-directional fluid-structure interaction, using the computational fluid dynamics and particle image velocimetry methods, combined with the sliding and moving mesh technologies, we theoretically and experimentally study the hemodynamic indicators such as blood flow lines, blood pressure, equivalent stress, deformation and wall shear stress of blood vessels when the robot precesses, rotates or does not intervene in the pulsating blood flow. The results show that the intervention of the robot increase the blood flow rate, blood pressure, equivalent stress and deformation of the vessels by 76.4%, 55.4%, 76.5%, and 346%, respectively. The operating mode of the robot during low-speed operation has little impact on the hemodynamic indicators. Using the methyl silicone oil as the experimental fluid, the elastic silicone pipe as the experimental pipe, and the intervention robot having a bioplastic outer shell, the velocity of the fluid around the robot is measured on the developed experimental device for fluid flow field in a pulsating flow when the robot runs. The experimental results are similar to the numerical results. Our work provides an important reference for the hemodynamic study and optimization of the mobile interventional devices.

Fluid-structure interaction modeling of compliant aortic valves using the lattice Boltzmann CFD and FEM methods

The lattice Boltzmann method (LBM) has been increasingly used as a stand-alone CFD solver in various biomechanical applications. This study proposes a new fluid-structure interaction (FSI) co-modeling framework for the hemodynamic-structural analysis of compliant aortic valves. Toward that goal, two commercial software packages are integrated using the lattice Boltzmann (LBM) and finite element (FE) methods. The suitability of the LBM-FE hemodynamic FSI is examined in modeling healthy tricuspid and bicuspid aortic valves (TAV and BAV), respectively. In addition, a multi-scale structural approach that has been employed explicitly recognizes the heterogeneous leaflet tissues and differentiates between the collagen fiber network (CFN) embedded within the elastin matrix of the leaflets. The CFN multi-scale tissue model is inspired by monitoring the distribution of the collagen in 15 porcine leaflets. Different simulations have been examined, and structural stresses and resulting hemodynamics are analyzed. We found that LBM-FE FSI approach can produce good predictions for the flow and structural behaviors of TAV and BAV and correlates well with those reported in the literature. The multi-scale heterogeneous CFN tissue structural model enhances our understanding of the mechanical roles of the CFN and the elastin matrix behaviors. The importance of LBM-FE FSI also emerges in its ability to resolve local hemodynamic and structural behaviors. In particular, the diastolic fluctuating velocity phenomenon near the leaflets is explicitly predicted, providing vital information on the flow transient nature. The full closure of the contacting leaflets in BAV is also demonstrated. Accordingly, good structural kinematics and deformations are captured for the entire cardiac cycle.

Fluid-structure interaction (FSI) simulation for studying the impact of atherosclerosis on hemodynamics, arterial tissue remodeling, and initiation risk of intracranial aneurysms

The biomechanical and hemodynamic effects of atherosclerosis on the initiation of intracranial aneurysms (IA) are not yet clearly discovered. Also, studies for the observation of hemodynamic variation due to atherosclerotic stenosis and its impact on arterial remodeling and aneurysm genesis remain a controversial field of vascular engineering. The majority of studies performed are relevant to computational fluid dynamic (CFD) simulations. CFD studies are limited in consideration of blood and arterial tissue interactions. In this work, the interaction of the blood and vessel tissue because of atherosclerotic occlusions is studied by developing a fluid and structure interaction (FSI) analysis for the first time. The FSI presents a semi-realistic simulation environment to observe how the blood and vessels' structural interactions can increase the accuracy of the biomechanical study results. In the first step, many different intracranial vessels are modeled for an investigation of the biomechanical and hemodynamic effects of atherosclerosis in arterial tissue remodeling. Three physiological conditions of an intact artery, the artery with intracranial atherosclerosis (ICAS), and an atherosclerotic aneurysm (ACA) are employed in the models with required assumptions. Finally, the obtained outputs are studied with comparative and statistical analyses according to the intact model in a normal physiological condition. The results show that existing occlusions in the cross-sectional area of the arteries play a determinative role in changing the hemodynamic behavior of the arterial segments. The undesirable variations in blood velocity and pressure throughout the vessels increase the risk of arterial tissue remodeling and aneurysm formation.

An integrated fluid-structure interaction and thrombosis model for type B aortic dissection

False lumen thrombosis (FLT) in type B aortic dissection has been associated with the progression of dissection and treatment outcome. Existing computational models mostly assume rigid wall behavior which ignores the effect of flap motion on flow and thrombus formation within the FL. In this study, we have combined a fully coupled fluid–structure interaction (FSI) approach with a shear-driven thrombosis model described by a series of convection–diffusion reaction equations. The integrated FSI-thrombosis model has been applied to an idealized dissection geometry to investigate the interaction between vessel wall motion and growing thrombus. Our simulation results show that wall compliance and flap motion can influence the progression of FLT. The main difference between the rigid and FSI models is the continuous development of vortices near the tears caused by drastic flap motion up to 4.45 mm. Flap-induced high shear stress and shear rates around tears help to transport activated platelets further to the neighboring region, thus speeding up thrombus formation during the accelerated phase in the FSI models. Reducing flap mobility by increasing the Young’s modulus of the flap slows down the thrombus growth. Compared to the rigid model, the predicted thrombus volume is 25% larger using the FSI-thrombosis model with a relatively mobile flap. Furthermore, our FSI-thrombosis model can capture the gradual effect of thrombus growth on the flow field, leading to flow obstruction in the FL, increased blood viscosity and reduced flap motion. This model is a step closer toward simulating realistic thrombus growth in aortic dissection, by taking into account the effect of intimal flap and vessel wall motion.

Evaluating the Impact of Calcification on Plaque Vulnerability from the Aspect of Mechanical Interaction Between Blood Flow and Artery Based on MRI

Acute cerebral ischemic events and thrombosis are associated with the rupture/erosion of carotid atherosclerotic plaques. The aim of the present study was to determine the impact of calcification deposition on the wall shear stress (WSS) and stresses within the plaques using 3D fluid-structure interaction (FSI) models. Six patients with calcified carotid atherosclerosis underwent multisequence magnetic resonance imaging (MRI) and were divided into three groups according to the calcification volume. To evaluate the role of the calcification deposition on the stresses, the calcification content was replaced by lipids and arterial tissue, respectively. By comparing the results from the simulation with calcification, and when changing it to lipids there was a significant increment in the stresses at the fibrous cap (p = 0.004). Instead, by changing it to arterial tissue, there was no significant difference (p = 0.07). The calcification shapes that presented the highest stresses were thin concave arc-shaped (AS1) and thin convex arc-shaped (AS3), with mean stress values of 107 ± 54.2 and 99.6 ± 23.4 kPa, respectively. It was also observed that, the calcification shape has more influence on the level of stress than its distance to the lumen. Higher WSS values were associated with the presence of calcification. Calcification shape plays an important role in producing high stresses in the plaque. This work further clarifies the impact of calcification on plaque vulnerability.

Numerical evaluation of transcatheter aortic valve performance during heart beating and its post-deployment fluid-structure interaction analysis

Transcatheter aortic valve replacement (TAVR) is a minimally invasive procedure that provides an effective alternative to open-heart surgical valve replacement for treating advanced calcific aortic valve disease patients. However, complications, such as valve durability, device migration, paravalvular leakage (PVL), and thrombogenicity may lead to increased overall post-TAVR morbidity and mortality. A series of numerical studies involving a self-expandable TAVR valve were performed to evaluate these complications. Structural studies were performed with finite element (FE) analysis, followed by computational fluid dynamics (CFD) simulations, and fluid-structure interaction (FSI) analysis. The FE analysis was utilized to study the effect of TAVR valve implantation depth on valve anchorage in the Living Heart Human Model, which is capable of simulating beating heart during repeated cardiac cycles. The TAVR deployment cases where no valve migration was observed were then used to calculate the post-deployment thrombogenic potential via CFD simulations. FSI analysis followed to further assess the post-deployment TAVR hemodynamic performance for different implantation depths. The deployed valves PVL, geometric and effective orifice areas, and the leaflets structural and flow stress magnitudes were compared to determine the device optimal landing zone. The combined structural and hemodynamic analysis indicated that with the TAVR valve deployed at an aft ventricle position an optimal performance was achieved in the specific anatomy studied. Given the TAVR's rapid expansion to younger lower-risk patients, the comprehensive numerical methodology proposed here can potentially be used as a predictive tool for both procedural planning and valve design optimization to minimize the reported complications.

Search for an Optimal Design of a Bioprosthetic Venous Valve: In silico and in vitro Studies

OBJECTIVE/BACKGROUND

Valve incompetence is a progressive disease of the venous system that may eventually lead to venous hypertension, pain, and ulcers. There is a need for a venous valve prosthesis to replace incompetent valves. Computational and experimental investigations on venous valve design and associated haemodynamics will undoubtedly advance prosthesis design and treatments. Here, the objective is to investigate the effect of venous valve on the fluid and solid mechanics. The hypothesis is that there exists a valve geometry that maximises leaflet shear stress (LSS) but minimises leaflet intramural stress (LIS; i.e., minimise stress ratio = LIS/LSS).

METHODS

To address the hypothesis, fully dynamic fluid-structure interaction (FSI) models were developed. The entire cycle of valve opening and closure was simulated. The flow validation experiments were conducted using a stented venous valve prosthesis and a pulse duplicator flow loop.

RESULTS

Agreement between the output of FSI simulations and output of pulse duplicator was confirmed. The maximum flow rates were within 6% difference, and the total flow during the cycle was within 10% difference. The simulated high stress ratio region at the leaflet base (five times the leaflet average) predicted the disease location of the vast majority of explanted venous valves reported in clinical literature. The study found that the reduced valve height and leaflet dome shape resulted in optimal performance to provide the lowest stress ratio.

CONCLUSION

This study proposes an effective design of venous prostheses and elaborates on the correlations of venous valve with clinical observations.

CardioFAN: open source platform for noninvasive assessment of pulse transit time and pulsatile flow in hyperelastic vascular networks

A profound analysis of pressure and flow wave propagation in cardiovascular systems is the key in noninvasive assessment of hemodynamic parameters. Pulse transit time (PTT), which closely relates to the physical properties of the cardiovascular system, can be linked to variations of blood pressure and stroke volume to provide information for patient-specific clinical diagnostics. In this work, we present mathematical and numerical tools, capable of accurately predicting the PTT, local pulse wave velocity, vessel compliance, and pressure/flow waveforms, in a viscous hyperelastic cardiovascular network. A new one-dimensional framework, entitled cardiovascular flow analysis (CardioFAN), is presented to describe the pulsatile fluid-structure interaction in the hyperelastic arteries, where pertaining hyperbolic equations are solved using a high-resolution total variation diminishing Lax-Wendroff method. The computational algorithm is validated against well-known numerical, in vitro and in vivo data for networks of main human arteries with 55, 37 and 26 segments, respectively. PTT prediction is improved by accounting for hyperelastic nonlinear waves between two arbitrary sections of the arterial tree. Consequently, arterial compliance assignments at each segment are improved in a personalized model of the human aorta and supra-aortic branches with 26 segments, where prior in vivo data were available for comparison. This resulted in a 1.5% improvement in overall predictions of the waveforms, or average relative errors of 5.5% in predicting flow, luminal area and pressure waveforms compared to prior in vivo measurements. The open source software, CardioFAN, can be calibrated for arbitrary patient-specific vascular networks to conduct noninvasive diagnostics.

Numerical simulation of osteocyte cell in response to directional mechanical loadings and mechanotransduction analysis: Considering lacunar-canalicular interstitial fluid flow

The osteocyte cell is a bone cell that also functions as a bone mechanosensor. In this work, a three-dimensional (3D) fluid-structure interaction (FSI) model of an osteocyte cell under different mechanical loading conditions was used to obtain a better understanding of osteocyte cell behavior under different physiological conditions. In the current study, both fluid and solid parts of osteocyte cell were considered in order to allow for more accurate results. Five different loading conditions have been applied to the osteocyte cell, and consequently the different interstitial fluid flow velocities and shear stresses have been investigated. Furthermore, using a mathematical model, the change in the stimulus function value with shear stress and NO enzyme was revealed. This work suggests that changes in osteocyte morphology and direction of loadings affect cell stimulation. It was found that cell is mostly stimulated and expanded in the direction experiencing the most shear stress. Finally, the amount of cell stimulation was shown quantitatively and there was strong dependency between stimulus function, shear stress, calcium, and NO concentration.

A fluid-structure interaction study using patient-specific ruptured and unruptured aneurysm: the effect of aneurysm morphology, hypertension and elasticity

Fluid-structure interaction (FSI) simulations using five patient-specific aneurysm geometries are carried out to investigate the difference between ruptured and unruptured aneurysms. Two different blood pressure conditions (normal and hypertension, for all cases), and two different values of elastic modulus (1 and 2MPa, for two cases) are tested. Ruptured aneurysms (RA) generally displayed larger displacement at the dome, lower area-average WSS and higher von Mises stress than unruptured aneurysms (URA) regardless of elasticity or blood pressure condition. RAs had a longitudinal expansion whereas URAs had a radial expansion, which was the key difference between the two types. The difference in expansion pattern may be one of the keys to explaining aneurysm rupture, and further analysis is required in the future to confirm this theory.

The influence of elastic upstream artery length on fluid-structure interaction modeling: a comparative study using patient-specific cerebral aneurysm

Fluid-structure interaction (FSI) simulations using a patient-specific geometry are carried out to investigate the influence the length of elastic parent artery and the position of constraints in the solid domain on the accuracy of patient-specific FSI simulations. Three models are tested: Long, Moderate, and Short, based on the length of the elastic parent artery. All three models use same wall thickness (0.5 mm) and the elastic modulus (5 MPa). The maximum mesh displacement is the largest for the Long model (0.491 mm) compared to other models (0.3 mm for Moderate, and 0.132 mm for Short). The differences of hemodynamic and mechanical variables, aneurysm volume and cross-sectional area between three models are all found to be minor. In addition, the Short model takes the least amount of computing time of the three models (11h compared to 21 h for Long and 19 h for Moderate). The present results indicate that the use of short elastic upstream artery can shorten the time required for pati ent-specific FSI simulations without impacting the overall accuracy of the results.