时域内数值水槽的阻尼消波与Ｓｔｏｋｅｓ波数模计算.pdfVIP

时域内数值水槽的阻尼消波和 Stokes波数模计算 孙大鹏李玉成 (大连理工大学海岸和近海工程国家重点实验室， 116024) 摘要本文采用时域内对波面运动位置追踪的边界元方法，建土丁一种非线性波浪变形计算的三维数值 模式，并借助Sponge 波的数值造波水槽，进而对水槽内的Slokes波进行了波浪变形计算． 关键词数值造浊水槽·阻尼消波，{!婴es躜，边界元 1前言 无论是线性波浪变形计算还是非线性遮浪变形计算．数值模型中都将涉及到界面上边 界条件怎样设定的问题。通常情况F，不考虑底面上的渗流影响．我们可以给出底部速度 势函数法向导数为零的边界条件：但是如果波场的l司体岸壁有波能吸收，或波场中存在着 透水建筑物．那么如何给出波浪的部分反射边界条件?另外．发生在计算域内的这种反射 波以及建筑物周边的散乱波将向外传播和扩散．这时如何确定开敞边羿(即入射边界)处势 函数或许函数法向导数边界条什?这些问题的解次不但关系剑数值模，弘的完备与否亦将 J．p 影响数值计镡精度。应该说，在线性理论下，关丁二边界处理的数值方法已日趋成熟ll J， 在开敝边界处通过内、外域波势函数的榴台匹配给出包括敝乱波在内的入射波边界条r{：， 从而可以较灵活地选取开敞边界的位置，而无需将计算域选取得很犬：在嗣体岸艇上，采 _ffj波浪在任意波高反射系数下的边界条栩。，使得数学模型既使对于复杂j犁壁形状、1I恒定 水深区域的波场变形问题，亦能够正确地给出其波场波势函数的数值解。 但是．在处理时域内的非线性波浪变形问题时，在固体岸壁上．适用于线性波浪理论 的部分反射条件己失去其成立的前提，如果允许反射波亦夹杂在域内波场中，那么将使得 我们几乎无法正确给出开敞边界处的波浪入射条件。通常的做法是仿效差分法，把固体岸 壁处理成Sommerfeld型的反射边界，即波势函数满足： 掣+c霉：o (I) 口 日 由于上式中，系数c是单一的波动位相速度，对于具有非线性变形特性的波浪，边界 的吸收、透过效率将会大大降低刚。在探讨、比较几种边界处理方式的基础上．本文首先 在反射边界前的一段距离内，布设具有消波功能的阻尼段(Spongelayer)P)衰减波能的大部 分，同时仍然利用(1)式的放射条件，使得波能未被衰减的部分能够透过边界传到域外． 即采用一种“阻尼+透过”的边界处理方式。计算结果表明：这种边界处理方式可以保证 在有限区段内非线性波浪的连续计算，适合于非线性数值波动水槽中的波浪变形数值模 拟。 ·233· 2水槽内的阻尼消波 采BJ“阻尼十透过”的消波处理方式，即在反射边界前设置一段长度为口的Sponge 动·在Spongelayer区域Ek内为牯性流陆运动。如图1。 —缸互 ．， 辟———工—————刊 嗣1 I哿er阻尼消波层的布设 Sponge 以往【j’对Spongelayer阻尼消波层的应用经验表明，对于波长大于Spongelayer层长 度口的波浪，Spongelayel的消波效率并非很高。而且，从计算机容量及占用计算机时的 方面考虑，单纯地加大Spongelayer层K度以谋求消波效果的做法无疑是很不经济的。为 此·在Spongelayer阻尼消波层的后边界面曼上，设置透过边界．即波势函数满足 Sommerfeld辐射条件。假设Sponge layer内的流体运动阻尼系数∥(x)在Spongelayer内取 撒线性分布，井认为Sponge lay