馬玉祥 艾叢芳 董國(guó)海

非靜壓水波模型研究綜述*

馬玉祥 艾叢芳①董國(guó)海

(大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室 遼寧大連 116024)

水波數(shù)值模擬一直是水利、海岸、海洋工程和物理海洋等學(xué)科關(guān)注的重要研究領(lǐng)域之一。在諸多的水波模型之中, 非靜壓水波模型以其計(jì)算精度和效率較為平衡的特點(diǎn)受到青睞。非靜壓水波模型經(jīng)過(guò)近三十年的發(fā)展, 為科學(xué)研究、工程設(shè)計(jì)分析和海洋資源開發(fā)等提供了重要的技術(shù)手段。然而, 如何開發(fā)更為高效的非靜壓模型, 拓寬模型的應(yīng)用領(lǐng)域一直是研究的主要研究方向。文章綜述了非靜壓模型概念的提出以及數(shù)值開發(fā)的難點(diǎn), 介紹了非靜壓模型在波浪傳播演化和波浪與結(jié)構(gòu)物相互作用這兩方面的應(yīng)用進(jìn)展, 將有助于研究者了解非靜壓水波模型的研究和應(yīng)用現(xiàn)狀。

數(shù)值模擬; 水波模型; 非靜壓模型; 波浪與結(jié)構(gòu)物相互作用

諸多的水波模型都是基于一定的假設(shè)條件下由Navier-Stokes方程簡(jiǎn)化而來(lái), 它們都有各自的應(yīng)用范圍。例如, 對(duì)Navier-Stokes方程沿水深積分平均, 在一定的假設(shè)條件下, 可以推導(dǎo)出二維淺水方程。二維淺水方程適用于水平尺度遠(yuǎn)大于垂向尺度, 表面波波長(zhǎng)大于水深的水流問(wèn)題(包括河道洪水演進(jìn)、湖泊環(huán)流、海洋潮汐、風(fēng)暴潮等)的模擬。對(duì)于洪水演進(jìn)這樣的水流問(wèn)題, 二維淺水方程甚至一維淺水方程(即圣維南方程)的計(jì)算結(jié)果完全可以滿足工程上的需要。然而, 自然界中存在的水流運(yùn)動(dòng)都是三維的, 一、二維的模型無(wú)法得到水平速度沿水深的變化, 因此, 不能模擬水流和輸運(yùn)物質(zhì)在垂向的升降運(yùn)動(dòng)。早期由于計(jì)算機(jī)計(jì)算能力的制約, 廣泛應(yīng)用的三維模型是基于Navier-Stokes方程引入靜壓假設(shè)得到的三維淺水模型。由于三維淺水方程便于數(shù)值離散求解、計(jì)算量小等優(yōu)點(diǎn), 已被廣泛應(yīng)用于天然河道、湖泊、河口, 海岸及近岸水流的模擬中(Casulli, 1992; Li, 1998; Miglio, 1999;劉樺等, 2000; 賴錫軍等, 2006; Xing, 2012)。然而, 在水平尺度相對(duì)垂向尺度較小的表面波問(wèn)題以及潮汐與海底山脊相互作用引起的內(nèi)波問(wèn)題中, 靜壓假設(shè)已不合理, 這時(shí)需要采用考慮到非靜壓效應(yīng)的數(shù)學(xué)模型來(lái)求解這些問(wèn)題。

考慮到非靜壓效應(yīng)的數(shù)學(xué)模型有很多, 包括Boussinesq方程(Mohapatra, 2004; Mignot, 2009; Kim, 2011; Chang, 2014; Cantero- Chinchilla, 2016; 林鵬程等, 2021; 饒永紅等, 2021; 吳亞楠等, 2021)和勢(shì)流理論(Ma, 2001a, 2001b; Hu, 2002; Wang, 2010)等。Boussinesq方程考慮到了水平和垂直流速沿水深的變化, 通過(guò)方程中的頻散項(xiàng)來(lái)考慮非靜壓效應(yīng)的影響。相對(duì)于上面提到的淺水方程, Boussinesq方程被廣泛地應(yīng)用于近岸波浪傳播演化的模擬。在MIKE 21和SMS等商業(yè)軟件中均包含了基于Boussinesq方程開發(fā)的水波模型。在不可壓縮流動(dòng)的條件下, 假設(shè)流動(dòng)是無(wú)旋的, 速度勢(shì)便滿足拉普拉斯方程, 即可應(yīng)用勢(shì)流理論來(lái)求解涉及非靜壓效應(yīng)的水波問(wèn)題?；趧?shì)流理論的數(shù)值模型關(guān)鍵是求解速度勢(shì)滿足的拉普拉斯方程。速度勢(shì)包含了非靜壓效應(yīng)的影響, 求解拉普拉斯方程得到速度勢(shì)之后, 通過(guò)速度勢(shì)與流場(chǎng)的關(guān)系就可以直接計(jì)算出流場(chǎng)。與基于Euler方程或Navier-Stokes方程的數(shù)值模型相比, 基于勢(shì)流理論的數(shù)值模型求解的拉普拉斯方程與壓力Poisson方程類似, 在數(shù)學(xué)上都屬于橢圓型方程。然而, 基于勢(shì)流理論的數(shù)值模型無(wú)需求解動(dòng)量方程來(lái)計(jì)算流場(chǎng), 這是由流動(dòng)無(wú)旋假設(shè)導(dǎo)致的, 也是基于勢(shì)流理論的數(shù)值模型與基于Euler方程或Navier-Stokes方程的數(shù)值模型之間的最大區(qū)別。

從廣義上來(lái)講, 凡是已考慮到非靜壓效應(yīng)的數(shù)學(xué)模型, 并具備模擬頻散影響的數(shù)值模型皆可稱之為非靜壓模型。從狹義上來(lái)講, 以Euler方程或Navier-Stokes方程作為控制方程, 采用水位函數(shù)法計(jì)算自由表面的運(yùn)動(dòng)的數(shù)值模型才稱之為非靜壓模型(Lin, 2002; Stelling, 2003; Yuan, 2004; Bradford, 2005; Zijlema, 2005; Anthonio, 2006; Wu, 2010; Young, 2010; Ai, 2011, 2019a; Ma, 2012)。本文綜述的非靜壓水波模型即為狹義上來(lái)說(shuō)的非靜壓模型(Kirby, 2017)。狹義上所謂的“非靜壓模型”是相對(duì)淺水模型而來(lái)的。在淺水模型中, 只包括了靜水壓強(qiáng)的影響, 而非靜壓模型中,通常把壓力項(xiàng)分解為靜壓項(xiàng)和非靜壓項(xiàng)之和。相對(duì)于淺水模型, 非靜壓模型剛好多了非靜壓項(xiàng), 此即“非靜壓模型”的名稱由來(lái)。值得注意的是, 非靜壓模型的頻散效應(yīng)與非靜壓項(xiàng)直接相關(guān), 非靜壓項(xiàng)的求解是非靜壓模型最為耗時(shí)和關(guān)鍵的一步。

1 非靜壓模型在波浪運(yùn)動(dòng)模擬方面的進(jìn)展

非靜壓模型從20世紀(jì)末誕生之日起即被廣泛地應(yīng)用于模擬波浪運(yùn)動(dòng)。非靜壓模型的開發(fā)需要解決的一個(gè)關(guān)鍵問(wèn)題是如何準(zhǔn)確施加“零壓力自由表面邊界條件”。Casulli(1999)較早地將三維淺水模型(Casulli, 1992)擴(kuò)展, 建立了三維非靜壓自由表面流動(dòng)模型。然而, 此類非靜壓模型并不是“完全”非靜壓的模型, 它將自由表面的零壓力邊界條件直接施加在整個(gè)表層上, 即所謂的“表層靜壓假設(shè)”, 從而簡(jiǎn)化了基于標(biāo)準(zhǔn)交錯(cuò)網(wǎng)格非靜壓模型的離散求解, 構(gòu)造了對(duì)稱正定的壓力Poisson方程, 提高了非靜壓模型的求解效率。這類基于“表層靜壓假設(shè)”構(gòu)建的非靜壓模型(Casulli, 1999; Casulli, 2002)雖然具有較高的計(jì)算效率, 但模擬波浪運(yùn)動(dòng)需要增加垂向分層數(shù)來(lái)提高模型的計(jì)算精度。盡管如此, 在波浪傳播幾個(gè)周期之后仍會(huì)產(chǎn)生明顯的相位誤差, 而且波高也會(huì)有明顯的衰減。

Ai等(2011)和Ai等(2019a)開發(fā)的非靜壓表面流模型拋棄了傳統(tǒng)的交錯(cuò)定義變量的方式, 將垂向流速由分層網(wǎng)格界面移到分層網(wǎng)格中心定義, 同時(shí)水平流速的仍然按照交錯(cuò)的方式定義(圖1), 同樣構(gòu)建了模擬波浪運(yùn)動(dòng)的“完全”非靜壓模型。相對(duì)于其他模型而言(Stelling, 2003; Yuan, 2004; Zijlema, 2005; Anthonio, 2006; Cea, 2009; Wu, 2010; Young, 2010; Ma, 2012), 采用這種變量定義方式構(gòu)建的非靜壓模型的最大特點(diǎn)是最終求解的壓力Possion方程是對(duì)稱正定的, 可以采用預(yù)條件共軛梯度法高效求解, 這樣極大地提高了非靜壓模型的計(jì)算效率。Ai等(2011)基于垂向邊界擬合坐標(biāo)系統(tǒng)建立了模擬波浪運(yùn)動(dòng)的非靜壓模型(Ai, 2011)。該模型采用投影法, 即壓力修正法, 求解不可壓縮Euler方程。在水平笛卡爾網(wǎng)格框架和垂向邊界擬合坐標(biāo)系下, 提出了上述這種新穎的網(wǎng)格變量定義方式。通過(guò)與解析解和試驗(yàn)數(shù)據(jù)對(duì)比, 表明此非靜壓模型采用兩個(gè)垂向分層就能夠準(zhǔn)確高效地模擬波浪淺化、非線性、色散、折射和繞射現(xiàn)象。對(duì)于線性色散關(guān)系的模擬, 采用兩個(gè)垂向分層即能準(zhǔn)確模擬(誤差小于1%)=π的深水波浪(圖2)。為實(shí)現(xiàn)近岸波浪破碎爬高和強(qiáng)非線性波群的準(zhǔn)確模擬, Ai等(2012)對(duì)模型進(jìn)行了擴(kuò)展, 采用動(dòng)量守恒的計(jì)算格式離散動(dòng)量方程中的水平對(duì)流項(xiàng), 這使得模型能夠有效地模擬包括破碎波和水躍在內(nèi)的間斷流動(dòng)問(wèn)題。同時(shí), 再引入干濕動(dòng)邊界處理方法即實(shí)現(xiàn)了模型模擬波浪爬高的能力。數(shù)值模擬結(jié)果表明(圖3), 這樣構(gòu)建的非靜壓模型可以有效準(zhǔn)確地模擬近岸水深變化導(dǎo)致的波浪破碎和爬高過(guò)程。為降低非靜壓模型在模擬強(qiáng)非線性波群時(shí)由于壓力梯度項(xiàng)離散導(dǎo)致的數(shù)值誤差, 保證了非靜壓模型可以準(zhǔn)確地模擬強(qiáng)非線性波群的演化, Ai等(2014)首次在非靜壓模型中引入了廣義垂向邊界擬合坐標(biāo)系統(tǒng)。數(shù)值模擬結(jié)果表明(Ai, 2014), 基于廣義垂向邊界擬合坐標(biāo)系統(tǒng)建立的非靜壓模型可以準(zhǔn)確地模擬畸形波浪的產(chǎn)生(圖4)。最近, Ai等(2019a)提出了兩種半隱的非靜壓模型來(lái)模擬波浪運(yùn)動(dòng)。這兩種半隱模型同樣采用上述變量定義方式構(gòu)建, 一種是非迭代模型, 另一種需要迭代求解。這類半隱的非靜壓模型與之前的模型相對(duì)比, 具有時(shí)間步長(zhǎng)不受制于表面波波速的優(yōu)點(diǎn), 在模擬深水波浪問(wèn)題方面更具優(yōu)勢(shì)。研究評(píng)估和討論了這兩個(gè)半隱模型求解線性色散關(guān)系的精度和執(zhí)行效率。通過(guò)從淺水到深水的多個(gè)數(shù)值算例對(duì)這兩個(gè)模型進(jìn)行了驗(yàn)證。驗(yàn)證結(jié)果表明(Ai, 2019a), 這兩個(gè)半隱模型的結(jié)果非常相似, 均與實(shí)驗(yàn)數(shù)據(jù)吻合良好(圖5)。然而, 迭代模型的執(zhí)行效率低于非迭代模型, 求解迭代模型花費(fèi)的時(shí)間是非迭代模型的1.1～2.5倍。

圖1 非靜壓模型新穎的變量定義方式(Ai et al, 2011)

注:,,分別為,,方向的網(wǎng)格索引;,分別為水平流速和垂向流速;為非靜壓項(xiàng)

圖2 線性波波速的計(jì)算結(jié)果與解析解的對(duì)比 (Ai et al, 2011)

圖3 波浪爬高破碎過(guò)程(紅色實(shí)線: 模型結(jié)果; 圓點(diǎn): 實(shí)測(cè)值)(Ai et al, 2012)

圖4 畸形波的計(jì)算結(jié)果與實(shí)測(cè)數(shù)據(jù)的對(duì)比(紅色實(shí)線: 模型結(jié)果; 圓形: 實(shí)測(cè)值)(Ai et al, 2014)

圖5 聚焦點(diǎn)處三維深水聚焦波的計(jì)算結(jié)果與實(shí)測(cè)數(shù)據(jù)的對(duì)比(紅色虛線: 非迭代模型結(jié)果; 藍(lán)色實(shí)線: 迭代模型結(jié)果; 圓形: 實(shí)測(cè)值)(Ai et al, 2019a)

針對(duì)模擬強(qiáng)非線性波浪存在沿程波面衰減劇烈, 處理波浪破碎、強(qiáng)間斷等問(wèn)題上性能不足, 壓力分裂方法沒(méi)有及時(shí)更新變量信息以及簡(jiǎn)化的表面應(yīng)力條件不適用于強(qiáng)非線性波面等問(wèn)題, Dong等(2019)和He等(2020)提出了三階中心加權(quán)基本無(wú)振蕩格式(CWENO)和多級(jí)黎曼近似解MUSTA的組合方式提高重構(gòu)精度, 并構(gòu)建了動(dòng)壓預(yù)估校正兩步法, 同時(shí)完整考慮表面垂向和切向應(yīng)力, 實(shí)現(xiàn)了以較少的垂向分層對(duì)深水波浪傳播演化(Dong, 2019)、深水波浪破碎的準(zhǔn)確捕捉和模擬(He, 2020)。

波浪在孔隙介質(zhì)的傳播以及與孔隙結(jié)構(gòu)的相互作用涉及非常復(fù)雜的多尺度流動(dòng)耦合, 針對(duì)孔隙介質(zhì)大小形狀不一, 內(nèi)部湍流運(yùn)動(dòng)復(fù)雜, 傳統(tǒng)摩擦阻力系數(shù)難以率定以及實(shí)尺度、大范圍模擬計(jì)算成本高等問(wèn)題, He等(2022)提出了將體積平均的雷諾時(shí)均方程與非靜壓模型相結(jié)合的方式, 對(duì)體積平均后的-epsilon湍流模型進(jìn)行封閉, 統(tǒng)一了孔隙介質(zhì)內(nèi)外流動(dòng)的控制方程, 實(shí)現(xiàn)了波浪在包含孔隙結(jié)構(gòu)地形上傳播演化的高效準(zhǔn)確模擬。

2 非靜壓模型在波浪-結(jié)構(gòu)物相互作用模擬方面的進(jìn)展

在實(shí)際工程中, 常見的各類海洋結(jié)構(gòu)物尤其是目前發(fā)展較快的深水海洋結(jié)構(gòu)物當(dāng)中, 圓柱結(jié)構(gòu)十分常見, 如為海上船舶運(yùn)輸服務(wù)的各類透空式碼頭下部結(jié)構(gòu)中、利用海上風(fēng)能發(fā)電的各類固定式海上風(fēng)機(jī)下部結(jié)構(gòu)中都會(huì)見到圓柱結(jié)構(gòu)。波浪與直立圓柱的相互作用問(wèn)題一直是海岸和海洋工程界關(guān)注的問(wèn)題。在波浪與直立圓柱相互作用的過(guò)程中, 波浪的爬高和沖擊力可能會(huì)對(duì)圓柱造成嚴(yán)重的破壞, 因此諸如直立圓柱式海洋結(jié)構(gòu)物的設(shè)計(jì)必須考慮波浪的作用。目前, 研究者開發(fā)了許多數(shù)值模型來(lái)研究波浪與直立圓柱的相互作用(Li, 1999, 2001; Zhao, 2007; Zhong, 2009)。大多數(shù)非靜壓模型都是在矩形網(wǎng)格基礎(chǔ)上開發(fā)的。要實(shí)現(xiàn)波浪與直立圓柱或更為復(fù)雜外形結(jié)構(gòu)物相互作用的模擬, 應(yīng)構(gòu)建適應(yīng)于復(fù)雜邊界的非靜壓數(shù)值模型。Choi等(2011)基于水平邊界適應(yīng)的曲線坐標(biāo)系統(tǒng)建立了模擬波浪與復(fù)雜邊界相互作用的非靜壓模型。該模型采用基于高階積分方法的積分表層動(dòng)量方程以提高表層非靜壓項(xiàng)的處理精度。數(shù)值模擬結(jié)果表明, 模型可以較為準(zhǔn)確地模擬橢圓余弦波在直立圓柱上的爬高。此后, Choi等(2012)在上述模型的基礎(chǔ)上, 又引入了二階Adams- Bashforth格式和三階迎風(fēng)格式來(lái)離散對(duì)流項(xiàng)以提高模型的精度。在模擬非線性波浪在直立圓柱上的爬高問(wèn)題時(shí), 模型計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)和二階繞射理論結(jié)果吻合較好。除了上述兩個(gè)基于水平邊界適應(yīng)的曲線坐標(biāo)系統(tǒng)建立的非靜壓模型具備模擬波浪與復(fù)雜邊界相互作用的能力之外, 其他基于非結(jié)構(gòu)三角形網(wǎng)格構(gòu)建的非靜壓模型同樣具備這樣的能力(Casulli, 2002; Cea, 2009)。

圖6 三角形與四邊形混合網(wǎng)格(Ai et al, 2017)

大連理工大學(xué)在國(guó)內(nèi)較早地建立了模擬波浪與復(fù)雜邊界結(jié)構(gòu)物相互作用的非結(jié)構(gòu)網(wǎng)格非靜壓模型(Ai, 2010, 2017)。三維計(jì)算網(wǎng)格采用覆蓋水平計(jì)算域的二維計(jì)算網(wǎng)格和垂向分層式的結(jié)構(gòu)化網(wǎng)格構(gòu)建。水平二維計(jì)算網(wǎng)格可以采用三角形網(wǎng)格(Ai, 2010)或三角形與四邊形相結(jié)合的混合網(wǎng)格(Ai, 2017)。水平對(duì)流項(xiàng)的離散采用Perot格式實(shí)現(xiàn)(Perot, 2000)。該非靜壓模型可以模擬波浪與單直立圓柱或多直立圓柱的相互作用。通過(guò)與試驗(yàn)數(shù)據(jù)對(duì)比, 證實(shí)了模型可以準(zhǔn)確地模擬非線性波浪在直立圓柱周圍的爬高現(xiàn)象以及作用在直立圓柱上的波浪力。

除了非靜壓模型在波浪與直立圓柱相互作用方面的模擬應(yīng)用, 非靜壓模型也被應(yīng)用于波浪與水下結(jié)構(gòu)物或浮式結(jié)構(gòu)物相互作用的模擬。Lin(2006)開發(fā)了三維多層坐標(biāo)非靜壓模型, 以模擬波浪與結(jié)構(gòu)物的相互作用, 該模型克服了之前非靜壓模型(Lin, 2002)無(wú)法準(zhǔn)確模擬涉及間斷水深的問(wèn)題, 僅采用三個(gè)垂向分層, 該模型即可準(zhǔn)確地模擬波浪與浮式或水下結(jié)構(gòu)物的相互作用。Kang等(2015)對(duì)該多層坐標(biāo)模型進(jìn)行了擴(kuò)展, 以模擬涉及曲面問(wèn)題的波浪與結(jié)構(gòu)物相互作用問(wèn)題, 并成功地將模型應(yīng)用于波浪與不同淹沒(méi)深度的垂直圓柱相互作用的模擬研究。Ma等(2016)同樣提出了模擬波浪與結(jié)構(gòu)物相互作用的坐標(biāo)非靜壓模型。對(duì)于波浪與底部固定結(jié)構(gòu)物的相互作用, 模型結(jié)果與實(shí)驗(yàn)數(shù)據(jù)之間取得了良好的一致性, 但模型在解決波浪與浮式或水下結(jié)構(gòu)物相互作用方面的能力尚未得到證實(shí)。Rijnsdorp等(2016)提出了一個(gè)基于SWASH (Zijlema, 2011)的非靜壓數(shù)值模型。他們的模型能夠預(yù)測(cè)波浪的非線性演變及其對(duì)受約束船舶的影響。然而, 該模型不能直接用于模擬波浪與浮式結(jié)構(gòu)物的相互作用。

近年來(lái), Ai等(2018), Ai等(2019b), Ma等(2019)和Ai等(2022b)在其開發(fā)的非靜壓波浪數(shù)值模型的基礎(chǔ)上, 通過(guò)引入浸入邊界法實(shí)現(xiàn)了波浪與浮式或水下結(jié)構(gòu)物相互作用的模擬。浸入邊界法的引入使得非靜壓模型實(shí)現(xiàn)了基于結(jié)構(gòu)化網(wǎng)格模擬波浪與任意形狀結(jié)構(gòu)物相互作用的能力。基于廣義垂向邊界擬合坐標(biāo)系統(tǒng), 分別采用半隱分部算法求解Navier-Stokes方程(Ai, 2018)和顯式投影法求解Euler方程(Ma, 2019), 實(shí)現(xiàn)了波浪與浮式或水下任意形狀結(jié)構(gòu)物相互作用的立面二維模擬。數(shù)值模擬結(jié)果表明, 所建立的非靜壓模型可以準(zhǔn)確地模擬出浮式或水下結(jié)構(gòu)物對(duì)波浪傳播演化的影響(圖7)以及波浪作用下結(jié)構(gòu)物的受力情況(Ai, 2018; Ma, 2019)。另外, 模型還可以模擬波浪與結(jié)構(gòu)物相互作用過(guò)程中旋渦的產(chǎn)生以及結(jié)構(gòu)物附近旋渦的演化過(guò)程(圖8)。由于上述兩個(gè)非靜壓模型均為立面二維模型, 為實(shí)現(xiàn)波浪與三維結(jié)構(gòu)物相互作用的模擬, 對(duì)上述模型進(jìn)行擴(kuò)展建立了三維非靜壓模型(Ai, 2019b)。這樣建立的三維模型, 由于浸入邊界法的引入破壞了壓力Poisson方程的對(duì)稱性, 降低了模型的求解效率。為改進(jìn)模型的執(zhí)行效率, Ai等(2022b)首次將浸入邊界法與全局連續(xù)性方程相結(jié)合建立了高效的三維非靜壓模型, 在該模型中全局連續(xù)性方程用來(lái)替代固液邊界處的Newmann邊界條件, 避免了浸入邊界法對(duì)壓力Poisson方程的不利影響, 保證了壓力Poisson方程對(duì)稱正定的性質(zhì), 有效提高了模型的計(jì)算效率。在三維波浪與浮式結(jié)構(gòu)物相互作用的模擬中(圖9), 計(jì)算效率最高可提高7.6倍(Ai, 2022b)。

3 結(jié)論

非靜壓水波模型的應(yīng)用十分廣泛, 凡是涉及非靜壓效應(yīng)的自由表面流動(dòng)現(xiàn)象都可采用非靜壓模型模擬。本文綜述了非靜壓模型概念的提出以及數(shù)值開發(fā)的難點(diǎn), 主要介紹了非靜壓模型在波浪傳播演化和波浪與結(jié)構(gòu)物相互作用這兩方面的應(yīng)用進(jìn)展。本文介紹的非靜壓模型是基于不可壓縮Euler方程或Navier-Stokes方程建立的, 在這類模型中, 壓力Poisson方程的求解是不可避免的。壓力Poisson方程的性質(zhì)直接關(guān)系整個(gè)模型的執(zhí)行效率。高效的非靜壓水波模型一直是研究者不懈追求的目標(biāo)。作者希望本文能引起更多的研究者加入到非靜壓水波模型的開發(fā)應(yīng)用中, 拓寬非靜壓模型的應(yīng)用領(lǐng)域。除了本文綜述的內(nèi)容之外, 非靜壓模型在內(nèi)波產(chǎn)生和傳播演化(Lai, 2010; Ai, 2016)、海嘯波傳播(Oishi, 2013)、滑坡涌浪(Ai, 2021b)、波動(dòng)涌潮(Ai, 2021a)、潰壩波傳播演化(Ai, 2022a)等方面均有應(yīng)用。另外, 為提高非靜壓模型的執(zhí)行效率, 也可以采用GPU技術(shù)并行求解模型(Ai, 2019a)。

圖7 受水下結(jié)構(gòu)物影響的波面時(shí)間序列的計(jì)算結(jié)果與實(shí)測(cè)數(shù)據(jù)的對(duì)比(藍(lán)色實(shí)線: 模型結(jié)果; 圓形: 實(shí)測(cè)值)(Ma et al, 2019)

圖8 結(jié)構(gòu)物附近的流場(chǎng)

圖9 波浪與浮式沉箱的相互作用

劉樺, 何友聲, 2000. 河口三維流動(dòng)數(shù)學(xué)模型研究進(jìn)展[J]. 海洋工程, 18(2): 87-93.

吳亞楠, 武賀, 周慶偉, 等, 2021. 基于Boussinesq波浪模型的港池波浪數(shù)值模擬與泊穩(wěn)分析[J]. 海洋通報(bào), 40(3): 301-308.

林鵬程, 劉忠波, 劉勇, 2021. 基于Boussinesq數(shù)值模型的波浪速度垂向分布模擬研究[J]. 海洋湖沼通報(bào), 43(4): 7-15.

饒永紅, 劉忠波, 梁書秀, 等, 2021. 雙層Boussinesq模型非線性波浪模擬研究[J]. 水道港口, 42(5): 614-622.

賴錫軍, 曲卓杰, 周杰, 等, 2006. 非結(jié)構(gòu)網(wǎng)格上的三維淺水流動(dòng)數(shù)值模型[J]. 水科學(xué)進(jìn)展, 17(5): 693-699.

AI C F, DING W Y, 2016. A 3D unstructured non-hydrostatic ocean model for internal waves [J]. Ocean Dynamics, 66(10): 1253-1270.

AI C F, DING W Y, JIN S, 2014. A general boundary-fitted 3D non-hydrostatic model for nonlinear focusing wave groups [J]. Ocean Engineering, 89: 134-145.

AI C F, DING W Y, JIN S, 2017. A hybrid-grid 3D model for regular waves interacting with cylinders [J]. Journal of Hydraulic Research, 55(1): 129-134.

AI C F, JIN S, 2010. Non-hydrostatic finite volume model for non-linear waves interacting with structures [J]. Computers & Fluids, 39(10): 2090-2100.

AI C F, JIN S, 2012. A multi-layer non-hydrostatic model for wave breaking and run-up [J]. Coastal Engineering, 62: 1-8.

AI C F, JIN S, LV B, 2011. A new fully non-hydrostatic 3D free surface flow model for water wave motions [J]. International Journal for Numerical Methods in Fluids, 66(11): 1354-1370.

AI C F, MA Y X, DING W Y,, 2021a. An efficient three-dimensional non-hydrostatic model for undular bores in open channels [J]. Physics of Fluids, 33(12): 127111.

AI C F, MA Y X, DING W Y,, 2022a. Three-dimensional non-hydrostatic model for dam-break flows [J]. Physics of Fluids, 34(2): 022105.

AI C F, MA Y X, YUAN C F,, 2018. Semi-implicit non-hydrostatic model for 2D nonlinear wave interaction with a floating/suspended structure [J]. European Journal of Mechanics-B/Fluids, 72: 545-560.

AI C F, MA Y X, YUAN C F,, 2019a. Development and assessment of semi-implicit nonhydrostatic models for surface water waves [J]. Ocean Modelling, 144: 101489.

AI C F, MA Y X, YUAN C F,, 2019b. A 3D non-hydrostatic model for wave interactions with structures using immersed boundary method [J]. Computers & Fluids, 186: 24-37.

AI C F, MA Y X, YUAN C F,, 2021b. A three-dimensional non-hydrostatic model for tsunami waves generated by submarine landslides [J]. Applied Mathematical Modelling, 96: 1-19.

AI C F, MA Y X, YUAN C F,, 2022b. An efficient 3D non-hydrostatic model for predicting nonlinear wave interactions with fixed floating structures [J]. Ocean Engineering, 248: 110810.

ANTHONIO S L, HALL K R, 2006. High-order compact numerical schemes for non-hydrostatic free surface flows [J]. International Journal for Numerical Methods in Fluids, 52(12): 1315-1337.

BRADFORD S F, 2005. Godunov-based model for nonhydrostatic wave dynamics [J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 131(5): 226-238.

CANTERO-CHINCHILLA F N, CASTRO-ORGAZ O, DEY S,, 2016. Nonhydrostatic dam break flows. I: physical equations and numerical schemes [J]. Journal of Hydraulic Engineering, 142(12): 04016068.

CASULLI V, 1999. A semi-implicit finite difference method for non-hydrostatic, free-surface flows [J]. International Journal for Numerical Methods in Fluids, 30(4): 425-440.

CASULLI V, CHENG R T, 1992. Semi-implicit finite difference methods for three-dimensional shallow water flow [J]. International Journal for Numerical Methods in Fluids, 15(6): 629-648.

CASULLI V, ZANOLLI P, 2002. Semi-implicit numerical modeling of nonhydrostatic free-surface flows for environmental problems [J]. Mathematical and Computer Modelling, 36(9/10): 1131-1149.

CEA L, STELLING G, ZIJLEMA M, 2009. Non-hydrostatic 3D free surface layer-structured finite volume model for short wave propagation [J]. International Journal for Numerical Methods in Fluids, 61(4): 382-410.

CHANG T J, CHANG K H, KAO H M, 2014. A new approach to model weakly nonhydrostatic shallow water flows in open channels with smoothed particle hydrodynamics [J]. Journal of Hydrology, 519: 1010-1019.

CHOI D Y, WU C H, YOUNG C C, 2011. An efficient curvilinear non-hydrostatic model for simulating surface water waves [J]. International Journal for Numerical Methods in Fluids, 66(9): 1093-1115.

CHOI D Y, YUAN H L, 2012. A horizontally curvilinear non-hydrostatic model for simulating nonlinear wave motion in curved boundaries [J]. International Journal for Numerical Methods in Fluids, 69(12): 1923-1938.

DONG G H, FU R L, MA Y X,, 2019. Simulation of unidirectional propagating wave trains in deep water using a fully non-hydrostatic model [J]. Ocean Engineering, 180: 254-266.

HE D B, MA Y X, DONG G H,, 2020. Predicting deep water wave breaking with a non-hydrostatic shock-capturing model [J]. Ocean Engineering, 216: 108041.

HE D B, MA Y X, DONG G H,, 2022. A numerical investigation of wave and current fields along bathymetry with porous media [J]. Ocean Engineering, 244: 110333.

HU P X, WU G X, MA Q W, 2002. Numerical simulation of nonlinear wave radiation by a moving vertical cylinder [J]. Ocean Engineering, 29(14): 1733-1750.

KANG A Z, LIN P Z, LEE Y J,, 2015. Numerical simulation of wave interaction with vertical circular cylinders of different submergences using immersed boundary method [J]. Computers & Fluids, 106: 41-53.

KIM D H, LYNETT P J, 2011. Dispersive and nonhydrostatic pressure effects at the front of surge [J]. Journal of Hydraulic Engineering, 137(7): 754-765.

KIRBY J T, 2017. Recent advances in nearshore wave, circulation, and sediment transport modeling [J]. Journal of Marine Research, 75(3): 263-300.

LAI Z G, CHEN C S, COWLES G W,, 2010. A nonhydrostatic version of FVCOM: 2. Mechanistic study of tidally generated nonlinear internal waves in Massachusetts Bay [J]. Journal of Geophysical Research, 115(C12): C12049.

LI Y S, LIU S X, YU Y X,, 1999. Numerical modeling of Boussinesq equations by finite element method [J]. Coastal Engineering, 37(2): 97-122.

LI Y S, ZHAN J M, 1998. Three-dimensional finite-element model for stratified coastal seas [J]. Journal of Hydraulic Engineering, 124(7): 699-703.

LI Y S, ZHAN J M, 2001. Boussinesq-type model with boundary-fitted coordinate system [J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 127(3): 152-160.

LIN P Z, 2006. A multiple-layer-coordinate model for simulation of wave-structure interaction [J]. Computers & Fluids, 35(2): 147-167.

LIN P Z, LI C W, 2002. A-coordinate three-dimensional numerical model for surface wave propagation [J]. International Journal for Numerical Methods in Fluids, 38(11): 1045-1068.

MA G F, FARAHANI A A, KIRBY J T,, 2016. Modeling wave-structure interactions by an immersed boundary method in a-coordinate model [J]. Ocean Engineering, 125: 238-247.

MA G F, SHI F Y, KIRBY J T, 2012. Shock-capturing non-hydrostatic model for fully dispersive surface wave processes [J]. Ocean Modelling, 43/44: 22-35.

MA Q W, WU G X, TAYLOR R E, 2001a. Finite element simulation of fully non-linear interaction between vertical cylinders and steep waves. Part 1: methodology and numerical procedure [J]. International Journal for Numerical Methods in Fluids, 36(3): 265-285.

MA Q W, WU G X, TAYLOR R E, 2001b. Finite element simulations of fully non-linear interaction between vertical cylinders and steep waves. Part 2: numerical results and validation [J]. International Journal for Numerical Methods in Fluids, 36(3): 287-308.

MA Y X, YUAN C F, AI C F,, 2019. Comparison between a non-hydrostatic model and OpenFOAM for 2D wave-structure interactions [J]. Ocean Engineering, 183: 419-425.

MIGLIO E, QUARTERONI A, SALERI F, 1999. Finite element approximation of quasi-3D shallow water equations [J]. Computer Methods in Applied Mechanics and Engineering, 174(3/4): 355-369.

MIGNOT E, CIENFUEGOS R, 2009. On the application of a Boussinesq model to river flows including shocks [J]. Coastal Engineering, 56(1): 23-31.

MOHAPATRA P K, CHAUDHRY M H, 2004. Numerical solution of Boussinesq equations to simulate dam-break flows [J]. Journal of Hydraulic Engineering, 130(2): 156-159.

OISHI Y, PIGGOTT M D, MAEDA T,, 2013. Three-dimensional tsunami propagation simulations using an unstructured mesh finite element model [J]. Journal of Geophysical Research: Solid Earth, 118(6): 2998-3018.

PEROT B, 2000. Conservation properties of unstructured staggered mesh schemes [J]. Journal of Computational Physics, 159(1): 58-89.

RIJNSDORP D P, ZIJLEMA M, 2016. Simulating waves and their interactions with a restrained ship using a non-hydrostatic wave-flow model [J]. Coastal Engineering, 114: 119-136.

STELLING G, ZIJLEMA M, 2003. An accurate and efficient finite-difference algorithm for non-hydrostatic free-surface flow with application to wave propagation [J]. International Journal for Numerical Methods in Fluids, 43(1): 1-23.

WANG C Z, WU G X, 2010. Interactions between fully nonlinear water waves and cylinder arrays in a wave tank [J]. Ocean Engineering, 37(4): 400-417.

WU C H, YOUNG C C, CHEN Q,, 2010. Efficient nonhydrostatic modeling of surface waves from deep to shallow water [J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 136(2): 104-118.

XING Y, AI C F, JIN S, 2013. A three-dimensional hydrodynamic and salinity transport model of estuarine circulation with an application to a macrotidal estuary [J]. Applied Ocean Research, 39: 53-71.

YOUNG C C, WU C H, 2010. Nonhydrostatic modeling of nonlinear deep-water wave groups [J]. Journal of Engineering Mechanics, 136(2): 155-167.

YUAN H L, WU C H, 2004. An implicit three-dimensional fully non-hydrostatic model for free-surface flows [J]. International Journal for Numerical Methods in Fluids, 46(7): 709-733.

ZHAO M, CHENG L, TENG B, 2007. Numerical simulation of solitary wave scattering by a circular cylinder array [J]. Ocean Engineering, 34(3/4): 489-499.

ZHONG Z Y, WANG K H, 2009. Modeling fully nonlinear shallow-water waves and their interactions with cylindrical structures [J]. Computers & Fluids, 38(5): 1018-1025.

ZIJLEMA M, STELLING G S, 2005. Further experiences with computing non-hydrostatic free-surface flows involving water waves [J]. International Journal for Numerical Methods in Fluids, 48(2): 169-197.

ZIJLEMA M, STELLING G, SMIT P, 2011. SWASH: an operational public domain code for simulating wave fields and rapidly varied flows in coastal waters [J]. Coastal Engineering, 58(10): 992-1012.

A REVIEW ON NON-HYDROSTATIC WATER WAVE MODELS

MA Yu-Xiang, AI Cong-Fang, DONG Guo-Hai

(Dalian University of Technology, State Key Laboratory of Coastal and Offshore Engineering, Dalian 116024, China)

Water wave numerical simulation has always been one of the important research fields of hydraulic engineering, coastal engineering, marine engineering, and physical oceanography. Among many water wave models, non-hydrostatic water wave models are favored by researchers because of its balance between computational accuracy and efficiency. After nearly three decades of development, non-hydrostatic water wave models provide important technical means for scientific research, engineering design and analysis, and marine resources development. However, developing more-efficient non-hydrostatic models and broadening the application field of the models have always been the pursuit of non-hydrostatic model researchers. In this paper, we first introduce the concept of “non-hydrostatic pressure”, and then review the application of non-hydrostatic models in wave propagations and evolutions, and wave-structure interactions from the perspective of model development and application.

numerical simulation; water wave model; non-hydrostatic model; wave-structure interactions

* 國(guó)家自然科學(xué)基金項(xiàng)目, 51979029號(hào); 國(guó)家自然科學(xué)基金項(xiàng)目, 52171248號(hào); 國(guó)家自然科學(xué)基金項(xiàng)目, 51720105010號(hào); 中央高?；究蒲袠I(yè)務(wù)費(fèi), DUT21LK01號(hào)。馬玉祥, 教授, E-mail: yuxma@dlut.edu.cn

艾叢芳, 碩士生導(dǎo)師, 副教授, E-mail: aicongfang@dlut.edu.cn

2022-02-28,

2022-03-23

P753

10.11693/hyhz20220200041