三维分形絮团沉降的格子Boltzmann模拟.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 三维分形絮团沉降的格子 模拟 张金凤 张庆河 林列 天津大学建筑工程学院 天津 南开大学信息技术科学学院 天津 摘要 黏性泥沙絮团的形成机理和沉降特性对河口海岸细颗粒泥沙运动规律研究具有重要作用 为了从微观结构 出发研究絮团运动机理 由扩散受限絮凝体聚集模型生成不同大小的分形絮团 引入格子 方法模拟三维 分形絮团的静水沉降 获得了絮团沉速的变化过程 比较格子 算法模拟结果与采用盒子计数分维数的 絮团沉速公式 发现二者有较好的一致性 表明格子 模拟方法可应用于河口海岸环境中的黏性 泥沙运动研究 关键词 黏性泥沙 分形絮团 格子 方法 静水沉降 中图分类号 文献标识码 絮凝现象是河口海岸细颗粒黏性泥沙研究的关键问题 分形理论应用到絮凝过程为定量描述絮团 的复杂结构形态和絮凝发育的随机聚集过程提供了全新方法 可以揭示絮团内在特性 利 用絮凝器研究了由流动剪切和不等速沉降形成的絮团分维特性 和 通过现场絮团的照 片 采用面积和周长的关系计算了絮团分维数 和 利用絮团分形的特点描述 絮团生长和破裂 同时 利用分维数建立了沉速与絮团大小的关系 和 总结了 现场水下摄影系统得到的实际河口海岸絮团特性 并利用 沉速表达式计算了絮团分 维数 等在研究河口泥沙絮凝时 考虑絮团内部颗粒排列的分形结构 提出了聚结分形球体模 型 但到目前为止 人们对分形絮团沉降特性的研究还主要是通过理论分析和现场或实验室观测进行 关于絮团形成 破裂乃至沉降的直接数值模拟的工作很少 其原因在于求解宏观方程的方法 如 方程 二相流方程等很难描述具有复杂微观结构的絮团运动 而格子 即 简称 模拟方法能够从微观尺度对颗粒的运动进行模拟 因此 利用 方法研究絮团沉降 及其相关问题有可能为黏性泥沙研究开辟新的途径 的开创性工作为悬浮颗粒运动的 模拟奠定了基础 等在 模型的基础上 模拟了不同形状的颗粒在各种边界条件下的悬浮 和 对粒子悬浮的 方法进行发展 引 入润滑 理论 成功地模拟了相互接触颗粒的沉降 吴锤结等 将各固体颗粒的牛顿力学和 方法相结合 研究了二维不规则形状悬浮颗粒的运动 模拟了静水中移动颗粒周围的流 场 和 采用插值方法处理边界条件 提高了颗粒问题的计算精度 可以用较少的网格实现对 颗粒运动的模拟 研究了二维分形物体受到的拖曳力及其影响因素 最近 还研究了链 状聚合物的水动力特性 然而 目前关于黏性泥沙三维絮团沉降的 模拟的工作尚未见到 为了探 索 方法在黏性细颗粒泥沙絮凝研究中的应用 本文将基于 模型分析三维分形絮团的沉降特性 收稿日期 基金项目 国家自然科学基金资助项目 作者简介 张金凤 女 河北人 博士生 主要从事海岸河口水动力及泥沙研究 方法的理论基础 方法是近十几年来发展起来的一种新的数值计算方法 传统流动问题的数值计算主要是直接 对以连续介质假定