Mechanical Engineering - Theses

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 175
  • Item
  • Item
    Thumbnail Image
    Analysis of numerical error on unstructured meshes: with applications to fluid dynamics
    Chenoweth, Samuel Keith Maxwell. (University of Melbourne, 2010)
  • Item
    Thumbnail Image
    Turbulent plumes in confined spaces
    George, Nitheesh ( 2019)
    Studies pertaining to turbulent plumes in confined spaces are of utmost interest due to its relevance in practical flows that are associated with, but not restricted to, the propagation of smoke and hot gases generated by fires in buildings, road and railway tunnels, etc. In this dissertation, direct numerical simulations (DNS) of the governing equations are carried out to analyze such flows with the focus on (i) free turbulent line plumes, and (ii) wall attached turbulent line plumes, in confined spaces. In all cases, the computation domain is rectangular with no-slip and adiabatic boundary conditions at the top, bottom, and lateral side walls. In free turbulent line plume simulations, the plume originates from a line heat source of length, L, located at the centre of the bottom wall and rises until it impinges on the top wall and eventually spreading out laterally thereby producing a buoyant fluid layer at the top wall. Since the region is confined, the continuous heat source forces the top layer to move downwards, until it reaches the bottom wall, when the flow is said to be at the asymptotic state (Baines and Turner 1969). DNS data at three Reynolds numbers (ReH), 1800, 3600 and 7200, based on box height H and the buoyant velocity scale, F_1/3 0 , where F_0 is buoyancy flux per unit length, are presented for plume lengths, L/H = 1, 2 and 4 and box aspect ratio, R/H = 1. Here, R is the box half-width. Following the initial transient dynamics, a flapping motion of the plume is observed, where the plume oscillates around the centre plane of the box. The DNS results reveal that the long-term behavior of the flow consists of a meandering, flapping plume with a counter-rotating vortex pair on either side of the plume. Additionally, the plume volume, momentum, and buoyancy fluxes obtained from the simulations are compared to the theoretical models proposed by Baines and Turner (1969) and Barnett (1991). Further, simulations of turbulent line plumes are carried out at increased box aspect ratios R/H = 1, 2, 4, 8 and 16, to study the horizontal outflow of the buoyant fluid layer after the plume impinges on the top wall. Following the axisymmetric plume model of Kaye and Hunt (2007), a theoretical model to compute the horizontal outflow properties is developed for turbulent line plumes. In the case of wall attached thermal plumes, the plume originates from a local line heat source placed at the bottom left corner of the box. The plume develops along the vertical side wall while remaining attached to it before spreading across the top wall forming a buoyant fluid layer and eventually moving downwards and filling the whole box. The simulations are carried out at ReH = 14530 and L/H = 0.5, and a parametric study is conducted for boxes of aspect ratios R/H = 1 and 2. Furthermore, the original filling box model of Baines and Turner Baines and Turner (1969) is modified to incorporate the wall shear stress and are compared against the results obtained from the DNS. A reasonable agreement is observed for the volume and momentum fluxes in the quiescent uniform environment and for the time-dependent buoyancy profiles calculated further away from the plume. Finally, the entrainment processes in both free and wall attached line plumes are assessed, using the DNS data. Both cases show similar contributions to entrainment due to net buoyancy. However, a deficit in the entrainment coefficient is observed for wall plumes due to the effect of the wall, which in turn suppressed the turbulent kinetic energy production.
  • Item
    Thumbnail Image
    The Application of Surface Effect Ships to Freight Transport
    Minty, Paul John ( 2019)
    The research project was conducted in two phases. The first phase, a policy-based investigation of the application of a suitable technology to commercial freight transport in Southeast Asia, concluded that the candidate technology that best met the sought objectives was a hybrid catamaran hovercraft configuration, or Surface Effect Ship (SES). In the second phase, a specific physical issue associated with the Surface Effect Ship (SES) reported in the available literature was investigated, specifically, the control of water-wave induced ship motions caused by air pressure resonance in the SES air cushion. Air pressure resonance was investigated experimentally using a custom designed and fabricated air chamber apparatus and associated numerical simulation validation studies. The novel phenomenon of increase in time-wise mean air cushion pressure at higher frequency oscillations of an air cushion volume was predicted using numerical simulation. Experimental observations were in agreement with the predicted behaviour. It was concluded that some dynamic characteristics of the SES that could inhibit the wider use of this type of craft were able to be replicated in laboratory apparatus and predicted with a numerical model, thereby facilitating future design tools for practical SESs.
  • Item
    Thumbnail Image
    Microstructures of TiAl additively manufactured by EBM and LMD
    Jacob, Johnson ( 2019)
    Titanium aluminide (mainly TiAl) based intermetallic alloys have superior mechanical properties up to 800 C and are much lighter than the currently used nickel superalloys. They have been developed for turbine engine applications in the aerospace, automobile and energy industries. Because of the inherent brittleness of TiAl, conventional manufacturing requires specially designed post processing, which is not cost and time effective. Additive manufacturing (AM) is considered as an alternative production technique, that can replace the conventional process. AM facilitates printing of complex TiAl components directly from alloy powders. Solid-state transformations (SST) associated with the unique layer-by-layer manufacturing and in situ cyclic heating in AM offer the potential to tailor the microstructure of components, creating opportunity for improving mechanical properties. The aim of this study is to systematically investigate the microstructural evolution at different processing conditions and the effect of in situ heat treatment on different solid-state transformations of gamma-TiAl alloy during additive manufacturing by Electron Beam Melting (EBM) and Laser Engineered Net Shaping (LENS). Nearly fully dense Ti-48Al-2Cr-2Nb samples with minimum defects were printed by EBM and LENS, at high, medium and low energy input conditions. Detailed microstructural characterisation along the build direction of samples printed by EBM at a high energy density (ED) of 6 J/mm2 was conducted. X-Ray Diffraction (XRD) studies and analytical electron microscopic analysis at different locations along the build direction of the sample revealed microstructural instabilities such as discontinuous coarsening (DC) and alpha2 decomposition, caused by the in situ cyclic heating with the peak temperature of each cycle approaching close to the melting point of TiAl and gradually reducing in intensity. Discontinuous coarsening with associated grain boundary migration was identified to be the dominant solid-state transformation taking place in alpha2+gamma primary lamellar structure (PLS), resulting in high coherence lamellar interphases in the TiAl alloy produced by both EBM and LENS. It was shown that DC replaces fine PLS by coarsened lamellae with reduced colony size, along with significant change in chemical composition and volume fraction of phases. Reduction in interphase energy achieved by the coarsening of lamellae, with an average coarsening ratio around 18.6, was identified as a driving force for DC. High chemical free energy in PLS mainly caused by supersaturated alpha2 phase with a non-equilibrium phase fraction of 26%, can be considered as another factor inducing DC. Ascertaining this assumption, alpha2 phase fraction approaches its equilibrium composition after 11 cycles of in-situ heating. Colony size also considerably reduced as a consequence of DC reaction. beta phase formation, which is identified as another side outcome of these driving forces, shows a variation in phase fraction with build temperature, from 3.5% at 950 C to 8% at 1050 C. DC, along with beta phase formation greatly influences microstructural morphology and mechanical properties of TiAl. A similar trend was observed in LENS fabricated samples. The major variation was the absence of beta formation in LENS, which is attributed to the lack of pre-heating during LENS. Alloy modifications such as addition of 0.5 Si (at.%) as experimented in this study, result in precipitate pinning at interphase, potentially reducing excessive coarsening. This study summarises the process of DC and beta phase formation in TiAl alloys during EBM and LENS processes and discusses the strategies to regulate DC and beta formation through AM of TiAl and utilize the intrinsic solid-state transformations in AM for controlling the microstructure to achieve application-specific material properties.
  • Item
    Thumbnail Image
    Investigation of Direct Combustion Noise in Turbulent Premixed Jet Flames Using Direct Numerical Simulations
    Brouzet, Davy Joao-Etienne ( 2019)
    Direct combustion noise plays a key role in initiating thermoacoustic instability in lean, premixed gas turbines. Moreover, many combustion devices produce a high level of noise while being subjected to increasingly stringent noise regulations. Therefore, achieving a better understanding of sound generation by premixed flames is crucial for designing safer and quieter combustion devices. With the advancement in High-Performance Computing (HPC), high fidelity simulations such as Direct Numerical Simulation (DNS) have received increasing attention as a means to improve our fundamental understanding of turbulent flames. This thesis aims to study the mechanism of sound generation in turbulent premixed jet flames using DNS and state-of-the-art post-processing methods. A DNS dataset featuring sound generation by turbulent premixed flames with simple chemistry is first analysed. Using Spectral Proper Orthogonal Decomposition (SPOD), two types of flame coherent structures responsible for combustion noise are identified. The first type arises in the jet's shear layer, originating from the Kelvin-Helmholtz (K-H) instability and is indirectly producing sound through the deformation of the flame front. The second type is found near the jet centreline and is linked to small, non-linear flame dynamics. Even though their energy content is lower than that of K-H structures, they are an important feature to explain the broadband nature of combustion noise. Then, a framework to identify the location and topology of annihilation events, and to study their generated sound, is presented. This study reveals that different topologies are similar in terms of the generated sound. In addition, a spectral analysis shows that flame annihilation is the physical mechanism by which air-fuel ratio affects the radiated sound amplitude at high frequencies. Finally, DNS datasets of turbulent premixed jet flames with a semi-global and a skeletal chemical mechanism are produced and analysed, to investigate the impact of chemical modelling on combustion noise. Large differences between the two cases are observed in the OASPL value and on the high-frequency side of the acoustic spectrum. Analysis of the acoustic source term resulting from the heat release rate fluctuations demonstrates that the post-flame region has minimal contribution in terms of sound generation. Furthermore, the most exothermic reaction in each mechanism is by far the dominant source of heat release rate fluctuations, and hence sound generation. It is observed that the OASPL discrepancy between the two chemical mechanisms arises from the differences in the peak amplitude of the heat release rate. Then, a modelling approach shows that the acoustic spectrum in the high frequency range results from highly curved flamelets and can be estimated from the flame curvature statistics. This approach demonstrates that the high-frequency acoustic discrepancy arises from a more wrinkled flame in the case featuring the more complex chemistry. Overall, to accurately predict the sound generated by turbulent premixed flames, a reduced chemical mechanism needs to correctly capture the flame response to turbulence.
  • Item
    Thumbnail Image
    Optimal Performance of a Wind Farm with and without Battery Energy Storage
    Naemi, Mostafa ( 2019)
    The increasing penetration of wind in power systems can result in various challenges for system security and reliability, as well as wind farm investability. These challenges require further research to develop an understanding of the operation of a wind farm with and without the use of complementary technologies, such as energy storage, particularly, when the wind turbines have Frequency Control Ancillary Services (FCAS) capability. This thesis therefore first presents a hierarchical, data-driven, reduced order model of a wind farm, which accounts for the correlation between wind turbines' power output. It finds that the cross-correlation between each turbines' power generation is dependent on the frequency of the wind disturbance and the distance between wind turbines, so that the cross-correlation can be related to the convective length scale in the incoming wind. An investigation of the number of wind turbines required to simulate the wind farm's power generation then indicates that there is an inherent trade-off between model accuracy and complexity. An optimisation model is then developed to investigate the optimal performance of the wind farm participating in the energy and FCAS markets with and without a battery storage system. This requires modelling of the battery costs along with the revenues in the energy and FCAS markets, as well as a simplified version of the Causer Pays Method used in the Australian National Electricity Market (NEM). The optimal performance of the Mt Mercer wind farm, located in Victoria, Australia, is then examined. It is found that the wind farm participation in the FCAS markets can improve its financial performance. This analysis shows that the wind farm mainly tends to participate in the FCAS lower regulation market due to the higher prices of this service and no requirement for curtailment before starting a dispatch interval (precurtailment). An investigation of the impact of wind generation forecast accuracy on the system performance also finds that with a better forecasting system, the total performance of the wind farm improves. Finally, the optimal integration of a lithium-ion battery into a wind farm is examined. Assessments find that the battery is not investable without substantial subsidies when the battery participates only in the energy market. However, it is also found that participation in the FCAS markets can significantly improve the battery's investability, but its financial viability is highly sensitive to FCAS prices. In addition, it is found that the introduction of wind farm frequency control capability reduces the optimal value and capacity of battery storage. Investigation of different wind generation forecasting systems identifies that the improvement of forecast accuracy also reduces the optimal battery value and capacity, as there are fewer opportunities for the battery to reduce the wind farm regulation payments.
  • Item
    Thumbnail Image
    Dynamic clamp analysis of ion channel function
    Deerasooriya, Yadeesha Hasalanka ( 2019)
    Ion channels regulate neuronal excitability by controlling the ion flow through the neuronal membrane. Therefore, neuronal ion channel dysfunction can directly impact the function of neurons in the brain and lead to numerous disorders such as epilepsy and autism. As such, characterizing ion channels can provide valuable insights into these disorders and facilitate treatment. Current ion channel characterization methods predominantly use the voltage clamp (VC) approach which typically involves multiple time-consuming step protocols to capture ion channel dynamics corresponding to different membrane voltages. As such, these approaches neither recapitulate the natural behavior of the ion channel in the brain nor provide the means to directly investigate the relationship between these dynamics and neuronal excitability. This thesis addresses the limitations in current VC methods and implements the dynamic clamp (DC) approach to characterize ion channels. DC represents a real-time closed-loop system that incorporates computational models of neuronal systems with real ion channels. It provides biologically more natural recordings and is capable of directly determining the impact of ion channel dynamics on neuronal excitability. I investigate two applications to characterize ion channels using DC in this thesis. First, ion channel kinetics are mathematically modelled based on DC recordings. Second, the ion channels are functionally characterized using features extracted from DC recordings. My approach to mathematically model fast kinetics of ion channels requires only DC recordings of short time durations. It utilizes efficient global optimization algorithms to estimate the model that best matches the recorded DC data. To further enhance the performance of the approach, I have identified an optimal DC stimulation strategy, and extended optimization method is proposed when recorded DC data is noisy. This approach was more accurate than two existing VC based methods. The model derived from this approach could also predict firing patterns of experimental data with high accuracy. The DC-based approach was next extended to model slow kinetics together with fast kinetics. This involved a new DC stimulation strategy to record DC data and a three-step post-hoc optimization to estimate model parameters. The extended approach could estimate models that accurately predict AP firing patterns that occur during longer/sustained stimulation and could recreate the AP firing seen in experimental data. This thesis also proposes a workflow to functionally characterize ion channel variants such as mutations using DC data. The workflow creates a two-dimensional map of the variants where their positionings correspond to their functional characteristics. Multiple variants of two major neuronal sodium channels, Nav1.1 and Nav1.2, were investigated and their two-dimensional mappings were determined. The results suggested a clear functional separation between variants, not only corroborating the findings of previous conventional functional studies but also providing new insights into variant functionality. The two applications of the DC presented in this thesis demonstrate the potential of DC for ion channel characterization. Collectively, they provide crucial information on ion channel dynamics that will assist the development of effective treatments for neurological disorders involving mutant channels and will enable assessing the direct impact of pharmacological interventions on neuronal excitability.
  • Item
    Thumbnail Image
    Understanding loss mechanisms in turbomachinery to increase efficiency
    Akolekar, Harshal Deepak ( 2016)
    The improvement in the design of any mechanical device is always carried out to increase its efficiency. Turbomachinery, more specifically, the low pressure turbine (LPT) in aircraft engines, is no exception. In order to enhance its efficiency, the specific fuel consumption has to be reduced, which implies different loss mechanisms (mechanisms of entropy generation) in the LPT have to be studied. Conventional loss calculations consider a control volume around the blade and find the total loss in that region. However, there are a number of sources of losses around LPT blades, prominently the mixing out of incident wakes from the stator to the rotor. In order to design more efficient blades, knowing the total loss in the blade control volume is not enough. It is important to quantify the various sources of loss due to wake mixing. Denton has derived an analytical expression for losses due to mixing, which consists of 3 terms - losses from the trailing edge region, losses due to boundary layer effects, and losses due to blade blockage effects. This equation, however, is built upon a number of assumptions such as steady, incompressible flow conditions which are not realised for real flows. Denton does not define a distinctive trailing edge region, methods to calculate boundary layer thickness and base pressures. In order to enhance the applicability of Denton’s mixed out loss equation, it is important to identify the dependence of Denton’s equation on these parameters; which is an important objective of this work. This work aims to improve the robustness of Denton’s equation by analysing the effect of not defining the trailing edge region properly on the total loss; and then proposes 4 criteria to define the trailing edge region. This analysis was conducted on a LPT blade with steady flow conditions and it was found that the losses calculated with the aid of the trailing edge criteria lie within 5% of the loss in the blade control volume. Sensitivity analyses on the total loss from Denton’s equation has been conducted using boundary layer thickness criteria, base pressure and non-uniform flow averaging techniques as input parameters. It was found that there are large variances in the total loss, if the boundary layer thickness is not defined correctly, and therefore researchers need to be very careful and consistent in their selection of boundary layer thickness. Visual observation, along with the boundary layer thickness criteria identified in this work, will serve as a good method to determine the boundary layer thickness. Denton’s equation does not deal with unsteady dissipation effects. In order to analyse those effects, an analysis has been conducted on an unsteady flow field consisting of incoming wakes; which has been broken up into 20 quasi-steady phases. Denton’s losses were calculated for each of these phases in an attempt to understand the effects of unsteady flow on Denton’s equation. Based on the analysis conducted, it was concluded that better averaging techniques are required to average the losses from the 20 phases to quantify the losses due to unsteady effects. Overall, Denton’s analysis, which is dependent on a number of criteria, has the potential to give a good rough estimate of the loss source. Due to the absence of more accurate methods, it is helpful for designers; and if conducted by the same user, allows for good qualitative comparison of different flow configurations. It is however difficult to quantify the exact amounts of losses based on this method. For a quantitative comparison with total losses in the control volume, it is probably necessary to come up with new ways of quantifying the different sources. The report provides ideas for future analyses to further improve the understanding of loss mechanisms via Denton’s equation and other methods.
  • Item
    Thumbnail Image
    Computational coronary arterial fluid dynamics: from stenotic to rough-wall flow
    Li, Shuang ( 2018)
    Coronary heart disease (CAD) is a leading cause of mortality. CAD is usually caused by the built-up plaque, clinically known as stenosis, that narrows the arteries and hence limits blood flow to the heart. It has been found the development of stenosis is closely associate with local haemodynamics, which can be altered by the changes in arterial shape, such as plaque formation or stenting. The objective of this study is to better understand the post-narrowing and in-stent haemodynamic environments to gain insights into the flow physics associated with the stenosis development. Computational fluid dynamics (CFD) technique is used to firstly simulate pulsatile blood flow in full-scale straight and curved stenotic coronary arteries under a physiological inlet velocity waveform. Flow-pressure relation in simplified multiple sequential stenotic flow is then discussed. Lastly, investigations are carried out in simplified smooth models superposed with the egg carton type surface to mimic in-stent coronary arterial flow. The findings of flow characteristics can contribute towards the future prediction and diagnosis of coronary-related complications. Straight and curved arterial flows with three different degrees of stenosis are studied in both Newtonian and non-Newtonian fluids. The time-dependent inertial momentum is found to contribute to reverse flow development in the proximity of post-stenotic region and negatively correlate with the reverse flow size. Flow velocity and WSS in Newtonian and non-Newtonian fluids exhibit larger difference in response to the increase of stenosis degree. In the presence of curvature, low WSS is found to concentrate at the inner wall after the stenosis collocated with the reverse flow region. With the progressive stenosis severity, the secondary flow morphology distal to the stenosis also evolves into double-paired vortices structure and promotes the growth of reverse flow size. For non-Newtonian flow, smaller reverse flow bubble distal to the stenosis are observed and the difference in Newtonian and non-Newtonian fluids is more profound in higher degree stenosis cases. Overall, the haemodynamic behaviours downstream of stenosis are affected simultaneously by stenosis degree, the instant inertial momentum and secondary flow morphology if curvature presented. The relative location of low WSS (and reverse flow) becomes the potential trigger for the growth of stenosis. The correlation between the haemodynamics in post-stenotic region and potential clinical complications implies the necessity to determine the severity of stenosis. Virtual Fractional Flow Reserve (vFFR), a computational technique calculating pressure drop across stenosis, is considered as an adjunct application to invasive determined FFR, a current standard of clinical practice. Flow-pressure relation across multiple stenoses are analysed using both experimental and numerical approaches. Linear correlation between pressure drop and flow rate irrespective the number of stenosis is found. Negligible difference between steady and pulsatile flows is also observed. The conclusion may improve the clinical applicability of vFFR. In the in-stent coronary arterial flow, both increasing the roughness height and decreasing spacing reduces the shear rates (due to the increased proportion of pressure drag) near the trough of roughness, and hence encourage reverse flow formation. In non-Newtonian fluids, elevated relative viscosities are pronounced near the trough of the roughness while low viscosities are found around the peak of the roughness. This trend becomes more profound by increasing roughness height or decreasing wavelength. As a result, reverse flow is less likely to occur near the trough of roughness in non-Newtonian fluid. By comparing time-averaged velocity and WSS using different blood rheology models, the results show consistency in both qualitative and quantitative perspectives and suggest an interchangeable aspect of rheology models in simulations.