Optimal scheduling of multiple chlorine sources in water distribution systems

The specified range of free chlorine residual (between minimum and maximum) in water distribution systems needs to be maintained to avoid deterioration of the microbial quality of water, control taste and/or odor problems, and hinder formation of carcino-genic disinfection by-products. Multiple water quality sources for providing chlorine input are needed to maintain the chlorine residuals within a specified range throughout the distribution system. The determination of source dosage (i.e., chlorine concentrations/chlorine mass rates) at water quality sources to satisfy the above objective under dynamic conditions is a complex process. A nonlinear optimization problem is formulated to determine the chlorine dosage at the water quality sources subjected to minimum and maximum constraints on chlorine concentrations at all monitoring nodes. A genetic algorithm (GA) approach in which decision variables (chlorine dosage) are coded as binary strings is used to solve this highly nonlinear optimization problem, with nonlinearities arising due to set-point sources and non-first-order reactions. Application of the model is illustrated using three sample water distribution systems, and it indicates that the GA,is a useful tool for evaluating optimal water quality source chlorine schedules.

Water quality parameter estimation in steady-state distribution system

The maintenance of chlorine residual is needed at all the points in the distribution system supplied with chlorine as a disinfectant. The propagation and level of chlorine in a distribution system is affected by both bulk and pipe wall reactions. It is well known that the field determination of wall reaction parameter is difficult. The source strength of chlorine to maintain a specified chlorine residual at a target node is also an important parameter. The inverse model presented in the paper determines these water quality parameters, which are associated with different reaction kinetics, either in single or in groups of pipes. The weighted-least-squares method based on the Gauss-Newton minimization technique is used for the estimation of these parameters. The validation and application of the inverse model is illustrated with an example pipe distribution system under steady state. A generalized procedure to handle noisy and bad (abnormal) data is suggested, which can be used to estimate these parameters more accurately. The developed inverse model is useful for water supply agencies to calibrate their water distribution system and to improve their operational strategies to maintain water quality.

Estimation of groundwater recharge of shallow aquifer on humid environment in Yaounde, Cameroon using hybrid water-fluctuation and hydrochemistry methods

A study of environmental chloride and groundwater balance has been carried out in order to estimate their relative value for measuring average groundwater recharge under a humid climatic environment with a relatively shallow water table. The hybrid water fluctuation method allowed the split of the hydrologic year into two seasons of recharge (wet season) and no recharge (dry season) to appraise specific yield during the dry season and, second, to estimate recharge from the water table rise during the wet season. This well elaborated and suitable method has then been used as a standard to assess the effectiveness of the chloride method under forest humid climatic environment. Effective specific yield of 0.08 was obtained for the study area. It reflects an effective basin-wide process and is insensitive to local heterogeneities in the aquifer system. The hybrid water fluctuation method gives an average recharge value of 87.14 mm/year at the basin scale, which represents 5.7% of the annual rainfall. Recharge value estimated based on the chloride method varies between 16.24 and 236.95 mm/year with an average value of 108.45 mm/year. It represents 7% of the mean annual precipitation. The discrepancy observed between recharge value estimated by the hybrid water fluctuation and the chloride mass balance methods appears to be very important, which could imply the ineffectiveness of the chloride mass balance method for this present humid environment.

Theory of Gas Hydrates: Effect of the Approximation of Rigid Water Lattice

One of the assumptions of the van der Waals and Platteeuw theory for gas hydrates is that the host water lattice is rigid and not distorted by the presence of guest molecules. In this work, we study the effect of this approximation on the triple-point lines of the gas hydrates. We calculate the triple-point lines of methane and ethane hydrates via Monte Carlo molecular simulations and compare the simulation results with the predictions of van der Waals and Platteeuw theory. Our study shows that even if the exact intermolecular potential between the guest molecules and water is known, the dissociation temperatures predicted by the theory are significantly higher. This has serious implications to the modeling of gas hydrate thermodynamics, and in spite of the several impressive efforts made toward obtaining an accurate description of intermolecular interactions in gas hydrates, the theory will suffer from the problem of robustness if the issue of movement of water molecules is not adequately addressed.

String-like propagation of the 5-coordinated defect state in supercooled water: molecular origin of dynamic and thermodynamic anomalies

We find that at low temperature water, large amplitude (similar to 60 degrees) rotational jumps propagate like a string, with the length of propagation increasing with lowering temperature. The strings are formed by mobile 5-coordinated water molecules which move like a Glarum defect (J. Chem. Phys., 1960, 33, 1371), causing water molecules on the path to change from 4-coordinated to 5-coordinated and again back to 4-coordinated water, and in the process cause the tagged water molecule to jump, by following essentially the Laage-Hynes mechanism (Science, 2006, 311, 832-835). The effects on relaxation of the propagating defect causing large amplitude jumps are manifested most dramatically in the mean square displacement (MSD) and also in the rotational time correlation function of the O-H bond of the molecule that is visited by the defect (transient transition to the 5-coordinated state). The MSD and the decay of rotational time correlation function, both remain quenched in the absence of any visit by the defect, as postulated by Glarum long time ago. We establish a direct connection between these propagating events and the known thermodynamic and dynamic anomalies in supercooled water. These strings are found largely in the regions that surround the relatively rigid domains of 4-coordinated water molecules. The propagating strings give rise to a noticeable dynamical heterogeneity, quantified here by a sharp rise in the peak of the four-point density response function, chi(4)(t). This dynamics heterogeneity is also responsible for the breakdown of the Stokes-Einstein relation.

Modeling Indian Ocean circulation: Bay of Bengal fresh water plume and the Arabian Sea Mini Warm Pool

The Indian subcontinent divides the north Indian Ocean into two tropical basins, namely the Arabian Sea and the Bay of Bengal. The Arabian Sea has high salinity whereas the salinity of the Bay of Bengal is much lower due to the contrast in freshwater forcing of the two basins. The freshwater received by the Bay in large amounts during the summer monsoon through river discharge is flushed out annually by ocean circulation. After the withdrawal of the summer monsoon, the Ganga – Brahmaputra river plume flows first along the Indian coast and then around Sri Lanka into the Arabian Sea creating a low salinity pool in the southeastern Arabian Sea (SEAS). In the same region, during the pre-monsoon months of February – April, a warm pool, known as the Arabian Sea Mini Warm Pool (ASMWP), which is distinctly warmer than the rest of the Indian Ocean, takes shape. In fact, this is the warmest region in the world oceans during this period. Simulation of the river plume and its movement as well as its implications to thermodynamics has been a challenging problem for models of Indian Ocean. Here we address these issues using an ocean general circulation model – first we show that the model is capable of reproducing fresh plumes in the Bay of Bengal as well as its movement and then we use the model to determine the processes that lead to formation of the ASMWP. Hydrographic observations from the western Bay of Bengal have shown the presence of a fresh plume along the northern part of the Indian coast during summer monsoon. The Indian Ocean model when forced by realistic winds and climatological river discharge reproduces the fresh plume with reasonable accuracy. The fresh plume does not advect along the Indian coast until the end of summer monsoon. The North Bay Monsoon Current, which flows eastward in the northern Bay, separates the low salinity water from the more saline southern parts of the bay and thus plays an important role in the fresh water budget of the Bay of Bengal. The model also reproduces the surge of the fresh-plume along the Indian coast, into the Arabian Sea during northeast monsoon. Mechanisms that lead to the formation of the Arabian Sea Mini Warm Pool are investigated using several numerical experiments. Contrary to the existing theories, we find that salinity effects are not necessary for the formation of the ASMWP. The orographic effects of the Sahyadris (Western Ghats) and resulting reduction in wind speed leads to the formation of the ASMWP. During November – April, the SEAS behave as a low-wind heatdominated regime where the evolution of sea surface temperature is solely determined by atmospheric forcing. In such regions the evolution of surface layer temperature is not dependent on the characteristics of the subsurface ocean such as the barrier layer and temperature inversion.

The Behaviour of Two-Phase Flow of DNAPL and Water through a Fractured Rock under Confining Pressure

This study presents the characterization of DNAPL and water flow in a fracture under confining pressure. A comprehensive mathematical model and the conditions under which DNAPL will enter an initially water-saturated deforming rock fracture are discussed. A numerical model with which to predict the quantity of each phase in terms of their saturations in deforming rock joint is developed. The effect of varying confining stresses on the traverse time of DNAPL across a fractured aquitard is studied. The sensitivity analysis for physical and hydraulic properties like initial fracture apertures, fracture dips, equivalent fracture aperture and confining pressures are performed and discussed.

Exploring use of climate information for optimal design of water resources infrastructure

We conducted surveys of fire and fuels managers at local, regional, and national levels to gain insights into decision processes and information flows in wildfire management. Survey results in the form of fire managers’ decision calendars show how climate information needs vary seasonally, over space, and through the organizational network, and help determine optimal points for introducing climate information and forecasts into decision processes. We identified opportunities to use climate information in fire management, including seasonal to interannual climate forecasts at all organizational levels, to improve the targeting of fuels treatments and prescribed burns, the positioning and movement of initial attack resources, and staffing and budgeting decisions. Longer-term (5–10 years) outlooks also could be useful at the national level in setting budget and research priorities. We discuss these opportunities and examine the kinds of organizational changes that could facilitate effective use of existing climate information and climate forecast capabilities.

Spin Probe Signature of Freezing in Water Slow Cooling Experiment: Our Observation is Global or Local?

First systematic spin probe ESR study of water freezing has been conducted using TEMPOL and TEMPO as the probes. The spin probe signature of the water freezing has been described in terms of the collapse of narrow triplet spectrum into a single broad line. This spin probe signature of freezing has been observed at an anomalously low temperature when a milimoler solution of TEMPOL is slowly cooled from room temperature. A systematic observation has revealed a spin probe concentration dependence of these freezing and respective melting points. These results can be explained in terms of localization of spin probe and liquid water,most probably in the interstices of ice grains, in an ice matrix. The lowering of spin probe freezing point, along with the secondary evidences, like spin probe concentration dependence of peak-to-peak width in frozen limit signal, indicates a possible size dependence of these localizations/entrapments with spin probe concentration. A weak concentration dependence of spin probe assisted freezing and melting points, which has been observed for TEMPO in comparison to TEMPOL, indicates different natures of interactions with water of these two probes. This view is also supported by the relaxation behavior of the two probes.