Researcher profile

Sean McNamara

Sean McNamara contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2016arXiv

Eshelby inclusions in granular matter: theory and simulations

We present a numerical implementation of an active inclusion in a granular material submitted to a biaxial test. We discuss the dependence of the response to this perturbation on two parameters: the intra-granular friction coefficient on one hand, the degree of the loading on the other hand. We compare the numerical results to theoretical predictions taking into account the change of volume of the inclusion as well as the anisotropy of the elastic matrix.

preprint2012arXiv

On the relevance of numerical simulations to booming sand

We have performed a simulation study of 3D cohesionless granular flows down an inclined chute. We find that the oscillations observed in [L.E. Silbert, Phys. Rev. Lett., 94, 098002 (2005)] near the angle of repose are harmonic vibrations of the lowest normal mode. Their frequencies depend on the contact stiffness as well as on the depth of the flow. Could these oscillations account for the phenomena of "booming sand"? We estimate an effective contact stiffness from the Hertz law, but this leads to frequencies several times higher than observed. However, the Hertz law also predicts interpenetrations of a few nanometers, indicating that the oscillations frequencies are governed by the surface stiffness, which can be much lower than the bulk one. This is in agreement with previous studies ascribing the ability to sing to the presence of a soft coating on the grain surface.

preprint2010arXiv

Aspects of weakening in granular assemblies

We study numerical simulations of large (N~10^4) two-dimensional quasi-static granular assemblies subjected to a slowly increasing deviator stress. We report some peculiarities in the behavior of these packings that have not yet been adressed. The number of sliding contacts is not necessarily related to stability: first the number of sliding contacts rises linearly and smoothly with the applied stress. Then, at approximately half the peak stress, the increase slows down, a plateau develops, and a decrease follows.

preprint2008arXiv

Avalanches in anisotropic sheared granular media

We study the influence of particle shape anisotropy on the occurrence of avalanches in sheared granular media. We use molecular dynamic simulations to calculate the relative movement of two tectonic plates. % with transform boundaries. Our model considers irregular polygonal particles constituting the material within the shear zone. We find that the magnitude of the avalanches is approximately independent on particle shape and in good agreement with the Gutenberg-Richter law, but the aftershock sequences are strongly influenced by the particle anisotropy yielding variations on the exponent characterizing the empirical Omori's law. Our findings enable one to identify the presence of anisotropic particles at the macro-mechanical level only by observing the avalanche sequences of real faults. In addition, we calculate the probability of occurrence of an avalanche for given values of stiffness or frictional strength and observe also a significant influence of the particle anisotropy.

preprint2007arXiv

Numerical Improvement of the Discrete Element Method applied to Shear of Granular Media

We present a detailed analysis of the bounds on the integration step in Discrete Element Method (DEM) for simulating collisions and shearing of granular assemblies. We show that, in the numerical scheme, the upper limit for the integration step, usually taken from the average time $t_c$ of one contact, is in fact not sufficiently small to guarantee numerical convergence of the system during relaxation. In particular, we study in detail how the kinetic energy decays during the relaxation stage and compute the correct upper limits for the integration step, which are significantly smaller than the ones commonly used. In addition, we introduce an alternative approach, based on simple relations to compute the frictional forces, that converges even for integration steps above the upper limit.

preprint2006arXiv

Computer Simulation of Particle Suspensions

Particle suspensions are ubiquitous in our daily life, but are not well understood due to their complexity. During the last twenty years, various simulation methods have been developed in order to model these systems. Due to varying properties of the solved particles and the solvents, one has to choose the simulation method properly in order to use the available compute resources most effectively with resolving the system as well as needed. Various techniques for the simulation of particle suspensions have been implemented at the Institute for Computational Physics allowing us to study the properties of clay-like systems, where Brownian motion is important, more macroscopic particles like glass spheres or fibers solved in liquids, or even the pneumatic transport of powders in pipes. In this paper we will present the various methods we applied and developed and discuss their individual advantages.

preprint2004arXiv

Fragmentation of a Circular Disc by Impact on a Frictionless Plate

The break-up of a two-dimensional circular disc by normal and oblique impact on a hard frictionless plate is investigated by molecular dynamics simulations. The disc is composed of numerous unbreakable randomly shaped convex polygons connected together by simple elastic beams that break when bent or stretched beyond a certain limit. It is found that for both normal and oblique impacts the crack patterns are the same and depend solely on the normal component of the impact velocity. Analysing the pattern of breakage, amount of damage, fragment masses and velocities, we show the existence of a critical velocity which separates two regimes of the impact process: below the critical point only a damage cone is formed at the impact site (damage), cleaving of the particle occurs at the critical point, while above the critical velocity the disc breaks into several pieces (fragmentation). In the limit of very high impact velocities the disc suffers complete disintegration (shattering) into many small fragments. In agreement with experimental results, fragment masses are found to follow the Gates-Gaudin-Schuhmann distribution (power law) with an exponent independent of the velocity and angle of impact. The velocity distribution of fragments exhibit an interesting anomalous scaling behavior when changing the impact velocity and the size of the disc.

preprint2003arXiv

Vibrated granular media as experimentally realizable Granular Gases

We report numerical simulations of strongly vibrated granular materials designed to mimic recent experiments performed both in presence [1] or absence [2] of gravity. We show that a model with impact velocity dependent restitution coefficient is necessary to bring the simulations into agreement with experiments. We measure the scaling exponents of the granular temperature, collision frequency, impulse and pressure with the vibrating piston velocity. As the system changes from a homogeneous gas state at low density to a clustered state at high density, these exponents are all found to decrease continuously with the particle number. In absence of gravity, a loose cluster appears near the upper wall, opposite the piston, and acts as a buffer for fastest particles leading to unexpected non-extensive scaling exponents ; whereas in presence of gravity, the cluster bounces as a single inelastic body. All these results differ significantly from classical inelastic hard sphere kinetic theory and previous simulations, both based on a constant restitution coefficient.