Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.

Experimental investigation of contaminant transport in coastal groundwater

Contaminant transport in coastal aquifers is of increasing interest since, with the development of coastal areas, contaminants from surface sources may enter coastal aquifers and pollute the groundwater flow. Coastal groundwater flow is complicated because of the presence of a freshwater-saltwater diffusion zone and the tidal variation of sea level at the seaward end. This paper investigates experimentally the behaviour of contaminant plumes with different densities in an unconfined coastal aquifer. Experiments were performed in a flow tank filled with glass beads as the porous medium. Results show that the dense contaminant has a more diffusive front than the less dense one in the seaward direction towards the coastline. The plume becomes more diffusive when it travels closer to the saltwater interface. On the contrary, the less dense contaminant presents a relatively sharp outline. It tends to migrate in the upper portion of the aquifer and exits in a concentrated manner over a small discharge area at the coastline, not further seaward under the sea. Non-dimensional parameters show that instabilities occur in our experiments for a density difference of 1.2% or larger between the contaminant and the ambient water. The experimental results provide guidance for field monitoring and numerical modelling. (C) 2002 Elsevier Science Ltd. All rights reserved.

Numerical modelling of contaminant transport in coastal aquifers

This paper employs a two-dimensional variable density flow and transport model to investigate the transport of a dense contaminant plume in an unconfined coastal aquifer. Experimental results are also presented to show the contaminant plume in a freshwater-seawater flow system. Both the numerical and experimental results suggest that the neglect of the seawater interface does not noticeably affect the horizontal migration rate of the plume before it reaches the interface. However, the contaminant will travel further seaward and part of the solute mass will exit under the sea if the higher seawater density is not included. If the seawater density is included, the contaminant will travel upwards towards the beach along the freshwater-saltwater interface as shown experimentally. Neglect of seawater density, therefore, will result in an underestimate of solute mass rate exiting around the coastline. (C) 2002 IMACS. Published by Elsevier Science B.V. All rights reserved.

Confined-unconfined flow in a horizontal confined aquifer

Mixed confined and unconfined groundwater flow occurs in a bounded initially dry aquifer when the hydraulic head at the side boundary suddenly rises above the elevation of the aquifer's top boundary. The flow problem as modelled by the Boussinesq equation is non-trivial because of the involvement of two moving boundaries. The transformed equation (based on a similarity transformation) can, however, be dealt with more easily. Here, we present an approximate analytical solution for this flow problem. The approximate solution is compared with an 'exact' numerical solution and found to be a very accurate description for describing the mixed confined and unconfined flow in the confined aquifer. (C) 2002 Elsevier Science B.V. All rights reserved.

Influence of hysteresis on tidal capillary fringe dynamics in a well-sorted sand

Nielsen and Perrochet [Adv. Water Resour. 23 (2000) 503] presented experimental data for cyclic water movement in the vadose zone above an oscillating watertable. The response of the watertable to cyclic forcing was characterised by the ratios of the forcing head to watertable amplitudes and their associated phase lag. They found that their non-hysteretic Richards' equation model failed to represent the observed behaviour of these parameters. This paper explores the effect on the simulated capillary fringe dynamics (in terms of these parameters) of including varying degrees of hysteresis in the moisture retention curve used in a numerical model of their experiment. It is clear that hysteresis can indeed account for observed discrepancies between simulation and experiment and that the effect of hysteresis varies with the frequency of oscillation. The use of a single-valued mean retention curve, as advocated by some authors, fails to provide a match between the simulated and observed behaviour of the Nielsen and Perrochet parameters, but is shown to be adequate for predicting time-averaged soil moisture profiles. (C) 2003 Elsevier Ltd. All rights reserved.

The application of wave induced forces to a two-dimensional finite element long wave hydrodynamic model

This paper describes the modification of a two-dimensional finite element long wave hydrodynamic model in order to predict the net current and water levels attributable to the influences of waves. Tests examine the effects of the application of wave induced forces, including comparisons to a physical experiment. An example of a real river system is presented with comparisons to measured data, which demonstrate the importance of simulating the combined effects of tides and waves upon hydrodynamic behavior. (C) 2002 Elsevier Science Ltd. All rights reserved.

Modelling the effect of shear history on activated sludge flocculation

The aim of this paper is to investigate the effect of shear history on activated sludge flocculation dynamics and to model the observed relationships using population balances. Activated sludge flocs are exposed to dramatic changes in the shear rate within the treatment process, as they pass through localised high and low mixing intensities within the aeration basin and are cycled through the different unit operations of the treatment process. We will show that shear history is a key factor in determining floc size, and that the floc size varies irreversibly with changes in shear rate. A population balance model of the flocculation process is also introduced and evaluated.

Biotechnological production of lactic acid integrated with potato wastewater treatment by Rhizopus arrhizus

This paper describes a feasibility study of a for lactic acid production integrated with are treatment of wastewater from an industrial starch plant. Rhizopus oryzae two strains, Rhizopus arrhizus and Rhizopus oligosporus were tested with respect to their capability to carry out simultaneous saccharification and fermentation to lactic acid using potato wastewater. Rhizopus arrhizus DAR 36017 was identified as a suitable strain that demonstrated a high capacity for starch saccharification and lactic acid synthesis. The optimal conditions, in terms of pH, temperature and starch concentration, for lactic acid production were determined. The selected fungal strain grew well in a pH range from 3.0 to 7.0. The addition of CaCO(3)10 g dm(-3) maintained the pH at 5.0-6.0 and significantly enhanced lactic acid production. Kinetic study revealed that almost complete starch saccharification and a lactic acid yield of 450g kg(-1) could be achieved in 20 h and 28 h cultivation, respectively. The maximum lactic acid production 21 g dm(-3) and mycelial biomass (1.7 g dm(-3)) were obtained at 30degreesC. Besides the multiple bioproducts, total removal of suspended solids and 90% reduction of COD were achieved in a single no-aseptic operation. (C) 2003 Society of Chemical Industry.

Settling velocity of sediments at high concentrations

New data on the settling velocity of artificial sediments and natural sands at high concentrations are presented. The data are compared with a widely used semiempirical Richardson and Zaki equation (Trans. Inst. Chem. Eng. 32 (1954) 35), which gives an accurate measure of the reduction in velocity as a function of concentration and an experimentally determined empirical power n. Here, a simple method of determining n is presented using standard equations for the clear water settling velocity and the seepage flow within fixed sediment beds. The resulting values for n are compared against values derived from new and existing laboratory data for beach and filter sands. For sands, the appropriate values of n are found to differ significantly from those suggested by Richardson and Zaki for spheres, and are typically larger, corresponding to a greater reduction in settling velocity at high concentrations. For fine and medium sands at concentrations of order 0.4, the hindered settling velocity reduces to about 70% of that expected using values of n derived for spheres. At concentrations of order 0.15, the hindered settling velocity reduces to less than half of the settling velocity in clear water. These reduced settling velocities have important implications for sediment transport modelling close to, and within, sheet flow layers and in the swash zone.

Undular tidal bore dynamics in the Daly Estuary, Northern Australia

Measurements in the macro-tidal Daly Estuary show that the presence of an undular tidal bore contributed negligibly to the dissipation of tidal energy. No recirculation bubble was observed between a trough and the following wave crest in the lee waves following the undular bore. This differs to stationary undular bores in laboratory experiments at larger Froude numbers where a recirculation bubble exists. Secondary motions and the turbulence generated by the undular bore had no measurable influence on the sediment transport. This situation contrasts with the intense sediment resuspension observed in breaking tidal bores. The tidally averaged sediment budget in the Daly Estuary was controlled by the asymmetry of tidal currents. The undular bore may widen the river by breaking along the banks that it undercuts, leading to bank slippage. A patch of river-wide macro-turbulence of 3-min duration occurred about 20 min after the passage of the bore during accelerating tidal currents. (C) 2004 Elsevier Ltd. All rights reserved.

Effect of initial disturbance amplitude in gravity affected jet break-up

The initial disturbance amplitude has an effect on stretching jets that is not observed for capillary jet instability where gravitational acceleration is not significant. For inviscid and viscous fluids, gravity diminishes the effect that the initial amplitude has on jet length and its ability to prevent satellite formation. In stretching jets, not only the dimensionless frequency of the disturbance but also its initial amplitude must be known to properly study their satellite forming nature. Indirect methods of relating the applied disturbance energy to an initial velocity perturbation are not simple when the gravity parameter G is changing. When G A 0, the optimum disturbance frequency Omega(opt) and the initial disturbance amplitude are related, with Omega(opt) proportional to f (G) x In(1 /epsilon(nu)). Results from numerical simulations and experiments are presented here. (c) 2005 Elsevier Ltd. All rights reserved.

Structure of a cellulose degrading bacterial community during anaerobic digestion

It is widely accepted that cellulose is the rate-limiting substrate in the anaerobic digestion of organic solid wastes and that cellulose solubilisation is largely mediated by surface attached bacteria. However, little is known about the identity or the ecophysiology of cellulolytic microorganisms from landfills and anaerobic digesters. The aim of this study was to investigate an enriched cellulolytic microbial community from an anaerobic batch reactor. Chemical oxygen demand balancing was used to calculate the cellulose solubilisation rate and the degree of cellulose solubilisation. Fluorescence in situ hybridisation (FISH) was used to assess the relative abundance and physical location of three groups of bacteria belonging to the Clostridium lineage of the Firmicutes that have been implicated as the dominant cellulose degraders in this system. Quantitation of the relative abundance using FISH showed that there were changes in the microbial community structure throughout the digestion. However, comparison of these results to the process data reveals that these changes had no impact on the cellulose solubilisation in the reactor. The rate of cellulose solubilisation was approximately stable for much of the digestion despite changes in the cellulolytic population. The solubilisation rate appears to be most strongly affected by the rate of surface area colonisation and the biofilm architecture with the accepted model of first order kinetics due to surface area limitation applying only when the cellulose particles are fully covered with a thin layer of cells. (c) 2005 Wiley Periodicals, Inc.