离散颗粒层被横向推移过程中的力学行为研究

研究发现,颗粒物质层被匀速推移挤压过程中,所需推移力先以线性规律增加,在某一确定点后,则会以指数规律增加.而颗粒物质是由众多离散颗粒组成的软凝聚态物质,其宏观上反映的是离散颗粒的个体性质和凝聚态物质的集体效应.颗粒与颗粒之间以及颗粒与边界之间的细观尺度接触力链的构成以及演变规律将会直接影响各种宏观受力情况,其摩擦力与挤压力便是力链的主要构成形式.围绕着定量描述细观力链特征,从而揭示力的变化规律这一目标,采用计算机模拟的方法,依照球形颗粒Hertz法向接触理论和Mindlin-Deresiewicz切向接触理论,对重力作用下不同数目的三维等径球体颗粒层的推移情况进行了离散元仿真模拟,量化分析了推移力变化规律、各摩擦力变化规律以及力链分布规律,发现摩擦力与挤压力在颗粒堆积的不同阶段对力链的构成起到了不同的主次作用,使得力链发生强弱演变,从而发现了推移颗粒物质层时推移力的变化规律及原因.这些结果有利于从力链角度揭示颗粒内部和颗粒与各边界之间的受力情况.     

光诱导准一维光子晶格中离散空间光孤子的相互作用

基于分步束传播法数值分析了离散空间光孤子在准一维光诱 导光子晶格中的相干与非相干相互作用过程. 结果表明: 对于相干孤子, 同相时相互吸引, 反相时相互排斥. 然而, 由于非线性响应的各向异性, 横向排布的非相干孤子会因间隔波导数目的增加而由相互吸引变为相互排斥. 并且, 沿对角方向排布的两个非相干孤子在孤子相 互作用力和布拉格反射的共同影响下, 会呈现出"钟摆式"振荡传输现象. 研究结果有助于进一步理解非线性各向异性对离散孤子相互作用的影响机制, 并为后续实验研究提供理论参考.     

The effects of cubic potentials on discrete breathers in a mixed Klein－Gordon /Fermi－Pasta－Ulam chain

Nonlinearity has a crucial impact on the symmetry properties of dynamical systems. This paper studies a one-dimensional mixed Klein--Gordon/Fermi--Pasta--Ulam diatomic chain using the expanded rotating plane-wave approximation and numerical calculations to determine the effect of cubic potentials on the symmetry properties of discrete breathers in this system. The results will be very useful to researchers in the field of numerical calculations on discrete breathers.     

离散元与壳体有限元结合的多尺度方法及其应用

在深入研究复杂结构和非均质材料冲击响应和破坏机理的过程中,往往遇到多尺度计算问题。本文尝试建立三维离散元与壳体有限元结合的多尺度方法用于处理圆柱壳问题,该方法采用三维离散元对感兴趣的局域进行局部模拟,利用平板壳体有限元进行整体模拟,采用一种特殊的过渡层使离散元区和有限元区能很好的衔接。我们将这一方法应用于激光辐照下充压柱壳的热/力耦合冲击破坏响应,得到的模拟结果与文献报道有较好的吻合。     

二维多群辐射输运程序LARED-R-1的并行化

利用有向图描述数据依赖关系,应用已有的并行流水线通量扫描算法,实现基于非协调网格的二维辐射输运程序LARED-R-1的并行化.同时,采用消息缓冲技术提高并行程序的性能.经测试,对于典型的问题规模(100群、3800个网格单元、40个方向),在某并行机的64个和128个处理器上,并行程序分别获得80%和53%的并行效率.     

Influence of the tangential velocity on the compressible Kelvin−Helmholtz instability with nonequilibrium effects

Kelvin−Helmholtz (KH) instability is a fundamental fluid instability that widely exists in nature and engineering. To better understand the dynamic process of the KH instability, the influence of the tangential velocity on the compressible KH instability is investigated by using the discrete Boltzmann method based on the nonequilibrium statistical physics. Both hydrodynamic and thermodynamic nonequilibrium (TNE) effects are probed and analyzed. It is found that, on the whole, the global density gradients, the TNE strength and area firstly increase and decrease afterwards. Both the global density gradient and heat flux intensity in the vertical direction are almost constant in the initial stage before a vortex forms. Moreover, with the increase of the tangential velocity, the KH instability evolves faster, hence the global density gradients, the TNE strength and area increase in the initial stage and achieve their peak earlier, and their maxima are higher for a larger tangential velocity. Physically, there are several competitive mechanisms in the evolution of the KH instability. (i) The physical gradients increase and the TNE effects are strengthened as the interface is elongated. The local physical gradients decrease and the local TNE intensity is weakened on account of the dissipation and/or diffusion. (ii) The global heat flux intensity is promoted when the physical gradients increase. As the contact area expands, the heat exchange is enhanced and the global heat flux intensity increases. (iii) The global TNE intensity reduces with the decreasing of physical gradients and increase with the increasing of TNE area. (iv) The nonequilibrium area increases as the fluid interface is elongated and is widened because of the dissipation and/or diffusion.     

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model(DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation.It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model.Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model.To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified.     

一种易于硬件实现的图像有损压缩算法

原始图像数据经过小波变换后,根据子带系数之间在空间上和方向上所呈现出的带间相似性,将像素集合重新定义为集合树。通过对直流系数量化编码和交流系数比特深度编码中的参数K和熵编码中长度分别为3bit码字（choicel）和4bit码字（choice2）的编码方式的选择实现数据的渐进性编码,使产生的码流具有内嵌特性,可以对任意码流进行截取和译码,实现对码率的有效控制。提出了一种具有复杂度低、易于硬件实现的图像有损压缩算法。经过系统仿真证明此算法有利于提高重建图像的质量。     

回转窑多粒径颗粒流动与偏析的DEM模拟与实验研究

为明晰回转窑内颗粒的运动行为及偏析机理,以绿豆、黄豆和黑豆为颗粒介质,依次对3种装填顺序下的颗粒流动过程进行离散元模拟与实验研究,以颗粒质量分数和平均粒度为判据,对颗粒偏析进行评价。结果表明,回转窑内颗粒流动区可分为自由滚落区、渗流呆滞区以及窑壁携带区,自由滚落区颗粒流速最大,而渗流呆滞区流速最小。窑内颗粒沿轴向输运过程发生径向偏析,形成夹层结构,小颗粒受渗流作用在渗流呆滞区中心形成内核,大粒径和中等粒径颗粒集中在自由滚落区和窑壁携带区。窑内颗粒力链分布不均匀,强力链分布于近窑壁区,弱力链分布于自由滚落区和渗流呆滞区,且渗流呆滞区力链细而密集。当窑头附近不同粒径颗粒存在轴向速度差时,颗粒在轴向发生掺混,并产生径向偏析。     

混凝土搅拌车搅拌筒磨损数值分析

混凝土运输过程中搅拌筒的磨损一直是一种常见潜在危害,磨损严重时会导致叶片失效,对搅拌的质量和出料匀质性产生影响。通过实验获取搅拌筒内部的磨损费时费力,因此有必要采取一种数值分析方法对搅拌筒的磨损进行预测并提出改进。本文采用摩擦磨损实验的方法来标定颗粒与搅拌筒之间的Archard磨损常数,采用JKR接触模型表征混凝土的流动性能,采用离散元方法(DEM)对搅拌筒筒体及叶片磨损进行预测分析。通过法向接触能量与切向接触能量的对比,证明搅拌筒中的磨损主要为伴有冲击作用的磨粒磨损,搅拌车搅拌筒中搅拌叶片顶部的磨损较为严重。针对磨损比较严重的叶片顶部进行改进,采用T型耐磨结构等改进叶片结构,叶片顶部结构改进后搅拌筒使用寿命能得到明显提升。