粘性细颗粒泥沙静水絮凝沉降生长的计算机模拟.pdfVIP

泥 沙 研 究 年 月 第 期 粘性细颗粒泥沙静水絮凝沉降 生长的计算机模拟 杨铁笙 李富根 梁朝皇 清华大学水利水电工程系河流海洋研究所 北京 摘要 利用分形生长模型对粘性细颗粒泥沙的絮凝沉降过程进行了计算机模拟 在扩散受限的絮团聚集模型 模型 的基础上 初步研究了粘性细颗粒泥沙的絮凝沉降分形生长过程 并对最终的絮凝体进行了图 像处理 得到了分形维数和颗粒浓度的关系 关键词 粘性细颗粒泥沙 絮凝沉降 分形维数 计算机模拟 中图分类号 文献标识码 文章编号 前言 在一定条件下 粘性细颗粒泥沙在水体中会凝聚在一起形成絮团 粘性细颗粒泥沙的絮凝沉降 和压缩密实问题是泥沙运动力学理论研究中的一个重要而复杂的课题 高含沙水流的许多特点 诸如 紊动的遏制 流核的出现 非牛顿流体的形成以及巨大的挟沙力等 都跟粘性细颗粒泥沙的絮凝有密不 可分的关系 另外 阵流现象 浆河现象和揭河底现象等 也和粘性细颗粒泥沙的絮凝有密切关系 目前对粘性细颗粒絮凝的研究主要还是停留在沉降特性 和影响因素的定性分析 阶段 另外 杨美卿 和 等人 对粘性泥沙和高岭土的絮凝体进行了电镜观测 对粘性细颗粒泥沙絮 凝体的结构有了初步的认识 但是 对于粘性细颗粒絮凝沉降缺乏观测手段 絮凝过程的研究遇到了很 大困难 随着分形理论的发展及普遍应用 使用计算机模拟成为研究絮凝体生长过程的有效辅助手段 可在 一定程度上弥补泥沙絮凝发育过程研究中的不足 粘性细颗粒泥沙的絮凝沉降也是一种分形生长现 象 年 和 提出扩散受限聚集模型 即 模型 模型假设有一个不动的中心微粒 这与实际情况相去甚远 此外 它只刻划了单个聚集体的生长 而没有涉及多个聚集体同时随机游走 相互碰撞的情况 后来 在 模型基础上 进一步发展出 模型 即扩散受限的絮团聚集模型 它取消了研究视场内 存在中心不动颗粒的限制 允许场内所有颗粒 聚集体同时作二维随机游走并互相结合 形成颗粒更多 尺度更大 内部结构更复杂的絮团 这与粘性细颗粒泥沙的实际絮凝沉降过程非常接近 文献 也指 出了可以对粘性细颗粒泥沙应用 模型进行絮凝过程模拟 粘性细颗粒泥沙絮凝沉降过程的计算机模拟 模拟的基本假设和简化 液体中分散介质在一定条件下 在各种原因 布朗运动 水流紊动 差速沉降等 引起的碰撞 接触 收稿日期 基金项目 国家自然科学基金资助课题 作者简介 杨铁笙 男 北京人 教授 主要从事河流动力学的教学与研究工作 中 由于微观作用力引起的颗粒粘结现象 称为絮凝 长程微观力 或凝聚 短程微观力 在本数值模拟 中 作如下简化假设 悬浮颗粒在静水中沉降 体系中的泥沙颗粒是大小均一的球体 颗粒的运动 由布朗运动和重力沉降运动所控制 尽管布朗运动的影响很小 颗粒由于相互接触碰撞而发生永久 性粘结 在数值模拟过程中 圆形颗粒是否发生粘结作用是根据其外切矩形是否有接触来判断的 因此 颗 粒的接触实际上是圆形颗粒外切矩形的相互接触 实际模拟时 圆形颗粒外切矩形的相互接触又可以 分为四种情况进行判断 如图 所示 图 颗粒粘结的几种假设情况 模拟的基本参数和公式选取 基本参数 泥沙颗粒质量密度 泥沙颗粒粒径 模拟时表示一个粒径的尺度为 个 象素 泥沙浓度 根据实验室的试验情况 浓度可选 采用人机对话方式输入 时间步长 按 照单颗粒沉降速度和模拟区域的空间尺度