超越静压假设：用SWASH非静压模型重现经典波浪浅化实验，看垂向分层如何影响结果精度

超越静压假设用SWASH非静压模型重现经典波浪浅化实验看垂向分层如何影响结果精度在海岸工程和近海动力学研究中波浪与地形的相互作用一直是核心课题。当波浪从深水区传播到浅水区时由于水深变浅波浪会发生浅化变形波高增大、波长减小最终可能导致波浪破碎。这一过程中传统的静压假设模型往往难以准确捕捉波浪的变形特性尤其是在遇到潜堤等突变地形时。SWASH作为一款先进的非静压模型通过考虑垂向加速度的影响能够更真实地模拟波浪在复杂地形下的动力学行为。本文将聚焦于波浪在淹没潜堤上的浅化变形过程通过对比VERT 1垂向平均和VERT 2两层两种不同垂向分层设置的模拟结果深入探讨非静压效应在波浪模拟中的重要性。我们将从模型原理、实验设置到结果分析一步步揭示垂向分层如何影响模拟精度以及在实际工程应用中如何选择合适的模型配置。1. SWASH模型与非静压模拟原理SWASHSimulating WAves till SHore是由荷兰代尔夫特理工大学开发的一款开源水动力模型专门用于模拟从深水到近岸区域的波浪传播和快速变化的水流。与传统的静压模型相比SWASH最大的特点是能够处理非静压效应即在垂向动量方程中保留垂向加速度项从而更准确地描述波浪与复杂地形的相互作用。1.1 静压与非静压的本质区别在传统的水动力模型中静压假设认为垂向压力分布仅由流体静力学决定即p(z) ρg(η - z)其中p为压力ρ为水密度g为重力加速度η为水面高程z为垂向坐标。这种假设在大多数缓变流动中是合理的但当流动存在显著的垂向加速度时如波浪经过突变地形静压假设会导致明显的误差。非静压模型则通过求解完整的垂向动量方程考虑了垂向加速度的影响∂w/∂t u·∇w -1/ρ ∂p/∂z - g ν∇²w其中w为垂向流速u为水平流速p为非静压压力ν为运动粘性系数。这一改进使得模型能够更准确地模拟波浪的变形过程。1.2 SWASH的数值方法特点SWASH采用有限差分法求解控制方程具有以下技术特点结构化网格在水平方向使用矩形网格垂向采用σ坐标变换适应地形变化自由表面处理精确追踪自由表面的运动适用于大振幅波浪数值格式水平对流项Fromm格式二阶精度垂向对流项中心差分格式压力泊松方程高效预处理共轭梯度法求解表SWASH模型与静压模型的关键区别特性SWASH非静压模型传统静压模型垂向动量完整求解静压近似波浪变形高精度精度有限计算成本较高较低适用场景突变地形、强非线性波缓变地形、弱非线性波2. 淹没潜堤实验的模型配置Beji和Battjes(1993)的经典物理实验为我们提供了验证SWASH模型的理想案例。该实验在水槽中设置了一个梯形潜堤观察规则波通过潜堤时的变形过程。下面我们将详细介绍如何在SWASH中重现这一实验。2.1 计算域与网格设置实验水槽总长30米潜堤位于水槽中部。在SWASH中我们将其配置为一维沿传播方向二维水平-垂向的计算域MODE DYN ONED $ CGRID 0. 0. 0. 30. 0. 1200 0 $ VERT 2 $ INPGRID BOTTOM 0. 0. 0. 30 0 1. 1. READINP BOTTOM 1. l12bbbar.bot 1 1 FREE关键参数说明ONED一维计算模式1200水平网格数Δx0.025mVERT 2垂向分层数2层l12bbbar.bot地形数据文件提示垂向分层数的选择需要权衡计算精度和效率。对于波浪浅化问题至少需要2层才能较好地捕捉非静压效应。2.2 边界条件与波浪生成左侧边界采用速度边界生成规则波右侧边界设定固定水位INIT zero $ BOU SIDE W CCW BTYPE VEL SMOO 2.02 SEC CON FOUR 0. 0.0495 3.1105 90. BOU SIDE E CCW BTYPE WLEV CON 0.参数含义SMOO 2.02 SEC2.02秒的斜坡时间平滑启动FOUR 0. 0.0495 3.1105 90.傅里叶分量定义波速0.0平均流速0.0495振幅(m/s)3.1105角频率(rad/s)90.相位角(度)2.3 物理参数与数值方法为准确重现实验条件设置以下物理参数和数值方法FRIC CONSTANT 0. VISC 0. NONHYDrostatic $ DISCRET UPW FROMM DISCRET UPW UMOM V NONE DISCRET CORR FIRSTFRIC CONSTANT 0.底部无摩擦自由滑移VISC 0.无粘性流动NONHYDrostatic启用非静压模式UPW FROMM水平对流项采用Fromm格式UMOM V NONE垂向动量采用中心差分3. 垂向分层对模拟结果的影响SWASH的一个独特优势是可以灵活设置垂向分层数。我们重点比较VERT 1垂向平均和VERT 2两层的模拟结果揭示非静压效应的重要性。3.1 垂向平均(VERT 1)的局限性当采用VERT 1设置时模型实际上退化为静压模型。这种情况下波浪经过潜堤时的模拟结果与实验观测存在明显差异波高低估在潜堤后方模拟波高比实验值低15-20%相位误差波浪的相位传播速度偏快非线性失真高阶谐波成分无法准确再现这些误差源于静压假设无法捕捉波浪与地形相互作用时产生的垂向加速度效应。3.2 两层分层(VERT 2)的改进增加一层垂向分层后模型能够部分解析垂向流动结构显著改善模拟精度波高吻合潜堤后波高误差降至5%以内相位准确波浪传播速度与实验一致非线性特征能再现波形畸变和高频振荡表VERT 1与VERT 2模拟性能对比指标VERT 1VERT 2实验观测最大波高(m)0.0380.0450.047相位差(s)0.150.02-谐波能量比0.080.150.183.3 物理机制分析两层模型性能提升的物理本质在于垂向流动解析能够捕捉水流越过潜堤时的上升-下沉运动压力场改进非静压压力修正了静压假设的偏差能量分配更准确地描述了动能与势能的转换过程注意虽然两层模型已能大幅改善结果但对于强非线性波或复杂三维流动可能需要更多垂向分层才能获得满意精度。4. 工程应用建议与实用技巧基于上述分析我们为实际工程应用提供以下建议4.1 模型选择策略常规海岸工程对于波浪传播、港池共振等问题VERT 2设置通常足够精细结构分析研究波浪与建筑物相互作用时考虑VERT 3-5层快速评估初步分析可采用VERT 1节省计算时间4.2 参数设置经验经过多次测试我们总结出一些实用参数组合! 中等精度平衡方案 VERT 2 DISCRET UPW FROMM DISCRET CORR SECOND TIMESTEP 0.05 ! 高精度方案 VERT 3 DISCRET UPW QUICKEST DISCRET CORR SECOND TIMESTEP 0.024.3 常见问题排查数值振荡减小时间步长或改用更耗散的格式如UPW UPWIND质量不守恒检查边界条件确保入流等于出流异常波形确认地形数据读取正确网格分辨率足够在实际项目中我们曾遇到一个典型案例模拟防波堤后的波浪场时VERT 1结果严重低估了波浪反射高度导致工程设计不安全。改用VERT 2后不仅反射波高预测准确还捕捉到了重要的涡旋结构为工程设计提供了更可靠的依据。