相互碰撞的圆粒子在竖直通道中沉降的数值研究

采用格子Boltzmann方法,数值模拟了考虑相互碰撞的两圆粒子的沉降,分析了Re数、初始相对位移及通道宽度对粒子沉降的影响．结果表明,在0.1＜Re＜20范围内, 粒子沉降具有周期性．Re数越大,两圆粒子的相互作用越强,粒子横向位移的幅度也越大．在大Re数时,沉降的过程是两个粒子交替领先;在较小Re数时,当后面的粒子接近的时候,领先的粒子被向右侧推了一段后仍会继续领先;对中等Re数,原先在后面的粒子在第一次加速后将取得领先位置并一直保持下去．粒子的初始分布位置对沉降的形态影响不大．管道宽度变化时,粒子总的沉降特性不变,而周期改变,管道越宽,周期越长．     

纳米粒子在弯管中的输运和沉降特性

用摄动法对不同Reynolds数和Schmidt数的纳米粒子在圆截面弯管中的运输和沉积进行了求解.结果表明,当悬浮纳米粒子在直管中流动时,粒子输运模式不依赖于粒子的大小和其它参数.在弯管中运动时,管道外弯侧具有最多的沉积粒子而内弯侧的沉积粒子最少.在管道的上部和下部,不同Schmidt数粒子的沉降特性一致.管道曲率、Reynolds数和Schmidt数对粒子相对沉积效率的影响具有二阶,四阶和一阶的作用.     

三维矩形槽道中颗粒沉降的数值模拟

采用三维格子Boltzmann方法对矩形通道中的颗粒沉降进行了模拟研究．单颗粒沉降的模拟结果表明,颗粒最终的稳定沉降位置沿槽道中心线,不受颗粒初始位置和直径的影响．颗粒和壁面之间的两体相互效应可以用无因次沉降速度定量描述,无因次沉降速度的模拟结果和实验结果定量上吻合一致．模拟分析了双颗粒沉降的DKT（drafting, kissing and tumbling）过程,探讨了颗粒直径比以及壁面效应对DKT过程的影响．模拟发现当颗粒直径相同时,双颗粒的沉降过程为周期性的DKT过程,从而形成双螺旋形式的沉降轨迹,此螺旋沉降轨迹的频率和振幅受颗粒初始位置影响．从模拟结果中还得到颗粒群的最终稳定构型,并进行了构型对比分析．最后对包含49个颗粒的颗粒群沉降行为进行了模拟,说明多体相互作用在对称性的情况下可以简化．     

月尘静电浮扬现象的理论模拟

对月球上的月尘静电浮扬现象进行了理论研究,对比分析了影响月球表面静电浮扬强弱的主要因素．研究过程分为两个步骤:首先采用一维PIC（particle in cell）模拟计算了月尘和月球表面的充电过程,然后基于这一结果引进试验粒子,对月尘的静电浮扬现象进行了分析研究．结果表明月尘的静电浮扬主要受两个因素的影响:太阳角和月尘颗粒的大小．月尘静电浮扬现象在日出日落时分更容易发生,即太阳角越小越容易引起剧烈的月尘浮扬现象;并且月尘粒径越小,其浮扬高度越高．     

球面介质作用下的测度值分枝过程

我们考虑空间上一粒子系统，当其受到分布于求面上的介质作进行粒子分枝和衍生，产生新了体，而新粒子仍按原粒子的运动规则继续空间运动。通过合理的假设和极限过程，粒子在空间的散布一测度值分枝过程来刻划。     

曲率的形状梯度和经典梯度:微纳米曲面上的驱动力

近期的实验和分子动力学模拟均表明:圆锥面上粘附液滴能自发地定向运动,且自发定向运动的方向与粘附面的亲水、疏水性质无关．针对这一重要现象,拟从曲面微纳米力学几何化的角度,提供一般性的理论解释．借助于粒子对势,研究了孤立粒子与微纳米硬曲面之间的相互作用,分析了粒子/硬曲面相互作用的几何学基础．可以证实:（a） 粒子/硬曲面的作用势均具有统一的曲率化形式,均可以统一地表达成曲面平均曲率和Gauss曲率的函数;（b） 基于曲率化的作用势,能够实现曲面微纳米力学的几何化;（c） 曲率与曲率的内蕴梯度构成卷曲空间上的驱动力;（d） 驱动力方向与曲面的亲水、疏水性质无关,解释了自发定向运动实验．     

梢涡流场中气泡轨迹及形态数值模拟

针对气泡在舰船尾迹涡流场运动特性,根据其是否为尾涡所捕获,将数值模拟过程分为两个阶段:准球状运动阶段和非球状运动阶段．分别应用单向耦合质点粒子追踪法（PTM）和边界元法（BEM）模拟这两个阶段,将第1阶段结束的物理量作为第2阶段的初始条件,从而完成整个数值模拟过程．在已有数值研究结果和实验数据基础上,探讨空化发生条件,追踪尾迹空泡运动轨迹,模拟尾迹气泡的运动、变形、溃灭等,以及被尾涡捕获后的撕裂等运动特性,旨在为优化设计尾流场提供参考．     

纳米单晶与多晶铜薄膜力学行为的数值模拟研究

通过对不同温度下单晶薄膜的拉伸性能的分子动力学模拟，从微观角度揭示了温度效应对材料性能的影响. 结果表明温度效应对材料的变形机理影响很大.0K温度下由于缺乏热激活软化的影响, 粒子运动所受到的阻碍较大, 薄膜的强度较高, 塑性变形主要来自于粒子的短程滑移.温度升高，粒子的热运动加剧，屈服强度降低, 塑性变形将主要来自于大范围的位错长程扩展.多晶薄膜的模拟结果表明, 虽然其晶粒形状较为特殊, 但是它仍然遵循反Hall-Petch关系.在模拟过程中，侧向应力最大值比拉伸方向应力的最大值滞后出现.位错只会从晶界产生并向晶粒内部传播，晶粒间界滑移是多晶薄膜塑性变形的主要来源.     

考虑浅埋破坏的拖曳锚在黏土中安装运动特性分析

拖曳锚是海洋工程中一种常见的系泊基础,因造价低廉和高承载特性而得到广泛应用.其在海床中的安装轨迹和运动特性受到锚与土体之间复杂相互作用的影响,使得精确定位仍存在挑战.目前已有的塑性屈服面方法被广泛用于计算拖曳锚的运动特性,即假设整个拖曳过程为深埋板在不同深度的破坏过程.实际上,拖曳锚的安装是从浅埋到深埋的连续贯入,因此该方法不能考虑浅埋破坏对拖曳锚运动特性的影响,从而可能导致预测的轨迹不准确.通过有限元分析研究了锚板方位角及埋深比对单向承载和复合荷载下屈服面的影响,确定了锚板浅埋破坏时的屈服面,补充了塑性屈服面法对浅埋破坏效应的考虑;进一步地,考察了锚板方位角、承载系数、浅埋区域大小对拖曳锚轨迹预测和运动特性的影响,并与传统仅假设深埋破坏情况对比分析.结果表明:浅埋破坏时锚板方位角与浅埋区域大小决定了锚板的运动特性和轨迹;合理考虑浅埋破坏后,与纯假设深埋破坏情况比,锚板在达到稳定状态之前的预测埋深和锚链力要小,但极限嵌入深度一致.     

湍动尺度的模糊聚类分析

湍流运动可看成是大小不同尺度的涡体运动的叠加.定量地确定湍动尺度的分类,对于更好地描述不同尺度的涡体运动,探讨不同尺度涡体之间的相互作用,建立较好的湍流模式都具有重要的意义. 对事物按一定要求进行分类的数学方法,叫做聚类分析.由于湍动尺度的分类具有一定程度的模糊性,因而本文采用模糊聚类的方法,对壁面光滑及租糙两种边界条件下的湍动尺度进行了分类,并对各类结构的特性及其相互作用进行了探讨.     

Numerical simulation of dynamic sand dunes

The physics of granular materials is interesting from many points of view because they exhibit a wealth of phenomena that have both fluid and solid aspects [C.S. Campbell, Annu. Rev. Fluid. Mech. 22 (1990) 57, H.M. Jaeger, S.R. Nagel, R.P. Behringer, Phys. Today 494 (1996) 32]. Recently a difficult pattern was observed if sand falls in the space between two plates and passes an obstacle [Y. Amarouchene, J.F. Boudet, H. Kellay, Phys. Rev. Lett. 86 (2001) 4286]. The interesting behaviour occurs on top of the obstacle where a dynamic dune with a parabolic tip is formed. Inside this parabola, a triangular region of non- or very slow flowing sand is observed. Using factor analysis it is possible to extract latent parameters from a dynamic process. Applying a three factor model we can clearly identify the inner triangle (1st factor) and the outer parabolic pattern (3rd factor). The second factor we interpret as shock wave. Most interactions between particles take place in a relatively small region. We show that the pattern formation process depends on the restitution coefficients (particle–particle and particle–obstacle) and also on the particle size. These findings cannot be observed if standard velocity profiles are used to analyse the data. Our findings show, that most interactions take place in a relatively small area correlating with the particle size. If the interactions between different particles and particle–obstacle are elastic the formation of a non-flowing triangular region is more difficult as if inelastic collisions are used. The factor curves also clearly show that a pattern formation process has to be finished, before the next pattern can be formed.     

Computation of Screening Phenonema in a Vertical Tumbling Cylinder

Granular materials are an integral part of our environment. Due to their wide variety of applications in industrial and technological processes, they have captured a great interest in the recent research, see [1] and [2]. The related studies are often based on numerical simulations and it is considered as challenging to investigate computational phenomena of dense granular systems. Particle screening is an essential technology in many industrial fields and important in granular studies. The particular problem of interest is the separation of round shape particles of different geometrical sizes using a rotating tumbling vertical cylinder. The concept of discrete element method (DEM) that considers the motion of each single particle individually is applied in this study. Particle-to-particle and particle-to-wall collisions will appear under the tumbling motion of the rotating structure. The normal and frictional forces between particles themselves and particles and surrounding walls of the structure are calculated according to the rules of a penalty method, which employs spring-damper models for this purpose. As a result of collisions, the particles will dissipate kinetic energy due to the normal and frictional contact losses. Particle distribution and sifting rate of the separated particles have been studied taking into consideration different rotational speeds of the machine, various damping and frictional coefficients and different sizes of holes in the sifting plates at different levels of the structure. In an attempt to better understand the mechanism of the particle transport between the different layers of the sifting system, different computational studies have been performed. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling and simulation of the motion of nanoparticles in cylindrical capillaries allowing particle‐to‐wall interactions

The process of transporting nanoparticles at the blood vessels level stumbles upon various physical and physiological obstacles; therefore, a Mathematical modeling will provide a valuable means through which to understand better this complexity. In this paper, we consider the motion of nanoparticles in capillaries having cylindrical shapes (i.e., tubes of finite size). Under the assumption that these particles have spherical shapes, the motion of these particles reduces to the motion of their centers. Under these conditions, we derive the mathematical model, to describe the motion of these centers, from the equilibrium of the gravitational force, the hemodynamic force and the van der Waals interaction forces. We distinguish between the interaction between the particles and the interaction between each particle and the walls of the tube. Assuming that the minimum distance between the particles is large compared with the maximum radius R of the particles and hence neglecting the interactions between the particles, we derive simpler models for each particle taking into account the particles‐to‐wall interactions. At an error of order O(R) or O(R³)(depending if the particles are 'near' or 'very near' to the walls), we show that the horizontal component of each particle's displacement is solution of a nonlinear integral equation that we can solve via the fixed point theory. The vertical components of the displacement are computable in a straightforward manner as soon as the horizontal components are estimated. Finally, we support this theory with several numerical tests. Copyright © 2016 John Wiley & Sons, Ltd.     

Direct numerical simulation of bi-disperse particle-laden gravity currents in the channel configuration

We present a numerical investigation of bi-disperse particle-laden gravity currents in the lock-exchange configuration. Previous results, based on numerical simulation and laboratory experiments, are used to establish comparisons. Our discussion focuses on explaining how the presence of more than one particle diameter influences the main features of the flow, such as deposit profile, the evolution of the front location and suspended mass. We develop the complete energy budget equation for bi-disperse flows. A set of two and three-dimensional direct numerical simulations (DNS), with different initial compositions of coarse and fine particles, are carried out for Reynolds number equal to 4000. Such simulations show that the energy terms are strongly affected by varying the initial particle fractions. The addition of a small amount of fine particles into a current predominantly composed of coarse particles increases its run-out distance. In particular, it is shown that higher amounts of coarse particles have a dumping effect on the current development. Comparisons show that the two-dimensional simulation does not reproduce the intense turbulence generated in 3D cases accurately, which results in a significant difference in the suspended mass, front position as well as the dissipation term due to the advective motion.     

Paramagnetic Meissner effect of highT c superconductors II

A range of powdered Bi:2 212 samples exhibiting the paramagnetic Meissner effect (PME) are systematically examined. Interpretation of the results is made in terms of a phenomenological model in which there is a concentration within the material of small local moments that can be polarized during a field cooling. Information about the magnitudes of these local m0oments and their distribution are deduced. Relations between the local moments and the particle sizes, the weak link, oxygen content and the interactions between the local moments are also discussed. Comparison of the results from small particles and bulk samples shows that conclusions obtained from small particle experiments are reliable and universal.     

Nonlinear limit for a system of diffusing particles which alternate between two states

We study a system of particles and the nonlinear McKean-Vlasov diffusion that is its limit for weak interactions. Each particle switches between two states, both with their own diffusion dynamics. There is interaction, in particular, in the rates of the switches. We show existence and uniqueness for the system of particles by stopping-time techniques. For the nonlinear martingale problem, we use a time-change that allows us to return to a strong pathwise representation, and then we use a contraction argument for an appropriate metric. Finally, we show propagation of chaos.     

Particle transport in a bottom-feed separation vessel

A two-dimensional, axisymmetric numerical model of particle separation in a bottom-feed separation vessel is presented. The model includes six separate particle classes and assumes that the settling velocity of each particle class is sufficiently small when compared to the high inflow turbulence levels that the effect of the particles on turbulence can be neglected. Low particle settling velocities coupled with low particle volume fractions allows application of a drift-flux multi-phase model. The comparison between numerical results and measured plant data is in good agreement for overflow of all particle classes. Results of simulations show that bottom feeding results in a negatively buoyant, particle-laden jet being formed in the core of the vessel. The fraction of large particles that is carried out through the overflow is found to be critically dependent on the inlet velocity. The most effective way to reduce carry-over of large particles at the same time as maintaining through-put is to increase the diameter of the inlet feed pipe.