Ityp = 0

ブロックフォーマットのキーワードこの材料則を使用すれば、停滞点からデータを指定することにより、気体流入条件をモデル化することができます。気体は理想気体とします。入力カードは/MAT/LAW11 (BOUND)と似ていますが、乱流パラメータを定義するための2つの新しい行を導入します。

フォーマット

(1)	(3)	(5)
/MAT/B-K-EPS/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}^{stagnation}$	$ρ_{0}^{stagnation}$
`Ityp`	`P`_sh	`Fscale`_T

Ityp =0: 気体流入（停滞点データから）

(1)	(2)	(3)	(5)	(6)	(7)
`node_ID`_v		$γ$			`C`_d
`fct_ID`_$ρ$
`fct_ID`_p		$P_{0}^{stagnation}$
空白のフォーマット
$ρ_{0} κ_{0}$		$ρ_{0} ε_{0}$	`fct_ID`_k	`fct_ID`_$ε$
$C_{μ}$		$σ_{κ}$	$σ_{ε}$		$P_{r} / P_{r t}$
`fct_ID`_T	`fct_ID`_v

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
$ρ_{i}^{stagnation}$	初期停滞密度（実数）	$[\frac{kg}{m^{3}}]$
$ρ_{0}^{stagnation}$	E.O.S（状態方程式）で使用される基準密度デフォルト $ρ_{0}^{stagnation} = ρ_{i}^{stagnation}$ （実数）	$[\frac{kg}{m^{3}}]$
`Ityp`	境界条件タイプ 1 = 0 気体流入（停滞点データから） = 1 液体流入（停滞点データから） = 2 一般流入 / 流出 = 3 サイレント境界（整数）
`P`_sh	圧力シフト 2 （実数）	$[Pa]$
`Fscale`_T	時間スケールファクター 3 （実数）
`node_ID`_v	速度計算用の節点識別子 3 = 0 $v_{i n} = \min_{n o d e \in f a c e} (v_{n o d e} \cdot n)$ > 0 $v_{i n} = ‖ v_{n o d e_I D} ‖$ （整数）
$γ$	理想気体定数（実数）
`C`_d	流量係数 5 デフォルト = 0.0（実数）
`fct_ID` $ρ$	停滞密度用の関数 $f_{ρ} (t)$ 識別子 3 = 0 $ρ^{stagnation} (t) = ρ_{i}^{stagnation}$ > 0 $ρ^{stagnation} (t) = ρ_{i}^{stagnation} \cdot f_{ρ} (t)$ （整数）
`fct_ID`_p	停滞圧力用の関数 $f_{P} (t)$ 識別子 3 = 0 $P^{stagnation} (t) = P_{0}^{stagnation}$ > 0 $P^{s t a g n a t i o n} (t) = P_{0}^{s t a g n a t i o n} \cdot f_{P} (t)$ （整数）
$P_{0}^{stagnation}$	初期停滞圧力 3 （実数）	$[Pa]$
$ρ_{0} κ_{0}$	初期乱流エネルギー（実数）	$[J]$
$ρ_{0} ε_{0}$	初期乱流散逸（実数）	$[J]$
`fct_ID`_k	乱流モデリングの関数 $f_{κ} (t)$ 識別子 = 0 $κ = κ_{adjacent}$ > 0 $κ = κ_{0} \cdot f_{κ} (t)$ （整数）
`fct_ID`_$ε$	乱流モデリングの関数 $f_{ε} (t)$ 識別子 = 0 $ε = ε_{adjacent}$ > 0 $ε = ε_{0} \cdot f_{ε} (t)$ （整数）
$C_{μ}$	乱流粘性係数デフォルト = 0.09（実数）
$σ_{κ}$	`k`パラメータの拡散係数デフォルト = 1.00（実数）
$σ_{ε}$	$ε$ パラメータの拡散係数デフォルト = 1.30（実数）
$P_{r} / P_{r t}$	層流プラントル数（デフォルトは0.7）と乱流プラントル数（デフォルトは0.9）の比率（実数）
`fct_ID`_T	流入温度用の関数 $f_{T} (t)$ 識別子 3 6 = 0 `T = T`_v = n $T = T_{0} \cdot f_{T} (t)$ （整数）
`fct_ID`_Q	流入熱流束用の関数 $f_{Q} (t)$ 識別子 3 6 = 0 強制流束なし = n $Q = f_{Q} (t)$ （整数）

停滞点から指定された気体状態 $(ρ_{stagnation}, P_{stagnation})$ が流入気体状態の計算に使用されます。
全エンタルピー定式化、断熱則、および状態方程式を含む一連の式により流入状態の完全な定義が可能となります：(1)
$ρ_{in} = ρ_{stagnation} \cdot {[1 - \frac{γ - 1}{2 γ} \cdot \frac{ρ_{stagnation}}{P_{stagnation}} \cdot (1 + C_{d}) \cdot ν_{in}^{2}]}^{\frac{1}{γ - 1}}$
(2)
$P_{in} = P_{stagnation} {(\frac{ρ_{in}}{ρ_{stagnation}})}^{γ}$
(3)
${(ρ e)}_{in} = \frac{P_{in}}{γ - 1}$
P_SHパラメータを使用すれば、P-P_SHにもなる出力圧力をシフトさせることができます。 $P_{s h} = P (t = 0)$ を使用している場合は、出力圧力が $Δ P$ （初期値は0.0）になります。
関数が定義されていない場合は、関連する量 $P_{s t a g n a t i o n}, ρ_{s t a g n a t i o n}, T$ またはQが一定になり、初期値に設定されます。ただし、すべての入力量 $P_{s t a g n a t i o n}, ρ_{s t a g n a t i o n}, T$ およびQは、指定された関数識別子を使用して時間依存関数として定義できます。横軸関数は、 $f (t)$ ではなく $f (F s c a l e_{t}, t)$ の使用につながるFscale_Tパラメータを使用してスケーリングすることもできます。
流入速度 $v_{i n}^{}$ はBernoulliの定理で使用されます（固定速度）。
流量係数は、入力損失が考慮されており、形状開口部によって異なります。

図 1.
熱モデリングを使用すれば、すべての熱データ（ $T_{0}, ρ_{0} C_{p}$ , …）を/HEAT/MATを使って定義できます。
この境界材料則は、多相材料ALE /MAT/LAW37 (BIPHAS)および/MAT/LAW51 (MULTIMAT)と一緒に使用することができません。

Ityp = 0

フォーマット

定義

コメント