Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Evolution of plant water relations research in Australia

Ecophysiological research in Australia has focussed, at different times, on the fundamental similarities in function between all plant species, and on the peculiarity of Australian species with respect to their survival in stressful environments. Early work on plant water relations emphasised the differences between species, and indicated that diverse structural and functional attributes occurred in species from the same water-limited environment. Most recent research has emphasised processes that optimise rates of carbon dioxide exchange, but the understanding of functioning in plants with different morphological arrangements is incomplete. Variation in functions between individual plants and geographic populations in wild species has been examined to a lesser extent. The great variety within and between populations of wild plant species warrants further study for both understanding and more effective management of this biological resource.

Effects of nitrogen supply on xylem cytokinin delivery, transpiration and leaf expansion of pea genotypes differing in xylem-cytokinin concentration

The rms2 and rms4 pea ( Pisum sativum L.) branching mutants have higher and lower xylem-cytokinin concentration, respectively, relative to wild type (WT) plants. These genotypes were grown at two levels of nitrogen (N) supply for 18 - 20 d to determine whether or not xylem-cytokinin concentration (X-CK) or delivery altered the transpiration and leaf growth responses to N deprivation. Xylem sap was collected by pressurising de-topped root systems. As sap-flow rate increased, X-CK declined in WT and rms2, but did not change in rms4. When grown at 5.0 mM N, X-CKs of rms2 and rms4 were 36% higher and 6-fold lower, respectively, than WT at sap-flow rates equivalent to whole-plant transpiration. Photoperiod cytokinin (CK) delivery rates ( the product of transpiration and X-CK) decreased more than 6-fold in rms4. Growth of plants at 0.5 mM N had negligible (< 10%) effects on transpiration rates expressed on a leaf area basis in WT and rms4, but decreased transpiration rates of rms2. The low-N treatment decreased leaf expansion by 20 - 25% and expanding leaflet N concentration by 15%. These changes were similar in all genotypes. At sap-flow rates equivalent to whole-plant transpiration, the low N treatment decreased X-CK in rms2 but had no discernible effect in WT and rms4. Since the low N treatment decreased transpiration of all genotypes, photoperiod CK delivery rates also decreased in all genotypes. The similar leaf growth response of all genotypes to N deprivation despite differences in both absolute and relative X-CKs and deliveries suggests that shoot N status is more important in regulating leaf expansion than xylem-supplied cytokinins. The decreased X-CK and transpiration rate of rms2 following N deprivation suggests that changes in xylem-supplied CKs may modify water use.

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.

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.

The effect of different levels of dietary crude protein on urea metabolism of rusa deer (Cervus timorensis)

Two experiments were conducted to measure urea recycling and rumen flow dynamics in young rusa deer fed low (LP) or high (HP) protein diets. Pool size and flux rate of labelled urea. into and out of the blood pool were measured using single intravenous (i.v.) injection solutions containing [C-14] - and [N-15]-urea. A curve peeling technique was used to fit the enrichment of N-15 or specific radioactivity (SRA) of C-14 to exponential equations. Body urea-N pool size was significantly greater (P < 0.05) when a HP, compared to a LP diet, was fed. Urea space, expressed as a percent of live weight, total flux rate of urea through the blood pool and the irreversible loss of urea was similar for both diets. The mean (+/- S.E.M.) concentration of plasma urea-N was greater when animals were fed the HP diet compared to the LP diet (2 1.1 +/- 0.3 versus 14.4 +/- 1.4 mg/100 ml, respectively). Voluntary feed intake and digestibility of dietary components were also measured. Daily dry matter intakes were not affected by the crude protein (CP) content of the diet, although apparent DM digestibility was significantly greater for HP diet fed in both experiments. An intraruminal infusion of CrEDTA was used to determine rumen flow dynamics. Ruminal mean retention time, relative net outflow rate of water and passage rate constant (k(w)) were significantly greater (P < 0.05) when the HP diet was fed compared to the LP diet. The extent of urea metabolism and flux rates of urea between the blood and secondary pools appear similar to those previously reported for other ruminants fed diets contrasting in CP content. (c) 2005 Elsevier B.V. All rights reserved.

Tidal impacts on riparian salinities near estuaries

Groundwater-dependent riparian biota is known to be sensitive to changes in soil and groundwater salinity in estuarine systems. The groundwater flow and salinity behaviour in a phreatic aquifer adjoining a partially penetrating, tidal. estuary is investigated through two-dimensional numerical experiments for a lateral cross-section, which explore the influence of factors, such as aquifer and soil materials, tidal amptitudes, and regional groundwater hydraulic gradients. The density contrast between estuarine water and the fresh groundwater drives saltwater penetration of the aquifer even in the case of a marked groundwater hydraulic gradient towards the estuary. We show that tidal fluctuations in estuaries can significantly affect the groundwater salinity distribution in adjacent density-stratified phreatic aquifers. This has consequences for the expected distribution of salinity-sensitive biota in the hyporheic zone as well as vegetation and fauna dependent on water in the riparian soil and aquifer. The shape of the dense saltwater wedge propagating into the adjacent groundwater system is also modified by the estuarine tidal signal, although this effect appears to have only minor influence on the maximum distance penetrated into the aquifer (i.e., location of the 'toe' of the wedge). Tide-induced changes to riparian groundwater salinity are advection-driven, as evidenced by the modified time-averaged groundwater flow dynamics. (c) 2006 Elsevier B.V. All rights reserved.

Prediction of steady-state flux through variably saturated zones within a septic absorption trench

This paper describes effluent flow dynamics within a septic absorption system and the prediction of flow through the biomat and sub-biomat zone. Using soil hydraulic properties in a one dimensional model we demonstrate how soil hydraulic properties interact with biomat resistances to determine long-term acceptance rate (LTAR). The LTAR is a key parameter used in the Australian and New Zealand Standard AS1547:2000 to calculate the area of trench required to ensure trenches are not overloaded. Results show that several orders of magnitude variation in saturated hydraulic conductivity (Ks) collapse to a one order of magnitude variation in LTAR. These results are calculated from a model using basic flow theory, allowing LTAR to be estimated for any combination of biomat resistance and soil hydraulic properties. To increase the reliability of prediction of septic trench hydrology, HYDRUS 2D was used to model two dimensional flow. For more permeable soils, the exfiltration zone above sidewall biomat growth is shown to be a key pathway for excess effluent flow.

Study of flow behaviour in a three-product cyclone using computational fluid dynamics

Simplicity in design and minimal floor space requirements render the hydrocyclone the preferred classifier in mineral processing plants. Empirical models have been developed for design and process optimisation but due to the complexity of the flow behaviour in the hydrocyclone these do not provide information on the internal separation mechanisms. To study the interaction of design variables, the flow behaviour needs to be considered, especially when modelling the new three-product cyclone. Computational fluid dynamics (CFD) was used to model the three-product cyclone, in particular the influence of the dual vortex finder arrangement on flow behaviour. From experimental work performed on the UG2 platinum ore, significant differences in the classification performance of the three-product cyclone were noticed with variations in the inner vortex finder length. Because of this simulations were performed for a range of inner vortex finder lengths. Simulations were also conducted on a conventional hydrocyclone of the same size to enable a direct comparison of the flow behaviour between the two cyclone designs. Significantly, high velocities were observed for the three-product cyclone with an inner vortex finder extended deep into the conical section of the cyclone. CFD studies revealed that in the three-product cyclone, a cylindrical shaped air-core is observed similar to conventional hydrocyclones. A constant diameter air-core was observed throughout the inner vortex finder length, while no air-core was present in the annulus. (c) 2006 Elsevier Ltd. All rights reserved.

Population dynamics and gene flow of Helicoverpa armigera (Lepidoptera : Noctuidae) on cotton and grain crops in the Murrumbidgee Valley, Australia

The population dynamics of Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) in the Murrumbidgee Valley, Australia, has been characterized using five highly variable microsatellite loci. In the 2001-2002 growing season, there were very high levels of migration into the Murrumbidgee Valley with no detectable genetic structuring, consistent with previous analyses on a national scale. By contrast, there was significant genetic structuring over the 2002-2003 growing season, with three distinct genetic types detected. The first type corresponded to the first two generations and was derived from local individuals emerging from diapause and their progeny. The second genetic type corresponded to generation 3 and resulted from substantial immigration into the region. There was another genetic shift in generation 4, which accounts for the third genetic type of the season. This genetic shift occurred despite low levels of immigration. During the third generation of the 2002-2003 growing season, different population dynamics was characterized for H. armigera on maize, Zea mays L., and cotton Gossipium hirsutum L. Populations on cotton tended to cycle independently with very little immigration from outside the region or from maize within the region. Maize acted as a major sink for immigrants from cotton and from outside the region. If resistance were to develop on cotton under these circumstances, susceptible individuals from maize or from other regions would not dilute this resistance. In addition, resistance is likely to be transferred to maize and be perpetuated until diapause, from where it may reemerge next season. If low levels of immigration were to occur on transgenic cotton, this may undermine the effectiveness of refugia, especially noncotton refugia.

Measuring and predicting the dynamics of linear monodisperse entangled polymers in rapid flow through an abrupt contraction. A small angle neutron scattering study

Small-angle neutron scattering measurements on a series of monodisperse linear entangled polystyrene melts in nonlinear flow through an abrupt 4:1 contraction have been made. Clear signatures of melt deformation and subsequent relaxation can be observed in the scattering patterns, which were taken along the centerline. These data are compared with the predictions of a recently derived molecular theory. Two levels of molecular theory are used: a detailed equation describing the evolution of molecular structure over all length scales relevant to the scattering data and a simplified version of the model, which is suitable for finite element computations. The velocity field for the complex melt flow is computed using the simplified model and scattering predictions are made by feeding these flow histories into the detailed model. The modeling quantitatively captures the full scattering intensity patterns over a broad range of data with independent variation of position within the contraction geometry, bulk flow rate and melt molecular weight. The study provides a strong, quantitative validation of current theoretical ideas concerning the microscopic dynamics of entangled polymers which builds upon existing comparisons with nonlinear mechanical stress data. Furthermore, we are able to confirm the appreciable length scale dependence of relaxation in polymer melts and highlight some wider implications of this phenomenon.

Numerical Study of Hypersonic Leeward Flow over a Blunt Nosed Delta Wing

A kinetic theory based Navier-Stokes solver has been implemented on a parallel supercomputer (Intel iPSC Touchstone Delta) to study the leeward flowfield of a blunt nosed delta wing at 30-deg incidence at hypersonic speeds (similar to the proposed HERMES aerospace plane). Computational results are presented for a series of grids for both inviscid and laminar viscous flows at Reynolds numbers of 225,000 and 2.25 million. In addition, comparisons are made between the present and two independent calculations of the some flows (by L. LeToullec and P. Guillen, and S. Menne) which were presented at the Workshop on Hypersonic Flows for Re-entry Problems, Antibes, France, 1991.