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.

Web crippling tests of hollow flange channel beams – end one flange and interior one flange load cases

LiteSteel beam (LSB) is a new cold-formed steel hollow flange channel section produced using a patented manufacturing process involving simultaneous cold-forming and dual electric resistance welding. The LSBs were commonly used as floor joists and bearers with web openings in residential, industrial and commercial buildings. Due to the unique geometry of LSBs, as well as its unique residual stress characteristics and initial geometric imperfections resultant of manufacturing processes, much of the existing research for common cold-formed steel sections is not directly applicable to LSBs. Many research studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions, predominant shear and combined actions. However, to date, no investigation has been conducted into the web crippling behaviour and strength of LSB sections. Hence detailed experimental studies were conducted to investigate the web crippling behaviour and strengths of LSBs under EOF (End One Flange) and IOF (Interior One Flange) load cases. A total of 26 web crippling tests was conducted and the results were compared with current AS/NZS 4600 design rules. This comparison showed that AS/NZS 4600 (SA, 2005) design rules are very conservative for LSB sections under EOF and IOF load cases. Suitable design equations have been proposed to determine the web crippling capacity of LSBs based on experimental results. This paper presents the details of this experimental study on the web crippling behaviour and strengths of LiteSteel beams under EOF and IOF load cases.

Web crippling tests of hollow flange channel beams : ETF AND ITF load cases

This paper presents the details of an experimental study of a cold-formed steel hollow flange channel beam known as LiteSteel Beam (LSB) subject to web crippling actions (ETF and ITF). Due to the geometry of the LSB, as well as its unique residual stress characteristics and initial geometric imperfections resultant of manufacturing processes, much of the existing research for common cold-formed steel sections is not directly applicable to LSB. 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 were conducted to assess the web crippling behaviour and strengths of LSBs. Twenty eight web crippling tests were conducted and the results were compared with the current AS/NZS 4600[1] and AISI S100 [2]design equations. Comparison of the ultimate web crippling capacities from tests showed that AS/NZS 4600[1] and AISI S100 [2] design equations are unconservative for LSB sections under ETF and ITF load cases. Hence new equations were proposed to determine the web crippling capacities of LSBs. Suitable design rules were also developed under the DSM format.

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.

The crystal structure and vibrational spectroscopy of jarosite and alunite minerals

The alunite supergroup of minerals is a large hydroxy-sulfate mineral group, which has seen renewed interest following their discovery on Mars. Numerous reviews exist concerning nomenclature, formation, and natural occurrence of this mineral group. Sulfate minerals in general are widely studied and their vibrational spectra are well characterized. However, no specific review concerning alunite and jarosite spectroscopy and crystal structure has been forthcoming. This review focuses on the controversial aspects of the crystal structure and vibrational spectroscopy of jarosite and alunite minerals. Inconsistencies regarding band assignments especially in the 1000–400 cm−1 region plague these two mineral groups and result in different band assignments among the various spectroscopic studies. There are significant crystallographic and magnetic structure ambiguities with regards to ammonium and hydronium end-members, namely, the geometry these two ions assume in the structure and the fact that hydronium jarosite is a spin glass. It was also found that the synthetic causes for the super cell in plumbojarosite, minamiite, huangite, and walthierite are not known.

Generation of two isomers of C5H from the corresponding anions. A theoretically motivated mass spectrometric study

Molecular orbital calculations have predicted the stability of a range of connectivities for the radical C5H potential surface. The most energetically favorable of these include the linear C4CH geometry and two ring-chain structures HC2C3 and C2C3H The corresponding anions are also shown to be theoretically stable, and furthermore, a fourth isomer, C2CHC2, is predicted to be the most stable anion connectivity. These results have motivated experimental efforts. Methodologies for the generation of the non-ring-containing isomeric anions C4CH and C2CHC2 have been developed utilizing negative ion mass spectrometry. The absolute connectivities of the anions have been established using deuterium labeling, charge reversal, and neutralization reionization techniques. The success of the latter experiment confirms theoretical predictions of stability of the corresponding neutral species. This is the first reported observation of the neutral C2CHC2 species that calculations predict to be substantially less stable than the C4CH connectivity but still bound relative to isomerization processes.

Energy-based positioning control of underactuated vehicles using manifold regulation

In this paper, we consider the problem of position regulation of a class of underactuated rigid-body vehicles that operate within a gravitational field and have fully-actuated attitude. The control objective is to regulate the vehicle position to a manifold of dimension equal to the underactuation degree. We address the problem using Port-Hamiltonian theory, and reduce the associated matching PDEs to a set of algebraic equations using a kinematic identity. The resulting method for control design is constructive. The point within the manifold to which the position is regulated is determined by the action of the potential field and the geometry of the manifold. We illustrate the performance of the controller for an unmanned aerial vehicle with underactuation degree two-a quadrotor helicopter.

Monte Carlo simulations of small field output factors in solid phantoms

Introduction Total scatter factor (or output factor) in megavoltage photon dosimetry is a measure of relative dose relating a certain field size to a reference field size. The use of solid phantoms has been well established for output factor measurements, however to date these phantoms have not been tested with small fields. In this work, we evaluate the water equivalency of a number of solid phantoms for small field output factor measurements using the EGSnrc Monte Carlo code. Methods The following small square field sizes were simulated using BEAMnrc: 5, 6, 7, 8, 10 and 30 mm. Each simulated phantom geometry was created in DOSXYZnrc and consisted of a silicon diode (of length and width 1.5 mm and depth 0.5 mm) submersed in the phantom at a depth of 5 g/cm2. The source-to-detector distance was 100 cm for all simulations. The dose was scored in a single voxel at the location of the diode. Interaction probabilities and radiation transport parameters for each material were created using custom PEGS4 files. Results A comparison of the resultant output factors in the solid phantoms, compared to the same factors in a water phantom are shown in Fig. 1. The statistical uncertainty in each point was less than or equal to 0.4 %. The results in Fig. 1 show that the density of the phantoms affected the output factor results, with higher density materials (such as PMMA) resulting in higher output factors. Additionally, it was also calculated that scaling the depth for equivalent path length had negligible effect on the output factor results at these field sizes. Discussion and conclusions Electron stopping power and photon mass energy absorption change minimally with small field size [1]. Also, it can be seen from Fig. 1 that the difference from water decreases with increasing field size. Therefore, the most likely cause for the observed discrepancies in output factors is differing electron disequilibrium as a function of phantom density. When measuring small field output factors in a solid phantom, it is important that the density is very close to that of water.

Design of light gauge steel roofing systems subject to high wind forces

Light gauge steel roofing systems made of thin profiled roof sheeting and battens are used commonly in residential, industrial and commercial buildings. Their critical design load combination is that due to wind uplift forces that occur during high wind events such as tropical cyclones and thunderstorms. However, premature local failures at their screw connections have been a concern for many decades since cyclone Tracy that devastated Darwin in 1974. Extensive research that followed cyclone Tracy on the pull-through and pull-out failures of roof sheeting to batten connections has significantly improved the safety of roof sheeting. However, this has made the batten to rafter/truss connection the weakest, and recent wind damage investigations have shown the failures of these connections and the resulting loss of entire roof structures. Therefore an experimental research program using both small scale and full scale air-box tests is currently under way to investigate the pull-through failures of thin-walled steel battens under high wind uplift forces. Tests have demonstrated that occurrence of pull-through failures in the bottom flanges of steel batttens and the need to develop simple test and design methods as a function of many critical parameters such as steel batten geometry, thickness and grade, screw fastener sizes and other fastening details. This paper presents the details of local failures that occur in light fauge roofing systems, a review of the current design and test methods for steel battens and associated short comings, and the test results obtained to date on pull-through failures of battens from small scale and full scale tests. Finally, it proposes the use of suitable small scale test methods that can be used by both researchers and manufacturers of such screw-fastened light gauge steel batten systems.

Photoelectron spectroscopy of HCCN- and HCNC- reveals the quasilinear triplet carbenes, HCCN and HCNC

Negative ion photoelectron spectroscopy has been used to study the HCCN- and HCNC- ions. The electron affinities (EA) of cyanocarbene have been measured to be EA(HCCN (X) over tilde (3)Sigma(-)=2.003+/-0.014 eV and EA(DCCN (X) over tilde (3)Sigma(-))=2.009+/-0.020 eV. Photodetachment of HCCN- to HCCN (X) over tilde (3)Sigma(-) shows a 0.4 eV long vibrational progression in nu(5), the H-CCN bending mode; the HCCN- photoelectron spectra reveal excitations up to 10 quanta in nu(5). The term energies for the excited singlet state are found to be T-0(HCCN (a) over tilde (1)A('))=0.515+/-0.016 eV and T-0(DCCN (a) over tilde (1)A('))=0.518+/-0.027 eV. For the isocyanocarbene, the two lowest states switch and HCNC has a singlet ground state and an excited triplet state. The electron affinities are EA(HCNC (X) over tilde (1)A('))=1.883+/-0.013 eV and EA((X) over tilde (1)A(') DCNC)=1.877+/-0.010 eV. The term energy for the excited triplet state is T-0(HCNC (a) over tilde (3)A("))=0.050+/-0.028 eV and T-0(DCNC (a) over tilde (3)A("))=0.063+/-0.030 eV. Proton transfer kinetics in a flowing afterglow apparatus were used to re-measure the enthalpy of deprotonation of CH3NC to be Delta(acid)H(298)(CH3NC)=383.6+/-0.6 kcal mol(-1). The acidity/EA thermodynamic cycle was used to deduce D-0(H-CHCN)=104+/-2 kcal mol(-1) [Delta(f)H(0)(HCCN)=110+/-4 kcal mol(-1)] and D-0(H-CHNC)=106+/-4 kcal mol(-1) [Delta(f)H(0)(HCNC)=133+/-5 kcal mol(-1)]. (C) 2002 American Institute of Physics.