在本子計畫㆗的研究結果可包含兩部份，㆒部份是㆟機界面的實現，另㆒ 部份是模擬分析的結果與影響參數的探討。在程式流程㆖，是先以均質流模 型分析，並計算出噴嘴入口速度（Ui）、潛體入口直徑（da）、擴流器長度（Ld）、
需加入之氣體質量流率（
m
g）、噴嘴出口直徑與噴嘴喉部直徑比值（de/ d*）、
出口馬赫數（Me）及出口比推力（Isp），得到這些結果後，使用者可以選擇是 否使用雙流體模型進行修正，如果選擇「否」則程式將儲存結果並結束，若 選擇「是」程式則進入雙流體模型分析程式，程式會自動修正
d
a、L
d及 m
g 等 值，修正完後儲存結果並結束程式。均質流模型程式的操作條件及環境參數 輸入表單如㆘圖十。圖十、均質流模型程式的環境參數輸入表單
計算之後結果顯示如圖十㆒。
圖十㆒、均質流模型程式的計算結果
對於噴嘴之設計，在均質流模型計算㆗僅需知道潛體之操作參數及噴嘴的 入口直徑（di）和喉部直徑（d*），就可以計算出當噴嘴出口壓力（Pe）等於環 境壓力（Pa）時的出口直徑（de）以及出口速度（Ue），進而計算出口馬赫數
（Me）與出口比推力（Isp）等結果，並不需要噴嘴流道尺寸變化之詳細資料。
但是在使用雙流體模型計算時，必須要輸入噴嘴流道的幾何尺寸、噴嘴入口 之氣泡半徑等資料。其他參數則直接採用均質流模型計算之結果結果。雙流 體模型之目的在於修正兩相流體間，動力不平衡（dynamic non-equilibrium）
所造成之總質量流率的誤差，例如：滑動（slip）及氣泡體積因流場壓力改變 所呈現之非線性變化效應。使用雙流體模型程式所需額外增加的輸入條件及 參數如㆘圖十㆓、十㆔。
20
圖十㆓、雙流體模型程式所需額外增加的輸入條件:噴嘴資料
圖十㆔、雙流體模型程式所需額外增加的輸入條件:噴嘴入口參數
最後計算結果的輸出表單則如㆘圖十㆕。
圖十㆕、雙流體模型程式最後計算結果的輸出表單
至於在探討影響推進性能之設計參數部份，由於均質流模型遠比如Albagli and Gany [3]或 Wang and Chen[2]所用之雙流體模型要簡單許多，並且當噴嘴 流道變化幅度不太劇烈時，氣泡流可視為均質狀態，此時，均質流分析也有
㆒定的準確度，因此相當適合作為噴嘴之性能參數探討工具。故在本子計畫
㆗，皆以均質流模型為基礎，來探討設計參數對潛體推進性能之影響。而此 探討的目的將僅專注在如何設計㆒噴嘴截面積，以獲得最大的推力。
綜言之，影響潛體推進性能之參數有：
(1) 混合室之出口壓力：
混合室之出口壓力亦即噴嘴(Nozzle)的入口壓力。㆒般而言，噴嘴的入口 壓力越高，則引擎推力也越高；但在水㆘噴射推進之技術㆗，混合室的設計 是㆒個重要的關鍵。除了應盡可能保持由擴流器轉換之高壓外，加入的氣泡
22
比較引擎性能的最重要指標是比推力(Specific Impulse)，或是比衝值。由
㆘圖十五可以看出，(1) 加入氣泡之推進器其比推力皆比純水推進要高出許
圖十五、噴嘴出口比推力與出入口壓力比值之關係圖。
圖十六、噴嘴出口馬赫數與出入口壓力比值之關係圖。
圖十七則是對圖十六作㆒更精確的描繪，氣/液質量流率比（μ）及噴嘴出 入口壓力比值（Pa / Pi）對出口狀態（超音速或次音速）之關係。由圖十七可 看出，要使水動力推進達到效用，則入口壓力與氣液質量比的關係(μ, Pa / Pi)
24
必須落在虛線的㆘方。
圖十七、氣/液質量流率比（μ）及噴嘴出入口壓力比值（Pa / Pi）對出口狀態
（超音速或次音速）之關係圖。
(3) 壁面摩擦力的影響：
在本子計劃㆗，藉由若設定相同的入口條件，但分別求解包含與省略壁面 摩擦力的解，由分析結果可知只有在（Pa / Pi）較小時，㆓者間才會有些許的 差異，因為在出入口壓力比值（Pa / Pi）較小時，噴嘴之平均流速較大，因此 壁面摩擦效應的影響也較大，但㆓者相差仍在5 %以內。由於在此分析的效應 不大，故詳細的分析在此不贅述。
Reference:
1. Wang, Y.C. and Chen, E., Effects of phase relative motion on critical bubbly flows through a converging-diverging nozzle, Physics of Fluids A., vol. 14(9), pp.3215-3223.,（2002）.
2. Wang, Y.C. and Brennen, C.E., One-dimensional bubbly cavitating flows through a converging-diverging nozzle, ASME Journal of Fluids Engineering, vol. 120, pp.166-170., （1998）.
3. Albagli, D. and Gany, A., Analysis of a high speed bubbly water jet, AIAA/SAE/ASME/ASEE Joint Propulsion Conference., AIAA Paper 92-3380, 28th Nashville, TN, (1992).
26