# Tank Test

By using data from a tank test output, it is possible to obtain the temperature and the mass flow of the gas supplied, which can be used as input to Radioss.

With a tank test it is possible to measure the pressure at the injection point or in the middle of the tank, the two values are equal so the pressure variation is well known. Also, the amount of gas supplied and the characteristics of the gas in the test are known.

The following cases take into account that you may or may not know the injected temperature and the temperature in the tank.

- Initial and injected gas composition
Knowing the molecular weight ( $M{W}_{i}$ ) of each elements and the molar fraction ( ${X}_{i}$ ), it is possible to define the molar weight of the gas ( $MW$ ):

$$MW={\displaystyle \sum {X}_{i}\cdot M{W}_{i}}$$The average heat capacity, per mass unit of a mixture of gases is given by the Amagat-Leduc equation:

$${C}_{p}\left(T\right)=\frac{{\displaystyle \sum {m}_{i}\cdot {C}_{pa}}}{{\displaystyle \sum {m}_{i}}}+\frac{{\displaystyle \sum {m}_{i}\cdot {C}_{pb}}}{{\displaystyle \sum {m}_{i}}}T+\frac{{\displaystyle \sum {m}_{i}\cdot {C}_{pc}}}{{\displaystyle \sum {m}_{i}}}{T}^{2}$$

Using the previous equation, the heat capacity coefficients ( ${C}_{p}(T)$ ) of the initial and injected mixture may be defined.

- The perfect gas equation of state is:$$PV=nRT,\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=\frac{m}{MW}$$with $R=8.314\frac{J}{mole\cdot K}$
- And the adiabatic equation:$$H=const.$$with $H$ being the total enthalpy of the system (inflator + tank).

From the conservation of the energy, the basic energy equation of the tank test can be written as:

Here $d{H}_{out}$ =0, since the tank test is adiabatic. The constant volume of the tank test means $dV$ =0.

Therefore, Equation 5 is summarized as:

$\iff {\displaystyle \underset{0}{\overset{{T}_{mix}}{\int}}{m}_{\left(in+air\right)}}\cdot {C}_{V\left(in+air\right)}dT-{\displaystyle \underset{0}{\overset{{T}_{air}}{\int}}{m}_{\left(air\right)}}\cdot {C}_{V\left(air\right)}dT={\displaystyle \underset{0}{\overset{{T}_{in}}{\int}}{m}_{\left(in\right)}}\cdot {C}_{P\left(in\right)}dT$

$\begin{array}{l}\iff {m}_{\left(mix\right)}{T}_{mix}\left({C}_{pa\left(mix\right)}+{C}_{pb\left(mix\right)}\frac{{T}_{mix}}{2}+{C}_{pc\left(mix\right)}\frac{{T}_{mix}{}^{2}}{3}-\frac{R}{M{W}_{mix}}\right)-{m}_{\left(air\right)}{T}_{0}\\ \left({C}_{pa\left(air\right)}+{C}_{pb\left(air\right)}\frac{{T}_{0}}{2}+{C}_{pc\left(air\right)}\frac{{T}_{0}{}^{2}}{3}-\frac{R}{M{W}_{air}}\right)={m}_{\left(in\right)}{T}_{in}\left({C}_{pa\left(in\right)}+{C}_{pb\left(in\right)}\frac{{T}_{in}}{2}+{C}_{pc\left(in\right)}\frac{{T}_{in}{}^{2}}{3}\right)\end{array}$

## Inflator Temperature

In Equation 6, the unknown variable is only ${T}_{ini}$ .

- ${m}_{\left(in\right)}{m}_{\left(air\right)}=M{W}_{air}\frac{{P}_{0}V}{R{T}_{0}}$ and ${n}_{\left(mix\right)}=\frac{{m}_{\left(in\right)}}{M{W}_{in}}+\frac{{m}_{\left(air\right)}}{M{W}_{air}}$
- for i = a to c: ${C}_{pi}{}_{\left(in\right)},{C}_{pi}{}_{\left(air\right)},{C}_{pi}{}_{\left(mix\right)}$ are calculated with Equation 2, and
- ${T}_{mix}=\frac{{P}_{tank}{V}_{tank}}{{n}_{(mix)}R}$

Therefore, Equation 6 finds the temperature at the injector ${T}_{in}$ of the injected gas by iterating on ${T}_{in}$ . First, the temperature is guessed and six iterations are sufficient to converge to the solution.

## Mass Flow

Knowing the evolution of the pressure versus time at the top of the tank test, it is possible to determine the mass flow rate with:

- $\text{\Delta}P$
- Total pressure variation during the experiment
- $\text{\Delta}M$
- Total injected mass

## Inflator Gas Velocity

Since the pressure is quickly uniform, the following equation may be written: ${P}_{in}=P(t)$ , knowing ${T}_{in}$ the density may be expressed as a function of ${T}_{in}$ and $P(t)$ :

Besides, ${V}_{in}(t)=\frac{\raisebox{1ex}{$dm(t)$}\!\left/ \!\raisebox{-1ex}{$dt$}\right.}{S{\rho}_{in}(t)}$ , the sound speed into the gas is equal to:

- If the injector temperature is known:
$\frac{dm}{dt}=\frac{v}{\gamma \left(\gamma -1\right){C}_{V}{T}_{in}}\frac{dP}{dt}$

- If the temperature in the tank is known:
$\frac{dm}{dt}=\frac{V}{\gamma \left(\gamma -1\right){C}_{V}{T}^{2}}\left(T\frac{dP}{dt}-P\frac{dT}{dt}\right)$ or $m=\frac{PV}{\left(\gamma -1\right){C}_{V}T}$

${T}_{in}=\frac{{T}^{2}}{\gamma}\left(\frac{\frac{dP}{dt}}{T\frac{dP}{dt}-P\frac{dT}{dt}}\right)\text{\hspace{1em}}\text{or}\text{\hspace{1em}}{T}_{in}=\frac{T}{\gamma}+\frac{m}{\gamma}\frac{\frac{dT}{dt}}{\frac{dm}{dt}}$

- If the temperature in the tank is constant:
${T}_{in}=\frac{T}{\gamma}$