论文部分内容阅读
本文建立了一个用于波—流共同作用下长江口平面二维悬沙计算的数值模型系统,主要由四部分组成,即考虑风、浪、潮、径流的二维复合流场模型、SWAN海浪模型、波—流相互作用的底边界层模型以及二维悬沙输运模型。采用三个不同大小的区域进行嵌套计算,三个计算区域分别为东中国海区域、长江口杭州湾区域和长江口区域。 模型系统采用广义曲线坐标系下的形式,使用高精度的自适应网格拟合复杂岸界。复合流场学模型中考虑了波浪辐射应力的作用,底摩擦由波—流相互作用边界层模型提供。为了较高精度的模拟天文潮,在开边界加上了七个主要分潮。台风天气下,在东中国海区域计算由风产生的水位变化,插值到长江口杭州湾区域的外海开边界上作为余水位,这样就同时考虑了局地风和非局地风的影响。 悬沙输运模型利用切应力方法来确定对流扩散方程中的泥沙源函数,其中的临界起动速度利用经典的泥沙起动流速公式前面增加一个局地系数得到,这个系数能反映河床底质结构及固结程度,通过系列数值试验和实测资料的统计分析确定。计算中的底部剪切应力由波流边界层模型给出波—流共同作用下的形式。泥沙颗粒絮凝沉降速度考虑了流速、盐度、含沙浓度的影响。 将模型系统应用到长江口区域,经过了大量的实测资料的验证。其中复合流场模型中的天文潮利用众多测站的实测调和常数进行验证;波浪模型首先进行了水槽实验的检验,然后在台风过程中,进行了测站有效波高过程线的比较。复合流场模型的水位过程、流速过程,以及悬沙模型的含沙量过程则经过了洪、枯季及大、中、小潮的多个站点的实测过程验证。分析流场、波浪场以及正常和台风天气下的悬沙场的计算结果,表明该模型系统能合理地反映长江河口区域的水动力场和泥沙场的分布。