Low-frequency oscillations in narrow vibrated granular systems

We present simulations and a theoretical treatment of vertically vibrated granular media. The systems considered are confined in narrow quasi-two-dimensional and quasi-one-dimensional (column) geometries, where the vertical extension of the container is much larger than both horizontal lengths. The additional geometric constraint present in the column setup frustrates the convection state that is normally observed in wider geometries. This makes it possible to study collective oscillations of the grains with a characteristic frequency that is much lower than the frequency of energy injection. The frequency and amplitude of these oscillations are studied as a function of the energy input parameters and the size of the container. We observe that, in the quasi-two-dimensional setup, low-frequency oscillations are present even in the convective regime. This suggests that they may play a significant role in the transition from a density inverted state to convection. Two models are also presented; the first one, based on Cauchy's equations, is able to predict with high accuracy the frequency of the particles' collective motion. This first principles model requires a single input parameter, i.e. the centre of mass of the system. The model shows that a sufficient condition for the existence of the low-frequency mode is an inverted density profile with distinct low and high density regions, a condition that may apply to other systems too. The second, simpler model just assumes an harmonic oscillator like behaviour and, using thermodynamic arguments, is also able to reproduce the observed frequencies with high accuracy

Mercury-DPM: Fast particle simulations in complex geometries

Mercury-DPM is a code for performing discrete particle simulations. That is to say, it simulates the motion of particles, or atoms, by applying forces and torques that stem either from external body forces, (e.g. gravity, magnetic fields, etc…) or from particle interactions. For granular particles, these are typically contact forces (elastic, viscous, frictional, plastic, cohesive), while for molecular simulations, forces typically stem from interaction potentials (e.g. Lennard-Jones). Often the method used in these packages is referred to as the discrete element method (DEM), which was originally designed for geotechnical applications. However, as Mercury-DPM is designed for simulating particles with emphasis on contact models, optimized contact detection for highly different particle sizes, and in-code coarse graining (in contrast to post-processing), we prefer the more general name discrete particle simulation. The code was originally developed for granular chute flows, and has since been extended to many other granular applications, including the geophysical modeling of cinder cone creation. Despite its granular heritage it is designed in a flexible way so it can be adapted to include other features such as long-range interactions and non-spherical particles, etc

Solar Magnetic Carpet I: Simulation of Synthetic Magnetograms

This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve through the processes of emergence, cancellation, coalescence and fragmentation. Model parameters for the emergence of bipoles are based upon the results of observational studies. Using this model, several simulations are considered, where the range of flux with which bipoles may emerge is varied. In all cases the model quickly reaches a steady state where the rates of emergence and cancellation balance. Analysis of the resulting magnetic field shows that we reproduce observed quantities such as the flux distribution, mean field, cancellation rates, photospheric recycle time and a magnetic network. As expected, the simulation matches observations more closely when a larger, and consequently more realistic, range of emerging flux values is allowed (4e16 - 1e19 Mx). The model best reproduces the current observed properties of the magnetic carpet when we take the minimum absolute flux for emerging bipoles to be 4e16 Mx. In future, this 2D model will be used as an evolving photospheric boundary condition for 3D non-potential modeling.Comment: 33 pages, 16 figures, 5 gif movies included: movies may be viewed at http://www-solar.mcs.st-and.ac.uk/~karen/movies_paper1

Discrete particle simulation of the spreading process in additive manufacturing

Selective Laser Sintering/Melting (SLS/SLM) are additive manufacturing (AM) technologies. Objects are produced by spreading successive layers of powder material and solidifying selected parts by sintering/melting them with a laser. The focus of this study is the powder spreading process for which the powder characteristics plays a major role for the powder layer quality, that in turn, influences the final product properties. The spreading process of a characteristic, frequently used, Ti-6Al-4V powder is simulated in MercuryDPM, using a discrete particle model. A parameter study varying cohesion, sliding and rolling friction allows us to quantify the influence of these powder properties on the layer characteristics, such as density and uniformity. The layer characteristics were obtained by coarse-graining, which generates grid-free continuum fields, e.g., density from discrete data. The density and homogeneity of the powder layer decreased with the increase of interparticle friction, leading to non-uniform layer, higher porosity, and dragged particles causing defects in the powder bed. However, the larger interparticle friction led to a rather good bed. In addition, the sliding friction had a little effect on the layer uniformity, but a large effect on particle segregation, whereas the rolling friction had a larger effect on layer uniformity. Further investigations will focus on additional parametric studies, experimental validation, the effect of humidity and spreading tool design evaluation

Modeling the break-up of nano-particle clusters in aluminum- and magnesium-based metal matrix nano-composites

Aluminum- and magnesium-based metal matrix nano-composites with ceramic nano-reinforcements promise low weight with high durability and superior strength, desirable properties in aerospace, automobile, and other applications. However, nano-particle agglomerations lead to adverse effects on final properties: large-size clusters no longer act as dislocation anchors, but instead become defects; the resulting particle distribution will be uneven, leading to inconsistent properties. To prevent agglomeration and to break-up clusters, ultrasonic processing is used via an immersed sonotrode, or alternatively via electromagnetic vibration. A study of the interaction forces holding the nano-particles together shows that the choice of adhesion model significantly affects estimates of break-up force and that simple Stokes drag due to stirring is insufficient to break-up the clusters. The complex interaction of flow and co-joint particles under a high frequency external field (ultrasonic, electromagnetic) is addressed in detail using a discrete-element method code to demonstrate the effect of these fields on de-agglomeration

Shape matters: Competing mechanisms of particle shape segregation

It is well-known that granular mixtures that differ in size or shape segregate when sheared. In the past, two mechanisms have been proposed to describe this effect, and it is unclear if both exist. To settle this question, we consider a bidisperse mixture of spheroids of equal volume in a rotating drum, where the two mechanisms are predicted to act in opposite directions. We present the first evidence that there are two distinct segregation mechanisms driven by relative over-stress. Additionally, we showed that for non-spherical particles, these two mechanisms can act in different directions leading to a competition between the effects of the two. As a result, the segregation intensity varies nonmonotonically as a function of AR, and at specific points, the segregation direction changes for both prolate and oblate spheroids, explaining the surprising segregation reversal previously reported. Consistent with previous results, we found that the kinetic mechanism is dominant for (almost) spherical particles. Furthermore, for moderate aspect ratios, the kinetic mechanism is responsible for the spherical particles segregation to the periphery of the drum, and the gravity mechanism plays only a minor role. Whereas, at the extreme values of AR, the gravity mechanism notably increases and overtakes its kinetic counterpart