非线性海岸低频波浪地研究.pdfVIP

非线性海岸低频波浪的研究 董国海 (大连理工大学．大连116024) 王晓琳 (大连大学，大连116011) j囊要t本文研究了海岸低频波浪的计算方法．波浪控制方程由欧拉方程在短波周期上平均而得到．为了 项采用了时间算子分裂法进行了处理，并对由此产生的误差进行了修正．波浪能量方程由Lax-Wendroff 差分格式求解．计算了双色波垂直入射海岸的情况，并与实验结果进行了对比． 关键词：碎波拍．破碎波浪．低频波浪 1引言 wBves)而产生的低频波浪 碎波拍(surfbeat)是海岸处伴随入射波群(incident (10wfrequencywaves)(周期为30秒以上达几分钟)，大量研究表明其对海岸泥沙输运， 港口波浪和系泊船运动都有重要作用。Roelvink!”求解以流量表达的一维低频波浪方程， 研究了碎波拍对海岸变形的作用。通过坐标变换将计算域变换为固定域，处理了海岸处动 Flux)方法…用来求解一维长 边界。本文将用于求解短波方程的WAF(WeightedAverage 波方程。对方程中与辐射应力和底摩擦力有关的非齐次项，采用了时间算子分裂法进行处 理，并对由此产生的误差进行了修正。计算对双色波进行，并与实验结果进行了对比。 2波浪方程 假定沿海岸方向地形变化是均匀的，这样方程沿水深平均后所研究的问题是一维的。 对欧拉方程沿水深和短波周期平均，得 妄E础+j8一fh(pu2+p)嘞dh托=o (1) 式中，p，h和“分别为流体密度、静水水深和流体速度，f为自由表面升高．p，p6和 “分别为压力、水底面上压力和海底摩擦力。“～”表示在短波周期上时平均。将速度“ 分解为低频波浪速度u和短波速度“’，即“=u+”’。由流体流量为Qf=Ez础定义 质量传输速度 矿=Q／d=U+玑／d 包 (2) 一59一 其中， c) (3) 玑=臣“r出=E／(p 为短波平均流量。而E和c分别为波浪能量和短波相速度。将以上表达式代入方程(1) 及连续方程，得由短波相互作用产生的长波的控制方程 要d+兰归)=0 (4) 砉怕，+丢(矿2d+；鲥2]=∥罢一耋[鲁]+妄(譬)2一吾cs， 式中，d：7+h，手为平均自由表面升高，g为重力加速度，s。为辐射应力 ＆=酾一；倒k【z詈一捧 ㈤ 波能演化方程 -国D-+羞陋+咯瑚—一岛 (，) 式中，Db表示由波浪破碎引起的波能耗散，通常由“水跃”能量耗散确定。本文采用Roelvink 的方法f11计算 斗唧㈣“胁 cs， 式中，取口=1,0和H=10。正为特征频率，对双色波取平均频率。y为波浪破碎指标，H 为波高日；√F，昭‘，D^在短波破碎前为零，波浪破碎点由关系式H=一确定。 3方程求解方法 由于低频波浪方程(4)和(5)是非齐次的，本文采用时间算子分裂法处理非齐次项， 即首先求解以下齐次方程得到中间值d’和V’。 警+妄咐)=。 (9)