/EOS/NOBLE-ABEL

ブロックフォーマットのキーワード補容積状態方程式 $P (v - b) = R T$ を記述します。

フォーマット

(1)	(3)	(5)	(7)	(9)
/EOS/NOBLE-ABEL/`mat_ID`/`unit_ID`
`eos_title`
`b`	$γ$	`E`₀	`P`_sh	$ρ_{0}$

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`eos_title`	EOSのタイトル（文字、最大100文字）
`b`	補容積（実数）	$[\frac{m^{3}}{k g}]$
$γ$	熱容量の比 $γ = \frac{C_{p}}{C_{v}}$ （実数）
`E`₀	単位参照体積あたりの初期エネルギー（実数）	$[\frac{J}{m^{3}}]$
`P`_sh	圧力シフト（実数）	$[Pa]$
$ρ_{0}$	参照密度デフォルト = 材料密度（実数）	$[\frac{k g}{m^{3}}]$

例

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/HYDPLA/7
NOBLE_ABEL
#              RHO_I               RHO_0
             1.22e-6                   0
#                  E                  nu
                   0                   0
#                  a                   b                   n             eps_max           sigma_max
                1E30                   0                   0                   0                   0
#               Pmin                 Psh
                   0
/EOS/NOBLE-ABEL/7
NOBLE-ABEL EOS 
#                  b               GAMMA                  E0                 PSH                RHO0          
                1E-3                 1.4        0.2499999997                   0             1.22e-6
/ALE/MAT/7

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

理想気体の熱状態方程式 $P v = R T$ の一般化は、右記の補容積状態方程式です： $P (v - b) = R T$
ここで、

$v$

比容積

b

補容積

R

比例気体定数

T

温度

以前のフォーム $P = P (v, T)$ は、 $P = P (μ, E)$ フォームで書き出すことができます。

ここで、

$µ = \frac{ρ}{ρ_{0}} - 1$

$E = \frac{E_{i n t}}{V_{0}}$

$P (μ, E) = \frac{(γ - 1) (1 + μ) E}{1 - b ρ_{0} (1 + μ)}$

ここで、 $γ = \frac{C_{p}}{C_{v}}$
この状態方程式は、分子自体で占められる体積が無視できなくなっている高圧における高密度気体に適用されます。
補容積bは通常、右記のの範囲にあります； $[0.9 \times 10^{- 3}, 1.1 \times 10^{- 3}]$ $[\frac{m^{3}}{k g}]$

理想気体状態方程式との比較

	理想ガス	NOBLE-ABEL
$P (v, T)$	$P v = R T$	$P (v - b) = R T$
$P (μ, E)$	$(γ - 1) (1 + μ) E$	$\frac{(γ - 1) (1 + μ) E}{1 - b ρ_{0} (1 + μ)}$
音速 $c$	$c = \sqrt{\frac{γ P}{ρ}}$	$c = \sqrt{\frac{γ P}{(1 - b ρ) ρ}}$
$E_{0} = E (0)$	$\frac{P_{0}}{γ - 1}$	$\frac{P_{0} (1 - b ρ_{0})}{γ - 1}$

Radiossにより流体力学的圧力の計算に用いられ、右記の材料則と適合性のある状態方程式。
- /MAT/LAW3 (HYDPLA)
- /MAT/LAW4 (HYD_JCOOK)
- /MAT/LAW6 (HYDROまたはHYD_VISC)
- /MAT/LAW10 (DPRAG1)
- /MAT/LAW12 (3D_COMP)
- /MAT/LAW49 (STEINB)
- /MAT/LAW102 (DPRAG2)
- /MAT/LAW103 (HENSEL-SPITTEL)

/EOS/NOBLE-ABEL

フォーマット

定義

例

コメント