显式迎流FEM在河口近海研究中应用.pdfVIP

显式迎流FEM在河口近海研究中的应用 倪勇强耿兆铨朱军政虞开森潘存鸿林炳尧楼越平 (浙江省河口海岸研究所，杭州310016) 擅叠· 河口近海特殊的地理形态和复杂的动力因幕，使FEM在其数值研究中应用越来越广。 本文介绍的显式迎流FEM，}|!f：具有显式FDM省内存、省机时的优点，又具有边界适应矬强、节点布设 灵活且便于局部加密网格的特色。模式与计弊前厝处理技术配套而成的软件包已广泛成用于潮汐河口、 近海、港湾、湖泊及诸多岛汉的复杂水域的计算，除流场外，还成功应用于河口规划整治工程方案研 究和污染物浓度场、温度场、盐度、堀沙及事故溢油等计算模拟。文章给出了应用实例，并讨论了有 关问胚。 关t葡： 迎流显式有限元}i|『Lj敬值模拟前处理后处理 1前言 潮汐河口的规划整治和歼发利用在国民经济中举足轻重。随着改革、开放的深入进 展，新建、改建或扩建的I：矿企业如雨后春笋，既要发展经济，又要保护好环境己成共 识。河口近海的环境水力学问蹶也越来越引起人们的重视。研究方法除了现场观测分析 外，物理模型试验、数值模拟计算也是行之有效的手段。然而，现场观测往往受制于经 费和天气，物理模型也常囿丁时间和场地设各。因此，采用数学模型越来越广泛。随着 计算机功能的飞速发展和现代计算技术的不断提高，数值模拟水平越来越高，研究成果 也层出不穷。就计算网格而言，二维数模通常有四边形和三角形之分。由于河口常具喇 叭形平面形态和起伏的水下地形，且兼受山、潮水作用，规则网格较难适应，而无序网 橹(如三角形)则比正交曲线网格有更好的灵活性。加之潮动力因素的较快变化．取较短 时间步长的显格式也是恰当的。这正是显式迎风三角形FEM能在河口研究中获广’泛应用 的主要原因。 本模式先后应用于钱塘江、甬江、瓯江、IIs云江、鸭绿江等河口，西湖、太湖等湖 泊，杭州湾、三门湾、深圳湾、三亚湾等海湾及舟山群岛等复杂水域的流场和物质浓度 式以及二维与一维耦合的计算模斌等。对于层化明显的河口，三维分层流模式正在耦台 凋试中。 2数学模型简介 2．1水流模式‘n 以2DH模式为例，流场计算的基本方程选用沿垂线积分的二维浅水潮波方程(含=阶 涡动粘滞项以反映有效应力的作用)： 塑+型+她：o (j) at ax 0v ～572—— 等+u詈+VT抽y+g詈一fv～鲁+gu≤≥}i”瓦”+丙一～苜+“{■r =砉ce。等，+专ce，雾， U 十 十 √u2+v1 堕孔 虬『 5”1可 ^一 h 十 E ，言一l蓦。～抓 扑一曲I若一h 学争 旷争 式中：z为潮位(米)；u，v分别为x，Y方向上的垂线平均流速分量(米／秒)；H为水深f米)； 系数；宕x，8，分别为x，Y方向的涡动扩散系数；w。，W，分别为x，Y方向的风应力分量：x,y直角坐 标；t为时间。 方程(1)为水流连续方程，方程(2)、(3)为x，Y方向的动量守恒方程。 上述方程组的初始条件： Z《x’y)l。。2 z。x，y)；u(x，列。=u。(x，y)：v(x，y】。=vo(x，y) 边界条件：水边界z(x，y，t-)=z+(x，y，t) 一”表示已知值： 陆边界寸：o 法向流速为零。 本文选用显式迎流有限元格式。’，即对式(1)=～(3)应用加权余量法伽辽金逼近， 其中二阶项运用分部积分。对流项采用迎流格式，时闻导数采用前差公式。采用三节点 三角形单元面积坐标，利用阶梯函数特性和集中质量概念(图1)．极易求得其彤参数， 代入并移项化简，可得局部有限元方程适宜编制计算机程序的计算公式： z∥=z；一△t盯丽“+丽“派^) (4) 叮‘=u-At叫(驯．Ⅵ、剖yu：卜g鲁_fv?