近日，中國(guó)科學(xué)院南海海洋研究所熱帶海洋環(huán)境國(guó)家重點(diǎn)實(shí)驗(yàn)室（LTO）謝培煒博士構(gòu)建了基于納維-斯托克斯求解器的強(qiáng)非線性海嘯波浪生成模型。該模型的構(gòu)建彌補(bǔ)了海嘯災(zāi)害領(lǐng)域中強(qiáng)非線性波浪生成模型的缺失，并在國(guó)際范圍內(nèi)首次實(shí)現(xiàn)對(duì)典型海嘯波浪記錄的模擬重現(xiàn)。相關(guān)成果由謝培煒博士為第一作者，杜巖研究員為共同作者，發(fā)表在國(guó)際期刊Coastal Engineering上。

　　2004年的印度洋海嘯記錄和2011年的日本海嘯記錄，都揭示了海嘯波浪除具有長(zhǎng)波的特征之外，還擁有明顯的強(qiáng)非線性特征，即在近岸海域具有接近于水深的波高，通常伴隨著一個(gè)前置波谷。在海嘯災(zāi)害領(lǐng)域中，運(yùn)用求解納維-斯托克斯方程組的計(jì)算流體力學(xué)求解器生成海嘯波浪是研究海嘯的工程影響及岸灘沖刷效應(yīng)的基礎(chǔ)。在現(xiàn)有的基于納維-斯托克斯求解器的海嘯波浪生成模型中，傳統(tǒng)的孤立波模型已被證實(shí)并不適用于海嘯波浪的模擬研究，而采用線性淺水方程的模型亦無(wú)法生成具有強(qiáng)非線性特征的波浪。為滿足海嘯研究領(lǐng)域中模擬生成強(qiáng)非線性波浪的需求，該研究采用非線性淺水方程推導(dǎo)求得海嘯波速邊界條件的控制方程，并將該控制方程編譯入計(jì)算納維-斯托克斯方程組的開(kāi)源求解器interFoam中，實(shí)現(xiàn)對(duì)海嘯波浪生成模型的構(gòu)建。通過(guò)與多個(gè)物理實(shí)驗(yàn)的比對(duì)驗(yàn)證得知，該模型具備了模擬海嘯波浪生成、傳播、陡化、頻散和上岸破碎過(guò)程的能力。運(yùn)用此模型實(shí)現(xiàn)了對(duì)2004年印度洋海嘯和2011年日本海嘯兩組典型的強(qiáng)非線性海嘯波浪的模擬重現(xiàn)，這在國(guó)際范圍內(nèi)為首例。

　　該文章還進(jìn)一步研究了波谷前置型海嘯波浪的爬升及上岸過(guò)程，展示了海表高度和流速場(chǎng)的空間和時(shí)間變化過(guò)程。研究表明，更大波幅的前置波谷會(huì)導(dǎo)致波浪上岸之前的岸線后退得更遠(yuǎn)，而后在上岸時(shí)產(chǎn)生更高的波浪爬高。隨著前置波谷波幅的增加，上岸波浪的海向和陸向流速都得到加強(qiáng)。前置波谷的存在使得上岸波浪的波前更加陡峭，并可能引發(fā)波浪破碎，造成劇烈的岸灘侵蝕。

　　該研究由國(guó)家自然科學(xué)基金、中國(guó)博士后基金和LTO實(shí)驗(yàn)室自主項(xiàng)目等共同資助完成。

　　相關(guān)論文信息：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

　

　　圖1：(a), (b)和(c)為海嘯波浪上岸至岸灘中段時(shí)的波形及流速場(chǎng)圖像，其中(a)無(wú)前置波谷；(b)前置波谷波幅較小；(c)前置波谷波幅較大。(d), (e)和(f)為波前抵達(dá)岸線位置且波浪破碎時(shí)的圖像。

　　圖2：數(shù)值模擬結(jié)果（紅線）與 (a) 2004年印度洋海嘯記錄和 (b) 2011年日本海嘯記錄的對(duì)比（星形）。