狭长矩形湾的风潮计算研究.pdfVIP

狭长矩形湾的风潮计算+ 曾一非 (上海交通大学船舶海洋与建筑工程学院，上海200030) 摘要应用多层模型和有限差分法计算三维N-S方程，采用交错网格技术，分别对一狭 长等深封闭海湾和半封闭海湾在常风应力作用下的风潮运动做了数值模拟计算。矩形湾水深 10m，假设整个流场上方均匀作用有一定常风应力，其风速恒为lom／s，垂向分为5层，水平方 向的网格间距为1000m，得到了其线性和非线性、定常和非定常的结果。并与对应的浅海风暴 潮线性理论和非线性理论模型进行了计算比较和讨论，给出了定常状态下与沿湾长风暴潮位 线性理论与非线性理论的结果比较，以及非定常模型下风暴潮位时历曲线的结果比较。结果 显示它们之间有很好的一致，但浅海风暴潮理论某些因子的选取值得探讨。 关键词风潮，风应力，封闭湾，半封闭湾，数值计算 1前言 风暴潮是一种灾害性的海面异常升高现象，可由热带气旋和温带气旋，以及寒潮引起，在 我国一年四季都有可能发生，因此尤其深受其害。导致风暴潮产生的两个重要动力因子是海 面上作用有强风和气压骤变，对这两者的计算也就直接影响着风暴潮的计算与预报精度。但 在浅海区域，相比气压变化，风应力起着主要的动力作用，这时可不计气压的变化作用，这种只 考虑纯风作用下的风暴潮有的学者称其为“风潮”，并在某些近似假设的前提下针对简单模型 提出了浅海风暴潮理论，得到了一些理想情况下的理论解。实际情况的风暴潮由于地形复杂 以及气压场和风场复杂等而难于使用以上的近似理论解来描述，但它的物理概念清晰，是认识 实际风暴潮的基础，在难于得到实际风暴潮数据的情况下，利用它与数值计算做些对比计算可 相互检验模型和计算方法的有效性，进而提高计算的精度。描述风暴潮运动的基本方程是N— S方程，这是一个三维方程组，有许多方面的计算将它简化为二维方程组进行，但为了描述它 的三维流场结构和探求它的发生发展机制，则有必要做三维流场计算。本文应用多层模型和 有限差分法计算该三维N—S方程，采用交错网格技术做差分计算，为与浅海风暴潮理论做比 较，计算对象分别取为一狭长等深封闭海湾和半封闭海湾，计算它们在常风应力作用下的风潮 运动。 + 国家教育部留学回国人员科研启动基金资助课题。 ·182· 数值模式和解析解 2．1 基本方程组 平均海面上取直角坐标系0一划2，2轴取向上为正。针对其垂向尺度相比其水平尺度小 许多的特点，采用静水压假定。不计入科氏力影响，其控制方程式参见文献[2]。 海面风应力同大多数文献采用二次律的经验公式，为％=C驴。叫I硼f，叫为风速，R为空 气密度，G为风拖曳力系数。有文献[3]说明G的取值与海面粗糙度有关，但在本次计算中取 为常数。 2．2 数值计算 水平差分采用交错网格技术，垂向按多层模型进行分层[2]，共分为5层。初始条件取善一 “=才=硼=o进行计算，壁面边界按FrerSlip边界条件确定，取U·竹=O。计算网格间距△z =匈=1 ooom，计算对象为一长215km，宽15km，深10m的狭长矩形湾，封闭海湾和半封闭 海湾分别对应左端是闭边界和开边界。取风速硼=10m／s，风拖曳力系数已一o．0026，空气 密度一1．29kg／m，海水密度P=1 025kg／m，水平及垂向涡粘系数A^=50m，A。=o． 001 m2／s，时间计算间隔为△=10s。 2．3 解析解‘1] 对于湾长L的封闭海湾，在常风应力％作用下定常风暴潮位的线性解析解为 o r=警(z一专) 2—矿oz—i) 非定常模型的线性解析解分为忽略惯性效应与考虑惯性力作用的两个方面，即 f一≮笋曼面‰·‘2一—矿。幺…西丽。 ”eXp[逝铲邙cos学， 9。一—矿。囊…西丽。善一一％萨。囊丽‰· hc。洲]}c。s堕产 卜去eXp[罢s‰t 式中如彩西乎匿巧。