Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Improvement of the Discrete Element Method applied to Shear of Granular Media
95
0
0.0
(
0
)
Added by
Andr\\'es Alfonso Pe\\~na Olarte
Publication date
2007
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
The shear viscosity in the dilute regime of a model for confined granular matter is studied by simulations and kinetic theory. The model consists on projecting into two dimensions the motion of vibrofluidized granular matter in shallow boxes by modifying the collision rule: besides the restitution coefficient that accounts for the energy dissipation, there is a separation velocity that is added in each collision in the normal direction. The two mechanisms balance on average, producing stationary homogeneous states. Molecular dynamics simulations show that in the steady state the distribution function departs from a Maxwellian, with cumulants that remain small in the whole range of inelasticities. The shear viscosity normalized with stationary temperature presents a clear dependence with the inelasticity, taking smaller values compared to the elastic case. A Boltzmann-like equation is built and analyzed using linear response theory. It is found that the predictions show an excellent agreement with the simulations when the correct stationary distribution is used but a Maxwellian approximation fails in predicting the inelasticity dependence of the viscosity. These results confirm that transport coefficients depend strongly on the mechanisms that drive them to stationary states.
We use the first Betti number of a complex to characterize the morphological structure of granular samples in mechanical equilibrium. We analyze two-dimensional granular packings after a tapping process by means of both simulations and experiments. States with equal packing fraction obtained with different tapping intensities are distinguished after the introduction of a filtration parameter which determines the particles (nodes in the network) that are joined by an edge. We first use numerical simulations to characterize the effect of the precision in the particles localization by artificially adding different levels of noise in this magnitude. The outcomes obtained for the simulations are then compared with the experimental results allowing a clear distinction of experimental packings that have the same density. This is accomplished by just using the position of the particles and no other information about the possible contacts, or magnitude of forces.
A nonlinear analysis of high-frequency thickness-shear vibrations of AT-cut quartz crystal plates is presented with the two-dimensional finite element method. We expanded both kinematic and constitutive nonlinear Mindlin plate equations and then truncated them to the first-order equations as an approximation, which is used later for the formulation of nonlinear finite element analysis with all zeroth- and first-order displacements and electric potentials. The matrix equation of motion is established with the first-order harmonic approximation and the generalized nonlinear eigensystem is solved by a direct iterative procedure. A backbone curve and corresponding mode shapes are obtained and analyzed. The nonlinear finite element program is developed based on earlier linear edition and can be utilized to predict nonlinear characteristics of miniaturized quartz crystal resonators in the design process.
We apply a recently developed quasiparticle self-consistent $GW$ method (QSGW) to Gd, Er, EuN, GdN, ErAs, YbN and GdAs. We show that QSGW combines advantages separately found in conventional $GW$ and LDA+$U$ theory, in a simple and fully emph{ab initio} way. qsgw reproduces the experimental occupied $4f$ levels well, though unoccupied levels are systematically overestimated. Properties of the Fermi surface responsible for electronic properties are in good agreement with available experimental data. GdN is predicted to be very near a critical point of a first-order metal-insulator transition.
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 Omoris 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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
