Ityp = 1

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

フォーマット

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

Ityp = 1：液体流入（停滞点データから）

(1)	(2)	(3)	(5)	(6)	(7)
`node_ID`_v		`C`₁			`C`_d
`fct_ID`_$ρ$
`fct_ID`_p		$P_{0}^{stagnation}$
`fct_ID`_E		$E_{0}^{stagnation}$
$ρ_{0} κ_{0}$		$ρ_{0} ε_{0}$	`fct_ID`_k	`fct_ID`_$ε$
$C_{μ}$		$σ_{κ}$	$σ_{ε}$		$P_{r} / P_{r t}$
`fct_ID`_T	`fct_ID`_Q

定義

フィールド	内容	SI 単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子. （整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
$ρ_{i}^{stagnation}$	初期停滞密度 3 （実数）	$[\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 （実数）	$[s]$
`C`₁`node_ID`_v	速度計算用の節点識別子 4 = 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`₁	液体の弾性率 9 （実数）
`C`_d	流量係数 5 デフォルト = 0.0（実数）
`fct_ID`_$ρ$	停滞密度用の関数 $f_{ρ} (t)$ 識別子 3 = 0 $ρ^{stagnation} (t) = ρ_{i}^{stagnation}$ > 0 $ρ^{s t a g n a t i o n} (t) = ρ_{i}^{s t a g n a t i o n} \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]$
$E_{0}^{stagnation}$	停滞点での初期比体積エネルギー 3 8 （実数）	$[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 $E^{s t a g n a t i o n} (t) = E_{0}^{s t a g n a t i o n} \cdot f_{E} (t)$ （整数）
$C_{μ}$	乱流粘性係数デフォルト = 0.09（実数）
$σ_{κ}$	`k`パラメータの拡散係数デフォルト = 1.00（実数）
$σ_{ε}$	$ε$ パラメータの拡散係数デフォルト = 1.30（実数）
$P_{r} / P_{r t}$	層流プラントル数（デフォルトは0.7）と乱流プラントル数（デフォルトは0.9）の比率（実数）
`fct_ID`_T	流入温度用の関数 $f_{T} (t)$ 識別子 = 0 `T = T`_neighbor = n $T = T_{0} \cdot f_{T} (t)$ （整数）
`fct_ID`_Q	流入熱流束用の関数 $f_{Q} (t)$ 識別子 = 0 強制流束なし = n $Q = f_{Q} (t)$ （整数）

停滞点から指定された気体状態 $(ρ_{stagnation}, P_{stagnation})$ が流入気体状態の計算に使用されます。その後で、Bernoulliが適用されます。
(1)
$P_{stagnation} = P_{in} + \frac{ρ_{in} ν_{in}^{2}}{2}$

これが次の流入状態につながります：(2)
$ρ_{in} = \frac{C_{1} \cdot ρ_{stagnation}}{C_{1} + \frac{ρ_{stagnation} ν_{in}^{2}}{2} (1 + C_{d})}$
(3)
$P_{in} = P_{stagnation} - \frac{ρ_{stagnation} ν_{in}^{2}}{2} (1 + C_{d})$
(4)
${(ρ e)}_{in} = (1 - \frac{ρ_{in}}{ρ_{stagnation}}) P_{in} + E_{stagnation}$
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パラメータを使用してスケーリングすることもできます。
流入速度 $ν_{in}$ は、Bernoulliの定理で使用されます。
流量係数は、入力損失が考慮されており、形状開口部によって異なります。

図 2.
熱モデリングを使用すれば、すべての熱データ（ $T_{0}, ρ_{0} C_{p}$ , …）を/HEAT/MATを使って定義できます。
この境界材料則は、多相材料ALE /MAT/LAW37 (BIPHAS)および/MAT/LAW51 (MULTIMAT)と一緒に使用することができません。
停滞エネルギーの定義は省略可能です。次のデフォルト値が推奨されます： $E_{0}^{stagnation} = 0.0$ 。線形EOS $Δ P = C_{1} μ$ はエネルギーに依存しないため、圧力は影響を受けず、初期エネルギーもユーザーによって設定されます。
固有体積エネルギーEは $= E_{i n t} / V$ として定義されます。

ここで、

$E_{i n t}$

内部エネルギー。これは、/TH/BRICを使用して出力できます。

固有質量エネルギーeは $e = E_{i n t} / m$ として定義されます。これが $ρ e = E$ につながります。固有質量エネルギーeは、/ANIM/ELEM/ENERを使用して出力できます。これは、ユーザーモデリングによって相対エネルギーになる場合があります。
液体体積弾性率は、通常、 $C_{1} = ρ_{0} c_{0}^{2}$ に設定されます。
ここで、

$c_{0}$

音速

Ityp = 1

フォーマット

定義

コメント