Steepness expansion for free surface flows in coastal aquifers

Free surface flow of groundwater in aquifers has been studied since the early 1960s. Previous investigations have been based on the Boussinesq equation, derived from the non-linear kinematic boundary condition. In fact, the Boussinesq equation is the zeroth-order equation in the shallow-water expansion. A key assumption in this expansion is that the mean thickness of the aquifer is small compared with a reference length, normally taken to be the linear decay length. In this study, we re-examine the expansion scheme for free surface groundwater flows, and propose a new expansion wherein the shallow-water assumption is replaced by a steepness assumption. A comparison with experimental data shows that the new model provides a better prediction of water table levels than the conventional shallow-water expansion. The applicable ranges of the two expansions are exhibited. (c) 2004 Elsevier B.V. All rights reserved.

Swash overtopping and sediment overwash on a truncated beach

New experimental laboratory data are presented on swash overtopping and sediment overwash on a truncated beach, approximating the conditions at the crest of a beach berm or inter-tidal ridge-runnel. The experiments provide a measure of the uprush sediment transport rate in the swash zone that is unaffected by the difficulties inherent in deploying instrumentation or sediment trapping techniques at laboratory scale. Overtopping flow volumes are compared with an analytical solution for swash flows as well as a simple numerical model, both of which are restricted to individual swash events. The analytical solution underestimates the overtopping volume by an order of magnitude while the model provides good overall agreement with the data and the reason for this difference is discussed. Modelled flow velocities are input to simple sediment transport formulae appropriate to the swash zone in order to predict the overwash sediment transport rates. Calculations performed with traditional expressions for the wave friction factor tend to underestimate the measured transport. Additional sediment transport calculations using standard total load equations are used to derive an optimum constant wave friction factor of f(w)=0.024. This is in good agreement with a broad range of published field and laboratory data. However, the influence of long waves and irregular wave run-up on the overtopping and overwash remains to be assessed. The good agreement between modelled and measured sediment transport rates suggests that the model provides accurate predictions of the uprush sediment transport rates in the swash zone, which has application in predicting the growth and height of beach berms. (c) 2005 Elsevier B.V. All rights reserved.

Discussion of "Hydrodynamic loading on river bridges" by Stefano Malavasi and Alberto Guadagnini

Abstract not available.

Process optimization of an aerobic upflow sludge blanket reactor system

Response of an aerobic upflow sludge blanket (AUSB) reactor system to the changes in operating conditions was investigated by varying two principle operating variables: the oxygenation pressure and the flow recirculation rate. The oxygenation pressure was varied between 0 and 25 psig (relative), while flow recirculation rates were between 1,300 and 600% correspondingly. The AUSB reactor system was able to handle a volumetric loading of as high as 3.8 kg total organic carbon (TOC)/m(3) day, with a removal efficiency of 92%. The rate of TOC removal by AUSB was highest at a pressure of 20 psig and it decreased when the pressure was increased to 25 psig and the flow recirculation rate was reduced to 600%. The TOC removal rate also decreased when the operating pressure was reduced to 0 and 15 psig, with corresponding increase in flow recirculation rates to 1,300 and 1,000%, respectively. Maintenance of a high dissolved oxygen level and a high flow recirculation rate was found to improve the substrate removal capacity of the AUSB system. The AUSB system was extremely effective in retaining the produced biomass despite a high upflow velocity and the overall sludge yield was only 0.24-0.32 g VSS/g TOC removed. However, the effluent TOC was relatively high due to the system's operation at a high organic loading.

Wave-generated flow on coral reefs - an analysis for two-dimensional horizontal reef-tops with steep faces

Waves breaking on the seaward rim of a coral reef generate a flow of water from the exposed side of the reef to the sheltered side and/or to either channels through the reef-rim or lower sections of the latter. This wave-generated flow is driven by the water surface gradient resulting from the wave set-up created by the breaking waves. This paper reviews previous approaches to modelling wave-generated flows across coral reefs and discusses the influence of reef morphology and roughness upon these flows. Laboratory measurements upon a two-dimensional horizontal reef platform with a steep reef face provide the basis for extending a previous theoretical analysis for wave set-up on a reef in the absence of a flow [Gourlay, M.R., 1996b. Wave set-up on coral reefs. 2. Set-up on reefs with various profiles. Coastal Engineering 28, 1755] to include the interaction between a unidirectional flow and the wave set-up. The laboratory model results are then used to demonstrate that there are two basic reef-top flow regimes-reef-top control and reef-rim control. Using open channel flow theory, analytical relationships are derived for the reef-top current velocity in terms of the offreef wave conditions, the reef-top water depth and the physical characteristics of the reef-top topography. The wave set-up and wave-generated flow relationships are found to predict experimental values with reasonable accuracy in most cases. The analytical relationships are used to investigate wave-generated flows into a boat harbour channel on Heron Reef in the southern Great Barrier Reef. (c) 2005 Elsevier B.V. All rights reserved.

A knowledge-based system for liquid retaining structure design with blackboard architecture

Owing to the high degree of vulnerability of liquid retaining structures to corrosion problems, there are stringent requirements in its design against cracking. In this paper, a prototype knowledge-based system is developed and implemented for the design of liquid retaining structures based on the blackboard architecture. A commercially available expert system shell VISUAL RULE STUDIO working as an ActiveX Designer under the VISUAL BASIC programming environment is employed. Hybrid knowledge representation approach with production rules and procedural methods under object-oriented programming are used to represent the engineering heuristics and design knowledge of this domain. It is demonstrated that the blackboard architecture is capable of integrating different knowledge together in an effective manner. The system is tailored to give advice to users regarding preliminary design, loading specification and optimized configuration selection of this type of structure. An example of application is given to illustrate the capabilities of the prototype system in transferring knowledge on liquid retaining structure to novice engineers. (C) 2004 Elsevier Ltd. All rights reserved.

Sources of hydrogen sulfide in groundwater on reclaimed land

Fisherman Islands is an area of reclaimed land at the mouth of the Brisbane River in Queensland, Australia. Ongoing groundwater monitoring has found elevated concentrations of hydrogen sulfide (H2S) in the groundwater on the island. The presence of H2S on Fisherman Islands is of concern because of its toxic nature, the potential for acid sulfate soil formation, and its noxious odor. The aim of this study was,to, identify the sources of H2S within the groundwater on Fisherman Islands. It was hypothesized that the H2S is being formed by sulfate reducing bacteria acting on sulfate from seawater, rather than the introduction of sulfide with the dredge sediments. Groundwater and soil samples were collected and analyzed for sulfide, sulfate, and organic carbon contents. Elevated concentrations of sulfides coincide with,elevated concentrations of sulfate in the groundwater and elevated concentrations of organic carbon in the sediments, supporting the hypothesis that sulfide formation is the result of heterotrophic, sulfate reducing organisms.

Numerical solutions of the sediment conservation law: A review and improved formulation for coastal morphological modelling

Numerical solutions of the sediment conservation law are reviewed in terms of their application to bed update schemes in coastal morphological models. It is demonstrated that inadequately formulated numerical techniques lead to the introduction of diffusion, dispersion and the bed elevation oscillations previously reported in the literature. Four different bed update schemes are then reviewed and tested against benchmark analytical solutions. These include a first order upwind scheme, two Lax-Wendroff schemes and a non-oscillating centred scheme (NOCS) recently applied to morphological modelling by Saint-Cast [Saint-Cast, F., 2002. Modelisation de la morphodynamique des corps sableux en milieu littoral (Modelling of coastal sand banks morphodynamics), University Bordeaux 1, Bordeaux, 245 pp.]. It is shown that NOCS limits and controls numerical errors while including all the sediment flux gradients that control morphological change. Further, no post solution filtering is required, which avoids difficulties with selecting filter strength. Finally, NOCS is compared to a recent Lax-Wendroff scheme with post-solution filtering for a longer term simulation of the morphological evolution around a trained river entrance. (C) 2006 Elsevier B.V. All rights reserved.

Hydraulics of skimming flows on stepped chutes: The effects of inflow conditions?

Modern stepped spillways are typically designed for large discharge capacities corresponding to a skimming flow regime for which flow resistance is predominantly form drag. The writer demonstrates that the inflow conditions have some effect on the skimming flow properties. Boundary layer calculations show that the flow properties at inception of free-surface aeration are substantially different with pressurized intake. The re-analysis of experimental results highlights that the equivalent Darcy friction factor is f similar to 0.2 in average on uncontrolled stepped Chute and f similar to 0.1 on stepped chute with pressurized intake. A simple design chart is presented to estimate the residual flow velocity, and the agreement of the calculations with experimental results is deemed satisfactory for preliminary design.

Dam break wave of thixotropic fluid

Thixotropy is the characteristic of a fluid to form a gelled structure over time when it is not subjected to shearing, and to liquefy when agitated. Thixotropic fluids are commonly used in the construction industry (e.g., liquid concrete and drilling fluids), and related applications include some forms of mud flows and debris flows. This paper describes a basic study of dam break wave with thixotropic fluid. Theoretical considerations were developed based upon a kinematic wave approximation of the Saint-Venant equations down a prismatic sloping channel. A very simple thixotropic model, which predicts the basic theological trends of such fluids, was used. It describes the instantaneous state of fluid structure by a single parameter. The analytical solution of the basic flow motion and theology equations predicts three basic flow regimes depending upon the fluid properties and flow conditions, including the initial degree of jamming of the fluid (related to its time of restructuration at rest). These findings were successfully compared with systematic bentonite suspension experiments. The present work is the first theoretical analysis combining the basic principles of unsteady flow motion with a thixotropic fluid model and systematic laboratory experiments.

Estimation of water retention curve from mercury intrusion porosimetry and van Genuchten model

The water retention curve (WRC) is a hydraulic characteristic of concrete required for advanced modeling of water (and thus solute) transport in variably saturated, heterogeneous concrete. Unfortunately, determination by a direct experimental method (for example, measuring equilibrium moisture levels of large samples stored in constant humidity cells) is a lengthy process, taking over 2 years for large samples. A surrogate approach is presented in which the WRC is conveniently estimated from mercury intrusion porosimetry (MIP) and validated by water sorption isotherms: The well-known Barrett, Joyner and Halenda (BJH) method of estimating the pore size distribution (PSD) from the water sorption isotherm is shown to complement the PSD derived from conventional MIP. This provides a basis for predicting the complete WRC from MIP data alone. The van Genuchten equation is used to model the combined water sorption and MIP results. It is a convenient tool for describing water retention characteristics over the full moisture content range. The van Genuchten parameter estimation based solely on MIP is shown to give a satisfactory approximation to the WRC, with a simple restriction on one. of the parameters.

An advanced regularization methodology for use in watershed model calibration

A calibration methodology based on an efficient and stable mathematical regularization scheme is described. This scheme is a variant of so-called Tikhonov regularization in which the parameter estimation process is formulated as a constrained minimization problem. Use of the methodology eliminates the need for a modeler to formulate a parsimonious inverse problem in which a handful of parameters are designated for estimation prior to initiating the calibration process. Instead, the level of parameter parsimony required to achieve a stable solution to the inverse problem is determined by the inversion algorithm itself. Where parameters, or combinations of parameters, cannot be uniquely estimated, they are provided with values, or assigned relationships with other parameters, that are decreed to be realistic by the modeler. Conversely, where the information content of a calibration dataset is sufficient to allow estimates to be made of the values of many parameters, the making of such estimates is not precluded by preemptive parsimonizing ahead of the calibration process. White Tikhonov schemes are very attractive and hence widely used, problems with numerical stability can sometimes arise because the strength with which regularization constraints are applied throughout the regularized inversion process cannot be guaranteed to exactly complement inadequacies in the information content of a given calibration dataset. A new technique overcomes this problem by allowing relative regularization weights to be estimated as parameters through the calibration process itself. The technique is applied to the simultaneous calibration of five subwatershed models, and it is demonstrated that the new scheme results in a more efficient inversion, and better enforcement of regularization constraints than traditional Tikhonov regularization methodologies. Moreover, it is argued that a joint calibration exercise of this type results in a more meaningful set of parameters than can be achieved by individual subwatershed model calibration. (c) 2005 Elsevier B.V. All rights reserved.