01 研究背景

  泽布鲁日港（Zeebrugge）是比利时最重要的海港之一。泽布鲁日港附近区域经常发生强烈的潮汐流，在涨潮时潮汐振幅可达5米。由于潮汐的影响，在港口测量到不断变化的复杂涡流。研究在港口处产生的涡旋，特别是悬浮泥沙流经港口后的重新分配，不论是出于航海目的还是为了加深对港口淤积的了解，都是很有必要的。

  模型采用Smagorinsky格式计算水平涡旋粘度；采用Nezu和Nakagaw混合长度参数化模型计算垂直涡旋粘度。工程师通过对比实测的水位和流速数据来进行曼宁摩擦系数的校准，得出曼宁摩擦系数的值为0.02s/m1/3，并在靠近港口入口的地方适当增加摩擦系数来表示由堤岸的大型混凝土构件的摩擦带来的影响。IMDC的工程师使用三位水动力模块，以20s的时间步长，进行了时长15天的仿真计算。

  在模型的计算结果中，高水位的出现时间在整个潮汐周期中更晚。结果重现了ADCP测量结果中的强射流，（如图3c所示），而射流的位置更靠东。伴随着强射流的出现，同样形成了强烈的顺时针旋转的主涡旋。在主涡旋出现后1h（如图3d所示），同样地，射流也促使了较弱的逆时针旋转的次级涡旋的产生。

  根据ADCP的测量结果，在高水位到达时（如图4a所示），不再有射流流入港口，涡旋与海上流动解耦，开始遵循自身的动力学，主涡旋的大小和量级开始减少。直至高水位到达的1h后（如图4b所示），港口内只剩下逆时针旋转的次级涡旋，且次级涡旋缓缓向海港入口移动，并逐渐减弱。

  而在高水位到达约1h后，模型计算结果中仍然观察到了主涡旋和次级涡旋同时存在（如图4c所示），但随着时间的推移，主涡旋的大小迅速减弱，次级涡旋的强度较小更慢；最后主涡旋完全消失（如图4d所示），这与ADCP的测量结果是吻合的。然而在模型的计算结果中，在主涡旋的北部可以观察到一个与次级涡旋旋转方向相同的，非常微弱的三级涡旋。该三级涡旋没有被ADCP测量到，其原因可能是由于ADCP测量提供的空间分辨率有限，不足以观测到如此微小的涡旋。

  为了评估河床粗糙度对涡旋形成的敏感性，IMDC的工程师不再对港口沿岸的摩擦系数进行增大处理，采用全场恒定的曼宁粗糙度系数0.02 s/m1/3进行仿真。与基础工况相比，模拟的涡流模式发生了显著变化(图5)。首先，可以看到相同的射流发展(图5a)。然而，它的位置略有不同。因此初级涡旋的环流略微增强，而形成的次级环流则略弱于基础工况的次级环流(图5b)。因此，在使用全场统一的曼宁粗糙度系数后，主涡旋不会随着时间的推移而消失(图5c)，最终在退潮阶段在港口仍存在有两个涡流，而不是一个(图5d)，这与ADCP的测量结果是不一致的。

  B．湍流模型/涡流粘度

  为了测试水平湍流模型设置对港口模型的敏感性，IMDC的工程师采用了水平涡流粘度恒定为1 m2/s的模型来进行仿真计算。计算表明，采用恒定水平涡流粘度模型时结果比水平涡流粘度约为0.01 m2/s的Smagorinsky格式略高。高粘度系数的使用削弱了港口入口处的射流(图6a)。随着时间演进，次级涡旋产生，但仍弱于主涡旋(图6b)。因此在这种情况下，在退潮开始时，主涡旋仍留在港口(图6c)，且比次级涡旋略强。然而，随着时间演进，主涡旋最终消失，流动方向最终与ADCP的测量结果相同(图6d)。