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.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

On the microwave millisecond spike emission and its associated phenomena during the impulsive phase of large solar flare

We propose a model of two acceleration regions, which can explain, on the basis of microwave maser caused by a “hollow-beam” distribution of electrons, the presence of millisecond spikes in the event of 1981 May 16 and their absence in the event of 1981 October 12, and the enhanced continuous emission in the latter. We have also uncovered relations among the features, e.g. the Type III_G, Type IV_DCIM and hard x-ray bursts, that accompany the microwave millisecond spikes during the impulsive phase of a large flare.     

Simulating methane emission from a Chinese rice field as influenced by fertilizer and water level

Natural and agricultural wetlands are considered to be the major sources of global atmospheric methane (CH₄). A one‐dimensional model was developed to simulate methane emission and used to examine the influence of various physical processes on the rate of methane emission. Three processes involved in the methane emission are implemented in the model: production, reoxidation and transport. Three transport pathways were considered: diffusion across water–air or soil–air interfaces, ebullition and diffusion through plants. These pathways are influenced by soil properties, plant growth, water‐table conditions, temperature and external inputs (e.g. fertilizer). The model was used to examine the seasonal variation of the methane emission at a rice field in Hunan, China, which was observed during a field experiment for consecutive (early and late) rice seasons in 1992. The observed seasonal variations of methane emission, and role of plants in transporting methane to the atmosphere, are captured by the model simulation. Further model applications were conducted to simulate effects of fertilizer and water‐level condition on the methane emission. The results indicate that unfermented organic fertilizer produces a higher methane emission rate than mineral fertilizer. The simulations with treatments of a deep‐water covering and constant moisture reduced the methane emission. The rice field study provides a framework for further development of the model towards simulations based on spatially distributed variables (e.g. water table, soil temperature and vegetation) at a regional scale. Copyright © 2003 John Wiley & Sons, Ltd.     

Evolution of short period meteoroid streams

The dynamical evolution of meteoroid streams associated with cornets Encke, Halley, Machholz 1986 VIII and asteroid Phaethon is discussed. It is shown that the planetary perturbations can greatly increase the streams thickness and each stream may produce several couples of meteor showers active in different seasons of the year. The theoretical and observed data are in a satisfactory accordance.     

Spine–sheath polarization structures in four active galactic nuclei jets

We present the results of multifrequency (15 + 8 + 5 GHz) polarization Very Long Baseline Array (VLBA) observations of the three BL Lacertae objects 0745+241, 1418+546 and 1652+398 together with 5-GHz VLBI Space Observatory Programme (VSOP) observations of 1418+546 and 1.6- and 5-GHz VSOP observations of the blazar 1055+018. The jets of all these sources have polarization structure transverse to the jet axis, with the polarization E vectors aligned with the jet along the jet spine and 'sheaths' of orthogonal E vectors at one or both edges of the jet. The presence of polarization aligned with the jet near the 'spine' may indicate that the jets are associated with helical B fields that propogate outward with the jet flow; the presence of orthogonal polarization near the edges of the jet may likewise be a consequence of a helical jet B field, or may be owing to an interaction with the ambient medium on parsec scales. We have tentatively detected interknot polarization in 1055+018 with E aligned with the local jet direction, consistent with the possibility that the jet of this source is associated with a helical B field.     

A kinetic model of gas flow in a porous cometary mantle

A numerical study of gas flow through a porous cometary mantle is presented. A kinetic model based on the well-known Test Particle Monte Carlo Method for the solution of rarefied gas dynamics problems is proposed. The physical model consists of two spatial plane regions: the condensed ice phase and a porous dust mantle. The structure of the porous dust layer is described as a bundle of cylindrical inclined channels not crossing each other. A vertical temperature gradient may exist across the dust mantle. The aim is to investigate how the characteristics of molecular flow depend on the capillary length, inclination angle, and temperature gradient. Examples illustrating a significant deviation of some results from equilibrium values are shown. In particular, the gas velocity distribution at both ends of the pore is strongly non-Maxwellian if there is an important temperature contrast across the pore. The emergent gas flow rate is found to vary with the pore length/radius ratio in excellent agreement with Clausing's empirical formula. The degree of collimation of the flow is quantitatively studied as a function of the length/radius ratio, and consequences for the jet force of outgassing through a dust mantle or, indeed, a rough surface are estimated.     

The production of light elements in unconventional cosmologies with VMOs

The abundances of light elements (deuterium, helium and lithium) are calculated in unconventional cosmologies in which these elements are produced by interaction of fast -particles ejected by massive pregalactic stars with a purely hydrogenous cosmological substrate. The number density and energetics of such stars needed for suitable production of light elements are estimated. It is shown that the models always lead to overproduction of lithium. The intensity of ultraviolet, X-ray, and -ray emission caused by explosions of pregalactic massive stars is calculated. The distortions of cosmic background radiation by Compton losses of electrons heated by explosive stars are also calculated.     

Hydrodynamical models of type II supernovae

A series of hydrodynamical models of type-II supernova outbursts (SNII) has been calculated. Approximate relations connecting the total outburst energy ε, the mass of envelope ejectedM, the presupernova radiusR, and the amount of ionizing quanta radiated by the supernovaeN _H with such values as the duration of the light curve plateau Δt, and absolute magnitude in the wavelength bandV and photospheric velocityU _PH observed near the middle of the plateau have been established. Advantage has been taken of the relations to obtain a preliminary evaluation for the characteristics of the average SN II: ε=7×10⁵⁰ erg,M=6M _⊙,R=500R _⊙,N _H=2×10⁵⁸. The SNIIs with plateau-like light curves seem to be accounted for by thermonuclear explosions of degenerate cores of red giant stars and result in a total disruption of the star without any stellar remnant. To the contrary, SNIIs with linear light curves have substantially different properties (in particular, they throw considerably less massive envelopes off). These SNII must signify the birth of collapsed objects—neutron stars (pulsars) or black holes.