Planning and operation of large-scale water distribution systems with preemptive priorities

This paper describes the development of an optimization model for the management and operation of a large-scale, multireservoir water supply distribution system with preemptive priorities. The model considers multiobjectives and hedging rules. During periods of drought, when water supply is insufficient to meet the planned demand, appropriate rationing factors are applied to reduce water supply. In this paper, a water distribution system is formulated as a network and solved by the GAMS modeling system for mathematical programming and optimization. A user-friendly interface is developed to facilitate the manipulation of data and to generate graphs and tables for decision makers. The optimization model and its interface form a decision support system (DSS), which can be used to configure a water distribution system to facilitate capacity expansion and reliability studies. Several examples are presented to demonstrate the utility and versatility of the developed DSS under different supply and demand scenarios, including applications to one of the largest water supply systems in the world, the Sao Paulo Metropolitan Area Water Supply Distribution System in Brazil.

Residual stresses in TiN, DLC and MoS(2) coated surfaces with regard to their tribological fracture behaviour

Thin hard coatings on components and tools are used increasingly due to the rapid development in deposition techniques, tribological performance and application skills. The residual stresses in a coated surface are crucial for its tribological performance. Compressive residual stresses in PVD deposited TiN and DLC coatings were measured to be in the range of 0.03-4 GPa on steel substrate and 0.1-1.3 GPa on silicon. MoS(2) coatings had tensional stresses in the range of 0.8-1.3 on steel and 0.16 GPa compressive stresses on silicon. The fracture pattern of coatings deposited on steel substrate were analysed both in bend testing and scratch testing. A micro-scale finite element method (FEM) modelling and stress simulation of a 2 mu m TiN-coated steel surface was carried out and showed a reduction of the generated tensile buckling stresses in front of the sliding tip when compressive residual stresses of 1 GPa were included in the model. However, this reduction is not similarly observed in the scratch groove behind the tip, possibly due to sliding contact-induced stress relaxation. Scratch and bending tests allowed calculation of the fracture toughness of the three coated surfaces, based on both empirical crack pattern observations and FEM stress calculation, which resulted in highest values for TiN coating followed by MoS(2) and DLC coatings, being K(C) = 4-11, about 2, and 1-2 MPa M(1/2), respectively. Higher compressive residual stresses in the coating and higher elastic modulus of the coating correlated to increased fracture toughness of the coated surface. (C) 2009 Elsevier B.V. All rights reserved.

A comparison of the critical impeller speed for solids suspension in a bench-scale and a pilot-scale mechanical flotation cell

This paper compares the critical impeller speed results for 6 L Denver and Wemco bench-scale flotation cells with findings from a study by Van der Westhuizen and Deglon [Van der Westhuizen, A.P., Deglon, D.A., 2007. Evaluation of solids suspension in a pilot-scale mechanical flotation cell: the critical impeller speed. Minerals Engineering 20,233-240; Van der Westhuizen, A.P., Deglon, D.A., 2008. Solids suspension in a pilot scale mechanical flotation cell: a critical impeller speed correlation. Minerals Engineering 21, 621-629] conducted in a 125 L Batequip flotation cell. Understanding solids suspension has become increasingly important due to dramatic increases in flotation cell sizes. The critical impeller speed is commonly used to indicate the effectiveness of solids suspension. The minerals used in this study were apatite, quartz and hematite. The critical impeller speed was found to be strongly dependent on particle size, solids density and air flow rate, with solids concentration having a lesser influence. Liquid viscosity was found to have a negligible effect. The general Zwietering-type critical impeller speed correlation developed by Van der Westhuizen and Deglon [Van der Westhuizen, A.P., Deglon, D.A., 2008. Solids suspension in a pilot scale mechanical flotation cell: a critical impeller speed correlation. Minerals Engineering 21, 621-629] was found to be applicable to all three flotation machines. The exponents for particle size, solids concentration and liquid viscosity were equivalent for all three cells. The exponent for solids density was found to be less significant than that obtained by the previous authors, and to be consistent with values reported in the general literature for stirred tanks. Finally, a new dimensionless critical impeller speed correlation is proposed where the particle size is divided by the impeller diameter. This modified equation generally predicts the experimental measurements well, with most predictions within 10% of the experimental. (C) 2009 Elsevier Ltd. All rights reserved.

Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

MODELLING AND CONTROL OF ERSW PROCESSES BY NEURONAL NETWORK

A study on the use of artificial intelligence (AI) techniques for the modelling and subsequent control of an electric resistance spot welding process (ERSW) is presented. The ERSW process is characterized by the coupling of thermal, electrical, mechanical, and metallurgical phenomena. For this reason, early attempts to model it using computational methods established as the methods of finite differences, finite element, and finite volumes, ask for simplifications that lead the model obtained far from reality or very costly in terms of computational costs, to be used in a real-time control system. In this sense, the authors have developed an ERSW controller that uses fuzzy logic to adjust the energy transferred to the weld nugget. The proposed control strategies differ in the speed with which it reaches convergence. Moreover, their application for a quality control of spot weld through artificial neural networks (ANN) is discussed.

Modelling grain boundary migration during geometric dynamic recrystallization

High-angle grain boundary migration is predicted during geometric dynamic recrystallization (GDRX) by two types of mathematical models. Both models consider the driving pressure due to curvature and a sinusoidal driving pressure owing to subgrain walls connected to the grain boundary. One model is based on the finite difference solution of a kinetic equation, and the other, on a numerical technique in which the boundary is subdivided into linear segments. The models show that an initially flat boundary becomes serrated, with the peak and valley migrating into both adjacent grains, as observed during GDRX. When the sinusoidal driving pressure amplitude is smaller than 2 pi, the boundary stops migrating, reaching an equilibrium shape. Otherwise, when the amplitude is larger than 2 pi, equilibrium is never reached and the boundary migrates indefinitely, which would cause the protrusions of two serrated parallel boundaries to impinge on each other, creating smaller equiaxed grains.

Model for facilities or vendors location in a global scale considering several echelons in the Chain

The facilities location problem for companies with global operations is very complex and not well explored in the literature. This work proposes a MILP model that solves the problem through minimization of the total logistic cost. Main contributions of the model are the pioneer carrying cost calculation, the treatment given to the take-or-pay costs and to the international tax benefits such as drawback and added value taxes in Brazil. The model was successfully applied to a real case of a chemical industry with industrial plants and sales all over the world. The model application recommended a totally new sourcing model for the company.

Thermodynamic modelling of an Al(2)O(3)-MnO system using the ionic model

The thermodynamic assessment of an Al(2)O(3)-MnO pseudo-binary system has been carried out with the use of an ionic model. The use of the electro-neutrality principles in addition to the constitutive relations, between site fractions of the species on each sub-lattice, the thermodynamics descriptions of each solid phase has been determined to make possible the solubility description. Based on the thermodynamics descriptions of each phase in addition to thermo-chemical data obtained from the literature, the Gibbs energy functions were optimized for each phase of the Al(2)O(3)-MnO system with the support of PARROT(R) module from ThemoCalc(R) package. A thermodynamic database was obtained, in agreement with the thermo-chemical data extracted from the literature, to describe the Al(2)O(3)-MnO system including the solubility description of solid phases. (C) 2009 Elsevier Ltd. All rights reserved.

Emerging Threats and Opportunities for Large-Scale Ecological Restoration in the Atlantic Forest of Brazil

Over the past 150 years, Brazil has played a pioneering role in developing environmental policies and pursuing forest conservation and ecological restoration of degraded ecosystems. In particular, the Brazilian Forest Act, first drafted in 1934, has been fundamental in reducing deforestation and engaging private land owners in forest restoration initiatives. At the time of writing (December 2010), however, a proposal for major revision of the Brazilian Forest Act is under intense debate in the National Assembly, and we are deeply concerned about the outcome. On the basis of the analysis of detailed vegetation and hydrographic maps, we estimate that the proposed changes may reduce the total amount of potential areas for restoration in the Atlantic Forest by approximately 6 million hectares. As a radically different policy model, we present the Atlantic Forest Restoration Pact (AFRP), which is a group of more than 160 members that represents one of the most important and ambitious ecological restoration programs in the world. The AFRP aims to restore 15 million hectares of degraded lands in the Brazilian Atlantic Forest biome by 2050 and increase the current forest cover of the biome from 17% to at least 30%. We argue that not only should Brazilian lawmakers refrain from revising the existing Forest Law, but also greatly step up investments in the science, business, and practice of ecological restoration throughout the country, including the Atlantic Forest. The AFRP provides a template that could be adapted to other forest biomes in Brazil and to other megadiversity countries around the world.

Watershed scale temporal stability of soil water content

The recognition of temporally stable locations with respect to soil water content is of importance for soil water management decisions, especially in sloping land of watersheds. Neutron probe soil water content (0 to 0.8 m), evaluated at 20 dates during a year in the Loess Plateau of China, in a 20 ha watershed dominated by Ust-Sandiic Entisols and Aeolian sandy soils, were used to define their temporal stability through two indices: the standard deviation of relative difference (SDRD) and the mean absolute bias error (MABE). Specific concerns were (a) the relationship of temporal stability with soil depth, (b) the effects of soil texture and land use on temporal stability, and (c) the spatial pattern of the temporal stability. Results showed that temporal stability of soil water content at 0.2 m was significantly weaker than those at the soil depths of 0.6 and 0.8 m. Soil texture can significantly (P<0.05) affect the stability of soil water content except for the existence of an insignificant difference between sandy loam and silt loam textures, while temporal stability of areas covered by bunge needlegrass land was not significantly different from those covered by korshinsk peashrub. Geostatistical analysis showed that the temporal stability was spatially variable in an organized way as inferred by the degree of spatial dependence index. With increasing soil depth, the range of both temporal stability indices showed an increasing trend, being 65.8-120.5 m for SDRD and 148.8-214.1 m for MABE, respectively. This study provides a valuable support for soil water content measurements for soil water management and hydrological applications on sloping land areas. (C) 2010 Elsevier B.V. All rights reserved.

Impacts of land-use pattern on soil water-content variability on the Loess Plateau of China

Long-term vegetation restoration carried out on the slopes of the Loess Plateau of China employed different spatial and temporal land-use patterns but very little is known about the effects of these patterns on soil water-content variability. For this study the small Donggou catchment was selected to investigate soil water-content distributions for three spatial scales, including the entire catchment area, sampling transects, and land-use systems. Gravimetric soil water contents were determined incrementally to a soil depth of 1.20 m, on 10 occasions from April to October, 2007, at approximately 20-day intervals. Results indicated that soil water contents were affected by the six land-use types, resulting in four distinct patterns of vertical distribution of soil moisture (uniform, increasing, decreasing, and fluctuating with soil depth). The soil water content and its variation were also influenced in a complex manner by five land-use patterns distributed along transects following the gradients of five similar slopes. These patterns with contrasting hydrological responses in different components, such as forage land (alfalfa)-cropland-shrubland or shrubland-grassland (bunge needlegrass)-cropland-grassland, showed the highest soil water-content variability. Soil water at the catchment scale exhibited a moderate variability for each measurement date, and the variability of soil water content decreased exponentially with increasing soil water content. The minimum sample size for accurate data for use in a hydrological model for the catchment, for example, required many more samples for drier (69) than for wet (10) conditions. To enhance erosion and runoff control, this study suggested two strategies for land management: (i) to create a mosaic pattern by land-use arrangement that located units with higher infiltration capacities downslope from those with lower soil infiltrabilities; and (ii) raising the soil-infiltration capacity of units within the spatial mosaic pattern where possible.

On the use of correlated beta random variables with animal population modelling

We give reasons why demographic parameters such as survival and reproduction rates are often modelled well in stochastic population simulation using beta distributions. In practice, it is frequently expected that these parameters will be correlated, for example with survival rates for all age classes tending to be high or low in the same year. We therefore discuss a method for producing correlated beta random variables by transforming correlated normal random variables, and show how it can be applied in practice by means of a simple example. We also note how the same approach can be used to produce correlated uniform triangular, and exponential random variables. (C) 2008 Elsevier B.V. All rights reserved.

Modelling the potential distribution of the invasive tomato red spider mite, Tetranychus evansi (Acari: Tetranychidae)

Predicting the potential geographical distribution of a species is particularly important for pests with strong invasive abilities. Tetranychus evansi Baker & Pritchard, possibly native to South America, is a spider mite pest of solanaceous crops. This mite is considered an invasive species in Africa and Europe. A CLIMEX model was developed to predict its global distribution. The model results fitted the known records of T. evansi except for some records in dry locations. Dryness as well as excess moisture stresses play important roles in limiting the spread of the mite in the tropics. In North America and Eurasia its potential distribution appears to be essentially limited by cold stress. Detailed potential distribution maps are provided for T. evansi in the Mediterranean Basin and in Japan. These two regions correspond to climatic borders for the species. Mite establishment in these areas can be explained by their relatively mild winters. The Mediterranean region is also the main area where tomato is grown in open fields in Europe and where the pest represents a threat. According to the model, the whole Mediterranean region has the potential to be extensively colonized by the mite. Wide expansion of the mite to new areas in Africa is also predicted. Agricultural issues highlighted by the modelled distribution of the pest are discussed.

Parametric and nonparametric statistical modelling of crop yield: implications for pricing crop insurance contracts

This article considers alternative methods to calculate the fair premium rate of crop insurance contracts based on county yields. The premium rate was calculated using parametric and nonparametric approaches to estimate the conditional agricultural yield density. These methods were applied to a data set of county yield provided by the Statistical and Geography Brazilian Institute (IBGE), for the period of 1990 through 2002, for soybean, corn and wheat, in the State of Paran. In this article, we propose methodological alternatives to pricing crop insurance contracts resulting in more accurate premium rates in a situation of limited data.

Modelling the progress of Asiatic citrus canker on Tahiti lime in relation to temperature and leaf wetness

The combined effect of temperature (15A degrees C, 20A degrees C, 25A degrees C, 30A degrees C, 35A degrees C, 40A degrees C and 42A degrees C) and leaf wetness duration (0, 4, 8 12, 16, 20 and 24 h) on infection and development of Asiatic citrus canker (Xanthomonas citri subsp. citri) on Tahiti lime plant was examined in growth chambers. No disease developed at 42A degrees C and zero hours of leaf wetness. Periods of leaf wetness as short as 4 h were sufficient for citrus canker infection. However, a longer leaf duration wetness (24 h) did not result in much increase in the incidence of citrus canker, but led to twice the number of lesions and four times the disease severity. Temperature was the greatest factor influencing disease development. At optimum temperatures (25-35A degrees C), there was 100% disease incidence. Maximum disease development was observed at 30-35A degrees C, with up to a 12-fold increase in lesion density, a 10-fold increase in lesion size and a 60-fold increase in disease severity.