SWASH模型实战：从线性波模拟到参数文件深度解析

1. SWASH模型入门：从理论到实战的第一课

第一次接触SWASH模型时，我被它强大的非静压模拟能力吸引，但面对密密麻麻的参数文件却一头雾水。经过多次实践，我发现理解这个模型的关键在于掌握参数文件的编写逻辑。SWASH作为Delft大学开发的水动力模拟工具，特别擅长处理波浪传播、近岸水动力等复杂场景。与常见的静压模型不同，它能更精确地模拟波浪的非线性特征。

对于初学者来说，线性波水槽算例是最佳的入门选择。这个经典案例能帮助我们理解模型的基本工作原理，就像学编程时先写"Hello World"一样。我建议从官网下载最新版本，目前稳定版是4.0.1。安装过程很简单，Windows用户直接运行exe文件，Linux用户需要编译源码。记得把可执行文件路径加入系统环境变量，这样在任何目录都能调用swashrun命令。

2. 参数文件深度解析：手把手教你配置.sws文件

2.1 文件头与基础设置

每个SWASH模拟都从参数文件(.sws)开始。打开文本编辑器，第一行必须是PROJ指令。我习惯用单引号包裹项目名称，比如：

PROJ 'MyFirstWave' '001' 'Linear wave test'

名称限制在16字符内，运行编号用于区分相同项目的不同算例。接下来是MODE指令，决定模拟类型：

MODE DYN ONED

DYN表示动态模拟，ONED指一维水槽。如果是三维模拟就用TWOD。

2.2 网格设置的艺术

网格配置直接影响计算精度和效率。CGRID指令定义计算域：

CGRID 0. 0. 0. 240. 0. 360 0

这里设置了240米长的水槽，分为360个网格，Δx=0.67米。垂向分层用VERT指令：

VERT 10 10. 10. 10. 10. 7.5 7.5 6. 6. 5. 5.

表示10层σ坐标，从水面到底部逐渐加密。地形数据通过INPGRID和READINP读取：

INPGRID BOTTOM 0. 0. 0. 1 0 240. 1. READINP BOTTOM 10. 'flat_bottom.bot' 1 0 FREE

3. 边界条件设置：让波浪听话的关键

3.1 波浪入射边界

左边界设置波浪入射条件：

BOU SIDE W CCW BTYPE WEAK SMOO 1.7896 SEC CON REG 0.02 3.5791

WEAK表示弱反射边界，SMOO是1.79秒的渐变时间，CON REG定义波高0.02米、周期3.58秒的规则波。实测发现渐变时间太短会导致初始扰动。

3.2 吸收边界处理

右边界需要消除反射波：

BOU SIDE E CCW BTYPE RADIATION SPON RI 40.

RADIATION表示辐射边界，SPON设置40米长的海绵层。我曾对比过不同厚度，40米对波长36λ的波吸收效果最佳。

4. 数值方法选择：精度与效率的平衡

4.1 非静压模式核心设置

启用非静压模拟必须声明：

NONHYDROSTATIC BOX PREC ILU

BOX表示采用Keller-box格式，ILU是压力求解器。对于线性波，这些默认设置就足够。

4.2 时间步长控制

时间积分参数很关键：

TIMEI 0.2 0.5 COMPUTE 000000.000 0.025 SEC 000112.000

CFL数范围0.2-0.5，实际步长0.025秒，总时长112秒。建议先用大时间步测试，再逐步缩小。

5. 结果输出与可视化技巧

5.1 测点布置策略

输出水位时间序列：

GROUP 'GAUGE' 120 120 1 1 TABLE 'GAUGE' NOHEAD 'output.tbl' TSEC WATL OUTPUT 000000.000 0.0125 SEC

在x=120米处设置测点，输出间隔0.0125秒。多位置对比能验证波传播准确性。

5.2 MATLAB后处理实战

读取.tbl文件的MATLAB脚本：

data = load('output.tbl'); time = data(:,1); eta = data(:,2); plot(time, eta); xlabel('Time (s)'); ylabel('Surface elevation (m)');

无量纲化处理能更好展示波形：

T = 3.5791; % 波周期 a0 = 0.01; % 波幅 plot(time/T, eta/a0);

6. 常见问题排查指南

第一次运行常遇到路径错误，确保.sws和.bot文件与swash.exe同目录。如果出现"NaN"结果，检查时间步是否过大。我遇到过边界条件设置错误导致波形畸变的情况，这时需要逐步检查每个参数。模型手册是最好帮手，按F1键能快速查询指令说明。

记得保存不同版本的参数文件，我用"v1_","v2_"前缀区分。当模拟复杂场景时，可以先用粗网格快速测试，再逐步细化。SWASH对硬件要求不高，我的笔记本(i5处理器)跑这个算例只需2分钟。