數(shù)值波浪水槽及其應用研究.pdf

1、波浪是近岸海洋環(huán)境中最重要的動力因素之一，威脅工程結構物的安全，影響工程結構物的規(guī)劃、設計以及淺海區(qū)的生產(chǎn)作業(yè)。因此，近岸波浪與結構物相互作用以及數(shù)值波浪水槽的研究一直是水動力學的研究熱點之一。本文旨在開發(fā)一套適用于求解帶自由表面流動的兩相流計算程序，建立高效數(shù)值波浪水槽以用于波浪傳播以及與結構物相互作用的精細研究。
　　首先，建立了處理大變形自由液面問題的兩相流數(shù)學模型并對該模型的數(shù)值求解算法進行了系統(tǒng)的研究?？刂品匠贪臻g濾

2、波的NS方程、連續(xù)性方程、求解體積比率F的輸運方程以及使方程組封閉的附加方程如湍流模型方程等。自由液面的求解采用單相流數(shù)學模型進行描述，自由表面追蹤采用VOF方法，湍流模型采用大渦模擬中的Smagorinsky模型。
　　對兩相流數(shù)值模擬算法的研究包括以下幾個方面:(1)采用兩相流模型對空氣和水進行同時求解。系統(tǒng)研究了動量方程空間離散格式的選取，研究發(fā)現(xiàn)密度的處理不當會引起較大的數(shù)值耗散和不穩(wěn)定，因此，NS方程求解中涉及密度的地方

3、給予了謹慎處理，選擇了帶限制函數(shù)的TVD格式。(2)對自由面劇烈變形、破碎、空氣卷吸和不規(guī)則波浪傳播過程的模擬需要高效、穩(wěn)定和高精度的求解器。詳細研究了NS方程的求解算法，基于SIMPLE算法提出了對速度和壓力皆具有兩階時間精度的求解算法。(3)自由面捕捉精度是兩相流計算的核心部分，根據(jù)求解問題的特點研究了VOF方法所涉及的界面重構算法和界面輸運算法。為了解決自由表面追蹤方法中質量不守恒的現(xiàn)象，提出了質量重新分布技術。(4)考慮到笛卡爾

4、網(wǎng)格簡潔高效且有利于VOF方法的實現(xiàn)，采用浸入邊界法中的離散力法來處理結構物邊界，這種隱性處理壁面邊界條件的方法不僅降低了結構化網(wǎng)格系統(tǒng)中處理復雜幾何邊界的難度，且與VOF方法有很好的契合度。(5)為研究波浪對結構物拍擊問題中空氣可壓縮性的影響，將不可壓縮兩相流動的算法擴展到考慮空氣可壓縮的兩相流求解算法，給出了具體的模型推導和數(shù)值方法實現(xiàn)過程。(6)針對笛卡爾網(wǎng)格簡單易行的特點以及自由面求解需要較高的網(wǎng)格分辨率需求，初步研究了局部加密

5、的笛卡爾網(wǎng)格技術。該種網(wǎng)格結構可根據(jù)流場特點進行局部加密，增強了求解復雜問題的有效性。
　　針對所建立的數(shù)值求解器，采用了不同算例對所采用的數(shù)值方法及處理技術進行了驗證，主要考察了NS方程求解算法，處理自由表面的VOF算法以及VOF方法與NS方程求解耦合的兩相流求解器的正確性以及有效性。通過算例驗證了本文所開發(fā)的數(shù)值求解器可有效地、穩(wěn)定地求解帶自由表面的兩相復雜流動現(xiàn)象。
　　為了分析空氣可壓縮性對晃蕩問題結構受力的影響，分

6、別采用了不可壓縮流體的兩相流模型和視空氣為可壓縮氣體的兩相流模型模擬了晃蕩問題。數(shù)值結果表明空氣不可壓縮的兩相流模型不能反映包裹氣泡的影響，包裹氣泡對拍擊力的影響表現(xiàn)為壓力以一定頻率的振蕩。模擬了高充水晃蕩問題并與實驗結果進行了對比分析，從數(shù)值模擬角度驗證了采用空氣可壓縮模型對準確計算兩相流抨擊作用力是必要的。
　　基于所開發(fā)的求解器，建立了數(shù)值波浪水槽。將解析松弛法應用于求解NS方程的流場體系中，實現(xiàn)了耦合造波、消波功能，解決了

7、消除計算域中的二次反射問題。模擬規(guī)則波行波以及駐波驗證了造波、消波技術的有效性。為了進一步驗證數(shù)值波浪水槽的性能，分別模擬了一系列不同波參數(shù)的規(guī)則波和不規(guī)則波，并和理論結果進行了比較，通過大量的數(shù)值實驗給出了模擬非線性規(guī)則波和不規(guī)則波的基本網(wǎng)格尺度，對數(shù)值求解波浪傳播問題的網(wǎng)格收斂性進行了分析。
　　系統(tǒng)地研究了斜坡堤上的越浪問題。在規(guī)則波越浪的研究中，采用一系列不同堤頂高程和來波要素組合工況對越浪進行了數(shù)值模擬，將越浪量的計算結

8、果與已有的實驗結果和數(shù)值結果進行了比較，并擬合了越浪量和堤頂越浪流最大速度的關系。針對堤頂和背坡的沖刷問題，定量研究了后坡越浪流的流態(tài)和最大沖刷速度。通過系統(tǒng)數(shù)值模擬分析了非破碎和破碎波越浪后堤頂流深和最大流速沿程變化規(guī)律。堤頂流深沿程呈指數(shù)衰減，最大流速沿堤頂呈一定程度增大趨勢。
　　對于不規(guī)則波越浪的研究，針對不規(guī)則波作為平穩(wěn)各態(tài)歷經(jīng)隨機過程的特性，采用了一個工況多個波序列(5-7個)的方式進行越浪量統(tǒng)計計算。通過該途徑得到了