Computation of Turbulent Flow around Submarine Under Oblique Towing Condition

本研究以一艘模仿德國S209-1400型潛艦外型為對象，利用CFD方法進行斜航狀態時的流場計算，計算方法採用有限體積法，以結構化網格，計算穩態的黏性流場。在計算潛艦周圍流場之前，先對KVLCC2這艘油輪進行斜航時的穩態流場計算，因該油輪已有實驗資料，藉由與實驗資料比對，以便先針對斜航計算所採用之方法進行驗證。油輪的計算條件為弗勞得數(Froud number) 0.142以及雷諾數(Reynold number) 3.945x106，斜航的角度為平擺(yawing)0、3、6、9與12度。計算結果和實驗結果比較，在阻力方面誤差約在10%以內，在流場局部數值(船殼表面壓力分布、流場速度分布)方面則和實驗相當接近。 接著進行潛艦斜航狀況下的流場計算，潛艦的計算條件為雷諾數9.86x107。其中裸船殼的部分：縱搖方向斜航，計算-20~20等18個斜航角度；平擺方向斜航，計算0、1、2、3、4、5、6、7、10、15、20等11個斜航角度；潛艦(含帆罩、控制翼)部分：縱搖方向斜航，計算-15~20等8個斜航角度。計算結果顯示，潛艦的附屬物在潛艦斜航時，對於船體周圍流場以及船體阻力皆會產生重大影響。另外將計算結果與實驗結果相比較，在力與力矩係數與斜航角度關係方面不僅趨勢上吻合，數值上的誤差都在可接受的範圍。The turbulent flow field around the KVLCC2 tanker and the submarine are numerically simulated by solving the Reynolds Averaged Navier Stokes equations. Finite volume method is implemented in the RANS solver and the structured grids are used. For the KVLCC2 tanker in steady drift motion without the free surface, computations are carried out at angles of yaw from 0 to 12 degrees in model scale. The computational results are validated against the experimental data provided by the “CFD Workshop Tokyo 2005” in terms of various global and local quantities. All the results are in good agreement with experimental data. For the submarine in steady drift motion, the hull form of submarine is imitation of German S209-1400 submarine. Computations for bare hull are carried out at angles of pitch from -20 to 20 degrees and angles of yaw from 0 through 20 degrees. Additionally, simulations for the submarine with appendages including sail and control fins are performed at angle of pitch from -15 to 20 degrees. The simulated results show that the appendages (sail, control fins) have great influence on flow field and resistance of the submarine. The simulated results also compare with experimental data. The comparison result shows that the present computational approach predicts the outstanding force coefficients and moment coefficients well.摘 要 i ABSTRACT ii 一、 緒 論 1 1.1 求解RANS方程式在斜航流場計算之應用 2 1.2 紊流模型的選擇 3 1.3 網格的數目與分布 4 二、 控制方程式與數值方法 6 2.1 控制方程式 6 2.2 紊流模型 6 2.3 數值方法 8 2.4 船殼所受合力計算 9 三、 計算方法的驗證 10 3.1 實驗狀態 10 3.2 計算實例 12 3.3 數值網格建構 12 3.4 計算結果與討論 15 3.4.1 不同網格數目計算結果比較 15 3.4.2 網格分佈的調整 17 3.4.3 不同角度斜航時船體受力情形 19 3.4.4 不同角度斜航時船殼表面壓力係數分布 21 3.4.5 不同角度斜航時螺槳位置軸向速度分佈 25 3.4.6 採用不同紊流模型計算結果之比較 27 四、 潛艦斜航計算 28 4.1 潛艦幾何外型 28 4.2 計算實例 29 4.3 數值網格建構 30 4.3.1 裸船殼縱搖方向斜航計算之數值網格 30 4.3.2 裸船殼平擺方向斜航計算之數值網格 31 4.3.3 潛艦(含附屬物)縱搖方向斜航計算之數值網格 32 4.4 計算結果與討論 34 4.4.1 裸船殼縱搖斜航計算結果 34 4.4.2 裸船殼平擺方向斜航計算結果 38 4.4.3 潛艦(含附屬物)縱搖方向斜航計算結果 40 五、 結論 47 參考文獻 49 附錄A 51 附錄B 5