Modelling of arc welding : the importance of including the arc plasma in the computational domain

We present results of computational simulations of tungsten-inert-gas and metal-inert-gas welding. The arc plasma and the electrodes (including the molten weld pool when necessary) are included self-consistently in the computational domain. It is shown, using three examples, that it would be impossible to accurately estimate the boundary conditions on the weld-pool surface without including the arc plasma in the computational domain. First, we show that the shielding gas composition strongly affects the properties of the arc that influence the weld pool: heat flux density, current density, shear stress and arc pressure at the weld-pool surface. Demixing is found to be important in some cases. Second, the vaporization of the weld-pool metal and the diffusion of the metal vapour into the arc plasma are found to decrease the heat flux density and current density to the weld pool. Finally, we show that the shape of the wire electrode in metal-inert-gas welding has a strong influence on flow velocities in the arc and the pressure and shear stress at the weld-pool surface. In each case, we present evidence that the geometry and depth of the weld pool depend strongly on the properties of the arc.

Computational fluid dynamics simulation and turbomachinery code validation of a high pressure ratio radial-inflow turbine

The present study explores reproducing the closest geometry of a high pressure ratio single stage radial-inflow turbine applied in the Sundstrans Power Systems T-100 Multipurpose Small Power Unit. The commercial software ANSYS-Vista RTD along with a built in module, BladeGen, is used to conduct a meanline design and create 3D geometry of one flow passage. Carefully examining the proposed design against the geometrical and experimental data, ANSYS-TurboGrid is applied to generate computational mesh. CFD simulations are performed with ANSYS-CFX in which three-dimensional Reynolds-Averaged Navier-Stokes equations are solved subject to appropriate boundary conditions. Results are compared with numerical and experimental data published in the literature in order to generate the exact geometry of the existing turbine and validate the numerical results against the experimental ones.

CFD simulations of flow and heat transfer through the porous interface of a metal foam heat exchanger

This paper offers numerical modelling of a waste heat recovery system. A thin layer of metal foam is attached to a cold plate to absorb heat from hot gases leaving the system. The heat transferred from the exhaust gas is then transferred to a cold liquid flowing in a secondary loop. Two different foam PPI (Pores Per Inch) values are examined over a range of fluid velocities. Numerical results are then compared to both experimental data and theoretical results available in the literature. Challenges in getting the simulation results to match those of the experiments are addressed and discussed in detail. In particular, interface boundary conditions specified between a porous layer and a fluid layer are investigated. While physically one expects much lower fluid velocity in the pores compared to that of free flow, capturing this sharp gradient at the interface can add to the difficulties of numerical simulation. The existing models in the literature are modified by considering the pressure gradient inside and outside the foam. Comparisons against the numerical modelling are presented. Finally, based on experimentally-validated numerical results, thermo-hydraulic performance of foam heat exchangers as waste heat recovery units is discussed with the main goal of reducing the excess pressure drop and maximising the amount of heat that can be recovered from the hot gas stream.

On the effect of topology on cellular automata rule spaces

This thesis presents an empirical study of the effects of topology on cellular automata rule spaces. The classical definition of a cellular automaton is restricted to that of a regular lattice, often with periodic boundary conditions. This definition is extended to allow for arbitrary topologies. The dynamics of cellular automata within the triangular tessellation were analysed when transformed to 2-manifolds of topological genus 0, genus 1 and genus 2. Cellular automata dynamics were analysed from a statistical mechanics perspective. The sample sizes required to obtain accurate entropy calculations were determined by an entropy error analysis which observed the error in the computed entropy against increasing sample sizes. Each cellular automata rule space was sampled repeatedly and the selected cellular automata were simulated over many thousands of trials for each topology. This resulted in an entropy distribution for each rule space. The computed entropy distributions are indicative of the cellular automata dynamical class distribution. Through the comparison of these dynamical class distributions using the E-statistic, it was identified that such topological changes cause these distributions to alter. This is a significant result which implies that both global structure and local dynamics play a important role in defining long term behaviour of cellular automata.

A numerical method for the fractional Fitzhugh–Nagumo monodomain model

A fractional FitzHugh–Nagumo monodomain model with zero Dirichlet boundary conditions is presented, generalising the standard monodomain model that describes the propagation of the electrical potential in heterogeneous cardiac tissue. The model consists of a coupled fractional Riesz space nonlinear reaction-diffusion model and a system of ordinary differential equations, describing the ionic fluxes as a function of the membrane potential. We solve this model by decoupling the space-fractional partial differential equation and the system of ordinary differential equations at each time step. Thus, this means treating the fractional Riesz space nonlinear reaction-diffusion model as if the nonlinear source term is only locally Lipschitz. The fractional Riesz space nonlinear reaction-diffusion model is solved using an implicit numerical method with the shifted Grunwald–Letnikov approximation, and the stability and convergence are discussed in detail in the context of the local Lipschitz property. Some numerical examples are given to show the consistency of our computational approach.

Comparison of eight published lumbar spine finite element models

Due to its ability to represent intricate systems with material nonlinearities as well as irregular loading, boundary, geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in spinal biomechanics. Current FE models generally account for a single distinct spinal geometry with one set of material properties despite inherently large inter-subject variability. The uncertainty and high variability in tissue material properties, geometry, loading and boundary conditions has cast doubt on the reliability of their predictions and comparability with reported in vitro and in vivo values. A multicenter study was undertaken to compare the results of eight well-established models of the lumbar spine that have been developed, validated and applied for many years. Models were subjected to pure and combined loading modes and their predictions were compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges; their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with previously published median in vitro values. However, the ranges of predictions were larger and exceeded the in vitro ranges, especially for facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with in vivo values. The simulations yielded median facet joint forces of 0 N in flexion, 38 N in extension, 14 N in lateral bending and 60 N in axial rotation that could not be validated due to the paucity of in vivo facet joint forces. In light of high inter-subject variability, one must be cautious when generalizing predictions obtained from one deterministic model. This study demonstrates however that the predictive power increases when FE models are combined together. The median of individual numerical results can hence be used as an improved tool in order to estimate the response of the lumbar spine.

Numerical analysis for the time distributed-order and Riesz space fractional diffusions on bounded domains

Subdiffusion equations with distributed-order fractional derivatives describe some important physical phenomena. In this paper, we consider the time distributed-order and Riesz space fractional diffusions on bounded domains with Dirichlet boundary conditions. Here, the time derivative is defined as the distributed-order fractional derivative in the Caputo sense, and the space derivative is defined as the Riesz fractional derivative. First, we discretize the integral term in the time distributed-order and Riesz space fractional diffusions using numerical approximation. Then the given equation can be written as a multi-term time–space fractional diffusion. Secondly, we propose an implicit difference method for the multi-term time–space fractional diffusion. Thirdly, using mathematical induction, we prove the implicit difference method is unconditionally stable and convergent. Also, the solvability for our method is discussed. Finally, two numerical examples are given to show that the numerical results are in good agreement with our theoretical analysis.

Fractional models in space for diffusive processes in heterogeneous media with applications in cell motility and electrical signal propagation

This work addresses fundamental issues in the mathematical modelling of the diffusive motion of particles in biological and physiological settings. New mathematical results are proved and implemented in computer models for the colonisation of the embryonic gut by neural cells and the propagation of electrical waves in the heart, offering new insights into the relationships between structure and function. In particular, the thesis focuses on the use of non-local differential operators of non-integer order to capture the main features of diffusion processes occurring in complex spatial structures characterised by high levels of heterogeneity.

Space Station Design Report 2007

Design Proposal for the Blue Lunar Support Hub The conceptual design of a space station is one of the most challenging tasks in aerospace engineering. The history of the space station Mir and the assembly of the International Space Station demonstrate that even within the assembly phase quick solutions have to be found to cope with budget and technical problems or changing objectives. This report is the outcome of the conceptual design of the Space Station Design Workshop (SSDW) 2007, which took place as an international design project from the 16th to the 21st of July 2007 at the Australian Centre for Field Robotics (ACFR), University of Sydney, Australia. The participants were tasked to design a human-tended space station in low lunar orbit (LLO) focusing on supporting future missions to the moon in a programmatic context of space exploration beyond low Earth orbit (LEO). The design included incorporating elements from systems engineering to interior architecture. The customised, intuitive, rapid-turnaround software tools enabled the team to successfully tackle the complex problem of conceptual design of crewed space systems. A strong emphasis was put on improving the integration of the human crew, as it is the major contributor to mission success, while always respecting the boundary conditions imposed by the challenging environment of space. This report documents the methodology, tools and outcomes of the Space Station Design Workshop during the SSDW 2007. The design results produced by Team Blue are presented.

Model updating of real structures with ambient vibration data

It is important to develop reliable finite element models for real structures not only in the design phase but also for the structural health monitoring and structural maintenance purposes. This paper describes the experience of the authors in using ambient vibration model identification techniques together with model updating tools to develop reliable finite element models of real civil engineering structures. Case studies of two real structures are presented in this paper. One is a 10 storey concrete building which is considered as a non-slender structure with complex boundary conditions. The other is a single span concrete foot bridge which is also a relatively inflexible planar structure with complex boundary conditions. Both structures are located at the Queensland University of Technology (QUT) and equipped with continuous structural health monitoring systems.

Real-time traffic state estimation in urban corridors from heterogeneous data

In recent years, rapid advances in information technology have led to various data collection systems which are enriching the sources of empirical data for use in transport systems. Currently, traffic data are collected through various sensors including loop detectors, probe vehicles, cell-phones, Bluetooth, video cameras, remote sensing and public transport smart cards. It has been argued that combining the complementary information from multiple sources will generally result in better accuracy, increased robustness and reduced ambiguity. Despite the fact that there have been substantial advances in data assimilation techniques to reconstruct and predict the traffic state from multiple data sources, such methods are generally data-driven and do not fully utilize the power of traffic models. Furthermore, the existing methods are still limited to freeway networks and are not yet applicable in the urban context due to the enhanced complexity of the flow behavior. The main traffic phenomena on urban links are generally caused by the boundary conditions at intersections, un-signalized or signalized, at which the switching of the traffic lights and the turning maneuvers of the road users lead to shock-wave phenomena that propagate upstream of the intersections. This paper develops a new model-based methodology to build up a real-time traffic prediction model for arterial corridors using data from multiple sources, particularly from loop detectors and partial observations from Bluetooth and GPS devices.

Bending and torsion of hollow flange channel beams

The LiteSteel beam (LSB) is a cold-formed high strength steel channel section made of two torsionally rigid closed flanges and a slender web. Due to its mono-symmetric characteristics, its centroid and shear centre do not coincide. The LSBs can be used in floor systems as joists or bearers and in these applications they are often subjected to transverse loads that are applied away from the shear centre. Hence they are often subjected to combined bending and torsion actions. Previous researches on LSBs have concentrated on their bending or shear behaviour and strengths, and only limited research has been undertaken on their combined bending and torsion behaviour. Therefore in this research a series of nine experiments was first conducted on LSBs subject to combined bending and torsion. Three LSB sections were tested to failure under eccentric loading at mid-span, and appropriate results were obtained from seven tests. A special test rig was used to simulate two different eccentricities and to provide accurate simple boundary conditions at the supports. Finite element models of tested LSBs were developed using ANSYS, and the ultimate strengths, failure modes, and load–displacement curves were obtained and compared with corresponding test results. Finite element analyses agreed well with test results and hence the developed models were used in a parametric study to investigate the effects of load locations, eccentricities, and spans on the combined bending and torsion behaviour of LSBs. The interaction between the ultimate bending and torsional moment capacities was studied and a simple design rule was proposed. This paper presents the details of the tests, finite element analyses, and parametric study of LSBs subject to combined bending and torsion, and the results.

Development of a new design expression for the in-plane shear capacity of the partially grouted concrete masonry walls

Partially grouted masonry walls subjected to in-plane shear exhibit a complex behaviour because of the influence of the aspect ratio, the pre-compression, the grouting pattern, the ratios of the horizontal and the vertical reinforcements, the boundary conditions and the characteristics of the constituent materials. The existing in-plane shear expressions for the partially grouted masonry are formulated as sum of strength of three parameters, namely, the masonry, the reinforcement and the axial force. The parameter ‘masonry’ includes the wall aspect ratio and the masonry compressive strength; the aspect ratio of the unreinforced panel inscribed into the grouted cores and bond beams are not considered, although failure is often dominated by these unreinforced masonry panels. This paper describes the dominance of these panels, particularly those that are squat, to the shear capacity of whole of shear walls. Further, the current design formulae are shown highly un-conservative by many researchers; this paper provides a potential reason for this un-conservativeness. It is shown that by including an additional term of the unreinforced panel aspect ratio a rational design formula could be established. This new expression is validated with independent test results reported in the literature – both Australian and overseas; the predictions are shown to be conservative.

Micro paddy lysimeter for monitoring solute transport in paddy environment

The micro paddy lysimeter (MPL) was developed and evaluated for its performance to simulate solute transport in paddy environment under laboratory conditions. MPLs were constructed using soil collected from Field Museum Honmachi of Tokyo University of Agriculture and Technology, Japan. For the physical characteristics of the hardpan layer, parameters such as thickness, and soil aggregate size, affecting the percolation rate were studied. For the plow layer, two types of plow soils, sieved and un-sieved soils were compared. The sieved soil plow layer was produced by mixing air-dried soils of different aggregate sizes of D > 9.50, 9.50 ≥ D > 4.75, 4.75 ≥ D > 2.0 mm and D ≤ 2.0 mm at 47.1, 19.5, 20.6, and 12.8%, respectively. The un-sieved plow layer soil was directly used after collecting from the field. Inert tracer was applied to ponding water with controlled boundary conditions to evaluate the reproducibility of the soil hydraulic characteristics. HYDRUS-1D was used to evaluate the movement of bromide tracer in the MPL. The proposed conditions of the MPL were that the hardpan layer can be made from soil aggregates smaller than 0.425 mm with 2 cm thickness and that the plow layer can be prepared with sieved or un-sieved soils. With these conditions, the obtained results proved that MPLs can be a useful tool to simulate solute transport in paddy environment.

Modelling of the evaporation of a droplet suspended in a binary atmosphere

The process of spray drying is applied in a number of contexts. One such application is the production of a synthetic rock used for storage of nuclear waste. To establish a framework for a model of the spray drying process for this application, we here develop a model describing evaporation from droplets of pure water, such that the model may be extended to account for the presence of colloid within the droplet. We develop a spherically-symmetric model and formulate continuum equations describing mass, momentum, and energy balance in both the liquid and gas phases from first principles. We establish appropriate boundary conditions at the surface of the droplet, including a generalised Clapeyron equation that accurately describes the temperature at the surface of the droplet. To account for experiment design, we introduce a simplified platinum ball and wire model into the system using a thin wire problem. The resulting system of equations is transformed in order to simplify a finite volume solution scheme. The results from numerical simulation are compared with data collected for validation, and the sensitivity of the model to variations in key parameters, and to the use of Clausius–Clapeyron and generalised Clapeyron equations, is investigated. Good agreement is found between the model and experimental data, despite the simplicity of the platinum phase model.