中等Péclet数具势悬浮系统的沉降.pdfVIP

20卷 化工冶金增刊 (1999) 537 中等Peclet数具势悬浮系统的沉降* 王浩 温景篙 (南开人学物理科学学院 天津 300071) 摘 要 尝试解决中等Peclet数情7X下悬浮系统中受势力影响的粒 子沉降问题。采用流守恒 式的离散化方法将对分布方程离散化，分析了稳定系统的对分布函数的近场和远场行为，然 后通过迭代法得到了对分布函数的分布并计算了具势系统的沉降系数。 关键词 Peclet数，沉降系数，粒子间势 1 前 言 粒子沉降在理论上和实际应用上的重要性导致了许多相关的理论、实验工作的出 现。稀释悬浮系统的沉降问题。由Batchelor(,]以及 Batchelor和Wen[31的系列工作做了 系统的处理，并己经被广泛接受。但是两个难题的存在一直在阻碍对沉降问题的深入研 究:其 一是如何描述稳定系统的粒子间势，其二是具势的对分布方程的求解。我们先期 的系列工作提出了稳定悬浮系统中的 个势函数模型，解决了第 一个难题4[,5I第三个难 题的解决较第一个更为困难，我们先期的系列工作也仅处理介 些极端情况，例如单分 散、低Peclet数和高Peclet数等。据我们所知，关于中等Peclet数沉降问题的研究目前 还不存在。本文旨在尝试数值求解中等 Peclet数情况下的对分布方程，以求得到对分布 函数的分布，计算体系的沉降系数，从而与己发表的研究结果作为一个整体，给稳定系 统中的粒子沉降以一个全面的结论 2 对分布方程的离散化及边界条件的处理 稀释稳定悬浮系统中的粒子分布由对分布函数P。描述，决定P;的方程即对分布方 程，无因次化后它可以表示为石 V. f 及 中_飞 众 Pe0(.:,Pr)一Vv.·i弓ppi。,·..n,1V·(了1)一v·(下斋 V·Ar)=0. (I) 、畔 ‘’“ {‘“理 、kT 、衅」 由于是变系数的Fokker-Planck型方程，一般情况 卜中〔等Peclet数)对分布方程的求解 极为困难，目前只在一些极端条件有些结果[6J。处理中等Peclet数条件下对分布方程的 求解问题，用差分法离散化方程进行数值解是一个可选择的方案。如何离散化方程有很 多种选择，考虑到问题本身的一个重要特点是稳定性 (即J粒子在跨越以i粒子中心为 半径的球面上的通量为零)，我们选择了守恒型差分格式以确保这一特性在求解过程中 国家自然科学基金资助项目(项目编号 王浩:男，26岁，博十研究生，理论物理专业 王浩等:中等Peclet数具势悬浮系统的沉降 不被改变。为清楚起见，我们用图1来说明守恒型的差分方法。我们知道，图1中小体 积元Y内对分布的变化是由重力对流、势力和布朗扩散三种输送作用造成的，其中社向 和法向上的输送分量分别为 _Gd(cp,lkT) _OP，、 (2) p 一。一又~ 片5 {一PeLcosOp,,ds “ ] P(eMsin9p。一Hc豁、 (3) 因为无源无汇，对分布的变化纯粹由边界面上的输送造成，所以以上两分量在四个边界 上输送通量的和应该为零·将其中的app/as和中plae用中心差分代替后，即可将对 分布方程离散化成五点格式的差分形式，并可以用迭代方法求解。 具体求解离散化后的对分布方程时， 还需要找出对分布函数的内外边界值进行 迭代。对内外边界问题处理如下: (1)外边界 在远场，粒子间对流相互作用肯定超 过势力相互作用和布朗效应而居于主导地 位 作为一级近似，我们取Batchelor和 Wen在无穷高Peclet数情况下得到的对分 布方程的解析解9作为外边界条件121 (2)内边