A novel numerical approximation for the space fractional advection-dispersion equation

In this paper, we consider a space fractional advection–dispersion equation. The equation is obtained from the standard advection–diffusion equation by replacing the first- and second-order space derivatives by the Riesz fractional derivatives of order β1 ∈ (0, 1) and β2 ∈ (1, 2], respectively. The fractional advection and dispersion terms are approximated by using two fractional centred difference schemes. A new weighted Riesz fractional finite-difference approximation scheme is proposed. When the weighting factor θ equals 12, a second-order accuracy scheme is obtained. The stability, consistency and convergence of the numerical approximation scheme are discussed. A numerical example is given to show that the numerical results are in good agreement with our theoretical analysis.

A comparison of finite difference and finite volume methods for solving the space fractional advection-dispersion equation with variable coefficients

Transport processes within heterogeneous media may exhibit non-classical diffusion or dispersion; that is, not adequately described by the classical theory of Brownian motion and Fick's law. We consider a space fractional advection-dispersion equation based on a fractional Fick's law. The equation involves the Riemann-Liouville fractional derivative which arises from assuming that particles may make large jumps. Finite difference methods for solving this equation have been proposed by Meerschaert and Tadjeran. In the variable coefficient case, the product rule is first applied, and then the Riemann-Liouville fractional derivatives are discretised using standard and shifted Grunwald formulas, depending on the fractional order. In this work, we consider a finite volume method that deals directly with the equation in conservative form. Fractionally-shifted Grunwald formulas are used to discretise the fractional derivatives at control volume faces. We compare the two methods for several case studies from the literature, highlighting the convenience of the finite volume approach.

A catchment modelling approach integrating surface and groundwater processes, land use and distribution of nutrients : Elimbah Creek, southeast Queensland

As the world’s population is growing, so is the demand for agricultural products. However, natural nitrogen (N) fixation and phosphorus (P) availability cannot sustain the rising agricultural production, thus, the application of N and P fertilisers as additional nutrient sources is common. It is those anthropogenic activities that can contribute high amounts of organic and inorganic nutrients to both surface and groundwaters resulting in degradation of water quality and a possible reduction of aquatic life. In addition, runoff and sewage from urban and residential areas can contain high amounts of inorganic and organic nutrients which may also affect water quality. For example, blooms of the cyanobacterium Lyngbya majuscula along the coastline of southeast Queensland are an indicator of at least short term decreases of water quality. Although Australian catchments, including those with intensive forms of land use, show in general a low export of nutrients compared to North American and European catchments, certain land use practices may still have a detrimental effect on the coastal environment. Numerous studies are reported on nutrient cycling and associated processes on a catchment scale in the Northern Hemisphere. Comparable studies in Australia, in particular in subtropical regions are, however, limited and there is a paucity in the data, in particular for inorganic and organic forms of nitrogen and phosphorus; these nutrients are important limiting factors in surface waters to promote algal blooms. Therefore, the monitoring of N and P and understanding the sources and pathways of these nutrients within a catchment is important in coastal zone management. Although Australia is the driest continent, in subtropical regions such as southeast Queensland, rainfall patterns have a significant effect on runoff and thus the nutrient cycle at a catchment scale. Increasingly, these rainfall patterns are becoming variable. The monitoring of these climatic conditions and the hydrological response of agricultural catchments is therefore also important to reduce the anthropogenic effects on surface and groundwater quality. This study consists of an integrated hydrological–hydrochemical approach that assesses N and P in an environment with multiple land uses. The main aim is to determine the nutrient cycle within a representative coastal catchment in southeast Queensland, the Elimbah Creek catchment. In particular, the investigation confirms the influence associated with forestry and agriculture on N and P forms, sources, distribution and fate in the surface and groundwaters of this subtropical setting. In addition, the study determines whether N and P are subject to transport into the adjacent estuary and thus into the marine environment; also considered is the effect of local topography, soils and geology on N and P sources and distribution. The thesis is structured on four components individually reported. The first paper determines the controls of catchment settings and processes on stream water, riverbank sediment, and shallow groundwater N and P concentrations, in particular during the extended dry conditions that were encountered during the study. Temporal and spatial factors such as seasonal changes, soil character, land use and catchment morphology are considered as well as their effect on controls over distributions of N and P in surface waters and associated groundwater. A total number of 30 surface and 13 shallow groundwater sampling sites were established throughout the catchment to represent dominant soil types and the land use upstream of each sampling location. Sampling comprises five rounds and was conducted over one year between October 2008 and November 2009. Surface water and groundwater samples were analysed for all major dissolved inorganic forms of N and for total N. Phosphorus was determined in the form of dissolved reactive P (predominantly orthophosphate) and total P. In addition, extracts of stream bank sediments and soil grab samples were analysed for these N and P species. Findings show that major storm events, in particular after long periods of drought conditions, are the driving force of N cycling. This is expressed by higher inorganic N concentrations in the agricultural subcatchment compared to the forested subcatchment. Nitrate N is the dominant inorganic form of N in both the surface and groundwaters and values are significantly higher in the groundwaters. Concentrations in the surface water range from 0.03 to 0.34 mg N L..1; organic N concentrations are considerably higher (average range: 0.33 to 0.85 mg N L..1), in particular in the forested subcatchment. Average NO3-N in the groundwater has a range of 0.39 to 2.08 mg N L..1, and organic N averages between 0.07 and 0.3 mg N L..1. The stream bank sediments are dominated by organic N (range: 0.53 to 0.65 mg N L..1), and the dominant inorganic form of N is NH4-N with values ranging between 0.38 and 0.41 mg N L..1. Topography and soils, however, were not to have a significant effect on N and P concentrations in waters. Detectable phosphorus in the surface and groundwaters of the catchment is limited to several locations typically in the proximity of areas with intensive animal use; in soil and sediments, P is negligible. In the second paper, the stable isotopes of N (14N/15N) and H2O (16O/18O and 2H/H) in surface and groundwaters are used to identify sources of dissolved inorganic and organic N in these waters, and to determine their pathways within the catchment; specific emphasis is placed on the relation of forestry and agriculture. Forestry is predominantly concentrated in the northern subcatchment (Beerburrum Creek) while agriculture is mainly found in the southern subcatchment (Six Mile Creek). Results show that agriculture (horticulture, crops, grazing) is the main source of inorganic N in the surface waters of the agricultural subcatchment, and their isotopic signature shows a close link to evaporation processes that may occur during water storage in farm dams that are used for irrigation. Groundwaters are subject to denitrification processes that may result in reduced dissolved inorganic N concentrations. Soil organic matter delivers most of the inorganic N to the surface water in the forested subcatchment. Here, precipitation and subsequently runoff is the main source of the surface waters. Groundwater in this area is affected by agricultural processes. The findings also show that the catchment can attenuate the effects of anthropogenic land use on surface water quality. Riparian strips of natural remnant vegetation, commonly 50 to 100 m in width, act as buffer zones along the drainage lines in the catchment and remove inorganic N from the soil water before it enters the creek. These riparian buffer zones are common in most agricultural catchments of southeast Queensland and are indicated to reduce the impact of agriculture on stream water quality and subsequently on the estuary and marine environments. This reduction is expressed by a significant decrease in DIN concentrations from 1.6 mg N L..1 to 0.09 mg N L..1, and a decrease in the �15N signatures from upstream surface water locations downstream to the outlet of the agricultural subcatchment. Further testing is, however, necessary to confirm these processes. Most importantly, the amount of N that is transported to the adjacent estuary is shown to be negligible. The third and fourth components of the thesis use a hydrological catchment model approach to determine the water balance of the Elimbah Creek catchment. The model is then used to simulate the effects of land use on the water balance and nutrient loads of the study area. The tool that is used is the internationally widely applied Soil and Water Assessment Tool (SWAT). Knowledge about the water cycle of a catchment is imperative in nutrient studies as processes such as rainfall, surface runoff, soil infiltration and routing of water through the drainage system are the driving forces of the catchment nutrient cycle. Long-term information about discharge volumes of the creeks and rivers do, however, not exist for a number of agricultural catchments in southeast Queensland, and such information is necessary to calibrate and validate numerical models. Therefore, a two-step modelling approach was used to calibrate and validate parameters values from a near-by gauged reference catchment as starting values for the ungauged Elimbah Creek catchment. Transposing monthly calibrated and validated parameter values from the reference catchment to the ungauged catchment significantly improved model performance showing that the hydrological model of the catchment of interest is a strong predictor of the water water balance. The model efficiency coefficient EF shows that 94% of the simulated discharge matches the observed flow whereas only 54% of the observed streamflow was simulated by the SWAT model prior to using the validated values from the reference catchment. In addition, the hydrological model confirmed that total surface runoff contributes the majority of flow to the surface water in the catchment (65%). Only a small proportion of the water in the creek is contributed by total base-flow (35%). This finding supports the results of the stable isotopes 16O/18O and 2H/H, which show the main source of water in the creeks is either from local precipitation or irrigation waters delivered by surface runoff; a contribution from the groundwater (baseflow) to the creeks could not be identified using 16O/18O and 2H/H. In addition, the SWAT model calculated that around 68% of the rainfall occurring in the catchment is lost through evapotranspiration reflecting the prevailing long-term drought conditions that were observed prior and during the study. Stream discharge from the forested subcatchment was an order of magnitude lower than discharge from the agricultural Six Mile Creek subcatchment. A change in land use from forestry to agriculture did not significantly change the catchment water balance, however, nutrient loads increased considerably. Conversely, a simulated change from agriculture to forestry resulted in a significant decrease of nitrogen loads. The findings of the thesis and the approach used are shown to be of value to catchment water quality monitoring on a wider scale, in particular the implications of mixed land use on nutrient forms, distributions and concentrations. The study confirms that in the tropics and subtropics the water balance is affected by extended dry periods and seasonal rainfall with intensive storm events. In particular, the comprehensive data set of inorganic and organic N and P forms in the surface and groundwaters of this subtropical setting acquired during the one year sampling program may be used in similar catchment hydrological studies where these detailed information is missing. Also, the study concludes that riparian buffer zones along the catchment drainage system attenuate the transport of nitrogen from agricultural sources in the surface water. Concentrations of N decreased from upstream to downstream locations and were negligible at the outlet of the catchment.

Engineering the surface of a sugarcane roller mill roll for grip and durability

Roller mills are typically used to crush sugarcane to express the juice from which sugar is manufactured. The mill rolls need to provide sufficient grip to ensure minimal sliding of the sugarcane along the roll surface. The rolls are subject to pressures up to 55 MPa from the sugarcane bagasse (as the sugarcane is called after first being crushed between a pair of rolls). The insoluble component of sugarcane includes typically 10% ash that largely originates from soil that is harvested with the cane. The sugarcane juice is acidic with pH typically between 5.0 and 5.5. As a result of ash and juice, the mill rolls are subjected to a range of abrasive and corrosive wear mechanisms. Solutions to provide grip and resist wear involve the selection of an appropriate roll shell material and compatible hard facing to provide the desired grip and wear characteristics. This paper reviews the various solutions that have been adopted for grip and durability for mill rolls and highlights the advantages and disadvantages of each method.

Tragedy and assimilation : reciprocity between surface and subject

Here we unveil a tragic triptych of three Australian women painfully painted onto the walls of interior surfaces. The woman at the centre of the triptych is Florence Broadhurst whose tragic death still remains a mystery. To the right is Australian skin illustrator Emma Hack who recreates Broadhurst’s wallpapers, mimicking their colourful patterns onto live models. Hack perfectly assimilates the models’ body into the wallpaper, camouflaging bodies except for small hints at something more in the foreground. In the process of Hack’s images, the models become statues, standing painfully still holding their breath for minutes at a time. The third woman, to the left of the triptych, is the fictional character Candy from the 2006 Australian film Candy. Candy’s traumatic struggle with addiction ends with her conveying her pain in a poem she writes on the walls of her home; culminating her tragic story into a disturbed domestic wall surface. This research tries to understand this relationship with the surface through tragedy as a reciprocal agreement between surface and subject and not a permanent transference between one state and another. What the surface provides in times of personal struggle and turmoil is a method for us to come to terms with out material existence.

Ultra-trace detection of diagnostically important biomarkers using functionalised-Surface Enhanced Raman Spectroscopy (SERS)

Here we report an ultrasensitive method for detecting bio-active compounds in biological samples by means of functionalised nanoparticles interrogated by surface enhanced Raman spectroscopy (SERS). This method is applicable to the recovery and detection of many diagnostically important peptidyl analytes such as insulin, human growth hormone, growth factors (IGFs) and erythropoietin (EPO), as well as many small molecule analytes and metabolites. Our method, developed to detect EPO, demonstrates its utility in a complex yet well defined biological system. Recombinant human EPO (rhEPO) and EPO analogues have successfully been used to treat anaemia in end-stage renal failure, chronic disorders and infections, cancer and AIDS. Current methods for EPO testing are lengthy, laborious and relatively insensitive to low concentrations. In our rapid screening methodology, gold nanoparticles were functionalised with anti-EPO antibodies to provide very high selectivity towards the EPO protein in urine. These “smart sensor” nanoparticles interact with and trap EPO. Subsequent SERS screening allows for the detection and quantisation of ultra trace amounts (<<10-15 M) of EPO in urine samples with minimal sample preparation. We present data showing that the SERS spectrum differentiates between human endogenous EPO and rhEPO in unpurified urine, and potentially distinguishes between purified EPO isoforms. The elimination of sample preparation and direct screening in biological fluids significantly reduces the time required by current methods. Antibody recognition against a variety of biological targets and the availability of portable commercial SERS analysers for rapid onsite testing suggest broad diagnostic applicability in a flexible analytical platform.

Natural articular joints : model of lamellar-rollerbearing lubrication and the nature of the cartilage surface

The influence of pH on interfacial energy and wettability distributed over the phospholipid bilayer surface were studied, and the importance of cartilage hydrophobicity (wettability) on the coefficient of friction (f) was established. It is argued that the wettability of cartilage signifi antly depends on the number of phospholipid bilayers acting as solid lubricant; the hypothesis was proven by conducting friction tests with normal and lipid- depleted cartilage samples. A lamellar-roller-bearing lubrication model was devised involving two mechanisms: (i) lamellar frictionless movement of bilayers, and (ii) roller-bearing lubrication mode through structured synovial fluid, which operates when lamellar spheres, liposomes and macromolecules act like a roller-bearing situated between two cartilage surfaces in effective biological lubrication.

The influence of temperature and introduction point on the detection of Rhyzopertha dominica in stored grain

The presence of insect pests in grain storages throughout the supply chain is a significant problem for farmers, grain handlers, and distributors world-wide. Insect monitoring and sampling programmes are used in the stored grains industry for the detection and estimation of pest populations. At the low pest densities dictated by economic and commercial requirements, the accuracy of both detection and abundance estimates can be influenced by variations in the spatial structure of pest populations over short distances. Geostatistical analysis of Rhyzopertha dominica populations in 2 and 3 dimensions showed that insect numbers were positively correlated over short (0.5 cm) distances, and negatively correlated over longer (.10 cm) distances. At 35 C, insects were located significantly further from the grain surface than at 25 and 30 C. Dispersion metrics showed statistically significant aggregation in all cases. The observed heterogeneous spatial distribution of R. dominica may also be influenced by factors such as the site of initial infestation and disturbance during handling. To account for these additional factors, I significantly extended a simulation model that incorporates both pest growth and movement through a typical stored-grain supply chain. By incorporating the effects of abundance, initial infestation site, grain handling, and treatment on pest spatial distribution, I developed a supply chain model incorporating estimates of pest spatial distribution. This was used to examine several scenarios representative of grain movement through a supply chain, and determine the influence of infestation location and grain disturbance on the sampling intensity required to detect pest infestations at various infestation rates. This study has investigated the effects of temperature, infestation point, and grain handling on the spatial distribution and detection of R. dominica. The proportion of grain infested was found to be dependent upon abundance, initial pest location, and grain handling. Simulation modelling indicated that accounting for these factors when developing sampling strategies for stored grain has the potential to significantly reduce sampling costs while simultaneously improving detection rate, resulting in reduced storage and pest management cost while improving grain quality.

Interaction of water with the fluorine-covered anatase TiO2(001) surface

The interaction of water with the fluorine-covered (001) surface of anatase titanium dioxide (TiO2) has been studied within the framework of density functional theory (DFT). Our results show that water dissociation is unfavorable due to repulsive interactions between surface fluorine and oxygen. We also found that the reaction of hydrofluoric acid with a surface hydroxyl group to form a surface Ti–F bond is exothermic, while the removal of fluorine from the surface needs additional energy of about half an eV. Therefore, water molecules are predicted to remain intact at the interface with the F-terminated anatase (001).

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.

Role of lithium vacancies in accelerating the dehydrogenation kinetics on a LiBH4(010) surface : an Ab initio study

In this work, ab initio density functional calculations were performed to explore the effect of surface lithium vacancies on the initial dehydrogenation kinetics of lithium borohydride. We found that some B−H bonds in neighboring BH4-1 complexes around the vacancy became elongated (weakened). The activation barriers for the recombination of H atoms to form H2 were decreased from 3.64 eV for the stoichiometrically complete LiBH4(010) surface to 1.53 and 0.23 eV in the presence of mono- and di-vacancies, respectively. Our results indicate that the creation of Li vacancies may play a critical role in accelerating the dehydrogenation kinetics of LiBH4.

Vacancy mediated desorption of hydrogen from a sodium alanate surface : An ab initio spin-polarized study

Ab initio spin-polarized density functional theory calculations are performed to explore the effect of single Na vacancy on NaAlH4(001) surface on the initial dehydrogenation kinetics. The authors found that two Al–H bond lengths become elongated and weakened due to the presence of a Na vacancy on the NaAlH4(001) surface. Spontaneous recombination from the surface to form molecular hydrogen is observed in the spin-polarized ab initio molecular dynamics simulation. The authors’ results indicate that surface Na vacancies play a critical role in accelerating the dehydrogenation kinetics in sodium alanate. The understanding gained here will aid in the rational design and development of complex hydride materials for hydrogen storage

First-principle study of adsorption of hydrogen on ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H2 molecules, respectively). Additionally, a molecular adsorption state of H2 above the Ti atom is observed for the first time and is attributed to the polarization of the H2 molecule by the Ti cation. Our results parallel recent findings for H2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.

Ab initio studies of hydrogen desorption from low index magnesium hydride surface

The low index Magnesium hydride surfaces, MgH2(0 0 1) and MgH2(1 1 0), have been studied by ab intio Density Functional Theory (DFT) calculations. It was found that the MgH2(1 1 0) surface is more stable than MgH2(0 0 1) surface, which is in good agreement with the experimental observation. The H2 desorption barriers vary depending on the crystalline surfaces that are exposed and also the specific H atom sites involved – they are found to be generally high, due to the thermodynamic stability of the MgH2 system, and are larger for the MgH2(0 0 1) surface. The pathway for recombinative desorption of one in-plane and one bridging H atom from the MgH2(1 1 0) surface was found to be the lowest energy barrier amongst those computed (172 KJ/mol) and is in good agreement with the experimental estimates.

Catalytic effects of subsurface carbon in the chemisorption of hydrogen on a Mg(0001) surface: an ab-initio study

Ab initio density functional theory (DFT) calculations are performed to explore possible catalytic effects on the dissociative chemisorption of hydrogen on a Mg(0001) surface when carbon is incorporated into Mg materials. The computational results imply that a C atom located initially on a Mg(0001) surface can migrate into the subsurface and occupy an fcc interstitial site, with charge transfer to the C atom from neighboring Mg atoms. The effect of subsurface C on the dissociation of H2 on the Mg(0001) surface is found to be relatively marginal: a perfect sublayer of interstitial C is calculated to lower the barrier by 0.16 eV compared with that on a pure Mg(0001) surface. Further calculations reveal, however, that sublayer C may have a significant effect in enhancing the diffusion of atomic hydrogen into the sublayers through fcc channels. This contributes new physical understanding toward rationalizing the experimentally observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.