二维非恒定粘性泥石流运动堆积的数值模拟.pdfVIP

泥 沙 研 究 年 月 第 期 二维非恒定粘性泥石流运动堆积的数值模拟 胡健 匡尚富 徐永年 中国水利水电科学研究院 北京 摘要 以结构二相流模型为基础 建立了二维粘性泥石流运动堆积的数学模型 模型采用三角网格的有限差 分方法 具有计算效率高 边界处理精度高的优点 针对粘性泥石流运动速度快 龙头部位前缘坡度陡的特 性 简述了计算中的动边界处理及泥石流停积条件的处理方法 以云南东川蒋家沟泥石流堆积试验为例 模 拟了不同来流补给量下粘性泥石流的运动堆积过程 关键词 粘性泥石流 运动堆积 二维数值模拟 三角网格 中图分类号 文献标识码 文章编号 前言 泥石流运动模型大体可分为三类 一是宾汉体模型 这种模型以宾汉体的流变关系为基础 考虑粘 性力和剪切力的作用 以宾汉体本构关系模化阻力 用有限元 有限差分计算模拟泥石流的运动 适合于 描述泥流的运动 另一类是颗粒流模型 由日本的高桥堡提出 采用该模型对运动阻力进行概化处理 能 较好地模拟水石流的运动 第三类是结构二相流模型 将泥石流视为一种混合体 同时考虑细颗粒与水 组成的浆体和粗颗粒对流体运动的影响以及流体紊动作用 能描述各类泥石流的运动 由于两相流模 型方程中的固液两相间的相互作用力难以模化 因此在实际应用中还存在很大的困难 粘性泥石流中的固体颗粒含量较高 粒径分布范围广 流体中液相的浆体和粗颗粒固体都对泥石流 的运动产生影响 是一种浆体和固体粗颗粒组成的结构性二相流 针对这种泥石流的结构运动特性 近 似采用简化的单流体模型 模型采用有限节点法进行离散 这种方法是由赵士清教授 提出的 实质上是一种三角网格的有 限差分法 其思路上接近于积分形式的有限差分法 而在表现形式上则与有限元法更相似 泥石流的运动速度快 而且龙头前缘部分较陡 泥深变化大 计算时存在类似间断解的情况 要求算 法的计算时间空间步长小 精度高 并能处理动边界问题 该方法在离散过程中采用显示差分方法 不 需要求解大型方程组 计算效率高 计算时间短且动边界处理方便 并克服了一般差分方法边界模拟精 度差的缺点 能较正确地模拟复杂的天然边界 但对稳定性有较高的要求 实践证明 对于要求计算时 间空间步长小 存在运动边界而且龙头部分存在间断的粘性泥石流计算 该方法动边界处理方便而且计 算效率高 具有较强的适用性 基本方程 描述泥石流运动的结构二相流模型的简化单流体模型基本方程 为 收稿日期 作者简介 胡健 男 江西人 硕士 工程师 主要从事水力学及河流动力学研究 式中 为泥石流深度 分别为 方向的垂线平均流速 分别为 方向的床面坡降 为床面摩阻坡降在 方向的分量 为泥石流的综合粘度 为泥石流垂线平均密度 可 表达为 其中 为河床底面高程 床面摩阻坡降为液 固两相的摩阻坡降之和 为 式中 为泥石流的底面坡降 为宾汉应力 为泥石流固液混合物的容重 为泥石流浆体粘性系 数 为泥石流固相体积浓度 为固相颗粒容重 为液相的浆体容重 为固相摩擦角 为二维流速向量在 方向上的符号 方程 中阻力项由三部分构成 右侧第一项为粘性泥石 流液相浆体的宾汉应力引起的摩阻项 第二项与液相浆体粘性应力相关的摩阻坡降 等式右侧第三项来 源于粘性泥石流中的固相阻力 主要来自于固相颗粒的准静态摩擦 计算方法 将平面区域剖分为有限个三角单元的集合体 在每一个三角单元内 任意函数 可用三角单 元的节点坐标的线性函数来近似 其中 为三角形单元的面积 为所有以节点 为顶点的三角形单元的面积 系数 与单元节 点坐标有关 为 则三角单元节点 的一阶偏导数可表达为 二阶偏导数为 差分格式 将结构二相流的泥