移动边界的有限体积KFVS方法在一维溃坝波计算中的应用.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 移动边界的有限体积 方法在一维溃坝波计算中的应用 鲍远林 周晓阳 华中科技大学数学系 湖北武汉 摘要 本文从 方程出发 采用有限体积 方法 对微观粒子的速度分布进行 组合 给出了模拟河道洪水演进时干湿交界处移动边界的一种新方法 并成功地模拟了大坝下游为干河道的一维 溃坝模型 通过和精确解的比较得出 此方法在处理河道移动边界问题上是成功和有效的 关键词 方程 浅水方程 有限体积 方法 移动边界 中图分类号 文献标识码 在流体计算的研究中 流体与河岸交界线的移动边界的处理一直是难解问题 移动边界是水力计 算中常见的边界条件 如干河床上的瞬时溃坝波问题 干河床上的洪水演进问题 潮流涨落变化引起的 陆地动边界问题以及波浪爬高问题等等 目前处理这一问题最常用的方法是假设干地区域的初始条件 被一层水深很小的薄水层所覆盖 再利用变动的网格进行计算 这种处理方法虽然简单实用 但是 它要求计算域必须覆盖到洪水可能到达的所有区域 这在计算中会造成很多额外的计算网格和计算时 间 而且 对于某些问题如溃坝波问题 使用这种处理方法可能会产生完全不同的物理现象 文献 采用耗散有限元方法 并用网格变形技术及引入计算边界条件计算洪水演进前锋面 计算过程中需不断 调整计算网格 且新网格是通过旧网格上的信息插值而得 处理不当可能影响计算精度 文献 采用 变换技术给出拉格朗日坐标系下变动网格的一般处理方法 这些方法的共同特点是从描述连续介质宏 观状态的圣维南方程出发 采用某种离散化技术 得到方程的离散化近似解 为了研究分子系统的不平衡状态 在 年提出了著名的分子运动控制方程 后来称之 为 方程 其描述了分子速度分布函数的时空变化 近年来 以 方程为基础建立的 流体数值模型 在流体计算技术中受到很大关注 它本质上是追踪流体粒子 或粒子团 的运动规律 在模拟计算时引入实际粒子运动应遵循的物理规律作为限定条件 因此物理意义明确 构造方法简单 由于宏观物理量是大量微观粒子的相互作用导致的综合效应 因此宏观守衡律如 方程 都 可由 方程的矩来获得 显见 该方法的建模思路与经典方法有极大差异 等人最先利 用 形式的 方程建立了可压缩大气流动数值模型 邓家泉第一个应用 方程建立了适应水流 水质的数值模型 但尚未见到该方法应用于移动边界的文章 为此 本文将 从 方程出发 用有限体积 方法 开展对移动边界问题的计 算尝试 方程及 变量与宏观变量之间的关系 一维无碰撞 方程可写为如下形式 收稿日期 基金项目 国家自然科学基金资助项目 作者简介 鲍远林 男 安徽六安人 硕士生 研究方向为复杂系统建模与分析 式中 为粒子速度的 分布函数 为时间 为空间坐标 为粒子的随机速 度 为水深 为宏观平均流速 为作用在粒子上的合外力 为粒子质量 式 表明 粒子速度随时间的变化率 式 左边第一项 通过粒子对流 式 左边第二项 和作用 于粒子上的外力 式 左边第三项 进行改变 对于一维水流 粒子速度分布函数 与宏观变量水深 流速 等有以下基本关系 对式 求出粒子速度场的 阶和 阶矩 并由式 可得以下经典的一维浅水的连续方程 和动 量方程