01 研究背景

比利時的沿海地區(qū)有大量大型港口和風(fēng)電站，比利時的沿海地區(qū)對其經(jīng)濟(jì)發(fā)展至關(guān)重要。比利時海岸帶具有水位振幅大、水深小、懸浮沉積物濃度高的特點，因此非常需要可靠的數(shù)值模型來預(yù)測這一地區(qū)的水位和流速。

這個模型不僅要能夠?qū)Τr下由自北方進(jìn)入北海和從南方通過多佛海峽的潮汐波引起的水位和流速變化進(jìn)行準(zhǔn)確預(yù)測，還必須能夠在強(qiáng)風(fēng)極端天氣條件下給出準(zhǔn)確的預(yù)測結(jié)果。

02 案例展示

IMDC的工程師使用二維水動力通用仿真軟件建立了一個大陸架模型，對比利時海岸帶的水位進(jìn)行預(yù)測，利用該模型對2013年12月的Xaver氣旋的反演，證明該模型在極端氣候下也具有良好的預(yù)測能力。

03 模型搭建

模型研究的范圍囊括了西大西洋、北海、愛爾蘭海和波羅的海區(qū)域。網(wǎng)格共有49440個節(jié)點，海洋邊界處網(wǎng)格尺寸35 km，至近比利時海岸帶，瓦登群島和東丹麥等區(qū)域加密至500 m。

模型中的Scheldt河的河口部分建立了較密的河道網(wǎng)格，因此這個大陸架模型能夠很好地考慮到河流流動給海岸帶潮汐帶來的影響。由于河流上游來水流量的量級與河口潮汐流的流量相比非常小，因此忽略了河流上游來流的流量。

模型耦合OSU/TPXO潮汐數(shù)據(jù)庫對全場的水位進(jìn)行初始化，并在海洋邊界處施加水位邊界條件，工程師將科里奧利系數(shù)設(shè)為1.13 10-4定值，采用二維水動力通用仿真軟件以5min的時間步長進(jìn)行了仿真計算。使用myOcean.eu中水位數(shù)據(jù)對邊界潮汐水位和床面粗糙度進(jìn)行了校正后，將水位的校準(zhǔn)系數(shù)設(shè)置為1.0，曼寧粗糙度系數(shù)設(shè)置為0.0235 s/m1/3。

04 潮汐計算結(jié)果驗證

IMDC的工程師對該模型計算得到的一個月的數(shù)據(jù)和myOcean的測量數(shù)據(jù)進(jìn)行了對比分析，結(jié)果說明此模型很好地預(yù)測了潮汐的振幅譜分布。從這個分析可以看出，M2波最強(qiáng)，S2和N2具有相似的振幅，所有其他的潮汐波都較弱。

（藍(lán)色為二維水動力仿真模型結(jié)果數(shù)據(jù)；綠色為myOcean.eu測量數(shù)據(jù)；棕色為二者差異）

由于M2是潮汐波中最重要的組成成分，因此IMDC僅對M2成分的潮汐波的振幅和相位展開研究。為了驗證二維水動力仿真模型在不受氣象因素影響情況下的有效性，工程師利用此模型進(jìn)行了1年的模擬計算，得到了M2成分潮汐波的振幅與相位的時間過程，并與TOPEX的結(jié)果進(jìn)行了比較。

圖3展示了二維水動力通用仿真軟件計算得到的M2潮汐波振幅的空間分布，以及計算結(jié)果與TOPEX數(shù)據(jù)的差值。結(jié)果表明，在比利時海岸帶，法國海岸的部分地區(qū)（尤其在諾曼底）和英格蘭和威爾士的西海岸（特別是在賽弗斯河口）可以觀察到明顯的潮汐現(xiàn)象。

TOPEX中潮汐的振幅與二維水動力仿真計算結(jié)果的差異基本小于0.5m，而在淺水與靠近海岸的區(qū)域，振幅的差異往往更大。由于在這一部分區(qū)域，二維水動力模型的分辨率更高，計算結(jié)果會比TOPEX更加準(zhǔn)確。此外，模型的計算結(jié)果清晰表明，波羅的海附近的潮汐振幅非常小。事實上，波羅的海區(qū)域的水位變化主要是由氣象變化引起的。

05 Xaver氣旋的反演

Xaver氣旋是2013年12月發(fā)生的溫帶氣旋，氣旋的低壓和高風(fēng)速帶來的強(qiáng)暴風(fēng)雨在英格蘭、威爾士和比利時多個地區(qū)引起災(zāi)害，造成嚴(yán)重的損失。

為了對該氣旋進(jìn)行反演，IMDC的工程師使用來自GFS（全球預(yù)報系統(tǒng)）的空間分辨率為0.5度，時間分辨率為3h的氣壓與風(fēng)速數(shù)據(jù)，并線性插值到網(wǎng)格上，使用二維水動力仿真進(jìn)行了從2013年12月1日至2013年12月31日為期一個月的計算。

下圖展示了模型在Oostende和de Wandelaar兩個測點計算得到的水位結(jié)果與來自myOcean.eu的實測數(shù)據(jù)的對比。

（藍(lán)色為模型計算結(jié)果，橙色為站點實測水位）

從上圖可以看出，二維水動力仿真模型的計算結(jié)果是可靠的：模型很好地預(yù)測了12月5日的水位峰值，模擬結(jié)果與實際觀測水位僅相差0.1米；且在整個建模周期內(nèi)，計算水位的平均偏差誤差為0.02m，均方根誤差為0.37m。

06 研究結(jié)論

IMDC的工程師為了預(yù)測比利時海岸的水位和流速，建立了包含西大西洋、愛爾蘭海、波羅的海和北海的大陸架模型。通過模型計算結(jié)果與TOPEX數(shù)據(jù)的對比，對潮汐數(shù)據(jù)進(jìn)行了驗證：在北海和西大西洋的結(jié)果基本一致。

此外，工程師對2013年12月的熱帶氣旋Xaver進(jìn)行了反演分析，發(fā)現(xiàn)模型對比利時海岸帶的水位峰值預(yù)測非常準(zhǔn)確，且在氣旋期間，水位的偏差和均方根誤差僅分別為0.02米和0.27米。

07 小結(jié)

本文主要講述了IMDC的工程師利用二維水動力通用仿真軟件建立二維水動力模型，對比利時海岸帶的水位和流速進(jìn)行了模擬計算，并與TOPEX的實測數(shù)據(jù)與Xaver氣旋期間Oostende和de Wandelaar站點的實際測量結(jié)果進(jìn)行了對比。

IMDC的研究表明，使用二維水動力通用仿真軟件建立的大陸架模型，不僅可以很好地模擬常況下由潮汐波引起水位變化，更能夠很好地預(yù)測極端氣候條件下海岸帶的水位變化情況，具有相當(dāng)?shù)臏?zhǔn)確性和可靠性。

格物云CAE

一款國產(chǎn)可控云端仿真平臺，結(jié)構(gòu)、流體、水動力仿真軟件場景化模塊化，支持多格式網(wǎng)格導(dǎo)入(.med、.inp、.cdb、.cgns等)和高性能并行計算，降低CAE使用門檻，拓展CAE應(yīng)用范圍，加速工業(yè)企業(yè)研發(fā)制造數(shù)字化轉(zhuǎn)型。平臺支持云端CAE仿真生成工業(yè)APP，構(gòu)建完全交互式仿真社區(qū)，快速實現(xiàn)行業(yè)通用經(jīng)驗軟件化。

一鍵登錄，開啟仿真！

https://cae365.yuansuan.com

更多資訊可登錄格物CAE官方網(wǎng)站

https://cae.yuansuan.cn/

遠(yuǎn)算科技在bilibili、頭條、知乎、技術(shù)鄰定期發(fā)布課程視頻等內(nèi)容

或關(guān)注微信公眾號遠(yuǎn)算云仿真

敬請關(guān)注