A numerical computation of a non-dimensional form of stream water quality model with hydrodynamic advection–dispersion–reaction equations

Mathematical models of water quality assessment problems often arise in environmental science. The modelling often involves numerical methods to solve the equations. In this research, two mathematical models are used to simulate pollution due to sewage effluent in the nonuniform flow of water in a stream with varied current velocity. The first is a hydrodynamic model that provides the velocity field and elevation of the water flow. The second is a dispersion model, where the commonly used governing factor is the one-dimensional advection–dispersion–reaction equation that gives the pollutant concentration fields. In the simulation processes, we used the Crank–Nicolson method system of a hydrodynamic model and the backward time central space scheme for the dispersion model. Finally, we present a numerical simulation that confirms the results of the techniques.     

Compartmental model for nitrogen dynamics in citrus orchards

A simple model that represents the soil nitrogen dynamics in a citrus orchard is studied. This model consists of several compartments of organic and inorganic nitrogen representing the main processes that occur among the compartment of the soil column. Some of these processes are coupled to the carbon dynamics in the soil, thus, the evolution of organic carbon in the soil has been also described in the model. The dependence of these processes with soil moisture requires the coupling of a nitrogen model with a model of soil water dynamics. A compartmental model for the water has been used to simulate the dynamics of water in the root profile.The proposed model has been used to predict the soil mineral nitrogen content and the nitrate leaching in a citrus orchard placed in the area of Valencia and the obtained results have been compared with the results obtained with the nitrogen component of the widely used Leaching Estimation and Chemistry Model (LEACHM).     

The economics of water resources for the generation of electricity and other uses

An optimal control approach is used to analyze the tradeoff between the use of water resources for electricity generation versus other economic uses (irrigation, industry, etc.). For that purpose, a dynamic model is presented which establishes relationships between economic growth, water resources management, and energy policy in the context of the aforementioned tradeoff, in an economy whose energy matrix is heavily dependent upon hydroelectric power. Among other results, the analysis establishes that in the market, the price of water for non-energy uses should be twice the price of the energy goods, indicating the necessity of substituting other sources of energy for hydroelectric power.     

Nash bargaining and leader–follower models in water allocation: Application to the Zarrinehrud River basin,Iran

Mathematical models for conflict resolution are very important in integrated water resources and environmental management. This study proposes a new methodology to resolve conflicts among different water users and water suppliers while considering environmental requirements and the system’s constraints. A two-level leader–follower model is applied to maximize the net benefit with the Iran Water Resources Management Company as the leader and agricultural, domestic, and industrial users as followers subject to the system’s constraints. As a comparison, the Nash bargaining solution is also used to find a solution when simultaneous moves are assumed by the participants. The suggested method is then applied to the real case of the Zarrinehrud River basin that is one of the areas facing water shortages in Iran. For the actual optimization, Genetic Algorithm is used in order to avoid local optimum. As the contribution of this study, the results show that benefits for the leader in the leader–follower model increased in comparison with the Nash bargaining solutions.     

A kinetic formulation for a model coupling free surface and pressurised flows in closed pipes

The aim of this paper is to present a kinetic formulation of a model for the coupling of transient free surface and pressurised flows. Firstly, we revisit the system of Saint-Venant equations for free surface flow: we state some properties of Saint-Venant equations, we propose a kinetic formulation and we verify that this kinetic formulation leads to a Gibbs equilibrium that minimises (in some general case) an energy and preserves the still water steady state. Secondly, we propose a model for pressurised flows in a Saint-Venant-like conservative formulation. We then propose a kinetic formulation and we verify that this kinetic formulation leads to a Gibbs equilibrium that minimises in any case an energy and preserves the still water steady state. Finally, we propose a dual model that couples these two types of flow.     

Evaluating risk of water mains failure using a Bayesian belief network model

It has been reported that since year 2000, there have been an average 700 water main breaks per day only in Canada and the USA costing more than CAD 10 billions/year. Moreover, water main leaks affect other neighboring infrastructure that may lead to catastrophic failures. For this, municipality authorities or stakeholders are more concerned about preventive actions rather reacting to failure events. This paper presents a Bayesian Belief Network (BBN) model to evaluate the risk of failure of metallic water mains using structural integrity, hydraulic capacity, water quality, and consequence factors. BBN is a probabilistic graphical model that represents a set of variables and their probabilistic relationships, which also captures historical information about these dependencies. The proposed model is capable of ranking water mains within distribution network that can identify vulnerable and sensitive pipes to justify proper decision action for maintenance/rehabilitation/replacement (M/R/R). To demonstrate the application of proposed model, water distribution network of City of Kelowna has been studied. Result indicates that almost 9% of the total 259 metallic pipes are at high risk in both summer and winter.     

An equity approach to stream water quality management

A model for stream water quality management with particular emphasis on applicability in Turkey is developed via mathematical programming. A new criterion of equity is introduced and the model is formulated as a quadratic programming problem with two objectives, namely equity among dischargers and minimization of total cost. The model is also applied to a hypothetical stream and compared with other more common water quality management programs.     

Existence and qualitative theory for stratified solitary water waves

This paper considers two-dimensional gravity solitary waves moving through a body of density stratified water lying below vacuum. The fluid domain is assumed to lie above an impenetrable flat ocean bed, while the interface between the water and vacuum is a free boundary where the pressure is constant. We prove that, for any smooth choice of upstream velocity field and density function, there exists a continuous curve of such solutions that includes large-amplitude surface waves. Furthermore, following this solution curve, one encounters waves that come arbitrarily close to possessing points of horizontal stagnation.We also provide a number of results characterizing the qualitative features of solitary stratified waves. In part, these include bounds on the wave speed from above and below, some of which are new even for constant density flow; an a priori bound on the velocity field and lower bound on the pressure; a proof of the nonexistence of monotone bores in this physical regime; and a theorem ensuring that all supercritical solitary waves of elevation have an axis of even symmetry.     

Robust placement of sensors in dynamic water distribution systems

Designing a robust sensor network to detect accidental contaminants in water distribution systems is a challenge given the uncertain nature of the contamination events (what, how much, when, where and for how long) and the dynamic nature of water distribution systems (driven by the random consumption of consumers). We formulate a set of scenario-based minimax and minimax regret models in order to provide robust sensor-placement schemes that perform well under all realizable contamination scenarios, and thus protect water consumers. Single-and multi-objective versions of these models are then applied to a real water distribution system. A heuristic solution method is applied to solve the robust models. The concept of “sensitivity region” is used to visualize trade-offs between multiple objectives.     

Study of storm surge due to Typhoon Linda (1997) in the Gulf of Thailand using a three dimensional ocean model

Numerical integrations using the three dimensional ocean model based on the princeton ocean model (POM) were applied for the study of both sea level elevation and ocean circulation patterns forced by the wind fields during typhoons that moved over the Gulf of Thailand (GoT). The simulation concerned a case of Typhoon Linda which occurred during November 1-4, 1997. Typhoon Linda was one of the worst storms that passed the Gulf of Thailand and hit the southern coastal provinces of Thailand on November 3, 1997. It caused flooding and a strong wind covering large areas of agriculture and fisheries, which destroyed households, utilities and even human lives. The model is the time-dependent, primitive equation, Cartesian coordinates in a horizontal and sigma coordinate in the vertical. The model grid has 37 × 97 orthogonal curvilinear grid points in the horizontal, with variable spacing from 2 km near the head of the GoT to 55 km at the eastern boundary, with 10 sigma levels in the vertical conforming to a realistic bottom topography. Open boundary conditions are determined by using radiation conditions, and the sea surface elevation is prescribed from the archiving, validation and interpretation of satellite oceanographic data (AVISO). The initial condition is determined from the spin up phase of the first model run, which was executed by using wind stress calculated from climatological monthly mean wind, restoring-type surface heat and salt and climatological monthly mean freshwater flux. The model was run in spin up phase until an ocean model reached an equilibrium state under the applied force. A spatially variable wind field taken from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used to compute the wind stress directly from the velocity fluctuations. Comparison of tendency between the sea surface elevations from model and the observed significant wave heights of moored buoys in the Gulf of Thailand under Seawatch project is investigated. The model predicts the sea level elevation up to 68.5 cm at the Cha-Am area located in the north of where the typhoon strands to the shore. Results of sea level elevation show that there is an area of peak set-up in the upper gulf, particularly in the western coast, and the effects of the storm surge are small at the lower gulf. During the entire period of this study, the surge in the gulf was induced by the northeasterly wind blowing over it.     

Finite element modelling of surface waves in a channel

A variational formulation of the vertically-integrated differential equations for free surface wave motion is presented. A finite element model is derived for solving this nonlinear system of hydrodynamic equations. The time integration scheme employed is discussed and the results obtained demonstrate its good stability and accuracy.Several applications of the model are considered: the first problem is an open channel of uniform depth and the second an open channel of linearly varying depth. The ‘inflow’ boundary condition is prescribed in terms of the velocity which represents a wavemaker and/or a flow source, while the ‘outflow’ boundary condition is specified in terms of the water elevation. The outflow condition is adjusted for two cases, a reflecting boundary (finite channel) and a non-reflecting boundary (open-ended channel). The latter boundary condition is examined in some detail and the results obtained show that the numerical model can produce the non-reflecting boundary that is similar to the analytical radiation condition for waves. Computational results for a third problem, involving wave reflection from a submerged cylinder, are also presented and compared with both experimental data and analytical predictions.The simplicity and the performance of the computational model suggest that free surface waves can be simulated without excessively complicated numerical schemes. The ability of the model to simulate outflow boundary conditions properly is of special importance since these conditions present serious problems for many numerical algorithms.     

Numerical modeling of hydrodynamic and hydrothermal characteristics in subtropical alpine lake

A three-dimensional, time-dependent hydrodynamic and hydrothermal model was performed and applied to the subtropical alpine Yuan-Yang Lake (YYL) in northeastern region of Taiwan. The model was driven with discharge inflow, heat, and wind stress to simulate the hydrodynamic and hydrothermal in the lake. The model was validated with measured water surface elevation, current, and temperature in 2008. The overall model simulation results are in quantitative agreement with the available field data. The validated model was then used to investigate wind-driven current, mean circulation, and residence time in the YYL. The modeling results reveal that the velocity field along the wind axis present the variations over depth with return current where the velocity at the surface layer is along the wind direction while it is opposite near 1 m below water surface. The simulated mean current indicates that the surface currents flow towards the southwest direction and form a clock-wise rotation. The calculated residence time is strongly dependent on the inflows and wind effects. Regression analysis of model results reveals that an exponential regression equation can be employed to correlate the residence time to change of discharge input. The residence time without wind stress is higher than that with wind effect, indicating that wind plays an important role in lake mixing. The calculated residence time is approximately 2-2.5 days under low inflow with wind effect.     

A 3D hydrodynamic numerical model of the Río de la Plata and Montevideo’s coastal zone

The 3D hydrodynamic numerical model MOHID was applied in the Río de la Plata and Montevideo coastal zone in order to represent the main dynamics and to study its complex circulation pattern. The hydrodynamic model was calibrated and validated considering the following main forces: fresh water flow, astronomical and meteorological tides in the oceanic boundary, and wind acting on the water surface. A series of water levels measured at six coastal stations and vertical profiles of current velocity measured at four different locations in the estuarine zone of the Río de la Plata were used for calibrating and validating the hydrodynamic model. The calibration process was carried out in two steps. First the astronomical waves propagation was calibrated comparing harmonic constants of observed and computed sea surface elevation data. Next, both the astronomical and meteorological wave propagation was calibrated. Direct comparison of scatter plot and root-mean square errors of model results and field data were used when evaluating the calibration quality. The calibrated model shows good agreement with the measured water surface level in the entire domain with mean error values being minor than 20% of the measured data and correlation factors higher than 0.74. Also, the intensity and velocity direction observed in the currents data are well represented by the model in both bottom and surface levels with errors similar to 30% of the currents data components. Using the 3D calibrated model the bottom and surface residual circulation for a four month period of time was analyzed.     

Water waves diffraction by a circular plate

The interaction of water waves with circular plate within the framework of a linear theory is considered. The plate lies on the free surface in water of finite depth. The integral transform technique is used to solve this problem. The problem is reduced to a system of dual integral equations for a spectral function. The way to solve these equations consists in converting them into Fredholm integral equation of the second kind. The asymptotic solutions of this equation are obtained. Representations for diffraction field and for the forces on the plate are given.     

深圳市供水量的最优组合预测

城市供水系统是一个复杂的大系统,供水量受多种因素的共同影响。本文以深圳市最近20多年的供水量历史数据为基础,建立了深圳市供水量的最优组合预测模型。该模型具有较高的预测精度,组合预测的预测效果优于任意一种单一预测的预测效果,供水量预测结果对深圳市未来供水的短期或长期规划能起到重要的宏观指导作用。     

Superharmonic instability of stokes waves

A stability of nearly limiting Stokes waves to superharmonic perturbations is considered numerically in approximation of an infinite depth. Investigation of the stability properties can give one an insight into the evolution of the Stokes wave. The new, previously inaccessible branches of superharmonic instability were investigated. Our numerical simulations suggest that eigenvalues of linearized dynamical equations, corresponding to the unstable modes, appear as a result of a collision of a pair of purely imaginary eigenvalues at the origin, and a subsequent appearance of a pair of purely real eigenvalues: a positive and a negative one that are symmetric with respect to zero. Complex conjugate pairs of purely imaginary eigenvalues correspond to stable modes, and as the steepness of the underlying Stokes wave grows, the pairs move toward the origin along the imaginary axis. Moreover, when studying the eigenvalues of linearized dynamical equations we find that as the steepness of the Stokes wave grows, the real eigenvalues follow a universal scaling law, that can be approximated by a power law. The asymptotic power law behavior of this dependence for instability of Stokes waves close to the limiting one is proposed. Surface elevation profiles for several unstable eigenmodes are made available through  http://stokeswave.org website.     

A multi-level Taguchi-factorial two-stage stochastic programming approach for characterization of parameter uncertainties and their interactions: An application to water resources management

This paper presents a multi-level Taguchi-factorial two-stage stochastic programming (MTTSP) approach for supporting water resources management under parameter uncertainties and their interactions. MTTSP is capable of performing uncertainty analysis, policy analysis, factor screening, and interaction detection in a comprehensive and systematic way. A water resources management problem is used to demonstrate the applicability of the proposed approach. The results indicate that interval solutions can be generated for the objective function and decision variables, and a variety of decision alternatives can be obtained under different policy scenarios. The experimental data obtained from the Taguchi’s orthogonal array design are helpful in identifying the significant factors affecting the total net benefit. Then the findings from the multi-level factorial experiment reveal the latent interactions among those important factors and their curvature effects on the model response. Such a sequential strategy of experimental designs is useful in analyzing the interactions for a large number of factors in a computationally efficient manner.     

Dynamic model of methane–water diffusion

This paper presents a mathematical model to describe the time evolution of the diffusion process of methane exchange between liquid and gaseous phases. In order to reach the equilibrium, the distribution of gas in the liquid phase decreases to a constant value over time. Generalized model is analyzed under the assumption that both, reversible and irreversible form of absorption and desorption, occur at the same time. For the application in case of the real system of methane–water, model is developed under the assumption that the processes of absorption and desorption are irreversible. From the experimental data are determined the coefficients of model and their dependence on initial conditions, with constant end conditions. This paper tests the introduced theoretical model on the existing experimental data.     

Multi-decisions rating model: Establishing rescue policies for Regional Drinking Water Companies (PDAMs) in Indonesia

Fresh water is abundant in Indonesia, but only about 20% out of more than 200 million of its population have access to clean water supplied by 306 Regional Drinking Water Companies (PDAMs) around the country. The poor performance of most PDAMs is due to high level of debt, lack of investment, and inefficient operations. On one hand, the PDAMs need to increase its coverage capacity but on the other hand they have been operating under unfavorable regulated tariff policies. A government team involving the Indonesia Water Enterprise Association (PERPAMSI) has developed a multi-criteria model to rate the PDAMs, grouped them based on their scores, and formulated common policies for each group. Unfortunately, the recommended policies do not appear to be specific enough to make them useful. This paper proposes improvements to the original model which used the Analytic Hierarchy Process (AHP) while maintaining the original priority judgments. Model-1 is the AHP representation of the original model, Model-2 is similar to Model-1 with the principle of hierarchic composition satisfied, and Model-3 is a framework for multi-decisions decision making involving clustering the PDAMs based on their business performance. Model-2 did not prove useful for identifying specific group policies that fit the need of every PDAM in a given group. Model-3 enables one to identify a specific policy that is applicable to either all the PDAMs or those grouped as a certain cluster by rating them based on a set of criteria that is relevant to that policy. The results of the three models are evaluated and some examples of using the Model-3 framework to rate the PDAMs, each with a specific policy in mind, are provided.