Iform = 4

Blockフォーマットキーワードこの境界を使用すれば、多相材料ALE則（定式化： I_form = 0、1、10または11）の気体流入条件をシミュレートできます。

境界副材料状態は、ユーザーによって指定された停滞点での状態から計算されます。この機能の使用時は、速度が数値的なスキームによって計算される強制速度（/IMPVEL）を使用する必要はありません。

内容

ユーザーは、v=0の状態に相対する停滞点の状態 $α_{s t a g n a t i o n} = α_{0}$ 、 $ρ_{s t a g n a t i o n} = ρ_{0}$ および $E_{s t a g n a t i o n} = E_{0}$ を提供する必要があります。理想気体EOSより：

P_{0} = C_{0} + (1 + μ) \cdot (γ - 1) \cdot E_{0}

ここで、 $C_{4} = γ - 1$ 。ここから、 $P_{s t a g n a t i o n} = C_{0} + C_{4} \cdot E_{0}$ が導かれます。

各サイクルにおいて、Radiossは気体流入状態

ρ_{i n}, E_{i n}, P_{i n}

を、入力面における速度を用いてベルヌーイの定理が満足されるよう1計算します。

フォーマット


(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/MAT/LAW51/`mat_ID`/`unit_ID`
`mat_title`
空白のフォーマット
`I`_form

#グローバルパラメータ


(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`Scale`_time		`P`_EXT

#材料1パラメータ


(1)	(3)	(5)	(7)	(8)	(9)
$α_{0}^{mat_1}$	$ρ_{0}^{mat_1}$	$E_{0}^{m a t_ 1}$	`fct_ID`_{$α$ 1}	`fct_ID`_{$ρ$ 1}	`fct_ID`_E1
$C_{1}^{m a t_1}$			$C_{4}^{m a t_1}$
	$C_{0}^{m a t_1}$

#材料2パラメータ


(1)	(3)	(5)	(7)	(8)	(9)
$α_{0}^{mat_2}$	$ρ_{0}^{mat_2}$	$E_{0}^{m a t_ 2}$	`fct_ID`_{$α$ 2}	`fct_ID`_{$ρ$ 2}	`fct_ID`_E2
$C_{1}^{m a t_2}$			$C_{4}^{m a t_2}$
	$C_{0}^{m a t_2}$

#材料3パラメータ


(1)	(3)	(5)	(7)	(8)	(9)
$α_{0}^{mat_3}$	$ρ_{0}^{mat_3}$	$E_{0}^{m a t_ 3}$	`fct_ID`_{$α$ 3}	`fct_ID`_{$ρ$ 3}	`fct_ID`_E3
$C_{1}^{m a t_3}$			$C_{4}^{m a t_3}$
	$C_{0}^{m a t_3}$

定義


フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	Unit Identifier。（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
`I`_form	定式化フラグ = 4 気体流入（停滞点でのデータから計算される）（整数）
`Scale`_time	入力の関数の横軸のスケールファクタ 2 デフォルト = 1（実数）
`P`_EXT	外部（周囲）圧力 3 （実数）	$[Pa]$
$α_{0}^{mat_i}$	初期体積比率 4 （実数）
$ρ_{0}^{mat_i}$	停滞点での初期密度 1 （実数）	$[\frac{kg}{m^{3}}]$
$E_{0}^{m a t_ i}$	停滞点での初期エネルギー 5 （実数）	$[\frac{J}{m^{3}}]$
`fct_ID`_{$α$ i}	（オプション）体積比率スケーリング関数 $f_{α_{i}} (t)$ 識別子 6 = 0 $α^{m a t_{i}} (t) = α_{0}^{m a t_{i}}$ > 0 $α^{m a t_{i}} (t) = α_{0}^{m a t_{i}} f_{α_{i}} (t)$ （整数）
`fct_ID`_{$ρ$ i}	（オプション）密度比率スケーリング関数 $f_{ρ_{i}} (t)$ 識別子 = 0 $ρ^{m a t_{i}} (t) = ρ_{0}^{m a t_{i}}$ > 0 $ρ^{m a t_{i}} (t) = ρ_{0}^{m a t_{i}} . f_{ρ_{i}} (t)$ （整数）
`fct_ID`_Ei	（オプション）エネルギー比率スケーリング関数 $f_{E_{i}} (t)$ 識別子 = 0 $E^{m a t_{i}} (t) = E_{0}^{m a t_{i}}$ > 0 $E^{m a t_{i}} (t) = E_{0}^{m a t_{i}} f_{E_{i}} (t)$ （整数）
$C_{0}^{m a t_ 3}$	理想気体EOSの係数 5 （実数）	$[Pa]$
$C_{4}^{m a t_ i}$	理想気体（ $γ - 1$ ）定数 5 （実数）
$C_{0}^{m a t_ i}$	理想気体EOSの係数 5 （実数）	$[Pa]$

与えられた停滞点 $ρ_{stagnation}, P_{stagnation}$ が、気体流入状態の計算に用いられます。ベルヌーイの定理が適用されます。

$P_{stagnation} = P_{in} + \frac{ρ_{in} v_{in}^{2}}{2}$

これが次の流入状態につながります。
$ρ_{in} = ρ_{stagnation} {[1 - \frac{γ - 1}{2 γ} \cdot \frac{ρ_{stagnation}}{P_{stagnation}} \cdot (1 + C_{d}) \cdot v_{in}^{2}]}^{\frac{1}{γ - 1}}$
$P_{i n} = P_{s t a g n a t i o n} {(\frac{ρ_{i n}}{ρ_{s t a g n a t i o n}})}^{γ}$
${(ρ e)}_{i n} = \frac{P_{a}}{γ - 1} {(\frac{ρ_{i n}}{ρ_{stagnation}})}^{γ - 1}$
その後で、平均値を計算することにより、グローバル材料状態が決定されます。

圧力

$Δ P_{i n} = \sum_{i} α^{m a t_{i}} (t) Δ P_{i n}^{m a t_i}$

密度

$ρ_{i n} = \sum_{i} α^{m a t_{i}} (t) ρ_{i n}^{m a t_i}$

エネルギー

${(ρ e)}_{i n} = \sum_{i} α^{m a t_{i}} (t) E_{i n}^{m a t_i}$
オプションの比率関数は、体積、密度またはエネルギー比率をスケーリングするために使用できます。
パラメータ $P_{E X T}$ では、相対圧力 $Δ P_{min}^{m a t_i}$ を操作する場合に周囲圧力を考慮できます。このパラメータは、Radiossで各サイクルの正確なエネルギー統合に必要です。そうしないと、数値EOSの解が一般的に不正確になります。このパラメータは、全（物理）圧力を求めるためにEOS計算に含めなければならない圧力を表します。アニメーションファイル内の圧力コンターへの影響はありません。
線形EOSを使用した例：
全圧力： $P = P_{amb} + C_{1} μ$ また、 $P_{E X T} = 0$
相対圧力： $Δ P = C_{1} μ$ 、また $P_{E X T} = P_{a m b}$
体積比率によって、要素体積を3つの異なる材料で分け合うことができます。
材料毎に、 $α_{0}^{m a t_ i}$ を0と1の間に定義する必要があります。
初期体積比率の合計 $\sum_{i = 1}^{3} α_{0}^{m a t _ i}$ は1に等しい必要があります。
体積の自動初期比率については、/INIVOLをご参照ください。
理想気体EOSは $P (μ, E) = (γ - 1) (1 + μ)$ です。通常は、この一般形式 $P = C_{0} + C_{1} μ + C_{4} (1 + μ) E$ を使用して記述できます。ここで、 $C_{4} = (γ - 1)$ です。これは、圧力とエネルギーを合計にするか相対にするかにより、柔軟性が向上します：
$P (μ, E) = C_{4} (1 + μ) E$
ここで、 $C_{4} = (γ - 1)$ 、 $P_{E X T} = 0$
これにより、 $Δ P (μ, E) = C_{0} + C_{4} (1 + μ) E$ から一般形式が得られます。
$Δ P (μ, E) = C_{0} + C_{4} (1 + μ) E$
ここで、 $C_{4} = (γ - 1)$ 、 $C_{0} = - P_{0}$ および $P_{E X T} = P_{a m b}$
$Δ P (μ, Δ E) = C_{0} + C_{1} μ + C_{4} (1 + μ) Δ E$
ここで、 $C_{4} = (γ - 1)$ 、 $C_{1} = E_{0} (γ - 1)$ および $P_{E X T} = P_{a m b}$
$Δ P_{\min}^{m a t_ i}$ フラグは、計算される圧力の最小値です。
$P = Δ P + P_{E X T}$ のため、 $P_{E X T} = 0$ の定義は $Δ P \equiv P$ および $Δ P_{m i n} \equiv P_{m i n}$ を意味します。
液体材料圧力を正のままにして、引張り強度を避ける必要があります。これにより、 $P_{m i n} = 0$ から $Δ P_{m i n} = - P_{E X T}$ が得られます。
ソリッド材料については、 $Δ P_{\min}^{m a t_ i} = 10^{30}$ のデフォルト値が適切です
EOSパラメータは、隣接するMM-ALE領域からの気体EOSと一致する必要があります。

Iform = 4

内容

フォーマット

定義

コメント