/EOS/NOBLE-ABEL

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

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/EOS/NOBLE-ABEL/mat_ID/unit_ID
eos_title
b γ E0 Psh ρ 0

定義

フィールド 内容 SI 単位の例
mat_ID 材料識別子

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
eos_title 状態方程式のタイトル

(文字、最大100文字)

 
b 補容積

(実数)

[ m 3 k g ]
γ 熱容量の比 γ = C p C v

(実数)

 
E0 単位参照体積あたりの初期エネルギー

(実数)

[ J m 3 ]
Psh 圧力シフト

(実数)

[ Pa ]
ρ 0 参照密度

デフォルト = 材料密度(実数)

[ 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

コメント

  1. 理想気体の熱状態方程式 P v = R T の一般化は、右記の補容積状態方程式です: P ( v b ) = R T
    ここで、
    v
    比容積
    b
    補容積
    R
    比例気体定数
    T
    温度

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

    ここで、
    µ = ρ ρ 0 1
    E = E i n t V 0
    P ( μ , E ) = ( γ 1 ) ( 1 + μ ) E 1 b ρ 0 ( 1 + μ )

    ここで、 γ = C p C v

  2. この状態方程式は、分子自体で占められる体積が無視できなくなっている高圧における高密度気体に適用されます。
  3. 補容積bは通常、右記のの範囲にあります; [ 0.9 × 10 3 , 1.1 × 10 3 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaqaaaaaaaaa Wdbmaadmaapaqaa8qacaaIWaGaaiOlaiaaiMdacqGHxdaTcaaIXaGa aGima8aadaahaaWcbeqaa8qacqGHsislcaaIZaaaaOGaaiilaiaaig dacaGGUaGaaGymaiabgEna0kaaigdacaaIWaWaaWbaaSqabeaacqGH sislcaaIZaaaaaGccaGLBbGaayzxaaaaaa@4890@ [ m 3 k g ]
  4. 理想気体状態方程式との比較
      理想ガス NOBLE-ABEL
    P ( v , T ) P v = R T P ( v b ) = R T
    P ( μ , E ) ( γ 1 ) ( 1 + μ ) E ( γ 1 ) ( 1 + μ ) E 1 b ρ 0 ( 1 + μ )
    音速

    c

    c = γ P ρ c = γ P ( 1 b ρ ) ρ
    E 0 = E ( 0 ) P 0 γ 1 P 0 ( 1 b ρ 0 ) γ 1
  5. 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)