Pancreatic anastomosis training models: Current status and future directions.

Postoperative pancreatic fistula (POPF) is a major cause of morbidity and mortality after pancreatoduodenectomy (PD), and previous research has focused on patient-related risk factors and comparisons between anastomotic techniques. However, it is recognized that surgeon experience is an important factor in POPF outcomes, and that there is a significant learning curve for the pancreatic anastomosis. The aim of this study was to review the current literature on training models for the pancreatic anastomosis, and to explore areas for future research. It is concluded that research is needed to understand the mechanical properties of the human pancreas in an effort to develop a synthetic model that closely mimics its mechanical properties. Virtual reality (VR) is an attractive alternative to synthetic models for surgical training, and further work is needed to develop a VR pancreatic anastomosis training module that provides both high fidelity and haptic feedback.

Frequency and time dependent viscoelastic characterization of pediatric porcine brain tissue in compression.

Understanding the viscoelastic behavior of pediatric brain tissue is critical to interpret how external mechanical forces affect head injury in children. However, knowledge of the viscoelastic properties of pediatric brain tissue is limited, and this reduces the biofidelity of developed numeric simulations of the pediatric head in analysis of brain injury. Thus, it is essential to characterize the viscoelastic behavior of pediatric brain tissue in various loading conditions and to identify constitutive models. In this study, the pediatric porcine brain tissue was investigated in compression with determine the viscoelasticity under small and large strain, respectively. A range of frequencies between 0.1 and 40 Hz was applied to determine frequency-dependent viscoelastic behavior via dynamic mechanical analysis, while brain samples were divided into three strain rate groups of 0.01/s, 1/s and 10/s for compression up to 0.3 strain level and stress relaxation to obtain time-dependent viscoelastic properties. At frequencies above 20 Hz, the storage modulus did not increase, while the loss modulus increased continuously. With strain rate increasing from 0.01/s to 10/s, the mean stress at 0.1, 0.2 and 0.3 strain increased to approximate 6.8, 5.6 and 4.4 times, respectively. The brain compressive response was sensitive to strain rate and frequency. The characterization of brain tissue will be valuable for development of head protection systems and prediction of brain injury.

A geometry-based finite element tool for evaluating mitral valve biomechanics.

Mitral valve function depends on its complex geometry and tissue health, with alterations in shape and tissue response affecting the long-term restorarion of function. Previous computational frameworks for biomechanical assessment are mostly based on patient-specific geometries; however, these are not flexible enough to yield a variety of models and assess mitral closure for individually tuned morphological parameters or material property representations. This study details the finite element approach implemented in our previously developed toolbox to assess mitral valve biomechanics and showcases its flexibility through the generation and biomechanical evaluation of different models. A healthy valve geometry was generated and its computational predictions for biomechanics validated against data in the literature. Moreover, two mitral valve models including geometric alterations associated with disease were generated and analysed. The healthy mitral valve model yielded biomechanical predictions in terms of valve closure dynamics, leaflet stresses and papillary muscle and chordae forces comparable to previous computational and experimental studies. Mitral valve function was compromised in geometries representing disease, expressed by the presence of regurgitating areas, elevated stress on the leaflets and unbalanced subvalvular apparatus forces. This showcases the flexibility of the toolbox concerning the generation of a range of mitral valve models with varying geometric definitions and material properties and the evaluation of their biomechanics.

Development and experimental validation of a dynamic numerical model for human articular cartilage.

The purpose of this study was to create a preliminary set of experimentally validated Finite Element Analysis (FEA) models, in order to predict the dynamic mechanical behaviour of human articular cartilage (AC). Current models consider static loading with limited independent experimental validation, while the models for this study assess dynamic loading of AC, with direct comparison and validation to physical testing. Three different FEA models of AC were constructed, which considered both linear elastic and hyperelastic models; Neo-Hookean and Ogden. Models were validated using the data collected from compression testing of human femoral heads across 0-1.7 MPa (quasi-static tests and dynamic mechanical analysis). The linear elastic model was inadequate, with a 10-fold over prediction of the displacement dynamic amplitude. The Neo-Hookean model accurately predicted the dynamic amplitude but failed to predict the initial compression of the cartilage, with a 10 times overprediction. The Ogden model provided the best results, with both the initial compression lying within one standard deviation of that observed in the validation data set, and the dynamic amplitude of the same order of magnitude. In conclusion, this study has found that the fast dynamic response of human AC is best represented by a third order Ogden model.

Analysis of hydration and subchondral bone density on the viscoelastic properties of bovine articular cartilage.

BACKGROUND: Articular cartilage is known to be a viscoelastic material, however little research has explored the impact of cartilage water content and bone density on its viscoelasticity. This study aimed to isolate subchondral bone density and hydration of articular cartilage and analyse their effects on the viscoelastic properties of articular cartilage. METHODS: Dynamic mechanical analysis was used to test samples at frequencies of 1, 8, 12, 29, 49, 71, and 88 Hz. Synthetic bone material with densities of 663.7 kg/m3 and 156.8 kg/m3 were used to mimic the bone mineral density (BMD). Dehydration occurred in a stepwise manner at relative humidity (RH) levels of 100%, 30%, and 1%. These relative humidity levels led to water contents of approximately 76%, 8.5%, and ≈ 0% by mass, respectively. RESULTS: Samples from eight bovine femoral heads were tested under a sinusoidal load. Storage stiffness was lower on the lower substrate density. Storage stiffness, though, increased as cartilage samples were dehydrated from a water content of 76% to 8.5%; decreasing again as the water content was further reduced. Loss stiffness was lower on a lower density substrate and decreased as the water content decreased. CONCLUSIONS: In conclusions, a decrease in hydration decreases the loss stiffness, but a non-linear relationship between hydration and storage stiffness may exist. Additionally, higher BMD values led to greater storage and loss stiffnesses.

Numerical modelling of blood rheology and platelet activation through a stenosed left coronary artery bifurcation.

Coronary bifurcations are prone to atherosclerotic plaque growth, experiencing regions of reduced wall shear stress (WSS) and increased platelet adhesion. This study compares effects across different rheological approaches on hemodynamics, combined with a shear stress exposure history model of platelets within a stenosed porcine bifurcation. Simulations used both single/multiphase blood models to determine which approach best predicts phenomena associated with atherosclerosis and atherothrombosis. A novel Lagrangian platelet tracking model was used to evaluate residence time and shear history of platelets indicating likely regions of thrombus formation. Results show a decrease in area of regions with pathologically low time-averaged WSS with the use of multiphase models, particularly in a stenotic bifurcation. Significant non-Newtonian effects were observed due to low-shear and varying hematocrit levels found on the outer walls of the bifurcation and distal to the stenosis. Platelet residence time increased 11% in the stenosed artery, with exposure times to low-shear sufficient for red blood cell aggregation (>1.5 s). increasing the risk of thrombosis. This shows stenotic artery hemodynamics are inherently non-Newtonian and multiphase, with variations in hematocrit (0.163-0.617) and elevated vorticity distal to stenosis (+15%) impairing the function of the endothelium via reduced time-averaged WSS regions, rheological properties and platelet activation/adhesion.

Investigation of the Compressive Viscoelastic Properties of Brain Tissue Under Time and Frequency Dependent Loading Conditions.

The mechanical characterization of brain tissue has been generally analyzed in the frequency and time domain. It is crucial to understand the mechanics of the brain under realistic, dynamic conditions and convert it to enable mathematical modelling in a time domain. In this study, the compressive viscoelastic properties of brain tissue were investigated under time and frequency domains with the same physical conditions and the theory of viscoelasticity was applied to estimate the prediction of viscoelastic response in the time domain based on frequency-dependent mechanical moduli through Finite Element models. Storage and loss modulus were obtained from white and grey matter, of bovine brains, using dynamic mechanical analysis and time domain material functions were derived based on a Prony series representation. The material models were evaluated using brain testing data from stress relaxation and hysteresis in the time dependent analysis. The Finite Element models were able to represent the trend of viscoelastic characterization of brain tissue under both testing domains. The outcomes of this study contribute to a better understanding of brain tissue mechanical behaviour and demonstrate the feasibility of deriving time-domain viscoelastic parameters from frequency-dependent compressive data for biological tissue, as validated by comparing experimental tests with computational simulations.

A toolbox for generating scalable mitral valve morphometric models.

The mitral valve is a complex anatomical structure, whose shape is key to several traits of its function and disease, being crucial for the success of surgical repair and implantation of medical devices. The aim of this study was to develop a parametric, scalable, and clinically useful model of the mitral valve, enabling the biomechanical evaluation of mitral repair techniques through finite element simulations. MATLAB was used to parameterize the valve: the annular boundary was sampled from a porcine mitral valve mesh model and landmark points and relevant boundaries were selected for the parameterization of leaflets using polynomial fitting. Several geometric parameters describing the annulus, leaflet shape and papillary muscle position were implemented and used to scale the model according to patient dimensions. The developed model, available as a toolbox, allows for the generation of a population of models using patient-specific dimensions obtained from medical imaging or averaged dimensions evaluated from empirical equations based on the Golden Proportion. The average model developed using this framework accurately represents mitral valve shapes, associated with relative errors reaching less than 10% for annular and leaflet length dimensions, and less than 24% in comparison with clinical data. Moreover, model generation takes less than 5 min of computing time, and the toolbox can account for individual morphological variations and be employed to evaluate mitral valve biomechanics; following further development and validation, it will aid clinicians when choosing the best patient-specific clinical intervention and improve the design process of new medical devices.

The effect of injurious compression on the elastic, hyper-elastic and visco-elastic properties of porcine peripheral nerves.

The aim of this study was to characterise the viscoelastic and hyper-elastic properties of the ulnar nerve before and after compression has been induced, in order to aid the understanding of how the mechanical properties of nerves are altered during nerve compression, a contributing factor to cubital tunnel syndrome. Ulnar nerves were dissected from porcine legs and tensile tested to 10% strain. The Young's modulus and Yeoh hyper-elastic model were used to evaluate the materials elastic and hyper-elastic properties respectively. Dynamic mechanical analysis (DMA) was used to evaluate the viscoelastic properties over a range of frequencies between 0.5 Hz and 38 Hz. The nerves were then compressed to 40% for 60 s and the same tests were carried out after compression. The nerves were stiffer after compression, the mean Young's modulus before was 0.181 MPa and increased to 0.601 MPa after compression. The mean shear modulus calculated from the Yeoh hyper-elastic model was also higher after compression increasing from 5 kPa to 7 kPa. After compression, these properties had significantly increased (p < 0.05). The DMA results showed that the nerves exhibit frequency dependent viscoelastic behaviour across all tested frequencies. The median values of storage modulus before compression ranged between 0.605 and 0.757 MPa across the frequencies and after compression between 1.161 MPa and 1.381 MPa. There was a larger range of median values for loss modulus, before compression, median values ranged between 0.073 MPa and 0.216 MPa and after compression from 0.165 MPa to 0.410 MPa. There was a significant increase in both storage and loss modulus after compression (p < 0.05). The mechanical properties of the nerve change following compression, however the response to decompression in vivo requires further evaluation to determine whether the observed changes persist, which may have implications for clinical recovery after surgical decompression in entrapment neuropathy.

Computational fluid dynamics of the right atrium: Assessment of modelling criteria for the evaluation of dialysis catheters.

Central venous catheters are widely used in haemodialysis therapy, having to respect design requirements for appropriate performance. These are placed within the right atrium (RA); however, there is no prior computational study assessing different catheter designs while mimicking their native environment. Here, a computational fluid dynamics model of the RA, based on realistic geometry and transient physiological boundary conditions, was developed and validated. Symmetric, split and step catheter designs were virtually placed in the RA and their performance was evaluated by: assessing their interaction with the RA haemodynamic environment through prediction of flow vorticity and wall shear stress (WSS) magnitudes (1); and quantifying recirculation and tip shear stress (2). Haemodynamic predictions from our RA model showed good agreement with the literature. Catheter placement in the RA increased average vorticity, which could indicate alterations of normal blood flow, and altered WSS magnitudes and distribution, which could indicate changes in tissue mechanical properties. All designs had recirculation and elevated shear stress values, which can induce platelet activation and subsequently thrombosis. The symmetric design, however, had the lowest associated values (best performance), while step design catheters working in reverse mode were associated with worsened performance. Different tip placements also impacted on catheter performance. Our findings suggest that using a realistically anatomical RA model to study catheter performance and interaction with the haemodynamic environment is crucial, and that care needs to be given to correct tip placement within the RA for improved recirculation percentages and diminished shear stress values.

Dynamic mechanical characterization and viscoelastic modeling of bovine brain tissue.

Brain tissue is vulnerable and sensitive, predisposed to potential damage under various conditions of mechanical loading. Although its material properties have been investigated extensively, the frequency-dependent viscoelastic characterization is currently limited. Computational models can provide a non-invasive method by which to analyze brain injuries and predict the mechanical response of the tissue. The brain injuries are expected to be induced by dynamic loading, mostly in compression and measurement of dynamic viscoelastic properties are essential to improve the accuracy and variety of finite element simulations on brain tissue. Thus, the aim of this study was to investigate the compressive frequency-dependent properties of brain tissue and present a mathematical model in the frequency domain to capture the tissue behavior based on experimental results. Bovine brain specimens, obtained from four locations of corona radiata, corpus callosum, basal ganglia and cortex, were tested under compression using dynamic mechanical analysis over a range of frequencies between 0.5 and 35 Hz to characterize the regional and directional response of the tissue. The compressive dynamic properties of bovine brain tissue were heterogenous for regions but not sensitive to orientation showing frequency dependent statistical results, with viscoelastic properties increasing with frequency. The mean storage and loss modulus were found to be 12.41 kPa and 5.54 kPa, respectively. The material parameters were obtained using the linear viscoelastic model in the frequency domain and the numeric simulation can capture the compressive mechanical behavior of bovine brain tissue across a range of frequencies. The frequency-dependent viscoelastic characterization of brain tissue will improve the fidelity of the computational models of the head and provide essential information to the prediction and analysis of brain injuries in clinical treatments.

The role of strut chordae in mitral valve competence during annular dilation.

Strut chordae, on their own, are not typically thought to aid mitral valve competence. The aim of this study is to assess whether strut chordae aid mitral valve competence during acute annular dilation. Twelve porcine hearts were dissected and tested using an in vitro simulator, with the mitral annulus tested in either a 'normal' or a dilated configuration. The normal configuration included a diameter of 30 mm, a posterior leaflet 'radius' of 15 mm and a commissural corner 'radius' of 7.5 mm; the dilated annular template instead used dimensions of 50 mm, 25 mm and 12.5 mm, respectively. Each mitral valve underwent ten repeat tests with a target systolic pressure of 100 mmHg. No significant difference in the pressure was detected between the dilated and regular annuli for the mitral valves tested (95 ± 3 mmHg cf. 95 ± 2 mmHg). However, the volume of regurgitation for a dilated annulus was 28 ml greater than for a valve with a normal annulus. Following severing of strut chordae, there was a significant reduction in the systolic pressure withstood before regurgitation by mitral valves with dilated annuli (60 ± 29 mmHg cf. 95 ± 2 mmHg for normal annular dimensions; p < 0.05). In conclusion, strut chordae tendineae may play a role in aiding mitral valve competence during pathophysiology.

Mechanical testing of glutaraldehyde cross-linked mitral valves. Part one: In vitro mechanical behaviour.

The aim of this study was to perform an initial assessment, in vitro, of the feasibility of using a glutaraldehyde cross-linked porcine mitral valve to retain acute functionality, focusing on assessing mitral regurgitation. Six porcine hearts were tested using an in vitro simulator. Testing was repeated following cross-linking of mitral valves; where cross-linking was achieved by placing them in a glutaraldehyde solution. The simulator enabled systolic pressure on the ventricular side of the valve to be mimicked. Following testing, mitral valve leaflets underwent Scanning Electron Microscopy of the ventricular surface of both the anterior and posterior leaflets (1 cm2 samples). The peak pressure withstood by cross-linked valves was significantly lower than for untreated valves (108 mmHg cf. 128 mmHg for untreated valves; p < 0.05). The peak pressure was typically reached 0.5 s later than for the untreated valve. While both cross-linked and untreated valves exhibited endothelium denudation, the unfixed valve had less endothelial loss. Glutaraldehyde cross-linking of porcine mitral valves may be of potential value in assessing improved bioprosthetic mitral valve replacements. However, a more immobile valve exhibiting endothelial denudation (i.e. sclerosis) was a possible concerns identified following in vitro acute assessment.

Mechanical testing of glutaraldehyde cross-linked mitral valves. Part two: Elastic and viscoelastic properties of chordae tendineae.

The aim of this study was to assess whether the mechanical properties of mitral valve chordae tendineae are sensitive to being cross-linked under load. A total 64 chordae were extracted from eight porcine hearts. Two chordae (posterior basal) from each heart were subjected to uniaxial ramp testing and six chordae (two strut, two anterior basal and two posterior basal) were subjected to dynamic mechanical analysis over frequencies between 0.5 and 10 Hz. Chordae were either cross-linked in tension or cross-linked in the absence of loading. Chordae cross-linked under load transitioned from high to low extension at a lower strain than cross-linked unloaded chordae (0.07 cf. 0.22), with greater pre-transitional (30.8 MPa cf. 5.78 MPa) and post-transitional (139 MPa cf. 74.1 MPa) moduli. The mean storage modulus of anterior strut chordae ranged from 48 to 54 MPa for cross-linked unloaded chordae, as compared to 53-61 MPa cross-linked loaded chordae. The mean loss modulus of anterior strut chordae ranged from 2.3 to 2.9 MPa for cross-linked unloaded chordae, as compared to 3.8-4.8 MPa cross-linked loaded chordae. The elastic and viscoelastic properties of chordae following glutaraldehyde cross-linking are dependent on the inclusion/exclusion of loading during the cross-linking process; with loading increasing the magnitude of the material properties measured.

Impact of side-hole geometry on the performance of hemodialysis catheter tips: A computational fluid dynamics assessment.

Hemodialysis catheters are used to support blood filtration, yet there are multiple fundamentally different approaches to catheter tip design with no clear optimal solution. Side-holes have been shown to increase flow rates and decrease recirculation but have been associated with clotting/increased infection rates. This study investigates the impact of changing the shape, size and number of side-holes on a simple symmetric tip catheter by evaluating the velocity, shear stress and shear rate of inflowing blood. A platelet model is used to examine the residence time and shear history of inflowing platelets. The results show that side-holes improve the theoretical performance of the catheters, reducing the maximum velocity and shear stress occurring at the tip compared to non-side-hole catheters. Increasing the side-hole area improved performance up to a point, past which not all inflow through the hole was captured, and instead a small fraction slowly 'washed-out' through the remainder of the tip resulting in greater residence times and increasing the likelihood of platelet adhesion. An oval shaped hole presents a lower chance of external fibrin formation compared to a circular hole, although this would also be influenced by the catheter material surface topology which is dependent on the manufacturing process. Overall, whilst side-holes may be associated with increased clotting and infection, this can be reduced when side-hole geometry is correctly implemented though; a sufficient area for body diameter (minimising residence time) and utilising angle-cut, oval shaped holes (reducing shear stress and chances of fibrin formation partially occluding holes).

A technique for measuring the frictional torque of articular cartilage and replacement biomaterials.

Understanding the tribological behaviour of articular cartilage enables the development of effective replacement biomaterials. This study presents a technique for the investigation of the frictional torque of articular cartilage, for the assessment of replacement biomaterials. A calcium alginate hydrogel was used as the biomaterial for this study. Three different specimen types were examined to include articular cartilage, calcium alginate hydrogel, alone, and in combination with articular cartilage. An axial load, varying from 10 to 100 N, was applied to the specimen and the frictional torque measured whilst an indenter underwent axial rotation from 0° to 2° to 0° for 100 cycles. The resulting frictional torque magnitude was evaluated with a smooth curve fitting function. Linear regression identified a statistically significant relationship between torque magnitude and axial load (p < 0.05) for all specimen variations. From 10 to 100 N of applied load, mean torque magnitude ranged from 0.08±0.010 to 0.11±0.013 N m, 0.08±0.012 to 0.09±0.016 N m and 0.07±0.017 to 0.09±0.020 N m (mean ± standard deviation), for articular cartilage, calcium alginate separately and in combination with articular cartilage, respectively. This study has established a suitable frictional torque testing protocol for potential cartilage replacement biomaterials.

Assessment of surface roughness and blood rheology on local coronary haemodynamics: a multi-scale computational fluid dynamics study.

The surface roughness of the coronary artery is associated with the onset of atherosclerosis. The study applies, for the first time, the micro-scale variation of the artery surface to a 3D coronary model, investigating the impact on haemodynamic parameters which are indicators for atherosclerosis. The surface roughness of porcine coronary arteries have been detailed based on optical microscopy and implemented into a cylindrical section of coronary artery. Several approaches to rheology are compared to determine the benefits/limitations of both single and multiphase models for multi-scale geometry. Haemodynamic parameters averaged over the rough/smooth sections are similar; however, the rough surface experiences a much wider range, with maximum wall shear stress greater than 6 Pa compared to the approximately 3 Pa on the smooth segment. This suggests the smooth-walled assumption may neglect important near-wall haemodynamics. While rheological models lack sufficient definition to truly encompass the micro-scale effects occurring over the rough surface, single-phase models (Newtonian and non-Newtonian) provide numerically stable and comparable results to other coronary simulations. Multiphase models allow for phase interactions between plasma and red blood cells which is more suited to such multi-scale models. These models require additional physical laws to govern advection/aggregation of particulates in the near-wall region.

Biomechanical Assessment of Bicuspid Aortic Valve Phenotypes: A Fluid-Structure Interaction Modelling Approach.

PURPOSE: Bicuspid aortic valve (BAV) is a congenital heart malformation with phenotypic heterogeneity. There is no prior computational study that assesses the haemodynamic and valve mechanics associated with BAV type 2 against a healthy tricuspid aortic valve (TAV) and other BAV categories. METHODS: A proof-of-concept study incorporating three-dimensional fluid-structure interaction (FSI) models with idealised geometries (one TAV and six BAVs, namely type 0 with lateral and anterior-posterior orientations, type 1 with R-L, N-R and N-L leaflet fusion and type 2) has been developed. Transient physiological boundary conditions have been applied and simulations were run using an Arbitrary Lagrangian-Eulerian formulation. RESULTS: Our results showed the presence of abnormal haemodynamics in the aorta and abnormal valve mechanics: type 0 BAVs yielded the best haemodynamical and mechanical outcomes, but cusp stress distribution varied with valve orifice orientation, which can be linked to different cusp calcification location onset; type 1 BAVs gave rise to similar haemodynamics and valve mechanics, regardless of raphe position, but this position altered the location of abnormal haemodynamic features; finally, type 2 BAV constricted the majority of blood flow, exhibiting the most damaging haemodynamic and mechanical repercussions when compared to other BAV phenotypes. CONCLUSION: The findings of this proof-of-concept work suggest that there are specific differences across haemodynamics and valve mechanics associated with BAV phenotypes, which may be critical to subsequent processes associated with their pathophysiology processes.

A method for the assessment of the coefficient of friction of articular cartilage and a replacement biomaterial.

Replacement biomaterials for articular cartilage should encourage a coefficient of friction similar to the natural joint. Whilst the literature has assessed the coefficient of friction of articular cartilage against that of a potential biomaterial, it is unknown what the friction of articular cartilage in sliding against a surface defect, repaired with a biomaterial is. This evaluation is crucial to allow for the development of effective biomaterials to closely have the behaviour of articular cartilage. Thus, the aim of this study was to develop a novel technique for the assessment of the coefficient of friction of replacement biomaterials within articular cartilage, with this original testing configuration. For this study, a biomaterial was induced within an artificial defect perforated on the surface of bovine articular cartilage, whilst the material was assessed in sliding against articular cartilage itself. Calcium alginate was selected as the sample biomaterial for evaluation in this study. The tests were performed in sliding on a pin-on-disc tribometer in Ringer's solution. Two further tests were carried out, one as a benchmark comparison of a cartilage pin against a cartilage plate, as well as a cartilage pin against an aluminium plate. A constant induced stress of 0.06 MPa was applied at a frequency of 1 Hz. For the cartilage-cartilage, cartilage/hydrogel-cartilage and cartilage-aluminium test, the overall median coefficient of friction extracted across six repeats was of 0.36, 0.38 and 0.32, respectively. Statistical insignificance was identified across all three groups tested (p > 0.05). Similarity was observed in the coefficient of friction of cartilage-cartilage and cartilage/hydrogel-cartilage tests, however high-speed data identified the greatest wear for the cartilage/hydrogel-cartilage test.