粒间摩擦对含有点缺陷的二维颗粒堆积体底部力链的影响

力链对颗粒物质的宏观与微观力学性质起着决定性的作用。在离散元法的基础上,建立二维规则排列的颗粒物质系统,分别研究无缺陷颗粒系统在集中载荷变化与有缺陷颗粒系统在缺陷区域改变时,粒间摩擦系数对颗粒系统底部接触力分布规律的影响。结果表明:在无缺陷颗粒系统中,颗粒系统底部接触力的分布形式受摩擦系数和集中载荷的大小影响。随摩擦系数的增大,底部接触力由双峰形式经平台过渡,逐渐向单峰形式转变。在有缺陷颗粒系统中,摩擦系数和缺陷尺寸对底部接触力分布均有影响。同种载荷作用下,随缺陷尺寸的增大,底部接触力峰值显著增大;底部平均力被明显削弱,力向边界的传递增强。系统中轴线上缺陷的存在使底部中间区域受力削弱,当缺陷尺寸超过2层以上时,底部中间力随摩擦系数的变化特征由递增曲线演变为线形衰减曲线。     

EFFECT OF INTERGRANULAR FRICTION ON THE BOTTOM FORCE CHAIN OF TWO-DIMENSIONAL PARTICLE PILES WITH POINT DEFECTS1)

力链对颗粒物质的宏观与微观力学性质起着决定性的作用。在离散元法的基础上,建立二维规则排列的颗粒物质系统,分别研究无缺陷颗粒系统在集中载荷变化与有缺陷颗粒系统在缺陷区域改变时,粒间摩擦系数对颗粒系统底部接触力分布规律的影响。结果表明：在无缺陷颗粒系统中,颗粒系统底部接触力的分布形式受摩擦系数和集中载荷的大小影响。随摩擦系数的增大,底部接触力由双峰形式经平台过渡,逐渐向单峰形式转变。在有缺陷颗粒系统中,摩擦系数和缺陷尺寸对底部接触力分布均有影响。同种载荷作用下,随缺陷尺寸的增大,底部接触力峰值显著增大;底部平均力被明显削弱,力向边界的传递增强。系统中轴线上缺陷的存在使底部中间区域受力削弱,当缺陷尺寸超过2 层以上时,底部中间力随摩擦系数的变化特征由递增曲线演变为线形衰减曲线。     

颗粒物质在黏性流体中运动过程的数值分析

为了研究颗粒物质在黏性流体中整个下落过程的运动规律，对不同半径的颗粒物体运动过程进行了数值模拟分析。在给出颗粒物质下落过程模型的基础上，考虑黏性流体对颗粒物质的黏滞阻力，并通过受力分析建立了下落运动微分方程；利用计算机数值分析了不同半径的颗粒物质在同一种黏性流体中下落速度、入潜深度随时间变化的情况，并分析了不同半径的颗粒物质下落达到的稳态速度与所消耗时间的变化关系。计算结果表明：当颗粒物质半径 r<1mm 时，其在黏性流体中的下落距离与时间近似呈线性关系；随着颗粒物质尺寸的增加，其下落距离会呈现非线性增加；当颗粒物质半径 r>2mm 时，其达到稳态速度后的入潜深度与颗粒物质的半径呈非线性变化，且颗粒物质越大，达到稳态速度后的入潜深度越深；颗粒物质在黏性流体中下落后达到的稳态速度与颗粒物质的半径并非呈线性关系，且颗粒物质越小，达到的稳态速度越小，与黏性流体的深度无关，其相应入潜深度与所需时间呈非线性关系。     

飞行器船尾设计对亚声速底阻影响及风洞试验验证

减小阻力尤其是底部阻力,是飞行器增加航程的重要手段,而船尾布局则是减小底部阻力的有效措施.为研究船尾修型设计在亚声速段对飞行器阻力的影响规律,对某飞行器外形开展船尾修型设计并进行数值模拟,分析了不同船尾形状和船尾角度对飞行器阻力的影响情况,并开展了风洞试验验证.结果表明,船尾修型设计可以有效减小底部阻力;船尾采用曲线过...     

基于颗粒离散元的沙纹演变大涡模拟研究

应用大涡模拟、"点球" 浸入边界法和基于"事件驱动" 模型和"弹簧-阻尼" 模型的颗粒离散元法, 数值模拟了明渠湍流中沙纹的演变. 通过对不同谢尔兹数下无沙纹床面的推移质输沙率进行计算, 并与经典输沙率公式进行对比, 验证了模型的精度和可靠性. 随后, 对明渠湍流中沙纹床面的演变过程进行了数值模拟, 计算了有沙纹床面的推移质输沙率、沙纹长度和高度、等效床面高度的最大值、最小值和平均值、沙纹形状阻力、体积流速随时间的变化曲线. 研究发现:初始平整的床面在较短的时间内(tU_b=h≈100) 发展出数条沙纹, 随后沙纹逐渐发展, 在tU_b=h为1 600~2 000 时, 沙纹发生合并. 在沙纹数量不变的条件下, 沙纹高度随时间近似呈线性增长, 而沙纹的长度的平均值却保持恒定. 随着沙纹高度的增大, 输沙率和体积流速逐渐降低, 沙纹形状阻力则逐渐增大;当沙纹发生合并时, 沙纹高度快速增加, 输沙率、体积流速和沙纹形状阻力也出现了大幅跳跃. 在同等的水流强度条件下, 有沙纹床面的输沙率小于平整床面的输沙率.      

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象,该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径.为探究颗粒毛细效应的影响因素,采用离散元方法,模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性,考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律.结果表明,在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下,管径与粒径比D/d=3.33时,管内颗粒堵塞严重,使得颗粒上升缓慢,并造成颗粒柱中断; D/d=8.33时,起初毛细上升高度增加迅速,随后毛细上升高度的增大逐渐减缓,管内颗粒在管径方向几乎不存在速度梯度; D/d=15时,随着颗粒毛细上升高度的增大,管内颗粒柱分离为速度截然不同的两层,上层颗粒在管径方向几乎不存在速度梯度,而下层颗粒存在明显的速度梯度.研究还发现,在毛细效应能够发生的管径范围内,存在一个对应于颗粒最终毛细上升高度最大值的临界管径,当管径小于临界管径时,颗粒最终毛细上升高度随管径的增大而增大,当管径大于临界管径时,颗粒最终毛细上升高度随管径的增大而趋于减小;增大容器宽度,临界管径有所增大;增大振幅、适当提高频率能够有效促进临界管径的增大.     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象, 该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径. 为探究颗粒毛细效应的影响因素, 采用离散元方法, 模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性, 考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律. 结果表明, 在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下, 管径与粒径比$D/d = 3.33$时, 管内颗粒堵塞严重, 使得颗粒上升缓慢,并造成颗粒柱中断; $D/d = 8.33$时, 起初毛细上升高度增加迅速, 随后毛细上升高度的增大逐渐减缓, 管内颗粒在管径方向几乎不存在速度梯度; $D/d =15$时, 随着颗粒毛细上升高度的增大, 管内颗粒柱分离为速度截然不同的两层, 上层颗粒在管径方向几乎不存在速度梯度, 而下层颗粒存在明显的速度梯度.研究还发现, 在毛细效应能够发生的管径范围内, 存在一个对应于颗粒最终毛细上升高度最大值的临界管径, 当管径小于临界管径时, 颗粒最终毛细上升高度随管径的增大而增大, 当管径大于临界管径时, 颗粒最终毛细上升高度随管径的增大而趋于减小; 增大容器宽度,临界管径有所增大; 增大振幅、适当提高频率能够有效促进临界管径的增大.      

考虑墙趾阻力的加筋土挡墙的数值分析与筋材拉力计算方法

现行加筋土挡墙设计规范由于忽略了墙趾阻力对挡墙内部稳定性的贡献,计算的筋材拉力值远大于实测值,在设计上过于保守．本文对55座假定的模块式加筋土挡墙开展数值模拟,研究工作应力条件下筋材最大拉力随墙趾阻力的变化规律及其影响因素,在此基础上,采用响应面模型提出考虑墙趾阻力的筋材拉力计算方法,并用现场挡墙实测数据对该方法进行验证．结果表明,随着墙趾阻力增大,挡墙下部1/3高度范围内的筋材最大拉力值减小;墙趾阻力对筋材拉力的影响随着挡墙高度和填土内摩擦角的增加而减弱,随墙面仰角的增加而增强,受筋材刚度和间距影响不大;本文提出的方法可以合理预测加筋土挡墙在工作应力状态下的筋材拉力,并能反映出墙趾阻力使挡墙底部筋材拉力减小的现象．     

竖直振动管中颗粒毛细效应的离散元模拟

将一根细管插入填充有颗粒的静止容器中并对管施加竖直振动,颗粒将在管内发生上升运动,并最终稳定在一定高度,这一现象与液体毛细效应类似,被称为颗粒毛细效应.为探究颗粒毛细效应过程中伴随的颗粒尺度动力学行为及机理,基于离散元方法建立颗粒运动模型,对颗粒毛细效应动力学过程和特性开展数值模拟研究.模拟再现了文献中实验得到的颗粒毛细效应全过程,给出了管内颗粒柱高度随时间的演变规律,结果表明,受到颗粒系统参数的影响,本模拟条件下颗粒毛细效应过程呈现单周期上升、倍周期上升和倍周期稳定三个阶段,在倍周期上升阶段颗粒柱上升速度逐渐减小,平缓过渡到稳定阶段.在此基础上,分析了管内颗粒速度场和填充率分布随时间的演变特性,揭示了颗粒毛细效应过程中由容器传输到管内的颗粒的占比分布.研究发现,管内不同高度位置颗粒的运动并不同步,随着管的振动,管内出现速度波,速度波的传播引起管内颗粒出现膨胀和压缩交替的情况,从而管内颗粒填充率随时间发生周期性波动;在上升阶段,越接近管壁由容器传输到管内的颗粒占比越大,在稳定阶段,管内上层颗粒的对流引起容器传输到管内的颗粒占比发生反转.      

恒定大载荷划痕试验下紫铜的三维形貌及划痕硬度分析

采用圆锥形压头对紫铜进行划痕试验,并通过三维表面形貌仪获取划痕的三维形貌,研究正压力的变化对划痕沟槽所产生的影响. 结果表明:正压力的增大,使得划痕宽度和深度均线性增加,当正压力较大时,位错墙的形成使划痕深度出现周期性的波动,同时压头划刻过程伴有划痕两侧和前端的材料堆积现象,前端堆积高度和厚度、两侧堆积高度和宽度随着正压力的增加而线性增大. 通过“切削与塑性比”说明了压头对紫铜的刻划存在微犁耕和微切削两种变形机制,并且微切削机制在划刻过程中占主导地位. 磨损率随着载荷增加而线性增大,但划痕硬度不随载荷的变化而改变,约为0.77 GPa.      

Force and flow characteristics of an intruder immersed in 3D dense granular matter

The motion of a projectile impact onto a granular target results in both the resistance force exerted on the projectile and rheology of granular media. A horizontal arrangement of cylinder quasistatically and dynamically intruding into granular media under different velocities and angles is simulated using discrete element method. Three distinguished drag force regimes are exhibited, including hydrostatic-like force independent of velocity, viscous force related to velocity, and inertial drag force proportional to the square of velocity. Meanwhile, the influence of penetration angles on drag force is examined for these three regimes, and a force model, which is related to penetration depth and angle, is proposed for quasi-static penetration. Then, flow characteristics of the granular media, such as velocity field, pressure field, packing fraction etc., are traced, and a rheology model of packing fraction and inertial number is established.     

Discrete particle modeling of granular temperature distribution in a bubbling fluidized bed

The discrete hard sphere particle model (DPM) is applied in this work to study numerically the distributions of particle and bubble granular temperatures in a bubbling fluidized bed. The dimensions of the bed and other parameters are set to correspond to those of Müller et al. (2008). Various drag models and operational parameters are investigated to find their influence on particle and bubble granular temperatures. Various inlet superficial gas velocities are used in this work to obtain their effect on flow characteristics. It is found that the superficial gas velocity has the most important effect on granular temperatures including bubble granular temperature, particle translational granular temperature and particle rotational granular temperature. The drag force model affects more seriously the large scale variables such as the bubble granular temperature. Restitution coefficient influences all granular temperatures to some degree. Simulation results are compared with experimental results by Müller et al. (2008) showing reasonable agreement.     

Penetration into a granular bed in the presence of upward gas flows

We present a numerical study on the penetration of spherical projectiles into a granular bed in the presence of upward gas flows. Due to the presence of interstitial fluid, the force chains between particles in the granular bed are weakened significantly, and this distinguishes the penetration behavior from that in the absence of fluid. An interesting phenomenon, namely granular jet, is observed during the penetration, and the mechanism for its formation and growth is attributed to the merging of granular vortices generated by the interaction between the intruder and primary particles. Moreover, both the final penetration depth and the maximum diameter of the crater are found to follow a power-law dependence with the impact velocity, and the maximum height reached by the granular jet tends to increase linearly as the impact velocity increases, agreeing well with the experimental results reported in the literature.     

Investigation on effect of drag models on flow behavior of power-law fluid–solid two-phase flow in fluidized bed

In this study, a Eulerian-Eulerian two-fluid model combined with the kinetic theory of granular flow is adopted to simulate power-law fluid–solid two-phase flow in the fluidized bed. Two new power-law liquid–solid drag models are proposed based on the rheological equation of power-law fluid and pressure drop. One called model A is a modified drag model considering tortuosity of flow channel and ratio of the throat to pore, and the other called model B is a blending drag model combining drag coefficients of high and low particle concentrations. Predictions are compared with experimental data measured by Lali et al., where the computed porosities from model B are closer to the measured data than other models. Furthermore, the predicted pressure drop rises as liquid velocity increases, while it decreases with the increase of particle size. Simulation results indicate that the increases of consistency coefficient and flow behavior index lead to the decrease of drag coefficient, and particle concentration, granular temperature, granular pressure, and granular viscosity go down accordingly.     

Development of a filtered drag model considering effect of the solid shear rate

This study investigates a 2D gas–solid fast fluidized bed of typical Geldart A particles using highly resolved simulations with two-fluid model. The results show that the solid shear rate has a considerable impact on the orientation of the meso-scale structures and hence on the filtered drag force. On the basis of the correlation for the filtered drag force established in the literature using the traditional markers (such as filtered solid volume fraction, filtered slip velocity and filter scale), a correction correlated with the solid shear rate in the direction of gravity is proposed for better prediction of the filtered drag force. The corrected model is shown to produce improved results in posterior tests of flows in different fluidization regimes including bubbling, turbulent and fast fluidization.     

颗粒物质中的多尺度问题

颗粒物质是大量离散的固体颗粒相互作用而组成的复杂体系. 依据颗粒排布的稀疏程度, 体系可分为颗粒气体、颗粒流体和颗粒固体,它们有不同本质的动量传递和能量耗散机制. 后两者属于密集颗粒物质体系,内部形成了颗粒$\to $力链$ \to$体系的多尺度结 构,并涉及多个特征时间尺度,是典型的多尺度体系. 合理分割体系结构层 次、正确理解不同层次的物理过程、并确定它们之间的关联是密集颗粒物质研究的核心任务. 本文依次分析了密集颗粒物质的内在物理图像、多尺度结构层次和特征时间等,并介绍了多 尺度研究框架.     

Numerical modeling of fluid-particle interaction in granular media

Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method (CFD-DEM) approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier–Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.     

DEM simulation of polydisperse systems of particles in a fluidized bed

Numerical simulations based on three-dimensional discrete element model (DEM) are conducted for mono-disperse, binary and ternary systems of particles in a fluidized bed. Fluid drag force acting on each particle depending on its size and relative velocity is assigned. The drag coefficient corresponding to Ergun’s correlation is applied to the system of fluidized bed with particle size ratios of 1:1 for the mono-disperse system, 1:1.2, 1:1.4 and 1:2 for the binary system and 1:1.33:2 for the ternary system b...     

颗粒物质中最大静摩擦力随填充高度的变化规律

研究了圆棒在颗粒物质中受到的最大静摩擦力与颗粒填充高度的关系,并采用连续介质模型计算了圆棒在颗粒物质中受到的最大摩擦力.结果表明,当颗粒填充高度H小于圆筒直径的2倍时,最大静摩擦力F与H2成正比;当H大于2倍圆筒直径时,F与H成正比;理论计算结果与实验结果相符.