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.

Visualisation of groundwater flow using texture based visualisation techniques

Texture based techniques for visualisation of unsteady vector fields have been applied for the visualisation of a Finite volume model for variably saturated groundwater flow through porous media. This model has been developed by staff in the School of Mathematical Sciences QUT for the study of salt water intrusion into coastal aquifers. This presentation discusses the implementation and effectiveness of the IBFV algorithm in the context of visualisation of the groundwater simulation outputs.

Final report : a software diagnostic tool for evaluating distress mechanisms in bridges exposed to aggressive environments

Durability issues of reinforced concrete construction cost millions of dollars in repair or demolition. Identification of the causes of degradation and a prediction of service life based on experience, judgement and local knowledge has limitations in addressing all the associated issues. The objective of this CRC CI research project is to develop a tool that will assist in the interpretation of the symptoms of degradation of concrete structures, estimate residual capacity and recommend cost effective solutions. This report is a documentation of the research undertaken in connection with this project. The primary focus of this research is centred on the case studies provided by Queensland Department of Main Roads (QDMR) and Brisbane City Council (BCC). These organisations are endowed with the responsibility of managing a huge volume of bridge infrastructure in the state of Queensland, Australia. The main issue to be addressed in managing these structures is the deterioration of bridge stock leading to a reduction in service life. Other issues such as political backlash, public inconvenience, approach land acquisitions are crucial but are not within the scope of this project. It is to be noted that deterioration is accentuated by aggressive environments such as salt water, acidic or sodic soils. Carse, 2005, has noted that the road authorities need to invest their first dollars in understanding their local concretes and optimising the durability performance of structures and then look at potential remedial strategies.

Effect of height of bridge above water on salt deposition levels

Maintenance of bridge structures is a major issue for the Queensland Department of Main Roads. In the previous phase of this CRC project an initial approach was made towards the development of a program for lifetime prediction of metallic bridge components. This involved the analysis of five representative bridge structures with respect to salt deposition (a major contributor to metallic corrosion) to determine common elements to be used as “cases” - those defined for buildings are not applicable. The five bridges analysed included the Gladstone Port Access Road Overpass, Stewart Road Overpass, South Johnstone River Bridge, Johnson Creek Bridge and the Ward River Bridge.

Does sea-level rise have an impact on saltwater intrusion?

Climate change effects are expected to substantially raise the average sea level. It is widely assumed that this raise will have a severe adverse impact on saltwater intrusion processes in coastal aquifers. In this study we hypothesize that a natural mechanism, identified as the “lifting process” has the potential to mitigate or in some cases completely reverse the adverse intrusion effects induced by sea-level rise. A detailed numerical study using the MODFLOW-family computer code SEAWAT, was completed to test this hypothesis and to understand the effects of this lifting process in both confined and unconfined systems. Our conceptual simulation results show that if the ambient recharge remains constant, the sea-level rise will have no long-term impact (i.e., it will not affect the steady-state salt wedge) on confined aquifers. Our transient confined flow simulations show a self-reversal mechanism where the wedge which will initially intrude into the formation due to the sea-level rise would be naturally driven back to the original position. In unconfined systems, the lifting process would have a lesser influence due to changes in the value of effective transmissivity. A detailed sensitivity analysis was also completed to understand the sensitivity of this self-reversal effect to various aquifer parameters.

The structure of the copper(II) complex with 4,5-dichlorophthalic acid: The 1:1 complex 'adduct' salt trans-tetraaqua(2-carboxy-4,5-dichlorobenzoato)(4,5-dichlorobenzene-1,2-dicarboxylic acid)copper(II) 2-carboxy-4,5-dichlorobenzoate

The unusual (1:1) complex ‘adduct’ salt of copper(II) with 4,5-dichlorophthalic acid (H2DCPA), having formula [Cu(H2O)4(C8H3Cl2O4) (C8H4Cl2O4)] . (C8H3Cl2O4) has been synthesized and characterized using single-crystal X-ray diffraction. Crystals are monoclinic, space group P21/c, with Z = 4 in a cell with dimensions a = 20.1376(7), b =12.8408(4) c = 12.1910(4) Å, β = 105.509(4)o. The complex is based on discrete tetragonally distorted octahedral [CuO6] coordination centres with the four water ligands occupying the square planar sites [Cu-O, 1.962(4)-1.987(4) Å] and the monodentate carboxyl-O donors of two DCPA ligand species in the axial sites. The first of these bonds [Cu-O, 2.341(4) Å] is with an oxygen of a HDCPA monoanion, the second with an oxygen of a H2DCPA acid species [Cu-O, 2.418(4) Å]. The un-coordinated ‘adduct’ molecule is a HDCPA counter anion which is strongly hydrogen-bonded to the coordinated H2DCPA ligand [O… O, 2.503(6) Å] while a number of peripheral intra- and intermolecular hydrogen-bonding interactions give a two-dimensional network structure.

2,3-Dimethoxy-10-oxostrychnidinium 2-(2,4,6-trinitroanilino)benzoate monohydrate : a 1:1 proton-transfer salt of brucine with o-picraminobenzoic acid

In the structure of the 1:1 proton-transfer compound of brucine with 2-(2,4,6-trinitroanilino)benzoic acid C23H27N2O4+ . C13H7N4O8- . H~2~O, the brucinium cations form the classic undulating ribbon substructures through overlapping head-to-tail interactions while the anions and the three related partial water molecules of solvation (having occupancies of 0.73, 0.17 and 0.10) occupy the interstitial regions of the structure. The cations are linked to the anions directly through N-H...O(carboxyl) hydrogen bonds and indirectly by the three water molecules which form similar conjoint cyclic bridging units [graph set R2/4(8)] through O-H...O(carbonyl) and O(carboxyl) hydrogen bonds, giving a two-dimensional layered structure. Within the anion, intramolecular N-H...O(carboxyl) and N H...O(nitro) hydrogen bonds result in the benzoate and picrate rings being rotated slightly out of coplanarity inter-ring dihedral angle 32.50(14)\%]. This work provides another example of the molecular selectivity of brucine in forming stable crystal structures and also represents the first reported structure of any form of the guest compound 2-(2,4,6-trinitroanilino)benzoic acid.

Resolution of the chiral (1R,2S) enantiomer of cis-cyclohexane-1,2-dicarboxylic acid in the brucinium salt 2,3-dimethoxy-10-oxostrychnidinium (1R,2S)-2-carboxycyclohexane-1-carboxylate dihydrate

The structure of the 1:1 brucinium salt of cis-cyclohexane-1,2-dicarboxylic acid, 2,3-dimethoxy-10-oxostrychnidinium (1R,2S)-2-carboxycyclohexane-1-carboxylate dihydrate, has revealed the resolved (1R,2S) enantiomer of the acid. Crystals of the compound are orthorhombic, space group P212121, with unit cell dimensions a = 8.1955(3), b = 12.4034(3), c = 29.9073(9)Å, and Z = 4. The asymmetric unit comprises the brucinium cation, the hydrogen cis-cyclohexane-1,2-dicarboxylate cation, in which the carboxylate group is disordered over two sites (58, 42%), and two water molecules of solvation, one of which is occupies two 50% occupancy sites. The classic undulating brucinium cation substructures are present with the anion and the water molecules occupying the interstitial cavities and are hydrogen-bonded to them in a two-dimensional network structure.

Water absorption, permeability, and resistance to chloride-ion penetration of lightweight aggregate concrete

This paper presents an experimental study to evaluate effect of cumulative lightweight aggregate (LWA) content (including lightweight sand) in concrete [water/cement ratio (w/c) = 0.38] on its water absorption, water permeability, and resistance to chloride-ion penetration. Rapid chloride penetrability test (ASTM C 1202), rapid migration test (NT Build 492), and salt ponding test (AASHTO T 259) were conducted to evaluate the concrete resistance to chloride-ion penetration. The results were compared with those of a cement paste and a control normal weight aggregate concrete (NWAC) with the same w/c and a NWAC (w/c = 0.54) with 28-day compressive strength similar to some of the lightweight aggregate concrete (LWAC). Results indicate that although the total charge passed, migration coefficient, and diffusion coefficient of the LWAC were not significantly different from those of NWAC with the same w/c of 0.38, resistance of the LWAC to chloride penetration decreased with increase in the cumulative LWA content in the concretes. The water penetration depth under pressure and water sorptivity showed, in general, similar trends. The LWAC with only coarse LWA had similar water sorptivity, water permeability coefficient, and resistance to chloride-ion penetration compared to NWAC with similar w/c. The LWAC had lower water sorptivity, water permeability and higher resistance to chloride-ion penetration than the NWAC with similar 28-day strength but higher w/c. Both the NWAC and LWAC had lower sorptivity and higher resistance to chloride-ion penetration than the cement paste with similar w/c.

Physiological basis of salt stress tolerance in rice expressing the antiapoptotic gene SfIAP

Programmed cell death-associated genes, especially antiapoptosis-related genes have been reported to confer tolerance to a wide range of biotic and abiotic stresses in dicotyledonous plants such as tobacco (Nicotiana tabacum L.) and tomato (Solanum lycopersicum L.). This is the first time the antiapoptotic gene SfIAP was transformed into a monocotyledonous representative: rice (Oryza sativa L.). Transgenic rice strains expressing SfIAP were generated by the Agrobacterium-mediated transformation method and rice embryogenic calli, and assessed for their ability to confer tolerance to salt stress at both the seedling and reproductive stages using a combination of molecular, agronomical, physiological and biochemical techniques. The results show that plants expressing SfIAP have higher salt tolerance levels in comparison to the wild-type and vector controls. By preventing cell death at the onset of salt stress and maintaining the cell membrane’s integrity, SfIAP transgenic rice plants can retain plant water status, ion homeostasis, photosynthetic efficiency and growth to combat salinity successfully.