Patient specific finite element analysis results in more accurate prediction of stent fractures: Application to percutaneous pulmonary valve implantation

Stent fracture is a recognised complication following device implantation. Magnetic resonance data from a patient who underwent percutaneous pulmonary valve implantation (PPVI) and had subsequent stent fractures was used to create a finite element (FE) model of the patient's implantation site. Simulated expansion of the PPVI stent into this right ventricular outflow tract (RVOT) geometry was compared with free expansions of the PPVI stent up to a uniformly deployed configuration (conventional method employed in bench testing protocols), using FE analysis. PPVI biplane fluoroscopy images from the same patient were used to reconstruct the 3D shape and deformation of the stent in-situ and verify the FE geometrical results. Asymmetries were measured in all 3 orthogonal directions, in early systole and diastole.Although a simplified FE modelling of stent/implantation site interaction was adopted, this analysis gave useful information about the influence of the RVOT on the final geometry and mechanical performance of the stent. When deployed into the RVOT, the FE stent showed a non-uniform shape, similar to the geometry seen in the “real” fluoroscopy reconstructed stent, where the most expanded cells corresponded to the fracture locations. This asymmetrical geometry, when compared to the free-expanded stent, resulted in higher stresses in the portion of the stent where fractures occurred. Furthermore, fatigue fractures that were not predicted in the free-deployed stents, developed in the asymmetrically expanded device.In conclusion, the interaction between the PPVI device and the patient's RVOT is likely to be the crucial factor involved with this undesired event.     

Interim use of Jarvik-7 and Novacor artificial heart: blood rheology and transient ischemic attacks (TIA's)

The rheological properties of blood were studied in patients supported by both the Jarvik-7 total artificial heart (TAH) and Novacor left ventricular assist device (LVAD) as a bridge to cardiac transplantation. Both groups of patients had abnormalities in blood rheology which differed according to the type of device implanted as well as on the clinical state of the patient. The rheology of individual patients correlated well with their clinical status and outcome, with incidences of TIA's and/or stroke being accompanied by marked increases in relative blood viscosity, erythrocyte rigidity, fibrinogen concentration and platelet aggregation in varying combination. Observed abnormalities in blood rheology were also crucial to thrombus formation on artificial heart valves as well. Our results show that the therapeutic management of rheological parameters should prove to be a unique and clinically rewarding approach to these patients.     

Characterisation and simulation of an active microvalve for glaucoma

Glaucoma drainage device (GDD) has the potential to eliminate hypotony but still suffers from poor flow control and fibrosis. The ideal shunt should change its hydraulic resistance to achieve the desired intraocular pressure (IOP). In this study, the characterisation of a preliminary design of a new GDD is presented. This is activated by means of a diaphragm, which is actuated by conducting polymers. The valve can be manufactured employing microelectromechanical system technology by soft lithography. The characterisation process is performed by numerical simulation using the finite element method, considering the coupling between the fluid and the structure (diaphragm) obtaining the hydraulic resistance for several positions of the diaphragm. To analyse the hydraulic system of the microvalve implanted in a human eye, an equivalent circuit model was used. The parameters of the equivalent circuit model were obtained from numerical simulation. The hydraulic resistance of the designed GDD varies in the range of 13.08–0.36 mmHg min/μl compared with 3.38–0.43 mmHg min/μl for the Ahmed valve. The maximum displacement of the diaphragm in the vertical direction is 18.9 μm, and the strain in the plane is 2%. The proposed preliminary design allows to control the IOP by varying the hydraulic resistance in a greater range than the existing passive valves, and the numerical simulation facilitates the characterisation and the improvement of the design before its construction, reducing time and costs.     

Characterisation and simulation of an active microvalve for glaucoma

Glaucoma drainage device (GDD) has the potential to eliminate hypotony but still suffers from poor flow control and fibrosis. The ideal shunt should change its hydraulic resistance to achieve the desired intraocular pressure (IOP). In this study, the characterisation of a preliminary design of a new GDD is presented. This is activated by means of a diaphragm, which is actuated by conducting polymers. The valve can be manufactured employing microelectromechanical system technology by soft lithography. The characterisation process is performed by numerical simulation using the finite element method, considering the coupling between the fluid and the structure (diaphragm) obtaining the hydraulic resistance for several positions of the diaphragm. To analyse the hydraulic system of the microvalve implanted in a human eye, an equivalent circuit model was used. The parameters of the equivalent circuit model were obtained from numerical simulation. The hydraulic resistance of the designed GDD varies in the range of 13.08-0.36?mmHg?min/μl compared with 3.38-0.43?mmHg?min/μl for the Ahmed valve. The maximum displacement of the diaphragm in the vertical direction is 18.9?μm, and the strain in the plane is 2%. The proposed preliminary design allows to control the IOP by varying the hydraulic resistance in a greater range than the existing passive valves, and the numerical simulation facilitates the characterisation and the improvement of the design before its construction, reducing time and costs.     

Effect of Orientation of Fibers and Holes on the Radial Strain Amplification of Campaniform sensilla

##正## In this paper,a structural analysis is performed to gain insights on the synergistic mechanical amplification effect thatCampaniform sensilla have when combined in an array configuration.In order to simplify the analysis performed in this preliminaryinvestigation,an array of four holes in a single orthotropic lamina is considered.Firstly,a Finite Element Method(FEM) analysis is performed to discretely assess the influence that different geometrical parameters have on the mechanicalamplification properties of the array.Secondly,an artificial neural network is used to obtain an approximated multi-dimensionalcontinuous function,which models the relationship between the geometrical parameters and the amplification properties of thearray.Thirdly,an optimization is performed to identify the geometrical parameters yielding the maximum mechanical amplification.Finally,results are validated with an additional FEM simulation performed by varying geometrical parameters in theneighborhood of the identified optimal parameters.The method proposed in this paper can be fully automated and used to solvea wide range of optimization problems aimed at identifying optimal configurations of strain sensors inspired by Campaniformsensilla.     

A mechanical analysis of the closed Hancock heart valve prosthesis

In order to obtain mechanical specifications for the design of an artificial leaflet valve prosthesis, a geometrically non-linear numerical model is developed of a closed Hancock leaflet valve prosthesis. In this model, the fibre reinforcement of the leaflet and the viscoelastic properties of frame and leaflets are incorporated. The calculations are primarily restricted to 1/6 part of the valve and a time varying pressure load is applied. The calculations are verified experimentally by measuring the commissure displacements and leaflet centre displacement of a Hancock valve. The numerically obtained commissure displacements are found to be linearly dependent on the pressure load, and the slope of the curves is hardly dependent on loading type and loading velocity. Experimentally a difference is found between the three commissure displacements, which is also predicted numerically using a simplified asymmetric total valve model. Besides, experimentally a clear dependency of commissure displacements on frame size is found. For the leaflet centre displacement, a qualitative agreement exists between numerical prediction and experimental result, although the numerical predicted values are systematically higher. The numerically obtained stress distributions revealed that the maximum von Mises intensity in the membranes occurs in the vicinity of the commissure in the free leaflet area (0.2 N mm-2). Wrinkling of the membranes may occur in the coaptation area near the leaflet suspension. The maximum fibre stress is found near the aortic ring in the fibres which form the boundaries of the coaptation area (0.64 N mm-2). These locations seem to correlate with some common regions of tissue valve failure.     

Population-specific material properties of the implantation site for transcatheter aortic valve replacement finite element simulations

Patient-specific computational models are an established tool to support device development and test under clinically relevant boundary conditions. Potentially, such models could be used to aid the clinical decision-making process for percutaneous valve selection; however, their adoption in clinical practice is still limited to individual cases. To be fully informative, they should include patient-specific data on both anatomy and mechanics of the implantation site. In this work, fourteen patient-specific computational models for transcatheter aortic valve replacement (TAVR) with balloon-expandable Sapien XT devices were retrospectively developed to tune the material parameters of the implantation site mechanical model for the average TAVR population.Pre-procedural computed tomography (CT) images were post-processed to create the 3D patient-specific anatomy of the implantation site. Balloon valvuloplasty and device deployment were simulated with finite element (FE) analysis. Valve leaflets and aortic root were modelled as linear elastic materials, while calcification as elastoplastic. Material properties were initially selected from literature; then, a statistical analysis was designed to investigate the effect of each implantation site material parameter on the implanted stent diameter and thus identify the combination of material parameters for TAVR patients.These numerical models were validated against clinical data. The comparison between stent diameters measured from post-procedural fluoroscopy images and final computational results showed a mean difference of 2.5 ± 3.9%. Moreover, the numerical model detected the presence of paravalvular leakage (PVL) in 79% of cases, as assessed by post-TAVR echocardiographic examination.The final aim was to increase accuracy and reliability of such computational tools for prospective clinical applications.     

Fundamental mechanics of aortic heart valve closure

Stresses in a prosthetic heart valve at closure are determined by its geometrical and structural characteristics, by the mechanical support environment, and by the momentum of the valve leaflets or occluder and of the blood at the instant of closure. The mass of blood to be arrested is significantly greater than that of the leaflets or occluder, and is therefore likely to dominate the closure impulse. The kinetic energy of the blood must be transduced into potential energy in the structural components (valve leaflets, aortic root and aorta). This paper presents a methodology for computation and parameterisation of the blood momentum associated with a valve in the aortic position. It is suggested that the influence of physiological parameters, such as systolic waveform and systemic impedance, on the closure characteristics can be investigated based on the fluid dynamic implications. Detailed results are presented for a single leaflet mechanical valve (Bjork-Shiley 60 degrees Convexo-Concave). It is demonstrated that a simple analytical method can yield results that might be adequate for the purposes of valve design.     

The diastolic gradient of pressure and the effective area of section artifical mitral valve in aspect of its dysfunction

In article seven-year experience of treatment 126 sick mitral heart diseases to which have been implanted domestic monofolding the MIKS and folding MEDING-2 and ROSCARDIKS artificial mitral valve. At a comparative estimation of a haemodynamic efficiency and the analysis of frequency of occurrence of dysfunctions of the specified mechanical artificial valve it is revealed, that till three years the postimlantsperiod implant the MIKS and MEDING-2 possess advantage over ROSCARDIKS on haemodynamic properties, despite priority ROSCARDIKS on a standard size. It is shown, that initially low diastolic pressure gradient on mitral artificial valve and initially big area effective apertures mitral artificial valve have crucial importance in aspect of preventive maintenance of formation valve complications and reduction of number of repeated operations by open heart.     

Mechanics of the urethral duct: tissue constitutive formulation and structural modeling for the investigation of lumen occlusion

Urinary incontinence, often related to sphincter damage, is found in male patients, leading to a miserable quality of life and to huge costs for the healthcare system. The most effective surgical solution currently considered for men is the artificial urinary sphincter that exerts a pressure field on the urethra, occluding the duct. The evaluation of this device is currently based on clinical and surgical competences. The artificial sphincter design and mechanical action can be investigated by a biomechanical model of the urethra under occlusion, evaluating the interaction between tissues and prosthesis. A specific computational approach to urethral mechanics is here proposed, recalling the results of previous biomechanical experimental investigation. In this preliminary analysis, the horse urethra is considered, in the light of a significant correlation with human and in consideration of the relevant difficulty to get to human samples, which anyway represents the future advance. Histological data processing allow for the definition of a virtual and a subsequent finite element model of a urethral section. A specific hyperelastic formulation is developed to characterize the nonlinear mechanical behavior. The inverse analysis of tensile tests on urethra samples leads to the definition of preliminary constitutive parameters. The parameters are further refined by the computational analysis of inflation tests carried out on the entire urethral structure. The results obtained represent, in the light of the correlation reported, a valid preliminary support for the information to be assumed for prosthesis design, integrating surgical and biomechanical competences.     

Ventriculo-aortic junction in human root. A geometric approach

With advances in tissue engineering and improvement of surgical techniques, stentless biological valves and valve-sparing procedures have become alternatives to traditional aortic valve replacement with stented bioprostheses or mechanical valves. New surgical techniques preserve the advantages of native valves but require better understanding of the anatomical structure of the aortic root. Silicone rubber was injected in fresh aortic roots of nine human cadavers under the physiological closing pressure of 80 mmHg. The casts reproduced every detail of the aortic root anatomy and were used to digitize 27 leaflet attachment lines (LALs) of the aortic valves. LALs were normalized and described with a mathematical model. LALs were found to follow a pattern with the right coronary being the largest followed by the non-coronary and then the left coronary. During diastole, the aortic valve LAL can be described by an intersection between a created tube and an extruded parabolic surface. This geometrical definition of the LAL during end diastole gives a better understanding of the aortic root anatomy and could be useful for heart valve design and improvement of aortic valve reconstruction technique.     

A theoretical method to improve and optimize the design of bioartificial livers

Bioartificial livers (BALs) are a potentially effective countermeasure against liver failure, particularly in cases of acute or fulminant liver failure. It is hoped these devices can sustain a patient's liver function until recovery or transplant. However, no large‐scale clinical trial has yet proven that BALs are particularly effective and evidently design issues remain to be addressed. One aspect of BAL design that must be considered is the mass transfer of adequate oxygen to the hepatocytes within the device. We present here a mathematical modeling approach to oxygen mass transport in a BAL. A mathematical model based upon Krogh cylinders is outlined to describe a diffusion‐limited hollow fiber bioreactor. In addition, operating constraints are defined on the system—cells should not experience hypoxia and the cell population should be of adequate size. By combining modeling results with these operating constraints and presenting the results graphically, “operating region” charts can be constructed for the hollow fiber BAL (HF‐BAL). The effects of varying various operating parameters on the BAL are then established. It is found that smaller radii and short, thin walled fibers are generally advantageous while cell populations in excess of 10 billion could be supported in the BAL with a plasma flow rate of 200 mL/min. For fibers of intermediate length and lumen radius, the minimum number of fibers required to produce a viable design ranges approximately from 7,000–10,000. In theory, this may be enough to support patients with failing livers. Biotechnol. Bioeng. 2010;106: 980–988. © 2010 Wiley Periodicals, Inc.     

Geometrical stress-reducing factors in the anisotropic porcine heart valves

This study carries out a detailed parameter study based on a nonlinear anisotropic finite-element model published previously. The aim of this study is to identify the stress-reducing influences from geometrical parameters such as stent height, valve diameter, and the nonuniform thickness of porcine aortic valves under static loading condition. The anisotropy of the valve is considered to be transversely isotropic with fibers oriented along the circumferential directions, which enables us to use a simple anisotropic constitutive model using uniaxial experimental data. The results showed that in general, higher stent height and smaller diameter combined with nonuniform thickness give rise to a much more reduced overall stress level. Although the absolute values of the peak stresses may be influenced by the detailed orientations of fibers, the trends of the stress variation with the geometrical factors seem to be qualitatively consistent within the parameter ranges considered.     

The influence of the leaflets' curvature on the flow field in two bileaflet prosthetic heart valves

A successful mechanical prosthetic heart valve design is the bileaflet valve, which has been implanted for the first time more than 20 years ago. A key feature of bileaflet valves is the geometry of the two leaflets, which can be very important in determining the flow field. Laser Doppler anemometry (LDA) was used to perform an accurate study of the velocity and turbulence shear stress peak values (TSS(max)) fields at four distances from the valve plane. TSS(max) is a relevant parameter to assess the risk of hemolysis and platelet activation associated to the implantation of a prosthetic device, continuously interacting with blood. Two bileaflet valves were tested: the St. Jude HP and the Sorin Bicarbon, of the same nominal size (19mm). The former has flat leaflets, whereas the latter's leaflets have a cylindrical surface. A high regime (CO: 6l/min) was imposed, in order to test the two valves at maximum Reynolds number and consequent turbulence generation. The flat-leaflet design of the St. Jude generates a TSS field constant with distance; on the contrary, the Bicarbon's shear stress field undergoes an evident development, with an unexpected central peak at a distance comparable to the valve's dimensions (21mm). The two bileaflet valves tested, although very similar in design, behave very differently as for their turbulence properties. In particular, the concept of curved wake leads to conclude that the curvature of the leaflets' surface must be identified as an important parameter, which deserves careful attention in PHV design and development.     

A thin film nitinol heart valve

In order to create a less thrombogenic heart valve with improved longevity, a prosthetic heart valve was developed using thin film nitinol (NiTi). A "butterfly" valve was constructed using a single, elliptical piece of thin film NiTi and a scaffold made from Teflon tubing and NiTi wire. Flow tests and pressure readings across the valve were performed in vitro in a pulsatile flow loop. Bio-corrosion experiments were conducted on untreated and passivated thin film nitinol. To determine the material's in vivo biocompatibility, thin film nitinol was implanted in pigs using stents covered with thin film NiTi. Flow rates and pressure tracings across the valve were comparable to those through a commercially available 19 mm Perimount Edwards tissue valve. No signs of corrosion were present on thin film nitinol samples after immersion in Hank's solution for one month. Finally, organ and tissue samples explanted from four pigs at 2, 3, 4, and 6 weeks after thin film NiTi implantation appeared without disease, and the thin film nitinol itself was without thrombus formation. Although long term testing is still necessary, thin film NiTi may be very well suited for use in artificial heart valves.     

Lipid Bilayer Vesicle Generation Using Microfluidic Jetting

Bottom-up synthetic biology presents a novel approach for investigating and reconstituting biochemical systems and, potentially, minimal organisms. This emerging field engages engineers, chemists, biologists, and physicists to design and assemble basic biological components into complex, functioning systems from the bottom up. Such bottom-up systems could lead to the development of artificial cells for fundamental biological inquiries and innovative therapies^1,2. Giant unilamellar vesicles (GUVs) can serve as a model platform for synthetic biology due to their cell-like membrane structure and size. Microfluidic jetting, or microjetting, is a technique that allows for the generation of GUVs with controlled size, membrane composition, transmembrane protein incorporation, and encapsulation³. The basic principle of this method is the use of multiple, high-frequency fluid pulses generated by a piezo-actuated inkjet device to deform a suspended lipid bilayer into a GUV. The process is akin to blowing soap bubbles from a soap film. By varying the composition of the jetted solution, the composition of the encompassing solution, and/or the components included in the bilayer, researchers can apply this technique to create customized vesicles. This paper describes the procedure to generate simple vesicles from a droplet interface bilayer by microjetting.     

Dynamic modelling of prosthetic chorded mitral valves using the immersed boundary method

Current artificial heart valves either have limited lifespan or require the recipient to be on permanent anticoagulation therapy. In this paper, effort is made to assess a newly developed bileaflet valve prosthesis made of synthetic flexible leaflet materials, whose geometry and material properties are based on those of the native mitral valve, with a view to providing superior options for mitral valve replacement. Computational analysis is employed to evaluate the geometric and material design of the valve, by investigation of its mechanical behaviour and unsteady flow characteristics. The immersed boundary (IB) method is used for the dynamic modelling of the large deformation of the valve leaflets and the fluid-structure interactions. The IB simulation is first validated for the aortic prosthesis subjected to a hydrostatic loading. The predicted displacement fields by IB are compared with those obtained using ANSYS, as well as with experimental measurements. Good quantitative agreement is obtained. Moreover, known failure regions of aortic prostheses are identified. The dynamic behaviour of the valve designs is then simulated under four physiological pulsatile flows. Experimental pressure gradients for opening and closure of the valves are in good agreement with IB predictions for all flow rates for both aortic and mitral designs. Importantly, the simulations predicted improved physiological haemodynamics for the novel mitral design. Limitation of the current IB model is also discussed. We conclude that the IB model can be developed to be an extremely effective dynamic simulation tool to aid prosthesis design.     

Model reduction methodology for computational simulations of endovascular repair

Endovascular aneurysm repair (EVAR) is a current alternative treatment for thoracic and abdominal aortic aneurysms, but is still sometimes compromised by possible complications such as device migration or endoleaks. In order to assist clinicians in preventing these complications, finite element analysis (FEA) is a promising tool. However, the strong material and geometrical nonlinearities added to the complex multiple contacts result in costly finite-element models. To reduce this computational cost, we establish here an alternative and systematic methodology to simplify the computational simulations of stent-grafts (SG) based on FEA. The model reduction methodology relies on equivalent shell models with appropriate geometrical and mechanical parameters. It simplifies significantly the contact interactions but still shows very good agreement with a complete reference finite-element model. Finally, the computational time for EVAR simulations is reduced of a factor 6–10. An application is shown for the deployment of a SG during thoracic endovascular repair, showing that the developed methodology is both effective and accurate to determine the final position of the deployed SG inside the aneurysm.     

Numerical simulation of steady flow in a two-dimensional total artificial heart model.

In this paper, a numerical simulation of steady laminar and turbulent flow in a two-dimensional model for the total artificial heart is presented. A trileaflet polyurethane valve was simulated at the outflow orifice while the inflow orifice had a trileaflet or a flap valve. The finite analytic numerical method was employed to obtain solutions to the governing equations in the Cartesian coordinates. The closure for turbulence model was achieved by employing the k-epsilon-E model. The SIMPLER algorithm was used to solve the problem in primitive variables. The numerical solutions of the simulated model show that regions of relative stasis and trapped vortices were smaller within the ventricular chamber with the flap valve at the inflow orifice than that with the trileaflet valve. The predicted Reynolds stresses distal to the inflow valve within the ventricular chamber were also found to be smaller with the flap valve than with the trileaflet valve. These results also suggest a correlation between high turbulent stresses and the presence of thrombus in the vicinity of the valves in the total artificial hearts. The computed velocity vectors and turbulent stresses were comparable with previously reported in vitro measurements in artificial heart chambers. Analysis of the numerical solutions suggests that geometries similar to the flap valve (or a tilting disk valve) results in a better flow dynamics within the total artificial heart chamber compared to a trileaflet valve.     

Prediction of stenting related adverse events through patient-specific finite element modelling

Right ventricular outflow tract (RVOT) calcific obstruction is frequent after homograft conduit implantation to treat congenital heart disease. Stenting and percutaneous pulmonary valve implantation (PPVI) can relieve the obstruction and prolong the conduit lifespan, but require accurate pre-procedural evaluation to minimize the risk of coronary artery (CA) compression, stent fracture, conduit injury or arterial distortion.Herein, we test patient-specific finite element (FE) modeling as a tool to assess stenting feasibility and investigate clinically relevant risks associated to the percutaneous intervention.Three patients undergoing attempted PPVI due to calcific RVOT conduit failure were enrolled; the calcific RVOT, the aortic root and the proximal CA were segmented on CT scans for each patient. We numerically reproduced RVOT balloon angioplasty to test procedure feasibility and the subsequent RVOT pre-stenting expanding the stent through a balloon-in-balloon delivery system.Our FE framework predicted the occurrence of CA compression in the patient excluded from the real procedure. In the two patients undergoing RVOT stenting, numerical results were consistent with intraprocedural in-vivo fluoroscopic evidences. Furthermore, it quantified the stresses on the stent and on the relevant native structures, highlighting their marked dependence on the extent, shape and location of the calcific deposits. Stent deployment induced displacement and mechanical loading of the calcific deposits, also impacting on the adjacent anatomical structures.This novel workflow has the potential to tackle the analysis of complex RVOT clinical scenarios, pinpointing the procedure impact on the dysfunctional anatomy and elucidating potential periprocedural complications.