Water balance modelling in a tropical watershed under deciduous forest (Mule Hole, India): Regolith matric storage buffers the groundwater recharge process

Accurate estimations of water balance are needed in semi-arid and sub-humid tropical regions, where water resources are scarce compared to water demand. Evapotranspiration plays a major role in this context, and the difficulty to quantify it precisely leads to major uncertainties in the groundwater recharge assessment, especially in forested catchments. In this paper, we propose to assess the importance of deep unsaturated regolith and water uptake by deep tree roots on the groundwater recharge process by using a lumped conceptual model (COMFORT). The model is calibrated using a 5 year hydrological monitoring of an experimental watershed under dry deciduous forest in South India (Mule Hole watershed). The model was able to simulate the stream discharge as well as the contrasted behaviour of groundwater table along the hillslope. Water balance simulated for a 32 year climatic time series displayed a large year-to-year variability, with alternance of dry and wet phases with a time period of approximately 14 years. On an average, input by the rainfall was 1090 mm year(-1) and the evapotranspiration was about 900 mm year(-1) out of which 100 mm year(-1) was uptake from the deep saprolite horizons. The stream flow was 100 mm year(-1) while the groundwater underflow was 80 mm year(-1). The simulation results suggest that (i) deciduous trees can uptake a significant amount of water from the deep regolith, (ii) this uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers and (iii) water table response to recharge is buffered due to the long vertical travel time through the deep vadose zone, which constitutes a major water reservoir. This study stresses the importance of long term observations for the understanding of hydrological processes in tropical forested ecosystems. (C) 2009 Elsevier B.V. All rights reserved.

BDST modelling of sodium ion exchange column behaviour with strong acid cation resin in relation to coal seam water treatment

Reverse osmosis is the dominant technology utilized for desalination of saline water produced during the extraction of coal seam gas. Alternatively, ion exchange is of interest due to potential cost advantages. However, there is limited information regarding the column performance of strong acid cation resin for removal of sodium ions from both model and actual coal seam water samples. In particular, the impact of bed depth, flow rate, and regeneration was not clear. Consequently, this study applied Bed Depth Service Time (BDST) models to reveal that increasing sodium ion concentration and flow rates diminished the time required for breakthrough to occur. The loading of sodium ions on fresh resin was calculated to be ca. 71.1 g Na/kg resin. Difficulties in regeneration of the resin using hydrochloric acid solutions were discovered, with 86% recovery of exchange sites observed. The maximum concentration of sodium ions in the regenerant brine was found to be 47,400 mg/L under the conditions employed. The volume of regenerant waste formed was 6.2% of the total volume of water treated. A coal seam water sample was found to load the resin with only 53.5 g Na/kg resin, which was consistent with not only the co-presence of more favoured ions such as calcium, magnesium, barium and strontium, but also inefficient regeneration of the resin prior to the coal seam water test.

Numerical Simulation of Groundwater Flow with Gradually Increasing Heterogeneity due to Clogging

Well injection replenishes depleting water levels in a well field. Observation well water levels some distance away from the injection well are the indicators of the success of a well injection program. Simulation of the observation well response, located a few tens of meters from the injection well, is likely to be affected by the effects of nonhomogeneous medium, inclined initial water table, and aquifer clogging. Existing algorithms, such as the U.S. Geological Survey groundwater flow software MODFLOW, are capable of handling the first two conditions, whereas time-dependent clogging effects are yet to be introduced in the groundwater flow models. Elsewhere, aquifer clogging is extensively researched in theory of filtration; scope for its application in a well field is a potential research problem. In the present paper, coupling of one such filtration theory to MODFLOW is introduced. Simulation of clogging effects during “Hansol” well recharge in the parts of western India is found to be encouraging.

The high speed water tunnel facility at the Indian Institute of Science

The present article about the high speed water tunnel facility at the Indian Institute of Science, Bangalore, provides a general description of the tunnel circuit, and brief reports on the performance of the facility and some typical results from investigations carried out in it. A unique aspect of the facility is that it has a horizontal resorber in the form of a large cylindrical tank located in the lower leg of the circuit. The facility has been used, among other things, for flow visualization studies, and investigations on marine propeller hydrodynamics and “synthetic cavitation”. The last topic has been primarily developed at the Indian Institute of Science and shows considerable promise for basic work in cavitation inception and noise.

Study Of Fluid-Flow And Residence-Time Distribution In A Continuous Slab Casting Mold

The fluid-flow pattern and residence-time distribution (r.t.d.) of the fluid in a continuous casting mould have been studied using a water model. The two recirculating zones below the discharge ports have been found to be asymmetric. The effect of casting speed, discharge port diameter, shroud well depth and the immersion depth on r.t.d. have been investigated. The r.t.d. curve has been well represented by a model of two backmix cells of equal volume in series. The exist of the fluid has been found to be non-uniform across the cross-section of the mould. The fluid-flow pattern has been observed to change with time in a random fashion. Dead volume of upto 31.8% has been found with smaller discharge ports.

Influence of Osmotic Suction on the Soil-Water Characteristic Curves of Compacted Expansive Clay

Unsaturated clays are subject to osmotic suction gradients in geoenvironmental engineering applications and it therefore becomes important to understand the effect of these chemical concentration gradients on soil-water characteristic curves (SWCCs). This paper brings out the influence of induced osmotic suction gradient on the wetting SWCCs of compacted clay specimens inundated with sodium chloride solutions/distilled water at vertical stress of 6.25 kPa in oedometer cells. The experimental results illustrate that variations in initial osmotic suction difference induce different magnitudes of osmotic induced consolidation and osmotic consolidation strains thereby impacting the wetting SWCCs and equilibrium water contents of identically compacted clay specimens. Osmotic suction induced by chemical concentration gradients between reservoir salt solution and soil-water can be treated as an equivalent net stress component, (p(pi)) that decreases the swelling strains of unsaturated specimens from reduction in microstructural and macrostructural swelling components. The direction of osmotic flow affects the matric SWCCs. Unsaturated specimens experiencing osmotic induced consolidation and osmotic consolidation develop lower equilibrium water content than specimens experiencing osmotic swelling during the wetting path. The findings of the study illustrate the need to incorporate the influence of osmotic suction in determination of the matric SWCCs.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Analysis of flow near a dug well in an unconfined aquifer

A numerical analysis of flow to a dug well in an unconfined aquifer is made, taking into account well storage, elastic storage release, gravity drainage, anisotropy, partial penetration, vertical flow and seepage surface at the well face, and treating the water table in the aquifer and water level in the well as unknown boundaries. The pumped discharge is maintained constant. The solution is obtained by a two-level iterative scheme. The effects of governing parameters on the drawdown, development of seepage surface and contribution from aquifer flow to the total discharge are discussed. The degree of anisotropy and partial penetration are found to be the parameters which affect the flow characteristics most significantly. The effect of anisotropy on the development of seepage surface is very pronounced.

Blade manipulators in turbulent channel flow

We report here the results of a series of careful experiments in turbulent channel flow, using various configurations of blade manipulators suggested as optimal in earlier boundary layer studies. The mass flow in the channel could be held constant to better than 0.1%, and the uncertainties in pressure loss measurements were less than 0.1 mm of water; it was therefore possible to make accurate estimates of the global effects of blade manipulation of a kind that are difficult in boundary layer flows. The flow was fully developed at the station where the blades were mounted, and always relaxed to the same state sufficiently far downstream. It is found that, for a given mass flow, the pressure drop to any station downstream is always higher in the manipulated than in the unmanipulated flow, demonstrating that none of the blade manipulators tried reduces net duct losses. However the net increase in duct losses is less than the drag of the blade even in laminar flow, showing that there is a net reduction in the total skin friction drag experienced by the duct, but this relief is only about 20% of the manipulator drag at most.

Exact free-surface flows for shallow-water equations .1. - the incompressible case

A comprehensive exact treatment of free surface flows governed by shallow water equations (in sigma variables) is given. Several new families of exact solutions of the governing PDEs are found and are shown to embed the well-known self-similar or traveling wave solutions which themselves are governed by reduced ODEs. The classes of solutions found here are explicit in contrast to those found earlier in an implicit form. The height of the free surface for each family of solutions is found explicitly. For the traveling or simple wave, the free surface is governed by a nonlinear wave equation, but is arbitrary otherwise. For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed; in another case, the free surface is a horizontal plane while the flow underneath is a sine wave. The existence of simple waves on shear flows is analytically proved. The interaction of large amplitude progressive waves with shear flow is also studied.

Unsteady nonsimilar two-dimensional and axisymmetric water boundary layers with variable viscosity and prandtl number

The influence of temperature-dependent viscosity and Prandtl number on the unsteady laminar nonsimilar forced convection flow over two-dimensional and axisymmetric bodies has been examined where the unsteadiness and (or) nonsimilarity are (is) due to the free stream velocity, mass transfer, and transverse curvature. The partial differential equations governing the flow which involve three independent variables have been solved numerically using an implicit finite-difference scheme along with a quasilinearization technique. It is found that both the skin friction and heat transfer strongly respond to the unsteady free stream velocity distributions. The unsteadiness and injection cause the location of zero skin friction to move upstream. However, the effect of variable viscosity and Prandtl number is to move it downstream. The heat transfer is found to depend strongly on viscous dissipation, but the skin friction is little affected by it. In general, the results pertaining to variable fluid properties differ significantly, from those of constant fluid properties.

Water-modeling study of the surface disturbances in continuous slab caster

The present work is based on four static molds using nozzles of different port diameter, port angle, and immersion depth. It has been observed that the meniscus is wavy. The wave amplitude shows a parabolic variation with the nozzle exit velocity. The dimensionless amplitude is found to vary linearly with the Froude number. Vortex formation and bubble entrainment by the wave occurs at the meniscus beyond a critical flow rate, depending upon the nozzle configuration, immersion depth, and the mold aspect ratio.

Water Modeling Study Of The Jet Characteristics In A Continuous-Casting Mold

The jet characteristics and the fluid flow pattern in a continuous slab caster have been studied using a water model. The fluid jet is studied under free fall and submerged discharge conditions. In the latter case, the jet was followed by dye-injection technique and image analyser was used to find out the effect of nozzle parameters on jet-spread angle, jet-discharge angle and the volume entrainment by the jet. All free-fall jets with nozzle port angle zero and upward are found to be spinning. Some of the free-fall jets with downward nozzle-port angle are found to be spinning and rest are smooth. The spinning direction of the jets are found to change with time. The well depth, port diameter and the inner diameter of the nozzle have a clear effect on the free-fall jets with downward port angle. The jet-spread angle is found to be about 17-degrees for smooth jets. The spread angle for spinning jet increases as the nozzle-port angle is increased from downward 25 to upward 15-degrees. The jet-discharge angle is always downward even when the nozzle-discharge ports are angled upward. The extent of volume entrainment by the spinning jet is higher and it increases as the nozzle-port angle is increased from 25 downward to 15-degrees upward.

Levitation of a drop over a film flow

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 mi) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.