浅水波浪传播的数值模拟.pdfVIP

浅水波浪传播的数值模拟+ 柳淑笋俞聿修赖国璋 (大连理工大学海岸样垃海工程国家重点实验室，116024) 李毓湘 (香港理工大学土木及结构工程系) 摘 要 本文基于有限元方法，建立了求解Boussmesq方程的数值模型，模型中，采用节点周围单元 内的一阶导数的加权平均来确定单元节点上变量的一阶空间导教德．使模型可以采用线性单元，而对闯 的，可以模拟波浪的非线性传播过程。 ，关键词：有限元方法Boussinesq方程数值模型 前 言 近岸海浪条件是影响海岸建筑物设计的重要因素之一，由于海浪在浅水区传播时受海 岸建筑物及海底地形变化的影响，将发生波浪的绕射及折射，因此如何精确描述及模拟这 一过程成为历来重要的研究课题之一。 近岸浅水区波浪的传播过程是一复杂的非线性变化过程，Hmm等“1对有关的研究进行 了综台回顾，BoIl：ssinesq形式的方程中包含了弱的非线性及色散性，是描述海岸浅水速域 波浪传播较理想的数学模型，由于该方程沿水深方向进行积分，将三雏问题化为二维问题 求解，使大面积的波浪传播的数值模拟成为可能。 由于在数值计算中，有限差分法编程较简单，容易实现，而且计算速度较快，因此目 数值模型的求解精度问题，但众所周知，有限差分法具有～明显的弱点，不能对于海岸波 浪计算中经常遇到的复杂边界问题进行灵活处理。相反，尽管有限元法需要较复杂的编 程，但它可灵活处理各种较复杂的边界，因此该方法特别适用于海岸工程问题的计算，但 相对于差分法来说较困难，尤其是一般需模拟的海岸波浪传播区域一般较大，时域求解更 显困难，因此，目前应用有限元方法求解Boussin∞q方程较少，本文采用有限元方法，基 建立一求解波浪传播的数值模型，长距离孤立波传播的计算结果表明，本文提出的计算模 型是稳定的。 2数值计算模型 2．1有限元方程 ’奉文为田枣茸端事l学基盘{编号·$9479006)鬻助项目杆骨研觅成墨 仉+甲．b+_-】=。 VfV-oV胡 ”加w嘴V州-+芦善审p·帆)】+碰2 (2) d 一rt+口乒vp．0，)J～F誓v【V2 0 0 其中“=“一为二维沿水深平均的水质点速度，带为波面，^勘阮∥为水深，亭为重力 加速度，P为一常数，由方程的色散性，其值取为I／5。 将上述方程中的求导项展开，并将时间导数项归并～起，如果地形变化较缓，水深变 量h的二阶空间导数属于小量可以忽略，并设求解区域为口取，为方程的权函数，ⅢⅡ方程 的Galerkin加权余量式为 g卜却+昙晒)]．懒=。 (3) 以+，西西+g尘一型掣一型丝一pg．*堡+Vm 型一塑堡一堂封黜护。 射”矗旁增万一了矿一了面一万万一了i万一了万甜脚％ 矿 鼠^ 一 ^ 一 一 一 waxay=0 ]“j (6) 屈～，n一， 矿 虬，“， 反 冉 一 履一：岛一： ^ 葫～印葫一办 一 反一：反一： ^ A一毋西一乐 一 Waxdy=0 出心 如一分如～铆 办一埘办～矿 葫一【毽葫一旁 a一国a一砂 毋一盘a一印 葫一屠葫一砂 1叫J (7) 其中届=l帮，肛而切。 中的二阶空间导数项进行降阶处理。按有限元的求解方法，将求解区域锵4分成有限个单 元并将一典型单元记为口，可以得到其相应的单元有限元方程为 hfh奇。+Q嘞H。杜J+R嘶HlvJ+Q畸ufH]．+R畸v。Hi=0 (8) hfhP。+Q晦H。HJ+R畸v。MJ+sb