非定常二维水流有限元解及其应用.pdfVIP

非定常二维水流有限元解及应用 来军政 耿兆铨 0016 浙江省河口海岸研究所，杭州，31 摘要：本文应用一种基于阶梯函数特性和集中质量概念推导得到的显式迎风有限 元格式，求解水深平均的二维浅水流方程r数值模拟钱塘江尖山河湾实际流场·计 算结果与实际相符。 一、前言 计算水力学中，非定常二维水流的数值模拟发展了很多差分算法，如特征线法、 显式差分法、ADI差分法、分裂算子法等。由于水流复杂的边界几何形态．上述差分 法即使采用贴体坐标．实际应用中仍然很不方便。为此，基于三角形单元的各类算法 有着广泛的应用，如有限节点法、分步杂交有限元法、显式迎风有限元法等。 ● 本文运用阶梯函数特性(不连续有艰元)和集中质量概念推导了一种的迎风有限 元显格式。本文格式采用两种基函数，一为线性基函数．另～为不连续基函数(阶梯 函数)，因不涉及矩阵求逆问题．程序简单易编，内存需求小．计算迅速，同时对流 项粟用迎风处理技巧，使得格式在满足柯朗条件下．具有良好的收敛性。实例计算表 明，本文格式边界适应性强，计算稳定性好．计算省时且内存小的特色。 二、控制方程 流场计算的控制万程选用沿水深积分的二维泼水潮波方程： az．+掣+_a(Hv)：o (1) 吼 “ 口y 乱 孔 扎0z．“雨 a aJ a a1 瓦+“丙+”西+8夏～如+{；u—丽2磊(8x夏)+丙(8v≯(2) h h h 藐 j{雨。 弧 己 h i+“瓦+v_cry+8丙+m+gv—丽2夏(8x磊)+丙(8，丽)(3) 式中： z为潮位·即水面到某一基准面的距离：u、v分别为x、j，方向上的垂线平 均流速分量： H为水深，即水面到水底的距离；g为重力加速度；f为柯氏力参数 (f：2(o n 为糙率系数：￡、，￡，分别为x、y方向的涡动扩散系数：x，Y为直角坐标：t为时间。 ● 上述方程组是由雷诺方程组在下列假设下沿水深积分得到的：1)水流具有不可 ／ 压缩性：2)压强服从静压分布：3)忽略垂直加速度：4)粘性切应力在水平方向 上的变化率远较垂直方向的变化率小． 可忽略不汁。在推导过程中没有考虑水表面的 切应力。 初始条件： ● z(x枞，z。(Ⅵ) u(x硼：0=u。(x，y) v(Ky)l。。=v。(x，Y) 边界条件： z(x，Y，t)=z’(x，Y，t) 水边界 “}，，表示已知值 一V-=0 陆边界 法线方向流速为霉。 三、有限元求解格式 有了上述定解条件．就可用～定的离散格式求出