BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefan Luding

Stefan Luding appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.flu-dynphysics.comp-phcond-mat.dis-nncond-mat.mtrl-scinlin.AOnlin.PSphysics.data-an

Catalog footprint

What is connected

30works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

30 published item(s)

preprint2020arXiv

Stretching the limits of dynamic and quasi-static flow testing on limestone powders

Powders are a special class of granular matter due to the important role of cohesive forces. The flow behavior of powders depends on the flow states and stress and is therefore difficult to measure/quantify with only one experiment. In this study, the most commonly used characterization tests that cover a wide range of states are compared: (static, free surface) angle of repose, the (quasi-static, confined) ring shear steady state angle of internal friction, and the (dynamic, free surface) rotating drum flow angle are considered for free flowing, moderately and strongly cohesive limestone powders. The free flowing powder gives good agreement among all different situations (devices), while the moderately and strongly cohesive powders behave more interestingly. Starting from the flow angle in the rotating drum and going slower, one can extrapolate to the limit of zero rotation rate, but then observes that the angle of repose measured from the heap is considerably larger, possibly due to its special history. When we stretch the ring shear test to its lowest confining stress limit, the steady state angle of internal friction of the cohesive powder coincides with the flow angle (at free surface) in the zero rotation rate limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The influence of material properties and process parameters on the spreading process in additive manufacturing

Laser powder bed fusion (LPBF) is an additive manufacturing (AM) technology. To achieve high product quality, the powder is best spread as a uniform, dense layer. The challenge for LPBF manufacturers is to develop a spreading process that can produce a consistent layer quality for the many powders used, which show considerable differences in spreadability. Therefore, we investigate the influence of material properties, process parameters and the type of spreading tool on the layer quality. The discrete particle method is used to simulate the spreading process and to define metrics to evaluate the powder layer characteristics. We found that particle shape and surface roughness in terms of rolling resistance and interparticle sliding friction as well as particle cohesion all have a major (sometimes surprising) influence on the powder layer quality: more irregular shaped particles, rougher particle surfaces and/or higher interfacial cohesion usually, but not always, lead to worse spreadability. Our findings illustrate that there is a trade-off between material properties and process parameters. Increasing the spreading speed decreases layer quality for non- and weakly cohesive powders, but improves it for strongly cohesive ones. Using a counter-clockwise rotating roller as a spreading tool improves the powder layer quality compared to spreading with a blade. Finally, for both geometries, a unique correlation between the quality criteria uniformity and mass fraction is reported and some of the findings are related to size-segregation during spreading.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Memory of jamming - multiscale models for soft and granular matter

Soft, disordered, micro-structured materials are ubiquitous in nature and industry, and are different from ordinary fluids or solids, with unusual, interesting static and flow properties. The transition from fluid to solid -at the so-called jamming density- features a multitude of complex mechanisms, but there is no unified theoretical framework that explains them all. In this study, a simple yet quantitative and predictive model is presented, which allows for a variable, changing jamming density, encompassing the memory of the deformation history and explaining a multitude of phenomena at and around jamming. The jamming density, now introduced as a new state-variable, changes due to the deformation history and relates the system's macroscopic response to its microstructure. The packing efficiency can increase logarithmically slow under gentle repeated (isotropic) compression, leading to an increase of the jamming density. In contrast, shear deformations cause anisotropy, changing the packing efficiency exponentially fast with either dilatancy or compactancy. The memory of the system near jamming can be explained by a microstatistical model that involves a multiscale, fractal energy landscape and links the microscopic particle picture to the macroscopic continuum description, providing a unified explanation for the qualitatively different flow-behavior for different deformation modes. To complement our work, a recipe to extract the history-dependent jamming density from experimentally accessible data is proposed, and alternative state-variables are compared. The proposed simple macroscopic constitutive model is calibrated with the memory of microstructure. Such approach can help understanding predicting and mitigating failure of structures or geophysical hazards, and will bring forward industrial process design/optimization, and help solving scientific challenges in fundamental research.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Towards a general(ized) shear thickening rheology of wet granular materials under small pressure

We study the rheology of dry and wet granular materials in the steady quasistatic regime using the Discrete Element Method (DEM) in a split-bottom ring shear cell with focus on the macroscopic friction. The aim of our study is to understand the local rheology of bulk flow at various positions in the shear band, where the system is in critical state. The general(ized) rheology has four dimensionless control parameters that relate the time scales of five significant phenomena, namely, the time scales related to confining pressure $t_p$, shear rate $t_{\dotγ}$, particle stiffness $t_k$, gravity $t_g$ and cohesion $t_c$, respectively. We show that those phenomena collectively contribute to the rheology as multiplicative correction functions. While $t_{\dotγ}$ is large and thus little important for most of the data studied, it can increase the friction of flow in critical state, where the shear gradients are high. $t_g$ and $t_k$ are comparable to $t_p$ in the bulk, but become more or less dominant relative to $t_p$ at the extremes of the free surface and deep inside the bulk, respectively. We also measure the effect of strong wet cohesion on the flow rheology, as quantified by decreasing $t_c$. Furthermore, the proposed rheological model predicts well the shear thinning behavior both in the bulk and near the free surface; shear thinning develops to shear thickening near the free surface with increasing cohesion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tuning the bulk properties of bidisperse granular mixtures by small amount of fines

We study the bulk properties of isotropic bidisperse granular mixtures using discrete element simulations. The focus is on the influence of the size (radius) ratio of the two constituents and volume fraction on the mixture properties. We show that the effective bulk modulus of a dense granular (base) assembly can be enhanced by up to 20% by substituting as little as 5% of its volume with smaller sized particles. Particles of similar sizes barely affect the macroscopic properties of the mixture. On the other extreme, when a huge number of fine particles are included, most of them lie in the voids of the base material, acting as rattlers, leading to an overall weakening effect. In between the limits, an optimum size ratio that maximizes the bulk modulus of the mixture is found. For loose systems, the bulk modulus decreases monotonically with addition of fines regardless of the size ratio. Finally, we relate the mixture properties to the 'typical' pore size in a disordered structure as induced by the combined effect of operating volume fraction (consolidation) and size ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A contact model for sticking of adhesive meso-particles

The interaction between visco-elasto-plastic and adhesive particles is the subject of this study, where "meso-particles" are introduced, i.e., simplified particles, whose contact mechanics is not taken into account in all details. A few examples of meso-particles include agglomerates or groups of primary particles, or inhomogeneous particles with micro-structures of the scale of the contact deformation, such as core-shell materials. A simple, flexible contact model for meso-particles is proposed, which allows to model the bulk behavior of assemblies of many particles in both rapid and slow, quasi-static flow situations. An attempt is made to categorize existing contact models for the normal force, discuss all the essential mechanical ingredients that must enter the model (qualitatively) and finally solve it analytically. The model combines a short-ranged, non-contact part (resembling either dry or wet materials) with an elaborate, visco-elasto-plastic and adhesive contact law. Using energy conservation arguments, an analytical expression for the coefficient of restitution is derived in terms of the impact velocity (for pair interactions or, equivalently, without loss of generality, for quasi-static situations in terms of the maximum overlap or confining stress). Adhesive particles (or meso-particles) stick to each other at very low impact velocity, while they rebound less dissipatively with increasing velocity, in agreement with previous studies. For even higher impact velocities an interesting second sticking and rebound regime is reported. The low velocity sticking is due to non-contact adhesive forces, the first rebound regime is due to stronger elastic and kinetic energies with little dissipation, while the high velocity sticking is generated by the non-linearly increasing, history dependent plastic dissipation and adhesive contact force.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Master equation for force distributions in polydisperse frictional particles

An incremental evolution equation, i.e. a Master equation in statistical mechanics, is introduced for force distributions in polydisperse frictional particle packings. As basic ingredients of the Master equation, the conditional probability distributions of particle overlaps are determined by molecular dynamics simulations. Interestingly, tails of the distributions become much narrower in the case of frictional particles than frictionless particles, implying that correlations of overlaps are strongly reduced by microscopic friction. Comparing different size distributions, we find that the tails are wider for the wider distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A master equation for force distributions in soft particle packings - Irreversible mechanical responses to isotropic compression and decompression

Mechanical responses of soft particle packings to quasi-static deformations are determined by the microscopic restructuring of force-chain networks, where complex non-affine displacements of constituent particles cause the irreversible macroscopic behavior. Recently, we have proposed a master equation for the probability distribution functions of contact forces and interparticle gaps [K. Saitoh et al., Soft Matter 11, 1253 (2015)], where mutual exchanges of contacts and interparticle gaps, i.e. opening and closing contacts, are also involved in the stochastic description with the aid of Delaunay triangulations. We describe full details of the master equation and numerically investigate irreversible mechanical responses of soft particle packings to cyclic loading. The irreversibility observed in molecular dynamics simulations is well reproduced by the master equation if the system undergoes quasi-static deformations. We also confirm that the degree of irreversible responses is a decreasing function of the area fraction and the number of cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From the granular Leidenfrost state to buoyancy-driven convection

Grains inside a vertically vibrated box undergo a transition from a density inverted and horizontally homogeneous state, referred to as the granular Leidenfrost state, to a buoyancy-driven convective state. We perform a simulational study of the precursors of such a transition, and quantify their dynamics as the bed of grains is progressively fluidized. The transition is preceded by transient convective states, which increase their correlation time as the transition point is approached. Increasingly correlated convective flows lead to density fluctuations, as quantified by the structure factor, that also shows critical behaviour near the transition point. The amplitude of the modulations in the vertical velocity field are seen to be best described by a quintic supercritical amplitude equation with an additive noise term. The validity of such an amplitude equation, and previously observed collective semi-periodic oscillations of the bed of grains, suggests a new interpretation of the transition analogous to a coupled chain of vertically vibrated damped oscillators. Increasing the size of the container shows metastability of convective states, as well as an overall invariant critical behaviour close to the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Macroscopic model with anisotropy based on micro-macro informations

Physical experiments can characterize the elastic response of granular materials in terms of macroscopic state-variables, namely volume (packing) fraction and stress, while the microstructure is not accessible and thus neglected. Here, by means of numerical simulations, we analyze dense, frictionless, granular assemblies with the final goal to relate the elastic moduli to the fabric state, i.e., to micro-structural averaged contact network features as contact number density and anisotropy. The particle samples are first isotropically compressed and later quasi-statically sheared under constant volume (undrained conditions). From various static, relaxed configurations at different shear strains, now infinitesimal strain steps are applied to "measure" the effective elastic response; we quantify the strain needed so that plasticity in the sample develops as soon as contact and structure rearrangements happen. Because of the anisotropy induced by shear, volumetric and deviatoric stresses and strains are cross-coupled via a single anisotropy modulus, which is proportional to the product of deviatoric fabric and bulk modulus (i.e. the isotropic fabric). Interestingly, the shear modulus of the material depends also on the actual stress state, along with the contact configuration anisotropy. Finally, a constitutive model based on incremental evolution equations for stress and fabric is introduced. By using the previously measured dependence of the stiffness tensor (elastic moduli) on the microstructure, the theory is able to predict with good agreement the evolution of pressure, shear stress and deviatoric fabric (anisotropy) for an independent undrained cyclic shear test, including the response to reversal of strain.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Micro-macro transition and simplified contact models for wet granular materials

Wet granular materials in a quasi-static steady state shear flow have been studied with discrete particle simulations. Macroscopic quantities, consistent with the conservation laws of continuum theory, are obtained by time averaging and spatial coarse-graining. Initial studies involve understanding the effect of liquid content and liquid properties like the surface tension on the macroscopic quantities. Two parameters of the liquid bridge contact model have been studied as the constitutive parameters that define the structure of this model (i) the rupture distance of the liquid bridge model, which is proportional to the liquid content, and (ii) the maximum adhesive force, as controlled by the surface tension of the liquid. Subsequently a correlation is developed between these micro parameters and the steady state cohesion in the limit of zero confining pressure. Furthermore, as second result, the macroscopic torque measured at the walls, which is an experimentally accessible parameter, is predicted from our simulation results as a dependence on the micro-parameters. Finally, the steady state cohesion of a realistic non-linear liquid bridge contact model scales well with the steady state cohesion for a simpler linearized irreversible contact model with the same maximum adhesive force and equal energy dissipated per contact.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collective motion of macroscopic spheres floating on capillary ripples: Dynamic heterogeneity and dynamic criticality

When a dense monolayer of macroscopic slightly polydisperse spheres floats on chaotic capillary Faraday waves, a coexistence of large scale convective motion and caging dynamics typical for jammed systems is observed. We subtract the convective mean flow using a coarse graining and reveal subdiffusion for the caging time scales followed by a diffusive regime at later times. To test the system in the light of dynamic criticality, we apply the methods of dynamic heterogeneity to obtain the power-law divergent time and length scales as the floater concentration approaches the jamming point. We find that these are independent of the application of the coarse graining procedure. The critical exponents are consistent with those found in dense suspensions of colloids indicating universal stochastic dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effect of cohesion on shear banding in quasi-static granular material

It is widely recognized in particle technology that adhesive powders show a wide range of different bulk behavior due to the peculiarity of the particle interaction. We use Discrete Element simulations to investigate the effect of contact cohesion on the steady state of dense powders in a slowly sheared split-bottom Couette cell, which imposes a wide stable shear band. The intensity of cohesive forces can be quantified by the {\em granular Bond number} ($Bo$), namely the ratio between maximum attractive force and average force due to external compression. We find that the shear banding phenomenon is almost independent of cohesion for Bond numbers $Bo<1$, but for $Bo \ge 1$ cohesive forces start to play an important role, as both width and center position of the band increase for $Bo > 1$. Inside the shear band, the mean normal contact force is always independent of cohesion and depends only on the confining stress. In contrast, when the behavior is analyzed focusing on the eigen-directions of the local strain rate tensor, a dependence on cohesion shows up. Forces carried by contacts along the compressive and tensile directions are symmetric about the mean force (larger and smaller respectively), while the force along the third, neutral direction follows the mean force. This anisotropy of the force network increases with cohesion, just like the heterogeneity in all (compressive, tensile and neutral) directions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Experiments and Discrete Element Simulation of the Dosing of Cohesive Powders in a Simplified Geometry

We perform experiments and discrete element simulations on the dosing of cohesive granular materials in a simplified geometry. The setup is a simplified canister box where the powder is dosed out of the box through the action of a constant-pitch screw feeder connected to a motor. A dose consists of a rotation step followed by a period of rest before the next dosage. From the experiments, we report on the operational performance of the dosing process through a variation of dosage time, coil pitch and initial powder mass. We find that the dosed mass shows an increasing linear dependence on the dosage time and rotation speed. In contrast, the mass output from the canister is not directly proportional to an increase/decrease in the number coils. By calibrating the interparticle friction and cohesion, we show that DEM simulation can quantitatively reproduce the experimental findings for smaller masses but also overestimate arching and blockage. With appropriate homogenization tools, further insights into microstructure and macroscopic fields can be obtained. This work shows that particle scaling and the adaptation of particle properties is a viable approach to overcome the untreatable number of particles inherent in experiments with fine, cohesive powders and opens the gateway to simulating their flow in more complex geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recent advances in the simulation of particle-laden flows

A substantial number of algorithms exists for the simulation of moving particles suspended in fluids. However, finding the best method to address a particular physical problem is often highly non-trivial and depends on the properties of the particles and the involved fluid(s) together. In this report we provide a short overview on a number of existing simulation methods and provide two state of the art examples in more detail. In both cases, the particles are described using a Discrete Element Method (DEM). The DEM solver is usually coupled to a fluid-solver, which can be classified as grid-based or mesh-free (one example for each is given). Fluid solvers feature different resolutions relative to the particle size and separation. First, a multicomponent lattice Boltzmann algorithm (mesh-based and with rather fine resolution) is presented to study the behavior of particle stabilized fluid interfaces and second, a Smoothed Particle Hydrodynamics implementation (mesh-free, meso-scale resolution, similar to the particle size) is introduced to highlight a new player in the field, which is expected to be particularly suited for flows including free surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rheology of weakly wetted granular materials - a comparison of experimental and numerical data

Shear cell simulations and experiments of weakly wetted particles (a few volume percent liquid binders) are compared, with the goal to understand their flow rheology. Application examples are cores for metal casting by core shooting made of sand and liquid binding materials. The experiments are carried out with a Couette-like rotating viscometer. The weakly wetted granular materials are made of quartz sand and small amounts of Newtonian liquids. For comparison, experiments on dry sand are also performed with a modified configuration of the viscometer. The numerical model involves spherical, monodisperse particles with contact forces and a simple liquid bridge model for individual capillary bridges between two particles. Different liquid content and properties lead to different flow rheology when measuring the shear stress-strain relations. In the experiments of the weakly wetted granular material, the apparent shear viscosity $η_g$ scales inversely proportional to the inertial number $I$, for all shear rates. On the contrary, in the dry case, an intermediate scaling regime inversely quadratic in $I$ is observed for moderate shear rates. In the simulations, both scaling regimes are found for dry and wet granular material as well.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Role of gravity or confining pressure and contact stiffness in granular rheology

The steady shear rheology of granular materials is investigated in slow quasi-static states and inertial flows. The effect of the gravity field and contact stiffness, which are conventionally trivialized is the focus of this study. Series of Discrete Element Method simulations are performed on a weakly frictional granular assembly in a split-bottom geometry considering various gravity fields and contact stiffnesses. While traditionally the inertial number, i.e., the ratio of stress to strain-rate timescales describes the flow rheology, we find that a second dimensionless number, the ratio of softness and stress timescales, must also be included to characterize the bulk flow behavior. For slow, quasi-static flows, the density increases while the macroscopic friction decreases with respective increase in particle softness and gravity. This trend is added to the $μ(I)$ rheology and can be traced back to the anisotropy in the contact network, displaying a linear correlation between macroscopic friction and deviatoric fabric in the steady state. Interestingly, the linear relation holds when the external rotation rate is increased for a given gravity field and contact stiffness.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rolling, sliding & torsion of micron-sized silica particles - Experimental, numerical and theoretical analysis

The contact mechanics of individual, very small particles with other particles and walls is studied using a nanoindenter setup that allows normal and lateral displacement control and measurement of the respective forces. The sliding, rolling and torsional forces and torques are tested with borosilicate microspheres, featuring radii of about 10$μ$m. The contacts are with flat silicon substrates of different roughness for pure sliding and rolling and with silicon based, ion-beam crafted rail systems for combined rolling and torsion. The experimental results are discussed and compared to various analytical predictions and contact models, allowing for two concurrent interpretations of the effects of surface roughness, plasticity and adhesion. This enables us to determine both rolling and torsion friction coefficients together with their associated length scales. Interestingly, even though normal contacts behave elastically (Hertzian), all other modes of motion display effects due to surface roughness and consequent plastic deformation. The influence of adhesion is interpreted in the framework of different models and is very different for different degrees of freedom, being largest for rolling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Slow relaxation behavior of cohesive powders

We present findings from uniaxial (oedometric) compression tests on two cohesive industrially relevant granular materials (cocoa and limestone powder). Experimental results are presented for the compressibility, tested with two devices -- the FT4 Powder Rheometer and the custom made lambdameter. We focus on the stress response and the slow relaxation behavior of the cohesive samples tested. After compression ends, at constant volume, the ongoing stress relaxation is found to follow a power law consistently for both cohesive powders and for the different testing equipments. A simple (incremental algebraic evolution) model is proposed for the stress relaxation in cohesive powders, which includes a response timescale along with a second, dimensionless relaxation parameter. The reported observations are useful for both the improvement of discrete element simulations and constitutive macroscopic models for cohesive granular materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structure characterization of hard sphere packings in amorphous and crystalline states

The channel size distribution in hard sphere systems, based on the local neighbor correlation of four particle positions, is investigated for all volume fractions up to jamming. For each particle, all three particle combinations of neighbors define channels, which are relevant for the concept of caging. The analysis of the channel size distribution is shown to be very useful in distinguishing between gaseous, liquid, partially and fully crystallized, and glassy (random) jammed states. A common microstructural feature of four coplanar particles is observed in crystalline and glassy jammed states, suggesting the presence of "hidden" two-dimensional order in three-dimensional random close packings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fluid-particle flow and validation using two-way-coupled mesoscale SPH-DEM

First, a meshless simulation method is presented for multiphase fluid-particle flows with a two-way coupled Smoothed Particle Hydrodynamics (SPH) for the fluid and the Discrete Element Method (DEM) for the solid phase. The unresolved fluid model, based on the locally averaged Navier Stokes equations, is expected to be considerably faster than fully resolved models. Furthermore, in contrast to similar mesh-based Discrete Particle Methods (DPMs), our purely particle-based method enjoys the flexibility that comes from the lack of a prescribed mesh. It is suitable for problems such as free surface flow or flow around complex, moving and/or intermeshed geometries and is applicable to both dilute and dense particle flows. Second, a comprehensive validation procedure for fluid-particle simulations is presented and applied here to the SPH-DEM method, using simulations of single and multiple particle sedimentation in a 3D fluid column and comparison with analytical models. Millimetre-sized particles are used along with three different test fluids: air, water and a water-glycerol solution. The velocity evolution for a single particle compares well (less than 1% error) with the analytical solution as long as the fluid resolution is coarser than two times the particle diameter. Two more complex multiple particle sedimentation problems (sedimentation of a homogeneous porous block and an inhomogeneous Rayleigh Taylor instability) are also reproduced well for porosities 0.6 <= ε<= 1.0, although care should be taken in the presence of high porosity gradients. Overall the SPH-DEM method successfully reproduces quantitatively the expected behaviour in the test cases, and promises to be a flexible and accurate tool for other, realistic fluid-particle system simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From creeping to inertial flow in porous media: a lattice Boltzmann - Finite Element study

The lattice Boltzmann method has been successfully applied for the simulation of flow through porous media in the creeping regime. Its technical properties, namely discretization, straightforward implementation and parallelization, are responsible for its popularity. However, flow through porous media is not restricted to near zero Reynolds numbers since inertial effects play a role in numerous natural and industrial processes. In this paper we investigate the capability of the lattice Boltzmann method to correctly describe flow in porous media at moderate Reynolds numbers. The selection of the lattice resolution, the collision kernel and the boundary conditions becomes increasingly important and the challenge is to keep artifacts due to compressibility effects at a minimum. The lattice Boltzmann results show an accurate quantitative agreement with Finite Element Method results and evidence the capability of the method to reproduce Darcy's law at low Reynolds numbers and Forchheimer's law at high Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Influence of polydispersity on micromechanics of granular materials

We study the effect of polydispersity on the macroscopic physical properties of granular packings in two and three dimensions. A mean-field approach is developed to approximate the macroscale quantities as functions of the microscopic ones. We show that the trace of the fabric and stress tensors are proportional to the mean packing properties (e.g. packing fraction, average coordination number, and average normal force) and dimensionless correction factors, which depend only on the moments of the particle-size distribution. Similar results are obtained for the elements of the stiffness tensor of isotropic packings in the linear affine response regime. Our theoretical predictions are in good agreement with the simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Slow Dynamics Near Jamming

Static and dynamic properties of two-dimensional bidisperse dissipative particles are numerically studied near the jamming transition. We investigate the dependency of the critical scaling on the ratio of the different diameters and find a new scaling of the maximum overlap (not consistent with the scaling of the mean overlap). The ratio of kinetic and potential energies drastically slows down near the jamming transition, i.e. the relaxation time diverges at the jamming point.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Closure Relations for Shallow Granular Flows from Particle Simulations

The Discrete Particle Method (DPM) is used to model granular flows down an inclined chute. We observe three major regimes: static piles, steady uniform flows and accelerating flows. For flows over a smooth base, other (quasi-steady) regimes are observed where the flow is either highly energetic and strongly layered in depth for small inclinations, or non-uniform and oscillating for larger inclinations. For steady uniform flows, depth profiles of density, velocity and stress have been obtained using an improved coarse-graining method, which allows accurate statistics even at the base of the flow. A shallow-layer model for granular flows is completed with macro-scale closure relations obtained from micro-scale DPM simulations of steady flows. We thus obtain relations for the effective basal friction, shape factor, mean density, and the normal stress anisotropy as functions of layer thickness, flow velocity and basal roughness. For collisional flows, the functional dependencies are well determined and have been obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

From discrete particles to continuum fields near a boundary

An expression for the stress tensor near an external boundary of a discrete mechanical system is derived explicitly in terms of the constituents' degrees of freedom and interaction forces. Starting point is the exact and general coarse graining formulation presented by Goldhirsch in [I. Goldhirsch, Gran. Mat., 12(3):239-252, 2010], which is consistent with the continuum equations everywhere but does not account for boundaries. Our extension accounts for the boundary interaction forces in a self-consistent way and thus allows the construction of continuous stress fields that obey the macroscopic conservation laws even within one coarse-graining width of the boundary. The resolution and shape of the coarse-graining function used in the formulation can be chosen freely, such that both microscopic and macroscopic effects can be studied. The method does not require temporal averaging and thus can be used to investigate time-dependent flows as well as static and steady situations. Finally, the fore-mentioned continuous field can be used to define 'fuzzy' (highly rough) boundaries. Two discrete particle method simulations are presented in which the novel boundary treatment is exemplified, including a chute flow over a base with roughness greater than a particle diameter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Behavior of pressure and viscosity at high densities for two-dimensional hard and soft granular materials

The pressure and the viscosity in two-dimensional sheared granular assemblies are investigated numerically. The behavior of both pressure and viscosity is smoothly changing qualitatively when starting from a mono-disperse hard-disk system without dissipation and moving towards a system of (i) poly-disperse, (ii) soft particles with (iii) considerable dissipation. In the rigid, elastic limit of mono-disperse systems, the viscosity is approximately inverse proportional to the area fraction difference from $ϕ_η \simeq 0.7$, but the pressure is still finite at $ϕ_η$. In moderately soft, dissipative and poly-disperse systems, on the other hand, we confirm the recent theoretical prediction that both scaled pressure (divided by the kinetic temperature $T$) and scaled viscosity (divided by $\sqrt{T}$) diverge at the same density, i.e., the jamming transition point $ϕ_J > ϕ_η$, with the exponents -2 and -3, respectively. Furthermore, we observe that the critical region of the jamming transition becomes invisible as the restitution coefficient approaches unity, i.e. for vanishing dissipation. In order to understand the conflict between these two different predictions on the divergence of the pressure and the viscosity, the transition from soft to hard particles is studied in detail and the dimensionless control parameters are defined as ratios of various time-scales. We introduce a dimensionless number, i.e. the ratio of dissipation rate and shear rate, that can identify the crossover from the scaling of very hard, i.e. rigid disks to the scaling in the soft, jamming regime.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Transport coefficients for dense hard-disk systems

A study of the transport coefficients of a system of elastic hard disks, based on the use of Helfand-Einstein expressions is reported. The self-diffusion, the viscosity, and the heat conductivity are examined with averaging techniques especially appropriate for the use in event-driven molecular dynamics algorithms with periodic boundary conditions. The density and size dependence of the results is analyzed, and comparison with the predictions from Enskog's theory is carried out. In particular, the behavior of the transport coefficients in the vicinity of the fluid-solid transition is investigated and a striking power law divergence of the viscosity in this region is obtained, while all other examined transport coefficients show a drop in that density range.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A driven two-dimensional granular gas with Coulomb friction

We study a homogeneously driven granular gas of inelastic hard particles with rough surfaces subject to Coulomb friction. The stationary state as well as the full dynamic evolution of the translational and rotational granular temperatures are investigated as a function of the three parameters of the friction model. Four levels of approximation to the (velocity-dependent) tangential restitution are introduced and used to calculate translational and rotational temperatures in a mean field theory. When comparing these theoretical results to numerical simulations of a randomly driven mono-layer of particles subject to Coulomb friction, we find that already the simplest model leads to qualitative agreement, but only the full Coulomb friction model is able to reproduce/predict the simulation results quantitatively for all magnitudes of friction. In addition, the theory predicts two relaxation times for the decay to the stationary state. One of them corresponds to the equilibration between the translational and rotational degrees of freedom. The other one, which is slower in most cases, is the inverse of the common relaxation rate of translational and rotational temperatures.

0 citations0 reviewsNo code linkedOpen access