案例4 - 黄海沿岸风生波仿真模拟
1. 案例介绍
本案例体验利用格物云CAE三维海洋波浪模块,搭建黄海海岸区域的三维波浪模型,进行风生波仿真模拟。
2. 前处理操作
2.1 矢量数据
将研究区域的坐标系设置为32651。
利用画轮廓工具,描绘【黄海】轮廓,命名保存。 利用标点工具在研究区域内选择感兴趣的点位,并为其命名。
完成后,如下图所示:
2.2 生成几何
点击几何,点击生成几何,在弹出的窗口中主网格轮廓选择研究区域,点击生成几何。
点击几何分组,点击添加,在分组位置下方的方框里选中想要的线段,或者用鼠标靠近想要选中的线段,被靠近的线段会变成绿色,点击后该线段被选中并变成红色。
点击图中线段,设置为wave。
点击图中的线段,设置为分组in。
2.3 生成网格
点击网格,点击生成网格,在弹出的窗口中配置主网格和子网格的尺寸,可参考图中尺寸,保存。配置完成后点击生成网格,稍等片刻后,可以看到图中生成了网格。
接下来对网格进行地形赋值。点击左侧导航栏的【网格赋值】,点击添加按钮,会弹出赋值窗口。
- 赋值类型:高程
- 名称:高程
- 赋值边距:50m
- 赋值范围:插值区域
- 数值数据类型: 常量
- 赋值数据: -10
点击应用,赋值后网格如下图所示:
2.4 网格设置
垂向网格层数设置为5层
3. 水动力模型参数设置
点击水动力设置-模型菜单。
3.1 初始条件
点击初始条件,设置为均匀初始场,初始水位0m。
3.2 边界条件
对于边界条件,设定入海口的流量边界条件。流量设置为1700m3/s,施加位置设置在边界【in】上。
3.3 底部摩擦
假设整个计算域的底部摩擦均匀,选择曼宁摩擦系数0.03。
3.4 湍流
水平和垂向选择均一湍流模型。
3.5 计算时间
在计算时间选项卡中设置计算步长和计算时间。为保证计算稳定,这里设置时间步长为10s,模拟总时间14400s。
3.6 水动力结果配置
在结果界面可以勾选水深、水位、速度等。并设定结果输出频率为6时间步/帧,即每1分钟一个输出结果。
4. 海洋波浪模型参数设置
点击海洋波浪设置中的模型菜单。
4.1 初始条件
点击初始条件,设置模型初始条件,设置初始波浪场状况 。具体设置如下⾯两张图所⽰。
4.2 边界条件
设置波浪边界的输⼊形式,参照下图设置。 其中,波浪频率,以及波浪主⽅向,数值敏感性较⾼。将设置好的边界条件施加到【wave】边界上。
4.3 源汇项
点击源汇项菜单,开启考虑风的影响。
- X方向风速为5m/s,Y方向风速为-10m/s
- 风生波:WAM cycle3 模型,风生波系数填写1.2
- 线性波生长选择Cacaleri&Malanotte-Rizzoli(1981)
- 强逆流耗散:WESTHUYSEN公式
其余参数保持默认
4.3 波浪设置
- 方向数:12
- 最小频率:0.1
- 频率比率:2.2
- 频率数:20
其余模型参数设置保持默认即可。
4.4 计算时间
在计算时间选项卡中设置计算步长和计算时间。为保证计算稳定,这里设置时间步长为10s,计算步数为1440步,即模拟总时间14400s。
4.5 结果配置
在结果界面可以勾选有效波高、海底面高程、波浪能量、平均波浪方向、平均波向离散度。并设定结果输出频率为6时间步/帧,即每1分钟一个输出结果。
5. 仿真计算
全部配置完成后,点击仿真计算菜单,选择耦合计算,打开波浪驱动作用,耦合周期为1。选择32核,点击开始计算,等待完成计算。
6. 结果展示
6.1 波浪结果云图
可以查看波浪结果如有效波高、平均波浪方向(如下图)。
可以观察到平均波浪方向逐渐在风作用下发生转变。
6.2 点位结果
点击点位结果可以看到在不同的点平均波浪方向的不同。