A Dual-VENC Four-Dimensional Flow MRI Framework for Analysis of Subject-Specific Heterogeneous Nonlinear Vessel Deformation.

Advancement of subject-specific in silico medicine requires new imaging protocols tailored to specific anatomical features, paired with new constitutive model development based on structure/function relationships. In this study, we develop a new dual-velocity encoding coefficient (VENC) 4D flow MRI protocol that provides unprecedented spatial and temporal resolution of in vivo aortic deformation. All previous dual-VENC 4D flow MRI studies in the literature focus on an isolated segment of the aorta, which fail to capture the full spectrum of aortic heterogeneity that exists along the vessel length. The imaging protocol developed provides high sensitivity to all blood flow velocities throughout the entire cardiac cycle, overcoming the challenge of accurately measuring the highly unsteady nonuniform flow field in the aorta. Cross-sectional area change, volumetric flow rate, and compliance are observed to decrease with distance from the heart, while pulse wave velocity (PWV) is observed to increase. A nonlinear aortic lumen pressure-area relationship is observed throughout the aorta such that a high vessel compliance occurs during diastole, and a low vessel compliance occurs during systole. This suggests that a single value of compliance may not accurately represent vessel behavior during a cardiac cycle in vivo. This high-resolution MRI data provide key information on the spatial variation in nonlinear aortic compliance, which can significantly advance the state-of-the-art of in-silico diagnostic techniques for the human aorta.

Quantification of the regional bioarchitecture in the human aorta.

Regional variance in human aortic bioarchitecture responsible for the elasticity of the vessel is poorly understood. The current study quantifies the elements responsible for aortic compliance, namely, elastin, collagen and smooth muscle cells, using histological and stereological techniques on human tissue with a focus on regional heterogeneity. Using donated cadaveric tissue, a series of samples were excised between the proximal ascending aorta and the distal abdominal aorta, for five cadavers, each of which underwent various staining procedures to enhance specific constituents of the wall. Using polarised light microscopy techniques, the orientation of collagen fibres was studied for each location and each tunical layer of the aorta. Significant transmural and longitudinal heterogeneity in collagen fibre orientations were uncovered throughout the vessel. It is shown that a von Mises mixture model is required accurately to fit the complex collagen fibre distributions that exist along the aorta. Additionally, collagen and smooth muscle cell density was observed to increase with increasing distance from the heart, whereas elastin density decreased. Evidence clearly demonstrates that the aorta is a highly heterogeneous vessel which cannot be simplistically represented by a single compliance value. The quantification and fitting of the regional aortic bioarchitectural data, although not without its limitations, including mean cohort age of 77.6 years, facilitates the development of next-generation finite element models that can potentially simulate the influence of regional aortic composition and microstructure on vessel biomechanics.

Review of Mechanical Testing and Modelling of Thrombus Material for Vascular Implant and Device Design.

A thrombus or blood clot is a solid mass, made up of a network of fibrin, platelets and other blood components. Blood clots can form through various pathways, for example as a result of exposed tissue factor from vascular injury, as a result of low flow/stasis, or in very high shear flow conditions. Embolization of cardiac or vascular originating blood clots, causing an occlusion of the neurovasculature, is the major cause of stroke and accounts for 85% of all stroke. With mechanical thrombectomy emerging as the new standard of care in the treatment of acute ischemic stroke (AIS), the need to generate a better understanding of the biomechanical properties and material behaviour of thrombus material has never been greater, as it could have many potential benefits for the analysis and performance of these treatment devices. Defining the material properties of a thrombus has obvious implications for the development of these treatment devices. However, to-date this definition has not been adequately established. While some experimentation has been performed, model development has been extremely limited. This paper reviews the previous literature on mechanical testing of thrombus material. It also explores the use of various constitutive and computational models to model thrombus formation and material behaviour.

Numerical Simulation of Stent Angioplasty with Predilation: An Investigation into Lesion Constitutive Representation and Calcification Influence.

It is acceptable clinical practice to predilate a severely occluded vessel to allow better positioning of endovascular stents, and while the impact of this intervention has been examined for aggregate response in animals there has been no means to examine whether there are specific vessels that might benefit. Finite element methods offer the singular ability to explore the mechanical response of arteries with specific pathologic alterations in mechanics to stenting and predilation. We examined varying representations of atherosclerotic tissue including homogeneous and heterogeneous dispersion of calcified particles, and elastic, pseudo-elastic, and elastic-plastic constitutive representations of bulk atherosclerotic tissue. The constitutive representations of the bulk atherosclerotic tissue were derived from experimental test data and highlight the importance of accounting for testing mode of loading. The impact of arterial predilation is presented and, in particular, its effect on intimal predicted damage, atherosclerotic tissue von Mises and maximum principal stresses, and luminal deformation was dependent on the type of constitutive representation of diseased tissue, particularly in the presence of calcifications.

Evaluation of cover effects on bare stent mechanical response.

Covered tracheobronchial stents are used to prevent tumour growth from reoccluding the airways. In the present work a combination of experimental and computational methods are used to present the mechanical effects that adhered covers can have on stent performance. A prototype tracheobronchial stent is characterised in bare and covered configurations using radial force, flat plate and a novel non-uniform radial force test, while computational modelling is performed in parallel to extensively inform the physical testing. Results of the study show that cover configuration can have a significant structural effect on stent performance, and that stent response (bare or covered) is especially loading specific, highlighting that the loading configuration that a stent is about to be subjected to should be considered before stent implantation.

Webbing and Delamination of Drug Eluting Stent Coatings.

The advancement of the drug-eluting stent technology raises the significant challenge of safe mechanical design of polymer coated stent systems. Experimental images of stent coatings undergoing significant damage during deployment have been reported; such coating damage and delamination can lead to complications such as restenosis and increased thrombogenicity. In the current study a cohesive zone modeling framework is developed to predict coating delamination and buckling due to hinge deformation during stent deployment. Models are then extended to analyze, for the first time, stent-coating damage due to webbing defects. Webbing defects occur when a bond forms between coating layers on adjacent struts, resulting in extensive delamination of the coating from the strut surfaces. The analyzes presented in this paper uncover the mechanical factors that govern webbing induced coating damage. Finally, an experimental fracture test of a commercially available stent coating material is performed and results demonstrate that the high cohesive strength of the coating material will prevent web fracture, resulting in significant coating delamination during stent deployment.

An Experimental and Computational Investigation of Bone Formation in Mechanically Loaded Trabecular Bone Explants.

Understanding how bone marrow multipotent stromal cells (MSCs) contribute to new bone formation and remodeling in vivo is of principal importance for informing the development of effective bone tissue engineering strategies in vitro. However, the precise in situ stimuli that MSCs experience have not been fully established. The shear stress generated within the bone marrow of physiologically loaded samples has never been determined, but could be playing an important role in the generation of sufficient stimulus for MSCs to undergo osteogenic differentiation. In this study fluid structure interaction (FSI) computational models were used in conjunction with a bioreactor which physiologically compresses explanted trabecular bone samples to determine whether MSCs can be directly stimulated by mechanical cues within the bone marrow. Experimentally loaded samples were found to have greater osteogenic activity, as verified by bone histomorphometry, compared to control static samples. FSI models demonstrated a linear relationship between increasing shear stress and decreasing bone volume. The FSI models demonstrated that bone strain, not marrow shear stress, was likely the overall driving mechanical signal for new bone formation during compression. However, the shear stress generated in the models is within the range of values which has been shown previously to generate an osteogenic response in MSCs.

Computational Bench Testing to Evaluate the Short-Term Mechanical Performance of a Polymeric Stent.

Over the last decade, there has been a significant volume of research focussed on the utilization of biodegradable polymers such as poly-L-lactide-acid (PLLA) for applications associated with cardiovascular disease. More specifically, there has been an emphasis on upgrading current clinical shortfalls experienced with conventional bare metal stents and drug eluting stents. One such approach, the adaption of fully formed polymeric stents has led to a small number of products being commercialized. Unfortunately, these products are still in their market infancy, meaning there is a clear non-occurrence of long term data which can support their mechanical performance in vivo. Moreover, the load carry capacity and other mechanical properties essential to a fully optimized polymeric stent are difficult, timely and costly to establish. With the aim of compiling rapid and representative performance data for specific stent geometries, materials and designs, in addition to reducing experimental timeframes, Computational bench testing via finite element analysis (FEA) offers itself as a very powerful tool. On this basis, the research presented in this paper is concentrated on the finite element simulation of the mechanical performance of PLLA, which is a fully biodegradable polymer, in the stent application, using a non-linear viscous material model. Three physical stent geometries, typically used for fully polymeric stents, are selected, and a comparative study is performed in relation to their short-term mechanical performance, with the aid of experimental data. From the simulated output results, an informed understanding can be established in relation to radial strength, flexibility and longitudinal resistance, that can be compared with conventional permanent metal stent functionality, and the results show that it is indeed possible to generate a PLLA stent with comparable and sufficient mechanical performance. The paper also demonstrates the attractiveness of FEA as a tool for establishing fundamental mechanical characteristics of polymeric stent performance.

Micro-scale testing and micromechanical modelling for high cycle fatigue of CoCr stent material.

This paper presents a framework of experimental testing and crystal plasticity micromechanics for high cycle fatigue (HCF) of micro-scale L605 CoCr stent material. Micro-scale specimens, representative of stent struts, are manufactured via laser micro-machining and electro-polishing from biomedical grade CoCr alloy foil. Crystal plasticity models of the micro-specimens are developed using a length scale-dependent, strain-gradient constitutive model and a phenomenological (power-law) constitutive model, calibrated from monotonic and cyclic plasticity test data. Experimental microstructural characterisation of the grain morphology and precipitate distributions is used as input for the polycrystalline finite element (FE) morphologies. Two microstructure-sensitive fatigue indicator parameters are applied, using local and non-local (grain-averaged) implementations, for the phenomenological and length scale-dependent models, respectively, to predict fatigue crack initiation (FCI) in the HCF experiments.

Improving the finite element model accuracy of tissue engineering scaffolds produced by selective laser sintering.

In bone tissue engineering, both geometrical and mechanical properties of a scaffold play a major part in the success of the treatment. The mechanical stresses and strains that act on cells on a scaffold in a physiological environment are a determining factor on the subsequent tissue formation. Computational models are often used to simulate the effect of changes of internal architectures and external loads applied to the scaffold in order to optimise the scaffold geometry for the prospective implantation site. Finite element analysis (FEA) based on computer models of the scaffold is a common technique, but would not take into account actual inaccuracies due to the manufacturing process. Image based FEA using CT scans of fabricated scaffolds can provide a more accurate analysis of the scaffold, and was used in this work in order to accurately simulate and predict the mechanical performance of bone tissue engineering scaffolds, fabricated using selective laser sintering (SLS), with a view to generating a methodology that could be used to optimise scaffold design. The present work revealed that an approach that assumes isotropic properties of SLS fabricated scaffolds will lead to inaccurate predictions of the FE model. However, a dependency of the grey value of the CT scans and the mechanical properties was discovered, which may ultimately lead to accurate FE models without the need of experimental validation.

Thermo-electrical equivalents for simulating the electro-mechanical behavior of biological tissue.

Equivalence is one of most popular techniques to simulate the behavior of systems governed by the same type of differential equation. In this case, a thermo-electrical equivalence is considered as a method for modelling the inter-dependence of electrical and mechanical phenomena in biological tissue. We seek to assess this approach for multi-scale models (from micro-structure to tissue scale) of biological media, such as nerve cells and cardiac tissue, in which the electrical charge distribution is modelled as a heat distribution in an equivalent thermal system. This procedure allows for the reduction in problem complexity and it facilitates the coupling of electrical and mechanical phenomena in an efficient and practical way. Although the findings of this analysis are mainly addressed towards the electro-mechanics of tissue within the biomedical domain, the same approach could be used in other studies in which a coupled finite element analysis is required.

Mechanical stimulation of bone marrow in situ induces bone formation in trabecular explants.

Low magnitude high frequency (LMHF) loading has been shown to have an anabolic effect on trabecular bone in vivo. However, the precise mechanical signal imposed on the bone marrow cells by LMHF loading, which induces a cellular response, remains unclear. This study investigates the influence of LMHF loading, applied using a custom designed bioreactor, on bone adaptation in an explanted trabecular bone model, which isolated the bone and marrow. Bone adaptation was investigated by performing micro CT scans pre and post experimental LMHF loading, using image registration techniques. Computational fluids dynamic models were generated using the pre-experiment scans to characterise the mechanical stimuli imposed by the loading regime prior to adaptation. Results here demonstrate a significant increase in bone formation in the LMHF loaded group compared to static controls and media flow groups. The calculated shear stress in the marrow was between 0.575 and 0.7 Pa, which is within the range of stimuli known to induce osteogenesis by bone marrow mesenchymal stem cells in vitro. Interestingly, a correlation was found between the bone formation balance (bone formation/resorption), trabecular number, trabecular spacing, mineral resorption rate, bone resorption rate and mean shear stresses. The results of this study suggest that the magnitude of the shear stresses generated due to LMHF loading in the explanted bone cores has a contributory role in the formation of trabecular bone and improvement in bone architecture parameters.

Nitinol stent design - understanding axial buckling.

Nitinol׳s superelastic properties permit self-expanding stents to be crimped without plastic deformation, but its nonlinear properties can contribute towards stent buckling. This study investigates the axial buckling of a prototype tracheobronchial nitinol stent design during crimping, with the objective of eliminating buckling from the design. To capture the stent buckling mechanism a computational model of a radial force test is simulated, where small geometric defects are introduced to remove symmetry and allow buckling to occur. With the buckling mechanism ascertained, a sensitivity study is carried out to examine the effect that the transitional plateau region of the nitinol loading curve has on stent stability. Results of this analysis are then used to redesign the stent and remove buckling. It is found that the transitional plateau region can have a significant effect on the stability of a stent during crimping, and by reducing the amount of transitional material within the stent hinges during loading the stability of a nitinol stent can be increased.

Modelling of atherosclerotic plaque for use in a computational test-bed for stent angioplasty.

A thorough understanding of the diseased tissue state is necessary for the successful treatment of a blocked arterial vessel using stent angioplasty. The constitutive representation of atherosclerotic tissue is of great interest to researchers and engineers using computational models to analyse stents, as it is this in silico environment that allows extensive exploration of tissue response to device implantation. This paper presents an in silico evaluation of the effects of variation of atherosclerotic tissue constitutive representation on tissue mechanical response during stent implantation. The motivation behind this work is to investigate the level of detail that is required when modelling atherosclerotic tissue in a stenting simulation, and to give recommendations to the FDA for their guideline document on coronary stent evaluation, and specifically the current requirements for computational stress analyses. This paper explores the effects of variation of the material model for the atherosclerotic tissue matrix, the effects of inclusion of calcifications and a lipid pool, and finally the effects of inclusion of the Mullins effect in the atherosclerotic tissue matrix, on tissue response in stenting simulations. Results indicate that the inclusion of the Mullins effect in a direct stenting simulation does not have a significant effect on the deformed shape of the tissue or the stress state of the tissue. The inclusion of a lipid pool induces a local redistribution of lesion deformation for a soft surrounding matrix and the inclusion of a small volume of calcifications dramatically alters the local results for a soft surrounding matrix. One of the key findings from this work is that the underlying constitutive model (elasticity model) used for the atherosclerotic tissue is the dominant feature of the tissue representation in predicting tissue response in a stenting simulation.

Computational micromechanics of bioabsorbable magnesium stents.

Magnesium alloys are a promising candidate material for an emerging generation of absorbable metal stents. Due to its hexagonal-close-packed lattice structure and tendency to undergo twinning, the deformation behaviour of magnesium is quite different to that of conventional stent materials, such as stainless steel 316L and cobalt chromium L605. In particular, magnesium exhibits asymmetric plastic behaviour (i.e. different yield behaviours in tension and compression) and has lower ductility than these conventional alloys. In the on-going development of absorbable metal stents it is important to assess how the unique behaviour of magnesium affects device performance. The mechanical behaviour of magnesium stent struts is investigated in this study using computational micromechanics, based on finite element analysis and crystal plasticity theory. The plastic deformation in tension and bending of textured and non-textured magnesium stent struts with different numbers of grains through the strut dimension is investigated. It is predicted that, unlike 316L and L605, the failure risk and load bearing capacity of magnesium stent struts during expansion is not strongly affected by the number of grains across the strut dimensions; however texturing, which may be introduced and controlled in the manufacturing process, is predicted to have a significant influence on these measures of strut performance.

A physical corrosion model for bioabsorbable metal stents.

Absorbable metal stents (AMSs) are an emerging technology in the treatment of heart disease. Computational modelling of AMS performance will facilitate the development of this technology. In this study a physical corrosion model is developed for AMSs based on the finite element method and adaptive meshing. The model addresses a gap between currently available phenomenological corrosion models for AMSs and physical corrosion models that have been developed for more simple geometries than those of a stent. The model developed in this study captures the changing surface of a corroding three-dimensional AMS structure for the case of diffusion-controlled corrosion. Comparisons are made between model predictions and those of previously developed phenomenological corrosion models for AMSs in terms of predicted device geometry and mechanical performance during corrosion. Relationships between alloy solubility and diffusivity in the corrosion environment and device performance during corrosion are also investigated.

A novel flow chamber for biodegradable alloy assessment in physiologically realistic environments.

In order to better understand the in vivo corrosion of biodegradable alloys, it is necessary to replicate the physiological environment as closely as possible. In this study, a novel flow chamber system is developed that allows the investigation of biodegradable alloy corrosion in a simulated physiological environment. The system is designed to reproduce flow conditions encountered in coronary arteries using a parallel plate setup and to allow the culturing of cells. Computational fluid dynamics and analytical methods are used as part of the design process to ensure that suitable flow conditions are maintained in the test region. The system is used to investigate the corrosion behavior of AZ31 alloy foils of different thickness, in test media with and without proteins and in static and dynamic solutions. It is observed that pulsatile flows, similar to those in the coronary arteries, significantly increase corrosion rates and lead to a different corrosion surface morphologies relative to static immersion tests.

Influence of statistical size effects on the plastic deformation of coronary stents.

The dimensions of coronary stent struts are similar to those of the metallic grains of their constituent alloys. This means that statistical size effects (SSEs), which are evident in polycrystals with few grains through their dimensions, can have detrimental effects on the mechanical performance of stent struts undergoing large plastic deformation. Current trends in coronary stent design are towards thinner struts, potentially increasing the influence of SSEs. In order to maintain adequate device performance with decreasing strut thickness, it is therefore important to assess the role of SSEs in the plastic deformation of stents. In this study, finite element modelling and crystal plasticity theory are used to investigate SSEs in the deformation of struts in tension and bending. The relationships between SSEs and microstructure morphology, alloy strain hardening behaviour and secondary phases are also investigated. It is predicted that reducing the number of grains through the strut cross section and increasing the number of grains along the strut length have detrimental effects on mechanical performance. The magnitudes of these effects are predicted to be independent of the uniformity of the studied microstructures, but dependent on alloy strain hardening behaviour. It is believed that model predictions will aid in identifying a lower bound on suitable strut thicknesses in coronary stents for a range of alloys and microstructures.

Micromotion and friction evaluation of a novel surface architecture for improved primary fixation of cementless orthopaedic implants.

A new surface architecture (OsteoAnchor) for orthopaedic stem components has been developed, which incorporates a multitude of tiny anchor features for embedding into the bone during implantation. It was tested for its ability to provide improved primary fixation compared to existing surface coatings. Friction testing was performed on bovine trabecular bone. It was found that OsteoAnchor provided up to 76% greater resistance to transverse motion under simultaneous normal loading compared to the porous tantalum. Micromotion testing was performed on stem components implanted in cadaver ovine femurs. The micromotion amplitudes for the OsteoAnchor stem were significantly lower than for a corresponding plasma sprayed stem. These results demonstrate that OsteoAnchor has the potential to provide improved primary fixation for stem components in joint replacement operations.

Computational modelling of the mechanics of trabecular bone and marrow using fluid structure interaction techniques.

Bone marrow found within the porous structure of trabecular bone provides a specialized environment for numerous cell types, including mesenchymal stem cells (MSCs). Studies have sought to characterize the mechanical environment imposed on MSCs, however, a particular challenge is that marrow displays the characteristics of a fluid, while surrounded by bone that is subject to deformation, and previous experimental and computational studies have been unable to fully capture the resulting complex mechanical environment. The objective of this study was to develop a fluid structure interaction (FSI) model of trabecular bone and marrow to predict the mechanical environment of MSCs in vivo and to examine how this environment changes during osteoporosis. An idealized repeating unit was used to compare FSI techniques to a computational fluid dynamics only approach. These techniques were used to determine the effect of lower bone mass and different marrow viscosities, representative of osteoporosis, on the shear stress generated within bone marrow. Results report that shear stresses generated within bone marrow under physiological loading conditions are within the range known to stimulate a mechanobiological response in MSCs in vitro. Additionally, lower bone mass leads to an increase in the shear stress generated within the marrow, while a decrease in bone marrow viscosity reduces this generated shear stress.