Evaluation of dentinal fluid flow behaviours: a fluid-structure interaction simulation

This study uses the fluid-structure interaction (FSI) method to investigate the fluid flow in dental pulp. First, the FSI method is used for the biomechanical simulation of dental intrapulpal responses during force loading (50, 100 and 150 N) on a tooth. The results are validated by comparison with experimental outcomes. Second, the FSI method is used to investigate an intact tooth subjected to a mechanical stimulus during loading at various loading rates. Force loading (0–100 N) is applied gradually to an intact tooth surface with loading rates of 125, 62.5, 25 and 12.5 N/s, respectively, and the fluid flow changes in the pulp are evaluated. FSI analysis is found to be suitable for examining intrapulpal biomechanics. An external force applied to a tooth with a low loading rate leads to a low fluid flow velocity in the pulp chamber, thus avoiding tooth pain.     

Numerical simulation of two-phase non-Newtonian blood flow with fluid-structure interaction in aortic dissection

The behavior of blood cells and vessel compliance significantly influence hemodynamic parameters, which are closely related to the development of aortic dissection. Here the two-phase non-Newtonian model and the fluid-structure interaction (FSI) method are coupled to simulate blood flow in a patient-specific dissected aorta. Moreover, three-element Windkessel model is applied to reproduce physiological pressure waves. Important hemodynamic indicators, such as the spatial distribution of red blood cells (RBCs) and vessel wall displacement, which greatly influence the hemodynamic characteristics are analyzed. Results show that the proximal false lumen near the entry tear appears to be a vortex zone with a relatively lower volume fraction of RBCs, a low time-averaged wall shear stress (TAWSS) and a high oscillatory shear index (OSI), providing a suitable physical environment for the formation of atherosclerosis. The highest TAWSS is located in the narrow area of the distal true lumen which might cause further dilation. TAWSS distributions in the FSI model and the rigid wall model show similar trend, while there is a significant difference for the OSI distributions. We suggest that an integrated model is essential to simulate blood flow in a more realistic physiological environment with the ultimate aim of guiding clinical treatment.     

Mechanical instability of normal and aneurysmal arteries

Tortuous arteries associated with aneurysms have been observed in aged patients with atherosclerosis and hypertension. However, the underlying mechanism is poorly understood. The objective of this study was to determine the effect of aneurysms on arterial buckling instability and the effect of buckling on aneurysm wall stress. We investigated the mechanical buckling and post-buckling behavior of normal and aneurysmal carotid arteries and aorta’s using computational simulations and experimental measurements to elucidate the interrelationship between artery buckling and aneurysms. Buckling tests were done in porcine carotid arteries with small aneurysms created using elastase treatment. Parametric studies were done for model aneurysms with orthotropic nonlinear elastic walls using finite element simulations. Our results demonstrated that arteries buckled at a critical buckling pressure and the post-buckling deflection increased nonlinearly with increasing pressure. The presence of an aneurysm can reduce the critical buckling pressure of arteries, although the effect depends on the aneurysm’s dimensions. Buckled aneurysms demonstrated a higher peak wall stress compared to unbuckled aneurysms under the same lumen pressure. We conclude that aneurysmal arteries are vulnerable to mechanical buckling and mechanical buckling could lead to high stresses in the aneurysm wall. Buckling could be a possible mechanism for the development of tortuous aneurysmal arteries such as in the Loeys–Dietz syndrome.     

Smart platform for the analysis of cupula deformation caused by otoconia presence within SCCs

In this paper, we first briefly describe the mechanical model of cupula deformation with the appropriate analytical solution. Then, we present the numerical solution of the same problem and compare it with the analytical one. Besides, we provide another numerical solution based on the Finite Element Method procedure, in an effort to encompass a more realistic approach to the problem such as considering the real geometry of the SCCs and the obstruction of the fluid flow during head movement due to the presence of otoconia. As a result, we obtain fifty solutions for a head rotation angle in a range from 0° to 120°, taking into account that such a manoeuvre lasts exactly 3?seconds. In the end, we present a mobile client–server application for visualising the finite element solutions in a way that is convenient for the clinical practice.     

Influence of flow rate and scaffold pore size on cell behavior during mechanical stimulation in a flow perfusion bioreactor

Mechanically stimulating cell-seeded scaffolds by flow-perfusion is one approach utilized for developing clinically applicable bone graft substitutes. A key challenge is determining the magnitude of stimuli to apply that enhances cell differentiation but minimizes cell detachment from the scaffold. In this study, we employed a combined computational modeling and experimental approach to examine how the scaffold mean pore size influences cell attachment morphology and subsequently impacts upon cell deformation and detachment when subjected to fluid-flow. Cell detachment from osteoblast-seeded collagen-GAG scaffolds was evaluated experimentally across a range of scaffold pore sizes subjected to different flow rates and exposure times in a perfusion bioreactor. Cell detachment was found to be proportional to flow rate and inversely proportional to pore size. Using this data, a theoretical model was derived that accurately predicted cell detachment as a function of mean shear stress, mean pore size, and time. Computational modeling of cell deformation in response to fluid flow showed the percentage of cells exceeding a critical threshold of deformation correlated with cell detachment experimentally and the majority of these cells were of a bridging morphology (cells stretched across pores). These findings will help researchers optimize the mean pore size of scaffolds and perfusion bioreactor operating conditions to manage cell detachment when mechanically simulating cells via flow perfusion.     

Development and Experimental Validationof a Flmd/Structure-Interaction Finite Element Modelof a Vacuum-Driven Cell Culture Mechanostimulus System

A new fluid/structure-interaction finite element formulation is reported, by means of which reactive fluid stresses can be determined for what is currently the most widely used laboratory apparatus (; the Flexercell Strain Unit ) for delivering controlled in vitro mechanical stimuli to cultured cells. The apparatus functions by means of cyclic vacuum application to the undersurface of a membrane-like circular rubber substrate. When operated in Us original embodiment ( i.e., without axial constraint to substrate motion), the pulsatile vacuum causes appreciable pulsatile excursions ( often several millimeters) of the substrate. The mechanical stimuli experienced by cells attached atop the substrate include not only substrate distention, but also potentially confounding reactive fluid stresses due to coupled motions of the overlying liquid culture nutrient medium. Since it is impractical to directly measure reactive fluid stress in such environments, a corresponding mathematical model has been developed. The formulation involves transient continuum finite element solutions for the nutrient medium flow field and for the deformation of the substrate, coupled at their mutual interface ( the substrate culture surface) Besides the nonlinearities inherent in the flow field and substrate treatments per se, the numerical problem is complicated by the presence of moving boundaries at the nutrient free surface and at the nutrient/substrate interface, as well as by the need to enforce fluid/structure interaction throughout the duty cycle. Algorithmic considerations appropriate to achieving physically realistic numerical performance are reported, and a confirmatory laboratory validation experiment is described.     

Stratification of a population of intracranial aneurysms using blood flow metrics

Indices of the intra-aneurysm hemodynamic environment have been proposed as potentially indicative of their longitudinal outcome. To be useful, the indices need to be used to stratify large study populations and tested against known outcomes. The first objective was to compile the diverse hemodynamic indices reported in the literature. Furthermore, as morphology is often the only patient-specific information available in large population studies, the second objective was to assess how the ranking of aneurysms in a population is affected by the use of steady flow simulation as an approximation to pulsatile flow simulation, even though the former is clearly non-physiological. Sixteen indices of aneurysmal hemodynamics reported in the literature were compiled and refined where needed. It was noted that, in the literature, these global indices of flow were always time-averaged over the cardiac cycle. Steady and pulsatile flow simulations were performed on a population of 198 patient-specific and 30 idealised aneurysm models. All proposed hemodynamic indices were estimated and compared between the two simulations. It was found that steady and pulsatile flow simulations had a strong linear dependence (r ≥ 0.99 for 14 indices; r ≥ 0.97 for 2 others) and rank the aneurysms in an almost identical fashion (ρ ≥ 0.99 for 14 indices; ρ ≥ 0.96 for other 2). When geometry is the only measured piece of information available, stratification of aneurysms based on hemodynamic indices reduces to being a physically grounded substitute for stratification of aneurysms based on morphology. Under such circumstances, steady flow simulations may be just as effective as pulsatile flow simulation for estimating most key indices currently reported in the literature.     

A numerical analysis of the aortic blood flow pattern during pulsed cardiopulmonary bypass

In the modern era, stroke remains a main cause of morbidity after cardiac surgery despite continuing improvements in the cardiopulmonary bypass (CPB) techniques. The aim of the current work was to numerically investigate the blood flow in aorta and epiaortic vessels during standard and pulsed CPB, obtained with the intra-aortic balloon pump (IABP). A multi-scale model, realized coupling a 3D computational fluid dynamics study with a 0D model, was developed and validated with in vivo data. The presence of IABP improved the flow pattern directed towards the epiaortic vessels with a mean flow increase of 6.3% and reduced flow vorticity.     

Influence of cortical canal architecture on lacunocanalicular pore pressure and fluid flow

Bone is a dynamic tissue that undergoes structural modification in response to its mechanical environment, but how bone cells sense and respond to loading conditions remains incompletely understood. Current theories focus on strain-induced fluid flow for the primary means of mechanotransduction. To examine the influence of age-related cortical rarefaction on lacunocanalicular fluid characteristics, coupled fluid flow and mechanical computational models of bone specimens representing young, mid-age and aged samples were derived artificially from the same original micro-computed tomography image data. Simulated mechanical loading was applied to the bone models to induce pressure-driven interstitial fluid flow. Results demonstrated a decrease in pore pressure and fluid velocity magnitudes with age as a result of increased cortical porosity. Mean canal separation, as opposed to canal size, was implicated as a primary factor affecting age-related fluid dynamics. Future investigations through refinement of the model may implicate fluid stasis or inadequate nutrient transport experienced by osteocytes as a key factor in the initiation of cortical remodelling events.     

Effectiveness of pedicle screw inclusion at the fracture level in short-segment fixation constructs for the treatment of thoracolumbar burst fractures: a computational biomechanics analysis

When treating thoracolumbar burst fractures (BF), short-segment posterior fixation (SSPF) represents a less invasive alternative to the traditional long-segment posterior fixation (LSPF) approach. However, hardware failure and loss of sagittal alignment have been reported in patients treated with SSPF. Including pedicle screws at the fracture level in SSPF constructs has been proposed to improve stiffness and reliability of the construct. Accordingly, the biomechanical performance of the proposed construct was compared to LSPF via a computational analysis. Pedicle screws at fracture level improved the performance of the short-segment construct. However, LSPF still represent a biomechanically superior option for treating thoracolumbar BF.     

Computational fluid dynamics in abdominal aorta bifurcation: non-Newtonian versus Newtonian blood flow in a real case study

Hemodynamic in abdominal aorta bifurcation was investigated in a real case using computational fluid dynamics. A Newtonian and non-Newtonian (Walburn-Schneck) viscosity models were compared. The geometrical model was obtained by 3D reconstruction from CT-scan and hemodynamic parameters obtained by laser-Doppler. Blood was assumed incompressible fluid, laminar flow in transient regime and rigid vessel wall. Finite volume-based was used to study the velocity, pressure, wall shear stress (WSS) and viscosity throughout cardiac cycle. Results obtained with Walburn-Schneck’s model, during systole, present lower viscosity due to shear thinning behavior. Furthermore, there is a significant difference between the results obtained by the two models for a specific patient. During the systole, differences are more pronounced and are preferably located in the tortuous regions of the artery. Throughout the cardiac cycle, the WSS amplitude between the systole and diastole is greater for the Walburn-Schneck’s model than for the Newtonian model. However, the average viscosity along the artery is always greater for the non-Newtonian model, except in the systolic peak. The hemodynamic model is crucial to validate results obtained with CFD and to explore clinical potential.     

Transient blood flow in elastic coronary arteries with varying degrees of stenosis and dilatations: CFD modelling and parametric study

In this paper, we have analysed pulsatile flow through partially occluded elastic arteries, to determine the haemodynamic parameters of wall shear stress (WSS), wall pressure gradient and pressure drops (ΔP), contributing to enhanced flow resistance and myocardial ischaemic regions which impair cardiac contractility and cause increased work load on the heart. In summary, it can be observed that stenoses in an artery significantly influence the haemodynamic parameters of wall shear stress and pressure drop in contrast to dilatations case. This deduces that stenosis plays a more critical role in plaque growth and vulnerability in contrast to dilatation, and should be the key element in cardiovascular pathology and diagnosis. Through quantitative analysis of WSS and ΔP, we have provided a clearer insight into the haemodynamics of atherosclerotic arteries. Determination of these parameters can be helpful to cardiologists, because it is directly implicated in the genesis and development of atherosclerosis.     

Theoretical analysis of magnetic field interactions with aortic blood flow

The flow of blood in the presence of a magnetic field gives rise to induced voltages in the major arteries of the central circulatory system. Under certain simplifying conditions, such as the assumption that the length of major arteries (e.g., the aorta) is infinite and that the vessel walls are not electrically conductive, the distribution of induced voltages and currents within these blood vessels can be calculated with reasonable precision. However, the propagation of magnetically induced voltages and currents from the aorta into neighboring tissue structures such as the sinuatrial node of the heart has not been previously determined by any experimental or theoretical technique. In the analysis presented in this paper, a solution of the complete Navier-Stokes equation was obtained by the finite element technique for blood flow through the ascending and descending aortic vessels in the presence of a uniform static magnetic field. Spatial distributions of the magnetically induced voltage and current were obtained for the aortic vessel and surrounding tissues under the assumption that the wall of the aorta is electrically conductive. Results are presented for the calculated values of magnetically induced voltages and current densities in the aorta and surrounding tissue structures, including the sinuatrial node, and for their field-strength dependence. In addition, an analysis is presented of magnetohydrodynamic interactions that lead to a small reduction of blood volume flow at high field levels above approximately 10 tesla (T). Quantitative results are presented on the offsetting effects of oppositely directed blood flows in the ascending and descending aortic segments, and a quantitative estimate is made of the effects of assuming an infinite vs. a finite length of the aortic vessel in calculating the magnetically induced voltage and current density distribution in tissue. © 1996 Wiley-Liss, Inc.     

Numerical analysis of pulsatile blood flow and vessel wall mechanics in different degrees of stenoses

Hemodynamics factors and biomechanical forces play key roles in atherogenesis, plaque development and final rupture. In this paper, we investigated the flow field and stress field for different degrees of stenoses under physiological conditions. Disease is modelled as axisymmetric cosine shape stenoses with varying diameter reductions of 30%, 50% and 70%, respectively. A simulation model which incorporates fluid-structure interaction, a turbulence model and realistic boundary conditions has been developed. The results show that wall motion is constrained at the throat by 60% for the 30% stenosis and 85% for the 50% stenosis; while for the 70% stenosis, wall motion at the throat is negligible through the whole cycle. Peak velocity at the throat varies from 1.47 m/s in the 30% stenosis to 3.2m/s in the 70% stenosis against a value of 0.78 m/s in healthy arteries. Peak wall shear stress values greater than 100 Pa were found for > or =50% stenoses, which in vivo could lead to endothelial stripping. Maximum circumferential stress was found at the shoulders of plaques. The results from this investigation suggest that severe stenoses inhibit wall motion, resulting in higher blood velocities and higher peak wall shear stress, and localization of hoop stress. These factors may contribute to further development and rupture of plaques.     

Computational biomechanical analysis of postoperative inferior tibiofibular syndesmosis: a modified modeling method

To analyze the biomechanical effect of syndesmotic screw through three and four cortices, a total of 12 finite element models simulating healthy ankles, tibiofibular syndesmosis injured ankles, and post-operative ankles by screw fixations through three or four cortices were built. A set of biomechanical data were obtained to find that screw fixation methods for inferior tibiofibular syndesmosis can help recover most of the biomechanical relations of the ankle especially the tricortical fixation, while the screw of quadricortical fixation bear more stress than the tricortical fixation. The modeling method for finite element models was also modified for saving more time and realizing personalized modeling for clinical application.     

Comparative study of flow in right-sided and left-sided aortas: numerical simulations in patient-based models

A right-sided aorta is a rare malformation which may be associated with other various types of congenital heart disease. We utilised haemodynamic, echocardiographic measurements, computerised tomography and image reconstruction software packages that were integrated in a computational fluid dynamics model to determine blood flow patterns in patient-based aortas. In the left-sided aorta, a systolic clockwise rotational component was present, while helical flow was depicted in the aortic arch that was converted in the descending aorta as counter-rotating vortices with accompanying retrograde flow. The right-sided configuration has not altered the orientation of the three-dimensional vortex, but intensification of polymorphic flow patterns, alterations in wall shear stress distribution and development of a lateral pressure gradient at the area of an aneurysmal anomaly was observed. Moreover, increments of Reynolds, Womersley and Dean numbers were evident. These phenomena along with the formation of the aneurysm might influence cardiovascular risk in patients with right-sided aortas.