基于非静压模型的波浪破碎模拟-水科学进展.PDFVIP

第２８卷 第３期 水 科 学 进 展 Ｖｏｌ．２８，Ｎｏ．３ ２０１７年５月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｍａｙ，２０１７ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１７０３０１５ 基于非静压模型的波浪破碎模拟 张景新１，２ （１． 上海交通大学水动力学教育部重点实验室，上海　 ２００２４０； ２ 上海交通大学船舶海洋与建筑工程学院工程力学系， 上海　 ２００２４０） 摘要： 波浪破碎是海岸工程所关注的关键水动力学问题之一， 而波浪破碎的数值模拟技术的研究方兴未艾。 高效 的浅水方程基于静压假定， 而通过引入动压项建立的完全非静压模型， 可成功应用于色散水波的模拟。 自由表面 的捕捉采用的垂向坐标变换技术， 较之ＶＯＦ（Ｖｏｌｕｍｅ ｏｆ Ｆｌｕｉｄ）模型， 计算效率较高。 但垂向坐标变换不能模拟大曲 率自由表面变形， 即波浪破碎过程。 对于破碎波的模拟， 一种高效的模式分裂法应用至非静压模型， 即在波浪破 碎局部水域将模型分裂为静压模型和非静压模型， 破碎波波峰附近退化至静压模型， 并持续至波浪破碎结束再恢 复为动压模型。 通过典型算例， 验证了模式分裂法的适用性及对于波浪破碎过程的模拟精度， 鉴于模式分裂法对 于波浪破碎过程的模拟未引入新的概化模型， 计算效率较高， 可应用于大尺度的海岸带波浪的变形、 破碎及越浪 的数值模拟。 关键词： 浅水方程； 非静压模型； 色散水波； 波浪破碎 中图分类号： ＴＶ１３５　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１７）０３⁃０４３８⁃０７ 静压假定模型广泛应用于地表水流运动， 因控制方程的简化， 使得模型可以应用于诸如潮汐流、 明渠流 等大尺度水流运动， 已开发出了大型商用软件， 如荷兰的Ｄｅｌｆｔ３Ｄ和丹麦的Ｍｉｋｅ系列软件。 然而对于陡变地 ［１］ 形条件下的流动或色散水波的模拟， 该假定条件已然不适用 。 这类水流运动的模拟需要突破静压假定， ［２⁃４］ 非静压模型即在此基础上发展而成 。 对于自由表面水流运动的数值模拟， 自由表面捕捉方法的精度及效率至关重要。 通过垂向坐标变换， 可精 确拟合自由表面，σ坐标变换即属于此类。 较之ＶＯＦ（Ｖｏｌｕｍｅ ｏｆ Ｆｌｕｉｄ）等精确捕捉自由表面时空变化特征的数 值技术而言， 坐标变换的计算量显著降低，从而使得模型应用于大尺度自由表面水流运动的模拟成为可能。 对 于缓变的自由表面， 该方法是适用的， 但对于波陡较大， 发生波浪破碎、 翻卷等情况， 坐标变换方法不再适 用。 基于垂向σ坐标变换模拟自由表面水流运动的数值模型对于破碎波的模拟需要发展新的模拟技术。 利用浅水方程模拟自由表面水流运动， 对于破碎波的模拟， 已经开发了若干模式法。 通过动量方程中引 入人工涡黏模型， 从而模式化破波过程中的能量耗散， 这类方法已成功应用于浅水方程、 Ｂｏｕｓｓｉｎｅｓｑ类方 ［５⁃１０］ 程 。 近年来， 有别于动量方程引入能量耗散机制， 研究发现波浪破碎形成了间断， 可类比于水跃， 而非 ［１１］ 线性浅水方程对于这类间断流动的模拟， 若数值格式具有守恒性， 则可获得满意的模拟结果 。 当波浪演 化至破波点时， 在波前局部水域将Ｂｏｕｓｓｉｎｅｓｑ方程的色散项关闭， 模型自动退化至非线性浅水方程， 对于破 碎波演化的模拟获得了成功［１２⁃１３］。 该方法不需要额外的参数化模型， 只需在局部区域略去高阶色散项的求 解， 计算量没有丝毫增加， 效率较高。 非静压模型可借鉴这一技术， 于波浪破碎局部水域关闭动压项， 从而 退化至静压模型， 即模式的局部分裂法。 本文简要介绍了非静压模型的建立、 数值格式及波浪破碎的模式分裂法， 通过典型算例验证了模型的有 效性及模拟精