空蝕模型理論

Creo Flow Analysis > 預處理 > 定義實體 > 氣蝕 > 物理 > 空蝕模型理論

空蝕模型理論

在蒸汽傳輸方程式 2.169 中，Re 與 Rc 是與空蝕流中蒸汽氣泡的增長與塌陷相關的質量傳遞源項。這些項用於解釋空蝕過程中汽相與液相之間的質量交換。Re 與 Rc 會基於用來描述單一蒸汽氣泡在液體中的增長的 Rayleigh-Plesset 方程式進行建模。

液汽質量傳遞

欲推導空蝕中淨相變率的運算式，請考慮使用滑動速度為零的不可壓縮液汽二相流 (不考慮不凝結氣體)。如果引入 R 表示從液體到蒸汽的淨質量傳遞率，則液體與蒸汽體積分數方程式與總質量連續性方程式如下所示：

• 液相：

方程式 2.179

• 汽相：

方程式 2.180

• 混合物 (總質量連續性)

方程式 2.181

在液汽二相系統中，可透過蒸汽體積分數與相密度表示混合物密度 ρ：

方程式 2.182

由於假設液體與蒸汽密度恆定 (不可壓縮)，因此可透過方程式 2.181 和方程式 2.182 推導速度梯度與蒸汽體積分數之間的關係︰

方程式 2.183

合併方程式 2.179 與方程式 2.183，所得淨質量源項的運算式如下：

方程式 2.184

將方程式 2.184 代入方程式 2.180，蒸汽體積分數方程式可以改寫為一般式：

方程式 2.185

透過套用蒸汽質量分數與體積分數之間的關係方程式 2.173，會取得以蒸汽質量分數表示的方程式 2.185︰

方程式 2.186

根據方程式 2.185 與方程式 2.186，可瞭解到，當液相與汽相之間的滑動速度為零時，可透過額外的蒸汽質量分數傳輸方程式將空蝕建模為單相流；或透過液汽質量傳遞建模為尤拉多相混合流。如果不考慮擴散與相速度差的影響，這兩個方法在數學運算方面相同。Creo Flow Analysis 使用單相方法對空蝕流進行建模。

氣泡動力學考量

在大多數自然現象及工程系統中，液體中存在足夠數量的核 (例如氣泡、不凝結氣體等)，從而為空蝕初生提供了條件。因此，欲對空蝕過程進行建模，應主要關注正確解釋氣泡增長與塌陷過程。假設流動液體中液體與蒸汽氣泡之間的滑動速度為零，可透過描述液體中氣體氣泡增長的廣義 Rayleigh Plesset 方程式來推導氣泡動力學方程式：

方程式 2.187

其中，

RB	氣泡半徑
pB	氣泡內部壓力 (不存在其他氣體時，假設為液體溫度下的蒸汽壓力)
p	液體中氣泡周圍的壓力
σ	液體與蒸汽之間的表面張力係數

此方程式透過機械平衡 (沒有阻礙氣泡增長的熱障礙) 推導得出。

如果忽略二階時間導數 (適用於低振盪頻率)、黏性阻尼項與表面張力，則會取得適用於漸進狀態之方程式 2.187 的簡化運算式：

方程式 2.188

此簡化 Rayleigh Plesset 方程式提供用於在空蝕模型中引入氣泡動力學效果的物理方法。氣泡半徑可增加或減小，具體取決於 (pB–p) 的符號︰當 p<(pB 時，氣泡會增長；當 p>pB 時，氣泡會塌陷。因此，方程式 2.188 可改寫為如下所示：

方程式 2.189

如果 η 為液體中的蒸汽氣泡數密度 (單位體積內存在的氣泡數)，並且所有蒸汽氣泡都是具有相同半徑 RB 的完美球體，則汽相體積分數如下：

方程式 2.190

假設蒸汽氣泡無法在液體中形成或破裂，但氣泡可以在空蝕過程中增長 (蒸發) 和塌陷 (凝結)。在方程式 2.190 中，蒸汽氣泡數密度 (η) 保持恆定，但氣泡半徑 (RB) 將增大或減小。蒸汽體積分數的時間導數計算如下

方程式 2.191

將方程式 2.191 代入方程式 2.186，取得用於控制汽相質量分數的傳輸方程式：

方程式 2.192

透過套用方程式 2.189，單位體積內液體與蒸汽之間的淨質量傳遞率形式如下：

方程式 2.193

其中方程式 2.193表示在空蝕中，單位體積質量傳遞率 (R) 是汽相與液相密度的函數 (成正比)，並且還與混合物密度成反比。由於方程式 2.192 是直接從相與混合流質量連續性推導得到的，因此它是精確的並且應能夠準確地表示空蝕中液相與汽相之間的質量傳遞。引入氣泡動力學之後，方程式 2.193 會使用類似方法對兩個相反 (即在物理上不同) 的質量傳遞過程進行建模：從液體到蒸汽 (氣泡增長或蒸發)，以及從蒸汽到液體 (氣泡塌陷或凝結)。對於蒸汽質量分數傳輸方程式 2.192，氣泡增長為源項，而氣泡塌陷則視為匯項。

在實際空蝕模型中，局部遠場壓力 p 通常視為與单元中心壓力相同。如果不存在溶解氣體、質量傳遞與黏性阻尼，則氣泡壓力 pB 等於飽和蒸汽壓力 (psat，材料屬性)，即 pB=psat

如果將方程式 2.192 和方程式 2.193 與一般蒸汽質量分數方程式 2.169 相比較，源項 Re 與 Rc 如下：

• 如果局部流動壓力低於飽和蒸汽壓力，即 p<psat，則僅會發生蒸發，因此：

方程式 2.194

• 如果局部流動壓力大於飽和蒸汽壓力，即 p>psat，則僅存在凝結：

方程式 2.195

幾乎所有可用機械二相空蝕模型都基於方程式 2.194 與方程式 2.195。Creo Flow Analysis 使用了 Singhal 等人開發的建模方法。

氣體吸收/溶解與釋放

工作流體中通常存在不凝結氣體，這些氣體會對空蝕造成重大影響。有時，不凝結氣體會隨著流體流動自由傳輸，也會在液體中溶解或釋放，本質上是為了在液相和氣相之間實現質量濃度動態平衡。液體中的氣體吸收或溶解與釋放也是一種液氣質量傳遞現象，這種現象由質量濃度差和梯度驅動。為了對空蝕流進行建模，還需要考慮混合流中不凝結氣體及可能的液氣質量傳遞的影響。

假設在液氣二相流中，液相 (溶解氣體) 與氣相 (自由氣體) 中均存在不凝結氣體 (例如空氣或氧氣)，每個相中的氣體質量分數傳輸方程式如下：

• 自由氣體 (氣相)

方程式 2.196

• 溶解氣體 (液相)

方程式 2.197

其中，

fg,f	自由氣體的質量分數
fg,d	自由氣體與溶解氣體的質量分數
Sg,f, Sg,d	外部源或使用者定義的源。
Dg,f	自由氣體與溶解氣體的擴散係數
Dg,d	溶解氣體的擴散係數

如果不凝結氣體的質量分數預先規定為 fg，則可得：

fg=fg,f

或者，可透過求解下列方程式取得其在空間和時間上的傳輸：

方程式 2.199

請注意，對於方程式 2.196、方程式 2.197 與方程式 2.199，應僅直接求解其中的兩個方程式。

在方程式 2.196 與方程式 2.197 中，來源 Rd→f 表示溶解氣體釋放率，Rf→d 表示自由氣體吸收或溶解率。

當兩個相接觸時，易發生自由氣體 f 與溶解氣體 d 從一個相傳輸至另一相，以實現二相之間的動態平衡。平衡模型假設體積質量傳遞率取決於質量濃度差或梯度：

方程式 2.200

方程式 2.201

其中，

AI	液氣界面面積
kf,d(=kd,f)	整體體積質量傳遞係數
ρg,d(=ρfg,d)	溶解氣體的局部質量濃度
ρg,f(=ρfg,f)	自由氣體的局部質量濃度
	溶解氣體在其主相中的平衡質量濃度
	自由氣體在其主相中的平衡質量濃度