Development of stochastic dynamic Nash game model for reservoir operation. I. The symmetric stochastic model with perfect information

Increasing water demands, higher standards of living, depletion of resources of acceptable quality and excessive water pollution due to agricultural and industrial expansions have caused intense social and political predicaments, and conflicting issues among water consumers. The available techniques commonly used in reservoir optimization/operation do not consider interaction, behavior and preferences of water users, reservoir operator and their associated modeling procedures, within the stochastic modeling framework. In this paper, game theory is used to present the associated conflicts among different consumers due to limited water. Considering the game theory fundamentals, the Stochastic Dynamic Nash Game with perfect information (PSDNG) model is developed, which assumes that the decision maker has sufficient (perfect) information regarding the associated randomness of reservoir operation parameters. The simulated annealing approach (SA) is applied as a part of the proposed stochastic framework, which makes the PSDNG solution conceivable. As a case study, the proposed model is applied to the Zayandeh-Rud river basin in Iran with conflicting demands. The results are compared with alternative reservoir operation models, i.e., Bayesian stochastic dynamic programming (BSDP), sequential genetic algorithm (SGA) and classical dynamic programming regression (DPR). Results show that the proposed model has the ability to generate reservoir operating policies, considering interactions of water users, reservoir operator and their preferences.     

Stability of buoyancy opposed mixed convection in a vertical channel and its dependence on permeability

Non-Darcy mixed convective flow of water due to external pressure gradient and buoyancy opposed forces are considered in a vertical channel filled with porous medium, which can be either isotropic or anisotropic. The linear theory of stability analysis has been used to numerically investigate the dependence of the transition behavior of the fully developed basic flow on the permeability of the medium. Numerical experiments indicate that mainly two main instability modes appear: Rayleigh–Taylor (R–T) and buoyant instability. For Darcy numbers (Da) ?10⁻⁹, R–T instability dominates within the entire Reynolds number (Re) range considered here. It was also found that for the same Re, the fully developed base flow is highly unstable (stable) for porous media with high (low) permeability. Further, it was seen that the disturbance isotherm cells migrate from the channel walls toward the centerline when permeability is reduced. Reducing the permeability by one order of magnitude (corresponding to a decrease of Darcy number from 10⁻⁶ to 10⁻⁷) increases base flow stability approximately 20-fold. For higher Reynolds numbers, buoyant, mixed and shear instability of the basic flow were found when Da was increased from 10⁻⁷ to 10⁻³. However, for cases in which permeability and porosity behaved as suggested by Carman–Kozeny relation (CKR), buoyant stability was the only mode of instability. Critical values of the Rayleigh (Ra) and Darcy (Da) numbers in the R–T mode of instability were related to each other by the hyperbolic function RaDa = −2.465.     

Mineralogical and geochemical investigation of some pegmatites from the Mashad area,northeastern Iran

Summary Pegmatitic dikes and masses of various sizes occur in a granite intrusion in the Mashad area, northeastern Iran. Their geochemistry, as well as that of the parent granite was investigated using some 40 specimens. Major constituents and trace elements were determined by means of X-ray fluorescence and spectrographic methods. Anomalously high concentrations of Be and Nb and indications of their minerals beryl and columbite were observed. Therefore, the pegmatites of this area may prove to be a potential source for Be, Nb, Ta and possibly Li.

---

Mineralogische und geochemische Untersuchung einiger Pegmatite aus der Gegend von Mashad, Nordost-Iran

Zusammenfassung In einer Granitintrusion in der Gegend von Mashad, Nordost-Iran, kommen pegmatitische Gänge und Massen verschiedener Größe vor. Unter Verwendung von etwa 40 Proben wurde ihre Geochemie und die ihres Muttergranites untersucht. Die Hauptbestandteile und Spurenelemente wurden mit Röntgenfluoreszenz und spektrographischen Methoden bestimmt. Außergewöhnlich hohe Konzentrationen an Be und Nb und Hinweise auf ihre Minerale Beryll und Columbit wurden beobachtet. Die Pegmatite dieser Gegend könnten sich folglich als potentielle Lieferanten für Be, Nb, Ta und möglicherweise Li herausstellen.

---

---

With 1 Figure     

Structural Control on the Salmas Geothermal Region, Northwest Iran, from Fractal Analysis and Paleostress Data

The Salmas geothermal field is located in NW Iran. Subduction of Neo-Tethys oceanic crust beneath the Iranian microcontinent caused to propagation of the magmatic-Arc. Fractures and faults in the convergent zone have created path-ways for the circulation of geothermal fluid. Fracture concentration in the Salmas geothermal field has been characterized using of the fractal method and creation of a fracture density map that shows the highest concentration in the central part of the study area. The permeability of fractures has been evaluated by analyzing their orientation in respect to the paleostress axes. Also, the fractal analyzing result indicates the maximum fractal dimension (1.96) is around the thermal spring outlet. Paleostress analyzing revealed that in the central part of the study area, σ1 axes orientation is S90°W/10° and the σ2 dip is near to the vertical in this stress field, where strike slip faults can be propagated. In the SE part near the recharge of the thermal springs, the σ3 plunge increases to 70? and σ1 orientation is N15°E/20°, in this local tectonic regime thrust fault developed. Fractures have an important role in the circulation of fluid and the fractal dimension increases near the thermal springs in the Salmas geothermal field. Regarding the paleostress data fracture with N-S direction such as the F1 fault zone (parallel to the σ1 axes), a suitable pathway for deep circulation of geothermal fluid flow has been created.     

Design of raft–pile foundation using combined optimization and finite element approach

This paper describes the application of a structural optimization approach combined with the finite element method for the optimal design of a raft–pile foundation system. The analysis takes into account the non-linear behaviour of the soil medium and the piles. For the optimization process, the sensitivity analysis is carried out using the approximate semi-analytical method while the constraint approximation is obtained from the combination of extended Bi-point and Lagrangian polynomial approximation methods. The objective function of the problem is the cost of the foundation. The design variables are the raft thickness, cross-section, length and number of piles. The maximum displacement and differential displacement are selected as the constraints. The proposed method is shown to achieve an optimum design of raft–pile foundation efficiently and accurately. Copyright © 1999 John Wiley & Sons, Ltd.     

A peridynamics model for strain localization analysis of geomaterials

Geomaterials such as sand and clay are highly heterogeneous multiphase materials. Nonlocality (or a characteristic length scale) in modeling geomaterials based on the continuum theory can be associated with several factors, for instance, the physical interactions of material points within finite distance, the homogenization or smoothing process of material heterogeneity, and the particle or problem size-dependent mechanical behavior (eg, the thickness of shear bands) of geomaterials. In this article, we formulate a nonlocal elastoplastic constitutive model for geomaterials by adapting a local elastoplastic model for geomaterials at a constant suction through the constitutive correspondence principle of the state-based peridynamics theory. We numerically implement this nonlocal constitutive model via the classical return-mapping algorithm of computational plasticity. We first conduct a one-dimensional compression test of a soil sample at a constant suction through the numerical model with three different values of the nonlocal variable (horizon) δ. We then present a strain localization analysis of a soil sample under the constant suction and plane strain conditions with different nonlocal variables. The numerical results show that the proposed nonlocal model can be used to simulate the inception and propagation of shear banding as well as to capture the thickness of shear bands in geomaterials at a constant suction.     

An evolutionary approach to modelling the thermomechanical behaviour of unsaturated soils

A new data‐mining approach is presented for modelling of the stress–strain and volume change behaviour of unsaturated soils considering temperature effects. The proposed approach is based on the evolutionary polynomial regression (EPR), which unlike some other data‐mining techniques, generates a transparent and structured representation of the behaviour of systems directly from raw experimental (or field) data. The proposed methodology can operate on large quantities of data in order to capture nonlinear and complex relationships between contributing variables. The developed models allow the user to gain a clear insight into the behaviour of the system. Unsaturated triaxial test data from the literature were used for development and verification of EPR models. The developed models were also used (in a coupled manner) to produce the entire stress path of triaxial tests. Comparison of the EPR model predictions with the experimental data revealed the robustness and capability of the proposed methodology in capturing and reproducing the constitutive thermomechanical behaviour of unsaturated soils. More importantly, the capability of the developed models in accurately generalizing the predictions to unseen data cases was illustrated. The results of a sensitivity analysis showed that the models developed from data are able to capture and represent the physical aspects of the unsaturated soil behaviour accurately. The merits and advantages of the proposed methodology are also discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

Hot springs mediate spatial exchange of heat and mass in the enclosed sediment domain: A stability perspective

In many natural environments, such as in underwater hot springs and hydrothermal vents, thermal gradients are accompanied with changes in the concentration of chemical compounds transported to the seawater, causing the so-called double-diffusive, mixed convection. To study the physical scenarios in such systems, a vertical channel filled with a porous medium saturated with saline water is considered. The motion in the sediment-filled channel is induced by two buoyancy forces and an external pressure gradient, similar to the situation in a vent with an upward flow direction. The fluid flow has been modeled by an extended Darcy model, and the flow instability mechanisms have been studied numerically. The linear stability analysis is performed considering a wide range of Darcy number (Da = 10⁻⁵ -10⁻⁸). The instability boundary curve showed three distinct dynamic regimes: (i) Rayleigh-Taylor (R-T), (ii) log-log non-linear variation, and (iii) log-log linear variation. The domain of different regimes were sensitive to external pressure gradient as well as permeability. Similar to cross-diffusive natural convection in pure viscous fluids, a linear relationship between logarithmic absolute values of critical thermal Rayleigh number (∣Ra_T∣) and solute Rayleigh number (Ra_C) is found in the third regime. Based on the permeability, for any solute Rayleigh number (Ra_C), there existed a minimum value of Reynolds number (Re), below which R-T type of instability appeared. Above this minimum value, the instability was due to two buoyancy forces, known as buoyant instability. Simulations of secondary flow via energy analysis demonstrated the development of complex dynamics at the critical state in all three regimes characterized by transition of multi to uni-cellular structures and vice verse.     

Enhancement of Carbon Dioxide Capture by Amine‐Functionalized Multi‐Walled Carbon Nanotube

In this study, crude multi‐walled carbon nanotubes (MWCNT) was functionalized by a two‐step process; first using strong mixed acids (H₂SO₄/HNO₃) and then treatment with 1,3‐phenylenediamine (mPDA). The equilibrium adsorption of CO₂ on pristine MWCNT and amine functionalized MWCNT (MWCNT‐NH₂) were investigated. Experiments were preformed via application of volumetric method in a dual sorption vessel at temperature range of 298–318 K and pressures up to 40 bars. The results obtained indicated that the equilibrium uptake of CO₂ increased after functionalizing of MWCNT. The increase in CO₂ capture by MWCNT‐NH₂ was attributed to the existence of great affinity between CO₂ molecules and amine sites on this adsorbent at low pressures. The experimental data were analyzed by means of Freundlich and Langmuir adsorption isotherm models. The data obtained revealed a fast kinetics for the adsorption of CO₂ in which most of adsorption occurred at initial period of adsorption experiments. This renders MWCNT as a suitable adsorbent for practical applications. Values of isosteric heat of adsorption were evaluated based on Clausius–Clapeyron equation. The results demonstrated that both chemisorption and physisorption played important role in CO₂ adsorption on MWCNT‐NH₂, whereas the physisorption process was dominant for CO₂ adsorption on MWCNT.