水动力模型软件Delft3D操作过程与实例分析.docx
- 文档编号:371764
- 上传时间:2022-10-09
- 格式:DOCX
- 页数:12
- 大小:438.35KB
水动力模型软件Delft3D操作过程与实例分析.docx
《水动力模型软件Delft3D操作过程与实例分析.docx》由会员分享,可在线阅读,更多相关《水动力模型软件Delft3D操作过程与实例分析.docx(12页珍藏版)》请在冰豆网上搜索。
水动力模型软件Delft3D操作过程与实例分析
水动力模型软件Delft3D操作过程与实例分析
基于Delft3D-Flow模块的三维模型垂向在网格使用的是笛卡尔Z坐标,是基于二维浅水方程对不可压缩粘性流体求解NavierStokes方程,在浅水和方程假设下。
在垂直动量方程中,忽略了垂直加速度,导致了静水压方程。
在3D模型中,垂直速度是从连续性方程计算出来的。
在有限差分网格上求解与一组合适的初始边界条件相结合的偏微分方程组
连续方程:
垂向二维平均连续性方程如下
其中,Q表示单位面积上的源、汇通量,以及降水和蒸发的影响。
水平方向上的动量方程
在ξ和η方向上的三维动量方程如下:
和
为垂向紊动粘性系数,与水体的动粘度和三维湍流闭合模型计算得到的三维紊动粘性系数有关。
Pξ和Pη代表两个方向的压力梯度,动量方程中Fξ和Fη代表水平雷诺应力的非平衡量,Mξ和Mη代表源汇项。
泥沙输运方程
对于河口海岸中溶解性物质、盐分、热量的输运,delft3d-flow解算了三维平流-扩散方程,并考虑了一阶衰减过程。
输运方程在水平曲线正交网格和垂向坐标系以守恒形式表示。
其中DH是水平扩散系数,DV是垂向扩散系数,d代表一阶衰减过程。
1.水动力模型构建过程
1.1网格搭建
1.1.1设置工作空间
工作空间(workingdirectory)是Delft3d运行的工作文件夹,选择正确的工作空间可以提高计算效率,减少失误。
该课程论文以Delft3D自带的实例操作为主,因此在这里选择文件夹:
D:
\Delft3D4.02.01\tutorial\flow
1.1.2导入陆地边界
陆地边界相当于海岸线,该处的作用主要为规定划分网格和模拟区域。
在菜单栏选择File-attribute-openlandboundary,打开
图1陆地边界示意图
1.1.3绘制曲线(spline)
曲线是沿着陆地边界绘制而成,曲线的主要功能是生成网格,因此在绘制曲线的时候,曲线的范围要比陆地边界的范围大。
操作步骤:
在菜单栏点击Edit-spline-edit,首先沿陆地边界拖动鼠标,在陆地边界外围绘制曲线,绘制完曲线右击鼠标完成一条曲线;然后在网格内部绘制等距曲线,右击结束。
结果如图所示
图2spline曲线示意图
1.1.4根据曲线生成网格
在菜单栏单击operation-growsplineintogrid,这是网格会按照曲线自动生成。
需注意的是,此处生成的网格并不能直接用于模型计算,还需要对网格进行调整。
自动生成的网格如图所示。
1.1.5网格后处理
1.1.5.1网格移动、扩大与缩小
按住ctrl键然后移动鼠标可以自由移动网格;鼠标滑轮可以对网格进行自由缩放。
1.1.5.2网格加密与稀疏
在菜单栏单击operation-regulargridcoarseness-refine为加密网格,derefine为网格稀疏。
加密的行列数可以在setting里设置
1.1.5.3网格正交化
网格正交化主要针对的是网格与网格之间的垂直程度,单击菜单栏operation-orthogonalise,可以直接对自动网格进行正交化处理。
1.1.5.4网格平滑度处理
网格平滑度是指网格与网格之间的长宽比,最完美的情况是网格内部所有小网格大小完全一致。
操作步骤:
Edit-regulargrid-blocksmooth,单击之后没有变化,此时,分别单击需要进行平滑处理网格的两端,然后右击鼠标,就会自动进行平滑处理。
图3初步生成网格示意图
1.1.6网格精度检验
1.1.6.1正交性
网格正交性以余弦值作为参考值,其最合适的取值范围为<0.02,但是最大不能超过0.05,否则计算结果会出错。
1.1.6.2平滑度
平滑度的取值范围为10%~20%。
经过以上步骤调整后,得到网格即为模型计算所需要的网格,完成之后结果如下图
1.1.7
导出网格
图5完成网格示意图
网格导出功能在file菜单栏,单击export,选择保存路径直接保存即可。
导出网格包含两个文件,分别是.grd和.enc。
1.2水底地形文件
水深数据的生成是在生成网格的基础上,将离散水深点通过插值离散的方式赋值到网格,进而以面的形式表现出来。
因此,要生成水深数据需要两种基础数据,分别是网格(前期以做好的网格)和离散水深点。
制作水底地形文件用Delft3D中的Quickin工具,首先打开Quickin。
1.2.1导入陆地边界、网格与水深点
陆地边界在这里作为一个辅助边界进行展示,仅作为可视化的形状。
导入方法:
File-attributes,导入陆地边界选择importlandboundary(选择文件mar_01.ldb),导入水深点选择imprtsample(选择文件mar_02.xyz)。
之后选择相应文件确定即可。
导入网格:
File-importgrid,选择以生成的网格即可。
(选择网格文件mar_02.grd)。
图6离散水深点示意图
1.2.2水深点插值形成水底地形文件
导入完成的陆地边界、网格和水深如图所示。
之后对水深点进行插值即可得到水底地形。
操作:
菜单栏operations-triangularinterpolation。
图7插值水深成品示意图
1.2.3导出水深
完成插值后需要对水深进行导出,以便于输入模型。
操作:
File菜单栏-export-exprtdepth。
导出的文件为.dep格式的文件。
需要注意的是,一个网格只对应一个水深数据。
1.3输入模型参数
将之前做好的网格与水深数据以特定方式输入模型中,作为一个模型运行的基础;输入模型的开边界作为整个模型运动的驱动因素;规定模型的计算间隔与模拟结果输出的时间间隔;除此之外,还需要对模型中的参数进行计算调整,以满足验证数据的需要,以验证模型的精确性。
模型文件及参数输入界面如下图所示。
1.3.1导入网格与水深(domain)
首先导入模型网格与水深数据。
图8Flow模块设置主界面
网格(grid)主要输入.grd与.enc两种格式的文件;除此之外还需定义所在纬度,以考虑科氏力对模型的影响;由于现在构建的是一个二维模型,因此在numberoflayers一栏输入“1”,代表二维模型。
见左图。
水深(bathymetry)则主要输入之前生成的水深文件。
直接在file菜单下选择相应的.dep水深文件即可。
完成见右图。
图9网格与水深导入位置及其使用
1.3.2
设置模型计算时间
图10模型时间参数设置
在模型计算时间方面,将计算间隔设为5分钟,区时采用默认为0。
计算长度为1天,为1990年8月5日至1990年8月6日,具体时间设置见下图。
1.3.3设置模型初始值
模型初始值代表模型在开始计算之前的状态,在这里设置水面初始值为0,代表模型为冷启动方式。
在冷启动方式下,模型在开始时的变化可能会比较剧烈,但是一段时间之后会逐渐趋于稳定。
1.3.4设置模型边界
在边界条件设置上采用动边界条件作为模型的驱动条件。
边界条件为根据已有的潮汐数据对水位变化进行调和而成。
具体边界条件类型以及设置条件见下图。
1.3.5
设置物理参数(physicalparameters)
图11模型开边界条件设置
模型运动需要设置的物理参数比较多,如重力加速度、水体密度、水底粗糙度、边界粗糙度、水体粘滞系数等。
因为该模型为Delft3D的示例,所以在物理参数方面不做调整,直接采用Delft3D给出的默认值即可。
1.3.6
设置观测点与观测剖面(monitoring)
图12模型观测点与观测剖面设置界面
观测点与观测观测剖面的设置主要针对的模型的结果,在对网格节点上添加观察点,可以对指定位置的模型结果进行定量观测,并将结果作为单独的文件进行存储,这在很大程度上可以在后期通过对特定点位置的结果进行分析,或者与模型验证数据进行比较,从而调整模型参数,提高模型结果的精度。
具体的观察点设置结果如下图所示。
1.3.7模型结果输出设置
模型结果的输出有两种形式,分别是mapresults与historyresults。
historyresults主要是基于上一步设置的观察点与观测剖面而言的结果,mapresults则是指在网格的基础上形成的面状动态显示的结果。
两种结果都是模型输出的重要形式,各有侧重点,都是结果输出的重要部分。
在输出设置方面,需要对输出结果的时间间隔与时间范围进行定量的规定,最终模型输出结果的设置界面如下图所示。
图13模型输出结果设置
下图为最终的模型结构设置,包括水位开边界、水深、网格、观测点与观测剖面等信息。
图14模型最终设置完成后的示意图
2.实例结果分析
第一部分主要是对模型所需数据文件的制作,所需参数以及如何将模型所需数据输入Delft3D进行计算。
模型在经过一段时间的运行计算之后,会生成相应我们所需要的模型输出结果,包括流速、水位变化等方面的数据,之后要做的就是对所输出的结果进行分析。
1.流速分析
流速(flowvelocity)是指液体单位时间内的位移。
质点流速是描述液体质点在某瞬时的运动方向和运动快慢的矢量。
其方向与质点轨迹的切线方向一致。
由于在该模型设置时使用的二维模型,所以在结果输出时只有水深平均速度(depthaveragevelocity)。
图1为在模拟时间内的流速分布图,整个模拟区域的流速分布图。
从图中可以看出,该模型模拟的是一个潮汐通道的水动力变化情况。
水流在从外海向近岸传输过程中,水流速度在潮汐通道处始终保持高流速运动,最大流速可以达到90cm/s,其最小速度也在40cm/s以上。
另外,由于大部分流速在潮汐通道被拦截下来,所以近岸处的流速普遍较低,基本除在10~20cm/s的范围内;而外海的流速变化具有周期性,其流速变化可以周期性的增大或减小,这与其水深与开边界处的参数设置有关。
图14流速变化的时空分布图
图2为对四个点进行定量观测后所获得的流速变化曲线,从四副图来看,其流速变化具有明显的周期性变化特点,其中变化幅度最大的是Obs4位置点,其最大流速可达100cm/s,最低可达10cm/s;Obs2的变化幅度也较大,但是小于Obs4;其次为Obs3,最后是Obs1。
通过对观测点位置的分析可以看出,该模拟区域该模拟区域东北侧的流速比东西侧的流速较大,潮汐通道处的流速比北侧的流速要大。
图15定点观测下的流速随时间变化图
2.水位分析
水位是指对水面高度周期性变化的描述,水位是反映水体水情最直观的因素,它的变化主要由于水体水量的增减变化引起的。
水位过程线是某处水位随时间变化的曲线,横座标为时间,纵座标为水位。
分别对水位的时空分布范围与定点观测水位的周期性变化进行变化分析。
图3为水位在模拟周期时间内的水位时空变化图,从图中可以看出,水位变化具有周期性的涨落趋势,水位最高是在12时左右,其最高值可达200cm;水位最低约在5时,最低值为-150cm。
其余时刻的水位变化均在-150~200之间来回摆动。
从范围来看,最高值与最低值也存在交替出现现象,当外海是最高水位时,近岸地区为低水位;当近岸地区为低水位时,外海地区为高水位。
图16模拟水位的周期变化
通过对四个地点的水位变化曲线分析发现,四地的水位变化具有同步性,在同一时间里其水位变化趋势大致相同,最高水位与最低水位的时间也基本相同,其之间具有的差异性主要在于水位的变化幅度。
Obs4点的水位变化相对较缓慢,起伏幅度比其他三个点要小,这可能与其所处位置——潮汐通道有关;另外三个外海位置点的水位变化没有太大区别。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 动力 模型 软件 Delft3D 操作过程 实例 分析