Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.     

An infinity of conservation laws of fKdV equation

An infinity of conservation laws of fKdV equation is derived in terms of the Miura and Gardner''s transform.The pseudo-mass and energy theorems are studied by the first two conservation laws.As a typical example,the theoretical mean wave resistance and the regional distribution of energy of the precursor soliton generation are determined by means of the first and the second conservation laws.     

Constructal view of scaling laws of river basins

River basins are examples of naturally organized flow architectures whose scaling properties have been noticed long ago. Based on data of geometric characteristics, Horton [Horton, R.E., 1932. Drainage basin characteristics. EOS Trans. AGU 13, 350–361.], Hack [Hack, J.T., 1957. Studies of longitudinal profiles in Virginia and Maryland. USGS Professional Papers 294-B, Washington DC, pp. 46–97.], and Melton [Melton, M.A, 1958. Correlation structure of morphometric properties of drainage systems and their controlling agents. J. of Geology 66, 35–56.] proposed scaling laws that are considered to describe rather accurately the actual river basins. What we show here is that these scaling laws can be anticipated based on Constructal Theory, which views the pathways by which drainage networks develop in a basin not as the result of chance but as flow architectures that originate naturally as the result of minimization of the overall resistance to flow (Constructal Law).     

Division of Gas Accumulation System and Laws Controlling Distribution of Natural Gas in Typical Petroliferous Basins of China

Considering the existing problems of the petroleum system, this paper brings forward the concept of natural gas accumulation system and presents the dividing principles. Then detailed statistics on the accumulation factors of the 32 typical natural gas accumulation systems in China and studies on the laws controlling distribution of gas are collected. The research shows that the petroleum accumulation system is the basic unit controlling petroleum generation, migration and accumulation. Generating intensity, generating amount, accumulating efficiency and migration distance plays an important role in the distribution of natural gas. Through analysis on results of resources evaluation, discovered reserves and residual reserves, potential areas in middle-scaled petroliferous basins in China are forecasted in this paper. Ordos, Sichuan, Tarim and Qaidam basins are found out to be the main basins developing and enriching gas accumulation systems.     

Scaling effects on modeled surface energy-balance components using the NOAH-OSU land surface model

As surface exchange processes are highly non-linear and heterogeneous in space and time, it is important to know the appropriate scale for the reasonable prediction of these exchange processes. For example, the explicit representation of surface variability has been vital in predicting mesoscale weather events such as late-afternoon thunderstorms initiated by latent heat exchanges in mid-latitude regions of the continental United States. This study was undertaken to examine the effects of different spatial scales of input data on modeled fluxes, so as to better understand the resolution needed for accurate modeling. A statistical procedure was followed to select two cells from the Southern Great Plains 1997 hydrology experiment region, each 20 km×20 km, representing the most homogeneous and the most heterogeneous surface conditions (based on soil and vegetation) within the study region. The NOAH-OSU (Oregon State University) Land Surface Model (LSM) was employed to estimate surface energy fluxes. Three scales of study (200 m, 2 and 20 km) were considered in order to investigate the impacts of the aggregation of input data, especially soil and vegetation inputs, on the model output. Model results of net radiation and latent, sensible and ground heat fluxes were compared for the three scales. For the heterogeneous area, the model output at the 20-km resolution showed some differences when compared with the 200-m and 2-km resolutions. This was more pronounced in latent heat (12% decrease), sensible heat (22% increase), and ground heat flux (44% increase) estimation than in net radiation. The scaling effects were much less for the relatively homogeneous land area with 5% increase in sensible heat and 4% decrease in ground heat flux estimation. All of the model outputs for the 2- and 20-km resolutions were in close agreement. The results suggested that, for this study region, soils and vegetation input resolution of about 2 km should be chosen for realistic modeling of surface exchange processes. This resolution was sufficient to capture the effects of sub-grid scale heterogeneity, while avoiding the data and computational difficulties associated with higher spatial resolutions.     

Measuring the k–S–p relations on DNAPLs migration

Migration characteristics of dense non-aqueous phase liquids (DNAPLs) in the subsurface can be expressed as a water–DNAPL two-phase system based on mainly the k–S–p relations, which describe the relations among relative permeability (k), degree of water saturation (S), and capillary pressure head (p). The aims of this research are to develop an experimental system with new type of probes, such as the electrical conductivity probe for measuring the degree of water saturation and the hydrophilic and hydrophobic tensiometer for measuring the pore water and DNAPL pressures, and also to estimate the results obtained from the developed experimental system. From these tests, the excellent k–S–p relations were obtained and the efficiency and validity of this developed experimental system have been confirmed in terms of the concept of the scaling coefficient calculated by interfacial tensions.     

Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows

We use high resolution Monte Carlo simulations to study the dispersive mixing in two-phase, immiscible, porous media flow that results from the interaction of the nonlinearities in the flow equations with geologic heterogeneity. Our numerical experiments show that distinct dispersive regimes occur depending on the relative strength of nonlinearity and heterogeneity. In particular, for a given degree of multiscale heterogeneity, controlled by the Hurst exponent which characterizes the underlying stochastic model for the heterogeneity, linear and nonlinear flows are essentially identical in their degree of dispersion, if the heterogeneity is strong enough. As the heterogeneity weakens, the dispersion rates cross over from those of linear heterogeneous flows to those typical of nonlinear homogeneous flows.     

中国北西部地区地浸砂岩型铀矿床的成矿条件与分布规律

中亚地区发现大量可地浸砂岩型铀矿床引起了我国的重视 ,目前正在积极开展此类型铀矿的找矿工作。我国北西部是寻找类似中亚地区层间氧化带后生地浸砂岩型铀矿床最理想的地区 ,本文在论述了区域大地构造特征、地壳演化与区域铀成矿带形成的基础上 ,着重阐明了区域中新生代盆地地浸砂岩型铀矿床形成的铀源条件、盆地类型与结构及富铀性对形成地浸砂岩型铀矿的影响 ;并以区内已有层间氧化带地浸砂岩型铀矿床为例 (伊犁盆地、二连盆地与测老庙盆地砂岩型铀矿床 ) ,论述了盆地砂岩型铀矿床形成的岩相古地理条件、水文地质与层间氧化带发育条件及对铀成矿的影响 ;并阐述了本区层间氧化带地浸砂岩型铀矿床的分布规律     

Hazard mapping based on macroseismic data considering the influence of geological conditions

The object of this study is to consider directly the influence of regional geological conditions on the assessment of seismic hazard. It is assumed that macroseismic data at individual locations contain, in an average way, the influence of geological conditions.A Data Base referring to 199 historical (5) and instrumental (194, in the 1947–1993 period) events with macroseismic information in 1195 locations of Portugal was built. For any given seismic event, whenever macroseismic information was available at a location (town, village, etc.), an EMS-92 intensity value was estimated. To each one of those locations a geological unit, representing the most common type of soil, was assigned, based on the Geological Portuguese Map at a scale 1:500 000; the geological units were grouped into three categories: soft, intermediate and hard soils.The Data Base was used to determine the attenuation laws in terms of macroseismic intensity for the three different geological site conditions, using multiple linear regression analysis. The reasonability of the laws was tested by (i) checking residual distributions and (ii) comparing the map of isoseismals of important earthquakes with the isoseismals generated by the attenuation curves derived for each one of the three different soil classes, taking into consideration the soil class of each site. The main results of attenuation modeling are: high dispersion on macroseismic intensity data; all the models predict intensity values, for short hypocentral distances, lower than the ones observed; and for some important analyzed earthquakes and for the observed range of distances, the models confirm the expectancy that macroseismic intensity increases from hard to soft soil.The approach to obtain the hazard assessment at each location consisted in the use of the attenuation law specifically derived for the class of soil of that particular location. This method, which considers the influence of the regional geology, was illustrated with the mapping of hazard for the country for several return periods. Comparison with previous maps not taking into consideration the regional geological conditions emphasizes the importance of this new parameter. It can be concluded that (i) soil segmentation is clearly the cause for hazard increase in the region to the north of Lisbon, especially at sites with soft and intermediate soils as the ones in lower Tagus valley; the maximum increase on hazard is, in any case, less than one degree; (ii) when geological conditions are disregarded in the attenuation regression analysis, hazard pattern is similar to the one obtained for the case of hard soil everywhere.