本站小編為你精心準備了二維箱體結構反射和透射的模擬研究參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《哈爾濱工程大學學報》2016年第三期

摘要：

基于兩層流體簡化模型，研究了水面箱型結構在兩層流體中對波浪的反射和透射特性。為研究復雜海洋結構，采用外域級數展開和內域邊界元數值方法相結合的手段建立了頻域數學模型，在外域和內域交界的輻射面上利用連續性條件建立聯立方程組。其中內域計算采用Rankine源作為格林函數，既避免了求解復雜的格林函數，又可直接獲得自由水面及內界面上的速度勢。經過與級數展開的解析方法對比，以及兩層流中波能流守恒的驗證，證明了本文建立的這一聯合數值模型的正確性。利用這一數值模型通過計算反射系數和透射系數的變化情況，研究了兩層流中存在箱型結構時各模態波浪的轉化情況，分析了兩層流中箱型結構的反射和透射特性。

關鍵詞：

兩層流體；邊界元方法；Rankine源；反射系數；透射系數

波浪經過結構物后，透射波浪的能量和波面變化情況是設計關心的主要問題之一。通常的研究都假定海水密度均一，但實際海洋環境中，海水密度在垂向上存在顯著的層化現象。這一特性導致在海洋內部也會產生波動，因此對結構的反射、透射問-2-哈爾濱工程大學學報第卷維情況下作用于物體上的激振力以及運動響應特性，并與物理模型試驗結果進行了對比。NguyenandYeung[7]推導出三維情況下有限水深的兩層流體中滿足自由水面條件和內界面條件的格林函數。計算結果表明，在一定頻率范圍內分層效應對結構物水動力特性有顯著影響，并且當入射波浪以內波模態入射時會在特定頻率范圍內引起方箱較大的轉動。勾瑩等[8]采用簡單格林函數建立了時域方法求解積分方程，對漂浮圓柱、方箱等的水動力特性進行了研究。石強[9,10]基于分離變量法，建立了大直徑圓柱體和方箱在兩層流體中繞射勢與輻射勢的解析計算方法。Zhangetal[11]運用特征函數匹配的方法，研究了兩層流中月池結構垂向振蕩的附加質量和輻射阻尼，以及在共振頻率時自由水面和內界面上的波動情況。本文采用數值方法和解析方法相結合的研究手段，將流域分為內域和外域兩部分，內域采用利用簡單格林函數的邊界元法求解，外域速度勢及導數利用解析展開式表示。這一模型既能求解任意物體的水動力特性，又避免了求解復雜的滿足自由水面條件的格林函數；同時外域采用解析方法也使無窮遠處的散射條件得到了很好的滿足。通過計算方箱的激振力并與解析方法結果對比，驗證了模型的正確性。根據單層流中波能流守恒原理推導出了兩層流中的波能流守恒關系，利用這一關系對模型的精確性進行了進一步的驗證。最后給出了兩層流體中各模態的反射系數及透射系數的定義，研究了某一模態入射波浪與物體作用后，自由水面及內界面上的波面變化情況。

1數學模型

考慮二維情況下的兩層理想流體，定義直角坐標系OXZ，坐標軸OZ垂直向上，OX與靜水面重合。

1.1入射速度勢兩層流模型中存在兩個模態的入射波浪，設射波的方向與x軸正向相同，波數為()0mk的入射波，速度勢可表示。

1.2反射和透射速度勢當入射波作用到結構物上后，會產生反射波和透射波，其速度勢用d表示，并且滿足以下定解條件。在自由水面、物面和內界面積分時，未知量為速度勢，而在輻射面積分時，通過輻射面上內、外域速度勢及其導數的連續性，將未知量轉化為外域速度勢展開式中的系數A和B。最后通過在自由水面、物面、內界面及輻射面劃分網格，對積分方程(11)、(12)分別進行離散可得兩組線性方程組。

2模型的驗證

為驗證所建立模型及求解過程的正確性，以水面固定箱型結構為例，計算了其波浪激振力以及波能流的守恒關系。

2.1網格劃分及收斂性驗證不同節點數的計算結果如圖2所示。圖中縱坐標表示表面波模態波浪從左側入射時，自由水面和內界面固定點x=8.5m處波面升高隨頻率的變化情況。經比較可知，節點數為1062時數值結果已經收斂。此時自由水面節點數為202，物面節點數為239，內界面節點數為343，上下輻射面分別為178和100。波浪激振力圖2至圖4分別給出了方箱在表面波模態和內波模態波浪入射時的無因次化激振力，其中水平波浪力和垂向波浪力均以()10mgah無因次化，繞y軸力矩以()10mgaBh無因次化。由數值計算結果與解析解[10]的對比可以看出，本文的計算結果與解析解吻合很好，驗證了本文求解模型的正確性。由圖可知，當入射波浪為表面波模態時，水平波浪力、垂向波浪力和繞y軸力矩都與單層流體中作用力基本重合，說明此時激振力基本不受分層效應的影響。當入射波浪為內波模態時，三個方向的作用力在低頻處都有多個極值和零點，并且遠小于表面波模態入射的激振力，隨著頻率的不斷增大，激振力最后都減小為0。

2.2能量守恒關系根據波能流守恒定理可得。圖5兩曲線分別為入射波為表面波模態和內波模態時，模型所計算出的能量P隨無因次化頻率的變化關系，可以看出所建立模型能很好地滿足能量守恒關系式(16)，進一步證明了模型的精確性。

3反射、透射系數的定義與計算結果

雖然自由水面和內界面上都存在兩種模態的波浪，但是由于自由水面上以表面波模態波浪成分為主，內界面上以內波模態波浪成分為主，因此在后續分析中，只給出自由水面上表面波模態波浪的反透射系數和內界面上內波模態波浪的反透射系數。其余系數可通過關系式(2)確定。圖6(a)給出了表面波模態波浪入射時，自由水面上表面波模態的反射和透射系數，其中L1為表面波模態波長。可以看出，隨波長的減小自由水面上表面波模態反射系數趨近于1，透射系數趨近于0，趨勢與單層流體一致。圖6(b)給出了表面波模態波浪入射時，內界面上內波模態波浪的反射和透射系數，其中L2為內波模態波長。由圖可知，表面波模態的入射波浪與結構物作用后將在內界面上產生較大幅度的波動，最高可接近自由表面入射波浪波幅的2倍左右。這是因為兩層流體間密度差很小，其恢復力很小，相當于將分層介質置于微重力場，流體受到很小的擾動就會偏離其平衡位置而產生“軒然大波”，但其對于能量的貢獻較小。從圖中還可以看出隨著入射波浪頻率的變化，反射、透射系數有多個零點，即在某些頻率處，內界面上沒有內波模態的反射波或透射波。此外，自由水面上表面波模態的反射系數和透射系數變化趨勢（圖6(a)），當結構位于上層流體時，內界面上內波模態的反射系數和透射系數的變化趨勢則相同。圖7(a)給出了內波模態波浪入射時，自由水面上表面波模態的反射和透射系數。由圖可知，當內波模態波浪入射時，入射波與結構物作用后，自由水面僅有非常微小的波動，所引起的波面升高基本可以忽略。圖7(b)給出了內波模態波浪入射時，內界面上內波模態波浪的反射和透射系數。由圖可知，由于本文所研究的結構物只位于上層流體，因此內界面上內波模態透射系數很大，接近于1，反射系數則很小，高頻時趨于0。根據圖6和圖7的結果可知，表面波模態入射時自由水面上表面波模態反透射特性與單層流體趨于一致；內波模態入射時自由水面上表面波模態反透射系數都很小接近于0，內界面上內波模態透射系數很大接近于1，反射系數則很小。因此下文僅研究表面波模態入射時，密度比、水深比及方箱寬度等參數的變化對內界面上內波模態反透射系數的影響規律。圖8為改變上下層密度比、入射波浪為表面波模態時，內界面上內波模態的反射和透射系數。由圖可以看出，隨著下層流體密度的增大，內界面反射和透射系數的最大值都減小，這是因為隨著密度差的增大，恢復力相應增大，內界面處抗擾動能力增強。同時，從結果中可以看出反射和透射系數零點的位置并不隨密度比的改變而改變，只與B/L2有關。反射系數零點出現在B/L2約為0.42、0.90和1.35的位置，透射系數零點出現在B/L2約為0.37、0.84和1.30的位置。圖9為改變上下層水深比、入射波浪為表面波模態時，內界面內波模態的反透射系數。由圖可知，結構物離內界面越近，則波浪入射后內界面的擾動越大，即反射和透射系數越大。內界面的位置同樣不影響反射和透射系數零點的位置。圖10為改變方箱寬度、入射波為表面波模態時，內界面內波模態的反射和透射系數。由圖可知隨結構物寬度的增加，內界面內波模態的反射和透射系數在低頻時變化不大，高頻時明顯增大。但反射系數和透射系數的零點位置都幾乎不發生變化。

4結論

本文基于海洋密度層化的兩層流體簡化模型，研究了兩層流體中箱型結構對波浪的反射和透射特性。在線性勢流理論下，采用外域級數展開和內域邊界元數值方法相結合的手段建立了頻域數學模型，在外域和內域交界的輻射面上利用連續性條件建立聯立方程組。在上下兩層流體區域分別選取Rankine源及Rankine源和它關于海底的鏡像作為格林函數。通過計算結構物在兩種模態波浪入射時的激振力以及能量守恒關系式，驗證所建立頻域數學模型的正確性。運用所建立模型計算了某一模態波浪入射時自由水面上箱型結構的反射及透射系數，從而得到了自由水面及內界面上的波動變化情況。根據計算結果可以得出結論如下：1）無論入射波為何種模態，隨著入射頻率的變化，在某些頻率處內界面上的擾動為零，零點的位置只與B/L2的值相關，而與上下層密度比、水深比及箱體寬度沒有明顯關系。并且表面波模態波浪入射時，內界面上內波模態波浪的反射和透射系數的變化趨勢基本相同；2）內波模態入射波浪與結構物作用后，自由水面將產生非常微小的擾動；表面波模態入射波浪與結構物作用后，在內界面上則會產生大幅擾動，在高頻區域擾動趨于零。

作者：尚玉超 勾瑩 滕斌 鄭艷娜 單位：大連理工大學 海岸和近海工程國家重點實驗室 大連海洋大學海洋與土木工程學院