Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Photoinhibition and photoprotection in symbiotic dinoflagellates from reef-building corals

Pulse-amplitude-modulation fluorometry and oxygen respirometry were used to investigate diel photosynthetic responses by symbiotic dinoflagellates to light levels in summer and winter on a high latitude coral reef. The symbiotic dinoflagellates from 2 species of reef-building coral (Porites cylindrica and Stylophora pistillata) showed photoinhibitory decreases in the ratio of variable (F-v) to maximal (F-m) fluorescence (F-v/F-m) as early as 09:00 h on both summer and winter days on the reefs associated with One Tree Island (23 degrees 30' S, 152 degrees 06' E; Great Barrier Reef, Australia). This was due to decreases in maximum, F-m, and to a smaller extent minimum, F-0, chlorophyll fluorescence. Complete recovery took 4 to 6 h and began to occur as soon as light levels fell each day. Chlorophyll fluorescence quenching analysis of corals measured during the early afternoon revealed classic regulation of photosystem II (PSII) efficiency through non-photochemical quenching (NPQ). These results appear to be similar to data collected for other algae and higher plants, suggesting involvement of the xanthophyll cycle of symbiotic dinoflagellates in regulating the quantum efficiency of PSII. The ability of symbiotic dinoflagellates to develop significant NPQ, however, depended strongly on when the symbiotic dinoflagellates were studied. Whereas symbiotic dinoflagellates from corals in the early afternoon showed a significant capacity to regulate the efficiency of PSII using NPQ, those sampled before sunrise had a slower and much reduced capacity, suggesting that elements of the xanthophyll cycle are suppressed prior to sunrise. A second major finding of this study is that the quantum efficiency of PSII in symbiotic dinoflagellates is strongly diurnal, and is as much as 50% lower just prior to sunrise than later in the day. When combined with oxygen flux data, these results indicate that a greater portion of the electron transport occurring later in the day is likely to be due to the increases in the rate of carbon fixation by Rubisco or to higher flutes through the Mehler-Ascorbate-Peroxidase (MAP) cycle.

Quantifying the influence of intra-aggregate concentration gradients on solute transport

The physical nonequilibrium of solute concentration resulting from preferential now of soil water has often led to models where the soil is partitioned into two regions: preferential flow paths, where solute transport occurs mainly by advection, and the remaining region, where significant solute transport occurs through diffusive exchange with the flow paths. These two-region models commonly ignore concentration gradients within the regions. Our objective was to develop a simple model to assess the influence of concentration gradients on solute transport and to compare model results with experiments conducted on structured materials. The model calculates the distribution of solutes in a single spherical aggregate surrounded by preferential now paths and subjected to alternating boundary conditions representing either an exchange of solutes between the two regions (a wet period) or no exchange but redistribution of solutes within the aggregate (a dry period). The key parameter in the model is the aggregate radius, which defines the diffusive time scales. We conducted intermittent leaching experiments on a column of packed porous spheres and on a large (300 mm long by 216 mm diameter) undisturbed field soil core to test the validity of the model and its application to field soils. Alternating wet and dry periods enhanced leaching by up to 20% for this soil, which was consistent with the model's prediction, given a fitted equivalent aggregate radius of 1.8 cm, If similar results are obtained for other soils, use of alternating wet and dry periods could improve management of solutes, for example in salinity control and in soil remediation.

Comparison of coherent and weakly incoherent transport models for the interlayer magnetoresistance of layered Fermi liquids

The interlayer magnetoresistance of layered metals in a tilted magnetic field is calculated for two distinct models for the interlayer transport. The first model involves coherent interlayer transport, and makes use of results of semiclassical or Bloch-Boltzmann transport theory. The second model involves weakly incoherent interlayer transport where the electron is scattered many times within a layer before tunneling into the next layer. The results are relevant to the interpretation of experiments on angular-dependent magnetoresistance oscillations (AMRO) in quasi-one- and quasi-two-dimensional organic metals. We find that the dependence of the magnetoresistance on the direction of the magnetic field is identical for both models except when the field is almost parallel to the layers. An important implication of this result is that a three-dimensional Fermi surface is not necessary for the observation of the Yamaji and Danner oscillations seen in quasi-two- and quasi-one-dimensional metals, respectively. A universal expression is given for the dependence of the resistance at AMRO maxima and minima on the magnetic field and scattering time (and thus the temperature). We point out three distinctive features of coherent interlayer transport: (i) a beat frequency in the magnetic oscillations of quasi-two-dimensional systems, (ii) a peak in the angular-dependent magnetoresistance when the field is sufficiently large and parallel to the layers, and (iii) a crossover from a linear to a quadratic field dependence for the magnetoresistance when the field is parallel to the layers. Properties (i) and (ii) are compared with published experimental data for a range of quasi-two-dimensional organic metals. [S0163-1829(99)02236-5].

CXTANNEAL: an improved program for estimating solute transport parameters

CXTANNEAL is a program for analysing contaminant transport in soils. The code, written in Fortran 77, is a modified version of CXTFIT, a commonly used package for estimating solute transport parameters in soils. The improvement of the present code is that it includes simulated annealing as the optimization technique for curve fitting. Tests with hypothetical data show that CXTANNEAL performs better than the original code in searching for optimal parameter estimates. To reduce the computational time, a parallel version of CXTANNEAL (CXTANNEAL_P) was also developed. (C) 1999 Elsevier Science Ltd. All rights reserved.

Which one of these is not like the others? An inter-hospital study of the viscosity of thickened fluids

Engineering This investigation examined the rheological (viscosity and yield stress) and material property (density) characteristics of the thickened meal-time and videofluorscopy fluids provided by 10 major metropolitan hospitals. Differences in the thickness of thickened fluids were considered as a source of variability and potential hazard for inter-hospital transfers of dysphagic patients. The results indicated considerable differences in the viscosity, density, and yield stress of both meal-time and videofluoroscopy fluids. In theory, the results suggest that dysphagic patients transferred between hospitals could be placed on inappropriate levels of fluid thickness because of inherent differences in the rheology and material property characteristics of the fluids provided by different hospitals. Slowed improvement or medical complications are potential worst-case scenarios for dysphagic patients if the difference between the thick fluids offered by 2 hospitals are extreme. The investigation outlines the most appropriate way to assess the rheological and material property characteristics of thickened fluids. In addition, it suggests a plan of quality improvement to reduce the variability of the thickness of fluids offered at different hospitals.

Transport properties of strongly correlated metals: A dynamical mean-field approach

The temperature dependence of the transport properties of the metallic phase of a frustrated Hubbard model on the hypercubic lattice at half-filling is calculated. Dynamical mean-held theory, which maps the Hubbard model onto a single impurity,Anderson model that is solved self-consistently, and becomes exact in the limit of large dimensionality, is used. As the temperature increases there is a smooth crossover from coherent Fermi liquid excitations at low temperatures to incoherent excitations at high temperatures. This crossover leads to a nonmonotonic temperature dependence for the resistance, thermopower, and Hall coefficient, unlike in conventional metals. The resistance smoothly increases from a quadratic temperature dependence at low temperatures to large values which can exceed the Mott-Ioffe-Regel value ha/e(2) (where a is a lattice constant) associated with mean free paths less than a lattice constant. Further signatures of the thermal destruction of quasiparticle excitations are a peak in the thermopower and the absence of a Drude peak in the optical conductivity. The results presented here are relevant to a wide range of strongly correlated metals, including transition metal oxides, strontium ruthenates, and organic metals.

Solute transport by surface runoff from low-angle slopes: theory and application

The removal of chemicals in solution by overland how from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface-mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff: to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright (C) 2000 John Wiley & Sons, Ltd.

Experimental examination of solute transport by surface runoff from low-angle slopes

The removal of chemicals in solution, by overland flow from agricultural land has the potential to be a significant source of chemical loss from zero-till and surface mulched farming systems. The objective of this study was to determine the magnitude of solute loss by surface runoff from agricultural systems. Previous experiments have enhanced the understanding of the exchange process, but the initial soil conditions together with the tracer application method in these experiments have meant that in some cases the results have limited applicability to field situations. In this study, two different sets of experiments were carried out to determine the magnitude of solute loss by surface runoff. These experiments entailed the surface application of bromide to (1) field scale plots 18 m long by 2 m wide and (2) repacked soil cores 236 mm in diameter; followed by the application of simulated rainfall in both cases. The most substantial finding of the field experiments was that the quantities of solute in surface runoff varied greatly with soil type and structure (0.07-14.9% of the applied bromide). Also, on some soils, large quantities of tracer were measured in the surface runoff even after several hours of infiltration. The experiments on soil cores showed that soil structure plays an important role in the quantity of chemical that may be transported in the surface runoff. These field results showed that, in certain systems, solute movement by overland flow is an important transport mechanism, which should be considered when budgeting for chemical loss. (C) 2000 Elsevier Science B.V. All rights reserved.

Identifying regional dust transport pathways: Application of kinematic trajectory modelling to a trans-Tasman case

PI kinematic trajectory model is used to investigate potential pathways of dust transport from Australia to New Zealand. Historically, these have been assumed to follow rather direct west-east trajectories spanning 2 to 3 days, often resulting in red snow events in the Southern Alps of New Zealand. However, results from the present study which examined the route taken by air parcels originating in southern Australia during dust storms on 24 and 25 May 1994, indicate that trans-Tasman dust transport trajectories are more diverse than previously thought, and display considerable variation during single events. These mon divergent pathways tie in more closely with aeolian dust sedimentation patterns identified by ocean coring in the Tasman Sea, and may account for the deposition of Australian dust on sub-Antarctic islands located well south of the Australian continent. Copyright 2000 John Wiley Sons, Ltd.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Labelled tableaux for non-normal modal logics

In this paper we show how to extend KEM, a tableau-like proof system for normal modal logic, in order to deal with classes of non-normal modal logics, such as monotonic and regular, in a uniform and modular way.

BeachWin: modelling groundwater effects on swash sediment transport and beach profile changes

Field and laboratory observations have shown that a relatively low beach groundwater table enhances beach accretion. These observations have led to the beach dewatering technique (artificially lowering the beach water table) for combating beach erosion. Here we present a process-based numerical model that simulates the interacting wave motion on the beach. coastal groundwater flow, swash sediment transport and beach profile changes. Results of model simulations demonstrate that the model replicates accretionary effects of a low beach water table on beach profile changes and has the potential to become a tool for assessing the effectiveness of beach dewatering systems. (C) 2002 Elsevier Science Ltd. All rights reserved.