Modelling sea water intrusion in coastal aquifers using heterogeneous computing

The objective of this PhD research program is to investigate numerical methods for simulating variably-saturated flow and sea water intrusion in coastal aquifers in a high-performance computing environment. The work is divided into three overlapping tasks: to develop an accurate and stable finite volume discretisation and numerical solution strategy for the variably-saturated flow and salt transport equations; to implement the chosen approach in a high performance computing environment that may have multiple GPUs or CPU cores; and to verify and test the implementation. The geological description of aquifers is often complex, with porous materials possessing highly variable properties, that are best described using unstructured meshes. The finite volume method is a popular method for the solution of the conservation laws that describe sea water intrusion, and is well-suited to unstructured meshes. In this work we apply a control volume-finite element (CV-FE) method to an extension of a recently proposed formulation (Kees and Miller, 2002) for variably saturated groundwater flow. The CV-FE method evaluates fluxes at points where material properties and gradients in pressure and concentration are consistently defined, making it both suitable for heterogeneous media and mass conservative. Using the method of lines, the CV-FE discretisation gives a set of differential algebraic equations (DAEs) amenable to solution using higher-order implicit solvers. Heterogeneous computer systems that use a combination of computational hardware such as CPUs and GPUs, are attractive for scientific computing due to the potential advantages offered by GPUs for accelerating data-parallel operations. We present a C++ library that implements data-parallel methods on both CPU and GPUs. The finite volume discretisation is expressed in terms of these data-parallel operations, which gives an efficient implementation of the nonlinear residual function. This makes the implicit solution of the DAE system possible on the GPU, because the inexact Newton-Krylov method used by the implicit time stepping scheme can approximate the action of a matrix on a vector using residual evaluations. We also propose preconditioning strategies that are amenable to GPU implementation, so that all computationally-intensive aspects of the implicit time stepping scheme are implemented on the GPU. Results are presented that demonstrate the efficiency and accuracy of the proposed numeric methods and formulation. The formulation offers excellent conservation of mass, and higher-order temporal integration increases both numeric efficiency and accuracy of the solutions. Flux limiting produces accurate, oscillation-free solutions on coarse meshes, where much finer meshes are required to obtain solutions with equivalent accuracy using upstream weighting. The computational efficiency of the software is investigated using CPUs and GPUs on a high-performance workstation. The GPU version offers considerable speedup over the CPU version, with one GPU giving speedup factor of 3 over the eight-core CPU implementation.

Role of land use and seasonal factors in water quality degradation

Surface water and groundwater are the most important water sources in the natural environment. Land use and seasonal factors play an important role in influencing the quality of these water sources. An in-depth understanding of the role of these two influential factors can help to implement an effective catchment management strategy for the protection of these water sources. This paper discusses the outcomes of an extensive research study which investigated the role of land use and seasonal factors on surface water and groundwater pollution in a mixed land use coastal catchment. The study confirmed that the influence exerted on the water environment by seasonal factors is secondary to that of land use. Furthermore, the influence of land use and seasonal factors on surface water and groundwater quality varies with the pollutant species. This highlights the need to specifically take into consideration the targeted pollutants and the key influential factors for the effective protection of vulnerable receiving water environments.

The use of an electrostatic lens to enhance the efficiency of the electrospinning process

Electrospun scaffolds manufactured using conventional electrospinning configurations have an intrinsic thickness limitation, due to a charge build-up at the collector. To overcome this limitation, an electrostatic lens has been developed that, at the same relative rate of deposition, focuses the polymer jet onto a smaller area of the collector, resulting in the fabrication of thick scaffolds within a shorter period of time. We also observed that a longer deposition time (up to 13 h, without the intervention of the operator) could be achieved when the electrostatic lens was utilised, compared to 9–10 h with a conventional processing set-up and also showed that fibre fusion was less likely to occur in the modified method. This had a significant impact on the mechanical properties, as the scaffolds obtained with the conventional process had a higher elastic modulus and ultimate stress and strain at short times. However, as the thickness of the scaffolds produced by the conventional electrospinning process increased, a 3-fold decrease in the mechanical properties was observed. This was in contrast to the modified method, which showed a continual increase in mechanical properties, with the properties of the scaffold finally having similar mechanical properties to the scaffolds obtained via the conventional process at longer times. This “focusing” device thus enabled the fabrication of thicker 3-dimensional electrospun scaffolds (of thicknesses up to 3.5 mm), representing an important step towards the production of scaffolds for tissue engineering large defect sites in a multitude of tissues.

Corporate sustainability reporting for water : water footprint, global reporting initiative and the water accounting framework

Three initiatives with respect to water reporting in the mining sector are compared in this paper to understand the quantities that are asked for by each initiative and the guidelines of those initiatives through means of a case study. The Global Reporting Initiative (GRI) was chosen because it has achieved widespread acceptance amongst mining companies and its water-related indicators are widely reported in corporate sustainability reporting. In contrast, the Water Footprint Network, which has been an important initiative in food and agricultural industries, has had low acceptance in the mining industry. The third initiative is the Water Accounting Framework, a collaboration between The Minerals Council of Australia and the Sustainable Minerals Institute of the University of Queensland. A water account had previously been created according to the Water Accounting Framework for the case study site, an open pit coal mine in the Bowen Basin. The resulting account provided consistent data for the Global Reporting Initiative (GRI) and the Water Footprint attributable to mining but in particular, a deficiency in the GRI indicator of EN10 reuse and recycling efficiency was illustrated quantitatively. This has far-reaching significance due to the widespread use of GRI indicators in mining corporate reports.

Water heating and dehumidification with heat pump utilizing solar, ambient and waste heat

The low temperature operation of a heat pump makes it an excellent match for the use of solar energy. At the National University of Singapore, a solar assisted heat pump system has been designed, fabricated and installed to provide water heating and drying. The system also utilizes the air con waste heat, which would normally be released to atmosphere adding to global warming. Experimental results show that the twophase unglazed solar evaporator-collector, instead of losing energy to the ambient, gained a significant amount due to low operating temperature of the collector. As a result, the collector efficiency attains a value greater than 1, when conventional collector equations are used. With this evaporator-collector, the system can be operated even in the absence of solar irradiation. The waste heat was collected from an air-con system, which maintained a room at 20-22 oC. In the condenser side, water at 60 oC was produced at a rate of 3 liter/minute and the drying capacity was 2.2kg/hour. Maximum COP of the system was found to be about 5.5.

Water-based nanoparticulate solar cells using a diketopyrrolopyrrole donor polymer

Organic photovoltaic devices with either bulk heterojunction (BHJ) or nanoparticulate (NP) active layers have been prepared from a 1:2 blend of (poly{3,6-dithiophene-2-yl-2,5-di(2-octyldodecyl)-pyrrolo[3,4-c]pyrrole-1, 4-dione-alt-naphthalene}) (PDPP-TNT) and the fullerene acceptor, ([6,6]-phenyl C71-butyric acid methyl ester) (PC70BM). Atomic force microscopy (AFM) and scanning electron microscopy (SEM) have been used to investigate the morphology of the active layers of the two approaches. Mild thermal treatment of the NP film is required to promote initial joining of the NPs in order for the devices to function, however the NP structure is retained. Consequently, whereas gross phase segregation of the active layer occurs in the BHJ device spin cast from chloroform, the nanoparticulate approach retains control of the material domain sizes on the length scale of exciton diffusion in the materials. As a result, NP devices are found to generate more than twice the current density of BHJ devices and have a substantially greater overall efficiency. The use of aqueous nanoparticulate dispersions offers a promising approach to control the donor acceptor morphology on the nanoscale with the benefit of environmentally- friendly, solution-based fabrication. © 2014 the Owner Societies.

Modeling chickpea growth and development: Nitrogen accumulation and use

Quantitative information regarding nitrogen (N) accumulation and its distribution to leaves, stems and grains under varying environmental and growth conditions are limited for chickpea (Cicer arietinum L.). The information is required for the development of crop growth models and also for assessment of the contribution of chickpea to N balances in cropping systems. Accordingly, these processes were quantified in chickpea under different environmental and growth conditions (still without water or N deficit) using four field experiments and 1325 N measurements. N concentration ([N]) in green leaves was 50 mg g-1 up to beginning of seed growth, and then it declined linearly to 30 mg g-1 at the end of seed growth phase. [N] in senesced leaves was 12 mg g-1. Stem [N] decreased from 30 mg g-1 early in the season to 8 mg g-1 in senesced stems at maturity. Pod [N] was constant (35 mg g-1), but grain [N] decreased from 60 mg g-1 early in seed growth to 43 mg g-1 at maturity. Total N accumulation ranged between 9 and 30 g m-2. N accumulation was closely linked to biomass accumulation until maturity. N accumulation efficiency (N accumulation relative to biomass accumulation) was 0.033 g g-1 where total biomass was -2 and during early growth period, but it decreased to 0.0176 g g-1 during the later growth period when total biomass was >218 g m-2. During vegetative growth (up to first-pod), 58% of N was partitioned to leaves and 42% to stems. Depending on growth conditions, 37-72% of leaf N and 12-56% of stem N was remobilized to the grains. The parameter estimates and functions obtained in this study can be used in chickpea simulation models to simulate N accumulation and distribution.

Remote sensing of nitrogen and water stress in wheat

Nitrogen (N) is the largest agricultural input in many Australian cropping systems and applying the right amount of N in the right place at the right physiological stage is a significant challenge for wheat growers. Optimizing N uptake could reduce input costs and minimize potential off-site movement. Since N uptake is dependent on soil and plant water status, ideally, N should be applied only to areas within paddocks with sufficient plant available water. To quantify N and water stress, spectral and thermal crop stress detection methods were explored using hyperspectral, multispectral and thermal remote sensing data collected at a research field site in Victoria, Australia. Wheat was grown over two seasons with two levels of water inputs (rainfall/irrigation) and either four levels (in 2004; 0, 17, 39 and 163 kg/ha) or two levels (in 2005; 0 and 39 kg/ha N) of nitrogen. The Canopy Chlorophyll Content Index (CCCI) and modified Spectral Ratio planar index (mSRpi), two indices designed to measure canopy-level N, were calculated from canopy-level hyperspectral data in 2005. They accounted for 76% and 74% of the variability of crop N status, respectively, just prior to stem elongation (Zadoks 24). The Normalised Difference Red Edge (NDRE) index and CCCI, calculated from airborne multispectral imagery, accounted for 41% and 37% of variability in crop N status, respectively. Greater scatter in the airborne data was attributable to the difference in scale of the ground and aerial measurements (i.e., small area plant samples against whole-plot means from imagery). Nevertheless, the analysis demonstrated that canopy-level theory can be transferred to airborne data, which could ultimately be of more use to growers. Thermal imagery showed that mean plot temperatures of rainfed treatments were 2.7 °C warmer than irrigated treatments (P < 0.001) at full cover. For partially vegetated fields, the two-Dimensional Crop Water Stress Index (2D CWSI) was calculated using the Vegetation Index-Temperature (VIT) trapezoid method to reduce the contribution of soil background to image temperature. Results showed rainfed plots were consistently more stressed than irrigated plots. Future work is needed to improve the ability of the CCCI and VIT methods to detect N and water stress and apply both indices simultaneously at the paddock scale to test whether N can be targeted based on water status. Use of these technologies has significant potential for maximising the spatial and temporal efficiency of N applications for wheat growers. ‘Ground–breaking Stuff’- Proceedings of the 13th Australian Society of Agronomy Conference, 10-14 September 2006, Perth, Western Australia.

Cooking in the ungra area: Fuel efficiency, energy losses, and opportunities for reducing firewood consumption

Cooking efficiency and related fuel economy issues have been studied in a particular rural area of India. Following a description of the cooking practices and conditions in this locale, cooking efficiency is examined. A cooking efficiency of only 6% was found. The use of aluminium rather than clay pots results in an increased efficiency. In addition, cooking efficiency correlates very well with specific fuel consumption. The latter parameter is much simpler to analyse than cooking efficiency. The energy losses during cooking are examined in the second part of this case study. The major energy losses are heating of excess air, heat carried away by the combustion products, heat transmitted to the stove body and floor, and the chemical energy in charcoal residue. The energy loss due to the evaporation of cooking water is also significant because it represents about one-third of the heat reaching the pots.

Using rowers’ perceptions of on-water stroke success to evaluate sculling catch efficiency variables via a boat instrumentation system

Aim An effective catch in sculling is a critical determinant of boat velocity. This study used rowers’ performance-based judgments to compare three measures of catch slip efficiency. Two questions were addressed: (1) would rower-judged Yes strokes be faster than No strokes? and (2) which method of quantifying catch slip best reflected these judgements? Methods Eight single scullers performed two 10-min blocks of sub maximal on-water rowing at 20 strokes per minute. Every 30 s, rowers reported either Yes or No about the quality of their stroke at the catch. Results It was found that Yes strokes identified by rowers had, on average, a moderate effect advantage over No strokes with a standardised effect size of 0.43. In addition, a quicker time to positive acceleration best reflected the change in performance; where the standardised mean difference score of 0.57 for time to positive acceleration was larger than the scores of 0.47 for time to PowerLine force, and 0.35 for time to 30% peak pin force catch slip measures. For all eight rowers, Yes strokes corresponded to time to positive acceleration occurring earlier than No strokes. Conclusion Rower judgements about successful strokes was linked to achieving a quicker time to positive acceleration, and may be of the most value in achieving a higher average boat velocity.

Water resources: Embed stormwater use in city planning [Correspondence]

Potable water resources are being depleted at an alarming rate worldwide. Storm water is a hugely under-utilized resource that could help as extreme weather events become more frequent...

Use of Abel integral equations in water wave scattering by two surface-piercing barriers

In this paper the classical problem of water wave scattering by two partially immersed plane vertical barriers submerged in deep water up to the same depth is investigated. This problem has an exact but complicated solution and an approximate solution in the literature of linearised theory of water waves. Using the Havelock expansion for the water wave potential, the problem is reduced here to solving Abel integral equations having exact solutions. Utilising these solutions,two sets of expressions for the reflection and transmission coefficients are obtained in closed forms in terms of computable integrals in contrast to the results given in the literature which,involved six complicated integrals in terms of elliptic functions. The two different expressions for each coefficient produce almost the same numerical results although it has not been possible to prove their equivalence analytically. The reflection coefficient is depicted against the wave number in a number of figures which almost coincide with the figures available in the literature wherein the problem was solved approximately by employing complementary approximations. (C) 2009 Elsevier B.V. All rights reserved.

Systems Analysis of Tank Irrigation: II. Delayed Start and Water Deficit

Tank irrigation systems in the semiarid regions of India are discussed in this paper. To optimize the grain yield of rice, it is essential to start the agricultural operations in the second week of July so that favorable climatic conditions will prevail during flowering and yield formation stages. Because of low inflow during the initial few weeks of the crop season, often farmers are forced to delay planting until sufficient sowing rain and inflow have occurred or to adopt deficit irrigation during this period. The delayed start affects the grain yield, but will lead to an improved irrigation efficiency. A delayed start of agricultural operations with increased irrigation efficiency leads to the energy resources becoming critical during the peak requirement week, particularly those of female labor and animal power. This necessitates augmenting these resources during weeks of their peak use, either by reorganizing the traditional methods of cultivation or by importing from outside the system.