Ityp = 2

ブロックフォーマットキーワードこの材料則を使用すれば、その状態を直接設定することにより、材料流入 / 流出をモデル化することができます。入力カードは/MAT/LAW11 (BOUND)と似ていますが、乱流パラメータを定義するための２つの新しい行を導入します。

フォーマット

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

Ityp = 2：一般流入 / 流出

(1)	(2)	(3)	(5)	(6)	(7)
空白のフォーマット
`fct_ID` $ρ$
`fct_ID`_p		`P`₀
`fct_ID`_E		`E`₀
$ρ_{0} κ_{0}$		$ρ_{0} ε_{0}$	`fct_ID`_k	`fct_ID`_e
$C_{μ}$		$σ_{κ}$	$σ_{ε}$		$P_{r} / P_{r t}$
`fct_ID`_T	`fct_ID`_Q

定義

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

例（気体）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/B-K-EPS/3
GAS INLET (unit: kg_m_s)
#              RHO_I
               .3828
#     ITYP                           Psh         Fscale_T
         2
#blank line
#  fct_RHO
         1
#    fct_P                           P_0
         0
#    fct_E                           E_0
         1                        253300
#             Rho0k0            Rho0Eps0     fct_k   fct_eps
                  20                   0         1         0
#                Cmu             Sigma-k       Sigma-epsilon              Pr/Prt
                   0                   0                   0                   0
# fct_T        fct_Q
/ALE/MAT/3
#     Modif. factor.
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1
CST
#                  X                   Y
                   0                   1
              1.0E20                   1
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

指定された状態が流入境界要素に直接強制されます。これが次の流入状態につながります。
(1)
$ρ_{i n} = ρ_{i} f_{ρ} (t)$
(2)
$P_{i n} = P_{0} f_{P} (t)$
(3)
$E_{i n} = {(ρ e)}_{i n} = E_{0} f_{E} (t)$

この定式化を使用すれば、境界節点上の速度を物理流入速度（/IMPVEL）と一致させることができます。/MAT/LAW11 - Ityp=0および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パラメータを使用してスケーリングすることもできます。
熱モデリングを使用すれば、すべての熱データ（ $T_{0}, ρ_{0} C_{p}$ , …）を/HEATを使って定義できます。
この境界材料則は、多相材料ALE /MAT/LAW37 (BIPHAS)および/MAT/LAW51 (MULTIMAT)と一緒に使用することができません。
固有体積エネルギーEは $E = \frac{E_{int}}{V}$ として定義されます。
ここで、

$E_{i n t}$

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

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

Ityp = 2

フォーマット

定義

例（気体）

コメント