数值模式发展与应用(一).PDF

會自由液面之三維不可壓縮黏性流 數值模式發展與應用(一) 鍾孟軒* 一、前言 TAO 在自然界中存在著許許多多含自由液面之流場，譬如河川之水流、海洋之波浪、湖泊之水面等。 在水利工程問題上，吾人亦常可發現含自由液面流之例子，如明渠變量流(水躍)之水面急遲變化， 水流經過水工結構物(橋驗、溢洪道、握、丁壩、堤防、間門等等)之流場、水庫水面波動、港池振 動、進水口捲氣漩渦、電廠冷卻水海洋放流等，比比皆是。因此自由液面流是一個相當重要而值得探 討之研究主題。 合自由液面的流場問題相當複雜而困難 ， 因此不管在數學或數值模式上都存在著許多簡化程度不 一之近似模式 ， 如波浪問題中的線性波 、 靜水壓假設下之明渠流或俗稱之淺水波方程(Shallow Water Equations) ，雖然這些簡化模式在某些工程應用上有其實際之價值及合理之適用性，但對於實際複雜之 工程問題(如第一段文章所舉之流場)則無法適用 。 然而一般工程師在應用這些簡化模式時，常有意或無 意地忽視了其理論與應用限制，而延伸或高估了這些簡化模式之功能，也常因此低估了原本物理問題 之複雜與困難度。基於此，本計畫之主要目的有二 : (l)發展一套求解三維(紊流) 自由液面流場之數 值模式; (2) 廣泛應用於各類具自由液面水利工程問題之模擬 。 三、研究方法 本研究計畫之最終目標在發展適用三維自由液面計算之數值模式，期能建立有效之水工結構物三 維水理分析技術與工具。研究計畫之進行將先從二維模式發展，以便更容易熟悉分析及評估所選用之 數學物理模式與數值技巧，二維模式發展完成後再進行三維模式之發展。同時模式之發展流程將由簡 入繁，先從定常性層流自由液面問題切入，再發展至時變性紊流自由液面問題。 本研究計畫在自由液面計算方法上，原則上將朝等位函數法 ( Level Set Method) 以及流體體積追 蹤法 ( Volume of Fluid , VOF) 兩種做法同時發展並彼此驗證 ， 但限於人力 ， 後者僅發展至可處理二維 簡單幾何邊界的問題，而在流場解子方面，本研究計畫採取目前廣泛使用之不可壓縮流解法， SIMPLE句pe 方法 。 紊流模式採用LES 方法中最省計算機時間且最簡單之定係數次網格尺度模型 ( Constant Coefficient Subgrid-Scale, CSGS) *中興工程顧問社土木工程研究中心研究員 • 93 • 本研究計畫致力於數值模式之發展與分析程式之建立，電腦程式以FORTRAN 90 語言撰寫 ， 而以 Compaq Visual Fortran Professional Edition 6. 1.0 編譯 。 至於相關之網格建立與計算結果繪圖顯像之前 、 後處理工作，直接借助現有可用之(商業)軟體，如 SEC-HY2 l及 TECPLOT。 三、研究成果 本研究計畫之成果主要以案例計算的方式呈現，並儘可能對流體體積追蹤法及等位函數法做相互 驗證，所執行的主要模擬案例包括二維乾床潰壩問題之模擬、簡單海嘯傳播問題之模擬、孤立波之撞 壁反射、水底有三角形障礙物之水槽漾 盪、具地面障礙物之二維潰壩問題、水平 Turbulent