Lattice Boltzmann simulation of behaviour of particles moving in blood vessels under the rolling massage

The rolling massage is one of the most important manipulations in Chinese massage, which is expected to eliminate many diseases. Here, the effect of the rolling massage on a pair of particles moving in blood vessels under rolling massage manipulation is studied by the lattice Boltzmann simulation. The simulated results show that the motion of each particle is considerably modified by the rolling massage, and it depends on the relative rolling velocity, the rolling depth, and the distance between particle position and rolling position. Both particles' translational average velocities increase almost linearly as the rolling velocity increases, and obey the same law. The increment of the average relative angular velocity for the leading particle is smaller than that of the trailing one. The result is helpful for understanding the mechanism of the massage and to further develop the rolling techniques.     

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.     

激光熔覆中金属粉末粒子与激光相互作用模型

为了对同轴激光熔覆过程中运动的金属粉末粒子的速度和温度进行理论分析,并研究各工艺参量的影响,建立了运动中金属粉末粒子的运动模型和热模型.模拟结果表明,粉嘴几何尺寸、粒子直径以及气/粉两相流初始速度是影响粒子运动行为的重要因素;粉嘴几何尺寸、激光焦点位置、激光发散角、激光功率、粒子直径以及气/粉两相流初始速度是影响粒子热行为的重要因素.在相同的工艺参量下(粉嘴出口内径r=2 mm,粉嘴倾角α=60°,初始气流速度v0=0.8 m/s),基于数字粒子图像测速(DPIV)技术,对316L不锈钢粉末粒子运动模型进行了实验验证.结果表明,运动理论模型是可靠的.该模型是掌握同轴激光熔覆过程中金属粉末粒子运动行为的有效工具;同时,热模型也是分析粉末粒子温度随不同参量变化的重要工具.     

The response of turbulent premixed flames to normal acoustic excitation

To model the thermo-acoustic excitation of flames in practical combustion systems, it is necessary to know how a turbulent flame front responds to an incident acoustic wave. This will depend partly on the way in which the burning velocity responds to the wave. In this investigation, the response of CH₄/air and CH₄/H₂/air mixtures has been observed in a novel flame stabilisation configuration, in which the premixture of fuel and air is made to decelerate under controlled conditions in a wide-angle diffuser. Control is provided by an annular wall-jet of air and by turbulence generators at the inlet. Ignition from the outlet of the diffuser allows an approximately flat flame to propagate downwards and stabilise at a height that depends on the turbulent burning velocity. When the flow is excited acoustically, the ensemble-averaged height oscillates. The fluctuations in flow velocity and flame height are monitored by phase-locked particle image velocimetry and OH-planar laser induced fluorescence, respectively. The flame stabilised against a lower incident velocity as the acoustic amplitude increased. In addition, at the lowest frequency of 52 Hz, the fluctuations in turbulent burning velocity (as represented by the displacement speed) were out-of-phase with the acoustic velocity. Thus, the rate of displacement of the flame front relative to the flow slowed as the flow accelerated, and so the flame movement was bigger than it would have been if the burning velocity had not responded to the acoustic fluctuation. With an increase in frequency to 119 Hz, the relative flame movement became even larger, although the phase-difference was reduced, so the effect on burning velocity was less dramatic. The addition of hydrogen to the methane, so as to maintain the laminar burning velocity at a lower equivalence ratio, suppressed the response at low amplitude, but at a higher amplitude, the effect was reversed.     

Simulation of triaxial response of granular materials by modified DEM

A modified discrete element method(DEM)with rolling effect taken into consideration is developed to examine macroscopic behavior of granular materials in this study.Dimensional analysis is firstly performed to establish the relationship between macroscopic mechanical behavior,mesoscale contact parameters at particle level and external loading rate.It is found that only four dimensionless parameters may govern the macroscopic mechanical behavior in bulk.The numerical triaxial apparatus was used to study their influence on the mechanical behavior of granular materials.The parametric study indicates that Poisson’s ratio only varies with stiffness ratio,while Young’s modulus is proportional to contact modulus and grows with stiffness ratio,both of which agree with the micromechanical model.The peak friction angle is dependent on both inter-particle friction angle and rolling resistance.The dilatancy angle relies on inter-particle friction angle if rolling stiffness coefficient is sufficiently large.Finally,we have recommended a calibration procedure for cohesionless soil,which was at once applied to the simulation of Chende sand using a series of triaxial compression tests.The responses of DEM model are shown in quantitative agreement with experiments.In addition,stress-strain response of triaxial extension was also obtained by numerical triaxial extension tests.     

Hydrodynamic model for a dynamical jammed-to-flowing transition in gravity driven granular media

Granular material on an inclined plane will flow like a fluid if the angle theta the plane makes with the horizontal is large enough. We study chute flow down a plane using a hydrodynamic model previously used to describe granular Couette flow. Our model predicts a jammed-to-flowing transition as theta is increased even though it does not include solid friction, which might seem necessary to stabilize a state without flow. The transition is driven by coupling between mean and fluctuating velocity. In agreement with experiments and simulations, it predicts flow for layers with a thickness H larger than a critical value H(stop)(theta) and absence of flow for H相似文献   

Solids Movement in Rotary Kilns in the Slumping Regime: Model Using a Control Plane Parallel to the Steepest Descent

A new approach is proposed to model the bulk movement of solids in rotary drums operating at low rotation speeds, in slumping and rolling regimes. The model yields an equation similar to Saeman's equation, but which is valid also for the slumping regime and active area in the rolling regime. The model was developed for constant depth using a control surface containing the steepest descent direction, so that any contribution from sliding particles to flow rate can be neglected. By considering an appropriate virtual kiln, the model is extended to the more general, variable depth situation.     

也谈“用绳拉船靠岸的船速问题”

指出解质点运动学问题中运用的所谓“速度合成与分解方法”有两种含义,第一种是只在一个参考系中求解,运用速度的矢量性,将速度矢量按平行四边形法则合成和分解,最常用的是在直角坐标系或其它坐标系中合成和分解;第二种是借助两个有相对运动的参考系,运用相对运动中的速度定理求解．同时澄清了有关文献中某些不当之处．     

受碾区域内颗粒轴向流动特性的离散元模拟

为探讨受碾状态颗粒的稳定流动, 在碾辊轴与筛筒组成的受碾区域内, 建立了轴向运动的颗粒流离散元物理模型. 研究结果表明: 受碾区域内各颗粒沿轴向运动能力的差异造成了颗粒流密度不均匀; 颗粒与筛筒间的静摩擦系数影响颗粒轴向流动的形态、速率及集散程度, 受碾区域内单层颗粒的轴向均方偏差与流动时间的平方正相关, 属于“super”扩散; 整体分析受碾区域发现, 颗粒的轴向平均速度沿轴向坐标逐渐增大, 而颗粒的三轴合成平均速度沿轴向坐标逐渐降低; 受碾区域内各轴向位置处颗粒运动的剧烈程度不同, 沿轴向坐标颗粒的波动速度平方呈现先增大后降低而后又增大的趋势; 单颗粒的碰撞总能量损失能谱也表明了颗粒运动程度不同, 即轴向流动时在受碾区域的前半段碰撞剧烈, 能量损失多, 在后半段碰撞程度弱, 能量损失较少. 通过对受碾区域内颗粒流动的数值模拟分析, 明晰了颗粒在受碾条件下稳定流动特性, 有益于碾磨工业对产品品质控制及设备参数优化的研究.     

Discrete Velocity Fields with Explicitly Computable Lagrangian Law

We introduce a class of random velocity fields on the periodic lattice and in discrete time having a certain hidden Markov structure. The generalized Lagrangian velocity (the velocity field as viewed from the location of a single moving particle) has similar hidden Markov structure, and its law is found explicitly. Its rate of convergence to equilibrium is studied in small numerical examples and in rigorous results giving absolute and relative bounds on the size of the second–largest eigenvalue modulus. The effect of molecular diffusion on the rate of convergence is also investigated; in some cases it slows convergence to equilibrium. After repeating the velocity field periodically throughout the integer lattice, it is shown that, with the usual diffusive rescaling, the single–particle motion converges to Brownian motion in both compressible and incompressible cases. An exact formula for the effective diffusivity is given and numerical examples are shown.     

对颗粒物质运动的一致性进行控制的随机力场

在颗粒流的研究中引入了一个正态分布的随机力场, 并通过计算机模拟研究了该力场对均匀颗粒流的影响. 结果发现: 随机力场基本上不改变均匀颗粒流的平均密度和速度, 对颗粒流密度的涨落也影响很小. 随机力场对均匀颗粒流的影响主要体现在它可提高速度的涨落, 它与颗粒体系的耗散性质相抗衡, 使颗粒流维持一定的波动能量. 研究结果还显示: 通过随机力场所获得的波动能量并没有均匀分布到各个自由度上, 由于颗粒体系的耗散性质颗粒体系难以达到能均分状态. 关键词： 颗粒物质 随机力 分子动力学模拟     

一种新的光子多普勒速度频谱分析方法

光子多普勒速度计可给出飞层表面某一速度带内颗粒群速度随时间演化的频谱数据, 在冲击动力学实验尤其是微喷射及其混合研究中得到广泛应用. 本文提出一种新的光子多普勒频谱数据分析方法, 可推断出混合区厚度变化和前端等效颗粒尺度. 利用该方法, 对一些典型状态下喷射混合速度频谱开展分析, 获得了不同冲击压力、气体条件下颗粒度数据, 证实了气体环境下喷射颗粒的气动破碎现象, 以及破碎后尺度与初始条件的依赖性, 为喷射混合物理规律研究提供了重要依据.     

The Thermal Conductivity Theory of Non-uniform Granular Flow and the Mechanism Analysis

According to the fractal characteristics appearing in non-uniform granular system, we found the fractal model to study the effective thermal conductivity in the mixed system. Considering the quasi-equilibrium, we bring forward the fractal velocity probability distribution function. The equipartition of energy is employed to the non-uniform granular system, and the granular temperature is derived. We investigate the thermal conductivity in granular flow due to the movement of the particles, namely the heat transfer induced by the streaming mode only. The thermal conductivity in the mixed system changes with the fractal parameters such as the solid fraction v, structural character parameter η, and fractal dimension D of size distribution. These parameters depict the characteristics of the thermal conductivity in the actual complex granular system. Comparing our conclusion with the correlative experimental data and the theoretical conclusion of binary mixture of granular materials, the results can qualitatively confirm the generality of our prediction on the granular system.     

Shear zones and wall slip in the capillary flow of concentrated colloidal suspensions

We image the flow of a nearly random close packed, hard-sphere colloidal suspension (a "paste") in a square capillary using confocal microscopy. The flow consists of a "plug" in the center while shear occurs localized adjacent to the channel walls, reminiscent of yield-stress fluid behavior. However, the observed scaling of the velocity profiles with the flow rate strongly contrasts yield-stress fluid predictions. Instead, the velocity profiles can be captured by a theory of stress fluctuations originally developed for chute flow of dry granular media. We verified this both for smooth and rough walls.     

稠密颗粒射流倾斜撞击颗粒膜特征

采用高速摄像仪对稠密颗粒射流倾斜撞击形成的类液体颗粒膜特征进行实验研究,考察了颗粒粒径、射流速度以及射流含固率等因素对颗粒膜形态及动态特征的影响.结果表明:随着颗粒粒径增大,稠密颗粒倾斜撞击流由颗粒膜向散射模式转变;随着射流速度增加,气固不稳定增强,射流流量出现脉动,正面与侧面分别表现为颗粒膜的非轴对称振荡和表面波纹结构;颗粒膜非轴对称振荡的振荡频率和振荡幅度随射流速度的增大而增大;表面波纹速度和波长沿传播方向增大,波纹间存在叠加现象.颗粒膜出现非轴对称振荡主要是因为喷嘴出口由气固不稳定性引起的射流流量脉动,射流流量脉动频率与撞击面振荡频率基本相当.     

Computing particle velocity and intensity in diffracting systems

In a previous paper the authors have shown how it is possible to compute the acoustic pressure due to an axi-symmetric radiator with different boundary conditions. The technique involves the conversion of a double integral to a single integral by means of a geometric transform, and may be applied with the Green's function appropriate to each boundary condition.The present work shows how the form of the single integrals makes it possible to derive and compute both the particle velocity and the wave intensity. It has become apparent that the particle velocity must be described as an elliptic vector. A new concept, the ‘diffraction front’, has been defined as a surface in which there is no net energy transfer, and this has been applied to the evaluation of the magnitude and direction of energy flow. Plots have been made of flow lines and diffraction fronts, and also three-dimensional representations of energy intensity have been made for those acoustic fields in which pressure has already been evaluated.     

稠密颗粒射流撞击壁面颗粒膜表面波纹特征

采用高速摄像仪对稠密颗粒射流撞击有限尺寸壁面的流动过程进行了实验研究,重点研究了颗粒膜及其表面波纹特征,考察了颗粒粒径、射流速度和固含率等因素对颗粒膜形态和表面波纹的影响.研究结果表明,随着颗粒粒径增大,稠密颗粒撞壁流由颗粒膜向散射模式转变.与液体射流撞壁液膜相比,颗粒膜扩展角较大,射流速度对其影响不显著.稠密颗粒射流撞壁颗粒膜表面波纹存在明显的叠加现象,颗粒膜表面波纹频率比液膜大约低一个数量级.颗粒膜表面波纹主要由射流脉动引起,表面波纹频率与射流脉动频率具有相同的数量级.     

Phase transitions in shear-induced segregation of granular materials

We computationally study shear-induced segregation of different-sized particles in vertical chute flow. We find that, for low solid fractions, large particles segregate toward regions of low shear rates where the granular temperature (velocity variance) is low. As the solid fraction increases, this trend reverses, and large particles segregate toward regions of high shear rates and temperatures. We find that this is a global phenomenon: local segregation trends reverse at high system solid fractions even where local solid fractions are small. The reversal corresponds to the growth of a single enduring cluster of 30%-60% of the particles that we propose changes the segregation dynamics.     

Rolling Friction

Doklady Physics - The dependence of rolling friction on velocity for various contact conditions is discussed. The principal difference between rolling and other types of relative motion (sliding...     

The Effect of the Relative Velocity on Traffic Flow

The optimal velocity model of traffic is extended to take the relative velocity into account. The traffic behavior is investigated numerically and analytically with this model. It is shown that the car interaction with the relative velocity can effect the stability of the traffic flow and raise critical density. The jamming transition between the freely moving and jamming phases is investigated with the linear stability analysis and nonlinear perturbation methods. The traffic jam is described by the kink solution of the modified Korteweg--de Vries equation. The theoretical result is in good agreement with the simulation.