Measurements of impact properties of small, nearly spherical particles

The authors report impact properties for collisions of small, nearly spherical particles that present interesting experimental challenges. They consider difficulties arising with surface reflectivity, slight asphericity, surface damage and collisions with particles affixed to a rigid plate. To measure these impact properties, the authors refine the experimental technique of Foersteret al. To permit straightforward incorporation in rapid granular theories, the impacts are described with three coefficients. The first is the Newtonian coefficient of normal restitution. The second represents the frictional properties of the contact surfaces. The last characterizes the restitution of the tangential component of the contact point velocity for impacts that involve negligible sliding.     

Collisions in a liquid fluidized bed

Collisional phenomena in a solid–liquid flow were studied in terms of two parameters: the collision frequency and the coefficient of restitution. Experimental measurements of these parameters were conducted inside a liquid fluidized bed by particle tracking in an index-matched array. Collision detection was based on the use of a peak acceleration threshold of the instantaneous speed of colored tracers. The measurements of collision frequency were compared with the theoretical expression derived from the kinetic theory for granular flow (KTGF). The normal and tangential restitution coefficients were measured from the trajectories before and after contact for both particle–particle and particle–wall collisions. A comparison with previous theoretical and experimental works is presented and discussed.     

Three-dimensional simulation of a dust lifting process with varying parameters

This article presents modelling considerations and simulation results for a dust lifting process in a three-dimensional domain. The Eulerian–Lagrangian modelling technique is used. Multiple simulations with different values for the number of particles were performed. The results of the simulations are shown as snapshots of particle position at certain points in time after the passage of a shock wave. Statistical data for the particle positions and collisions are presented. These are: the average height of the particles, the mean square displacement of the particles and the cumulative number of recorded collisions plotted as functions of time. The particle averaged kinetic energy and the mechanical energy lost by particles during collisions are recorded as functions of time in order to study the motion of particles. The results show that simulations of an increasing number of particles render a less intense lifting effect and, more importantly, that the inter-particle and particle–wall collisions represent essential phenomena and need to be included in this type of model. Also, a comparison between two-dimensional and three-dimensional simulations was performed. It was found that, although 2D simulations are still useful, they overestimate the lifting process and therefore a 3D model is preferable. The influence of the magnitude of the restitution and friction coefficients on the process was also studied.     

Inequality of the Coefficients of Restitution for Vertical and Horizontal Quasielastic Impacts of a Ball Agaist a Massive Plate

It is experimentally found out that there is a difference between the coefficients of restitution for horizontal and vertical quasielastic impacts of a steel ball with an initial velocity of no greater than 3 m/sec against a massive steel plate. Possible causes of the difference are considered.     

喷砂冲蚀实验中颗粒轨迹的数值预测

针对冲击磨损实验研究中磨粒群体的运动轨迹难以准确表征的问题,在负压喷射砂粒群冲击Q235钢板的实验中宏观测量了砂粒撞击的速度与位置分布,并使用数值方法模拟了实验砂粒与空气在喷嘴内外的双向耦合过程,以实现负压喷射砂粒群的轨迹预测。计算中提出了非球形粒子在相对马赫数接近1时的曳力模型,以反映空气可压缩引起砂粒表面流动分离的现象,并合理选择Magnus升力模型及壁面反射模型,最终数值预测的砂粒碰撞速度以及撞击位置与实验情况吻合良好。     

A discrete particle approach to simulate the combined effect of blast and sand impact loading of steel plates

The structural response of a stainless steel plate subjected to the combined blast and sand impact loading from a buried charge has been investigated using a fully coupled approach in which a discrete particle method is used to determine the load due to the high explosive detonation products, the air shock and the sand, and a finite element method predicts the plate deflection. The discrete particle method is based on rigid, spherical particles that transfer forces between each other during collisions. This method, which is based on a Lagrangian formulation, has several advantages over coupled Lagrangian-Eulerian approaches as both advection errors and severe contact problems are avoided. The method has been validated against experimental tests where spherical 150 g C-4 charges were detonated at various stand-off distances from square, edge-clamped 3.4 mm thick AL-6XN stainless steel plates. The experiments were carried out for a bare charge, a charge enclosed in dry sand and a charge enclosed in fully saturated wet sand. The particle-based method is able to describe the physical interactions between the explosive reaction products and soil particles leading to a realistic prediction of the sand ejecta speed and momentum. Good quantitative agreement between the experimental and predicted deformation response of the plates is also obtained.     

Plasma diffusion in a toroidal stellarator

Particle trajectories in a gently toroidal stellarator are investigated. The distribution function for the particles is determined in the absence of collisions.     

Modeling of softsphere normal collisions with characteristic of coefficient of restitution dependent on impact velocity

This letter presents a theoretical model of the normal (head-on) collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. After the contact force law between the contacted spheres during a collision is phenomenologically formulated in terms of the compression or overlap displacement under consideration of an elastic—plastic loading and a plastic unloading subprocesses, the coefficient of restitution is gained by the dynamic equation of the contact process once an initial impact velocity is input. It is found that the theoretical predictions of the coefficient of restitution varying with the impact velocity are well in agreement with the existing experimental characteristics which are fitted by the explicit formula.     

Efficient modelling of particle collisions using a non-linear viscoelastic contact force

In this paper the normal collision of spherical particles is investigated. The particle interaction is modelled in a macroscopic way using the Hertzian contact force with additional linear damping. The goal of the work is to develop an efficient approximate solution of sufficient accuracy for this problem which can be used in soft-sphere collision models for Discrete Element Methods and for particle transport in viscous fluids. First, by the choice of appropriate units, the number of governing parameters of the collision process is reduced to one, which is a simple combination of known material parameters as well as initial conditions. It provides a dimensionless parameter that characterizes all such collisions up to dynamic similitude. Next, a rigorous calculation of the collision time and restitution coefficient from the governing equations, in the form of a series expansion in this parameter is provided. Such a calculation based on first principles is particularly interesting from a theoretical perspective. Since the governing equations present some technical difficulties, the methods employed are also of interest from the point of view of the analytical technique. Using further approximations, compact expressions for the restitution coefficient and the collision time are then provided. These are used to implement an approximate algebraic rule for computing the desired stiffness and damping in the framework of the adaptive collision model (Kempe and Fröhlich, J. Fluid Mech. 709: 445–489, 2012). Numerical tests with binary as well as multiple particle collisions are reported to illustrate the accuracy of the proposed method and its superiority in terms of numerical efficiency.     

Precise measurement of soil deformation and fluctuation in drawbar pull for steel and rubber-coated rigid wheels

The travelling performance of rigid wheels on sand stratum is measured using two kinds of surface material, i.e. steel and steel coated with rubber. A new method for measuring the displacement of soil beneath the wheel has been developed using small polyester film markers. The trajectories of soil particles beneath the wheels are approximated by an exponential function and the fluctuations in the drawbar pull are represented by a sinusoidal function. The amplitude and basic wavelength of the fluctuation in the drawbar pull are discussed for both types of wheels.     

Analysis of repeated impacts on a steel rod with visco-plastic material behavior

In machine dynamics impacts are usually common phenomena, resulting from collisions of moving bodies. Even low velocity impacts might produce high stresses in the contact region, which result in inelastic deformation. Thereby, visco-plastic materials, such as steel, show a significant increase of the yield stress with the strain rate. In machine dynamics repeated collisions occur, resulting in repeated impacts on a previously deformed contact area. Then, inelastic deformation and the resulting residual stresses produced by previous impacts have an influence on the behavior of the following impacts. Thus, the impact behavior varies with the number of impacts. This paper presents a numerical and experimental evaluation of repeated impacts with identical impact velocity up to 3 m/s, whereby the deformation history of the contact area, due to previous impacts, is included. The approach is applied to longitudinal impacts of an elastic steel sphere on a steel rod with distinct visco-plastic material behavior which is identified by Split Hopkinson Pressure Bar tests. A Finite Element analysis and experimental verification using two Laser-Doppler-Vibrometers are performed. It is shown that for an accurate impact simulation the FE model must include the visco-plastic material behavior of the steel. Further it is found that the maximal contact force, the rebound velocity and the coefficient of restitution increase with the number of impacts, while the contact duration decreases with the number of impacts. After several impacts these quantities show saturation to a constant value, indicating no significant additional inelastic deformation in the later impacts. Further, the residual stress distribution, the maximal von Mises stress distribution and the local deformation at the contact point are evaluated and a characteristic force-deformation diagram is obtained. Finally, an analysis is performed to describe the relation between maximal force and remaining crater at the contact point.     

On the Direct and Radiated Components of the Collisional Particle Pressure in Liquid-Solid Flows

In a recent study the collisional particle pressure was measured for liquid fluidized beds and liquid-solid flows. The particle pressure was defined as the additional pressure generated by the presence of the particulate-solid phase in a liquid-solid mixture. The particle pressure generated by collisions of particles was found to be composed of two main contributions: one from pressure pulses generated by direct collisions of particles against the containing walls (direct component), and a second one from pressure pulses due to collisions between individual particles that are transmitted through the liquid (radiated component). This paper presents a summary of the technique to measure the particle pressure and the main results of that study.Additional experiments were performed to further study each one of the components of the particle pressure. The direct component was studied by impacting particles on the active face of the pressure transducer. The magnitude of the measured impulse was found to be related to the impact velocity, the mass and the size of the impacting particle. By comparing the measurements with the predictions from Hertzian theory, a quantification of the interstitial fluid effects could be obtained. The radiated component was investigated by generating binary collisions of particles in the vicinity of the transducer. The magnitude of the measured impulse was found to be a function of fluid density, particle size and impact velocity. Predictions based on impulse-pressure theory were obtained and compared with the experimental measurements. The model results showed good agreement with the experimental measurements.     

Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

Particle-wall collision in shear thinning fluids

The present study deals with the measurements of the impact w_i and rebound w_r velocities of steel particles in different fluids colliding with a rigid wall. The results are presented in terms of the coefficient of restitution e=w_r/w_i as a function of the Stokes number (ratio between the particle inertia and the viscous forces). We focus the attention on possible differences between rebounds that occur in Newtonian fluids and in non-Newtonian, shear thinning fluids. The measurements of wet coefficients of restitution in Newtonian fluid are in good agreement with the experimental data found by Gondret et al. (2002). In the range of Stokes number investigated, an increase of the coefficient of restitution with the shear thinning fluid is clearly observed with respect to the Newtonian data. Particular attention has been dedicated to techniques of image processing to perform an optimal estimation of the particle centroid in highly noisy images.     

Dynamic Deformation Behavior of a Golf Ball during Normal Impact

The normal impact between a golf ball and a rigid steel target was studied to examine ball deformation and the contact force during the impact. Using high-speed video images, the normal and tangential compression ratios of the ball were measured to analyze the ball deformation quantitatively. In addition, the inbound and rebound ball velocities, contact time, and coefficient of restitution were determined as basic parameters of the impact. As the inbound ball velocity increased, the maximum normal compression ratio increased while the maximum tangential compression ratio, contact time and coefficient of restitution decreased. The ball center displacements during the impact were measured to determine the ball center velocity and acceleration, and the contact force was calculated by the product of the mass and acceleration. The contact force increased almost linearly with the inbound ball velocity, and its relationship agreed well quantitatively with the results from a load-cell, and also agreed well qualitatively with Hertz contact theory.     

高浓度固-液两相流紊流的动理学模型

采用分子动理学方法,基于固-液两相流液相分子或颗粒相颗粒的Boltzmann方程,对Boltzmann方程分别取零矩和一次矩,则得到高浓度固-液两相流紊流的连续方程和动量方程,再和较成熟的低浓度两相流连续方程和动量方程比较,取低浓度两相流控制方程中较成熟合理的有关项和高浓度时由动理学方法推导出的颗粒间碰撞项,则得到高浓度固-液两相流紊流的最终控制方程:连续方程和动量方程.     

具有摩擦的刚体碰撞

刚体碰撞是力学上的一个经典问题，但目前大都采用给定恢复系数进行分析的方法，本文则直接从两刚体碰撞时Newton第二定律出发，建立了计及摩擦的两刚体碰撞基本理论，并指出了以往研究通过给定恢复系数方法的错误之处。文中通过利用两刚体接触相对变位加速度与碰撞力的关系，定义了法向等效质量，对二维平面碰撞，通过一些简单的比较和定义给出了各种碰撞状态的分类判别法，并给出了碰撞力和冲量的解析式。指出法向等效质量依赖于碰撞状态，在反向滑动和停止滑动这两状态下的法向等效质量是变化的，因此导致其恢复系数是变化的，这样，恢复系数不仅与接触力(接触点情况)有关，而且还与两刚体的运动和动力参数有关，同时给出各种碰撞状态下的恢复系数计算式。     

稠密气固两相流各向异性颗粒相矩方法

基于气体分子动力学和颗粒动理学方法,考虑颗粒速度脉动各向异性,建立颗粒相二阶矩模型.应用初等输运理论,对三阶关联项进行模化和封闭.考虑颗粒与壁面之间的能量传递和交换,建立颗粒相边界条件模型.采用Koch等计算方法模拟气固脉动速度关联矩.考虑气体-颗粒间相互作用,建立稠密气体-颗粒流动模型.数值模拟提升管内气固两相流动特性,模拟结果表明提升管内颗粒相湍流脉动具有明显的各向异性.预测颗粒速度、浓度和颗粒脉动速度二阶矩与Tartan等实测结果相吻合.模拟结果表明轴向颗粒速度脉动强度约为平均颗粒相脉动强度的1.5倍,轴向颗粒脉动能大约是径向颗粒脉动能3.0倍.     

粗糙颗粒动理学及稠密气固两相流动的数值模拟

考虑颗粒碰撞过程中摩擦作用,给出了粗糙颗粒碰撞动力学.引入颗粒相拟总温来表征颗粒平动和转动脉动能量的特征.基于气体分子运动论,建立颗粒碰撞中平动和旋转共同作用的粗糙颗粒动理学,给出了颗粒相压力和黏度等输运参数计算模型.运用基于颗粒动理学的欧拉-欧拉气固两相流模型,数值模拟了流化床内气体颗粒两相流动特性,分析了颗粒旋转流动对颗粒碰撞能量交换和耗散的影响.模拟得到的流化床内径向颗粒浓度和提升管内颗粒轴向速度与他人实验结果相吻合.模拟结果表明随着颗粒浓度的增加,颗粒相压力和能量耗散逐渐增加,而颗粒拟总温先增加后下降.随着颗粒粗糙度系数的增加,床内平均颗粒相拟总温和能量耗散增加,表明颗粒旋转产生的摩擦将导致颗粒旋转脉动能量的改变,影响床内气体-颗粒两相宏观流动特性.