Comparison of Models for Computing Drainage Discharge

The WAVE model describes the transport and transformations of matter and energy in the soil, crop, and vadose environment. A lateral field drainage subprogram was added to the WAVE model to simulate lateral subsurface drainage flow. The subsurface drainage is considered as the drainage provided by evenly spaced parallel drains with a free outlet: drain tubing or ditch. The rate of subsurface water movement into drain tubes or ditches depends on the hydraulic conductivity of the soil, drain or ditch spacing, hydraulic head in the drains, profile depth, and water table elevation. Hooghoudt's steady-state equation was selected for incorporation in the WAVE model. The subsurface drainage subprogram was calibrated and validated by comparison with the SWAP model (The Netherlands) and DRAINMOD (the United States) and partially by using 7 years of drain outflow data from an experimental field under fallow and cropped conditions. The comparative study revealed that the three models performed equally well and that the models were reliable and accurate tools for predicting the drainage flux as a function of rainfall-evapotranspiration and local conditions. The WAVE model, in comparison to the SWAP and DRAINMOD model, provided as good a prediction of the lateral subsurface drainage flow to drains. The statistical analysis between each model and observed data revealed that the three models were able to predict with sufficient accuracy the observed drainage discharge. The DRAINMOD model, however, has the advantage of giving a more accurate estimate of the discharge, resulting in a more precise modeling. The models were consistent in predicting water table levels, but they could not be verified against field data because of a lack of suitable measurements.     

Full-Flowing Collection Conduits with Nonuniform Inflow

Collection conduits flowing full with nonuniform inflow (variable rate of increase in flow with distance along the conduit length) include well screens (vertical and directionally drilled); submerged effluent collectors; and certain types of inboard weir configurations for settling tanks. Certain subsurface drains used in environmental engineering applications and civil engineering more generally may be inadvertently designed for full-flowing conditions. Formulation of the problem for such collection conduits is presented in terms of the applicable differential equation, slope invariance, the frictionless solution for a general cross section, a uniform-inflow solution, and the difference formulation. The importance of checking inflow uniformity is discussed and exemplified. Dimensional analysis is then employed to demonstrate the relationship between variables, leading to a new generalized numerical solution. That solution is presented in a graphical form which provides further useful display of the relationships between variables and quantitative information for design and analysis. The detrimental effects of flowing-full conditions in subsurface drains is demonstrated, and it is noted that existing design methods may unintentionally cause such conditions to occur.     

Field Study to Investigate WIDE Technology for TCE Extraction

A full-scale field study was conducted using well injection depth extraction (WIDE) technology to remove trichloroethylene (TCE) from subsurface profiles with fine-grained soils. WIDE incorporates the use of geosynthetic wick drains, attached to an aboveground PVC pipe network, designed to (1) extract contaminated fluids from a specific depth using vacuum pressure, and/or (2) inject flushing solutions. Both extraction-only and concurrent injection-extraction testing schemes were conducted, and fluid flow rates, TCE recovery rates, and groundwater elevations were monitored over a 9?month time interval. During extraction-only operational schemes, gas-phase TCE extraction rates were significantly higher than liquid-phase extraction rates due to the increased volume of air within the geosynthetic wells and PVC piping. TCE extraction rates were less than 2,000?mg/h for airflow rates less than 100,000?L/h, and increased significantly to 5,700?mg/h as the airflow rate approached 600,000?L/h. Long term testing and sampling is needed to quantify system performance and determine tailing and rebounding effects in comparison to conventional methods. The behaviors and trends observed during the field study are presented and discussed.     

Design and Management of Subsurface Horizontal Drainage to Reduce Salt Loads

Many irrigated areas have shallow water tables creating waterlogging and salinization problems. This has often been controlled by installation of subsurface horizontal pipe drainage; however, these systems export large amounts of salt off farm in the drainage effluent. Improved design and management of subsurface drainage systems to reduce drainage salt loads were tested in a replicated field experiment. Deep, widely spaced drains allowed to flow without control were compared to drains with management to reduce drain flow. These were also compared with shallow, closely spaced drains that protected the root zone only and an undrained control. The deep drains flowed continuously during the two irrigation seasons with an electrical conductivity of around 11 dS∕m resulting in a drainage salt load of 5,867 kg∕ha. The management measures reduced drainage volume and salinity resulting in a 50% reduction in salt load. The shallow drains only flowed directly after an irrigation or rainfall event with low salinity, around 2 dS∕m, resulting in a 95% reduction in salt load. This showed that by management there is great potential for reducing salt mobilization in existing drainage systems, and for new systems shallower drains will minimize salt loads.     

New Elastic Model of Pipe Flow for Stability Analysis of the Governor-Turbine-Hydraulic System

A higher-order elastic model of the flow in long pressurized pipelines is expected to be utilized for stability analysis of the governor-turbine-hydraulic system in hydropower stations. Because traditional elastic models are limited in lower order application because of their difficult decoupling in addition to the rigid model, a new linear elastic model of the flow in pressurized pipelines is derived on the basis of the equations of hydraulic vibration, in which each oscillatory flow with a different order has been obtained with ordinary differential equations in decoupling form. For water conveyance systems with branching pipes or parallel pipes in hydropower stations, the state equations to describe hydraulic characteristics of the governor-turbine-hydraulic system are established with the application of this new elastic model for diversion pipeline flow or tail tunnel flow. The influence of the elastic models with different order on a system’s stability are revealed in detail by two cases that illustrate that an elastic model with proper order should be used for the flow in pressurized pipelines of hydropower stations, according to their length, to improve the accuracy of stability analysis.     

Design of Two-Layered Porous Landscaping Detention Basin

Under the mandate of the Federal Clean Water Act, porous landscaping detention (PLD) has been widely used to increase on-site infiltration. A PLD system consists of a surface storage basin and subsurface filtering layers. The major design parameters for a PLD system are the infiltration rate on the land surface and the seepage rate through the subsurface medium. A low infiltration rate leads to a sizable storage basin while a high infiltration rate results in standing water if the subsurface seepage does not sustain the surface loading. In this study, the design procedure of a PLD basin is revised to take both detention flow hydrology and seepage flow hydraulics into consideration. The design procedure begins with the basin sizing according to the on-site water quality control volume. The ratio of design infiltration rate to sand-mix hydraulic conductivity is the key factor to select the thickness of sand-mix layer underneath a porous bed. The total filtering thickness for both sand-mix and gravel layers is found to be related to the drain time and infiltration rate. The recommended sand-mix and granite gravel layers underneath a PLD basin are reproduced in the laboratory for infiltration tests. The empirical decay curve for sand-mix infiltration rate was derived from the laboratory data and then used to maximize the hydraulic efficiency through the subsurface filtering layers. In this study, it is recommended that a PLD system be designed with the optimal performance to consume the hydraulic head available and then evaluated using the prolonged drain time for potential standing water problems under various clogging conditions.     

Recursive Design of Pressurized Branched Irrigation Networks

This paper presents a new method to design pressurized branched irrigation networks. This method is called recursive design and is based on application of the problem-solving technique known as backtracking to the problem of the optimum design of pressurized branched irrigation networks with a known delivery piezometric head (pipe-sizing). Recursive design is a heuristic optimizer, like genetic algorithms, and has been implemented in a fast, versatile computer application. After presenting and precisely defining the design problem, the writers review the theoretical foundations of some of the main existing design methods: maximum velocity, recommended velocity, Mougnie velocity, constant hydraulic slope, Lagrange multipliers, linear programming, Labye’s method, and genetic algorithms. Next, the writers explain what recursive design consists of and apply its methodology in detail to a simple network. In the results section, the solutions obtained by recursive design are compared with those obtained by the other design methods, giving satisfactory results. For example, in an analyzed standard network, genetic algorithms take more than 20?minutes to offer a solution, whereas recursive design offers a cheaper solution with less than 3?seconds of computation time.     

Improved Simulation of Flow Regime Transition in Sewers: Two-Component Pressure Approach

Operational problems and system damage have been linked to the flow regime transition between free surface and pressurized flow in rapidly filling stormwater and combined sewer systems. In response, emphasis has been placed on the development of numerical models to describe hydraulic bores and other flow phenomena that may occur in these systems. Current numerical models are based on rigid column analyses, shock-fitting techniques, or shock-capturing procedures employing the Preissmann slot concept. The latter approach is appealing due to the comparative simplicity, but suffers from the inability to realistically describe subatmospheric full-pipe flows. A new modeling framework is proposed for describing the flow regime transition utilizing a shock-capturing technique that decouples the hydrostatic pressure from surcharged pressures occurring only in pressurized conditions, effectively overcoming the cited Preissmann slot limitation. This new approach exploits the identity between the unsteady incompressible flow equations for elastic pipe walls and the unsteady open-channel flow equations, and the resulting numerical implementation is straightforward with only minor modifications to standard free surface flow models required. A comparison is made between the model predictions and experimental data; good agreement is achieved.     

Experimental Study of a Right-Angled End Junction between a Pipe and an Open Channel

This paper presents an experimental study of a junction between a closed conduit and an open channel. This study was undertaken to explore hydraulic properties of outlets of subsurface drainage or sewage networks into an open air stream during flood events. Experiments were conducted in a laboratory flume, with a main rectangular channel joined at right angle to a lateral circular pipe. Both branches were supplied with independent flow rates and downstream water level was controlled by an adjustable weir. Several flow patterns were identified, combining free-surface and pressurized flows. Transitions between these flow patterns, as well as changes in water level or energy, in response to the modifications of experimental variables, were studied and could be linked to known properties of single channels, single pipes, and homogeneous junctions. Transitions between free-surface and pressurized pipe flow appeared to be strongly dependent on the whole set of experimental variables and the pipe longitudinal slope. This work contributes to a better knowledge of hydraulic and hydrologic key processes for point source discharging.     

Drainage of Sloping Lands with Variable Recharge Model Validation and Applications

Semianalytical transient equations describing the behavior of water tables in subsurface drained soils when drains rest on a sloping impervious barrier have been derived previously and represent the bases of the SIDRA model. To validate (SIDRA), water table elevations and drainflow rates have been monitored for six years in the French Alps on fields with a slope of 8%. The predicted drain flow rates and water table shapes compare reasonably well with data of a drainage experiment site but the improvement provided by taking the slope into consideration was limited. Running the model with different slopes confirmed that high water table elevations and peak flows were not significantly changed with a slope of 8%. As could be predicted from the analysis in steady state, low water table elevations were the most affected by slope. With the soil parameters of the field experiment and from an analysis in tail recession conditions, it was shown that there is a clear threshold of 12% slope below which the slope has no significant influence and can be neglected in drainage design.     

Collection Conduits Including Subsurface Drains

Numerous types of pipes and channels with spatially increasing flows in environmental engineering applications are identified by type and function and referred to as collection conduits. An overview of methods for designing and analyzing collection conduits is provided. Full conduits with nonuniform and uniform inflow are first considered. Dimensional analysis is then employed to demonstrate the relationship between variables for open channels; that leads to the identification of possibilities for generalized numerical solutions. Prior collection conduit applications are discussed within the framework of the dimensional analysis (which also pertains to some constant-flow applications). A previously unpublished generalized numerical solution for rectangular collection conduits is presented. Subsurface drains are addressed with particular emphasis, including the use of numerical methods to develop a new generalized chart and relation to other design methods. Among the important conclusions for subsurface drains is that the somewhat common practice of using Manning’s equation alone for such problems is not generally adequate. Examples and practical design suggestions are included, and the use of computer-based numerical methods is discussed more generally.     

Coupled Mechanical and Hydraulic Modeling of Geosynthetic-Reinforced Column-Supported Embankments

Geosynthetic-reinforced column-supported (GRCS) embankments have increasingly been used in the recent years for accelerated construction. Numerical analyses have been conducted to improve understanding and knowledge of this complicated embankment system. However, most studies so far have been focused on its short-term or long-term behavior by assuming an undrained or drained condition, which does not consider water flow in saturated soft soil (i.e., consolidation). As a result, very limited attention has been paid to a settlement-time relationship especially postconstruction settlement, which is critical to performance of pavements on embankments or connection between approach embankments and bridge abutments. To investigate the time-dependent behavior, coupled two-dimensional mechanical and hydraulic numerical modeling was conducted in this study to analyze a well-instrumented geotextile-reinforced deep mixed column-supported embankment in Hertsby, Finland. In the mechanical modeling, soils and DM columns were modeled as elastic-plastic materials and a geotextile layer was modeled using cable elements. In the hydraulic modeling, water flow was modeled to simulate generation and dissipation of excess pore water pressures during and after the construction of the embankment. The numerical results with or without modeling water flow were compared with the field data. In addition, parametric studies were conducted to further examine the effects of geosynthetic stiffness, column modulus, and average staged construction rate on the postconstruction settlement and the tension in the geosynthetic reinforcement.     

Hydraulics of Tangential Vortex Intake for Urban Drainage

A tangential vortex intake is a compact structure that can convey storm water efficiently as a swirling flow down a vortex dropshaft. It has been studied in physical models and successfully employed in urban drainage and hydroelectric plant applications, but a comprehensive account of the key flow characteristics has not been reported and a theoretical design guideline of a tangential intake is not available. In this study the hydraulics of tangential slot vortex intakes is investigated via extensive experiments. It is found that the flow in the tapering and downward sloping vortex inlet channel is strongly dependent on the geometry of the inlet and dropshaft. Under some conditions, hydraulic instability and overflow can occur, rendering the design ineffective. It is shown that the hydraulic stability depends on the discharge at which flow control shifts from upstream to downstream (Qc), as well as the free drainage discharge (Qf). A theoretical design criterion for stable flow is developed in terms of Qf and Qc as a function of the vortex inlet geometry. For a “stable” design, the flow in the tapering inlet evolves from supercritical flow to subcritical flow smoothly as the discharge increases. Fifteen different tangential vortex intake models are tested. The experimental observations are in excellent agreement with the theoretical prediction. The present study provides a general guideline for designing a tangential vortex intake that can convey the flow smoothly without unstable fluctuating flow associated with a hydraulic jump.     

Mathematical Formulation and Validation of a Mixed Finite Element–Finite Difference Model for Simulating Phreatic Surfaces

The phreatic surface in an unconfined aquifer exists as a movable interface between the saturated and unsaturated zones. The movement of the phreatic surface depends on recharge, hydraulic conductivity, porosity, and horizontal and vertical flows. The location of the phreatic surface helps define the variably saturated flow domain in the subsurface. The variably saturated flow process in the subsurface is described by a parabolic partial differential equation. In this equation, the hydraulic conductivity and soil moisture capacity are used as the subsurface characteristics. The location of the phreatic surface is governed by a first-order partial differential equation. The governing parabolic partial differential equation is solved using a variational finite element formulation. The first order phreatic surface equation is then solved by loosely coupling with the governing parabolic partial differential equation describing the variably saturated flow. In the present study, a two-dimensional space is used to investigate the movement of the phreatic surface in a variably saturated unconfined flow domain. Based on the time-varying solutions of hydraulic heads, the location of the phreatic surface is simulated in a finite two-dimensional space.     

Evaluating the Long-Term Performance of Geosynthetic Clay Liners Exposed to Freeze-Thaw

An important consideration for landfill liners and covers constructed in the frost zone of cold climates is the possible deterioration in performance due to freeze-thaw cycling over the design life of the liner or cover system. Several examples in the literature show that geosynthetic clay liners can withstand a limited number of freeze-thaw events, but data on long-term freeze-thaw performance are lacking. The objective of this study was to examine the long-term performance of geosynthetic clay liners exposed to repeated freeze-thaw cycles, encompassing their application as a final cover as well as a bottom liner. Measurements of hydraulic conductivity were performed after as many as 150 freeze-thaw cycles, with no appreciable increases observed.     

Concepts for Lunar Outpost Development

This paper summarizes the results of a qualitative investigation to identify concepts for design and construction of near‐term lunar facilities. Accomplishing such construction will require an adaptation or transfer of current terrestrial technology and methods. Discussions on modularization, geosynthetic materials, aluminum materials, static load analysis, and dynamic load analysis provide illustrative examples of how terrestrial technologies can be adapted to lunar applications. These discussions provide support for the development of a phased lunar construction strategy. The initial stage of construction is characterized by small self‐supporting accomodation and laboratory modules. The assembly facility stage is characterized by the construction of a large pressurized module‐assembly facility. The module production stage is characterized by the fitting together of terrestrial or low earth‐orbit subassemblies into completed modules within the module assembly facility. The completed modules are also tested and moved to their final location in this stage. The lunar materials stage is characterized by the construction of facilities with maximum use of lunar materials.     

Undular Hydraulic Jumps in Circular Conduits

Undular hydraulic jumps in circular conduits are considered with an experimental approach. Based on previous findings in rectangular channels, this research indicates differences in terms of shape effects. All present results depend on the filling ratio of the upstream conduit flow in addition to the upstream Froude number. The results include information on the wave crests and troughs, wave lengths, and generalized axial surface profiles. The wall surface profile is shown to be similar to the axial wave profile, but with smaller wave extrema and a wave shift. The design of conduits containing undular jumps should be avoided because of unstable flow. It is also demonstrated that conduits may choke in the presence of undular jumps, with a previously established choking number relating to a design limit. For flows with choking numbers in excess of 1, choking occurs associated with a transition from the free surface to the pressurized conduit flow.     

Transit-Time Design for Diffusion through Composite Liners

Transit-time design methods are presented in this paper for determining the design thickness for composite liners consisting of a geomembrane and a compacted soil liner or geosynthetic clay liner. The design methods are based on a closed-form analytical solution for transient solute diffusion of volatile organic compounds in a composite liner and results from a numerical model. An analytical solution for diffusion in a two-layer soil profile, which is useful for transit-time design of composite liners, is also presented. The analytical solutions are used to develop graphical solution charts that can be used to design composite liners for which the effluent concentration and contaminant flux are less than a specified value. Design examples are included for a composite liner having a compacted soil liner and a composite liner having a geosynthetic clay liner. The method is relatively simple to apply and can be used for preliminary design of composite liners, evaluating experimental results, and verifying more complex numerical models.     

Design Method for Geogrid-Reinforced Unpaved Roads. I. Development of Design Method

A theoretically based design method for the thickness of the base course of unpaved roads is developed in this paper, which considers distribution of stress, strength of base course material, interlock between geosynthetic and base course material, and geosynthetic stiffness in addition to the conditions considered in earlier methods: traffic volume, wheel loads, tire pressure, subgrade strength, rut depth, and influence of the presence of a reinforcing geosynthetic (geotextile or geogrid) on the failure mode of the unpaved road or area. In this method, the required base course thickness for a reinforced unpaved road is calculated using a unique equation, whereas more than one equation was needed with earlier methods. This design method was developed for geogrid-reinforced unpaved roads. However, it can be used for geotextile-reinforced unpaved roads and for unreinforced roads with appropriate values of relevant parameters. The calibration of this design method using data from field wheel load tests and laboratory cyclic plate loading tests on unreinforced and reinforced base courses is presented in the companion paper by the authors.