Magnetic flux analysis in superconducting Fault Current Limiters under short-circuit conditions

In order to obtain a more compact Superconducting Fault Current limiter (SFCL), a special geometry of core and AC coil is required. This results in a unique magnetic flux pattern which differs from those associated with conventional round core arrangements. In this paper the magnetic flux density within a Fault Current Limiter (FCL) is described. Both experimental and analytical approaches are considered. A small scale prototype of an FCL was constructed in order to conduct the experiments. This prototype comprises a single phase. The analysis covers both the steady state and the short-circuit condition. Simulation results were obtained using commercial software based on the Finite Element Method (FEM). The magnetic flux saturating the cores, leakage magnetic flux giving rise to electromagnetic forces and leakage magnetic flux flowing in the enclosing tank are computed.

Patient-specific simulations of spinal deformity surgery

Adolescent idiopathic scoliosis (AIS) is a three-dimensional spinal deformity involving the side-to-side curvature of the spine in the coronal plane and axial rotation of the vertebrae in the transverse plane. For patients with a severe or rapidly progressing deformity, corrective instrumented fusion surgery is performed. The wide choice of implants and large variability between patients make it difficult for surgeons to choose optimal treatment strategies. This paper describes the patient specific finite element modelling techniques employed and the results of preliminary analyses predicting the surgical outcomes for a series of AIS patients. This report highlights the importance of not only patient-specific anatomy and material parameters, but also patient-specific data for the clinical and physiological loading conditions experienced by the patient who has corrective scoliosis surgery.

A coupled SPH/FEM analysis of portable water filled barriers

Portable water filled barriers (PWFB) are semi-rigid roadside barriers which have the potential to display good crash attenuation characteristics at low and moderate impact speeds. The traditional mesh based numerical methods alone fail to simulate this type of impact with precision, stability and efficiency. This paper proposes to develop an advanced simulation model based on the combination of Smoothed Particles Hydrodynamics (SPH), a meshless method, and finite element method (FEM) for fluid-structure analysis using the commercially available software package LS-Dyna. The interaction between SPH particles and FEA elements is studied in this paper. Two methods of element setup at the element boundary were investigated. The response of the impacted barrier and fluid inside were analysed and compared. The system response and lagging were observed and reported in this paper. It was demonstrated that coupled SPH/FEM can be used in full scale PWFB modelling application. This will aid the research in determining the best initial setup to couple FEA and SPH in road safety barrier for impact response and safety analysis in the future.

Saturated core high-temperature superconducting fault current limiters as an alternative to conventional series reactors in a distribution grid

Series reactors are used in distribution grids to reduce the short-circuit fault level. Some of the disadvantages of the application of these devices are the voltage drop produced across the reactor and the steep front rise of the transient recovery voltage (TRV), which generally exceeds the rating of the associated circuit breaker. Simulations were performed to compare the characteristics of a saturated core High-Temperature Superconducting Fault Current Limiter (HTS FCL) and a series reactor. The design of the HTS FCL was optimized using the evolutionary algorithm. The resulting Pareto frontier curve of optimum solution is presented in this paper. The results show that the steady-state impedance of an HTS FCL is significantly lower than that of a series reactor for the same level of fault current limiting. Tests performed on a prototype 11 kV HTS FCL confirm the theoretical results. The respective transient recovery voltages (TRV) of the HTS FCL and an air core reactor of comparable fault current limiting capability are also determined. The results show that the saturated core HTS FCL has a significantly lower effect on the rate of rise of the circuit breaker TRV as compared to the air core reactor. The simulations results are validated with shortcircuit test results.

Non-Newtonian natural convection flow along an isothermal horizontal circular cylinder using modified power-law model

Laminar two-dimensional natural convection boundary-layer flow of non-Newtonian fluids along an isothermal horizontal circular cylinder has been studied using a modified power-law viscosity model. In this model, there are no unrealistic limits of zero or infinite viscosity. Therefore, the boundary-layer equations can be solved numerically by using marching order implicit finite difference method with double sweep technique. Numerical results are presented for the case of shear-thinning as well as shear thickening fluids in terms of the fluid velocity and temperature distributions, shear stresses and rate of heat transfer in terms of the local skin-friction and local Nusselt number respectively.

Cohesive zone modelling of mineralized collagen fibril arrays in bending

Hard biological materials such as bone possess superior material properties of high stiffness and toughness. Two unique characteristics of bone microstructure are a large aspect ratio of mineralized collagen fibrils (MCF), and an extremely thin and large area of extrafibrillar protein matrix located between the MCF. The objective of this study is to investigate the effects of: (1) MCF aspect ratio, and (2) energy dissipation in extrafibrillar protein matrix on the mechanical behaviour of MCF arrays. In this study, notched specimens of MCF arrays in extrafibrillar protein matrix are subjected to bending. Cohesive zone model was implemented to simulate the failure of extrafibrillar protein matrix. The study reveals that the MCF array with a higher MCF aspect ratio and the MCF array with a higher protein energy dissipation in the interface direction are able to sustain a higher bending force and dissipate higher energy.

Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

Development of optimal designs of insulated rail joints

Proper functioning of Insulated Rail Joints (IRJs) is essential for the safe operation of the railway signalling systems and broken rail identification circuitries. The Conventional IRJ (CIRJ) resembles structural butt joints consisting of two pieces of rails connected together through two joint bars on either side of their web and the assembly is held together through pre-tensioned bolts. As the IRJs should maintain electrical insulation between the two rails, a gap between the rail ends must be retained at all times and all metal contacting surfaces should be electrically isolated from each other using non-conductive material. At the gap, the rail ends lose longitudinal continuity and hence the vertical sections of the rail ends are often severely damaged, especially at the railhead, due to the passage of wheels compared to other continuously welded rail sections. Fundamentally, the reason for the severe damage can be related to the singularities of the wheel-rail contact pressure and the railhead stress. No new generation designs that have emerged in the market to date have focussed on this fundamental; they only have provided attention to either the higher strength materials or the thickness of the sections of various components of the IRJs. In this thesis a novel method of shape optimisation of the railhead is developed to eliminate the pressure and stress singularities through changes to the original sharp corner shaped railhead into an arc profile in the longitudinal direction. The optimal shape of the longitudinal railhead profile has been determined using three nongradient methods in search of accuracy and efficiency: (1) Grid Search Method; (2) Genetic Algorithm Method and (3) Hybrid Genetic Algorithm Method. All these methods have been coupled with a parametric finite element formulation for the evaluation of the objective function for each iteration or generation depending on the search algorithm employed. The optimal shape derived from these optimisation methods is termed as Stress Minimised Railhead (SMRH) in this thesis. This optimal SMRH design has exhibited significantly reduced stress concentration that remains well below the yield strength of the head hardened rail steels and has shifted the stress concentration location away from the critical zone of the railhead end. The reduction in the magnitude and the relocation of the stress concentration in the SMRH design has been validated through a full scale wheel – railhead interaction test rig; Railhead strains under the loaded wheels have been recorded using a non-contact digital image correlation method. Experimental study has confirmed the accuracy of the numerical predications. Although the SMRH shaped IRJs eliminate stress singularities, they can still fail due to joint bar or bolt hole cracking; therefore, another conceptual design, termed as Embedded IRJ (EIRJ) in this thesis, with no joint bars and pre-tensioned bolts has been developed using a multi-objective optimisation formulation based on the coupled genetic algorithm – parametric finite element method. To achieve the required structural stiffness for the safe passage of the loaded wheels, the rails were embedded into the concrete of the post tensioned sleepers; the optimal solutions for the design of the EIRJ is shown to simplify the design through the elimination of the complex interactions and failure modes of the various structural components of the CIRJ. The practical applicability of the optimal shapes SMRH and EIRJ is demonstrated through two illustrative examples, termed as improved designs (IMD1 & IMD2) in this thesis; IMD1 is a combination of the CIRJ and the SMRH designs, whilst IMD2 is a combination of the EIRJ and SMRH designs. These two improved designs have been simulated for two key operating (speed and wagon load) and design (wheel diameter) parameters that affect the wheel-rail contact; the effect of these parameters has been found to be negligible to the performance of the two improved designs and the improved designs are in turn found far superior to the current designs of the CIRJs in terms of stress singularities and deformation under the passage of the loaded wheels. Therefore, these improved designs are expected to provide longer service life in relation to the CIRJs.

Modelling the failure of thin layered mortar joints in masonry

This paper deals with the failure of high adhesive, low compressive strength, thin layered polymer mortar joints in masonry through a contact modelling in finite element framework. Failure due to combined shear, tensile and compressive stresses are considered through a constitutive damaging contact model that incorporates traction–separation as a function of displacement discontinuity. The modelling method is verified using single and multiple contact analyses of thin mortar layered masonry specimens under shear, tensile and compressive stresses and their combinations. Using this verified method, the failure of thin mortar layered masonry under a range of shear to tension ratios and shear to compression ratios has been examined. Finally, this model is applied to thin bed masonry wallettes for their behaviour under biaxial tension–tension and compression–tension loadings perpendicular and parallel to the bed joints.

Transient air flow and heat transfer in a triangular enclosure with a conducting partition

Numerical investigation is carried out for natural convection heat transfer in an isosceles triangular enclosure partitioned in the centre by a vertical wall with infinite conductivity. A sudden temperature difference between two zones of the enclosure has been imposed to trigger the natural convection. As a result, heat is transferred between both sides of the enclosure through the conducting vertical wall with natural convection boundary layers forming adjacent to the middle partition and two inclined surfaces. The Finite Volume based software, Ansys 14.5 (Fluent) is used for the numerical simulations. The numerical results are obtained for different values of aspect ratio, A (0.2, 0.5 and 1.0) and Rayleigh number, Ra (10^5 <= Ra <= 10^8) for a fixed Prandtl number, Pr = 0.72 of air. It is anticipated from the numerical simulations that the coupled thermal boundary layers development adjacent to the partition undergoes several distinct stages including an initial stage, a transitional stage and a steady stage. Time dependent features of the coupled thermal boundary layers as well as the overall natural convection flow in the partitioned enclosure have been discussed in this study.

Steady natural convection of non-Newtonian power-law fluid in a trapezoidal enclosure

Numerically investigation of free convection heat transfer in a differentially heated trapezoidal cavity filled with non-Newtonian Power-law fluid has been performed in this study. The left inclined surface is uniformly heated whereas the right inclined surface is maintained as uniformly cooled. The top and bottom surfaces are kept adiabatic with initially quiescent fluid inside the enclosure. Finite volume based commercial software FLUENT 14.5 is used to solve the governing equations. Dependency of various flow parameters of fluid flow and heat transfer is analyzed including Rayleigh number, Ra ranging from 10^5 to 10^7, Prandtl number, Pr of 100 to 10,000 and power index, n of 0.6 to 1.4. Outcomes have been reported in terms of isotherms, streamline, and local Nusselt number for various Ra, Pr, n and inclined angles. Grid sensitivity analysis is performed and numerically obtained results have been compared with those results available in the literature and found good agreement.

How to increase the accuracy of analysis and reduce the computational time in ANSYS in the case of deformation study of orthopedic bone plates

Currently, finite element analyses are usually done by means of commercial software tools. Accuracy of analysis and computational time are two important factors in efficiency of these tools. This paper studies the effective parameters in computational time and accuracy of finite element analyses performed by ANSYS and provides the guidelines for the users of this software whenever they us this software for study on deformation of orthopedic bone plates or study on similar cases. It is not a fundamental scientific study and only shares the findings of the authors about structural analysis by means of ANSYS workbench. It gives an idea to the readers about improving the performance of the software and avoiding the traps. The solutions provided in this paper are not the only possible solutions of the problems and in similar cases there are other solutions which are not given in this paper. The parameters of solution method, material model, geometric model, mesh configuration, number of the analysis steps, program controlled parameters and computer settings are discussed through thoroughly in this paper.

Heat and mass transfer modeling of the osmo-convective drying of yacon roots (Smallanthus sonchifolius)

Osmotic treatments are often applied prior to convective drying of foods to impart sensory appeal aspects. During this process a multicomponent mass flow, composed mainly of water and osmotic agent, takes place. In this work, a heat and mass transfer model for the osmo-convective drying of yacon was developed and solved by the Finite Element Method using COMSOL Multiphysics®, considering a 2-D axisymmetric geometry and moisture dependent thermophysical properties. Yacon slices were osmotically dehydrated for 2 hours in a solution of sucralose and then dried in a tray dryer for 3 hours. The model was validated by experimental data of temperature, moisture content and sucralose uptake (R²> 0.90).

Effect of an adiabatic fin on natural convection heat transfer in a triangular enclosure

Natural convection thermal boundary layer adjacent to the heated inclined wall of a right angled triangle with an adiabatic fin attached to that surface is investigated by numerical simulations. The finite volume based unsteady numerical model is adopted for the simulation. It is revealed from the numerical results that the development of the boundary layer along the inclined surface is characterized by three distinct stages, i.e. a start-up stage, a transitional stage and a steady stage. These three stages can be clearly identified from the numerical simulations. Moreover, in presence of adiabatic fin, the thermal boundary layer adjacent to the inclined wall breaks initially. However, it is reattached with the downstream boundary layer next to the fin. More attention has been given to the boundary layer development near the fin area.

Behaviour and design of crest-fixed profiled steel roof claddings under wind uplift

Profiled steel roof claddings in Australia are commonly made of very thin high tensile steel and are crest-fixed with screw fasteners. At present the design of these claddings is entirely based on testing. In order to improve the understanding of the behaviour of these claddings under wind uplift, and thus the design methods, a detailed investigation consisting of a finite element analysis and laboratory experiments was carried out on two-span roofing assemblies of three common roofing profiles. It was found that the failure of the roof cladding system was due to a local failure (dimpling of crests/pull-through) at the fasteners. This paper presents the details of the investigation, the results and then proposes a design method based on the strength of the screwed connections, for which testing of small-scale roofing models and/or using a simple design formula is recommended.