タンクテスト

タンクテスト出力からのデータを使用し、Radioss入力で使用可能な、供給されるガスの温度とマスフローを得ることができます。

タンクテストで、インジェクションの点とタンク中央の点で圧力の計測が可能で、それらが等しければ、圧力変化を得ることができます。また、共有されるガスの量と、ガスの特性も得られます。

温度の場合、温度計が十分に正確でないことが多いため、正確なテストデータを取得することはより困難です。したがって、タンクテストの温度が正確でない可能性があります。

以下は、注入温度とタンク内の温度が既知であるかどうかを考慮に入れています。

インジェクターとタンクの温度が未知のケースでは：

初期と注入ガス組成
それぞれの要素のモル質量（ $M W_{i}$ ）とモル分率（ $X_{i}$ ）が既知であればガスのモル質量（ $M W$ ）を定義することが可能です：(1)
$M W = \sum X_{i} \cdot M W_{i}$

ガスの混合体の単位体積当たりの平均熱容量はAmagat-Leduc方程式で与えられます：(2)
$C_{p} (T) = \frac{\sum m_{i} \cdot C_{p a}}{\sum m_{i}} + \frac{\sum m_{i} \cdot C_{p b}}{\sum m_{i}} T + \frac{\sum m_{i} \cdot C_{p c}}{\sum m_{i}} T^{2}$

前の式を用いて、初期状態と注入された混合体の熱容量係数（ $C_{p} (T)$ ）の定義が可能な場合があります。

注入されたガスの特性が既知であれば、初期のガスと混合体のマスフローとインジェクターの温度を見つける事が可能です。以下の基礎方程式が解析の実行に用いられます。

理想気体の状態方程式は：(3)
$P V = n R T, n = \frac{m}{M W}$
ここで、 $R = 8.314 \frac{J}{m o l e \cdot K}$
断熱の方程式は：(4)
$H = c o n s t .$
ここで、 $H$ は系の全エンタルピー（インフレーター＋タンク）です。

エネルギー保存から、タンクテストの基礎のエネルギー方程式は次のように書くことができます：(5)

d E_{a i r b a g} = - P d V + d H_{i n} - d H_{o u t}

ここで、タンクテストは断熱のため、 $d H_{o u t}$ =0。タンクテストが一定体積であるため、 $d V$ =0。

したがって、式 5は下記のように要約することができます：(6)

U_{f} - U_{0} = d H_{i n}

$\Leftrightarrow \int_{0}^{T_{m i x}} m_{(i n + a i r)} \cdot C_{V (i n + a i r)} d T - \int_{0}^{T_{a i r}} m_{(a i r)} \cdot C_{V (a i r)} d T = \int_{0}^{T_{i n}} m_{(i n)} \cdot C_{P (i n)} d T$

$\begin{array}{l} \Leftrightarrow m_{(m i x)} T_{m i x} (C_{p a (m i x)} + C_{p b (m i x)} \frac{T_{m i x}}{2} + C_{p c (m i x)} \frac{T_{m i x}^{2}}{3} - \frac{R}{M W_{m i x}}) - m_{(a i r)} T_{0} \\ (C_{p a (a i r)} + C_{p b (a i r)} \frac{T_{0}}{2} + C_{p c (a i r)} \frac{T_{0}^{2}}{3} - \frac{R}{M W_{a i r}}) = m_{(i n)} T_{i n} (C_{p a (i n)} + C_{p b (i n)} \frac{T_{i n}}{2} + C_{p c (i n)} \frac{T_{i n}^{2}}{3}) \end{array}$

インフレーター温度

式 6で、未知の変数は $T_{i n i}$ のみです。

他の変数は既知であるか、または式 2と式 3を用いて決めることができます：

$m_{(i n)} m_{(a i r)} = M W_{a i r} \frac{P_{0} V}{R T_{0}}$ および $n_{(m i x)} = \frac{m_{(i n)}}{M W_{i n}} + \frac{m_{(a i r)}}{M W_{a i r}}$
i =aからcの場合： $C_{p i}_{(i n)}, C_{p i}_{(a i r)}, C_{p i}_{(m i x)}$ は式 2で計算され、
$T_{m i x} = \frac{P_{t a n k} V_{t a n k}}{n_{(m i x)} R}$

したがって、式 6で、注入されるガスのインジェクタ $T_{i n}$ における温度は、 $T_{i n}$ の反復で見つけることができます。最初に温度を推定し、収束解を得るのに６回の反復もあれば十分です。

マスフロー

タンクテストの上部での圧力の時間変化が既知であれば、マスフロー比率は次の式で決める事ができます：(7)

\dot{m} = \frac{\partial m}{\partial P} \frac{\partial P}{\partial t} ≅ \frac{Δ M}{Δ P} \frac{\partial P}{\partial t}

ここで、

$Δ P$: 実験の間の全圧力変化
$Δ M$: 全注入質量

式 7 は、質量変化と圧力変化の関係が完全に増加関数の場合に成り立ちますが、これはその場合です。

インフレーターガス速度

圧力は速やかに一定になるため、以下の式で表すことができるかもしれません：

P_{i n} = P (t)

、

T_{i n}

が既知であれば、密度は

T_{i n}

および

P (t)

の関数として表すことができます：(8)

ρ_{i n} = \frac{P (t)}{\frac{T_{i n} R}{M W_{i n}}}

V_{i n} (t) = \frac{\frac{d m (t)}{d t}}{S ρ_{i n} (t)}

の他に、ガス内の音速は：(9)

c_{i n}^{2} = γ \frac{R}{M W_{i n}} T_{i n}

また、

V_{i n} > c_{i n}

の場合、

V_{i n} = c_{i n}

となり、速度は音速となります。

インジェクタの温度が既知の場合：
$\frac{d m}{d t} = \frac{v}{γ (γ - 1) C_{V} T_{i n}} \frac{d P}{d t}$
タンク内の温度が既知の場合：
$\frac{d m}{d t} = \frac{V}{γ (γ - 1) C_{V} T^{2}} (T \frac{d P}{d t} - P \frac{d T}{d t})$ または $m = \frac{P V}{(γ - 1) C_{V} T}$

$T_{i n} = \frac{T^{2}}{γ} (\frac{\frac{d P}{d t}}{T \frac{d P}{d t} - P \frac{d T}{d t}}) or T_{i n} = \frac{T}{γ} + \frac{m}{γ} \frac{\frac{d T}{d t}}{\frac{d m}{d t}}$
タンク内の温度が一定の場合：
$T_{i n} = \frac{T}{γ}$