无网格数值方法在结构物水动力研究中的应用 - 水科学进展.PDF

第 15 卷 第 6 期 水 科 学 进 展 Vol15 ,No6 2004 年 11 月 　　 ADVANCES IN WATER SCIENCE 　 Nov. , 2004 　 无网格数值方法在结构物水动力研究中的应用 李绍武 , 尹振军 (天津大学建筑工程学院港口工程系 , 天津 　300072) 摘要 : 运用一种移动粒子半隐式法建立无网格水流数值模型 , 探讨该模型在结构物水动力特性研究中的可能性 。在 无网格条件下 , 粒子间的相互作用通过核子函数来表示 , 并通过建立物理量光滑模型 、梯度模型和二阶导数模型来 离散控制方程 。流体的不可压缩性是通过引入粒子数密度概念实现的。该方法不受流体变形程度和固边界形状的 限制 , 因而具有特定的适应性 。运用该模型模拟了水柱倒塌后有无透空块体两种情况下右端直墙上的压力过程 。计 算结果表明透空块体可使直墙压力减小 42 % 。 关 　键 　词 : 无网格数值模型 ; 移动粒子半隐式法 ; 核子函数 ; 透空块体 ; 结构物 ; 水动力学 ( ) 中图分类号 : TV1312 　　　文献标识码 : A 　　　文章编号 : 2004 ( ) ( ) 近年来 , MAC MarkAndCell 法和 VOF Volume Of Fluid 法等NS 方程求解方法在水动力学研究中已取得明 显进展[ 1 , 2 ] 。但这类数值模型在以下两种情况下产生困难 : ①当自由表面流体出现剧烈变形 ; ②边界几何形 状极不规则 。其原因在于上述两种情况下很难建立计算网格 , 在边界处不得不借助插值等近似处理方法 , 计算 过程烦琐 。由于对流项的存在 , 还不可避免地产生数值耗散 。 另一方面 , 近年来出现的新型防波堤型式由于其结构复杂 , 不便于用欧拉方法进行直接模拟计算 , 在破碎 波作用下尤为如此 。目前实际工程主要依赖实验和经验公式 。在此情况下 , 无网格数值方法将发挥其独特的优 势 。无网格数值方法也称为粒子法 , 其基本思想是将着眼点回归到研究流体运动的基本单元 ———质点上 , 数值 计算的过程即为追踪质点运动学和动力学性质变化的过程 。各质点之间通过核子函数相互作用 。控制方程的离 ( ) 散是通过建立物理量光滑模型、梯度模型和二阶导数模型 Laplace 模型 实现的。由于控制方程中不出现对流 项 , 自然就避免了数值耗散 。 ( ) [3 ] ( ) [4 ] 目前发展起来的无网格数值方法有自由元 Free Element