颗粒流在沙坡上运动的直接数值模拟.pdf

颗粒流在沙坡上运动的直接数值模拟1 江召兵 吴锤结2 (解放军理工大学理学院流体力学研究中心，南京，211101) 摘要：颗粒流是自然界的一个很普遍的现象。本文运用直接数值模拟方 法，把颗粒表面看成流体运动的内边界，合并流体运动方程和颗粒运动方程，得 到统一的弱解形式的有限元方程。我们对一个颗粒和大量颗粒在沙坡上运动进行 了数值模拟，对颗粒的速度和角速度与沙坡的角度的正弦值的关系用误差分析和 相关性分析，得出了一些非常新颖的有意义的结果；对颗粒的阻力系数和运动轨 迹也进行了分析，加深了我们对颗粒流在沙坡上运动的认识和理解，并为以后的 工作打下了基础。 关键词：颗粒流 直接数值模拟误差 相关性 1 引言 颗粒流是自然界的一个很普遍的现象，如花粉的输运与传播，河床的生成与迁移；它 也是工业应用值得关注的问题，如冶金工业金属粉尘的对人体的危害，石油工业微小颗粒 的沉淀；它与人们的生活也是密切相关的，如近几年来频繁发生的的沙尘暴现象[t]等。 本文采用有限元动网格技术，运用并行计算模拟了颗粒流在沙坡上的运动，数值结果 加深了我们对其运动机理的认识，对其运动规律进行了误差分析和相关性分析，并提出了 一些有益的新颖的见解。 2直接数值模拟方法 的技术求解离散方程。 2．1 控制方程 考虑N个固体颗粒在不可压缩流体中的运动[2J。设流体所占的区域为Q0(￡)，第i个颗 粒(i=1，2，…，N)所占的区域为Qi(￡)，Q0(￡)和Q(￡)的边界分别为3Q0(￡)和3Qi(￡)。 控制流体运动的质量守恒方程为： V·击一o (1) 控制流体运动的动量方程为： 一I ‘ LZJ(2) IDf·餐一IDf于+一v·～LaJIDf’了F—p，十 1 国家自然基金重点项目资助课题 2 E—mail：ciwu@jlonline．corn 20 苎!，去为流体速度矢量，p为压力，lDf为流体密度，[仃]为应力张量。控制颗粒运动牛 顿第二定律为： 一学一邑+置一6一k，M·丢ds (3) 控制颗粒运动的转动定律为： 学=}i：一，呱㈤c王～瓤①¨瑚 (4) 别为颗粒的平动速度和转动速度；上Fi和手i分别为作用在颗粒上的力和力矩。 流体和颗粒的初始条件分别为： ■■…1 乱一蛳，L却_J2L铴]，into(O) (5) 和 xi(o)一xo，@i(o)一@?，玑l-(o)=Wr，cuf(o)=∞?，fori=1，2，…，N (6) 鉴于颗粒流运动系统的复杂性，有限元方法是一个强有力的工具，合并方程(2)、 (3)和(4)可得到弱解形式的有限元方程。引入函数空问寺，假设： f占一(乏，Ⅵr，巩r，∞。，州)}磊∈Hlr(‰)。1 y一{Ⅵ∈缈，肼∈缈 l磊一彘+雠x(三一是) j 彗皇’H1(fro)3为流体速度场空间，卵为颗粒速度场空间，空间古是流体颗粒混合物的 自然空间。为得到方程的弱解形式，选择试函数： 卣％一(丢％，Ⅵr，古N叫1，∞N)∈砜r 把方程(2)两边乘以流体速度试函数击％，并积分得： fn Pf[甓+(乏吲乏一砰乏％蝌-f。。M：[讹肭。