【CAE案例】基于二維水動力仿真的大陸架建模

泡沫oO

2023年9月21日 11:14

瀏覽：2621 收藏：1

01 研究背景

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

這個模型不僅要能夠對常況下由自北方進入北海和從南方通過多佛海峽的潮汐波引起的水位和流速變化進行準確預測，還必須能夠在強風極端天氣條件下給出準確的預測結果。

02 案例展示

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

03 模型搭建

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

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

圖1. 模型的網格和底部高程及沿岸細化展示

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

04 潮汐計算結果驗證

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

圖2. Wandelaar的潮汐波振幅譜對比

（藍色為二維水動力仿真模型結果數據；

綠色為myOcean.eu測量數據；

棕色為二者差異）

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

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

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