直径任意分布球填充的数值模拟

提出球填充数值算法的新分类方法.改进原有的松弛算法,使其能够模拟直径任意分布的球填充问题,采用可变循环周期使不同球数情形下的填充率基本保持不变.算例数据表明,该算法的填充率和配位数均高于原算法.由于采用背景网格搜索和双向链表组数据结构,使得邻接球搜索效率有相当大的提高,算法的时间复杂度为O(N)(N为球数).在一台AMD Athlon 3200+PC上,对于10000个等径球的随机密排列,只需217s,填充率即可达到0.64.     

气体运动论统一算法在跨流域转动非平衡效应模拟中的应用

从考虑转动松弛变化特性的Rykov模型出发,基于Boltzmann模型方程求解跨流域气体运动论统一算法原理与计算规则,采用转动惯量描述气体分子自旋运动,研究考虑转动非平衡影响的Boltzmann模型方程数值求解方法.通过对氮气激波结构、二维钝头体和三维尖双锥跨流域绕流的模拟分析,验证该算法的跨流域一致适用性.     

气固鼓泡床的全尺度直接数值模拟

颗粒全尺度直接数值模拟由于不引入任何相间作用力封闭模型、同时又考虑了颗粒相与流体相之间的四向耦合,因而具有极高的计算精度和准度。本文采用内嵌边界多重直接力算法结合软球碰撞模型对实验室尺度的气固鼓泡流化床装置进行了全尺度直接数值模拟。计算结果真实还原了鼓泡流态化的流动情形,且能有效地捕捉装置内部流场的详细运动情况和涡结构。对流化床内颗粒在传统的离散元模型框架下进行曳力统计发现,平均曳力比全尺度模拟结果小约20%~30%。其具体数值因选取比较的曳力模型和统计网格而异。改进传统曳力模型需考虑颗粒群的非均匀性以及颗粒拟温度。     

同材质颗粒不同填充密度的随机介质中光场的空间分布

用时域有限差分法研究了同一材质的颗粒在不同填充密度下的随机介质中光局域化问题. 依据随机介质激光的实验参数,模拟了颗粒填充密度不同的随机介质中光场的空间分布. 结果表明：当散射颗粒的散射平均自由程与波长相当时,随机介质中的光场分布呈现局域化的特征,而且随着颗粒填充密度增大,光场的局域化程度增强. 因此,在同样的抽运激励下,颗粒填充密度越大的随机增益介质越容易产生激光辐射. 数值模拟结果与实验定性符合. 关键词： 随机激光 时域有限差分法 颗粒填充密度 局域化     

微靶装配中靶球几何量的高精度测量

为满足微靶装配过程中靶球和靶腔相对位置的严格的要求,保障惯性约束聚变实验的成功,针对机器视觉的几何量测量精度很难达到5 m以下的不足,设计了利用具有亚m级精度的激光共焦测头对装配中的靶球进行采点扫描以获取其尺寸和位置的非接触测量方法,借助于在线监测系统运动轴的高精度定位,通过上下两个激光共焦测头的组合测量,并采用一种快速而准确的基于点阵插值细分和目标两级提取的测量算法,使在线监测系统的球体几何量测量精度达到了2 m。     

含振动能激发Boltzmann模型方程气体动理论统一算法验证与分析

为模拟研究高温高马赫数下多原子气体内能激发对跨流域非平衡流动的影响,将转动能、振动能分别作为气体分子速度分布函数的自变量,把转动能和振动能处理为连续分布的能量模式,将Boltzmann方程的碰撞项分解成弹性碰撞项和非弹性碰撞项,同时将非弹性碰撞按一定松弛速率分解为平动-转动能松弛过程和平动-转动-振动能松弛过程,构造了一类考虑振动能激发的Boltzmann模型方程,并证明了其守恒性和H定理.基于内部能量变量对分布函数无穷积分,引入三个约化速度分布函数,得到一组考虑振动能激发的约化速度分布函数控制方程组,使用离散速度坐标法,基于LU-SGS隐式格式和有限体积法求解离散速度分布函数,建立含振动能激发的气体动理论统一算法.通过开展高稀薄流到连续流圆柱绕流问题统一算法与直接模拟蒙特卡罗法模拟结果对比分析,特别是过渡流区平动、转动、振动非平衡效应对绕流流场与物面力热特性的影响机制,证实了所建立的含振动能激发的Boltzmann模型方程及气体动理论统一算法的准确可靠性.     

一种光突发交换网络中降低填充开销的突发组装算法(英文)

针对光突发交换(OBS)网络现有组装算法在输入业务量较小时组装填充开销较大的性能缺陷,提出了一种新的突发组装算法,称之为BPRA算法.该组装算法通过改变突发生成时间,在不增加端到端时延的条件下利用偏置时间来延长突发组装时间,从而可以有效降低填充空闲比特的数目.同时提出一个近似分析模型来计算空闲比特填充概率.分析和仿真结果表明BPRA算法可以明显降低组装填充开销.     

高填充型复合材料的单胞模型有限元分析

通过有限元方法分析计算了高填充型复合材料面心立方和紧密堆积单胞模型的等效热导率。模拟研究了填料的填充量、填料的热导率和填料间接触面积的相对大小对复合材料等效热导率的影响。结果发现高填充量的面心立方复合材料可以获得较高等效热导率,大小颗粒混合填充的紧密堆积复合材料能够更大幅度地提高等效热导率。并且当填充颗粒之间出现接触时,等效热导率会随着填料间接触面积的增大而快速地增大。模拟研究表明,通过增加填料含量以及使填料之间相互接触形成热流通路是获得高热导率复合材料的有效途径,同时有助于指导相关的实验研究从而生产出高热导率的热界面材料。     

基于第一性原理分子动力学的填充方钴矿热输运性质及微观过程的研究

对于重要热电材料之一的填充方钴矿材料,其低热导率的成因存在两种观点:1)填充原子的局域振动引起共振散射降低热导率;2)填充原子的引入加强了三声子倒逆过程来降低热导率.本文采用含有限温度效应的第一性原理分子动力学方法模拟了YbFe_4Sb_(12)的动力学过程,并通过温度相关有效势场方法得到了充分包含非线性作用的等效非谐力常数,研究了微扰近似下的声子输运性质.结果显示,在填充原子振动全部参与三声子倒逆散射过程的近似下,相比于纯方钴矿体系,声子寿命大幅地降低,填充原子的振动是热阻的重要来源.但即便如此,理论计算结果与实验的晶格热导率之间仍存在明显偏离.不同填充原子振动之间的较弱关联性质也揭示其明显偏离经典的声子图像,表现为一种强烈的局域特征振动模式,并以此散射其他晶格声子,因而对热阻的贡献也超出了传统三声子的理论框架.通过将填充原子Yb振动模式的寿命进行共振散射形式的修正,可以使晶格热导率与实验结果符合较好.以上结果表明,YbFe_4Sb_(12)的低晶格热导率是由声子间相互作用以及具有局域振动特征的共振散射两方面因素导致.     

二维无摩擦颗粒体系中阻塞态的压缩过程

数值模拟无摩擦圆盘形颗粒体系的压缩过程,对单分散、双分散、线性作用、赫兹作用等进行对比.数值模拟表明:单分散体系容易晶格化,其临界填充率φC与临界平均接触数ZC都大于双分散的情形.以临界阻塞态为基准,阻塞态在压缩过程中的平均接触数增量ΔZ以及平均重叠量δ随填充率增量Δφ的变化与模拟条件几乎无关,压强P与Δφ的关系决定于作用力的形式,它们都能用幂律关系进行很好的拟合.     

Numerical Simulation of Random Close Packings in Particle Deformation from Spheres to Cubes

Variation of packing density in particle deforming from spheres to cubes is studied. A new model is presented to describe particle deformation between different particle shapes. Deformation is simulated by relative motion of component spheres in the sphere assembly model of a particle. Random close packings of particles in deformation form spheres to cubes are simulated with an improved relaxation algorithm. Packings in both 2D and 3D cases are simulated. With the simulations, we find that the packing density increases while the particle sphericity decreases in the deformation. Spheres and cubes give the minimum (0.6404) and maximum (0.7755) of packing density in the deformation respectively. In each deforming step, packings starting from a random configuration and from the final packing of last deforming step are both simulated. The packing density in the latter case is larger than the former in two dimensions, but is smaller in three dimensions. The deformation model can be applied to other particle shapes as well.     

Using level sets for creating virtual random packs of non-spherical convex shapes

Random packs of spheres have been used to model heterogeneous and porous material morphologies during simulations of physical processes such as burning of coal char, convective burning in porous explosives, and regression of solid rocket propellant. Sphere packs have also been used to predict thermo-mechanical properties, permeability, packing density, and dissolution characteristics of various materials. In this work, we have extended the Lubachevsky–Stillinger (LS) sphere packing algorithm to create polydisperse packs of non-spherical shapes for modeling heterogeneity in complex energetic materials such as HMX and pressed gun propellants. In the method, we represent the various particle shapes using level sets. The LS framework requires estimates of inter-particle collision times, and we predict these times by numerically solving a minimization problem. We have obtained results for dense random packs of various convex shapes such as cylinders, spherocylinders, and polyhedra, and we show results with these various particles packed together in a single pack to high packing fraction.     

形状渐变的颗粒与沉降液的相互作用

工业过程中存在大量的含非球形颗粒的多相流动,研究这类颗粒与流体的相互作用有重要的应用价值。本文选择底面为方形、形状随边长和高度规律变化的长方体颗粒,采用单颗粒沉降法测量了曳力系数,研究了形状变化、材料密度和沉降液黏度的影响。发现曳力系数在整个颗粒Re_p数范围(0.05～6300)内均高于球形颗粒,特别是采用阿基米德数Ar与Re_p数可揭示出所有数据的内在规律,提出了适用于这类长方体颗粒所有形状的统一曳力模型。     

Preliminary results of a non-spherical aerosol method for the retrieval of the atmospheric aerosol optical properties

There is experimental evidence that the non-sphericity of certain atmospheric particles can cause scattering properties different from those predicted by standard Mie theory. Numerous studies indicate the need to consider the presence of non-spherical particles in modeling the optical properties of atmospheric aerosols. On the other hand, natural aerosols show a great variety of shapes, making difficult a realistic choice of a particle shape (or shape mixture) model. In this paper, we test a parameterization of the particle shape in the retrieval of size distribution, phase function, single scattering albedo and asymmetry parameter from direct and sky-radiance measurements. For this purpose we have substituted the Kernel based on the Mie theory included in the model SKYRAD.PACK by one derived for non-spherical particles. The method is applied under different atmospheric conditions, including Saharan dust outbreak, polluted and local mineral episodes. We compare the results with those obtained by the well known spheroids algorithm used in the AERONET network.     

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

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

Particle Characterization by Projected Area Determination

The projected areas of non-spherical particles do not represent an unambiguous particle characteristic. Depending on the orientation towards a constant observational direction, different projected areas result. The spectrum of all projected area values of a particle, if determined representatively, gives the probability with which a certain value is obtained by a single measurement. In this work, the frequency distributions of different examples of test objects were both calculated and measured. The objects were a cube, a rectangular parallelepiped and also three model agglomerates consisting of spheres of the same size. Instead of just one projected area, during each measuring procedure three projected areas from three orthogonal directions can be obtained. A mean value is then calculated to reduce the ambiguity of the particle characteristic and enhance the resolution. A suitable measurement set-up is introduced. The results of calculation and measurement are compared for observation from just one direction and also simultaneous observation from three directions. The frequency distributions of the equivalent diameters of the particle projected areas show a characteristic trend of the total curve with remarkable properties. The simultaneous measurement of three values from mutually orthogonal directions and their mean value calculation result in a much narrower distribution. In this case, a non-sphericity factor can additionally be calculated, whose frequency distribution contains information in a characteristic manner about the degree to which the particle shape differs from a sphere.     

基于Bloch球面搜索的量子粒子群优化算法

通过分析量子势阱粒子群优化算法的设计过程,提出一种基于Bloch球面搜索的量子粒子群优化算法.首先用基于Bloch球面描述的量子位描述粒子,用泡利矩阵建立旋转轴,用Delta势阱模型计算旋转角度,用量子位在Bloch球面上的绕轴旋转实现搜索.然后用Hadamard门实现粒子变异,以避免早熟收敛.这种旋转可使当前量子位沿着Bloch球面上的大圆逼近目标量子位,从而可加速优化进程.仿真结果表明,该算法的优化能力优于原算法.     

非球形椭球粒子参数变化对光偏振特性的影响

针对自然界中多数沙尘、烟煤粒子的非球形问题, 在球形粒子偏振特性的基础上, 进一步研究非球形椭球粒子的折射率、有效半径、粒子形状等参数变化对光偏振特性的影响, 采用基于T矩阵的非球形粒子仿真方法, 模拟非偏振光经椭球粒子传输后光的偏振特性及其与球形粒子间的差异, 并以实际沙尘、海洋、烟煤三种气溶胶粒子为例说明结果的正确性. 结果表明: 当折射率实部越小, 虚部越大时, 球形粒子与非球形粒子的偏振差异越不明显; 当粒子有效半径增加时, 球形粒子偏振度的变化比非球形粒子更为明显, 且最大值分别出现在散射角为150°和120°的位置; 当粒子形状不同时, 不同形状椭球及球形粒子的差异在散射角小于60° 时并不明显, 且当椭球粒子纵横比互为倒数时, 两种粒子的偏振特性近似相同. 通过以上分析可知, 在光传输过程中, 椭球粒子多数情况下无法被近似为球形粒子进行计算.     

Numerical Simulation of Random Close Packing with Tetrahedra

The densest packing of tetrahedra is still an unsolved problem. Numerical simulations of random close packing of tetrahedra are carried out with a sphere assembly model and improved relaxation algorithm. The packing density and average contact number obtained for random close packing of regular tetrahedra is 0.6817 and 7.21 respectively, while the values of spheres are 0.6435 and 5.95. The simulation demonstrates that tetrahedra can be randomly packed denser than spheres. Random close packings of tetrahedra with a range of height are simulated as well. We find that the regular tetrahedron might be the optimal shape which gives the highest packing density of tetrahedra.     

Arbitrary direction incident Gaussian beam scattering by multispheres

Based on spherical vector wave functions and their coordinate rotation theory,the field of a Gaussian beam in terms of the spherical vector wave functions in an arbitrary unparallel Cartesian coordinate system is expanded.The beam shape coefficient and its convergence property are discussed in detail.Scattering of an arbitrary direction Gaussian beam by multiple homogeneous isotropic spheres is investigated.The effects of beam waist width,sphere separation distance,sphere number,beam centre positioning,and incident angle for a Gaussian beam with two polarization modes incident on various shaped sphere clusters are numerically studied.Moreover,the scattering characteristics of two kinds of shaped red blood cells illuminated by an arbitrary direction incident Gaussian beam with two polarization modes are investigated.Our results are expected to provide useful insights into particle sizing and the measurement of the scattering characteristics of blood corpuscle particles with laser diagnostic techniques.