Ityp = 3

ブロックフォーマットキーワードこの材料則を使用すれば、サイレント境界（NRF）をモデル化することができます。入力カードは/MAT/LAW11 (BOUND)と似ていますが、乱流パラメータを定義するための2つの新しい行を導入します。

フォーマット

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/MAT/B-K-EPS/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$		$ρ_{0}$
`Ityp`		`P`_sh

Ityp = 3：サイレント境界

(1)	(3)	(5)	(6)	(7)
	`c`	$l_{c}$
空白のフォーマット
空白のフォーマット
空白のフォーマット
$ρ_{0} κ_{0}$	$ρ_{0} ε_{0}$	`fct_ID`_k	`fct_ID`_e
$C_{μ}$	$σ_{κ}$	$σ_{ε}$		$P_{r} / P_{r t}$
空白のフォーマット

定義

フィールド	内容	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	圧力シフト 2 （実数）	$[Pa]$
`c`	流出音速 1 （実数）	$[\frac{m}{s}]$
$l_{c}$	特性長さ 1 （実数）	$[m^{3}]$
$ρ_{0} κ_{0}$	初期乱流エネルギー（実数）	$[J]$
$ρ_{0} ε_{0}$	初期乱流散逸（実数）	$[J]$
`fct_ID`_k	乱流モデリングの関数 $f_{κ} (t)$ 識別子 = 0 $κ = κ_{adjacent}$ > 0 $κ = κ_{0} \cdot f_{κ} (t)$ （整数）
`fct_ID`_e	関数 `f`_$ε$(`t`) エネルギーの識別子. = 0 $ε = ε_{adjacent}$ = n $ε = ε_{0} \cdot f_{ε} (t)$ （整数）
$C_{μ}$	乱流粘性係数デフォルト = 0.09（実数）
$σ_{κ}$	`k`パラメータの拡散係数デフォルト = 1.00（実数）
$σ_{ε}$	$ε$ パラメータの拡散係数デフォルト = 1.30（実数）
$P_{r} / P_{r t}$	層流プラントル数（デフォルトは0.7）と乱流プラントル数（デフォルトは0.9）の比率（実数）

例（気体）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/B-K-EPS/5
GAS OUTLET (unit: kg_m_s)
#              RHO_I               RHO_0         
               .3828                   0
#     ITYP                           Psh       
         3                           0.0 		 
#                                      c                  lc
                                     605                 0.3
#blank line

#blank line
                    
#blank line

#             Rho0k0            Rho0Eps0     fct_k   fct_eps
                  20                   0         0         0
#                Cmu             Sigma-k       Sigma-epsilon              Pr/Prt
                   0                   0                   0                   0  
#blank line

/ALE/MAT/5
#     Modif. factor.
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

サイレント境界定式化は、Bayliss & Turkelに基づきます。¹ この目的は、圧力と速度の急激な変化で変動する平均圧力を強制することです。
$\frac{\partial P}{\partial t} = ρ c (\frac{\partial}{\partial t} (V_{n}) - V_{n} \cdot div (V - V_{n} \cdot n)) + c \frac{(P_{\infty} - P)}{2 l_{c}}$
遠距離場の圧力 $P_{\infty}$ は時間関数を使用して設定されます。過渡圧力は、 $P_{\infty}$ 、局所速度場V、および出口面の法線から求めます：
- 密度、エネルギー、温度、乱流エネルギー、および散逸は、Ityp = 2のように時間関数で設定されます。
- 関数番号が0の場合は、連続性を尊重して隣接要素の値が使用されます。
- 音響インピーダンスは次のようになります： $ρ c$
- 一般的な長さ $l_{c}$ を使用して、設定された値に対する有効な圧力を緩和します。これは問題に関係する最高波長に比べて大きい値にする必要があります。緩和項がハイパスフィルタとして機能する場合の周波数カットオフは次のとおりです：(1)
  $f_{c} = \frac{c}{4 \cdot π \cdot l_{c}}$
ここで、音速cと特性長 $l_{c}$ は、2つの必須パラメータ（0以外）です。
P_shパラメータを使用すれば、P-P_shにもなる出力圧力をシフトさせることができます。 $P_{s h} = P (t = 0)$ を使用している場合は、出力圧力が $Δ P$ （初期値は0.0）になります。
熱モデリングを使用すれば、すべての熱データ（ $T_{0}, ρ_{0} C_{p}$ , …）を/HEAT/MATを使って定義できます。
この境界材料則は、多相材料ALE /MAT/LAW37 (BIPHAS)および/MAT/LAW51 (MULTIMAT)と一緒に使用することができません。

¹ A. Bayliss, E. Turkel, "Outflow Boundary Condition for Fluid Dynamics", NASA-CR-170367, Institute for Computer Application in Science and Engineering, August 7, 1980

Ityp = 3

フォーマット

定義

例（気体）

コメント