Quasi-static simulations of compact polydisperse particle systems

The Discrete Element Method (DEM) was originally devised by Cundall and Strack (1979), as a technique to examine the micromechanics of granular media with the anticipation that this would lead to more physically reliable continuum theories to describe the quasi-static deformation of granular material such as sand. However, the methodology models the evolution of a system of particles as a dynamic process. Consequently there have been numerous publications of the application of DEM to an increasingly wider v...     

Statistics of contact force network in dense granular matter

We investigate the distribution of contact forces in a static granular system and in annular shear granular flow, using the discrete element method, and considering the influences of both packing fraction and friction coefficient. We find the existence of a critical packing fraction. If the packing fraction is lower than this critical value, all contact forces in granular system vanish. For shear granular flow, the critical packing fraction is significantly smaller than that for static granular system. The ...     

Numerical solution to the shearing flow of granular materials between two plates

We study the shearing flow of granular materials between two horizontal flat plates where the top plate is moving with a constant speed. The constitutive relation used for the stress is based on the continuum model proposed by Rajagopal and Massoudi (DOE Report, DOE/PETC/TR-90/3, 1990). The material coefficients such as viscosity and normal stress coefficients are based on the model of Boyle and Massoudi (Int. J. Eng. Sci 28 (1990) 1261). The governing equations are non-dimensionalized and the resulting system of non-linear differential equations is solved numerically using finite difference technique.     

Particle dynamics in dense shear granular flow

The particle dynamics in an annular shear granular flow is studied using the discrete element method, and the influences of packing fraction, shear rate and friction coefficient are analyzed. We demonstrate the existence of a critical packing fraction exists in the shear granular flow. When the packing fraction is lower than this critical value, the mean tangential velocity profile exhibits a rate-independent feature. However, when the packing fraction exceeds this critical value, the tangential velocity profile becomes rate-dependent and varies gradually from linear to nonlinear with increasing shear rate. Furthermore, we find a continuous transition from the unjammed state to the jammed state in a shear granular flow as the packing fraction increases. In this transforming process, the force distribution varies distinctly and the contact force network also exhibits different features.     

Exact solution for entrance region flow between parallel plates

A method of obtaining an exact solution of the complete Navier-Stokes equations for incompressible flow has been developed to study the development of fluid flow between parallel plates. An analytical solution has been developed in the form of a convergent infinite series and a numerical solution obtained for Reynolds numbers of 10 and 500. These results have been compared with data available elsewhere. The method developed here is expected to have general application     

Boundary effects on behaviour of granular material during plane strain compression

The paper investigates the boundary effect on the behaviour of granular materials during plane strain compression using finite element method. A micro-polar hypoplastic constitutive model was used. The numerical calculations were carried out with different initial densities and boundary conditions. The behaviour of initially dense, medium dense and loose sand specimen with very smooth or very rough horizontal boundary was investigated. The formation of shear zones gave rise to different global and local stress and strain. Comparisons of the mobilized internal friction, dilatancy and non-coaxiality between global and local quantities were made.     

Numerical simulation of an isolated vortex and shear flow interaction

A numerical solution for the Navier-Stokes equations in the unbounded region is considered for the interaction of an isolated vortex and shear flow. A Chebyshey collocation method in space and finite-difference method for temporal discretization are used. The results of the numerical experiments for the interaction are discussed. It is shown that shear flow can both increase and decrease the vortex dissipation rate.     

Analytic solutions of incompressible laminar flow in the inlet region of a pipe

The laminar analytic solutions of velocities and pressure in the central zone of the inlet region of pipe flow are given under the condition of uniform inflow, based on the Navier-Stokes equations of incompressible viscous flow.     

Explicit adaptive calculations of wrinkled flame propagation

The aim of this work is to study the propagation of a curved premixed flame in an infinite two-dimensional tube. The numerical method combines some features of the finite-element and of the finite-difference methods, and uses a moving adaptive grid procedure in order to reduce the computational costs.     

Experimental study of gaseous detonation propagation over acoustically absorbing walls

The paper presents results from an experimental investigation of the propagation of gaseous detonation waves over tube sections lined with acoustically absorbent materials. The measurements were compared with results from control tests in a smooth wall section. The results show the increasing effectiveness of a perforated steel plate, wire mesh and steel wool in attenuating detonation. Received 25 October 2000 / Accepted 20 August 2001     

Mechanism of compressive stress formation during weak shock waves impact with granular materials

 Certain aspects of wave propagation and the dynamic reaction of a granular material when subjected to a long-duration impulse load are studied. In the majority of studies published on this subject the unsteady pressure behavior at the end-wall covered by a layer of granular material was observed and documented. However, up to now little attention was given to explaining the physical mechanism of this process. Experimental results, obtained in the course of this study, regarding the pressure fields inside granular layers of different materials, clearly show that the compaction effect strongly depends on the characteristics of the medium. This phenomenon manifests itself by changing the gas-particle interaction in the course of the gas filtration, and by variation in the contribution of the different forces and effective stress, σ, to the energy exchange between the gas, the particles and the shock-tube wall. The material permeability,  f, the relative density, ν, and the particle response time, τ _p , are the most important parameters affecting the stress formation at the end-wall covered by the granular layer. In addition to the effect of the material parameters, the effective stress, σ, was found to strongly depend on the granular layer height, h. Based on detailed pressure measurements a qualitative analysis regarding the role of the particle rearrangement in the formation of the unsteady peak at the end-wall was performed. The phenomenology of the particle–particle interaction includes rotation and consolidation of the granules and movement or sliding of particle planes within the layer over each other. Most of these processes are frictional in their nature. They are related to the energy losses and affect the profile and magnitude of the compressive stress as measured at the shock-tube end-wall covered by the granular layer. Received: 10 June 1996/Accepted: 15 October 1996     

A Numerical Study of the Applicability of the Shear Compression Specimen to Parabolic Hardening Materials

This paper addresses quasi-static loading of the shear compression specimen (SCS), that has been especially developed to investigate the shear response of materials at various strain rates. Previous work [4, 5] addressed bi-linear hardening materials, whereas the present work concerns parabolic hardening materials. The investigation is done numerically using three-dimensional elastoplastic finite element simulations. The analyses show that the averaged von Mises stress ( ) and strain ( ) on a mid-section of the gauge reflect accurately the prescribed parabolic hardening model. A method for finding the parabolic hardening coefficients and reducing the measured load, P, and displacement, d, into equivalent stress and strain is introduced and tested. A very good agreement is observed, thus confirming the potential of the technique for large strain testing of parabolic hardening materials.     

Shock and shear initiation of explosive

Three reaction-rate models, Forest Fire, Arrhenius and a combined Forest Fire and Arrhenius, are incorporated in the hydrocode DEFEL. Simulations with this code for PBX-9404 explosive for bare and covered explosives by fragment impact are made. Results confirm the well-known shock-initiation theory for bare explosives. For explosives with a thick cover plate, code results show macroscopic shear initiation, as proposed by Howe (1985) and Frey (1981).This article was processed using Springer-Verlag TEX Shock Waves macro package 1990.     

Effect of axial load on the propagation of elastic waves in helical beams

Helical structures are designed to support heavy loads, which can significantly affect the dynamic behaviour. This paper proposes a physical analysis of the effect of axial load on the propagation of elastic waves in helical beams. The model is based on the equations of motion of loaded helical Timoshenko beams. An eigensystem is obtained through a Fourier transform along the axis. The equations are made dimensionless for beams of circular cross-section and the number of parameters governing the problem is reduced to four (helix angle, helix index, Poisson coefficient, and axial strain). A parametric study is conducted. The effect of loading is quantified in high, medium and low-frequency ranges. Noting that the effect is significant in low frequencies, dispersion curves of stretched and compressed helical beams are presented for different helix angles and radii. This effect is greater as the helix angle increases. Both the effects of stress and geometry deformation are shown to be non-negligible on elastic wave propagation.     

The impact of non-DLVO forces on the onset of shear thickening of concentrated electrically stabilized suspensions

This paper exposes an extension of an activation model previously published by the authors. When particles arranged along the compression axis of a sheared suspension, they may overcome the electrostatic repulsion and form force chains associated with shear thickening. A percolation-based consideration allows an estimation of the impact of the force chains on a flowing suspension. It suggests that similar to mode coupling models, the suspension becomes unstable before the critical stress evaluated from the activation model is reached. The percolated force chains lead to discontinuous shear thickening. The model predictions are compared with results from two experimental studies on aqueous suspensions of inorganic oxides; in one of them, hydration repulsion and in the other hydrophobic attraction can be expected. It is shown that the incorporation of non-Derjaguin–Landau–Verwey–Overbeek forces greatly improve predictions of the shear thickening instability.     

散体系统冲击破碎的动力学分析

针对球形粒子组成的散体系统,基于离散单元法,将球形粒子离散成弹簧-球单元系统,给出了离散单元的运动方程,建立了离散单元之间的弹性力和接触力的计算模型,并用Mohr -Coulomb型破坏准则判断粒子的破碎。运用上述方法,对圆筒内由脆性材料组成的散体系统在冲击载荷下的挤压破碎过程进行了数值模拟;计算过程中,跟踪散体系统中每个粒子在不同时刻的破碎情况;分析了散体系统冲击破碎过程数值模拟结果的主要影响因素。结果显示：数值模拟过程中需综合考虑计算精度和计算时间之间的平衡;相同的计算条件下,颗粒的初始堆积方式不同,计算得到的散体系统的破碎程度不同。

     

Modeling of vibrational impact motion of mobile-based body

The paper deals with vibrational impact motion of a mobile-based body on an inclined plane, which is a characteristic for vibrational alignment of components subjected to an automated assembly. The alignment model is presented by a two-degree-of-freedom vibro-impact system in which the moving body is impacting the plane obliquely. Properties of the vibrational impact motion under kinematical excitement of the moving body and under the excitement of the plane in two perpendicular directions were investigated by numerical methods. It was determined that impact displacement occurs under transient regimes of motion from static to dynamic state of equilibrium. Dependencies of maximum displacement of the body and average displacement velocity on force of elastic resistance, amplitude of plane vibration in the direction of the displacement, and phase angle between the perpendicular members of vibration components were established. Zones of system and excitement parameter combination when the vibrational impact motion occurs were defined.     

Analytical investigation of the fluid flow in the interface region between a porous medium and a clear fluid in channels partially filled with a porous medium

In this paper analytical solutions for the steady fully developed laminar fluid flow in the parallel-plate and cylindrical channels partially filled with a porous medium and partially with a clear fluid are presented. The Brinkman-extended Darcy equation is utilized to model the flow in a porous region. The solutions account for the boundary effects and for the stress jump boundary condition at the interface recently suggested by Ochoa-Tapia and Whitaker. The dependence of the velocity on the Darcy number and on the adjustable coefficient in the stress jump boundary condition is investigated. It is shown that accounting for a jump in the shear stress at the interface essentially influences velocity profiles.     

The effect of an unsteady drag force on the structure of a non-equilibrium region behind a shock wave in a gas-particle mixture

For numerical analysis of shock wave propagation in gas-particle mixtures, drag coefficients of a sphere in steady flows are generally used. However, it is shown both experimentally and numerically that a shock loaded solid sphere experiences unsteady drag forces. The paper describes a model of unsteady drag force and its effect on the structure of the non-equilibrium region behind a shock front traveling in a dusty gas. The results are compared with those obtained by using a steady drag coefficient and are discussed. It is demonstrated that the large drag force at the early stage of the interaction between shock-wave induced flow and a solid particle affects the flow structure that is obtained with a steady drag force.