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.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

Radiative Transfer in Inhomogeneous Atmospheres. III

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.     

On cross-correlating weak lensing surveys

The present generation of weak lensing surveys will be superseded by surveys run from space with much better sky coverage and high level of signal-to-noise ratio, such as the Supernova/Acceleration Probe ( SNAP ). However, removal of any systematics or noise will remain a major cause of concern for any weak lensing survey. One of the best ways of spotting any undetected source of systematic noise is to compare surveys that probe the same part of the sky. In this paper we study various measures that are useful in cross-correlating weak lensing surveys with diverse survey strategies. Using two different statistics – the shear components and the aperture mass – we construct a class of estimators which encode such cross-correlations. These techniques will also be useful in studies where the entire source population from a specific survey can be divided into various redshift bins to study cross-correlations among them. We perform a detailed study of the angular size dependence and redshift dependence of these observables and of their sensitivity to the background cosmology. We find that one-point and two-point statistics provide complementary tools which allow one to constrain cosmological parameters and to obtain a simple estimate of the noise of the survey.     

Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

CATASTROPHIC ECO-ENVIRONMENTAL CHANGE IN THE SONGNEN PLAIN, NORTHEASTERN CHINA SINCE 1900S

Although the Songnen Plain in the northeastern China was developed relatively late in the temperate zone of the world, its eco-environment has changed greatly. This paper analyzes the changes of land cover and the rates and trends ofdesertification during the past 100 years in the Songnen Plain. According to the macroscopic analysis, we find that the eco-environment in the plain has reached to the threshold of catastrophic change since the 1950s. The Thorn Needle Catastrophic Model was used to determine and validate this conclusion. Human activities, including large-scale construction projects, such as huge dams and dikes, and excessive grazing were the primary factors contributing to regional eco-environmental catastrophe. And irrational reclamation of the wilderness also affected the eco-environmental change. The results reveal the complex human-land interactions.     

A MESO-β SCALE SIMULATION OF THE EFFECTS OF BOREAL FOREST ECOSYSTEM ON THE LOWER ATMOSPHERE

Based on the Intensive Field Campaign(IFC-1)data of Boreal Ecosystem-Atmosphere Study(BOREAS).a three-dimensional meso-β scale model is used to simulate the effect of boreal forests onthe lower atmosphere.A fine horizontal resolution of 2 km×2 km is used in order to distinguish thevegetative heterogeneity in the boreal region.A total of 20×25 grid points cover the entire sub-modeling area in BOREAS' South Study Area(SSA).The ecosystem types and their coverage ineach grid square are extracted from the North American Land Cover Characteristics Data Base(NALCCD)generated by the U.S.Geographical Survey(USGS)and the University of Nebraska-Lincoln(UNL).The topography of the study area is taken from the Digital Elevation Map(DEM)of USGS.The model outputs include the components of the energy balance budget within the canopyand at the ground.the turbulence parameters in the atmospheric boundary layer and the wind.temperature and humidity profiles extending up to a height of 1500 m.In addition to the fine timeand spatial step,the unique feature of the present model is the incorporation of both dynamic andbiological effects of the Boreal forest into the model parameterization scheme.The model resultscompare favorably with BOREAS' IFC-1 data in 1994 when the forest was in the luxuriant growingperiod.     

Forced oscillations in a hydrodynamical accretion disk and QPOs

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz and Klu’zniak (2001). In a first paper (P'etri, 2005a, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. After a discussion on the magnitude of this deformation applied to neutron stars, we show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to {three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric resonance}. In a second part, we focus on the linear response of a thin accretion disk in the 2D limit. {Waves are launched at the aforementioned resonance positions and propagate in some permitted regions inside the disk, according to the dispersion relation obtained by a WKB analysis}. In a last part, these results are confirmed and extended via non linear hydrodynamical numerical simulations performed with a pseudo-spectral code solving Euler's equations in a 2D cylindrical coordinate frame. {We found that for a weak potential perturbation, the Lindblad resonance is the only effective mechanism producing a significant density fluctuation}. In a last step, we replaced the Newtonian potential by the so called logarithmically modified pseudo-Newtonian potential in order to take into account some general-relativistic effects like the innermost stable circular orbit (ISCO). The latter potential is better suited to describe the close vicinity of a neutron star or a black hole. However, from a qualitative point of view, the resonance conditions remain the same. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.