非等温气固两相各向同性湍流的直接数值模拟

本文对非等温气固两相各向同性湍流中颗粒及颗粒所见流体温度的统计行为进行了直接数值模拟研究,系统地讨论了颗粒惯性的影响。对于Tp/Tk<1的颗粒,其温度脉动随惯性增大而减弱;对于Tp/Tk>1的颗粒,其趋势相反。颗粒及颗粒所见流体温度的自相关都随颗粒惯性的增加而减弱,而且随颗粒惯性的增加,颗粒所见流体温度的自相关比颗粒自身温度的自关联下降得快。     

颗粒分层过程三维离散元法模拟研究

采用软球干接触模型对球形及非球形颗粒的分层过程进行了三维离散元法模拟研究,从颗粒间作用力、转动力矩和能量变化的角度分析了颗粒分层机理,讨论了颗粒的粒度比对分层速度的影响规律.结果表明,分层过程中,大颗粒比小颗粒活跃,非球颗粒由于具有较高的动能而比球颗粒活跃,在一定程度上弥补了颗粒形状对分层过程的影响.大颗粒间的平均法向、切向作用力、平均力矩及平均动能均大于小颗粒.颗粒分层速度随着粒度比的增加而显著增大,当粒度比大于临界粒度比3时,分层速度的增幅减缓.     

椭圆截面非球形颗粒群的多重光散射

尽管非球形下一些特殊形状颗粒的单散射已被得到,球形多颗粒系（颗粒群）的多重散射也被研究,但至今仍未得到非球形颗粒群的多重散射.文中建立了一类椭圆截面非球形颗粒模型,求得其散射相位函数,借助于辐射传播方程,考虑形状及大小分布,得到了该类颗粒群的多重光散射.在两种特例情况下的结果能与已有的结果符合较好,说明了方法的可靠性.计算分析表明：非球形颗粒群的多重散射光强角分布要比球形颗粒平坦.椭圆截面颗粒的粒度或形状参数越大,多重散射光越集中于小的散射角;粒度分布或形状分布越宽,多重散射光强的角分布越平坦.随着光学厚 关键词： 多重光散射 颗粒 非球形 椭圆截面     

纳米一维双原子晶体颗粒的原子均方位移和均方速度

应用作者与合作者已导出的纳米一维双原子晶体颗粒的位移-位移格林函数,推导了纳米一维双原子晶体颗粒的速度-速度格林函数,并在此基础上,应用涨落-耗散定理,推导了纳米一维双原子晶体颗粒的原子均方位移公式和原子均方速度公式.数值计算表明,纳米晶体颗粒表面原子的均方位移大于内部的原子均方位移,而表面原子的均方速度小于内部原子的均方速度;当纳米颗粒尺寸减小时,表面原子的均方位移和内部原子的均方位移都减小,而表面原子的均方速度和内部原子的均方速度几乎不发生变化;温度增大时,纳米晶体颗粒中各原子的均方位移和均方速度都增大,在高温极限,各原子的均方位移和均方速度与温度成正比.     

基于超二次曲面的颗粒材料缓冲性能离散元分析

自然界或工业中普遍是由非球形颗粒组成的复杂体系,与球形颗粒相比,非球形颗粒间的高离散和咬合互锁可使冲击载荷引起的能量有效衰减实现缓冲作用.基于连续函数包络的超二次曲面单元能准确地描述非球形颗粒的几何形态,并可精确地计算单元间的接触碰撞作用.本文采用离散元方法对冲击载荷作用下非球形颗粒物质的缓冲性能进行数值分析,并与圆柱体冲击的理论结果和球体冲击的实验结果进行对比验证.在此基础之上,进一步研究了筒底作用力在不同颗粒层厚度和形状等因素影响下的变化规律.计算结果表明:不同颗粒形状都存在一个临界厚度H_c.当HH_c时,缓冲率随H的增加而增加;当HH_c时,缓冲率的变化不再显著并趋于稳定值.此外,减小颗粒表面尖锐度和增加或减小圆柱形和长方形颗粒的长宽比都会提高颗粒材料的缓冲效果.     

典型流化床内非球形颗粒流动行为模拟研究

气固流化床在能源、化工、制药、冶金等众多工业领域有着广泛的应用。本文通过耦合计算流体力学和离散单元模型的模拟来评估典型流化床中非球形颗粒的流动行为。结果表明，典型流化床内球形颗粒的旋转颗粒温度分布接近各向同性，非球形颗粒的旋转颗粒温度分布呈现明显的各向异性。随着颗粒长径比的增加，颗粒的平动颗粒温度和旋转颗粒温度增大，意味着颗粒平动和旋转速度的脉动加剧。     

基于有限元分析的高温气冷堆石墨颗粒碰壁过程研究

石墨颗粒碰壁过程的研究对定量计算高温气冷堆颗粒沉积率十分重要。本文基于有限元方法,利用ABAQUS软件对球形颗粒碰壁进行数值计算。仅考虑弹性变形时,计算结果与Hertz模型高度吻合.当考虑塑性变形导致的能量耗散时,计算得到的速度恢复系数变化趋势与实验数据吻合一致,但整体数值偏低,引入等效屈服强度概念后吻合较好.结果表明有限元方法数值计算颗粒碰壁准确可行,可为研究高温气冷堆石墨颗粒碰壁数值计算提供新的思路和方法。     

颗粒物质内自旋小球运动行为的数值模拟研究

针对大耳沙蜥在沙子中的运动行为,以球形物体为研究对象,利用3维离散元数值模拟程序LIGGGHTS模拟了旋转运动模式对颗粒物质中球形物体平动和上升/下降行为的影响,定量分析了旋转速度以及颗粒间摩擦系数等因素的具体影响.研究结果表明:球形物体与颗粒物质基底颗粒间的摩擦系数以及球形物体的自转角速度对球形物体的运动有明显影响,摩擦系数越大物体运动越明显,自转角速度越大物体运动越明显.该结果比较好地解释了沙漠生物外表具有鳞片的原因.     

用格子Boltzmann模型模拟垂直平板间的热对流

引入一个新的能量分布函数，利用该能量分布函数与粒子速度分布函数耦合来求解一个热流场. 因而，这一能量分布函数与粒子速度分布函数和Boltzmann方程构成了一个新的格子Boltzmann模型. 这一模型满足质量、动量和能量守恒的准则. 用该模型对垂直平板间的狭缝热对流进行了数值模拟，数值结果表明，在Prandtl数为1，Grashof数在1.3×10²—1×10⁶之间时，流场将出现多个旋涡结构的流型. 得出了与Lee相一致的结论. 关键词： 能量分布函数 Boltzmann方程 热对流     

火灾烟颗粒分形模型和球形模型光散射的比较研究

对烟颗粒的光散射进行模拟计算是研究火灾烟颗粒光散射特性的重要手段,目前对于火灾烟颗粒光散射的数值计算多采用球形或椭球模型.实际上,火灾烟颗粒的形貌与球形和椭球均存在着显著差异.扫描电子显微镜图像表明,烟颗粒具有近似分形的结构.本文利用离散偶极近似方法计算了随机取向的火灾烟颗粒分形凝团以及同体积的球形颗粒的光散射Muller矩阵,并对两者的归一化Muller矩阵元素随散射角的分布进行了比较.研究表明:火灾烟颗粒分形模型和球形模型的归一化矩阵元素F₁₁(θ)/     

Molecular rototranslation in condensed phases : Single particle theory

A multidimensional expansion of the Mori equation in terms of a chain of Markov equations is used to develop a theory of molecular rototranslation in condensed phases. The stochastic equations of motion are solved for transient and equilibrium averages of the relevant dynamical variables. The single particle rototranslational Langevin equations correspond to the first equation of the Markov chain and (with a rotational constraint) are solved using Wiener matrix algebra for a possible sixteen autocorrelation functions. The Einstein result for the mean-square velocity and angular velocity is generalized. The third dimension of the Markov chain corresponds mechanically to the (constrained) rototranslation of a molecule bound to a cage of nearest neighbours by a dissipative matrix γ. The cage is itself undergoing a rototranslational Brownian motion. The problem of evaluating the formal theory with experimental measurements is discussed in terms of the number of parameters associated with each approximant (or dimensionality of the Markov chain). It is possible to avoid using a least-mean-squares fitting procedure by using a broad enough range of data and simulator results.     

Dynamic Light Scattering Measurement of Nanometer Particles in Liquids

Dynamic light scattering (DLS) techniques for studying sizes and shapes of nanoparticles in liquids are reviewed. In photon correlation spectroscopy (PCS), the time fluctuations in the intensity of light scattered by the particle dispersion are monitored. For dilute dispersions of spherical nanoparticles, the decay rate of the time autocorrelation function of these intensity fluctuations is used to directly measure the particle translational diffusion coefficient, which is in turn related to the particle hydrodynamic radius. For a spherical particle, the hydrodynamic radius is essentially the same as the geometric particle radius (including any possible solvation layers). PCS is one of the most commonly used methods for measuring radii of submicron size particles in liquid dispersions. Depolarized Fabry-Perot interferometry (FPI) is a less common dynamic light scattering technique that is applicable to optically anisotropic nanoparticles. In FPI the frequency broadening of laser light scattered by the particles is analyzed. This broadening is proportional to the particle rotational diffusion coefficient, which is in turn related to the particle dimensions. The translational diffusion coefficient measured by PCS and the rotational diffusion coefficient measured by depolarized FPI may be combined to obtain the dimensions of non-spherical particles. DLS studies of liquid dispersions of nanometer-sized oligonucleotides in a water-based buffer are used as examples.     

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.     

周期场中非各态历经布朗运动

将自由状态下呈弹道扩散的非各态历经系统置于周期场中,进而将非各态历经布朗运动分为两类.第一类是阻尼核的Laplace变换的低频为零的系统,当温度远大于势垒高度时,系统平均能量的动能部分依赖粒子的初始速度分布;随温度降低,系统的各态历经性得到恢复.然后将第一类系统的稳定速度变量作为一个内部噪声,再去驱动一个自由布朗粒子,则阻尼核的Laplace变换在零频时为无穷大.结果发现,粒子扩散系数随温度的增加而趋于零,显示一种经典局域化特征,系统的渐进分布依赖于初始坐标分布.这是第二类非各态历经性运动,不能通过外加势而恢复. 关键词： 非各态历经 非Markov布朗运动 扩散系数 噪声谱     

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.     

非球体填充的组合球模型及松弛算法

现有的松弛算法由于仅用于球填充而只考虑颗粒的平动,故提出考虑非球体转动的改进松弛算法并采用组合球模型,使其能够模拟任意形状非球体的随机填充以及多种非球体的混合填充.用多个球体的外包络面近似一个非球体外形的组合球模型,将非球体之间的接触转化为球体之间的接触,从而简化并统一非球体接触判断算法.通过引入非球体的转矩和转角松弛机制,使改进松弛算法克服了"自锁"现象,并能生成非球体的随机密填充.算例表明,填充结果与现有的数值模拟及实验结果相符.     

Brownian dynamics simulations with spin Brownian motion on the negative magneto-rheological effect of a rod-like hematite particle suspension

We have investigated the behaviour of a suspension of magnetic rod-like hematite particles in a simple shear flow with the addition of an applied magnetic field. A significant feature of the present hematite particle suspension is the fact that the magnetic moment of the hematite particle lies normal to the particle-axis direction. From simulations, we have attempted to clarify the dependence of the negative magneto-rheological effect on the particle aggregation and orientational distribution of particles. The present Brownian dynamics method has a significant advantage in that it takes into account the spin rotational Brownian motion about the particle axis in addition to the ordinary translational and rotational Brownian motion. The net viscosity is decomposed into three components and discussed at a deeper level and in detail: these three viscosity components arise from (1) the torque due to the magnetic particle–field interaction, (2) the torque and (3) the force due to the interaction between particles. It is found that a slight change in the orientational distribution has a significant influence on the negative magneto-rheological effect. In a relatively dense suspension, the viscosity components arising from an applied magnetic field and the interaction between particles come to change rapidly for a certain strength of the magnetic particle–particle interaction, which is due to the onset of the formation of raft-like clusters.     

Lattice Boltzmann simulation of the two-dimensional flow of a magnetic suspension between two parallel walls under a non-uniform magnetic field

We have developed a lattice Boltzmann method based on fluctuation hydrodynamics that is applicable to the flow problem of a particle suspension. In this method, we have introduced the viscosity-modifying method, rather than the velocity-scaling method, in which a modified viscosity is used for generating random forces in lattice Boltzmann simulations. The viscosity-modifying method is found to be applicable to the simulation of a magnetic particle suspension. We have applied this method to the two-dimensional Poiseuille flow of a magnetic suspension between two parallel walls in order to investigate the behavior of magnetic particles in a non-uniform applied magnetic field. From the results of the snapshots, the pair correlation function between the magnetic pole and the magnetic particles and the averaged local particle velocity and magnetization distributions, it was observed that the behavior of the magnetic particles changes significantly depending upon which factor dominates the phenomenon in the balance between the magnetic particle–particle interaction, the non-uniform applied magnetic field and the translational and rotational Brownian motion.     

Spherical particle Brownian motion in viscous medium as non-Markovian random process

The Brownian motion of a spherical particle in an infinite medium is described by the conventional methods and integral transforms considering the entrainment of surrounding particles of the medium by the Brownian particle. It is demonstrated that fluctuations of the Brownian particle velocity represent a non-Markovian random process. The features of Brownian motion in short time intervals and in small displacements are considered.     

Mean-square displacement of a stochastic oscillator: Linear vs quadratic noise

Using the Langevin equations, we calculated the stationary second-order moment (mean-square displacement) of a stochastic harmonic oscillator subject to an additive random force (Brownian motion in a parabolic potential) and to different types of multiplicative noise (random frequency or random damping or random mass). The latter case describes Brownian motion with adhesion, where the particles of the surrounding medium may adhere to the oscillator for some random time after the collision. Since the mass of the Brownian particle is positive, one has to use quadratic (positive) noise. For all types of multiplicative noise considered, replacing linear noise by quadratic noise leads to an increase in stability.