案例3 - 长江入海口泥沙输运仿真

1. 案例介绍

本案例使用格物云CAE二维泥沙模块,搭建长江入海口的水沙耦合仿真模型,模拟分析长江入海口的泥沙输运现象。

2. 操作说明

2.1 矢量数据

将研究区域的坐标系设置为32651。

利用画轮廓工具,描绘长江及入海口轮廓,命名河道并保存。

利用画轮廓工具分别描绘三个主要岛屿,命名并保存。

编辑矢量,在河道上游、岛屿轮廓、岛屿附近的河道等位置适当添加点加密,保存。

再利用画轮廓工具,描绘一个插值区域,命名并保存。

利用矢量点工具,在河道中点几个后续想分析的点位,命名并保存。

或者本地上传准备好的矢量文件(点击下载)。

2.2 生成几何

点击几何,点击生成几何,在弹出的窗口中选择主网格轮廓为河道,分别添加三个岛屿为子网格轮廓,并把岛屿设为不可淹没区域。点击生成几何。

点击几何分组,点击添加,在分组位置下方的方框里选中想要的线段,或者用鼠标靠近想要选中的线段,被靠近的线段会变成绿色,点击后该线段被选中并变成红色。选中河流上游入口,并重命名为river,保存。

用相同的操作选中河流下游右侧的河道边界,并重命名为sea,保存。

2.3 生成网格

点击网格,点击生成网格,在弹出的窗口中配置主网格和子网格的尺寸,可参考图中尺寸。配置完成后点击生成网格。

接下来对网格进行地形赋值。点击网格赋值,点击添加按钮,会弹出赋值窗口。

对研究区域赋值:

  • 赋值类型:高程
  • 赋值边距:可选为1
  • 范围:选择 河道(main)
  • 数值类型: TIF/XYZ-value型数据
  • 上传地形数据:ytz_river_estuary.xyz
  • 点击应用

赋值后网格如下图所示。

2.4 水动力模型参数

点击水动力设置-模型菜单。

2.4.1 初始条件

点击初始条件,设置为均匀初始场,常数水位2m。

2.4.2 边界条件

对于边界条件,设定河道上游为流量边界,入海口为潮汐边界。 边界一:边界类型为流量边界,施加位置选择river,给出常值流量29000m^3^/s。

边界二:边界类型为潮汐边界,施加位置选择sea,上传水位边界和速度边界文件(点击下载文件:uv_DTv5.outh_DTv5.out)。时间可按个人喜好设置,比如直接点击日历和时钟栏内的今天和此刻,坐标系参数为UTM北半球51区:

2.4.3 底部摩擦

假设整个计算域的底部摩擦均匀,选择曼宁摩擦系数0.025。

2.4.4 潮汐

点击潮汐,查看潮汐设置中,潮差校准系数为1,海平面为0m。

计算格式如图设置,其余模型参数设置保持默认即可。

2.4.5 计算时间

在计算时间选项卡中设置计算步长和计算时间。为保证计算稳定,这里设置时间步长为10s,模拟总时间为10天,即864000s。

2.4.6 水动力结果配置

在结果界面可以勾选水深、水位、高程、速度、速度标量等。并设定结果输出频率为360时间步/帧,即每一小时一个输出结果。

2.5 泥沙输运模型设置

点击泥沙设置-模型菜单。

2.5.1 泥沙性质

添加一种粘性泥沙,颗粒半径为0.00002m,初始占比为0.33,设置悬沙初始浓度为0.15g/l,其他参数保持默认。

添加一种非粘性泥沙,颗粒半径为0.00003m,初始占比为0.67,悬沙初始浓度为0.35g/l,其他参数保持默认。

2.5.2 输运机制

激活悬移质输沙,选择Soulsby & van Rijin 公式。

2.5.3 床层模型

点击加号,设置1层1m的初始床层,床层孔隙率设为0.4。

2.5.4 粘性泥沙输运

点击加号,设置1层的侵蚀临界底床剪切应力为1N/m^2^,Krone-Partheniades侵蚀定律常数为0.00005。

2.5.4 泥沙边界条件

泥沙入流边界设为恒定悬沙浓度,对应添加两类泥沙的悬沙浓度,分别为0.15g/l和0.35g/l。

2.5.5 泥沙结果配置

泥沙的计算结果输出可以勾选底床演变、底床剪应力、每类泥沙的质量浓度等。

其他参数保持默认。

2.6 仿真计算

全部配置完成后,点击仿真计算菜单,选择水沙耦合求解,点击提交计算,开始仿真计算,计算运行时间约3~4分钟。

3. 结果展示

3.1 水动力结果云图

点击水动力结果云图,即可查看水深、水位、流速、流量、两类泥沙浓度等传播情况。如下动图为水位随时间变化的过程,设置色阶范围为-3~3。

3.2 泥沙结果云图

悬沙浓度的结果储存在水动力结果中,点击查看粘性和非粘性泥沙浓度随时间的变化,设置色阶范围为0~0.5,可以观察到两类泥沙在水体中迁移和分布的差异。

点击泥沙结果云图,可以看到底床剪应力随时间的变化,色阶范围为0~10。

点击查看底床演变,直接选择结果的最后一帧,配色方案可以选择erdc-icepeach,显示范围选-0.5~0.5,可见长江口河道冲淤趋势。