A dual-scale modelling approach for drying hygroscopic porous media

A new dualscale modelling approach is presented for simulating the drying of a wet hygroscopic porous material that couples the porous medium (macroscale) with the underlying pore structure (microscale). The proposed model is applied to the convective drying of wood at low temperatures and is valid in the so-called hygroscopic range, where hygroscopically held liquid water is present in the solid phase and water exits only as vapour in the pores. Coupling between scales is achieved by imposing the macroscopic gradients of moisture content and temperature on the microscopic field using suitably-defined periodic boundary conditions, which allows the macroscopic mass and thermal fluxes to be defined as averages of the microscopic fluxes over the unit cell. This novel formulation accounts for the intricate coupling of heat and mass transfer at the microscopic scale but reduces to a classical homogenisation approach if a linear relationship is assumed between the microscopic gradient and flux. Simulation results for a sample of spruce wood highlight the potential and flexibility of the new dual-scale approach. In particular, for a given unit cell configuration it is not necessary to propose the form of the macroscopic fluxes prior to the simulations because these are determined as a direct result of the dual-scale formulation.

Suspended sediment properties and suspended sediment flux estimates in an inundated urban environment during a major flood event

During a major flood event, the inundation of urban environments leads to some complicated flow motion most often associated with significant sediment fluxes. In the present study, a series of field measurements were conducted in an inundated section of the City of Brisbane (Australia) about the peak of a major flood in January 2011. Some experiments were performed to use ADV backscatter amplitude as a surrogate estimate of the suspended sediment concentration (SSC) during the flood event. The flood water deposit samples were predominantly silty material with a median particle size about 25 μm and they exhibited a non-Newtonian behavior under rheological testing. In the inundated urban environment during the flood, estimates of suspended sediment concentration presented a general trend with increasing SSC for decreasing water depth. The suspended sediment flux data showed some substantial sediment flux amplitudes consistent with the murky appearance of floodwaters. Altogether the results highlighted the large suspended sediment loads and fluctuations in the inundated urban setting associated possibly with a non-Newtonian behavior. During the receding flood, some unusual long-period oscillations were observed (periods about 18 min), although the cause of these oscillations remains unknown. The field deployment was conducted in challenging conditions highlighting a number of practical issues during a natural disaster.

Nitrous oxide emissions from subtropical horticultural soils : a time series analysis

Time series regression models were used to examine the influence of environmental factors (soil water content and soil temperature) on the emissions of nitrous oxide (N2O) from subtropical soils, by taking into account temporal lagged environmental factors, autoregressive processes, and seasonality for three horticultural crops in a subtropical region of Australia. Fluxes of N2O, soil water content, and soil temperature were determined simultaneously on a weekly basis over a 12-month period in South East Queensland. Annual N2O emissions for soils under mango, pineapple, and custard apple were 1590, 1156, and 2038 g N2O-N/ha, respectively, with most emissions attributed to nitrification. The N2O-N emitted from the pineapple and custard apple crops was equivalent to 0.26 and 2.22%, respectively, of the applied mineral N. The change in soil water content was the key variable for describing N2O emissions at the weekly time-scale, with soil temperature at a lag of 1 month having a significant influence on average N2O emissions (averaged) at the monthly time-scale across the three crops. After accounting for soil temperature and soil water content, both the weekly and monthly time series regression models exhibited significant autocorrelation at lags of 1–2 weeks and 1–2 months, and significant seasonality for weekly N2O emissions for mango crop and for monthly N2O emissions for mango and custard apple crops in this location over this time-frame. Time series regression models can explain a higher percentage of the temporal variation of N2O emission compared with simple regression models using soil temperature and soil water content as drivers. Taking into account seasonal variability and temporal persistence in N2O emissions associated with soil water content and soil temperature may lead to a reduction in the uncertainty surrounding estimates of N2O emissions based on limited sampling effort.

Soil N2O and CO2 emissions from cotton in Australia under varying irrigation management

Irrigation is known to stimulate soil microbial carbon and nitrogen turnover and potentially the emissions of nitrous oxide (N2O) and carbon dioxide (CO2). We conducted a study to evaluate the effect of three different irrigation intensities on soil N2O and CO2 fluxes and to determine if irrigation management can be used to mitigate N2O emissions from irrigated cotton on black vertisols in South-Eastern Queensland, Australia. Fluxes were measured over the entire 2009/2010 cotton growing season with a fully automated chamber system that measured emissions on a sub-daily basis. Irrigation intensity had a significant effect on CO2 emission. More frequent irrigation stimulated soil respiration and seasonal CO2 fluxes ranged from 2.7 to 4.1 Mg-C ha−1 for the treatments with the lowest and highest irrigation frequency, respectively. N2O emission happened episodic with highest emissions when heavy rainfall or irrigation coincided with elevated soil mineral N levels and seasonal emissions ranged from 0.80 to 1.07 kg N2O-N ha−1 for the different treatments. Emission factors (EF = proportion of N fertilizer emitted as N2O) over the cotton cropping season, uncorrected for background emissions, ranged from 0.40 to 0.53 % of total N applied for the different treatments. There was no significant effect of the different irrigation treatments on soil N2O fluxes because highest emission happened in all treatments following heavy rainfall caused by a series of summer thunderstorms which overrode the effect of the irrigation treatment. However, higher irrigation intensity increased the cotton yield and therefore reduced the N2O intensity (N2O emission per lint yield) of this cropping system. Our data suggest that there is only limited scope to reduce absolute N2O emissions by different irrigation intensities in irrigated cotton systems with summer dominated rainfall. However, the significant impact of the irrigation treatments on the N2O intensity clearly shows that irrigation can easily be used to optimize the N2O intensity of such a system.

Nitrous oxide emissions from irrigated wheat in Australia : impact of irrigation management

Background and Aims: Irrigation management affects soil water dynamics as well as the soil microbial carbon and nitrogen turnover and potentially the biosphere-atmosphere exchange of greenhouse gasses (GHG). We present a study on the effect of three irrigation treatments on the emissions of nitrous oxide (N2O) from irrigated wheat on black vertisols in South-Eastern Queensland, Australia. Methods: Soil N2O fluxes from wheat were monitored over one season with a fully automated system that measured emissions on a sub-daily basis. Measurements were taken from 3 subplots for each treatment within a randomized split-plot design. Results: Highest N2O emissions occurred after rainfall or irrigation and the amount of irrigation water applied was found to influence the magnitude of these “emission pulses”. Daily N2O emissions varied from -0.74 to 20.46 g N2O-N ha-1 day-1 resulting in seasonal losses ranging from 0.43 to 0.75 kg N2O N ha-1 season -1 for the different irrigation treatments. Emission factors (EF = proportion of N fertilizer emitted as N2O) over the wheat cropping season, uncorrected for background emissions, ranged from 0.2 to 0.4% of total N applied for the different treatments. Highest seasonal N2O emissions were observed in the treatment with the highest irrigation intensity; however, the N2O intensity (N2O emission per crop yield) was highest in the treatment with the lowest irrigation intensity. Conclusions: Our data suggest that timing and amount of irrigation can effectively be used to reduce N2O losses from irrigated agricultural systems; however, in order to develop sustainable mitigation strategies the N2O intensity of a cropping system is an important concept that needs to be taken into account.

Enhanced hydrogen separation by vertically-aligned carbon nanotube membranes with zeolite imidazolate frameworks as a selective layer

Vertically-aligned carbon nanotube (VACNT) membranes show very high permeation fluxes due to the inherent smooth and frictionless nature of the interior of the nanotubes. However, the hydrogen selectivities are all in the Knudsen range and are quite low. In this study we grew molecular sieve zeolite imidazolate frameworks (ZIFs) via secondary seeded growth on the VACNT membranes as a gas selective layer. The ZIF layer has a thickness of 5–6 μm and shows good contact with the VACNT membrane surface. The VACNT supported ZIF membrane shows much higher H2 selectivity than Ar (7.0); O2 (13.6); N2 (15.1) and CH4 (9.8). We conclude that tailoring metal–organic frameworks on the membrane surface can be an effective route to improve the gas separation performance of the VACNT membrane.

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.

Late Oligocene to Middle Miocene rifting and synextensional magmatism in the southwestern Sierra Madre Occidental, Mexico : the beginning of the Gulf of California rift

Although Basin and Range–style extension affected large areas of western Mexico after the Late Eocene, most consider that extension in the Gulf of California region began as subduction waned and ended ca. 14–12.5 Ma. A general consensus also exists in considering Early and Middle Miocene volcanism of the Sierra Madre Occidental and Comondú Group as subduction related, whereas volcanism after ca. 12.5 Ma is extension related. Here we present a new regional geologic study of the eastern Gulf of California margin in the states of Nayarit and Sinaloa, Mexico, backed by 43 new Ar-Ar and U-Pb mineral ages, and geochemical data that document an earlier widespread phase of extension. This extension across the southern and central Gulf Extensional Province began between Late Oligocene and Early Miocene time, but was focused in the region of the future Gulf of California in the Middle Miocene. Late Oligocene to Early Miocene rocks across northern Nayarit and southern Sinaloa were affected by major approximately north-south– to north-northwest– striking normal faults prior to ca. 21 Ma. Between ca. 21 and 11 Ma, a system of north-northwest–south-southeast high angle extensional faults continued extending the southwestern side of the Sierra Madre Occidental. Rhyolitic domes, shallow intrusive bodies, and lesser basalts were emplaced along this extensional belt at 20–17 Ma. Rhyolitic rocks, in particular the domes and lavas, often show strong antecrystic inheritance but only a few Mesozoic or older xenocrysts, suggesting silicic magma generation in the mid-upper crust triggered by an extension induced basaltic infl ux. In northern Sinaloa, large grabens were occupied by huge volcanic dome complexes ca. 21–17 Ma and filled by continental sediments with interlayered basalts dated as 15–14 Ma, a stratigraphy and timing very similar to those found in central Sonora (northeastern Gulf of California margin). Early to Middle Miocene volcanism occurred thus in rift basins, and was likely associated with decompression melting of upper mantle (inducing crustal partial melting) rather than with fluxing by fluids from the young and slow subducting microplates. Along the eastern side of the Gulf of California coast, from Farallón de San Ignacio island offshore Los Mochis, Sinaloa, to San Blas, Nayarit, a strike distance of >700 km, flat lying basaltic lavas dated as ca. 11.5–10 Ma are exposed just above the present sea level. Here crustal thickness is almost half that in the unextended core of the adjacent Sierra Madre Occidental, implying signifi cant lithosphere stretching before ca. 11 Ma. This mafic pulse, with subdued Nb-Ta negative spikes, may be related to the detachment of the lower part of the subducted slab, allowing an upward asthenospheric flow into an upper mantle previously modified by fluid fluxes related to past subduction. Widespread eruption of very uniform oceanic island basalt–like lavas occurred by the late Pliocene and Pleistocene, only 20 m.y. after the onset of rifting and ~9 m.y. after the end of subduction, implying that preexisting subduction-modified mantle had now become isolated from melt source regions. Our study shows that rifting across the southern-central Gulf Extensional Province began much earlier than the Late Miocene and provided a fundamental control on the style and composition of volcanism from at least 30 Ma. We envision a sustained period of lithospheric stretching and magmatism during which the pace and breadth of extension changed ca. 20–18 Ma to be narrower, and again after ca. 12.5 Ma, when the kinematics of rifting became more oblique.

Influence of nitrogen fertiliser application and timing on greenhouse gas emissions from a lychee (Litchi chinensis) orchard in humid subtropical Australia

Nitrous oxide emissions from intensive, fertilised agricultural systems have been identified as significant contributors to both Australia's and the global greenhouse gas (GHG) budget. This is expected to increase as rates of agriculture intensification and land use change accelerate to support population growth and food production. Limited data exists on N2O trace gas fluxes from subtropical or tropical tree cropping soils critical for the development of effective mitigation strategies.This study aimed to quantify GHG emissions over two consecutive years (March 2007 to March 2009) from a 30 year (lychee) orchard in the humid subtropical region of Australia. GHG fluxes were measured using a combination of high temporal resolution automated sampling and manually sampled chambers. No fertiliser was added to the plots during the 2007 measurement season. A split application of nitrogen fertiliser (urea) was added at the rate of 265kgNha-1 during the autumn and spring of 2008. Emissions of N2O were influenced by rainfall events and seasonal temperatures during 2007 and the fertilisation events in 2008. Annual N2O emissions from the lychee canopy increased from 1.7kgN2O-Nha-1yr-1 for 2007, to 7.6kgN2O-Nha-1yr-1 following fertiliser application in 2008. This represented an emission factor of 1.56%, corrected for background emissions. The timing of the split application was found to be critical to N2O emissions, with over twice as much lost following an application in spring (2.44%) compared to autumn (EF: 1.10%). This research suggests that avoiding fertiliser application during the hot and moist spring/summer period can reduce N2O losses without compromising yields.

Modeling nitrous oxide emissions from irrigated agriculture : testing DayCent with high-frequency measurements

A unique high temporal frequency dataset from an irrigated cotton-wheat rotation was used to test the agroecosystem model DayCent to simulate daily N2O emissions from sub-tropical vertisols under different irrigation intensities. DayCent was able to simulate the effect of different irrigation intensities on N2O fluxes and yield, although it tended to overestimate seasonal fluxes during the cotton season. DayCent accurately predicted soil moisture dynamics and the timing and magnitude of high fluxes associated with fertilizer additions and irrigation events. At the daily scale we found a good correlation of predicted vs. measured N2O fluxes (r2 = 0.52), confirming that DayCent can be used to test agricultural practices for mitigating N2O emission from irrigated cropping systems. A 25 year scenario analysis indicated that N2O losses from irrigated cotton-wheat rotations on black vertisols in Australia can be substantially reduced by an optimized fertilizer and irrigation management system (i.e. frequent irrigation, avoidance of excessive fertiliser application), while sustaining maximum yield potentials.

Crystalline Si nanoparticles below crystallization threshold : effects of collisional heating in non-thermal atmospheric-pressure microplasmas

Nucleation and growth of highly crystalline silicon nanoparticles in atmospheric-pressure low-temperature microplasmas at gas temperatures well below the Si crystallization threshold and within a short (100 μs) period of time are demonstrated and explained. The modeling reveals that collision-enhanced ion fluxes can effectively increase the heat flux on the nanoparticle surface and this heating is controlled by the ion density. It is shown that nanoparticles can be heated to temperatures above the crystallization threshold. These combined experimental and theoretical results confirm the effective heating and structure control of Si nanoparticles at atmospheric pressure and low gas temperatures.

Dense plasmas in magnetic traps : generation of focused ion beams with controlled ion-to-neutral flux ratios

Customized magnetic traps were developed to produce a domain of dense plasmas with a narrow ion beam directed to a particular area of the processed substrate. A planar magnetron coupled with an arc discharge source created the magnetic traps to confine the plasma electrons and generate the ion beam with the controlled ratio of ion-to-neutral fluxes. Images of the plasma jet patterns and numerical vizualizations help explaining the observed phenomena.

Ion impact distribution over plasma exposed nanocone arrays

The effect of an ordered array of nanocones on a conducting substrate immersed in the plasma on the transport of the plasma ions is investigated. The real conical shape of the cones is rigorously incorporated into the model. The movement of 10^5 CH3+ ions in the plasma sheath modified by the nanocone array is simulated. The ions are driven by the electric fields produced by the sheath and the nanostructures. The surface charge density and the total charge on the nanotips with different aspect ratios are computed. The ion transport simulation provides important characteristics of the displacement and velocity of the ions. The relative ion distribution along the lateral surfaces of the carbon nanotips is computed as well. It is shown that a rigorous account of the realistic nanostructure shape leads to very different distribution of the ion fluxes on the nanostructured surfaces compared to the previously reported works. The ion flux distribution is a critical factor in the nucleation process on the substrate and determines the nanostructure growth patterns.

Copper-capped carbon nanocones on silicon : plasma-enabled growth control

Controlled self-organized growth of vertically aligned carbon nanocone arrays in a radio frequency inductively coupled plasma-based process is studied. The experiments have demonstrated that the gaps between the nanocones, density of the nanocone array, and the shape of the nanocones can be effectively controlled by the process parameters such as gas composition (hydrogen content) and electrical bias applied to the substrate. Optical measurements have demonstrated lower reflectance of the nanocone array as compared with a bare Si wafer, thus evidencing their potential for the use in optical devices. The nanocone formation mechanism is explained in terms of redistribution of surface and volumetric fluxes of plasma-generated species in a developing nanocone array and passivation of carbon in narrow gaps where the access of plasma ions is hindered. Extensive numerical simulations were used to support the proposed growth mechanism.

Control of ion density distribution by magnetic traps for plasma electrons

The effect of a magnetic field of two magnetic coils on the ion current density distribution in the setup for low-temperature plasma deposition is investigated. The substrate of 400 mm diameter is placed at a distance of 325 mm from the plasma duct exit, with the two magnetic coils mounted symmetrically under the substrate at a distance of 140 mm relative to the substrate centre. A planar probe is used to measure the ion current density distribution along the plasma flux cross-sections at distances of 150, 230, and 325 mm from the plasma duct exit. It is shown that the magnetic field strongly affects the ion current density distribution. Transparent plastic films are used to investigate qualitatively the ion density distribution profiles and the effect of the magnetic field. A theoretical model is developed to describe the interaction of the ion fluxes with the negative space charge regions associated with the magnetic trapping of the plasmaelectrons. Theoretical results are compared with the experimental measurements, and a reasonable agreement is demonstrated.