SPH法在大坝表孔泄流数值模拟中的应用.pdfVIP

第 l9卷 第 6期 水 科 学 进 展 V01．19．No．6 2008年 11月 ADVANCESIN WAIER SCIENCE NOV．，2008 SPH法在大坝表孑L泄流数值模拟中的应用 李大鸣 ，刘江川 ，一，徐亚男 (1．天津大学建筑工程学院港 口与海洋工程教育部重点实验室，天津 300072；2．长春工程学院，吉林 长春 130012) 摘要 ：主要对光滑粒子流体动力学 (SPH)法进行研究，建立了大坝表孔泄流 的光滑粒子模型。将 SPH数学模型应用 于拉西瓦水电站表孔泄流中，提出了采用补水边界的方法来满足库区恒定水位条件，模拟 了表孔泄流的流场变化 及粒子运动过程。通过与物理模型实测的压力值 比较 ，堰表面压力变化基本一致。对模拟结果进行分析 ，表明光滑 粒子流体动力学可用于高速水流的计算模拟研究。 关 键 词：光滑粒子流体动力学；数学模型；大坝表孔泄流；补水边界 中图分类号：TV135．2 文献标识码 ：A 文章编号：1001．6791(2008)06．0841．05 光滑粒子流体动力学 (SPH)法是近年来发展起来的一种 Lagrange方法_1．2J，这种方法最初是用来解决天体物 理中涉及到流体质团无边界情况下在三维空间任意流动 问题的。近年来 SPH法涉及 的领域不断扩展，在水流 运动的计算模拟 中亦得到研究和应用-3．4J。SPH法把连续的物理量用多数粒子的集合来进行插值求解 ，连续体 的运动用有限数量的粒子运动来离散化，由于没有使用网格，因此不会发生界面变形大所引起的计算溢出问 题 。计算空间导数时不需要使用任何网格，而是通过一个称为 “核函数”的积分核进行 “核函数估值”近似， 将流体力学基本方程组转换成数值计算的SPH方程 ，整个流场离散成一系列 “粒子”，所有力学量 由这些粒子 负载，这些粒子可以按计算公式，亦即按流体力学流动的规律任意地运动 。 1 SPH基本方程 SPH法直接将流体看作为具有密度和压力等参量的粒子，通过在邻近粒子上的插值来计算密度、压力、温 度等变化，场中的任一宏观变量 A(r)都能借助于一组无序点上的值表示成积分插值计算得到。 A(r)=IA(r)W(r—r，h)dr (1) JD 式中 D为整个求解的区域 ；A(r)为预先设置点的宏观变量值 ；h为光滑长度 ；W(r—r，，h)称为核函数，它 有两个 自变量 ：粒子间距离 Ir—rI和光滑长度 h。 引入算子 ，可得到SPH基本方程如下 ： 质量守恒方程 dtoi ： ∑ (一 )v (ri一0，) 动量守恒方程 dui = 一 砉 考 ) 能量守恒方程 del = 筹VI(一)V(r一。，∽ 收稿El期 ：2007．10．24 作者简介：李大鸣(1957一)，男，天津人，教授，博士，博士生导师。主要从事水力学及河流动力学、海岸工程的研究。 E-mail：Lidaming@lju．edu．O13 842 水 科 学 进 展 第 19卷 质点运动位置 dr—A一+e砉tj／--tti附 ) 式中 lD为密度；m为质量；P为压力；u为速度；e为内能；W为核函数；e为0e1的常数； ：古(Pl+ )。 r SPH中核函数要满足一些基本条件[．6】，主要有以下几个方面：①满足归一化条件。I (r—r，h)dr=1；