关于缓速阻燃波的数值模拟.docxVIP

关于缓速阻燃波的数值模拟 规则波的定义 由于海滩非常缓慢，波浪变形非常强烈，而且底部的能耗非常显著。波浪变形和破碎有其特殊性。提出了以下问题：采用何种波浪变形计算方法？采用什么破碎指标？以及，由于大量土壤能耗，波传播在一定的条件下是不可能发生的。对于i1.50的背景，文献从不同角度进行了讨论，但对缓冲区河岸的研究较少。例如，nelson和bara获得了一些矛盾的结果，但作者也分析了这些结果。这些研究是基于规则波。对天然情况的波浪是属于随机性质的不规则波,却研究甚少。本文对此作一些分析与讨论。 1 规则波khb/nb 判别波浪是否破碎可以采用几何学指标(GSP)、运动学指标(KSP)及动力学指标(DSP)。(DSP)以水质点垂直加速度值为判据,该值在濒临破碎时变化剧烈,判断值难以选取;(KSP)以表面水质点速度u与波速C比值≥1.0作为准则,其物理含义明确,但u依赖于所采用的波浪理论,在不规则波中计算比较困难;(GSP)主要以当地破碎波高Hb与水深db之比Hb/db或波陡参数α=(Hb/Lb)/tanh(2πdb/Lb)作为判据,Lb为破碎波波长,在规则波条件下它们可由运动学指标导出,工程应用简便。文献在岸坡i≥1/50条件下,将这两个指标用于不规则波条件。在文献中,对于规则波Hb/db值可按照Goda方法确定: (Ηb/db)=A?{1-exp[-1.5πdbL0(1+1.5i4/3)]}/(dbL0),(1)(Hb/db)=A?{1?exp[?1.5πdbL0(1+1.5i4/3)]}/(dbL0),(1) 式中L0为深水波长,A为常数(0.17)。文献指出,对于不规则波,可采用式(1),但应做修正:a.A值应降为0.15(对发生破碎的大波)或为0.12(对波列中的有效波高Hs);b.不规则波的波长比规则波理论波长短,式中L0应为线性理论波长(L0=gT2/2π,g为重力加速度,T为波周期)的0.74。本文中采用Hb/db及α作为平缓岸坡上的不规则波的破碎指标,并通过实验与分析来确定底坡i1/100时不规则波的破碎指标。 1.1 波谱和波高值测定 实验在海岸和近海工程国家重点实验室波浪水槽内进行,水槽长69 m,宽2.0 m,高1.8 m。水槽一端配置液压驱动不规则波造波机,另一端设有1/6的消浪斜坡。模型底坡1/200,后方为一平台,它和浪高仪布置如图1。沿水槽设5个测站,1#站用于测初始波要素,距造波板约7 m,2#站在缓坡的前端,距1#站7 m,4#站设在缓坡末端,其后方为平台,3#站和2#站及4#站的距离分别为12 m和11 m;5#站设在4#站后3 m处。 实验时,设定波谱谱峰周期后,调节波谱能量值,使波浪在3#站与4#站区间发生破碎。另外在3#站和4#站前10 cm处分别加设一个浪高仪,用于测定当地的不规则波波长。5#测点用以测定破碎后的波况。同时目测3#站及4#站逐个波浪是否破碎和破碎情况,以确定各个破碎波的波高值等。 实验时取两种水深:4#站水深分别为15 cm及25 cm,对应的3#站水深分别为20.5 cm和30.5 cm,1#站水深为70 cm和80 cm。谱峰周期在浅水条件下为1.1 s,1.3 s和1.5 s,在深水条件下为1.2 s,1.4 s和1.6 s。 1.2 anh2d/l的对应值点 将实测破碎波的波长L实与同周期按线性理论计算的浅水波长L实=(gT2/2π)tanh 2πd/L的对应值点绘在图2上。以最小二乘法原理拟合二者关系,得L实=0.74L线,其相关系数r=0.962,标准差σ为0.07。此结果与Ochi的结果(L实=0.74L线)完全吻合,和文献的结果(L实=0.75L线)也十分相近。 1.3 有效值的验证 在本次实验的不规则破碎中,每个破碎波的(Hb/db)和db/L0的关系见图3。L0为修正后的波长。根据式(1),用最小二乘法拟合得A=0.15,标准差为0.031。此值与文献在1/50底坡上的结果完全相同。 图4为实测最大波高的Hmax/db和有效波高的Hs/db与db/L0相关图。同样,采用最小二乘法拟合相关数据,分别得相关的式(1)中A值各为0.162及0.11。其中有效值所对应的A=0.11略小于文献和文献分别在1/10～1/40和1/50岸坡上的结果(A=0.12)。如取本次实验点的外包线,则仍然可得A=0.12。 1.4 破碎波 根据各破碎波的实测Hb及L0可分别计算3#站和4#站各破碎波的α=(H/L)b/tanh(2πdb/Lb)值,其变化范围为0.114～0.131,平均值为0.120。对于3#站和4#站,该值无大差异。 1.5 谷深度t/t 如图5,在波面过程线上表示波形不对称性参数有三个:一是波峰前部历时T′与背部历时T″之比T′/T″;二是(T′+T″