Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Investigation on mechanical behaviour of AM60 magnesium alloys

Purpose: In this work, tension, impact, bend and fatigue tests were conducted in an AM60 magnesium alloy. The effects of environmental temperature and loading rates on impact and tension behavior of the alloy were also investigated. Design/methodology/approach: The tests were conducted using an Instron universal testing machine. The loading speed was changed from 1 mm/min to 300 mm/min to gain a better understanding of the effect of strain rate. To understand the failure behavior of this alloy at different environmental temperatures, Charpy impact test was conducted in a range of temperatures (-40~35°C). Plane strain fracture toughness (KIC) was evaluated using compact tension (CT) specimen. To gain a better understanding of the failure mechanisms, all fracture surfaces were observed using scanning electron microscopy (SEM). In addition, fatigue behavior of this alloy was estimated using tension test under tension-tension condition at 30 Hz. The stress amplitude was selected in the range of 20~50 MPa to obtain the S-N curve. Findings: The tensile test indicated that the mechanical properties were not sensitive to the strain rates applied (3.3x10-4~0.1) and the plastic deformation was dominated by twining mediated slip. The impact energy is not sensitive to the environmental temperature. The plane strain fracture toughness and fatigue limit were evaluated and the average values were 7.6 MPa.m1/2 and 25 MPa, respectively. Practical implications: Tested materials AM60 Mg alloy can be applied among others in automotive industry aerospace, communication and computer industry. Originality/value: Many investigations have been conducted to develop new Mg alloys with improved stiffness and ductility. On the other hand, relatively less attention has been paid to the failure mechanisms of Mg alloys, such as brittle fracture and fatigue, subjected to different environmental or loading conditions. In this work, tension, impact, bend and fatigue tests were conducted in an AM60 magnesium alloy.

Understanding chromosome disorders and their implications for special educators

More children are now being diagnosed with chromosome abnormalities. Some chromosome disorder syndromes are relatively well known; while others are so rare that there is only limited evidence about their likely impact on learning and development. For educators, a basic level of knowledge about chromosome abnormalities is important for understanding the literature and communicating with families and professionals. This paper describes chromosomes, and the numerical and structural anomalies that can occur, usually spontaneously during early cell division. Distinctive features of various chromosome syndromes are summarised before a discussion of the rare chromosome disorders that are labelled, not with a syndrome name, but simply by a description of the chromosome number, size and shape. Because of the potential within-group variability that characterises syndromes, and the scarcity of literature about the rare chromosome disorders, expectations for learning and development of individual students need to be based on the range of possible outcomes that may be achievable.

Numerical investigation of the influence of topography on simulated downburst wind fields

A, dry, non-hydrostatic sub-cloud model is used to simulate an isolated stationary downburst wind event to study the influence topographic features have on the near-ground wind structure of these storms. It was generally found that storm maximum wind speeds could be increased by up to 30% because of the presence of a topographic feature at the location of maximum wind speeds. Comparing predicted velocity profile amplification with that of a steady flow impinging jet, similar results were found despite the simplifications made in the impinging jet model. Comparison of these amplification profiles with those found in the simulated boundary layer winds reveal reductions of up to 30% in the downburst cases. Downburst and boundary layer amplification profiles were shown to become more similar as the topographic feature height was reduced with respect to the outflow depth.

Numerical simulation of downburst winds

The wind field of an intense idealised downburst wind storm has been studied using an axisymmetric, dry, non-hydrostatic numerical sub-cloud model. The downburst driving processes of evaporation and melting have been paramaterized by an imposed cooling source that triggers and sustains a downdraft. The simulated downburst exhibits many characteristics of observed full-scale downburst events, in particular the presence of a primary and counter rotating secondary ring vortex at the leading edge of the diverging front. The counter-rotating vortex is shown to significantly influence the development and structure of the outflow. Numerical forcing and environmental characteristics have been systematically varied to determine the influence on the outflow wind field. Normalised wind structure at the time of peak outflow intensity was generally shown to remain constant for all simulations. Enveloped velocity profiles considering the velocity structure throughout the entire storm event show much more scatter. Assessing the available kinetic energy within each simulated storm event, it is shown that the simulated downburst wind events had significantly less energy available for loading isolated structures when compared with atmospheric boundary layer winds. The discrepancy is shown to be particularly prevalent when wind speeds were integrated over heights representative of tall buildings. A similar analysis for available full scale measurements led to similar findings.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Numerical modelling and design of lipped channel beams subject to shear

Cold-formed steel members are increasingly used as primary structural elements in buildings due to the availability of thin and high strength steels and advanced cold-forming technologies. Cold-formed lipped channel beams (LCB) are commonly used as flexural members such as floor joists and bearers. Many research studies have been carried out to evaluate the behaviour and design of LCBs subject to pure bending actions. However, limited research has been undertaken on the shear behaviour and strength of LCBs. Hence a numerical study was undertaken to investigate the shear behaviour and strength of LCBs. Finite element models of simply supported LCBs with aspect ratios of 1.0 and 1.5 were considered under a mid-span load. They were then validated by comparing their results with test results and used in a detailed parametric study based on the validated finite element models. Numerical studies were conducted to investigate the shear buckling and post-buckling behaviour of LCBs. Experimental and numerical results showed that the current design rules in cold-formed steel structures design codes are very conservative for the shear design of LCBs. Improved design equations were therefore proposed for the shear strength of LCBs. This paper presents the details of this numerical study of LCBs and the results.

Numerical studies and design of hollow flange channel beams subject to combined bending and shear actions

LiteSteel beam (LSB) is a cold-formed steel hollow flange channel section produced using a patented manufacturing process involving simultaneous cold-forming and dual electric resistance welding. It is commonly used as floor joists and bearers in residential, industrial and commercial buildings. Design of the LSB is governed by the Australian cold-formed steel structures code, AS/NZS 4600. Due to the geometry of the LSB, as well as its unique residual stress characteristics and initial geometric imperfections resultant of manufacturing processes, currently available design equations for common cold-formed sections are not directly applicable to the LSB. Many research studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions and predominant shear actions. To date, however, no investigation has been conducted into the strength of LSB sections under combined bending and shear actions. Hence experimental and numerical studies were conducted to assess the combined bending and shear behaviour of LSBs. Finite element models of LSBs were developed to simulate their combined bending and shear behaviour and strength of LSBs. They were then validated by comparing the results with available experimental test results and used in a detailed parametric study. The results from experimental and finite element analyses were compared with current AS/NZS 4600 and AS 4100 design rules. Both experimental and numerical studies show that the AS/NZS 4600 design rule based on circular interaction equation is conservative in predicting the combined bending and shear capacities of LSBs. This paper presents the details of the numerical studies of LSBs and the results. In response to the inadequacies of current approaches to designing LSBs for combined bending and shear, two lower bound design equations are proposed in this paper.

Numerical estimation of 3D mechanical forces exerted by cells on non-linear materials

The exchange of physical forces in both cell-cell and cell-matrix interactions play a significant role in a variety of physiological and pathological processes, such as cell migration, cancer metastasis, inflammation and wound healing. Therefore, great interest exists in accurately quantifying the forces that cells exert on their substrate during migration. Traction Force Microscopy (TFM) is the most widely used method for measuring cell traction forces. Several mathematical techniques have been developed to estimate forces from TFM experiments. However, certain simplifications are commonly assumed, such as linear elasticity of the materials and/or free geometries, which in some cases may lead to inaccurate results. Here, cellular forces are numerically estimated by solving a minimization problem that combines multiple non-linear FEM solutions. Our simulations, free from constraints on the geometrical and the mechanical conditions, show that forces are predicted with higher accuracy than when using the standard approaches.

Numerical investigation of mechanical and thermal dynamic properties of the industrial transformer

Industrial transformer is one of the most critical assets in the power and heavy industry. Failures of transformers can cause enormous losses. The poor joints of the electrical circuit on transformers can cause overheating and results in stress concentration on the structure which is the major cause of catastrophic failure. Few researches have been focused on the mechanical properties of industrial transformers under overheating thermal conditions. In this paper, both mechanical and thermal properties of industrial transformers are jointly investigated using Finite Element Analysis (FEA). Dynamic response analysis is conducted on a modified transformer FEA model, and the computational results are compared with experimental results from literature to validate this simulation model. Based on the FEA model, thermal stress is calculated under different temperature conditions. These analysis results can provide insights to the understanding of the failure of transformers due to overheating, therefore are significant to assess winding fault, especially to the manufacturing and maintenance of large transformers.

Experimental and numerical analysis of the relationship between indoor and outdoor airborne particles in an operating room

This work investigated the impact of the HVAC filtration system and indoor particle sources on the relationship between indoor and outdoor airborne particle size and concentrations in an operating room. Filters with efficiency between 65% and 99.97% were used in the investigation and indoor and outdoor particle size and concentrations were measured. A balance mass model was used for the simulation of the impact of the surgical team, deposition rate, HVAC exhaust and air change rates on indoor particle concentration. The experimental results showed that high efficiency filters would not be expected to decrease the risk associated with indoor particles larger than approximately 1 µm in size because normal filters are relatively efficient for these large particles. A good fraction of outdoor particles were removed by deposition on the HVAC system surfaces and this deposition increased with particle size. For particles of 0.3-0.5 µm in diameter, particle reduction was about 23%, while for particles >10 µm the loss was about 78%. The modelling results showed that depending on the type of filter used, the surgical team generated between 93-99% of total particles, while the outdoor air contributed only 1-6%.

Experimental and numerical characteristics of portable water-filled road safety barrier system under different impact conditions

This research has developed an innovative road safety barrier system that will enhance roadside safety. In doing so, the research developed new knowledge in the field of road crash mitigation for high speed vehicle impact involving plastic road safety barriers. This road safety barrier system has the required feature to redirecting an errant vehicle with limited lateral displacement. Research was carried out using dynamic computer simulation technique support by experimental testing. Future road safety barrier designers may use the information in this research as a design guideline to improve the performance and redirectional capability of the road safety barrier system. This will lead to better safety conditions on the roadways and potentially save lives.

Unsteady natural convection within a porous enclosure of sinusoidal corrugated side walls

Numerically investigation of free convection within a porous cavity with differential heating has been performed using modified corrugated side walls. Sinusoidal hot left and cold right walls are assumed to receive sudden differentially heating where top and bottom walls are insulated. Air is considered as working fluid and is quiescent, initially. Numerical experiments reveal 3 distinct stages of developing pattern including initial stage, oscillatory intermediate and finally steady state condition. Implicit Finite Volume Method with TDMA solver is used to solve the governing equations. This study has been performed for the Rayleigh numbers ranging from 100 to 10,000. Outcomes have been reported in terms of isotherms, streamline, velocity and temperature plots and average Nusselt number for various Ra, corrugation frequency and corrugation amplitude. The effects of sudden differential heating and its resultant transient behavior on fluid flow and heat transfer characteristics have been shown for the range of governing parameters. The present results show that the transient phenomena are enormously influenced by the variation of the Rayleigh Number with corrugation amplitude and frequency.

Numerical investigation on stable crack growth in plane stress

Large deformation finite element analysis has been carried out to investigate the stress-strain fields ahead of a growing crack for compact tension .a=W D 0:5/ and three-point bend .a=W D 0:1 and 0:5/ specimens under plane stress condition. The crack growth is controlled by the experimental J -integral resistance curves measured by Sun et al. The results indicate that the distributions of opening stress, equivalent stress and equivalent strain ahead of a growing crack are not sensitive to specimen geometry. For both stationary and growing cracks, similar distributions of opening stress and triaxiality can be found along the ligament. During stable crack growth, the crack-tip opening displacement (CTOD) resistance curve and the cohesive fracture energy in the fracture process zone are independent of specimen geometry and may be suitable criteria for characterizing stable crack growth in plane stress.

Combined analysis of categorical and numerical descriptors of Australian groundnut accessions using nonlinear principal component analysis

For users of germplasm collections, the purpose of measuring characterization and evaluation descriptors, and subsequently using statistical methodology to summarize the data, is not only to interpret the relationships between the descriptors, but also to characterize the differences and similarities between accessions in relation to their phenotypic variability for each of the measured descriptors. The set of descriptors for the accessions of most germplasm collections consists of both numerical and categorical descriptors. This poses problems for a combined analysis of all descriptors because few statistical techniques deal with mixtures of measurement types. In this article, nonlinear principal component analysis was used to analyze the descriptors of the accessions in the Australian groundnut collection. It was demonstrated that the nonlinear variant of ordinary principal component analysis is an appropriate analytical tool because subspecies and botanical varieties could be identified on the basis of the analysis and characterized in terms of all descriptors. Moreover, outlying accessions could be easily spotted and their characteristics established. The statistical results and their interpretations provide users with a more efficient way to identify accessions of potential relevance for their plant improvement programs and encourage and improve the usefulness and utilization of germplasm collections.