Iform = 6

ブロックフォーマットキーワードこの境界材料は、多相材料則/MAT/LAW51（I_form = 0、1、10または11）についてサイレント流出境界のシミュレーションを可能にします。本材料は、無限領域のモデル化に有益です。この材料則を割り当てられた要素はサイレント要素境界のように機能します。

これは、I_form = 3と使用される定式化を改善したもので、波形の反射と逆流の問題を最小にするのに役立ちます。計算領域の外側限界上の波形反射を回避するための音響処理および体積比率用の“連続”モデルの両方が含まれます。

すべての入力を空白のままとすることが可能で、Radiossは隣接要素の値に基づいてデフォルト値を計算します。したがって、モデルのすべての遠距離場境界について1つのパートしか使用できません。

フォーマット

(1)	(4)	(6)
/MAT/LAW51/`mat_ID`/`unit_ID`
`mat_title`
空白行
`I`_form
`P`_ext	`T`_cp	`T`_ca

#材料1パラメータ

(1)	(3)	(5)	(7)	(9)
$α_{0}^{mat_1}$	$ρ_{0}^{mat_1}$	$E_{0}^{m a t_1}$	$P_{\min}^{m a t_1}$	$P_{0}^{m a t_1}$
$s s p_{0}^{m a t_1}$
空白行

#材料2パラメータ

(1)	(3)	(5)	(7)	(9)
$α_{0}^{mat_2}$	$ρ_{0}^{mat_2}$	$E_{0}^{m a t_2}$	$P_{\min}^{m a t_2}$	$P_{0}^{m a t_1}$
$s s p_{0}^{m a t_2}$
空白行

#材料3パラメータ

(1)	(3)	(5)	(7)	(9)
$α_{0}^{mat_3}$	$ρ_{0}^{mat_3}$	$E_{0}^{m a t_3}$	$P_{\min}^{m a t_3}$	$P_{0}^{m a t_1}$
$s s p_{0}^{m a t_3}$
空白行

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
`I`_form	=6 定式化フラグ（整数）
`P`_ext	外圧デフォルトはStarterによって計算されます（実数）	$[Pa]$
`T`_cp	圧力緩和のための特性時間デフォルトはStarterによって計算されます（実数）	$[s]$
`T`_ca	体積比率緩和のための特性時間デフォルトはStarterによって計算されます（実数）	$[s]$
$α_{0}^{mat_j}$	副材料#j=1...3の初期体積比率デフォルトはStarterによって計算されます（実数）
$ρ_{0}^{mat_j}$	副材料#j=1...3の初期密度デフォルトはStarterによって計算されます（実数）	$[\frac{kg}{m^{3}}]$
$E_{0}^{m a t_j}$	副材料#j=1...3の初期体積エネルギーデフォルトはStarterによって計算されます（実数）	$[\frac{J}{m^{3}}]$
$P_{\min}^{m a t_j}$	副材料#j=1...3の圧力カットオフデフォルトはStarterによって計算されます（実数）	$[Pa]$
$P_{0}^{m a t_j}$	副材料#j=1...3の初期圧力デフォルトはStarterによって計算されます（実数）	$[Pa]$
$s s p_{0}^{m a t_j}$	副材料#j=1...3の初期音速デフォルトはStarterによって計算されます（実数）	$[\frac{m}{s}]$

例

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW51/1
Boundary Element, AUTO-OUTLET                                                                                  

#    IFORM
         6
#---Global parameters------------------------------------------------------------------------------#
#               PEXT                 TCP            TC_ALPHA
                   0                   0                   0               
#---Material#1:------------------------------------------------------------------------------------#
#            ALPHA_1             RHO_0_1               E_0_1             P_MIN_1               P_0_1
                   0                   0                   0                   0                   0
#              SSP_1
                   0

#---Material#2:------------------------------------------------------------------------------------#
#            ALPHA_2             RHO_0_2               E_0_2             P_MIN_2               P_0_2
                   0                   0                   0                   0                   0
#              SSP_2
                   0

#---Material#3:------------------------------------------------------------------------------------#
#            ALPHA_3             RHO_0_3               E_0_3             P_MIN_3               P_0_3
                   0                   0                   0                   0                   0
#              SSP_3
                   0

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

この圧力定式化は、Bayliss and Turkelによる方程式から導出されています。 8 (1)
$\frac{\partial P}{\partial t} = ρ c (\frac{\partial V_{n}}{\partial t} - V_{n} \cdot div (V - V_{n} n)) + \frac{c}{T_{c p}} (P_{\infty} - P)$

ここで、

$ρ c$

音響インピーダンス

$P_{\infty}$

流束が領域の内側からである、すなわち法線速度 $V_{n}$ <0（停滞圧力が代わりにかけられる）でない限り入力として与えられる遠距離場圧力。
材料の移流は次のように定義されます：
- 流出流束： $α_{i} = α_{i}^{n e i g h b o r}$
局所変数は単に更新されます。
- 流入流束： $\frac{δ α_{i}}{δ t} = - \frac{| α_{i}^{0} - α_{i} |}{T_{c α}}$
これは領域内で起こっていることが遠距離場でも起こると仮定しているため、局所変数は上述のケースのようには更新されません。しかしながら、流れが逆流する以前にそうであった状態と材料が局所的に類似していると考えるのは理にかなっています。緩和係数が追加され、流入材料が入力で与えられている設定された値へと次第に戻っていくことが可能となります。
定義されていない場合、全体パラメータP_extと初期状態は、隣接要素から自動的に計算されます。
爆発副材料については、初期状態およびデータは常に、隣接要素からのデータを用いて自動的に計算されます。
各特性時間は、1桁の、圧力波（T_cp）または材料（T_ca）が領域の最も近い限界に達するために必要な最短時間として計算されます。非圧縮性の場合、T_cpとT_caは等しいとすることが可能です。爆轟または衝撃波のシミュレーションの場合、T_cpは1桁の波の伝播速度を用いて計算されます。
各流出要素は、ALE領域からの1つの要素と1つの面を共有することが推奨されます。

図 2.
流出境界材料は、境界に達する波は局所的に平面波または球状波に近似されるべきことを意味する音響単極でできています。したがって、遠距離場の減少を正しくモデル化するには、流出要素は、爆発または衝突領域から、爆発または衝突領域の1倍から10倍の大きさ分のところに位置している必要があります。
Alvin Bayliss and Eli Turkel, Outflow boundary condition for fluid dynamics, Report 80-21, Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, Virginia

Iform = 6

フォーマット

定義

例

コメント