Statistics of particle pair relative velocity in the homogeneous shear flow

Small scale clustering of inertial particles and relative velocity of particle pairs have been fully characterized for statistically steady homogeneous isotropic flows. Depending on the particle Stokes relaxation time, the spatial distribution of the disperse phase results in a multi-scale manifold characterized by local particle concentration and voids and, because of finite inertia, the two nearby particles have high probability to exhibit large relative velocities. Both effects might explain the speed-up of particle collision rate in turbulent flows. Recently it has been shown that the large scale geometry of the flow plays a crucial role in organizing small scale particle clusters. For instance, a mean shear preferentially orients particle patterns. In this case, depending on the Stokes time, anisotropic clustering may occur even in the inertial range of scales where the turbulent fluctuations which drive the particles have already recovered isotropy. Here we consider the statistics of particle pair relative velocity in the homogeneous shear flow, the prototypical flow which manifests anisotropic clustering at small scales. We show that the mean shear, by imprinting anisotropy on the large scale velocity fluctuations, dramatically affects the particle relative velocity distribution even in the range of small scales where the anisotropic mechanisms of turbulent kinetic energy production are sub-dominant with respect to the inertial energy transfer which drives the carrier fluid velocity towards isotropy. We find that the particles’ populations which manifest strong anisotropy in their relative velocities are the same which exhibit small scale clustering. In contrast to any Kolmogorov-like picture of turbulent transport these phenomena may persist even below the smallest dissipative scales where the residual level of anisotropy may eventually blow-up. The observed anisotropy of particle relative velocity and spatial configuration is suggested to influence the directionality of the collision probability, as inferred on the basis of the so-called “ghost collision” model.     

A probability density function of liftoff velocities in mixed-size wind sand flux

With the discrete element method (DEM), employing the diameter distribution of natural sands sampled from the Tengger Desert, a mixed-size sand bed was produced and the particle-bed collision was simulated in the mixed-size wind sand movement. In the simulation, the shear wind velocity, particle diameter, incident velocity and incident angle of the impact sand particle were given the same values as the experimental results. After the particle-bed collision, we collected all the initial velocities of rising sand particles, including the liftoff angular velocities, liftoff linear velocities and their horizontal and vertical components. By the statistical analysis on the velocity sample for each velocity component, its probability density functions were obtained, and they are the functions of the shear wind velocity. The liftoff velocities and their horizontal and vertical components are distributed as an exponential density function, while the angular velocities are distributed as a normal density function. Supported by the Key Project of the National Natural Science Foundation of China (Grant No. 10532040)     

Aeolian transport layer

We investigate the airborne transport of particles on a granular surface by the saltation mechanism through numerical simulation of particle motion coupled with turbulent flow. We determine the saturated flux q(s) and show that its behavior is consistent with classical empirical relations obtained from wind tunnel measurements. Our results also allow one to propose and explain a new relation valid for small fluxes, namely, q(s) = a(u*-u(t))alpha, where u* and u(t) are the shear and threshold velocities of the wind, respectively, and the scaling exponent is alpha approximately 2. We obtain an expression for the velocity profile of the wind distorted by the particle motion due to the feedback and discover a novel dynamical scaling relation. We also find a new expression for the dependence of the height of the saltation layer as a function of the wind velocity.     

Solid-fluid transition in a granular shear flow

The rheology of a granular shear flow is studied in a quasi-2D rotating cylinder. Measurements are carried out near the midpoint along the length of the surface flowing layer where the flow is steady and nonaccelerating. Streakline photography and image analysis are used to obtain particle velocities and positions. Different particle sizes and rotational speeds are considered. We find a sharp transition in the apparent viscosity (eta) variation with rms velocity (u). Below the transition depth we find that the rms velocity decreases with depth and eta proportional to u(-1.5) for all the different cases studied. The material approaches an amorphous solidlike state deep in the layer. The velocity distribution is Maxwellian above the transition point and a Poisson velocity distribution is obtained deep in the layer. The results indicate a sharp transition from a fluid to a fluid + solid state with decreasing rms velocity.     

Anomalous dispersion in correlated porous media: a coupled continuous time random walk approach

We study the causes of anomalous dispersion in Darcy-scale porous media characterized by spatially heterogeneous hydraulic properties. Spatial variability in hydraulic conductivity leads to spatial variability in the flow properties through Darcy’s law and thus impacts on solute and particle transport. We consider purely advective transport in heterogeneity scenarios characterized by broad distributions of heterogeneity length scales and point values. Particle transport is characterized in terms of the stochastic properties of equidistantly sampled Lagrangian velocities, which are determined by the flow and conductivity statistics. The persistence length scales of flow and transport velocities are imprinted in the spatial disorder and reflect the distribution of heterogeneity length scales. Particle transitions over the velocity length scales are kinematically coupled with the transition time through velocity. We show that the average particle motion follows a coupled continuous time random walk (CTRW), which is fully parameterized by the distribution of flow velocities and the medium geometry in terms of the heterogeneity length scales. The coupled CTRW provides a systematic framework for the investigation of the origins of anomalous dispersion in terms of heterogeneity correlation and the distribution of conductivity point values. We derive analytical expressions for the asymptotic scaling of the moments of the spatial particle distribution and first arrival time distribution (FATD), and perform numerical particle tracking simulations of the coupled CTRW to capture the full average transport behavior. Broad distributions of heterogeneity point values and lengths scales may lead to very similar dispersion behaviors in terms of the spatial variance. Their mechanisms, however are very different, which manifests in the distributions of particle positions and arrival times, which plays a central role for the prediction of the fate of dissolved substances in heterogeneous natural and engineered porous materials.     

Rigid and differential plasma crystal rotation induced by magnetic fields

Observations show that plasma crystals, suspended in the sheath of a radio-frequency discharge, rotate under the influence of a vertical magnetic field. Depending on the discharge conditions, two different cases are observed: a rigid-body rotation (all the particles move with a constant angular velocity) and sheared rotation (the angular velocity of particles has a radial distribution). When the discharge voltage is increased sufficiently, the particles may even reverse their direction of motion. A simple analytical model is used to explain qualitatively the mechanism of the observed particle motion and its dependence on the confining potential and discharge conditions. The model takes into account electrostatic, ion drag, neutral drag, and effective interparticle interaction forces. For the special case of rigid-body rotation, the confining potential is reconstructed. Using data for the radial dependence of particle rotation velocity, the shear stresses are estimated. The critical shear stress at which shear-induced melting occurs is used to roughly estimate the shear elastic modulus of the plasma crystal. The latter is also used to estimate the viscosity contribution due to elasticity in the plasma liquid. Further development is suggested in order to quantitatively implement these ideas.     

Suspension and fall of heavy particles in random two-dimensional flow

We investigate the settling of heavy particles in a steady, two-dimensional random velocity field, and find instances in which particle suspension occurs. This leads to a bimodal velocity distribution that may explain some apparently conflicting results reported in the literature. The bimodal distribution is typically smeared out by a time dependence of the ambient flow but, if the time variation is slow, the settling rates of some particles will be as well. The resulting broadbanded velocity distribution of the settling particles will have significance for processes such as rain drop formation, in which the spread of particle velocities affects the statistics of particle collisions.     

湿颗粒聚团碰撞解聚过程的离散元法模拟

基于线性接触模型、库仑滑移接触模型以及平行黏结三种接触模型的组合, 利用离散元法对包衣结构的湿颗粒聚团与壁面碰撞解聚的物理过程进行了数值模拟, 研究了碰撞过程中湿颗粒聚团解聚模式、解聚过程中聚团内各颗粒的速度变化以及颗粒间液桥断裂的规律, 分析了聚团的碰撞速度、黏附小颗粒的重力以及中心大颗粒的旋转对聚团解聚的影响. 研究发现: 聚团的碰撞解聚呈现出碰撞式、重力-碰撞式和剪切-碰撞式三种解聚模式. 湿颗粒聚团与壁面的碰撞打破了聚团内颗粒速度的一致性, 颗粒间出现相对运动而使颗粒间的液桥发生拉伸断裂. 液桥的断裂由聚团的碰撞点向外、由底部向上、由内层向外扩展. 聚团内液桥的断裂经历了缓慢断裂、快速断裂和完全断裂三个阶段. 碰撞速度越大, 黏附的小颗粒质量越大、大颗粒的转速越高, 湿颗粒聚团的缓慢断裂阶段越短暂且解聚程度越高. 模拟结果和实验符合.     

双横波声速法检测单向受压混凝土构件永存应力

该文力求寻找一种高效准确检测混凝土构件永存应力的方法.基于声弹性理论提出了一种双横波声速法检测单向受压混凝土构件永存应力的方法,该方法通过测试受力构件第一波速和第二波速,以第一波速和第二波速的平方差为基础构建综合声学参数来检测单向受压混凝土构件永存应力.所述第一波速为声波传播方向与应力方向垂直、质点振动方向与应力方向平...     

Directional motion of dust particles at different gear structures in a plasma

Directional motion of dust particles in a dusty plasma ratcnet is observed experimentally.The dusty plasma ratcnet consists of two concentric gears with asymmetric sawtooth.It is found that the sawtooth number affects the directional motion of dust particles along the saw channel.With the inerease of the sawtooth number,the particle velocity increases firstly and then decreases,and there is an optimum number of the sawtooth which could induce fast rotation of dust particles.The velocities of dust particles change as they are flowing along the saw channel.We also explore the force acting on the dust particle experimentally.     

CFD-DEM simulation of three-dimensional aeolian sand movement

A three-dimensional CFD-DEM model is proposed to investigate the aeolian sand movement.The results show that the mean particle horizontal velocity can be expressed by a power function of heights.The probability distribution of the impact and lift-off velocities of particles can be described by a log-normal function,and that of the impact and lift-off angles can be expressed by an exponential function.The probability distribution of particle horizontal velocity at different heights can be described as a lognormal function,while the probability distribution of longitudinal and vertical velocity can be described as a normal function.The comparison with previous two-dimensional calculations shows that the variations of mean particle horizontal velocity along the heights in two-dimensional and three-dimensional models are similar.However,the mean particle density of the two-dimensional model is larger than that in reality,which will result in the overestimation of sand transportation rate in the two-dimensional calculation.The study also shows that the predicted probability distributions of particle velocities are in good agreement with the experimental results.     

Nonlinear interaction between two components of motion upon precipitation of a heavy particle in a shear flow

When heavy particles move in a shear flow, the drag depends on the Reynolds number and, hence, on the magnitude of the particle velocity relatively to the medium. This leads to a nonlinear interaction between various components of motion. For example, when a particle precipitates in a horizontal air flow with vertical shear, it also acquires horizontal motion relative to air in addition to vertical motion. These two components of motion contribute to the hydrodynamic drag coefficient by affecting the Reynolds number and thereby influence each other. Steady motion of a particle in a flow with constant vertical shear is analyzed. Dimensionless criterion of significance of the nonlinear effect under consideration is determined. It is demonstrated that this effect can be significant in the near-surface atmospheric layer in the case of storm wind.     

Orientational correlation and velocity distributions in uniform shear flow of a dilute granular gas

Using particle simulations of the uniform shear flow of a rough dilute granular gas, we show that the translational and rotational velocities are strongly correlated in direction, but there is no orientational correlation-induced singularity at perfectly smooth (beta=-1) and rough (beta=1) limits for elastic collisions (e=1); both the translational and rotational velocity distribution functions remain close to a Gaussian for these two limiting cases. Away from these two limits, the orientational as well as spatial velocity correlations are responsible for the emergence of non-Gaussian high-velocity tails. The tails of both distribution functions follow stretched exponentials, with the exponents depending on normal (e) and tangential (beta) restitution coefficients.     

On aggregate structures in a rod-like haematite particle suspension by means of Brownian dynamics simulations

We have investigated aggregation phenomena in a suspension composed of rod-like haematite particles by means of Brownian dynamics simulations. The magnetic moment of the haematite particles lies normal to the particle axis direction and therefore the present Brownian dynamics method takes into account the spin rotational Brownian motion about the particle axis. We have investigated the influence of the magnetic particle–field and particle–particle interactions, the shear rate and the volumetric fraction of particles on the particle aggregation phenomena. Snapshots of aggregate structures are used for a qualitative discussion and the cluster size distribution, radial distribution function and the orientational correlation functions of the direction of particle axis and magnetic moment are the focus for a quantitative discussion. The significant formation of raft-like clusters is found to occur at a magnetic particle–particle interaction strength much larger than that required for a magnetic spherical particle suspension. This is because the rotational Brownian motion has a significant influence on the formation of clusters in a suspension of rod-like particles with a large aspect ratio. An applied magnetic field enhances the formation of raft-like clusters. A shear flow does not have a significant influence on the internal structure of the clusters, but influences the cluster size distribution of the raft-like clusters.     

Derivation of the probability distribution of velocities of a detector-recorded brownian particle

The probability distribution of velocities in the given space region (detector) is found for particles of a passive admixture in a stream of external gas. Since direct calculation of the above probability density involves significant difficulties, the solution is based on the classical problem of the probability distribution of coordinates and velocity of a Brownian particle at a fixed time. Analyzing dependence of the solution on the parameters of the initial problem, we obtain conditions under which the assumptions on the character of particle motion hold true. State University, Nizhny Novgorod, Russia. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 41, No. 10, pp. 1301–1313, October 1998.     

陶瓷材料在压剪联合冲击加载下动态响应的实验研究

 应用57 mm压剪炮和双磁场IMPS粒子测试系统对95%的Al₂O₃陶瓷材料在压剪联合加载条件下的力学行为进行了实验研究,测量得到了材料内部压缩波（P波）、剪切波（S波）的加载、卸载和传播规律。对剪切波的衰减随载荷的变化进行了研究,初步得到95% Al₂O₃陶瓷材料在压剪联合载荷冲击下发生损伤的纵向应力阀值是4.86 GPa。     

Influence of supraglottal structures on the glottal jet exiting a two-layer synthetic, self-oscillating vocal fold model

A synthetic two-layer, self-oscillating, life-size vocal fold model was used to study the influence of the vocal tract and false folds on the glottal jet. The model vibrated at frequencies, pressures, flow rates, and amplitudes consistent with human phonation, although some differences in behavior between the model and the human vocal folds are noted. High-speed images of model motion and flow visualization were acquired. Phase-locked ensemble-averaged glottal jet velocity measurements using particle image velocimetry (PIV) were acquired with and without an idealized vocal tract, with and without false folds. PIV data were obtained with varying degrees of lateral asymmetric model positioning. Glottal jet velocity magnitudes were consistent with those measured using excised larynges. A starting vortex was observed in all test cases. The false folds interfered with the starting vortex, and in some cases vortex shedding from the false folds was observed. In asymmetric cases without false folds, the glottal jet tended to skew toward the nearest wall; with the false folds, the opposite trend was observed. rms velocity calculations showed the jet shear layer and laminar core. The rms velocities were higher in the vocal tract cases compared to the open jet and false fold cases.     

Brownian motion and correlation in particle image velocimetry

In particle image velocimetry applications involving either low velocities or small seed particles, Brownian motion can be significant. This paper addresses the effects of Brownian motion. First, general equations describing cross-correlation particle image velocimetry are derived that include Brownian motion. When light-sheet illumination particle image velocimetry (PIV) is used Brownian motion diminishes the signal strength. A parameter describing this effect is introduced, and a weighting function describing the contribution to the measured velocity as a function of position is derived. The latter is unaffected by Brownian motion. Microscopic PIV Brownian motion also diminishes the signal strength. The weighting function for microscopic PIV is found to depend on Brownian motion, thus affecting an important experimental parameter, the depth of correlation. For both light-sheet illumination and microscopic PIV, a major consequence of Brownian motion is the spreading of the correlation signal peak. Because the magnitude of the spreading is dependent on temperature, PIV can, in principle, be used to simultaneously measure velocity and temperature. The location of the signal peak provides the velocity data, while the spreading of the peak yields temperature.     

On the orientational characteristics and transport coefficients of an oblate spheroidal hematite particle in a simple shear flow

We have developed the basic equation of the orientational distribution function of oblate spheroidal hematite particles with rotational Brownian motion in a simple shear flow under an applied magnetic field. An oblate spheroidal hematite particle has an important characteristic in that it is magnetized in a direction normal to the particle axis. Since a dilute dispersion is addressed in the present study, we have taken into account only the friction force (torque) whilst neglecting the hydrodynamic interactions among the particles. This basic equation has been solved numerically in order that we may investigate the dependence of the orientational distribution on the magnetic field strength, shear rate and rotational Brownian motion and the relationship between the orientational distribution and the transport coefficients such as viscosity and diffusion coefficient. We found that if the effect of the magnetic field is more dominant, the particle inclines in such a way that the oblate surface aligns in the magnetic field direction. If the Peclet number increases and the effect of the shear flow becomes more dominant, the particle inclines such that the oblate surface tilts in the shear flow direction. The viscosity due to the magnetic torque is shown to increase as the magnetic field increases, since the magnetic torque due to the applied magnetic field becomes the more dominant effect. Moreover, the viscosity increase is shown to be more significant for a larger aspect ratio or for a more oblate hematite particle. We have applied the analysis to the problem of particle sedimentation under gravity in the presence of a magnetic field applied in the sedimentation direction. The particles are found to sediment with the oblate surface aligning more significantly in the sedimentation direction as the applied magnetic field strength increases.     

A Mechanical Model of Brownian Motion for One Massive Particle Including Slow Light Particles

We provide a connection between Brownian motion and a classical mechanical system. Precisely, we consider a system of one massive particle interacting with an ideal gas, evolved according to non-random mechanical principles, via interaction potentials, without any assumption requiring that the initial velocities of the environmental particles should be restricted to be “fast enough”. We prove the convergence of the (position, velocity)-process of the massive particle under a certain scaling limit, such that the mass of the environmental particles converges to 0 while the density and the velocities of them go to infinity, and give the precise expression of the limiting process, a diffusion process.