Ityp = 1

ブロックフォーマットのキーワードこの材料則を使用すれば、停滞点からデータを指定することにより、液体流入条件をモデル化することができます。液体挙動は線形EOSを使用してモデル化されます。

フォーマット

(1)	(3)	(5)
/MAT/LAW11/`mat_ID`/`unit_ID`または/MAT/BOUND/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}^{stagnation}$	$ρ_{0}^{stagnation}$
`Ityp`	`P`_sh	`Fscale`_T

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

(1)	(2)	(3)	(7)
`node_ID`_V		`C`₁	`C`_d
$fct_{ID}_{ρ}$
`fct_ID`_p		$P_{0}^{stagnation}$
`fct_ID`_E		$E_{0}^{stagnation}$
空白のフォーマット
空白のフォーマット
`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]$
`node_ID`_V	速度計算用の節点識別子 4 = 0 $v_{in} = \min_{nod e ε face} {\vec{v}}_{n o d e} \cdot \vec{n}$ > 0 $v_{in} = ‖ v_{node_ID} ‖$ （整数）
`C`₁	液体の弾性率 9 （実数）
`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^{stagnation} (t) = P_{0}^{stagnation} \cdot f_{P} (t)$ （整数）
$P_{0}^{stagnation}$	初期停滞圧力 3 （実数）	$[Pa]$
`fct_ID`_E	停滞密度用の関数 $f_{E} (t)$ 識別子 3 = 0 $E^{stagnation} (t) = E_{0}^{stagnation}$ > 0 $E^{stagnation} (t) = E_{0}^{stagnation} \cdot f_{E} (t)$ （整数）
$E_{0}^{stagnation}$	停滞点での初期比体積エネルギー 3 8 （実数）	$[Pa]$
`fct_ID`_T	流入温度用の関数 $f_{T} (t)$ 識別子 3 6 = 0 $T = T_{adjacent}$ = 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})$ が流入気体状態の計算に使用されます。その後で、Bernoulliが適用されます。
(1)
$P_{stagnation} = P_{in} + \frac{ρ_{in} v_{in}^{2}}{2}$

これが次の流入状態につながります：(2)
$\begin{array}{l} ρ_{in} = \frac{C_{1} \cdot ρ_{stagnation}}{C_{1} + \frac{ρ_{stagnation} v_{in}^{2}}{2} (1 + C_{d})} \\ P_{in} = P_{stagnation} - \frac{ρ_{stagnation} v_{in}^{2}}{2} (1 + C_{d}) \\ {(ρ e)}_{in} = (1 - \frac{ρ_{in}}{ρ_{stagnation}}) P_{in} + E_{stagnation} \end{array}$
P_shパラメータを使用すれば、P-P_shにもなる出力圧力をシフトさせることができます。P_sh=P(t=0)を使用している場合は、出力圧力が $Δ P$ （初期値は0.0）になります。
関数が定義されていない場合は、関連する量 $(P_{stagnation}, ρ_{stagnation}, T, or Q)$ が一定になり、初期値に設定されます。ただし、すべての入力量 $(P_{stagnation}, ρ_{stagnation}, T, and Q)$ は、指定された関数識別子を使用して時間依存関数として定義できます。横軸関数は、f(t)ではなくf (Fscale_t* t)の使用につながるFscale_Tパラメータを使用してスケーリングすることもできます。
流入速度 $v_{in}$ は、Bernoulliの定理で使用されます。
流量係数は、入力損失が考慮されており、形状開口部によって異なります。

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

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

Ityp = 1

フォーマット

定義

コメント