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.     

Force-free magnetic fields on curvilinear coordinate surfaces

To seek nonlinear solutions of force-free magnetic fields, some symmetries or approximations are usually invoked. We consider magnetic fields lying on coordinate surfaces of an orthogonal curvilinear coordinate system. We conclude that only fields on parallel planes or spherical shells can be expressed in the form provided by Low in 1980s. These force-free fields are stable against small perturbations with rigid boundaries. Fields on cylindrical shells are also considered.     

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.     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

Constraining the origin and evolution of confined turbidite systems: southern Cretan margin,Eastern Mediterranean Sea (34°30–36°N)

Bathymetric, 9.5-kHz long-range sidescan sonar (OKEAN), seismic reflection and sediment-core data are used in the analysis of two tectonic troughs south of Crete, Eastern Mediterranean Sea. Here, up to 1.2 s two-way travel time (TWTT) of strata have accumulated since the Middle Miocene in association with extension in the South Aegean region. The study area comprises >100-km- long by >25-km-wide basins filled by sediments subdivided into two seismic units: (1) an upper Unit 1 deposited in sub-basins which follow the present-day configuration of the southern Cretan margin; (2) a basal Unit 2, more than 500 ms (TWTT) thick, accumulated in deeper half-graben/grabens distinct from the present-day depocentres. Both units overlap a locally stratified Unit 3 comprising the pre-Neogene core complex of Crete and Gavdos. In this work, the interpreted seismic units are correlated with the onshore stratigraphy, demonstrating that denudation processes occurring on Crete and Gavdos in response to major tectonic events have been responsible for high sedimentation rates along the proximal southern Cretan margin. Consequently, topographically confined sedimentary units have been deposited south of Crete in the last 12 Ma, including turbidites and other mass-flow deposits fed by evolving transverse and axial channel systems. Surface processes controlling facies distribution include the direct inflow of sediment from alluvial-fan systems and incising mountain rivers onto the Cretan slope, where significant sediment instability processes occur at present. In this setting, seismic profiles reveal eight different types of stratigraphic contacts on basin-margin highs, and basinal areas show evidence of halokinesis and/or fluid escape. The acquired data also show that significant changes to the margin’s configuration occurred in association with the post-Alpine tectonic and eustatic episodes affecting the Eastern Mediterranean.     

Geomorphological evidence for upslope canyon-forming processes on the northern KwaZulu-Natal shelf,SW Indian Ocean,South Africa

A geomorphological and statistical analysis of slope canyons from the northern KwaZulu-Natal continental margin is documented and compared with submarine canyons from the Atlantic margin of the USA. The northern KwaZulu-Natal margin is characterized by increasing upslope relief, concave slope-gradient profiles and features related to upslope growth of the canyon forms. Discounting slope-gradient profile, this morphology is strikingly similar to canyon systems of the New Jersey slope. Several phases of canyon incision indicate that downslope erosion is also an important factor in the evolution of the northern KwaZulu-Natal canyon systems. Despite the strong similarities between the northern KwaZulu-Natal and New Jersey slope-canyon systems, key differences are evident: (1) the concavity of the northern KwaZulu-Natal slope, contrasting with the ∼linear New Jersey slope; (2) the relative isolation of the northern KwaZulu-Natal canyons, rather than the dense clustering of the New Jersey canyons; and (3) the absence of strongly shelf-breaching canyons along the northern KwaZulu-Natal margin. In comparison with the New Jersey margin, we surmise a more youthful stage of canyon evolution, a result of either the canyons themselves being younger or the formative processes being less active. Less complicated patterns of erosion resulting from reduced sediment availability have developed in northern KwaZulu-Natal. The reduction in slope concavity on the New Jersey margin may be the result of grading of the upper slope by intensive headward erosion, a process more subdued—or less evident—on the KwaZulu-Natal margin.     

青岛文昌鱼的生活史──年龄、生长和死亡研究

于１９８９年８月－１９９０年７月对青岛文昌鱼的生活史进行研究。结果表明：（１）青岛文昌鱼种群由４个齿组组成，它至少生产４ａ以上，甚至可达６ａ左右。（２）在Ⅰ、Ⅱ齿时生长较快，平均年生长１３－１５ｍｍ，Ⅲ龄后生长速度减慢，即使同一年龄组的文晶鱼，在不同季节其生长速度也不同。（３）在冬春采集的文昌鱼中Ⅲ龄组的峰值较小。（４）繁殖季节过后，老龄文昌鱼因衰老而死亡。     

任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。