/EOS/IDEAL-GAS-VT

理想気体の状態方程式 $P = ρ r T$ を記述します。ここで、 $c_{p}$ は温度に依存し、次のようにJanafモデルを使用して記述されます： $c_{p} (T) = A_{0} + A_{1} T + A_{2} T^{2} + A_{3} T^{3} + A_{4} T^{4}$ (SI： $J . k g^{- 1} . K^{- 1}$ ）

フォーマット

(1)	(3)	(5)	(7)	(9)
/EOS/IDEAL-GAS-VT/`mat_ID`/`unit_ID`
`eos_title`
$r$	`T`₀	`P`₀	`P`_sh
`A`₀	`A`₁	`A`₂	`A`₃	`A`₄

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`eos_title`	状態方程式のタイトル。（文字、最大100文字）
$r$	理想気体定数。（実数）	$[J . k g^{- 1} . K^{- 1}]$
`T`₀	初期温度。`P`₀が定義されていない場合は、コメント3をご参照ください。（実数）	$[K]$
`P`₀	初期圧力`T`₀が定義されていない場合は、コメント3をご参照ください。（実数）	$[Pa]$
`P`_sh	圧力シフト（実数）	$[Pa]$
`A`₀	比熱容量パラメータ。（実数）	$[\frac{J}{k g . K}]$
`A`₁	比熱容量パラメータ。（実数）	$[\frac{J}{k g . K^{2}}]$
`A`₂	比熱容量パラメータ。（実数）	$[\frac{J}{k g . K^{3}}]$
`A`₃	比熱容量パラメータ。（実数）	$[\frac{J}{k g . K^{4}}]$
`A`₄	比熱容量パラメータ。（実数）	$[\frac{J}{k g . K^{5}}]$

例（空気）

#---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/LAW6/1/1
Air
#              RHO_I             
             1.22e-6     

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/EOS/IDEAL-GAS-VT/1/1
Air
#                  r                  T0                  P0                 PSH                
                 287                   0                 0.1                   0                               
#                 A0                  A1                  A2                  A3                  A4
              1004.5                   0                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

理想気体の状態方程式： $P V = R T$
ここで、

$V$

体積

$R$

一般理想気体定数（SI： $J . k g^{- 1} m o l^{- 1}$ ）

$T$

温度

これは、 $P v = r T$ または $P = ρ r T$ と記述することもできます。

ここで、

$v$

比容積（ $v = 1 / ρ$ ）

$r$

気体定数（SI： $J . k g^{- 1} K^{- 1}$ ）、 $r = R / M$ 。ここで、 $M$ は分子量
理想気体の代替形式は/EOS/IDEAL-GAS-VE（体積エネルギー）：(1)
$P (μ, E) = (γ - 1) (1 + μ) E$

ここで、

$γ = \frac{c_{p}}{c_{v}}$

定数

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

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

このケースでは、状態方程式/EOS/IDEAL-GAS-VT（体積温度）は $γ$ が定数でない場合に使用されるため、 $c_{p} (T)$ 関数が必要です。 $γ = \frac{c_{p}}{c_{v}} = \frac{c_{p}}{c_{p} - r}$ （ $c_{v} = c_{p} - r$ であるため）。

したがって、 $γ$ が定数の場合、 $A_{0} = c_{p} = \frac{γ r}{γ - 1}$ 、および $A_{1} = A_{2} = A_{3} = A_{4} = 0$ と定義すると、/EOS/IDEAL-GAS-VEの定式と同じ状態方程式が得られます。
T₀が設定されている場合、P₀は次のように計算されます：(2)
$P_{0} = r . ρ_{0} . T_{0}$

そうでない場合、P₀が設定されていれば、T₀は次のように計算されます：(3)
$T_{0} = P_{0} / r . ρ_{0}$
$c_{p} (T)$ 関数はJanafモデルを使用して次のように与えられます：(4)
$c_{p} (T) = A_{0} + A_{1} T + A_{2} T^{2} + A_{3} T^{3} + A_{4} T^{4}$

(SI： $J . k g^{- 1} . K^{- 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/IDEAL-GAS-VT

フォーマット

定義

例（空気）

コメント