应用体积力方法进行船舶於波浪中之自推模拟.PDF

第六屆海峽兩岸水動力學研討會 應用體積 力方法進行船舶 於波浪中 之自 推模 擬 1 1 1 1 2 2 黃斯寬 , 辛敬業, 呂玲 , 毛耀堂 , 陳紀川 , 林宗岳 (1. 國立臺灣海洋大學 系統工程暨造船學系 ; 2. CR Classification Society, Taipei, Taiwan) 摘要 “能源效率設計指標”(EEDI, Energy Efficiency Design Index) 中的氣候因子f w (weather factor)是一用來表示船舶因風浪失速的無因次化係數，而要進行船舶於波浪中的失速計算，就 必須進行船舶於靜水與波浪中的自推模擬，本文主要在探討應用體積力方法進行船舶在波浪 中自推之模擬。在2015 年的“第五屆海峽兩岸水動力學研討會”中，我們進行了“應用計算 方法預測船舶於波浪中之失速”的報告，進行了一些初步計算成果的報告，在本文中，我們 將進行更多結果的報告，以及不同體積力方法的應用等。體積力方法整合了黏性流計算程式 與勢流螺槳計算程式，用來模擬螺槳與船體互動的問題。本文在自推模擬上，螺槳的效應分 別以定常流與非定常流的體積力方法進行計算，船體阻力的計算方法分為三種。第一種方法， 靜水中船體的阻力直接應用黏性流 RANS 方法配合勢流方法進行興波阻力的計算，而波浪中 的波浪附加阻力則應用勢流截片理論計算；第二種方法，靜水中船體的阻力直接應用考慮自 由液面效應的黏性流 RANS 方法進行計算，而波浪附加阻力則仍應用勢流截片理論計算；第 三種方法，則是靜水與波浪中的船體阻力皆以考慮自由液面效應的黏性流 RANS 方法進行計 算，文中將比較這三種方法的計算時間與計算結果的準確度，以及對於流場的影響。 關鍵字 自推試驗，推進效率，船舶運動，失速，體積力方法，RANS 1 引言 “能源效率設計指標”(EEDI, Energy Efficiency Design Index) 中的氣候因子 f w (weather factor)是一用來表示船舶因風浪失速的無因次化係數，而要進行船舶於波浪中的失速計算，就 必須進行船舶於靜水與波浪中的自推模擬，本文主要在探討應用體積力方法進行船舶在波浪 中自推之模擬，而此模擬計算我們分為兩種計算方式，第一種方法應用黏性流對疊船模方法 計算船舶與螺槳互動之阻力，勢流興波阻力程式計算興波阻力，波浪附加阻力也應用勢流方 法進行計算；第二種方法則應用黏性流方法計算船舶與螺槳互動之阻力以及興波阻力，波浪 附加阻力仍然應用勢流方法進行計算，第三種方法，則是靜水與波浪中的船體阻力皆以考慮 自由液面效應的黏性流 RANS 方法進行計算。藉由靜水與波浪中二自推試驗模擬方法在三個 不同船速的模擬結果，即可進一步的計算船舶在波浪中航行時的失速。在 2015 年的“第五屆 海峽兩岸水動力學研討會”中，我們進行了“應用計算方法預測船舶於波浪中之失速”的報 告，進行了一些初步計算成果的報告，在本文中，我們將進行更多結果的報告。 - 1 - 第六屆海峽兩岸水動力學研討會 有關船舶失速的相關研究，Journée 在 1976 年即進行了研究[1][2] ，應用近似的計算方法考慮 船舶航行中船舶運動對於螺槳推進性能的影響，同時進行了實驗。Faltinsen[3]對於航行中阻力與推 進的整體變化進行了研究。日本的學者對於此問題進行了一系列的研究， Nakamura[4]與 Naito[5] 探討了螺槳於航行時負荷的變動， Nakatake[6]應用 source 與 sink 模擬船殼、螺槳與舵，並應用簡 單的黏性修正，以計算方法探討波浪對於螺槳、舵以及船殼交互作動的影響，Ando[7]隨後並進行 了實驗驗證數值計算的結果；並且應用了改進的以 Rankine source 為基礎的小板法推估波浪中船舶 的興波阻力與螺槳的推力[8] 。Kashiwagi[9]則由船舶運動理論，藉由所謂的 EUT(Enhanc