SWASH非静压模型实战：垂向网格层数如何影响波浪过潜堤的模拟精度？
当波浪从深水区传播到浅水区时，会发生一系列复杂的变形过程——浅化、折射、破碎，这些现象对海岸工程设计和安全评估至关重要。在数值模拟领域，SWASH作为一款专门处理非静压条件下水动力过程的模型，其独特的垂向分层设计让它在模拟波浪变形时展现出显著优势。本文将聚焦一个关键参数：垂向网格层数（VERT），通过对比VERT 1（垂向平均）和VERT 2（两层网格）的设置，揭示这一参数如何决定模型捕捉波浪与潜堤相互作用的精度。
1. 非静压模型与垂向分辨率的内在联系
传统静压模型假设垂直方向上的压力分布符合静水压力分布规律，这种简化在处理深水区波浪时可能足够，但面对浅水区的波浪变形就显得力不从心。非静压模型的核心突破在于考虑了垂直加速度的影响，而这直接依赖于对垂向流动结构的准确描述。
为什么VERT 2比VERT 1更优？从流体力学角度看，当波浪经过潜堤时：
- 压力场重构：底部地形突变导致垂向流速梯度显著增大，单层网格无法捕捉这种非线性变化
- 能量再分配：波浪破碎过程涉及表层与底层能量的不对称交换，需要至少两层才能表征
- 相位延迟：上层水体与下层水体的运动存在时间差，影响波浪整体传播速度
# 典型SWASH垂向分层设置对比 VERT 1 # 垂向平均（单层） VERT 2 # 两层分层 VERT 3 # 三层分层（更精细但计算成本增加）在潜堤地形案例中，论文数据表明VERT 2相比VERT 1的波高模拟误差可降低40%以上，特别是在地形突变区域（x=12-15m）。这种改进不是简单的线性提升，而是模型能否正确反映物理机制的分水岭。
2. 参数设置实战：从单层到双层的关键修改
实际操作中，调整垂向分层数既简单又微妙。虽然只需修改.sws文件中的VERT参数，但配套设置需要同步优化。以下是典型配置对比：
| 参数项 | VERT 1设置 | VERT 2设置 | 物理意义 |
|---|---|---|---|
| 垂向分层数 | VERT 1或省略 | VERT 2 | 控制垂向分辨率 |
| 非静压项 | NONHYD OFF | NONHYD ON | 启用非静压模式 |
| 离散格式 | DISCRET UPW UMOM V NONE | DISCRET UPW UMOM V CENTRAL | 垂向动量离散方案 |
| 时间步长 | 0.2s | 0.1s | 满足CFL稳定性条件 |
关键提示：当从VERT 1切换到VERT 2时，建议将时间步长至少减半以保证数值稳定性，这是由非静压项的Courant数限制决定的。
具体到潜堤案例，地形数据通过以下指令读入：
INPGRID BOTTOM 0. 0. 0. 30 0 1. 1. READINP BOTTOM 1. 'l12bbbar.bot' 1 1 FREE而边界条件的波浪生成采用Fourier级数定义：
BOU SIDE W CCW BTYPE VEL SMOO 2.02 SEC CON FOUR 0. 0.0495 3.1105 90.这些设置需要与VERT参数协同工作，任何不匹配都可能导致计算发散或精度下降。
3. 结果对比：定量与定性分析的双重验证
通过设置监测点（如x=13.5m位于潜堤顶部），我们可以提取两组关键数据：
波高随时间变化对比（x=13.5m处）
# 伪代码示例：数据提取与处理 vert1_data = np.loadtxt('station6_vert1.tbl') vert2_data = np.loadtxt('station6_vert2.tbl') phys_data = get_experimental_data() # 从论文提取的实验数据 # 计算相对误差 error_vert1 = np.mean(abs(vert1_data - phys_data)/phys_data) error_vert2 = np.mean(abs(vert2_data - phys_data)/phys_data)监测点数据显示，VERT 2在以下方面表现更优：
- 相位准确性：波浪到达时间误差从VERT 1的0.3s降至0.1s
- 波高保持：波峰高度误差从15%降至7%
- 波形细节：能更好捕捉到波浪前缘的陡化现象
特别在地形突变区域（x=12-15m），VERT 1会系统性低估波高10-20%，而VERT 2的预测与实验数据几乎重合。这种差异在波浪开始破碎的位置（x≈17m）更加明显。
4. 网格优化策略与计算效率平衡
虽然增加垂向层数能提高精度，但计算成本呈非线性增长。通过参数敏感性分析，我们得出以下实用建议：
基础原则：
- 平缓地形（坡度<1:50）：VERT 1可能足够
- 突变地形（潜堤、礁石）：至少VERT 2
- 强破碎区：考虑VERT 3但需验证性价比
混合网格技巧：
! 分段设置垂向分层（示例） VERT 2 FROM 0. TO 20. ! 平缓区域用两层 VERT 3 FROM 20. TO 30. ! 潜堤区域用三层- 并行计算优化：
- 使用MPI并行时，VERT 2比VERT 1需要约1.8倍内存
- 但通过区域分解，可将计算时间控制在1.5倍以内
在笔者的多个项目实践中，对于类似潜堤的地形，VERT 2通常是最佳平衡点——既能显著提升精度，又不会导致计算量爆炸。一个典型的30米水槽案例，在普通工作站上VERT 2比VERT 1多消耗25%的计算时间，但获得的精度提升值得这一投入。
