A micromechanical study of rolling and sliding contacts in assemblies of oval granules

The evolution of the microstructure of an assembly of cohesionless granular materials with associated pores, which carry the overall applied stresses through frictional contacts is a complex phenomenon. The macroscopic flow of such materials take place by the virtue of the relative rolling and sliding of the grains on the micro‐scale. A new discrete element method for biaxial compression simulations of random assemblies of oval particles with mixed sizes is introduced. During the course of deformation, the new positions of the grains are determined by employing the static equilibrium equations. A key aspect of the method is that, it is formulated for ellipse cross‐sectional particles, hence desirable inherent anisotropies are possible. A robust algorithm for the determination of the contact points between neighbouring grains is given. Employing the present methodology, many aspects of the behaviour of two‐dimensional assemblies of oval cross‐sectional rods have been successfully addressed. The effects of initial void ratio, interparticle friction angle, aspect ratio, and bedding angle on the rolling and sliding contacts are examined. The distribution of normals to the rolling and sliding contacts have different patterns and are concentrated along directions, which are approximately perpendicular to one another. On the other hand, the distribution of all contact normals (combined rolling and sliding) are close to that of rolling contacts, which confirm that rolling is the dominant mechanism. This phenomenon becomes more pronounced for higher intergranular friction angle. Characteristics of the rolling and sliding contacts are also discussed in the context of the force angle, which is the inclination of contact force with respect to the contact normal. Copyright © 2003 John Wiley & Sons, Ltd.     

Determining the groundwater potential recharge zone and karst springs catchment area: Saldoran region,western Iran

Assessing the groundwater recharge potential zone and differentiation of the spring catchment area are extremely important to effective management of groundwater systems and protection of water quality. The study area is located in the Saldoran karstic region, western Iran. It is characterized by a high rate of precipitation and recharge via highly permeable fractured karstic formations. Pire-Ghar, Sarabe-Babaheydar and Baghe-rostam are three major karstic springs which drain the Saldoran anticline. The mean discharge rate and electrical conductivity values for these springs were 3, 1.9 and 0.98 m³/s, and 475, 438 and 347 μS/cm, respectively. Geology, hydrogeology and geographical information system (GIS) methods were used to define the catchment areas of the major karstic springs and to map recharge zones in the Saldoran anticline. Seven major influencing factors on groundwater recharge rates (lithology, slope value and aspect, drainage, precipitation, fracture density and karstic domains) were integrated using GIS. Geology maps and field verification were used to determine the weights of factors. The final map was produced to reveal major zones of recharge potential. More than 80 % of the study area is terrain that has a recharge rate of 55–70 % (average 63 %). Evaluating the water budget of Saldoran Mountain showed that the total volume of karst water emerging from the Saldoran karst springs is equal to the total annual recharge on the anticline. Therefore, based on the geological and hydrogeological investigations, the catchment area of the mentioned karst springs includes the whole Saldoran anticline.     

Application of the Christensen Failure Criterion to Intact Rock

The Christensen criterion, originally introduced in materials science, has a simple mathematical form and uniaxial tensile and compressive strength as the only parameters, making it an attractive candidate for rock engineering purposes. In this study, the applicability of the criterion to rock materials is examined. Explicit equations for application of the criterion under biaxial, triaxial compression, triaxial extension, and polyaxial states of stresses are derived. A comprehensive strength data set including the results of tests on synthetic rock, chert dyke, Carrara marble and Westerly granite is utilized to examine the accuracy of the Christensen criterion to the failure of rock material. The two surprising findings about the Christensen criterion are the zero values of tensile strength and the very low slopes of the failure envelope obtained from fitting analyses for chert dyke and Westerly granite. It is shown that the two problems are interrelated and the values of tensile strength tend to zero to produce higher slopes. It is then mathematically proven that the maximum initial slope of the Christensen failure envelope is limited to 4 in triaxial compression and 2.5 in triaxial extension which is considerably lower than the slope of experimental data. The accuracy of the Christensen criterion was found to be significantly lower than the well-established Hoek–Brown criterion. The circular π-plane representations and brittle-to-ductile transition limits from the Christensen criterion are also inconsistent with the observed behavior of rocks.     

Simulation of Power Production from Dry Geothermal Well Using Down-hole Heat Exchanger in Sabalan Field,Northwest Iran

In this research, a simulation was performed for evaluating power production from an abandoned geothermal well as an enhanced geothermal system by injecting a secondary fluid. Abandoned wells, due to lack of fluid or very low transmissivity, are regarded among the low-to moderate-temperature resources that have the potential for heat production without any cost for deep drilling. Accordingly, they are taken as suitable sources of energy. In the present paper, an abandoned geothermal well at Meshkinshahr geothermal field in Sabalan district, northwestern Iran, with 3176 m depth was simulated. The bottom-hole temperature of 148 °C, as well as well casing size, and real thermal gradient for well were applied in the model. A 3D heat transfer simulation model was designed by considering a coaxial pipe as a down-hole heat exchanger between surrounding rocks of the well and injected fluid. Injected fluid to the well with specified pressure and temperature receives heat from rocks surrounding the well, until it reaches the bottom of the well and converts to vapor. The vapor returns to the surface from inner pipe with very low heat loss during its return. The inner pipe is isolated by a thin layer having a low heat conductivity to prevent heat loss from the returned fluid. It was observed that obtained heat in the well depends on temperature profile of the well, injection velocity, and fluid mass flow rate. The model results were optimized by selecting suitable parameters such as inlet injection speed and fluid flow rate to achieve the highest temperature of the fluid returned from the well. A binary power plant was also modeled to determine the extractable power using returned fluid as input using ammonia and isobutene, as working fluids in binary cycle. Finally, electric power of 270 kW was generated from well NWS3 using designed down-hole heat exchanger.     

Anisotropic behavior of the Bushehr carbonate sand in the Persian Gulf

The anisotropy effect is exhibited more prominently in sedimentary depositions, and it relates the soil’s mechanical specifications to the directions of imposed loads. Even though this phenomenon has been comprehensively explored in silica sands, few research has been conducted for studying the anisotropic behavior of marine carbonate sands. To bridge this gap, the present study investigates the anisotropy effect on the mechanical behavior of Bushehr carbonate sand acquired from the north shelf of the Persian Gulf in Iran. Toward this end, some undrained principal stress rotation tests are conducted using a hollow cylinder shear torsional apparatus in such a manner that the direction of the applied principal stresses are fixed along a desired orientation and the total mean stress and intermediate principal stress ratio are kept constant. Furthermore, prior to shearing, the samples are consolidated under three confining pressures and two isotropic and anisotropic states. The results show that dilative behavior is observed in all loading directions after initial contraction; this contradicts the response observed in silica sands. The anisotropy response of soil follows two different trends in the contractive and dilative phases. The relation of soil’s mechanical properties shows a descending trend with the angle of maximum principal stress in the contractive phase; on the other hand, the anisotropy behavior shows a dominant parabola trend in the dilative phase, where the maximum ultimate pore pressure and minimum soil strength occur in the stress direction with an angle of α?=?30°. By increasing the confining pressure in the soil element, the intensity of the anisotropy in some mechanical properties except the soil deformation is reduced. Furthermore, the deviatoric-to-effective mean stress ratio in the phase transformation state from contraction to dilation is independent of the loading direction and consolidation stress state, and it is considered one of the intrinsic properties of sand.     

Performance evaluation of the M8 algorithm to predict M7+ earthquakes in Turkey

The M8 algorithm is one of the most reliable intermediate-term middle-range earthquake prediction algorithms. The present study evaluates the ability of the M8 algorithm and its modified versions for predicting major events (M7+) in Turkey. Thirty different algorithms were developed by changing the radius of circle of investigation (CI) and the lower magnitude cutoff of the M8 algorithm. These modified algorithms were executed all over the territory of Turkey, and the results were evaluated using the error diagram. Each modified algorithm was executed for consecutive half-year intervals over a specified period of time. Subsequently, the seismic catalog was updated, and failures-to-predict ratio and the fraction of alarm were considered. Results showed that the location of areas of alarm change gradually over consecutive intervals, and no sudden changes can be observed. In addition, the annual changes of areas of alarm are not random and follow a pattern. This study also showed that the modified algorithm having a three to six annual average of events and a 393-km CI radius is an efficient algorithm for predicting the future seismic events in Turkey. This algorithm predicted six out of six target events, retrospectively, with a confidence level of 96.4 %. According to the obtained results, it will be possible to rely on this modified algorithm to predict near future earthquakes of Turkey. Furthermore, this study proves that it is possible to alter the M8 algorithm for being used in regional studies.     

Simulation of rainfall-runoff process using an artificial neural network (ANN) and field plots data

Theoretical and Applied Climatology - Rainfall-runoff modeling is necessary for many hydrological studies, such as estimating peak discharges and designing hydraulic structures. The intensity and...     

Rock Type Identification Using Analysis of the Acoustic Signal Frequency Contents Propagated While Drilling Operation

Geotechnical and Geological Engineering - One of the most important parameters specially in mining and oil drilling fields is the type of rocks. It is important to determine that rock structure is...