Passive control of a bi-ventricular assist device : an experimental and numerical investigation

For the last two decades heart disease has been the highest single cause of death for the human population. With an alarming number of patients requiring heart transplant, and donations not able to satisfy the demand, treatment looks to mechanical alternatives. Rotary Ventricular Assist Devices, VADs, are miniature pumps which can be implanted alongside the heart to assist its pumping function. These constant flow devices are smaller, more efficient and promise a longer operational life than more traditional pulsatile VADs. The development of rotary VADs has focused on single pumps assisting the left ventricle only to supply blood for the body. In many patients however, failure of both ventricles demands that an additional pulsatile device be used to support the failing right ventricle. This condition renders them hospital bound while they wait for an unlikely heart donation. Reported attempts to use two rotary pumps to support both ventricles concurrently have warned of inherent haemodynamic instability. Poor balancing of the pumps’ flow rates quickly leads to vascular congestion increasing the risk of oedema and ventricular ‘suckdown’ occluding the inlet to the pump. This thesis introduces a novel Bi-Ventricular Assist Device (BiVAD) configuration where the pump outputs are passively balanced by vascular pressure. The BiVAD consists of two rotary pumps straddling the mechanical passive controller. Fluctuations in vascular pressure induce small deflections within both pumps adjusting their outputs allowing them to maintain arterial pressure. To optimise the passive controller’s interaction with the circulation, the controller’s dynamic response is optimised with a spring, mass, damper arrangement. This two part study presents a comprehensive assessment of the prototype’s ‘viability’ as a support device. Its ‘viability’ was considered based on its sensitivity to pathogenic haemodynamics and the ability of the passive response to maintain healthy circulation. The first part of the study is an experimental investigation where a prototype device was designed and built, and then tested in a pulsatile mock circulation loop. The BiVAD was subjected to a range of haemodynamic imbalances as well as a dynamic analysis to assess the functionality of the mechanical damper. The second part introduces the development of a numerical program to simulate human circulation supported by the passively controlled BiVAD. Both investigations showed that the prototype was able to mimic the native baroreceptor response. Simulating hypertension, poor flow balancing and subsequent ventricular failure during BiVAD support allowed the passive controller’s response to be assessed. Triggered by the resulting pressure imbalance, the controller responded by passively adjusting the VAD outputs in order to maintain healthy arterial pressures. This baroreceptor-like response demonstrated the inherent stability of the auto regulating BiVAD prototype. Simulating pulmonary hypertension in the more observable numerical model, however, revealed a serious issue with the passive response. The subsequent decrease in venous return into the left heart went unnoticed by the passive controller. Meanwhile the coupled nature of the passive response not only decreased RVAD output to reduce pulmonary arterial pressure, but it also increased LVAD output. Consequently, the LVAD increased fluid evacuation from the left ventricle, LV, and so actually accelerated the onset of LV collapse. It was concluded that despite the inherently stable baroreceptor-like response of the passive controller, its lack of sensitivity to venous return made it unviable in its present configuration. The study revealed a number of other important findings. Perhaps the most significant was that the reduced pulse experienced during constant flow support unbalanced the ratio of effective resistances of both vascular circuits. Even during steady rotary support therefore, the resulting ventricle volume imbalance increased the likelihood of suckdown. Additionally, mechanical damping of the passive controller’s response successfully filtered out pressure fluctuations from residual ventricular function. Finally, the importance of recognising inertial contributions to blood flow in the atria and ventricles in a numerical simulation were highlighted. This thesis documents the first attempt to create a fully auto regulated rotary cardiac assist device. Initial results encourage development of an inlet configuration sensitive to low flow such as collapsible inlet cannulae. Combining this with the existing baroreceptor-like response of the passive controller will render a highly stable passively controlled BiVAD configuration. The prototype controller’s passive interaction with the vasculature is a significant step towards a highly stable new generation of artificial heart.

Theoretical and numerical investigation of plasmon nanofocusing in metallic tapered rods and grooves

Effective focusing of electromagnetic (EM) energy to nanoscale regions is one of the major challenges in nano-photonics and plasmonics. The strong localization of the optical energy into regions much smaller than allowed by the diffraction limit, also called nanofocusing, offers promising applications in nano-sensor technology, nanofabrication, near-field optics or spectroscopy. One of the most promising solutions to the problem of efficient nanofocusing is related to surface plasmon propagation in metallic structures. Metallic tapered rods, commonly used as probes in near field microscopy and spectroscopy, are of a particular interest. They can provide very strong EM field enhancement at the tip due to surface plasmons (SP’s) propagating towards the tip of the tapered metal rod. A large number of studies have been devoted to the manufacturing process of tapered rods or tapered fibers coated by a metal film. On the other hand, structures such as metallic V-grooves or metal wedges can also provide strong electric field enhancements but manufacturing of these structures is still a challenge. It has been shown, however, that the attainable electric field enhancement at the apex in the V-groove is higher than at the tip of a metal tapered rod when the dissipation level in the metal is strong. Metallic V-grooves also have very promising characteristics as plasmonic waveguides. This thesis will present a thorough theoretical and numerical investigation of nanofocusing during plasmon propagation along a metal tapered rod and into a metallic V-groove. Optimal structural parameters including optimal taper angle, taper length and shape of the taper are determined in order to achieve maximum field enhancement factors at the tip of the nanofocusing structure. An analytical investigation of plasmon nanofocusing by metal tapered rods is carried out by means of the geometric optics approximation (GOA), which is also called adiabatic nanofocusing. However, GOA is applicable only for analysing tapered structures with small taper angles and without considering a terminating tip structure in order to neglect reflections. Rigorous numerical methods are employed for analysing non-adiabatic nanofocusing, by tapered rod and V-grooves with larger taper angles and with a rounded tip. These structures cannot be studied by analytical methods due to the presence of reflected waves from the taper section, the tip and also from (artificial) computational boundaries. A new method is introduced to combine the advantages of GOA and rigorous numerical methods in order to reduce significantly the use of computational resources and yet achieve accurate results for the analysis of large tapered structures, within reasonable calculation time. Detailed comparison between GOA and rigorous numerical methods will be carried out in order to find the critical taper angle of the tapered structures at which GOA is still applicable. It will be demonstrated that optimal taper angles, at which maximum field enhancements occur, coincide with the critical angles, at which GOA is still applicable. It will be shown that the applicability of GOA can be substantially expanded to include structures which could be analysed previously by numerical methods only. The influence of the rounded tip, the taper angle and the role of dissipation onto the plasmon field distribution along the tapered rod and near the tip will be analysed analytically and numerically in detail. It will be demonstrated that electric field enhancement factors of up to ~ 2500 within nanoscale regions are predicted. These are sufficient, for instance, to detect single molecules using surface enhanced Raman spectroscopy (SERS) with the tip of a tapered rod, an approach also known as tip enhanced Raman spectroscopy or TERS. The results obtained in this project will be important for applications for which strong local field enhancement factors are crucial for the performance of devices such as near field microscopes or spectroscopy. The optimal design of nanofocusing structures, at which the delivery of electromagnetic energy to the nanometer region is most efficient, will lead to new applications in near field sensors, near field measuring technology, or generation of nanometer sized energy sources. This includes: applications in tip enhanced Raman spectroscopy (TERS); manipulation of nanoparticles and molecules; efficient coupling of optical energy into and out of plasmonic circuits; second harmonic generation in non-linear optics; or delivery of energy to quantum dots, for instance, for quantum computations.

Atomistic numerical investigation of single-crystal copper nanowire with surface defects

Based on the embedded atom method (EAM), a molecular dynamics (MD) simulation is performed to study the single-crystal copper nanowire with surface defects through tension. The tension simulations for nanowire without defect are first carried out under different temperatures, strain rates and time steps and then surface defect effects for nanowire are investigated. The stress-strain curves obtained by the MD simulations of various strain rates show a rate below 1 x 10(9) s-1 will exert less effect on the yield strength and yield point, and the Young's modulus is independent of strain rate. a time step below 5 fs is recommend for the atomic model during the MD simulation. It is observed that high temperature leads to low Young's modulus, as well as the yield strength. The surface defects on nanowires are systematically studied in considering different defect orientations. It is found that the surface defect serves as a dislocation source, and the yield strength shows 34.20% decresse with 45 degree surface defect. Both yield strength and yield point are significantly influenced by the surface defects, except the Young's modulus.

Theoretical and numerical investigation of bending properties of Cu nanowires

Based on the AFM-bending experiments, a molecular dynamics (MD) bending simulation model is established which could accurately account for the full spectrum of the mechanical properties of NWs in a double clamped beam configuration, ranging from elasticity to plasticity and failure. It is found that, loading rate exerts significant influence to the mechanical behaviours of nanowires (NWs). Specifically, a loading rate lower than 10 m/s is found reasonable for a homogonous bending deformation. Both loading rate and potential between the tip and the NW are found to play an important role in the adhesive phenomenon. The force versus displacement (F-d) curve from MD simulation is highly consistent in shapes with that from experiments. Symmetrical F-d curves during loading and unloading processes are observed, which reveal the linear-elastic and non-elastic bending deformation of NWs. The typical bending induced tensile-compressive features are observed. Meanwhile, the simulation results are excellently fitted by the classical Euler-Bernoulli beam theory with axial effect. It is concluded that, axial tensile force becomes crucial in bending deformation when the beam size is down to nanoscale for double clamped NWs. In addition, we find shorter NWs will have an earlier yielding and a larger yielding force. Mechanical properties (Young’s modulus & yield strength) obtained from both bending and tensile deformations are found comparable with each other. Specifically, the modulus is essentially similar under these two loading methods, while the yield strength during bending is observed larger than that during tension.

Numerical investigation of mechanical properties of nanowires : a review

Nanowires (NWs) have attracted intensive researches owing to the broad applications that arise from their remarkable properties. Over the last decade, immense numerical studies have been conducted for the numerical investigation of mechanical properties of NWs. Among these numerical simulations, the molecular dynamics (MD) plays a key role. Herein we present a brief review on the current state of the MD investigation of nanowires. Emphasis will be placed on the FCC metal NWs, especially the Cu NWs. MD investigations of perfect NWs’ mechanical properties under different deformation conditions including tension, compression, torsion and bending are firstly revisited. Following in succession, the studies for defected NWs including the defects of twin boundaries (TBs) and pre-existing defects are discussed. The different deformation mechanism incurred by the presentation of defects is explored and discussed. This review reveals that the numerical simulation is an important tool to investigate the properties of NWs. However, the substantial gaps between the experimental measurements and MD results suggest the urgent need of multi-scale simulation technique.

Numerical investigation of the influence of topography on simulated downburst wind fields

A, dry, non-hydrostatic sub-cloud model is used to simulate an isolated stationary downburst wind event to study the influence topographic features have on the near-ground wind structure of these storms. It was generally found that storm maximum wind speeds could be increased by up to 30% because of the presence of a topographic feature at the location of maximum wind speeds. Comparing predicted velocity profile amplification with that of a steady flow impinging jet, similar results were found despite the simplifications made in the impinging jet model. Comparison of these amplification profiles with those found in the simulated boundary layer winds reveal reductions of up to 30% in the downburst cases. Downburst and boundary layer amplification profiles were shown to become more similar as the topographic feature height was reduced with respect to the outflow depth.

Numerical investigation of mechanical and thermal dynamic properties of the industrial transformer

Industrial transformer is one of the most critical assets in the power and heavy industry. Failures of transformers can cause enormous losses. The poor joints of the electrical circuit on transformers can cause overheating and results in stress concentration on the structure which is the major cause of catastrophic failure. Few researches have been focused on the mechanical properties of industrial transformers under overheating thermal conditions. In this paper, both mechanical and thermal properties of industrial transformers are jointly investigated using Finite Element Analysis (FEA). Dynamic response analysis is conducted on a modified transformer FEA model, and the computational results are compared with experimental results from literature to validate this simulation model. Based on the FEA model, thermal stress is calculated under different temperature conditions. These analysis results can provide insights to the understanding of the failure of transformers due to overheating, therefore are significant to assess winding fault, especially to the manufacturing and maintenance of large transformers.

Numerical investigation on stable crack growth in plane stress

Large deformation finite element analysis has been carried out to investigate the stress-strain fields ahead of a growing crack for compact tension .a=W D 0:5/ and three-point bend .a=W D 0:1 and 0:5/ specimens under plane stress condition. The crack growth is controlled by the experimental J -integral resistance curves measured by Sun et al. The results indicate that the distributions of opening stress, equivalent stress and equivalent strain ahead of a growing crack are not sensitive to specimen geometry. For both stationary and growing cracks, similar distributions of opening stress and triaxiality can be found along the ligament. During stable crack growth, the crack-tip opening displacement (CTOD) resistance curve and the cohesive fracture energy in the fracture process zone are independent of specimen geometry and may be suitable criteria for characterizing stable crack growth in plane stress.

Numerical investigation into the buckling behavior of advanced grid stiffened composite cylindrical shell

Advanced grid stiffened composite cylindrical shell is widely adopted in advanced structures due to its exceptional mechanical properties. Buckling is a main failure mode of advanced grid stiffened structures in engineering, which calls for increasing attention. In this paper, the buckling response of advanced grid stiffened structure is investigated by three different means including equivalent stiffness model, finite element model and a hybrid model (H-model) that combines equivalent stiffness model with finite element model. Buckling experiment is carried out on an advanced grid stiffened structure to validate the efficiency of different modeling methods. Based on the comparison, the characteristics of different methods are independently evaluated. It is arguable that, by considering the defects of material, finite element model is a suitable numerical tool for the buckling analysis of advanced grid stiffened structures.

Numerical investigation of graphene and graphene-polymer nanocomposites

This thesis is a comprehensive and deep investigation on graphene and graphene-polymer nanocomposites. It explores the strong structure-property relationships in both graphene and graphene-based polymeric nanocomposites. A number of significant conclusions, including failure mechanism in graphene, interfacial load transfer and thermal transport mechanisms in graphene-polymer nanocomposites, have been drawn through both atomistic simulations and theoretical analysis. These results can provide direct guidelines for development of new graphene-based materials and devices.

Numerical investigation of deformation and failure mechanisms of osteopontin-hydroxyapatite interfaces and mineralized collagen fibril arrays

This thesis is a comprehensive study of deformation and failure mechanisms in bone at nano- and micro-scale levels. It explores the mechanical behaviour of osteopontin-hydroxyapatite interfaces and mineralized collagen fibril arrays, through atomistic molecular dynamics and finite element simulations. This thesis shows some main factors contributing to the excellent material properties of bone and provides some guidelines for development of new artificial biological materials and medical implants.

Numerical investigation of upstream fuel injection through porous media for scramjet engines via surrogate-assisted evolutionary algorithms

A multi-objective design optimization study has been conducted for upstream fuel injection through porous media applied to the first ramp of a two-dimensional scramjet intake. The optimization has been performed by coupling evolutionary algorithms assisted by surrogate modeling and computational fluid dynamics with respect to three design criteria, that is, the maximization of the absolute mixing quantity, total pressure saving, and fuel penetration. A distinct Pareto optimal front has been obtained, highlighting the counteracting behavior of the total pressure against the mixing efficiency and fuel penetration. The injector location and size have been identified as the key design parameters as a result of a sensitivity analysis, with negligible influence of the porous properties in the configurations and conditions considered in the present study. Flowfield visualization has revealed the underlying physics associated with the effects of these dominant parameters on the shock structure and intensity.

Numerical investigation of plant tissue porosity and its influence on cellular level shrinkage during drying

Dried plant food products are increasing in demand in the consumer market, leading to continuing research to develop better products and processing techniques. Plant materials are porous structures, which undergo large deformations during drying. For any given food material, porosity and other cellular parameters have a direct influence on the level of shrinkage and deformation characteristics during drying, which involve complex mechanisms. In order to better understand such mechanisms and their interrelationships, numerical modelling can be used as a tool. In contrast to conventional grid-based modelling techniques, it is considered that meshfree methods may have a higher potential for modelling large deformations of multiphase problem domains. This work uses a meshfree based microscale plant tissue drying model, which was recently developed by the authors. Here, the effects of porosity have been newly accounted for in the model with the objective of studying porosity development during drying and its influence on shrinkage at the cellular level. For simplicity, only open pores are modelled and in order to investigate the influence of different cellular parameters, both apple and grape tissues were used in the study. The simulation results indicated that the porosity negatively influences shrinkage during drying and the porosity decreases as the moisture content reduces (when open pores are considered). Also, there is a clear difference in the deformations of cells, tissues and pores, which is mainly influenced by the cell wall contraction effects during drying.

Experimental and numerical investigation of strain-rate dependent mechanical properties of single living cells

The objective of this project is to investigate the strain-rate dependent mechanical behaviour of single living cells using both experimental and numerical techniques. The results revealed that living cells behave as porohyperlastic materials and that both solid and fluid phases within the cells play important roles in their mechanical responses. The research reported in this thesis provides a better understanding of the mechanisms underlying the cellular responses to external mechanical loadings and of the process of mechanical signal transduction in living cells. It would help us to enhance knowledge of and insight into the role of mechanical forces in supporting tissue regeneration or degeneration.