本报讯 近日，中国科学院南海海洋研究所热带海洋环境国家重点实验室（LTO）谢培炜博士构建了基于纳维-斯托克斯求解器的强非线性海啸波浪生成模型。该模型的构建弥补了海啸灾害领域中强非线性波浪生成模型的缺失，并在国际范围内首次实现对典型海啸波浪记录的模拟重现。相关成果由谢培炜博士为第一作者，杜岩研究员为共同作者，发表在国际期刊Coastal Engineering上。

图1.(a), (b)和(c)为海啸波浪上岸至岸滩中段时的波形及流速场图像，其中(a)无前置波谷；(b)前置波谷波幅较小；(c)前置波谷波幅较大。(d), (e)和(f)为波前抵达岸线位置且波浪破碎时的图像

2004年的印度洋海啸记录和2011年的日本海啸记录，都揭示了海啸波浪除具有长波的特征之外，还拥有明显的强非线性特征，即在近岸海域具有接近于水深的波高，通常伴随着一个前置波谷。在海啸灾害领域中，运用求解纳维-斯托克斯方程组的计算流体力学求解器生成海啸波浪是研究海啸的工程影响及岸滩冲刷效应的基础。在现有的基于纳维-斯托克斯求解器的海啸波浪生成模型中，传统的孤立波模型已被证实并不适用于海啸波浪的模拟研究，而采用线性浅水方程的模型亦无法生成具有强非线性特征的波浪。为满足海啸研究领域中模拟生成强非线性波浪的需求，该研究采用非线性浅水方程推导求得海啸波速边界条件的控制方程，并将该控制方程编译入计算纳维-斯托克斯方程组的开源求解器interFoam中，实现对海啸波浪生成模型的构建。通过与多个物理实验的比对验证得知，该模型具备了模拟海啸波浪生成、传播、陡化、频散和上岸破碎过程的能力。运用此模型实现了对2004年印度洋海啸和2011年日本海啸两组典型的强非线性海啸波浪的模拟重现，这在国际范围内为首例。

图2.数值模拟结果（红线）与 (a) 2004年印度洋海啸记录和 (b) 2011年日本海啸记录的对比（星形）

该文章还进一步研究了波谷前置型海啸波浪的爬升及上岸过程，展示了海表高度和流速场的空间和时间变化过程。研究表明，更大波幅的前置波谷会导致波浪上岸之前的岸线后退得更远，而后在上岸时产生更高的波浪爬高。随着前置波谷波幅的增加，上岸波浪的海向和陆向流速都得到加强。前置波谷的存在使得上岸波浪的波前更加陡峭，并可能引发波浪破碎，造成剧烈的岸滩侵蚀。

该研究由国家自然科学基金、中国博士后基金和LTO实验室自主项目等共同资助完成。

相关论文信息：Xie, P., & Du, Y. (2023). Tsunami wave generation in Navier–Stokes solver and the effect of leading trough on wave run-up. Coastal Engineering, 104293.

https://doi.org/10.1016/j.coastaleng.2023.104293(完)

通讯员|付恬 徐晓璐