平面二维水流数学模型的建立与应用.docxVIP

平面二维水流数学模型的建立与应用 在计算二维水流运动的数值时，很难解决区域边界的复杂性。如何处理不规则边界，直接影响模型计算的可靠性和精度。传统的笛卡尔测量方法具有简单的网格校正原理，适合计算网格位移，但由于限制因素易于形成锯齿限制，边界值调整的精度不高。为了提高该方法对非规则区域的适应性，本工作介绍了一种新的不规则边界处理方法。该方法不仅保留了传统笛卡尔坐标方法的简单、正交、容易形成计算网格和计算效率的特点，而且改进了传统锯齿限制精度差的缺点。它靠近实际情况，使边境地区的洪水分布变得合理。在复杂边界分洪区的洪水计算中，具有很高的精度。 斜对角笛卡尔方法基本原则就是按就近逼近原则,把距实际水边线较近的网格节点连成线以拟合实际水边线.生成网格时,根据边界是否固定,可采取不同的判定方法.当边界为固定边界时,可把固定边界离散为若干离散节点,确定各离散节点坐标,按距离最小原则选定固定边界上各离散节点的逼近网格节点,然后依序连接各网格节点即形成拟合边界;对于动边界,由于边界线变化不定,则可借助网格节点间水深的干湿变化和水深值,参照相应网格节点间的地形高差,判定距离实际水边线较近的网格节点,合并纵横向的拟合网格节点,形成拟合瞬时水边线.根据斜对角笛卡尔方法基本原则,即以距交点最近为原则,选取对应的网格节点作为实际水边线的逼近点. 1 模型介绍 1.1 单宽流量n 平面二维水流数学模型的基本方程组由水流连续性方程和水流运动方程组成. 水流连续性方程: ?Z?t+?M?x+?N?y=0(1)?Ζ?t+?Μ?x+?Ν?y=0(1) 水流运动方程: ?M?t+?uM?x+?vM?y=?gh?Z?x?gn2uu2+v2√h1/3(2)?N?t+?uN?x+?vN?y=?gh?Z?x?gn2uu2+v2√h1/3(3)?Μ?t+?uΜ?x+?vΜ?y=-gh?Ζ?x-gn2uu2+v2h1/3(2)?Ν?t+?uΝ?x+?vΝ?y=-gh?Ζ?x-gn2uu2+v2h1/3(3) 式中:Z为水位,m;h为水深,m;u,v为垂线平均流速在x,y方向的分量,m/s;M,N为单宽流量在x,y方向的分量,M=hu,N=hv,m2/s;n为曼宁糙率系数;g为重力加速度,m/s2. 本模型采用有限体积法计算,离散的详细过程见文献.有限体积法兼有有限差分法物理概念清晰和有限元法适应不规则网格、复杂边界情况及计算精度高等方面的特点.有限体积法又称有限控制容积积分法,即把计算区域按划分的网格离散为若干点,以这些点为中心,把整个计算区域划分为若干互相连接但不重叠的控制体.计算物理量采用交错网格,其优越性在于能使水量保持守恒. 1.2 角形网格单元 采用斜对角笛卡尔法对实际水边界处理后,计算边界网格单元较传统笛卡尔锯齿法多4种类型,共8种基本类型(如图1所示). 由图1可知,边界网格单元不外乎四边形和三角形两种类型,计算时可通过网格标示数组判别.基本方程在网格上离散后,边界三角形网格单元对应的水流参数(单元形心)与原四边形网格单元不同,其分布在三角形形心处,而边界处流速点存储位置不变.为了求解计算方法的统一,需对三角形边界网格单元进行处理,若把三角形边界网格单元斜边中点视为一角点,则仍可视为四边形网格,与内部计算网格无异,可用统一方程求解. 由于三角形网格单元水深存储位置变化,而流速点存储位置不变,故对方程组离散参数影响主要体现在hi,j+1/2,hi+1,j+1/2,hi+1/2,j,hi+1/2,j+1.以hi,j+1/2为例,其值受不同边界类型组合影响,可通过加权平均的方法确定各类边界的hi,j+1/2值,基本类型有8种,如图2所示. (a),(b),(d),(e)边界类型有: hi,j+1/2=12(hi?1/2,j+1/2+hi+1/2,j+1/2)hi,j+1/2=12(hi-1/2,j+1/2+hi+1/2,j+1/2) (c),(f)边界类型有: hi,j+1/2=(35hi?1/2,j+1/2+25hi+1/2,j+1/2)hi,j+1/2=(35hi-1/2,j+1/2+25hi+1/2,j+1/2) (c),(f)边界类型有: hi,j+1/2=(25hi?1/2,j+1/2+35hi+1/2,j+1/2)hi,j+1/2=(25hi-1/2,j+1/2+35hi+1/2,j+1/2) 三角形网格边界处其他点处水深如hi+1,j+1/2,hi+1/2,j,hi+1/2,j+1处理与之相同. 1.3 大陆地区边界 为检验采用斜对角笛卡尔方法改进后的数学模型较传统笛卡尔方法优越性,本文选取一概化的具有复杂边界的分洪区为计算对象,如图3所示. 计算时,水流边界条件为进口10 m宽扒口,分洪流量恒为50 m3/s,外边界为四周封闭围堤.陆地边界的处理方法为