Modelling trickle irrigation: Comparison of analytical and numerical models for estimation of wetting front position with time

Irrigation practices that are profligate in their use of water have come under closer scrutiny by water managers and the public. Trickle irrigation has the propensity to increase water use efficiency but only if the system is designed to meet the soil and plant conditions. Recently we have provided a software tool, WetUp (http://www.clw.csiro.au/products/wetup/), to calculate the wetting patterns from trickle irrigation emitters. WetUp uses an analytical solution to calculate the wetted perimeter for both buried and surface emitters. This analytical solution has a number of assumptions, two of which are that the wetting front is defined by water content at which the hydraulic conductivity (K) is I mm day(-1) and that the flow occurs from a point source. Here we compare the wetting patterns calculated with a 2-dimensional numerical model, HYDRUS2D, for solving the water flow into typical soils with the analytical solution. The results show that the wetting patterns are similar, except when the soil properties result in the assumption of a point source no longer being a good description of the flow regime. Difficulties were also experienced with getting stable solutions with HYDRUS2D for soils with low hydraulic conductivities. (c) 2005 Elsevier Ltd. All rights reserved.

Stirring in 3-d spherical models of convection in the Earth's mantle

On a global scale basalts from mid-ocean ridges are strikingly more homogeneous than basalts from intraplate volcanism. The observed geochemical heterogeneity argues strongly for the existence of distinct reservoirs in the Earth's mantle. It is an unresolved problem of Geodynamics as to how these findings can be reconciled with large-scale convection. We review observational constraints, and investigate stirring properties of numerical models of mantle convection. Conditions in the early Earth may have supported layered convection with rapid stirring in the upper layers. Material that has been altered near the surface is transported downwards by small-scale convection. Thereby a layer of homogeneous depleted material develops above pristine mantle. As the mantle cools over Earth history, the effects leading to layering become reduced and models show the large-scale convection favoured for the Earth today. Laterally averaged, the upper mantle below the lithosphere is least affected by material that has experienced near-surface differentiation. The geochemical signature obtained during the previous episode of small-scale convection may be preserved there for the longest time. Additionally, stirring is less effective in the high viscosity layer of the central lower mantle [1, 2], supporting the survival of medium-scale heterogeneities there. These models are the first, using 3-d spherical geometry and mostly Earth-like parameters, to address the suggested change of convective style. Although the models are still far from reproducing our planet, we find that proposal might be helpful towards reconciling geochemical and geophysical constraints.

Pumice rafting and faunal dispersion during 2001-2002 in the Southwest Pacific: record of a dacitic submarine explosive eruption from Tonga

A new influx of sea-rafted pumice reached the eastern coast of Australia in October 2002, approximately 1 year after a felsic, shallow-marine explosive eruption at a previously unknown volcano (0403-091) along the Tofua volcanic arc (Tonga). The eruption produced floating pumice rafts that first became stranded in Fiji in November 2001, approximately I month after the eruption. Strandings of sea-rafted pumice along shorelines have been the only record of products from this submarine explosive eruption at the remote, submerged volcano. Computed drift trajectories of the sea-rafted pumice using numerical models of southwest Pacific surface wind fields and ocean currents indicate two cyclonic systems disturbed the drift of pumice to eastern Australia, as well as the importance of the combined wave and direct wind effect on pumice trajectory. Pumice became stranded along at least two-thirds (>2000 km) of the coastline of eastern Australia, being deposited on beaches during a sustained period of fresh onshore winds. Typical amounts of pumice initially stranded on beaches were 500-4000 individual clasts per in, and a minimum volume estimate of pumice that arrived to eastern Australia is 1.25 x 10(5) m(3). Pumice was beached below maximum tidal/storm surge levels and was quickly reworked back into the ocean, such that the concentration of beached pumice rapidly dissipated within weeks of the initial stranding, and little record of this stranding event now exists. Most stranded pumice clasts ranged in size from 2 to 5 cm in diameter; the largest measured clasts were 10 cm in Australia and 20 cm in Fiji. The pumice has a low phenocryst content (3500 km) and period of pumice floatation (greater than or equal to1 year), confirm the importance of sea-rafted pumice as a long-distance dispersal mechanism for marine organisms including marine pests and harmful invasive species. Billions of individual rafting pumice clasts can be generated in a single small-volume eruption, such as observed here, and the geological implications for the transport of sessile taxa over large distances are significant. An avenue for future research is to examine whether speciation events and volcanicity are linked; the periodic development of globalism for some taxa (e.g., corals, gastropods, bryozoa) may correlate in time and/or space with voluminous silicic igneous events capable of producing >10(6) km(3) of silicic pumice-rich pyroclastic material and emplaced into ocean basins. (C) 2004 Elsevier B.V. All rights reserved.

Mine waste beach profile and flow resistance equations

Four mine waste beach longitudinal profile equations are compared theoretically and in statistical analyses of profile data from 64 field and laboratory beaches formed by mine tailings, co-disposed coal mine wastes, and sand. All four equations fit the profile data well. The best performing equation both accounts for particle sorting and satisfies hydraulic constraints, and the combination of assumptions underlying it is considered to best represent the processes occurring on mine waste beaches. Combining these assumptions with the Lacey normal equation leads to a variant of the Manning resistance equation. Features that it is desirable to incorporate in theoretical and numerical models of mine waste beaches are listed.

Large scale experiments on gravel and mixed beaches: Experimental procedure, data documentation and initial results

This paper provides information on the experimental set-up, data collection methods and results to date for the project Large scale modelling of coarse grained beaches, undertaken at the Large Wave Channel (GWK) of FZK in Hannover by an international group of researchers in Spring 2002. The main objective of the experiments was to provide full scale measurements of cross-shore processes on gravel and mixed beaches for the verification and further development of cross-shore numerical models of gravel and mixed sediment beaches. Identical random and regular wave tests were undertaken for a gravel beach and a mixed sand/gravel beach set up in the flume. Measurements included profile development, water surface elevation along the flume, internal pressures in the swash zone, piezometric head levels within the beach, run-up, flow velocities in the surf-zone and sediment size distributions. The purpose of the paper is to present to the scientific community the experimental procedure, a summary of the data collected, some initial results, as well as a brief outline of the on-going research being carried out with the data by different research groups. The experimental data is available to all the scientific community following submission of a statement of objectives, specification of data requirements and an agreement to abide with the GWK and EU protocols. (C) 2005 Elsevier B.V. All rights reserved.

A review of flow modeling for dense medium cyclones

A critical assessment is presented for the existing fluid flow models used for dense medium cyclones (DMCs) and hydrocyclones. As the present discussion indicates, the understanding of dense medium cyclone flow is still far from the complete. However, its similarity to the hydrocyclone provides a basis for improved understanding of fluid flow in DMCs. The complexity of fluid flow in DMCs is basically due to the existence of medium as well as the dominance of turbulent particle size and density effects on separation. Both the theoretical and experimental analysis is done with respect to two-phase motions and solid phase flow in hydrocyclones or DMCs. A detailed discussion is presented on the empirical, semiempirical, and the numerical models based upon both the vorticity-stream function approach and Navier-Stokes equations in their primitive variables and in cylindrical coordinates available in literature. The existing equations describing turbulence and multiphase flows in cyclone are also critically reviewed.

Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.

Regional-scale, fully coupled modelling of stream-aquifer interaction in a tropical catchment

The planning and management of water resources in the Pioneer Valley, north-eastern Australia requires a tool for assessing the impact of groundwater and stream abstractions on water supply reliabilities and environmental flows in Sandy Creek (the main surface water system studied). Consequently, a fully coupled stream-aquifer model has been constructed using the code MODHMS, calibrated to near-stream observations of watertable behaviour and multiple components of gauged stream flow. This model has been tested using other methods of estimation, including stream depletion analysis and radon isotope tracer sampling. The coarseness of spatial discretisation, which is required for practical reasons of computational efficiency, limits the model's capacity to simulate small-scale processes (e.g., near-stream groundwater pumping, bank storage effects), and alternative approaches are required to complement the model's range of applicability. Model predictions of groundwater influx to Sandy Creek are compared with baseflow estimates from three different hydrograph separation techniques, which were found to be unable to reflect the dynamics of Sandy Creek stream-aquifer interactions. The model was also used to infer changes in the water balance of the system caused by historical land use change. This led to constraints on the recharge distribution which can be implemented to improve model calibration performance. (c) 2006 Elsevier B.V. All rights reserved.

Evaluation of parameters effective in FRP shear strengthening of RC beams using FE method

Shear strengthening is required when an RC beam is found deficient in shear, or when its shear capacity falls below its flexural capacity after flexural strengthening. A recent technique for the shear strengthening of RC beams is to provide additional FRP web reinforcement, commonly in the form of bonded external FRP strips/sheets. Over the last few years, several experimental studies have been conducted on this new strengthening technique, which has established its effectiveness. While experimental methods of investigation are extremely useful in obtaining information about the composite behaviour of FRP and reinforced concrete, the use of numerical models such as the one presented in this paper helps in developing a good understanding of the behaviour at lower costs. In the study presented in this paper, ANSYS finite element program is used to examine the response of beams strengthened in shear by FRPs. The FE model is calibrated against test results performed at the University of Kentucky. Once validated, the model is used to examine the influence of fibre orientation, compressive strength of concrete, area of tensile and compressive reinforcements, and amount and distance between stirrups on the strength and ductility of FRP strengthened beam.

Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

Numerical modelling of radiating superorbital flows

Experimental aerodynamic studies of the flows around new aerocapture spacecraft configurations are presently being done in the superorbital expansion tubes at The University of Queensland. Short duration flows at speeds of 10--13 km/s are produced in the expansion tube facility and are then applied to the model spacecraft. Although high-temperature effects, such as molecular dissociation, have long been a part of the computational modelling of the expansion tube flows for speeds below 10 km/s, radiation may now be a significant mechanism of energy transfer within the shock layer on the model. This paper will study the coupling of radiation energy transport for an optically thin gas to the flow dynamics in order to obtain accurate predictions of thermal loads on the spacecraft. The results show that the effect of radiation on the flowfields of subscale models for expansion tube experiments can be assessed by measurements of total heat transfer and radiative heat transfer.