フローモデル
質量および運動量保存のための「フロー」(Flow) モジュールの解法 (非定常 Navier-Stokes 方程式を使用) H.Ding、F.C. Visser、Y.Jiang、および M. Furmanczyk「Demonstration and Validation of a 3-D CFD Simulation Tool Predicting Pump Performance and Cavitation for Industrial Applications」FEDSM2009-78256、2009
Reynold's Averaged Navier-Stokes 方程式 (RANS) の積分形式 (保存) は次のとおりです。
連続性
運動量
応力テンソル
ここで、
τij
有効なせん断応力 (分子 + 乱流)
f
ボディのフォース
n
サーフェス法線
ρ
静的圧力 (Pa)
t
時間
v
流体速度
vσ
メッシュ速度
Ω(t)
時間の関数として制御ボリューム
r
ローカル流体密度の平均 (kg/m3)
σ
制御ボリュームのサーフェス
µ
動粘度 (ポアズまたは Pa-s)
µt
乱流の動粘度
δij
クロネッカーのデルタ (i=j では =1、i≠j では =0)
粘度モデル
「定数動的粘度」(Constant Dynamic Viscosity) - 選択した体積の流体の粘度を指定します。動粘度の単位は、Pa-s または N-s/m2 です。
動粘度の値は、「定数動的粘度」(Constant Dynamic Viscosity) 選択にあるボックスで指定します。
「定数キネマティック粘度」(Constant Kinematic Viscosity) - 選択した体積の流体の粘度を指定します。キネマティック粘度の単位は m2/s です。キネマティック粘度の値は、「Constant Kinematic Viscosity」選択にあるボックスで指定します。
「Sutherland Law」 - 選択した体積の流体の粘度を、動粘度 (Pa-s) で指定します。方程式と入力は次のとおりです。
ここで、
T
温度 (K)
µref
参照温度の粘度 (Pa-s)
S
サザランド温度 (K)
* 
T は、エネルギーモジュールがアクティブでない場合に、入力として必要とされる流体温度 (K) です。
サザランド則を使用して、理想気体の粘度を温度の関数として計算します。Sutherland, W. (1893), "The viscosity of gases and molecular force," Philosophical Magazine, S. 5, 36, pp. 507-531 (1893)。次の表に、選択した気体に対するサザランドの定数と参照温度を示します。参照: en.wikipedia.org/wiki/viscosity
気体
S (K)
Tref (K)
mref (Pa-s)
空気
120
291.15
18.27 e-6
窒素
111
300.55
17.81 e-6
酸素
127
292.25
20.81 e-6
二酸化炭素
240
293.15
14.8 e-6
一酸化炭素
118
288.15
17.2 e-6
水素
72
293.85
8.76 e-6
アンモニア
370
293.15
9.82 e-6
二酸化硫黄
416
293.65
12.54 e-6
ヘリウム
79.4
273
19 e-6
非ニュートン粘度モデル
非ニュートン粘度モデルは次のとおりです。
Herschel-Bulkley モデル
Bingham モデル
これらのモデルは、非ニュートン流れのプロパティを示すさまざまなタイプの流体に対して適切な粘度を提供します。Herschel-Bulkley モデルと Bingham モデルは、せん断応力とせん断率を次のように関連付けます。
ここで、
e0
重要なせん断率
k
一貫性指標
τ0
流体の降伏応力
n
べき乗則の指数Bingham モデルの場合、n =1
* 
せん断率 0 は、上のプロットのガンマ点と同じです。
抵抗モデル
「抵抗モデル」(Resistance Model) は、選択した体積で抵抗を設定するために使用できる流れモジュールオプションです。「抵抗モデル」(Resistance Model) には、次の 2 つのモデルが含まれています。
「圧力損失」(Pressure Loss) - 次の方程式に基づきます。
ここで、
Cl
直線ドラッグ係数 (Pa-s/m2)
Cd
2 次ドラッグ係数 (1/m)
β
多孔性
ρ
密度
「ダルシーの法則」(Darcy's Law) - 次の方程式に基づくモデル。
ここで、
β
多孔性
α
透磁率
µ
動粘度
V
速度
Cd
2 次ドラッグ係数 (1/m)
抵抗方程式で使用される速度は、ローカル速度です。この方程式の F は、単位 N/m3 (荷重/体積など)、圧力勾配 (Dp/Dx)、または rg で測定されます。インタフェース全体の圧力損失は、F と有限厚を掛け合わせることで求められます。多孔性は「共通」(Common) モジュールで設定します。
これは役に立ちましたか?