Numerical simulations of granular dynamics: I. Hard-sphere discrete element method and tests

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular “atmosphere” that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.     

A numerical model of cohesion in planetary rings

We present a numerical method that incorporates particle sticking in simulations using the N-body code pkdgrav to study motions in a local rotating frame, such as a patch of a planetary ring. Particles stick to form non-deformable but breakable aggregates that obey the (Eulerian) equations of rigid-body motion. Applications include local simulations of planetary ring dynamics and planet formation, which typically feature hundreds of thousands or more colliding bodies. Bonding and breaking thresholds are tunable parameters that can approximately mimic, for example, van der Waals forces or interlocking of surface frost layers. The bonding and breaking model does not incorporate a rigorous treatment of internal fracture; rather the method serves as motivation for first-order investigation of how semi-rigid bonding affects the evolution of particle assemblies in high-density environments.We apply the method to Saturn’s A ring, for which laboratory experiments suggest that interpenetration of thin, frost-coated surface layers may lead to weak cohesive bonding. These experiments show that frost-coated icy bodies can bond at the low impact speeds characteristic of the rings. Our investigation is further motivated by recent simulations that suggest a very low coefficient of restitution is needed to explain the amplitude of the azimuthal brightness asymmetry in Saturn’s A ring, and the hypothesis that fine structure in Saturn’s B ring may in part be caused by large-scale cohesion.This work presents the full implementation of our model in pkdgrav, as well as results from initial tests with a limited set of parameters explored. We find a combination of parameters that yields aggregate size distribution and maximum radius values in agreement with Voyager data for ring particles in Saturn’s outer A ring. We also find that the bonding and breaking parameters define two strength regimes in which fragmentation is dominated either by collisions or other stresses, such as tides. We conclude our study with a discussion of future applications of and refinements to our model.     

Moisture Movement Through Cracked Clay Soil Profiles

A continuum mechanics approach is used for the formulation of unsaturated hydraulic conductivity functions and the water storage functions for fractured or cracked clay soils in this parametric study. Suggested procedures are based on available research literature related to the behavior of cracked unsaturated porous media. The soil–water characteristic curve, hydraulic conductivity and water storage functions take on the character of bi-modal unsaturated soil property functions. The bimodal character arises out of the independent behavior of the cracks and the intact clay soil. Matric suction changes beneath a slab-on-grade foundation placed on a cracked clay soil profile are modeled for varied surface flux conditions using the proposed unsaturated hydraulic conductivity and water storage functions. The impact of various levels of surface cracking on soil suction distributions is discussed. Suction distribution patterns are dependent on the initial soil surface suction. In particular, the results are dependent upon whether the initial matric suction is less than or greater than the air entry of the cracked clay. There is an extremely wide range of possible conditions that could be modeled but the parametric study results presented in this paper are limited to a series of selected crack widths and densities for an exfiltration case and an infiltration case.     

The Big Circles in Jordan: First absolute ages using rock luminescence surface dating

In this study, we provide the first absolute ages for a Big Circle megalithic structure in Jordan, using rock surface luminescence dating of the buried surface of rocks collected from circle J4 in southern Jordan. Five rocks were used for this study. All rocks showed evidence of previous daylight exposure before being used in the construction of the circle. The exposure was sufficient to bleach the latent luminescence signal to a negligible level compared to the subsequent burial dose. Three rocks gave indistinguishable ages, and were last exposed to daylight in 1500 ± 100 BCE; this is very likely to be the date of circle construction. Two others gave younger ages, indicating later disturbance or reworking. These new results provide very strong evidence for construction during the Late Bronze Age, and refute the earlier hypothesis of construction during Umayyad period (661–750 CE) as a hunting trap.     

Brackish water faunas from the St Maughans Formation: The Old Red Sandstone section at Ammons Hill,Hereford and Worcester,UK, re-examined

The disused railway cutting at Ammons Hill, Hereford and Worcester, exposes a sequence of beds belonging to the Devonian St Maughans Formation of Lochkovian (Gedinnian) age. The beds are of Old Red Sandstone facies, but contain brackish water faunas. These faunas occur at a level generally considered to be above the level of marine influence that affected the older Raglan Mudstone Formation of mainly Přídolí Series age. The section, described by King in 1934, is now overgrown, but was excavated in 1986 by the British Geological Survey during its survey of the Worcester 1:50000 sheet. The evidence of the section calls for slight amendment of Allen's (1985) model of an interrupted transition from marine deposition in Ludlow time to freshwater deposition in Gedinnian time that was complete by the time of the formation of a regionally extensive calcrete palaeosol, the Psammosteus Limestone. Subsequent transgressive events took place before the establishment of apparently wholly fluvial and floodplain environments.