孤立波传播VOF数值模拟.pdf

孤立波传播VOF 数值模拟 倪培桐 江 洧 （广东省水利水电科学研究院，广州，510630） 摘 要：本文利用VOF 方法模拟自由水面，利用通度概念处理固-液界面，模拟了孤立波对直立墙的冲击过程，实 例研究表明数值模拟结果较好。 1 前言 对于直墙式建筑物/斜坡堤而言，破波作用是该类建筑物破坏的重要原因,我们需要了解冲击 荷载的详实数据及破波作用于建筑物时压力的时空分布。到目前为止，该类问题仍然缺乏理论解。 物理模型实验一直是解决该类问题的主要途径，近年来随着数值方法及计算机速度的不断提高， 数学模型的研究到学术界的重视。 早期的波浪数值模型多是基于势流理论来研究破波对建筑物的冲击荷载。由于自由表面附近 水流并非势流，这些方法在不能描述整个破波过程和冲击过程。随着 VOF （Volume of fluid, Nichols, etc,1980）方法的发展及改进，可用于描述自由表面破波及近岸斜坡破波作用和冲击作 用的波浪数值模型研究取得较大进展。王永学用VOF 方法分别对波浪在直墙前的破碎和数值波浪 槽的造波问题进行了研究，邹志利等研究了数值波浪槽生成非线性波（孤立波和椭圆余弦波）的 过程，万德成等应用VOF 技术研究了孤立波翻越直立方柱流动的全过程。基于VOF 技术，N.T.WU 等研究了二维N-S 方程研究破波对直墙的冲击过程，结果表明当破波过程中不存在掺气时，数模 结果在定性与定量上都与实验结果吻合较好，当存在掺气时，利用实验公式作适当修正仍然可以 得到较好的结果。M.R.A.Van Gent 等应用二维N-S 方程结合渗流模型研究了涌浪对渗透结构建筑 物的冲击作用。目前数值模拟的主要对象是规则波，而对于非规则波的研究有待于进一步的提高 （M.R.A.Van Gent，1994）。本文主要基于VOF 方法研究孤立波在堤前的传播及变形过程。 2 数值模型 2.1 控制方程 ∂ ∂ ∂ ∂ ∂2 ∂2  u + u + u − 1 p + +  u + u  （2-1） ∂t u ∂x v ∂y ρ ∂x g x υ∂x 2 ∂y 2   2 2  ∂ ∂ ∂ ∂ ∂ ∂ v + v + v − 1 p + +  v + v  （2-1） u v gy υ 2 2  ∂ ∂ ∂ ∂ ∂ ∂ t x y ρ y  x y  ∂u ∂v + 0 （2-3） ∂x ∂y 其中u,v ρ 为速度分量，g x , g y 是重力加速度分量，υ运动粘性系数。 为流体密度。 定义 F 函数代表流体占据网格单元的体积比值。1 代表网格充满流体，0 代表空网格。若 F 值在0 与１之间，则表示网格中含有一个自由表面。Ｆ随时间的变化由下控制方程描述： 50 ∂ ∂ ∂ F uF vF + + 0 （2-4） ∂