Experimental and numerical investigation of strain-rate-dependent behavior of kangaroo shoulder cartilage

This thesis introduces a new animal model, kangaroo, to biomechanical investigations of shoulder cartilage research. It examines the effect of cartilage structure and constituents on tissue behavior and its adaptation to mechanical loading. In doing so, the study explains the relationship of tissue's functional behaviors to its structure and constituents which has important implications for tissue engineering strategies catering joint specific cartilage tissue generation.

Effect of magnetic and electrical fields on dendritic freezing of aqueous solutions of sodium chloride

Aqueous solutions of sodium chloride were solidified under the influence of magnetic and electrical fields using two different freezing systems. In the droplet system, small droplets of the solution are introduced in an organic liquid column at −20°C which acts as the heat sink. In the unidirectional freezing system the solutions are poured into a tygon tube mounted on a copper chill, maintained at −70°C, from which the freezing initiates. Application of magnetic fields caused an increase in the spacing and promoted side branching of primary ice dendrites in the droplet freezing system, but had no measurable effect on the dendrites formed in the unidirectional freezing system. The range of electric fields applied in this investigation had no measurable effect on the dendritic structure. Possible interactions between external magnetic and electrical fields have been reviewed and it is suggested that the selective effect of magnetic fields on dendrite spacings in a droplet system could be due to a change in the nucleation behaviour of the solution in the presence of a magnetic field.

Solutions of Boundary Layer Equations Governing Free Convective and Magnetohydrodynamic Flows through Parametric Differentiation

In der vorliegenden Arbeit wird die Methode der parametrischen Differentiation angewendet, um ein System nichtlinearer Gleichungen zu lösen, das zwei- und dreidimensionale freie, konvektive Grenzschichströmungen bzw. eine zweidimensionale magnetohydrodynamische Grenzschichtströmung beherrscht. Der Hauptvorteil dieser Methode besteht darin, daß die nichlinearen Gleichungen auf lineare reduziert werden und die Nichtlinearität auf ein System von Gleichungen erster Ordnung beschränkt wird, das, verglichen mit den ursprünglichen Nichtlinearen Gleichungen, viel leichter gelöst werden kann. Ein anderer Vorzug der Methode ist, daß sie es ermöglicht, die Lösung von einer bekannten, zu einem bestimmten Parameterwert gehörigen Lösung aus durch schrittweises Vorgehen die Lösung für den gesamten Parameterbereich zu erhalten. Die mit dieser Methode gewonnenen Ergebnisse stimmen gut mit den entsprechenden, mit anderen numerischen Verfahren erzielten überein.

A note on non-separable solutions of linear partial differential equation

A method has been presented for constructing non-separable solutions of homogeneous linear partial differential equations of the type F(D, D′)W = 0, where D = ∂/∂x, D′ = ∂/∂y, Image where crs are constants and n stands for the order of the equation. The method has also been extended for equations of the form Φ(D, D′, D″)W = 0, where D = ∂/∂x, D′ = ∂/∂y, D″ = ∂/∂z and Image As illustration, the method has been applied to obtain nonseparable solutions of the two and three dimensional Helmholtz equations.

Regularized numerical integration of multibody dynamics with the generalized alpha method*

This paper discusses the consistent regularization property of the generalized α method when applied as an integrator to an initial value high index and singular differential-algebraic equation model of a multibody system. The regularization comes from within the discretization itself and the discretization remains consistent over the range of values the regularization parameter may take. The regularization involves increase of the smallest singular values of the ill-conditioned Jacobian of the discretization and is different from Baumgarte and similar techniques which tend to be inconsistent for poor choice of regularization parameter. This regularization also helps where pre-conditioning the Jacobian by scaling is of limited effect, for example, when the scleronomic constraints contain multiple closed loops or singular configuration or when high index path constraints are present. The feed-forward control in Kane's equation models is additionally considered in the numerical examples to illustrate the effect of regularization. The discretization presented in this work is adopted to the first order DAE system (unlike the original method which is intended for second order systems) for its A-stability and same order of accuracy for positions and velocities.

Numerical computation of the module of a quadrilateral

The module of a quadrilateral is a positive real number which divides quadrilaterals into conformal equivalence classes. This is an introductory text to the module of a quadrilateral with some historical background and some numerical aspects. This work discusses the following topics: 1. Preliminaries 2. The module of a quadrilateral 3. The Schwarz-Christoffel Mapping 4. Symmetry properties of the module 5. Computational results 6. Other numerical methods Appendices include: Numerical evaluation of the elliptic integrals of the first kind. Matlab programs and scripts and possible topics for future research. Numerical results section covers additive quadrilaterals and the module of a quadrilateral under the movement of one of its vertex.

Experimental and numerical investigation of jet injection in a wall bounded supersonic flow

Numerical and experimental studies of a supersonic jet (Helium) inclined at 45 degrees to a oncoming Mach 2 flow have been carried out. The numerical study has been used to arrive at a geometry that could reduce an oncoming Mach 5.75 flow to Mach 2 flow and in determining the jet parameters. Experiments are carried out in the IISc. hypersonic shock tunnel HST2 at similar conditions obtained from numerical studies. Flow visualization studies carried out using Schlieren technique clearly show the presence of the bow shock in front of the jet exposed to supersonic cross flow. The jet Mach number is experimentally found to be approximate to 3. Visual observations show that the jet has penetrated up to 60% of the total height of the chamber.

Mathematical modelling and numerical simulation of phosphine flow during grain fumigation in leaky cylindrical silos

The phosphine distribution in a cylindrical silo containing grain is predicted. A three-dimensional mathematical model, which accounts for multicomponent gas phase transport and the sorption of phosphine into the grain kernel is developed. In addition, a simple model is presented to describe the death of insects within the grain as a function of their exposure to phosphine gas. The proposed model is solved using the commercially available computational fluid dynamics (CFD) software, FLUENT, together with our own C code to customize the solver in order to incorporate the models for sorption and insect extinction. Two types of fumigation delivery are studied, namely, fan- forced from the base of the silo and tablet from the top of the silo. An analysis of the predicted phosphine distribution shows that during fan forced fumigation, the position of the leaky area is very important to the development of the gas flow field and the phosphine distribution in the silo. If the leak is in the lower section of the silo, insects that exist near the top of the silo may not be eradicated. However, the position of a leak does not affect phosphine distribution during tablet fumigation. For such fumigation in a typical silo configuration, phosphine concentrations remain low near the base of the silo. Furthermore, we find that half-life pressure test readings are not an indicator of phosphine distribution during tablet fumigation.

Numerical study of heat transfer in laminar film boiling by the finite-difference method

The finite-difference form of the basic conservation equations in laminar film boiling have been solved by the false-transient method. By a judicious choice of the coordinate system the vapour-liquid interface is fitted to the grid system. Central differencing is used for diffusion terms, upwind differencing for convection terms, and explicit differencing for transient terms. Since an explicit method is used the time step used in the false-transient method is constrained by numerical instability. In the present problem the limits on the time step are imposed by conditions in the vapour region. On the other hand the rate of convergence of finite-difference equations is dependent on the conditions in the liquid region. The rate of convergence was accelerated by using the over-relaxation technique in the liquid region. The results obtained compare well with previous work and experimental data available in the literature.

Experimental and numerical studies of hollow flange channel beams subject to web crippling under ETF and ITF load cases

This paper presents the details of experimental and numerical studies on the web crippling behaviour of hollow flange channel beams, known as LiteSteel beams (LSB). The LSB has a unique shape of a channel beam with two rectangular hollow flanges, made using a unique manufacturing process. Experimental and numerical studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions, predominant shear actions and combined actions. To date, however, no investigation has been conducted into the web crippling behaviour and strength of LSB sections under ETF and ITF load conditions. Hence experimental studies consisting of 28 tests were first conducted in this research to assess the web crippling behaviour and strengths of LSBs under two flange load cases (ETF and ITF). Experimental web crippling capacity results were then compared with the predictions from AS/NZS 4600 and AISI S100 design rules, which showed that AS/NZS 4600 and AISI S100 design equations are very unconservative for LSBs under ETF and ITF load cases. Hence improved equations were proposed to determine the web crippling capacities of LSBs. Finite element models of the tested LSBs were then developed, and used to determine the elastic buckling loads of LSBs under ETF and ITF load cases. New equations were proposed to determine the corresponding elastic buckling coefficients of LSBs. Finally suitable design rules were also developed under the Direct Strength Method format using the test results and buckling analysis results from finite element analyses.

Numerical modelling of rivet fastened rectangular hollow flange channel beams subject to shear

The intermittently rivet fastened Rectangular Hollow Flange Channel Beam (RHFCB) is a new cold-formed hollow section proposed as an alternative to welded hollow flange beams. Many experimental and numerical studies have been carried out in the past to investigate the shear behaviour of lipped channel beams. However, no research has been undertaken on the shear behaviour of rivet fastened RHFCBs. Therefore experimental and numerical studies were undertaken to investigate the shear behaviour and strength of rivet fastened RHFCBs. In this research finite element models of rivet fastened RHFCBs were developed to investigate their nonlinear shear behaviour including their buckling characteristics and ultimate shear strength. This paper presents the details of the finite element models of rivet fastened RHFCBs and the results. Both finite element analysis and experimental results showed that the current design rules are very conservative for the shear design of rivet fastened RHFCBs. Appropriate improvements have been proposed for the design rules of shear strength of rivet fastened RHFCBs within the Direct Strength Method format.

Material characterisation and numerical modelling of gypsum plasterboards in fire

The fire resistance characteristic of LSF wall systems mainly depends on the protective linings in use, commonly gypsum plasterboards. However, unclassified boards with varying composition and more notably with ambiguous thermal properties are increasingly becoming available in the market. Therefore a study was undertaken with an aim to set minimum standards for fire protective boards used in LSF wall applications. This paper presents the details of this study based on material characterisation and finite element thermal modelling of the most commonly used fire protective board, gypsum plasterboards, to address these critical issues related to fire safety design. In the material characterisation phase of this study, thermal properties of three different gypsum plasterboards manufactured in Australia were measured, analysed and compared. Subsequently, it proposes a thermal property based “k-factor” capable of giving an overall measure of the fire performance of boards, so that it can be used in appropriately classifying fire protective boards. As it is not known how this factor relates to the overall fire performance of LSF wall systems, numerical models were also developed and used to simulate the performance of LSF walls exposed to the standard fire. Finally, a correlation between time-temperature profiles from numerical analyses and calculated k-factors was established.

Entire Solutions of Hydrodynamical Equations with Exponential Dissipation

We consider a modification of the three-dimensional Navier-Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as e(vertical bar k vertical bar/kd) at high wavenumbers vertical bar k vertical bar. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than e-(C(k/kd) ln(vertical bar k vertical bar/kd)) for any C < 1/(2 ln 2). The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with C = C-* = 1/ ln 2. The same behavior with a universal constant C-* is conjectured for the Navier-Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier-Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed.

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.

A numerical study of projectile impact on thin aluminium plates

This article deals with a simulation-based Study of the impact of projectiles on thin aluminium plates using LS-DYNA by modelling plates with shell elements and projectiles with solid elements. In order to establish the required modelling criterion in terms of element size for aluminium plates, a convergence Study of residual velocity has been carried Out by varying mesh density in the impact zone. Using the preferred material and meshing criteria arrived at here, extremely good prediction of test residual velocities and ballistic limits given by Gupta et al. (2001) for thin aluminium plates has been obtained. The simulation-based pattern of failure with localized bulging and jagged edge of perforation is similar to the perforation with petalling seen in tests. A number Of simulation-based parametric studies have been carried out and results consistent with published test data have been obtained. Despite the robust correlation achieved against published experimental results, it would be prudent to conduct one's own experiments, for a final correlation via the present modelling procedure and analysis with the explicit LS-DYNTA 970 solver. Hence, a sophisticated ballistic impact testing facility and a high-speed camera have been used to conduct additional tests on grade 1100 aluminium plates of 1 mm thickness with projectiles Of four different nose shapes. Finally, using the developed numerical simulation procedure, an excellent correlation of residual velocity and failure modes with the corresponding test results has been obtained.