Iform = 5

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

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

ユーザーは、 $υ$ =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}$ （ $E$ はオプション）を提供する必要があります。1つの線形EOSより： $P_{0} = C_{0} + C_{1} μ$ ；したがって、 $P_{s t a g n a t i o n} = C_{0}$ 。

各サイクルにおいて、Radiossは液体流入状態を、入力面における速度を用いてベルヌーイの定理が満足されるよう計算します。

フォーマット

(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`₁	`fct_ID`_{$ρ$ 1}	`fct_ID`_E1
$C_{1}^{m a t_1}$
$Δ P_{min}^{mat_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`₂	`fct_ID`_{$ρ$ 2}	`fct_ID`_E2
$C_{1}^{m a t_2}$
$Δ P_{min}^{mat_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`₃	`fct_ID`_{$ρ$ 3}	`fct_ID`_E3
$C_{1}^{m a t_3}$
$Δ P_{min}^{mat_3}$	$C_{0}^{m a t_3}$

定義

フィールド	内容	SI 単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子. （整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
`I`_form	定式化フラグ = 5 液体流入（停滞点でのデータから計算される）（整数）
`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}} \cdot f_{α_{i}} (t)$ （整数）
`fct_ID`_{$ρ$ i}	（オプション）密度比率スケーリング関数 $f_{ρ_{i}} (t)$ 識別子。 = 0 $P^{m a t_{i}} (t) = P_{0}^{m a t_{i}}$ > 0 $P^{m a t_{i}} (t) = P_{0}^{m a t_{i}} \cdot f_{P_{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}} \cdot f_{E_{i}} (t)$ （整数）
$C_{1}^{m a t_ i}$	理想気体EOSの係数 5 （実数）	$[Pa]$
$Δ P_{min}^{mat_i}$	流体力学的キャビテーション圧力 6 デフォルト = -10^-30（実数）	$[Pa]$
$C_{0}^{m a t_ i}$	理想気体EOSの係数 5 （実数）	$[Pa]$

停滞点から指定された気体状態 $ρ_{stagnation}, P_{stagnation}$ が、液体流入状態の計算に使用されます。Bernoulliの定理が適用されます：(1)
$P_{stagnation} = P_{in} + \frac{ρ_{in} v_{in}^{2}}{2}$

これが、流入境界要素内の副材料状態につながります：(2)
$ρ_{i n} = \frac{C_{1} \cdot ρ_{s t a g n a t i o n}}{C_{1} + \frac{ρ_{s t a g n a t i o n} v_{i n}^{2}}{2} (1 + C_{d})}$

ここで、 $C_{d}$ はオプションの降下パラメータ。(3)
$P_{i n} = P_{s t a g n a t i o n} - \frac{ρ_{s t a g n a t i o n} v_{i n}^{2}}{2} (1 + C_{d})$
(4)
${(ρ e)}_{i n} = (1 - \frac{ρ_{i n}}{ρ_{s t a g n a t i o n}}) P_{i n} + E_{s t a g n a t i o n}$

その後で、平均値を計算することにより、グローバル材料状態が決定されます。

圧力

$Δ 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_EXTでは、相対圧力 $Δ P_{min}^{m a t_i}$ を操作する場合に周囲圧力を考慮できます。このパラメータは、Radiossで各サイクルの正確なエネルギー統合に必要です。そうしないと、数値EOSの解が一般的に不正確になります。このパラメータは、全（物理）圧力を求めるためにEOS計算に含めなければならない圧力を表します。アニメーションファイル内の圧力コンターへの影響はありません。
線形EOSを使用した例：

全圧力： $P = P_{a m b} + C_{1} μ$ かつP_EXT = 0

相対圧力： $Δ P = C_{1} μ$ かつ $P_{E X T} = P_{a m b}$
体積比率によって、要素体積を3つの異なる材料で分け合うことができます。
材料ごとに、 $α_{0}^{mat_i}$ を0と1の間に定義する必要があります。

初期体積比率の合計 $\sum_{i = 1}^{3} α_{0}^{mat_i}$ は1に等しい必要があります。

体積の自動初期比率については、/INIVOLをご参照ください。
線形EOSは：(5)
$P (μ, E) = C_{0} + C_{1} μ$

これは、圧力とエネルギーを合計にするか相対にするかにより、柔軟性を提供します：(6)
$P (μ, E) = C_{0} + C_{1} μ$

ここで、 $C_{0} = P_{amp}, C_{1} = ρ_{0} c_{0}^{2}$ 、 $P_{E X T} = 0$

これにより、 $P (μ, E) = C_{0} + C_{1} μ$ から一般形式が得られます。(7)
$Δ P (μ, E) = C_{1} μ$

ここで、 $C_{1} = ρ_{0} c_{0}^{2}$ 、 $P_{E X T} = P_{a m b}$ .
$Δ P_{min}^{mat_i}$ フラグは、計算される圧力の最小値です。
$P = Δ P + P_{E X T}$ のため、 $P_{E X T} = 0$ の定義は $Δ P = P$ および $Δ P_{\min} = P_{\min}$ を意味します。

流体材料圧力を正のままにして、引張り強度を避ける必要があります。これにより、 $P_{m i n} = 0$ から $Δ P_{m i n} = - P_{E X T}$ が得られます。

ソリッド材料については、 $Δ P_{min}^{mat_i}$ = 10³⁰ のデフォルト値が適切です。
停滞エネルギーは省略可能です。これは、このEOSがエネルギーに依存しないためです。エネルギー入力値は、アニメーションファイルと時刻歴内の出力値にのみ影響します。
EOSパラメータは、隣接するMM-ALE領域からの流体EOSと一致する必要があります。

Iform = 5

フォーマット

定義

コメント