溶解椭圆颗粒沉降的介观尺度数值模拟

在任意拉格朗日欧拉(ALE)算法模拟有热对流影响的颗粒两相流动的直接数值模拟基础上, 通过建立颗粒溶解速度和颗粒表面热流密度的关系, 对溶解的椭圆颗粒在垂直管道内牛顿流体中的沉降进行了直接数值模拟. 计算结果表明:与等温惰性椭圆颗粒沉降相比, 流体的对流运动、颗粒质量以及形状的变化等因素使溶解的椭圆颗粒在不同初始角度沉降时, 颗粒沉降动态尾迹、颗粒受力、颗粒沉降速度等都发了较大变化.     

热对流条件下固液密度比对颗粒沉降运动影响的直接数值模拟

本文应用任意拉格朗日-欧拉(ALE)算法对固液两相流流场中考虑热对流的非等温颗粒在竖直通道中的沉降运动进行了数值模拟.在牛顿流体中通过积分黏性应力和压力获得颗粒的受力跟踪颗粒运动,使用有限元方法数值求解流场的N-S方程和能量方程,模型不需经验假设.通过模拟来研究颗粒沉降的运动规律和热对流下固液密度比对固液两相流的影响作用.结果表明随着固液密度比的增加,颗粒经历了稳定沉降、周期性摆动,不规则摆动等过程;热对流使颗粒的摆动幅度和沉降速度发生变化;热对流对颗粒的影响作用随着固液密度比的增加而减小.     

微尺度稀薄效应下颗粒沉降的格子Boltzmann方法模拟

基于多松弛格子Boltzmann模型,对竖直细长微通道内颗粒自由沉降过程进行模拟,分析气体稀薄效应、初始位置以及颗粒间相互作用对微颗粒沉降特性的影响.研究表明：随Knudsen数增大,微通道内气体稀薄效应增强,颗粒表面气体滑移速度增大,气相流体有效粘度减小,颗粒相同运动状态下受到气体阻力相应减小,颗粒沉降平衡速度明显增大;不同初始位置颗粒沉降过程存在明显差异,初始位置偏离中心线颗粒将发生水平方向位移且呈振荡趋势,最终稳定于中心线平衡位置;在微尺度双颗粒沉降DKT现象过程中,气体稀薄效应影响颗粒运动特性,后颗粒跟随过程明显增长.     

溶解与热对流对固体颗粒运动影响的直接数值模拟

对牛顿流体内溶解与热对流对单颗粒在垂直管道中的沉降运动进行了直接数值模拟.流体运动由守恒方程计算,密度和黏性的变化考虑流场温度变化的影响,通过积分黏性应力和压力获得颗粒的受力跟踪颗粒运动,溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定.通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合.分别模拟了颗粒在等温流体、热流体、冷流体及颗粒溶解四种情况下的沉降过程.结果表明,在一定雷诺数内,热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒表面形态的变化引起了颗粒的横向摆动,并使颗粒沉降速度发生了变化. 关键词： 溶解 热对流 颗粒两相流 直接数值模拟     

热对流条件下颗粒沉降的格子Boltzmann方法模拟

采用格子Boltzmann方法模拟了在热对流条件下的颗粒沉降问题, 在研究单颗粒在等温流体、热流体和冷流体中运动的基础上, 进一步模拟了两个不同温度的颗粒在流体中的沉降.结果表明:两等温颗粒的沉降方式与雷诺数Re以及格拉晓夫数Gr密切相关, 而两不同温度的颗粒与两等温颗粒的沉降规律有显著不同.无论初始位置如何, 冷颗粒最终总位于热颗粒下方运动, Re较大时, 发生连续的拖曳、接触现象, 而Re较小时, 冷颗粒会以较大的沉降速度远离热颗粒. 关键词： 格子Boltzmann方法 颗粒沉降 热对流     

双颗粒在溶解条件下沉降的多相流动特性

对牛顿流体内溶解双颗粒在垂直管道中的沉降运动进行了直接数值模拟. 流体运动由守恒方程计算, 密度和黏性的变化考虑流场温度变化的影响, 通过积分粘性应力和压力获得颗粒的受力跟踪颗粒运动, 溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定. 通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合. 对双颗粒在等温流体无溶解条件和非等温流体溶解条件下的沉降过程进行了计算. 结果表明, 在一定雷诺数内, 热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒质量、颗粒表面形态的变化引起了颗粒的横向摆动, 并使颗粒沉降速度发生了变化.     

热对流条件下双颗粒沉降的直接数值模拟

在任意拉格朗日-欧拉(ALE)算法模拟等温惰性颗粒两相流的基础上,增加对能量方程的联立求解,对热对流条件下双颗粒的沉降进行了直接数值模拟.结果表明:双颗粒在等温流体、热流体和冷流体中运动形式是不同的,造成不同的主要因素是热对流引起了颗粒不同的动态尾迹,在热流体中颗粒尾部形成了涡的脱落,在冷流体中颗粒尾部形成了羽流.     

直接模拟中不同边界条件的实施及对沉降规律的影响

在牛顿流体中, 对颗粒在4种不同边界的垂直通道中的沉降运动进行了直接数值模拟. 计算结果表明:通过计算区域随颗粒运动而移动构建的无限长通道能准确模拟颗粒自由下落到稳定沉降的发展过程; 周期性边界条件由于流场变化, 对颗粒沉降产生了影响, 不能模拟颗粒的自由沉降过程; 底部封闭边界适合模拟封闭容器内颗粒与固壁的相互作用过程, 若颗粒达到稳定沉降, 也能模拟无限长通道内的沉降过程; 流化边界适合模拟流化床内气固两相流动. 计算结果有助于更好地理解和使用不同边界条件. 关键词： 直接数值模拟 边界条件 沉降 任意拉格朗日-欧拉方法     

粮仓内颗粒压力的测量:Janssen行为及其偏差

基于侧壁分为上下两部分的装置,测量了圆柱对称颗粒系统(竖直粮仓)的压力随高度的变化,并讨论了测得的仓内压力与Janssen理论模型的偏离情况.实验结果显示无论是否对粮仓进行Vanel等建议的让"侧壁摩擦力充分动员"的底座缓慢下沉操作处理,测得的仓内压力都有不同程度的偏离Janssen现象.另外如果没有沉降处理,在填充颗粒的荷载下上侧壁会发生微弱下沉,从而显著增大压力偏离的程度,但沉降处理可将其消除.这些结果表明,粮仓的颗粒重量在仓底和沿侧壁边界上的分配情况受制备方式和装置变形的影响明显,具体分布情况复杂多样.因此颗粒物质的静应力弹性理论分析应该取应力边界条件.     

嵌套椭圆时空管道的监控视频运动片段快速提取

针对现存运动片段提取方法大多步骤繁琐、计算量大,无法满足快速提取运动片段的现实需求的问题,提出一种嵌套椭圆时空管道的监控视频运动片段快速提取算法,以最大限度地减小计算量.该算法首先对监控视频进行椭圆时空采样形成时空管道,其次递进椭圆时空采样组成嵌套椭圆时空管道,然后展开嵌套椭圆时空管道生成时空平面图,最后对时空平面图进...     

不同振动模式下颗粒分离行为的数值模拟

利用三维离散元法对垂直方向上的直线、圆和椭圆振动模式颗粒分离过程进行了数值模拟研究,对直线振动时上层大颗粒的波动及圆和椭圆振动时出现的聚集、循环等现象的形成机理进行了分析,并讨论了振动强度对各振动模式下颗粒分离形态的影响规律. 研究表明,综合运用空隙填充、侧面驱动的颗粒运动和能量非均匀分布三种机理,并结合颗粒群的速度矢量分布情况能够较好地解释各振动模式下的颗粒分离行为. 振动强度对圆和椭圆振动模式的分离形态具有显著的影响,并在振动强度约为3时,各种振动模式均具有良好的颗粒分离效果和稳定的颗粒运动状态.     

柱状粒子沉降的数值模拟

通过对基于柱状粒子振荡频率ω而建立的粒子受力方程和力矩公式修正,得到粒子在静水中沉降的控制方程。用该方程求解长径比在0.1～10范围内沿轴向和垂直轴向运动的粒子的Stokes沉降速度,以及任意初始倾角高频振荡粒子沉降速度和倾角随时间的变化。所得结果与实验数据吻合较好。     

不同振动模式下颗粒分离行为的数值模拟

利用三维离散元法对垂直方向上的直线、圆和椭圆振动模式颗粒分离过程进行了数值模拟研究,对直线振动时上层大颗粒的波动及圆和椭圆振动时出现的聚集、循环等现象的形成机理进行了分析,并讨论了振动强度对各振动模式下颗粒分离形态的影响规律. 研究表明,综合运用空隙填充、侧面驱动的颗粒运动和能量非均匀分布三种机理,并结合颗粒群的速度矢量分布情况能够较好地解释各振动模式下的颗粒分离行为. 振动强度对圆和椭圆振动模式的分离形态具有显著的影响,并在振动强度约为3时,各种振动模式均具有良好的颗粒分离效果和稳定的颗粒运动状态. 关键词： 振动模式 颗粒分离 离散元法 数值模拟     

高温超导体中约瑟夫森涡旋流阻的振荡效应

通过数值计算耦合sine-Gordon方程组研究高温超导体中约瑟夫森涡旋的运动,得到约瑟夫森涡旋电压和流阻随平面磁场和驱动电流的变化规律.固定驱动电流,约瑟夫森涡旋电压和流阻随着磁场的增大出现周期性的振荡行为,振荡周期与每层约瑟夫森结中进入一个磁通量子相对应.分析和阐明了产生这种周期性振荡的原因. 关键词： 约瑟夫森涡旋 涡旋格子 高温超导     

莱顿弗罗斯特水滴振荡模式的影响因素及机理探究

探讨了莱顿弗罗斯特效应下水滴振荡的二维边界条件,在此基础上可求解液滴的本征频率.将振荡归结为由表面波叠加形成的驻波,并引入参量共振描述振荡的物理机制.探究了一定温度下振荡模式的演化规律,给出了振荡模式随半径的分布及振荡频率随时间的变化.验证了热对流为液滴与表面间主要的传热方式,并对真实液滴的几何形状做了定性探讨.     

液滴碰撞Janus颗粒球表面的行为特征

为研究液滴碰撞Janus颗粒(双亲性)球表面的独特行为特征,以粒径为5.0 mm铜球为材料制备了Janus颗粒,用直径为2.0 mm的液滴,在韦伯数(We)为2.7,10,20,30的测试情况下对Janus颗粒球表面进行了碰撞实验.结果表明:液滴碰撞Janus颗粒球表面后的运动可分为铺展、回缩、振荡和回弹4个过程.在不...     

不同条件下颗粒物质的粘弹性响应

研制了能承受100 MPa压强的测量装置,并用该装置对加压下制备的紧密沙堆进行粘弹性响应测量,得球形镍颗粒沙堆在周期性竖直振荡下储能模量G′与压强p 的关系. 实验结果表明,反映沙堆整体弹性特征的储能模G′与压强p 的关系为G′~p0.45;不同形状颗粒体系的粘弹性响应行为不同.     

神舟七号飞船小卫星伴飞原理分析

神舟七号伴飞小卫星和神舟七号飞船由于两个轨道的近地点和远地点不重合,导致它们在共同绕地球公转的同时,伴飞小卫星也在绕神舟七号飞船做周期性椭圆运动.伴飞周期等于飞船绕地球公转的周期,飞船位于椭圆轨道的几何中心.     

颗粒与顺列管束磨损的数值模拟

采用直接数值模拟方法模拟了在管道中直径80μm的颗粒对顺列10×10管束的碰撞和磨损.圆管与流场流体之间的相互作用采用内嵌边界方法进行计算并得到流场结构.颗粒的运动采用了拉格朗日跟踪方法计算并与流场之间进行了双向耦合.计算分析了沿流向以及垂直于流向各排管束受到颗粒的碰撞产生的磨损量.计算结果表明沿流向第三排管束以后的各排管束应采取防磨措旌.     

二维湍流热对流羽流运动路径对传热特性的影响

本文采用并行直接求解方法计算了0.7和4.3两个Pr数的二维湍流热对流,最高Ra=10~(13).高和极高Ra数的二维湍流热对流传热Nu数与反映羽流运动的大尺度环流路径随Ra数的变化具有很好的相关性,并具有两个Ra数转折点.第一转折点出现在大尺度环流由椭圆变为圆形,其周长C_(LSC)随Ra数变化突然减小,第二转折点出现在大尺度环流周长C_(LSC)最小值处,随后羽流变为旋涡团状而周长C_(LSC)随Ra数增高而增大.Pr数较小对应的转折Ra数较低.传热Nu数的变化规律用Ra~(0.3)补偿后发现, Nu的局部标度律在大尺度环流周长C_(LSC)较小的情况下减小,出现偏离GL理论预测倍数线的现象.当Ra数大于第二转折点后,极高Ra数二维湍流热对流Nu数随Ra数的变化与GL理论预测倍数线符合良好,走向保持一致.