什么是波浪衍射与港湾波浪计算
🙋
防波堤背后不是应该没有波浪吗？为什么有时候还会有，这个“衍射系数K_D”是什么？
🎓
简单来说，K_D就是衡量波浪能“绕”到防波堤背后多少的一个比例尺。K_D=1表示背后的波高和外面的入射波一样高，完全没挡住；K_D=0表示背后一点波浪都没有，完美遮蔽。在实际工程中，防波堤不是无限长的，波浪会从开口或者堤头绕射进去。比如你站在一个L形的防波堤背后，虽然你看不到外海，但波浪能量还是会传过来让你感觉到摇晃。试着在模拟器里把“开口宽度B”调小，你会看到港湾内部大部分区域的K_D值（颜色）立刻变深，表示波高变小了。
🙋
诶，真的吗？那波浪的周期长短也会有影响吗？感觉长波更容易绕进来？
🎓
你的直觉是对的！波长越长，衍射（绕射）能力越强，就越容易“拐弯”进入遮蔽区。这背后的关键物理量是波长$L$，它由周期$T$和水深$h$共同决定。在我们的模拟器里，你可以固定水深，单独拖动“波浪周期T”的滑块。你会发现，当周期变长（波长变长）时，港湾内部那些蓝色（低波高）的区域会明显缩小，橙红色（高波高）的区域会扩大，这意味着更长的波浪能更深入地传播到港湾内部。
🙋
原来如此！那“入射角”这个参数是干嘛的？如果波浪不是正面过来，而是斜着冲向防波堤开口，情况会怎样？
🎓
这个问题非常关键！在实际的港口设计中，主导波浪方向往往不是正对着开口的。改变“入射角θ”，就相当于改变了波浪攻击防波堤和开口的角度。比如，当波浪斜射时，衍射波的传播模式会变得不对称。你可以在模拟器里将入射角从0度（正对）调到30度或60度，然后观察彩色分布图的变化。你会看到高波区（红色）的分布会明显偏向一侧，而不是对称地分布在开口中轴线上。这直接影响到港口内哪些泊位会更颠簸，是设计时必须考虑的因素。
物理模型与关键公式
本模拟器的核心是基于Sommerfeld衍射理论的Penney-Price简化解，用于计算半无限长防波堤背后的波浪衍射系数 $K_D$。它描述了由于衍射作用，在几何阴影区内某点的波高与原始入射波高的比值。
$$K_D(\xi) = \frac{1}{2}\left|F(\xi_+) + F(\xi_-)\right|$$
其中，$F(\xi)$ 是复数的菲涅尔积分。$\xi_+$ 和 $\xi_-$ 是与计算点的位置坐标 $(x, y)$、波长 $L$、衍射角 $\beta$ 相关的无量纲参数。这个公式优美地将复杂的波动衍射问题，转化为对标准积分函数的计算。
波浪的能量传输能力用波能通量来描述，它是评估港口内波况和静稳度的另一个重要指标。单位宽度波峰线所传递的平均功率为：
$$P = \frac{\rho g^2 H^2 T}{32 \pi}$$
这里，$\rho$ 是海水密度（约1025 kg/m³），$g$ 是重力加速度，$H$ 是波高，$T$ 是波浪周期。在模拟器中，你看到的“波能通量”可视化，就是基于此公式，结合了各点的衍射波高 $H = K_D \cdot H_0$ 计算出来的，它能更直观地显示能量在港湾内的传播路径和集中区域。
现实世界中的应用
港口总体规划与防波堤设计：在港口设计的初始阶段，工程师利用此理论快速评估不同防波堤布局（如开口位置、宽度、走向）对港内静稳度的影响。比如，通过计算不同入射波向下的K_D分布，可以优化防波堤轴线，以最小的工程成本为关键泊位提供最佳掩护。
船舶泊稳与作业安全评估：港内波高（由K_D决定）直接关系到船舶的系泊安全、装卸作业效率以及乘客上下船的安全。例如，集装箱码头要求K_D通常低于0.3-0.4，以确保起重机能够安全高效地进行装卸作业。模拟结果可以帮助划定港内不同安全等级的作业区域。
数值模拟（CAE）的验证与边界条件设定：在使用OpenFOAM、MIKE21或BEM（边界元法）软件进行全港精细的数值波浪模拟时，Sommerfeld解析解常被用来验证复杂数值模型的正确性，或为这些模型提供可靠的入射波浪边界条件。
现有港口改造与升级：对于已建港口，当需要新增泊位或遭遇新的波浪气候模式时，可以通过此工具分析现有防波堤体系的掩护效能短板。比如，分析发现某方向长周期涌浪造成内港K_D过高，从而针对性提出加长防波堤或增设消浪设施的改造方案。
常见误解与注意事项
开始使用本工具时，有几个容易误解的地方。首先，“衍射系数K_D是波高之比，而非能量本身”。若K_D=0.5，波高虽减半，但波能量（波力与越浪量的参考指标）与波高的平方成正比，实际会降至0.25倍。设计中切勿因“波高较小”而掉以轻心，务必准确评估作用在构筑物上的力。
其次，水深h与波周期T并非独立参数。工具中虽可分别调整，但实际波浪通过色散关系 $L = (gT^2/2\pi) \tanh(2\pi h/L)$ 将水深与波长L关联。例如，周期10秒的波浪传播至水深5m的浅水区域时，波长会缩短至约50m（深水区约为150m）。这种“浅水变形”会影响衍射的扩展方式，因此使用实际数据时需注意水深与周期的匹配关系。
最后，需理解本计算基于“稳态单一波浪”的前提这一根本限制。实际海况为不规则波，且具有方向谱。因此，即使工具计算出K_D=0.3的静稳区域，仍可能有其他方向的波成分传入。在实际工程中，通常以此结果为基础，再通过考虑不规则波和方向分布的更高级数值模拟进行验证。