Limiting equilibrium method for slope/drilled shaft system

The use of drilled shafts to stabilize an unstable slope has been a widely accepted practice. There are two basic design and analysis issues involved: one is to determine the global factor of safety of the drilled shafts stabilized slope and the other one is to determine the design earth thrust on the drilled shafts for structural design of the shafts. In this paper, a limiting equilibrium method of slices based solution for calculating global factor of safety (FS) of a slope with the presence of a row of drilled shafts is developed. The arching mechanisms due to the presence of the drilled shafts on slope were taken into account by a load transfer factor. The method for calculating the net force applied to the drilled shaft from the soil mass was also developed. The interrelationships among the drilled shaft location on the slope, the load transfer factor, and the global FS of the slope/shaft system were derived utilizing the developed numerical closed‐form solution. An illustrative example is presented to elucidate the use of the solution in optimizing the location of the drilled shafts on slope to achieve the desired global factor of safety of the slope/shaft system. Copyright © 2009 John Wiley & Sons, Ltd.     

Mineral chemical composition and geodynamic significance of peridotites from Nain ophiolite,central Iran

The peridotites from north of the town of Nain in central Iran consist of clinopyroxene-bearing harzburgite and lherzolite with small lenses of dunite and chromitite pods. The lherzolite contains aluminous spinel with a Cr number (Cr^# = Cr/[Cr + Al]) of 0.17. The Cr number of spinels in harzburgite and chromitite is 0.38–0.42 and 0.62, respectively. This shows that the lherzolite and harzburgite resulted from <18% of partial melting of the source materials. The estimated temperature is 1100 ± 200 °C for peridotites, the estimated pressure is <15 ± 2.3 kbar for harzburgites and >16 ± 2.3 kbar for lherzolites and estimated fo₂ is 10^?1±0.5 for peridotites. Discriminant geochemical diagrams based on mineral chemistry of harzburgites indicate a supra-subduction zone (SSZ) to mid-oceanic ridge (MOR) setting for these rocks. On the basis of their Cr^#, the harzburgite and lherzolite spinels are analogous to those from abyssal peridotites and oceanic ophiolites, whereas the chromites in the chromitite (on the basis of Cr^# and boninitic nature of parental melts) resemble those from SSZ ophiolitic sequences. Therefore, the Nain ophiolite complex most likely originated in an oceanic crust related to supra-subduction zone, i.e. back arc basin. Field observations and mineral chemistry of the Nain peridotites, indicating the suture between the central Iran micro-continent (CIM) block and the Sanandaj–Sirjan zone, show that these peridotites mark the site of the Nain–Baft seaway, which opened with a slow rate of ocean-floor spreading behind the Mesozoic arc of the Sanandaj–Sirjan zone as a result of change of Neo Tethyan subduction régime during middle Cretaceous.     

Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.     

List of forthcoming papers

List of forthcoming papers     

Prediction of Land Subsidence Under Cyclic Pumping Based on Laboratory and Numerical Simulations

In this study subsidence due to groundwater withdrawal was investigated. Kerman Province in Iran is struggling with land subsidence problem due to extensive groundwater withdrawal mainly for farming. The rate and type of groundwater withdrawal has very important impact on settlement rate. In this research, effective parameters on land subsidence caused by groundwater withdrawal were determined by laboratory tests. Sampling had done up to depth of 300 m mainly with remolded specimens from Shams-abad, Nouq plain in Kerman province. Similar to the field preconsolidation pressure was applied on specimens in the laboratory. Rate of applied stress on prepared specimens was similar to effect of oscillation of groundwater level. In order to model the actual soil behavior in the laboratory, one-dimensional consolidation device (odometer) was adopted for testing. In these tests, the effect of loading caused by seasonal oscillation of groundwater table is considered by means of cyclic loading in the testing which has great effect on rate of settlements. The results of tests show that when the water table level periodically increases and decreases the amount of settlement decrease, comparing with the case when the groundwater table drop to a constant level. In order to predict the further effects of groundwater level oscillation and actual field condition on land subsidence, a finite element model based on Biots’ three-dimensional consolidation theory was developed. After calibration of finite element model with laboratory tests, this model was used for prediction the effect of groundwater level oscillation on actual field conditions.     

Development of a Darcy-Brinkman model to simulate water flow and tracer transport in a heterogeneous karstic aquifer (Val d’Orléans, France)

Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10^?6 to 5.5?×?10^?5?m² and the limestone permeability ranges from 8?×?10^?11 to 6?×?10^?10?m². The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated.     

Hydrochemical characteristics and quality assessment of shallow groundwater in a coastal area of Southwest Bangladesh

In this study, the hydrochemical characteristics of shallow groundwater in a coastal region (Khulna) of southwest Bangladesh have been evaluated based on different indices for drinking and irrigation uses. Water samples were collected from 26 boreholes and analyzed for major cations and anions. Other physico-chemical parameters like pH, electrical conductivity (EC), and total dissolved solids were also measured. Most groundwater is slightly alkaline and largely varies in chemical composition, e.g. EC ranges from 962 to 9,370 μs/cm. The abundance of the major ions is as follows: Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ = Cl⁻ > HCO₃ ⁻ > SO₄ ²⁻ > NO₃ ⁻. Interpretation of analytical data shows two major hydrochemical facies (Na⁺–K⁺–Cl⁻–SO₄ ²⁻ and Na⁺–K⁺–HCO₃ ⁻) in the study area. Salinity, total hardness, and sodium percentage (Na%) indicate that most of the groundwater samples are not suitable for irrigation as well as for domestic purposes and far from drinking water standard. Results suggest that the brackish nature in most of the groundwaters is due to the seawater influence and hydrogeochemical processes.     

Spatiotemporal variation and mechanisms of temperature inversion in the Bay of Bengal and the eastern equatorial Indian Ocean

In the northern Bay of Bengal, the existence of intense temperature inversion during winter is a widely accepted phenomenon. However, occurrences of temperature inversion during other seasons and the spatial distribution within and adjacent to the Bay of Bengal are not well understood. In this study, a higher resolution spatiotemporal variation of temperature inversion and its mechanisms are examined with mixed layer heat and salt budget analysis utilizing long-term Argo（2004 to 2020） and RAMA（2...     

Mixed precision parallel preconditioner and linear solver for compositional reservoir simulation

Reservoir simulation role in value creation and strategic management decisions cannot be over emphasized. Simulation of complex challenging reservoirs with millions of grid blocks especially in compositional mode is very time-consuming even with fast modern computers. On the other hand, high price of cluster supercomputers prevents them for being commonly used for fast simulation of such reservoirs. In recent years, the development of many-core processors like cell processors, DSPs, and graphical processing units (GPUs) has provided a very cost-effective hardware platform for fast computational operations. However, programming for such processors is much more difficult than conventional CPUs, and new parallel algorithm design and special parallel implementation methods are needed. Using the computational power of CPUs, GPUs, and/or any other processing unit, Open Computing Language (OpenCL) provides a framework for programming for heterogeneous platforms. In this paper, OpenCL is used to employ the computational power of a GPU to build a preconditioner and solve the linear system arising from compositional formulation of multiphase flow in porous media. The proposed parallel preconditioner is proved to be quite effective, even in heterogeneous porous media. Using data-parallel modules on GPU, the preconditioner/solver runtime reduced at least 1 order of magnitude compared to their serial implementation on CPU.     

Optimal water and waste load allocation in reservoir–river systems: a case study

In this paper, a new methodology is developed for optimization of water and waste load allocation in reservoir–river systems considering the existing uncertainties in reservoir inflow, waste loads and water demands. A stochastic dynamic programming (SDP) model is used to optimize reservoir operation considering the inflow uncertainty, and another model called PSO-SA is developed and linked with the SDP model for optimizing water and waste load allocation in downstream river. In the PSO-SA model, a particle swarm optimization technique with a dynamic penalty function for handling the constraints is used to optimize water and waste load allocation policies. Also, a simulated annealing technique is utilized for determining the upper and lower bounds of constraints and objective function considering the existing uncertainties. As the proposed water and waste load allocation model has a considerable run-time, some powerful soft computing techniques, namely, Regression tree Induction (named M5P), fuzzy K-nearest neighbor, Bayesian network, support vector regression and an adaptive neuro-fuzzy inference system, are trained and validated using the results of the proposed methodology to develop real-time water and waste load allocation rules. To examine the efficiency and applicability of the methodology, it is applied to the Dez reservoir–river system in the south-western part of Iran.