A One-Dimensional Integral Equation for a Porous Horizontal Disc under Water Waves

技术摘要：关于水波作用下水平多孔圆盘的一维积分方程

问题陈述
本研究探讨了淹没在深水中的薄型多孔圆形平板所产生的波散射与辐射现象。主要目标是建立并求解流体动力学问题，以确定作用在该结构上的附加质量与阻尼系数。分析重点在于复杂的多孔效应参数 $G$ 的影响，该参数表征了平板的阻力和惯性特性。该问题在数学上被建模为定义在单位圆盘上的二类超奇异 Fredholm 积分方程。

研究方法
研究采用了势流理论框架，假设流体运动为小振幅、无旋、不可压缩且无粘性的。边界值问题被简化为关于圆盘两侧速度势跳跃的超奇异积分方程。

关键方法步骤包括：

1. 降维：遵循 Farina & Martin (1997) 建立的数学框架，在轴对称运动的假设下，将三维问题简化为一维积分方程。这是通过将物理量展开为傅里叶级数并引入一个未知的辅助函数 $\psi$ 来实现的。
2. 控制方程：降维产生了一个具有奇异核的二类 Fredholm 积分方程（方程 21）。该方程包含了波数 $K$、多孔参数 $G$ 以及淹没深度 $d$。
3. 数值解法：由于方程的超奇异性质，无法直接应用带有高斯-勒让德求积法的标准 Nyström 方法。因此，作者利用边界元法（BEM）结合 NAG 库中的 D01GCF 程序进行求解。该程序采用 Korobov-Conroy 数论方法（NTM）来严格评估潜在积分。通过分解核函数并将积分区域划分为紧邻奇异点的子区间，实现了对奇异性的处理。
4. 验证：通过网格独立性测试以及与既有文献（特别是 De Freitas 等人 (2021) 以及 Farina & Martin (1997) 的工作）的结果对比，对数值实现进行了验证。

主要贡献与结果
本文提出了一种经过验证的数值公式，用于计算淹没多孔圆盘的流体动力学作用力。研究结果针对多孔参数 $G$（实数、虚数及复数）和淹没深度 $d$ 的五种特定场景进行了分析：

• 多孔性的影响：研究展示了基于 $G$ 性质的不同行为：
  • 实数 $G$（阻性）：代表高阻力多孔结构。随着实部增加，圆盘的行为趋向于固体平板，阻碍流体渗透。这导致与低阻力配置相比，附加质量和阻尼效应有所降低。
  • 虚数 $G$（惯性）：代表由孔隙内惯性效应主导的结构。在这种机制下，流体倾向于跟随圆盘运动，表现为有效的附加质量。这会导致附加质量和阻尼系数出现显著峰值，特别是在靠近自由表面的配置中，会诱发圆盘上方流体层的共振。
• 数值一致性：所提出的公式与 De Freitas 等人 (2021) 的参考数据表现出极佳的一致性。通过平均绝对误差（MAE）和均方根误差（RMSE）进行的定量验证，确认了所得附加质量与阻尼系数在各种波数（$Ka$）和深度下的连贯性。
• 与经典理论的联系：研究发现，当波数 $K=0$ 时，控制方程简化为类似于经典的 Love-Lieb 方程的形式。这建立了淹没多孔圆盘的流体动力学问题与同轴圆盘静电问题之间的理论联系，表明所提出的公式是 Love 方程在流固耦合背景下的推广。

意义与主张
作者声称，本研究成功地将最初为刚性平板开发的 Farina & Martin (1997) 的分析策略，应用到了多孔物体的情境中。其主要意义在于推导出了一个精简的一维积分公式，该公式在容纳多孔介质复杂性的同时，保留了原始模型的核心数学结构。

论文谦逊地指出，虽然数值一致性验证了该方法，但处理此类方程中的奇异性仍然是一个挑战，需要开发超越标准求积方法的替代数值策略。本研究局限于较小的淹没深度（$d=0.1$ 和 $d=0.2$），侧重于由靠近自由表面所强化的共振现象。作者得出结论，复数多孔参数 $G$ 是控制入射流体动力作用极端值的关键调节机制，其中实部起耗散作用，而虚部则会放大共振。