Adaption and use of a compressible flow code for turbomachinery design

The present paper presents and discusses the use of dierent codes regarding the numerical simulation of a radial-in ow turbine. A radial-in ow turbine test case was selected from published literature [1] and commercial codes (Fluent and CFX) were used to perform the steady-state numerical simulations. An in-house compressible- ow simulation code, Eilmer3 [2] was also adapted in order to make it suitable to perform turbomachinery simulations and preliminary results are presented and discussed. The code itself as well as its adaptation, comprising the addition of terms for the rotating frame of reference, programmable boundary conditions for periodic boundaries and a mixing plane interface between the rotating and non-rotating blocks are also discussed. Several cases with dierent orders of complexity in terms of geometry were considered and the results were compared across the dierent codes. The agreement between these results and published data is also discussed.

Groundwater visualisation system (GVS) : a software framework for integrated display and interrogation of conceptual hydrogeological models, data and time-series animation

Management of groundwater systems requires realistic conceptual hydrogeological models as a framework for numerical simulation modelling, but also for system understanding and communicating this to stakeholders and the broader community. To help overcome these challenges we developed GVS (Groundwater Visualisation System), a stand-alone desktop software package that uses interactive 3D visualisation and animation techniques. The goal was a user-friendly groundwater management tool that could support a range of existing real-world and pre-processed data, both surface and subsurface, including geology and various types of temporal hydrological information. GVS allows these data to be integrated into a single conceptual hydrogeological model. In addition, 3D geological models produced externally using other software packages, can readily be imported into GVS models, as can outputs of simulations (e.g. piezometric surfaces) produced by software such as MODFLOW or FEFLOW. Boreholes can be integrated, showing any down-hole data and properties, including screen information, intersected geology, water level data and water chemistry. Animation is used to display spatial and temporal changes, with time-series data such as rainfall, standing water levels and electrical conductivity, displaying dynamic processes. Time and space variations can be presented using a range of contouring and colour mapping techniques, in addition to interactive plots of time-series parameters. Other types of data, for example, demographics and cultural information, can also be readily incorporated. The GVS software can execute on a standard Windows or Linux-based PC with a minimum of 2 GB RAM, and the model output is easy and inexpensive to distribute, by download or via USB/DVD/CD. Example models are described here for three groundwater systems in Queensland, northeastern Australia: two unconfined alluvial groundwater systems with intensive irrigation, the Lockyer Valley and the upper Condamine Valley, and the Surat Basin, a large sedimentary basin of confined artesian aquifers. This latter example required more detail in the hydrostratigraphy, correlation of formations with drillholes and visualisation of simulation piezometric surfaces. Both alluvial system GVS models were developed during drought conditions to support government strategies to implement groundwater management. The Surat Basin model was industry sponsored research, for coal seam gas groundwater management and community information and consultation. The “virtual” groundwater systems in these 3D GVS models can be interactively interrogated by standard functions, plus production of 2D cross-sections, data selection from the 3D scene, rear end database and plot displays. A unique feature is that GVS allows investigation of time-series data across different display modes, both 2D and 3D. GVS has been used successfully as a tool to enhance community/stakeholder understanding and knowledge of groundwater systems and is of value for training and educational purposes. Projects completed confirm that GVS provides a powerful support to management and decision making, and as a tool for interpretation of groundwater system hydrological processes. A highly effective visualisation output is the production of short videos (e.g. 2–5 min) based on sequences of camera ‘fly-throughs’ and screen images. Further work involves developing support for multi-screen displays and touch-screen technologies, distributed rendering, gestural interaction systems. To highlight the visualisation and animation capability of the GVS software, links to related multimedia hosted online sites are included in the references.

Blast induced ground shock effects on pile foundations

Due to increased number of terrorist attacks in recent years, loads induced by explosions need to be incorporated in building designs. For safer performance of a structure, its foundation should have sufficient strength and stability. Therefore, prior to any reconstruction or rehabilitation of a building subjected to blast, it is important to examine adverse effects on the foundation caused by blast induced ground shocks. This paper evaluates the effects of a buried explosion on a pile foundation. It treats the dynamic response of the pile in saturated sand, using explicit dynamic nonlinear finite element software LS-DYNA. The blast induced wave propagation in the soil and the horizontal deformation of pile are presented and the results are discussed. Further, a parametric study is carried out to evaluate the effect of varying the explosive shape on the pile response. This information can be used to evaluate the vulnerability of piled foundations to credible blast events as well as develop guidance for their design.

Metal hexaborides with Sc, Ti or Mn

Comparison of well-determined single crystal data for stoichiometric, or near-stoichiometric, metal hexaborides con-firm previously identified lattice parameter trends using powder diffraction. Trends for both divalent and trivalent forms suggest that potential new forms for synthesis include Sc and Mn hexaborides. Density Functional Theory (DFT) calculations for KB6, CaB6, YB6, LaB6, boron octahedral clusters and Sc and Mn forms, show that the shapes of bonding orbitals are defined by the boron framework. Inclusion of metal into the boron framework induces a reduction in energy ranging from 1 eV to 6 eV increasing with ionic charge. For metals with d1 character, such a shift in energy brings a doubly degenerate band section along the G-M reciprocal space direction within the conduction bands tangential to the Fermi surface. ScB6 band structure and density of states calculations show directional and gap characteristics similar to those of YB6 and LaB6. These calculations for ScB6 suggest it may be possible to realize superconductivity in this compound if synthesized.

Lateral impact response and parametric studies of axially loaded square concrete filled steel tube columns

This paper presents a numerical study on the response of axially loaded slender square concrete filled steel tube (CFST) columns under low velocity lateral impact loading. A finite element analysis (FEA) model was developed using the explicit dynamic nonlinear finite element code LS -DYNA in which the strain rate effects of both steel and concrete, contact between steel tube and concrete and confinement effect provided by the steel tube for the concrete were considered. The model also benefited from a relatively recent feature of LS-DYNA for applying a pre-loading in the explicit solver. The developed numerical model was verified for its accuracy and adequacy by comparing the results with experimental results available in the literature. The verified model was then employed to conduct a parametric study to investigate the influence of axial load level, impact location, support conditions, and slenderness ratio on the response of the CFST columns. A good agreement between the numerical and experimental results was achieved. The model could reasonably predict the impact load-deflection history and deformed shape of the column at the end of the impact event. The results of the parametric study showed that whilst impact location, axial load level and slenderness ratio can have a significant effect on the peak impact force, residual lateral deflection and maximum lateral deflection, the influence of support fixity is minimal. With an increase of axial load to up to a certain level, the peak force increases; however, a further increase in the axial load causes a decrease in the peak force. Both residual lateral deflection and maximum lateral deflection increase as axial load level increases. Shifting the impact location towards the supports increases the peak force and reduces both residual and maximum lateral deflections. A rise in slenderness ratio decreases the peak force and increases the residual and maximum lateral deflections.

The effect of surface tension and kinetic undercooling on a radially-symmetric melting problem

The addition of surface tension to the classical Stefan problem for melting a sphere causes the solution to blow up at a finite time before complete melting takes place. This singular behaviour is characterised by the speed of the solid-melt interface and the flux of heat at the interface both becoming unbounded in the blow-up limit. In this paper, we use numerical simulation for a particular energy-conserving one-phase version of the problem to show that kinetic undercooling regularises this blow-up, so that the model with both surface tension and kinetic undercooling has solutions that are regular right up to complete melting. By examining the regime in which the dimensionless kinetic undercooling parameter is small, our results demonstrate how physically realistic solutions to this Stefan problem are consistent with observations of abrupt melting of nanoscaled particles.

Modal flexibility method for structural damage detection in suspension bridges

Suspension bridges meet the steadily growing demand for lighter and longer bridges in today’s infrastructure systems. These bridges are designed to have long life spans, but with age, their main cables and hangers could suffer from corrosion and fatigue. There is a need for a simple and reliable procedure to detect and locate such damage, so that appropriate retrofitting can be carried out to prevent bridge failure. Damage in a structure causes changes in its properties (mass, damping and stiffness) which in turn will cause changes in its vibration characteristics (natural frequencies, modal damping and mode shapes). Methods based on modal flexibility, which depends on both the natural frequencies and mode shapes, have the potential for damage detection. They have been applied successfully to beam and plate elements, trusses and simple structures in reinforced concrete and steel. However very limited applications for damage detection in suspension bridges have been identified to date. This paper examines the potential of modal flexibility methods for damage detection and localization of a suspension bridge under different damage scenarios in the main cables and hangers using numerical simulation techniques. Validated finite element model (FEM) of a suspension bridge is used to acquire mass normalized mode shape vectors and natural frequencies at intact and damaged states. Damage scenarios will be simulated in the validated FE models by varying stiffness of the damaged structural members. The capability of damage index based on modal flexibility to detect and locate damage is evaluated. Results confirm that modal flexibility based methods have the ability to successfully identify damage in suspension bridge main cables and hangers.

Comparison of functionals for metal hexaboride band structure calculations

Density functional calculations of the electronic band structure for superconducting and semi-conducting metal hexaborides are compared using a consistent suite of assumptions and with emphasis on the physical implications of computed models. Spin polarization enhances mathematical accuracy of the functional approximations and adds significant physical meaning to model interpretation. For YB6 and LaB6, differences in alpha and beta projections occur near the Fermi energy. These differences are pronounced for superconducting hexaborides but do not occur for other metal hexaborides.

Eigenstrain boundary integral equations with local Eshelby matrix for stress analysis of ellipsoidal particles

Aiming at the large scale numerical simulation of particle reinforced materials, the concept of local Eshelby matrix has been introduced into the computational model of the eigenstrain boundary integral equation (BIE) to solve the problem of interactions among particles. The local Eshelby matrix can be considered as an extension of the concepts of Eshelby tensor and the equivalent inclusion in numerical form. Taking the subdomain boundary element method as the control, three-dimensional stress analyses are carried out for some ellipsoidal particles in full space with the proposed computational model. Through the numerical examples, it is verified not only the correctness and feasibility but also the high efficiency of the present model with the corresponding solution procedure, showing the potential of solving the problem of large scale numerical simulation of particle reinforced materials.

Thermodynamic modeling of ethanol fumigation in a diesel engine

Due to rapidly diminishing international supplies of fossil fuels, such as petroleum and diesel, the cost of fuel is constantly increasing, leading to higher costs of living, as a result of the significant reliance of many industries on motor vehicles. Many technologies have been developed to replace part or all of a fossil fuel with bio-fuels. One of the dual fuel technologies is fumigation of ethanol in diesel engines, which injects ethanol into the intake air stream of the engine. The advantage of this is that it avoids any costly modification of the engine high pressure diesel injection system, while reducing the volume of diesel required and potentially increasing the power output and efficiency. This paper investigates the performance of a diesel engine, converted to implement ethanol fumigation. The project will use both existing experimental data, along with generating computer modeled results using the program AVL Boost. The data from both experiments and the numerical simulation indicate desirable results for the peak pressure and the indicated mean effective pressure (IMEP). Increase in ethanol substitution resulted in elevated combustion pressure and an increase in the IMEP, while the variation of ethanol injection location resulted in negligible change. These increases in cylinder pressure led to a higher work output and total efficiency in the engine as the ethanol substitution was increased. In comparing the numerical and experimental results, the simulation showed a slight elevation, due to the inaccuracies in the heat release models. Future work is required to improve the combustion model and investigate the effect of the variation of the location of ethanol injection.

Blast response and vulnerability assessment of piled foundations

This research treated the response of pile foundations to blast loads. The influence of important parameters was investigated. The research techniques and the results will enable safer design of pile foundations that are vulnerable to blast loads.

Vertical displacement measurement using fibre Bragg grating (FBG) sensors for structural health monitoring of bridges

This thesis developed a practical, cost effective, easy-to-use method for measuring the vertical displacements of bridges using fiber Bragg grating (FBG) sensors, which includes the curvature and inclination approaches. These approaches were validated by the numerical simulation tests on a full scale bridge and the laboratory-based tests. In doing so, a novel frictionless FBG inclination sensor with extremely high sensitivity and resolution has also been developed and validated.

Low-pressure planar magnetron discharge for surface deposition and nanofabrication

Current-voltage characteristics of the planar magnetron are studied experimentally and by numerical simulation. Based on the measured current-voltage characteristics, a model of the planar magnetron discharge is developed with the background gas pressure and magnetic field used as parameters. The discharge pressure was varied in a range of 0.7-1.7 Pa, the magnetic field of the magnetron was of 0.033-0.12 T near the cathode surface, the discharge current was from 1 to 25 A, and the magnetic field lines were tangential to the substrate surface in the region of the magnetron discharge ignition. The discharge model describes the motion of energetic secondary electrons that gain energy by passing the cathode sheath across the magnetic field, and the power required to sustain the plasma generation in the bulk. The plasma electrons, in turn, are accelerated in the electric field and ionize effectively the background gas species. The model is based on the assumption about the prevailing Bohm mechanism of electron conductivity across the magnetic field. A criterion of the self-sustained discharge ignition is used to establish the dependence of the discharge voltage on the discharge current. The dependence of the background gas density on the current is also observed from the experiment. The model is consistent with the experimental results. © 2010 American Institute of Physics.

Plasma Nanoscience : from nature's mastery to deterministic plasma-aided nanofabrication

This paper introduces the plasma-nanoscience research area and shows the way from Nature's mastery in assembling nanosized dust grains in the Universe to deterministic plasma-aided nanofabrication. The concept of deterministic nanoassembly is explained, and the multidisciplinary approach to bridge the spatial gap of nine orders of magnitude between the sizes of plasma reactors and atomic building units is discussed. Ongoing numerical simulation and experimental efforts on highly controlled synthesis of carbon nanotip and semiconducting quantum-dot structures show potential benefits of using ionized-gas environments in nanofabrication. © 2007 IEEE.

Surface fluxes of Si and C adatoms at initial growth stages of SiC quantum dots

Self-assembly of highly stoichiometric SiC quantum dots still remains a major challenge for the gas/plasma-based nanodot synthesis. By means of a multiscale hybrid numerical simulation of the initial stage (0.1-2.5 s into the process) of deposition of SiCSi (100) quantum dot nuclei, it is shown that equal Si and kst atom deposition fluxes result in strong nonstoichiometric nanodot composition due to very different surface fluxes of Si and C adatoms to the quantum dots. At this stage, the surface fluxes of Si and C adatoms to SiC nanodots can be effectively controlled by manipulating the SiC atom influx ratio and the Si (100) surface temperature. It is demonstrated that at a surface temperature of 800 K the surface fluxes can be equalized after only 0.05 s into the process; however, it takes more then 1 s at a surface temperature of 600 K. Based on the results of this study, effective strategies to maintain a stoichiometric ([Si] [C] =1:1) elemental ratio during the initial stages of deposition of SiCSi (100) quantum dot nuclei in a neutral/ionized gas-based process are proposed.