三棱柱形网格下自由表面流的三维数值模拟5x∫5y-水科学进展.pdfVIP

第 20 卷 第 4 期 水 科 学 进 展 Vol20 ,No4 2009 年 7 月 　　 ADVANCES IN WATER SCIENCE 　 Jul . , 2009 　 三棱柱形网格下自由表面流的三维数值模拟 1 ,2 1 ,2 1 ,2 1 ,2 王 　昆 , 金 　生 , 马志强 , 高述峰 ( 1 港口海岸及近海工程国家重点实验室 , 辽宁 大连 　116024 ; 2 大连理工大学土木水利学院 , 辽宁 大连 　116024) 摘要 : 为了模拟自由表面流动问题 , 建立了基于三棱柱形网格的三维 NavierStokes 方程的半隐式动压离散模型 。出 ε 流边界采用在动量方程右端附加线性衰减项的海绵层处理方式来消减波能 。引入 k 双方程紊流模型求解涡黏性系 数 , 使方程达到封闭。用规则周期波通过淹没障碍物的传播变形和复合明渠流动算例 , 验证了当遇到短波高频问题 时 , 所述模型模拟结果与测量值吻合良好 , 体现了模型的精确性和有效性 。 ε 关 　键 　词 : 三棱柱形网格 ; 自由表面流 ; NavierStokes 方程 ; k 紊流模型 ( ) 中图分类号 : TV1314 　　　文献标识码 : A 　　　文章编号 : 2009 自由表面流动 问题[ 12 ] 的数值模拟一直受到人们 的广泛关注 , 曾用过的方法主要有 marker and cell[3 ] 、 VOF[4 ] 、arbitrary LagrangianEulerian[5 ] 等 , 这几种方法各 自存在着不足 , 花费的计算代价难 以令人满意 。鉴于 此 , 建立了自由表面流动的静压模型[6 ] , 但三维静压控制方程中无垂向动量方程 , 无速度修正 , 仅在计算长波 问题时比较适用 , 当遇到如波浪运动的短波问题时 , 结果并不准确 , 一直应用比较广泛源于它可以节省计算成 本 。本文引入了自由面流动的动压模型 , 建立了基于三棱柱形网格的三维 NavierStokes 方程的半隐式离散格 式 , 可以很好地解决复杂计算域在边界处的精确模拟问题 。通过两个算例验证了在处理短波高频问题时 , 本模 型所得的结果与测量值吻合良好 , 体现了模型的精确性和有效性 。 1 　模型方程 不考虑斜压项及大气压力项影响的不可压流体的动量方程为[78 ] 2 2 d u η q h u u v u γ +