Issue 
Mechanics & Industry
Volume 20, Number 7, 2019
Application of Experimental and Numerical Methods in Fluid Mechanics and Energy



Article Number  706  
Number of page(s)  9  
DOI  https://doi.org/10.1051/meca/2019079  
Published online  12 December 2019 
Regular Article
Analysis of room overpressure after launching an integrated extinguishing system
^{1}
University of West Bohemia in Pilsen, New Technologies − Research Centre, Univerzitní 2732/8, 301 00 Pilsen, Czech Republic
^{2}
Safran Cabin CZ s.r.o., Univerzitní 1119/34, 301 00 Pilsen, Czech Republic
^{*} email: hurda@ntc.zcu.cz
Received:
31
March
2019
Accepted:
3
November
2019
The throttling of liquefied substances from high pressure vessels to an enclosed volume starting at atmospheric pressure is described in order to determine the thermodynamic state of gaseous contents of the extinguished room. A time dependent, 0D mathematical model is implemented describing the state inside the agent container, the isenthalpic throttling in the distribution system, agent vaporization and mixing with air. The agent is modelled as real gas, and solid–vapour two phase behaviour is approximated. Other influences on the process including the heat transfer from selected solid parts inside the room and the gas mixture leakage out of the room are taken into account. The main outcome is an MS Excel tool for the design optimization of integrated fire extinguishers. The optimization balances two contradictory requirements: volumetric concentration of the agent to sustain the fire extinguishing capabilities, and tolerable room overpressure. Agent fill weight and discharge time are adjusted. The discharge time is controlled by the distribution piping and design of the spray nozzles. System operation is checked with regards to various initial and boundary conditions.
Key words: Fire suppression / gas outflow from vessel / triple point
© L. Hurda et al., published by EDP Sciences 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The design tool for an integrated fire extinguishing system is described. Many aspects of the process are covered by standards but concerns about the emerging overpressure in the extinguished volume require careful analysis and a proposal of new techniques.
This system and process was already briefly discussed in [1].
The thermodynamic system consists of two main parts: an enclosed volume (a room) and extinguishing agent container (a pressure vessel).
The room is initially filled with clean, dry air at atmospheric or nearatmospheric pressure and ordinary temperature. The heat transfer between the walls and the surroundings is omitted completely as the reallife object being modelled is well insulated and relatively low temperature gradients are expected. The model is not completely closed regarding mass conservation. The simplest model of gas mixture leakage based on constant dimensions of a virtual nozzle i is discussed and sophisticated leakage model based on measurement data is implemented.
The vessel may have a different initial temperature than the room. The heat transfer through the vessel wall is also neglected as the conditions for its evaluation are unclear. The surfacetovolume ratio of pressure vessels of common types is relatively low, so the heat transfer should not have a large influence on the relatively fast process, even though the vessel materials are usually highly conductive and relatively large temperature gradients may occur. The vessel is installed in a separate compartment of the room, which prevents sudden changes in environment temperature. Also, the heat being transferred in the whole process has similar effect and the same magnitude whether it is conducted through the vessel wall, the walls of the distribution system piping or exchanged in direct contact with sprayed droplets.
The room is installed with a piping system ending with nozzle sets for agent distribution. This parts' heat capacity and heat transfer properties are neglected for similar reasons as for the vessel. The whole distribution system is modelled only in terms of pressure drop.
Real gas thermodynamic data for three substances from the NIST REFPROP database are used for various calculations of agent state change: CO_{2} and two hydrocarbon/halocarbon mixtures. These mixtures are used for approximating a pure compound used in reality that is not available in the properties database. The substitutes were chosen based on their interchangeability in refrigeration applications as stated in [2]. The agent is expected to be liquid initially.
The proper functioning of the fire extinguishing system is ensured by a sufficient volumetric concentration of agent in the extinguished volume, it is checked throughout the timedependent calculation. The unusual limit of the overpressure comes from the demand for integrity of the lightweight structural design of the room walls. Such walls are also highly pliable and can deflect while introducing room volume expansion by work expelled by gas and enlargement of the slots which mediate the agentair mixture leakage out of the room. The latter effect is discussed and addressed.
2 Thermodynamic state of extinguishing agent in the pressure vessel
The main reason for being concerned with the process going on during the emptying of the bottle is to find the starting point of isenthalpic expansion (throttling) which occurs in the piping and mainly in the nozzles (see Sect. 3.1), as this enables the enthalpy of the released agent to be determined.
The vessel is designed to release only the liquid phase of the agent, even though the agent enters the twophase region just after the start of release. This is ensured by a siphon tube or a bottom release socket. The presence of relatively well separated liquid and vapour phases is expected, despite the ongoing boiling process in the vessel.
The state of the released agent is then found on the liquid saturation curve. How the agent pressure and mass drops with time also needs to be examined, then the state of agent in the bottle can be determined in discrete time levels.
2.1 Initial release of pressure in fluid
In the first moment after the container socket is opened, the pressure drops quickly as the contents of the container are mostly incompressible. In reality, it is necessary for the vessel to be filled with a small fraction of inert gas to compensate for the temperature changes of the vessel contents, but this fraction occupies too small volume to have a significant effect on the pressure release.
The first time step of the calculation finds the point on the liquid saturation curve which is reached by isentropic expansion of fluid from the prescribed initial pressure and temperature, see (1). The isentropic nature is expected as the kinetic aspects of the process are minor.$${p}_{\mathrm{sat}}=f\left({s}_{\mathrm{sat}}\right);{s}_{\mathrm{sat}}={s}_{\text{init}}=f\left({p}_{\mathrm{init}},{T}_{\mathrm{init}}\right)\text{.}$$(1)
2.2 Mass flow estimation
An important assumption is made that the mass flow rate is constant during most of the agent discharge process. This conforms with [3,4]. With this assumption, the value of the mass flow rate can become an input quantity for the design of the distribution system, so it is not calculated in any way, but may be a subject of optimization.
An explanation why the two phase, boiling mix discharge rate does not vary with dropping pressure in the container may lie in the subsonic design of the spray nozzles which control the flow speed and the twophase state of the agent in front of the nozzles allowing the compressible flow to occur.
The fluid definitely becomes a twophase mixture before reaching the nozzles. The partial evaporation is caused by a pressure drop in the parts of the distribution system upstream from the nozzle (probably a minor drop compared to nozzle) and the heat transfer through the piping walls delivering certain amount of latent heat to the flowing agent. It is not known what determines the constant mass flow rate and if the expansion in the nozzles is dominated by liquid, gas or complex twophase behaviour. This would require tedious research. For this purpose, the experimental verification (see Sect. 2.3) and similar experience of other authors is considered satisfactory.
In case of gas dominated outflow, the pressure ratio is high enough to choke the flow through the convergent nozzles for most of the two phase fluid release. The critical pressure ratio ranges between 0.55 and 0.58 for all subjected gases.
The end of the discharge is determined by the point where there is only vapour and no liquid in the vessel. The period of subsonic outflow could possibly start near this point and is neglected due to the low total mass of the outflow.
2.3 Mass flow experimental verification
The mass flow estimation was verified experimentally. In Figure 1, a constant mass flow rate may be observed for most of the total mass of the agent discharged. The diminishing outflow on the right side of Figure 1 is the part of discharge after the “end of liquid” point. This point is approached extremely early in the CO_{2} test due to reaching the states under the triple point, where the vessel siphon tube freezes fully. This later discharge is neglected in the calculation and the discharge is considered to continue at the prescribed constant flow rate.
The linear parts of the curves from Figure 1 do confirm the mass flow rate is constant, however, for a given distribution system, an equal value of this constant should be observed. This discrepancy lies in the setup of the experiment (see Fig. 2), in which one to three pressure vessels of gas were used. These were standard industrial pressure containers equipped with valves, which are not designed to deplete the content in a matter of seconds. Therefore, the limiting sonic flow occurred first in the parallelmounted vessel valves. This hypothesis is confirmed by the fact that the mass flow rates for the three cases are approximately in the same ratio as the total initial mass and the number of vessels − 1:2:3.
Fig. 1 Discharge process for three experimental tests on a room prototype with CO_{2} used for verification of the constant mass flow rate. Thick grey lines indicate the linear approximation of mass, i.e. approximation of constant mass flow rate. 
Fig. 2 Simplified experiment layout. Standard integrated pressure vessel (grey) replaced with external pressure bottles of CO_{2} for parallel discharge controlled initiated by manual release valve (green). 
2.4 Model of state of the agent in the vessel during outflow
The state change of the agent in the pressure vessel is described with these assumptions: only liquid flows out of the vessel through the bottom outlet and all the vapour is trapped above the liquid level. The second assumption is that the process is adiabatic, which is an estimation since the surfacetovolume ratio of common types of pressure vessels is relatively low and therefore the heat transfer is negligible.
The total volume of agent fill is obtained from the total mass of agent fill and specific volume of the saturated liquid in a state after reaching saturation as described in Section 2.1. The known volume of the container is not considered as a valuable parameter because of the unknown volume of the inert gas fill occupying a portion of the vessel.
2.5 Isentropic expansion in twophase region
The first of the two processes used to describe the liquid release from the vessel mimics the expansion of the remaining contents of the vessel. The concurrent increase of twophase quality is described in the next subsection.
As the mass flow rate is known, the released mass of the liquid and the related specific volume for each time step is determined. This released mass has the volume of Δm_{e} ⋅ v_{eliq (n−1)}, which is taken as the total volume change for the analysed isentropic expansion. The specific volume after the assumed (or though) expansion can then be determined as the hypothetically enlarged total volume divided by the current remaining agent mass − as in (2):$${v}_{en}=\frac{V+\mathrm{\Delta}{m}_{e}\cdot {v}_{eliq\left(n1\right)}}{{m}_{en}}\text{.}$$(2)
The state after such expansion can then be expressed based on s_{en}, v_{en}. As the calculation is built in MS Excel, a tablebased numerical solution is made. The twophase quality is quantified for all available points of the saturation curve for both the entropy s_{en} and the specific volume v_{en}, see (3).$${x}_{ens}=\frac{{s}_{en}{s}_{liq}}{{s}_{vap}{s}_{liq}},\text{\hspace{1em}}{x}_{env}=\frac{{v}_{en}{v}_{liq}}{{v}_{vap}{v}_{liq}}\text{,}$$(3)
The pressuretemperature level which yields the smallest difference between these two values of quality (x_{ens} − x_{env} ≅ 0) is chosen as a resulting state after fictional isentropic expansion along with the quality x_{ens} = x_{env} = x_{enprerelease} .
The resulting temperature and pressure are those present in the vessel at the end of the time step. The fictional increase of volume compensates for the amount of liquid expelled from the bottle. This decrease of amount of liquid also causes an increase of the twophase quality which is described below.
2.6 State change due to liquid release
As the total volume increase discussed in the previous subsection is only assumed, the mass decrease in the pressure vessel is modelled further. An increase in the twophase quality at constant pressure and temperature is calculated. The siphon tube or bottom socket of the vessel ensures that only liquid is released.
The relation of the quality after assumed expansion (prerelease) to the real quality after the concurrent liquid release and expansion is expressed simply in (4).$${x}_{en}={x}_{enprerelease}\cdot \frac{{m}_{n1}}{{m}_{n}}\text{.}$$(4)
This practice is valid if the time step is short enough, so the increase of the quality is small. This may be measured by the difference of the achieved fluid state from the presented calculation and if the expansion is calculated after the liquid release − in the reverse order of that presented here. The proposed order was chosen to avoid the unstable behaviour of the quality increase calculation from saturated liquid state, where x = 0.
3 Thermodynamic state of gas mixture in the room
3.1 Released agent temperature and enthalpy
The isenthalpic throttling happening in the distribution system always starts at the saturated liquid state determined by the pressure in the vessel. This expansion process ends at the pressure level of the room contents, since the agent vaporizes partially and is therefore in the twophase region, the discharge temperature is constrained by this pressure. This temperature is lower than the room temperature in most cases.
The remaining latent heat is expected to be transferred to the sprayed agent droplets rapidly from the current gaseous contents of the room. This portion of the specific enthalpy of state change can be evaluated as follows (5) and the obtained value can be illustrated as the length of horizontal part of the green and cyan thick lines in Figure 3:$${q}_{evap}={h}_{dnvap}{h}_{dn}\text{.}$$(5)
The temperature of the room contents before the mixing calculation is initiated is lowered by the heat consumed by the phase change of the agent as described in equation (6):$${T}_{nrpremix}={T}_{r\left(n1\right)}\frac{{q}_{evap}\cdot \mathrm{\Delta}{m}_{n}}{{c}_{\left(n1\right)r}\cdot {m}_{\left(n1\right)r}}\text{.}$$(6)
Fig. 3 Schematic pressureenthalpy diagram of the extinguishing agent with annotated examples of phenomena described in (2) and (3). Red − saturation curves; blue − isotherms; magenta − isentropes, green − curves of constant twophase quality. Other curves annotated in the figure by text of corresponding colour. (Own work based on NIST REFPROP fluid thermophysical database program.) 
3.2 Carbon dioxide under triple point pressure
There is a serious issue regarding the approach of the twophase area of CO_{2} snow − between solid and vapour states. The process of calculation is the same whether the agent phase change is vaporization or desublimation, but as the thermodynamic properties database REFPROP focuses on fluids, the data in the area of solid and the 2phase solidvapour area are unavailable. There are only data of the saturated vapour curve as the righthand side limit of the snow area. Therefore, constant enthalpy and entropy of vaporization is assumed as found in [5].
Under this assumption, REFPROP data on CO_{2} can be expanded as in Figure 4. Just the values on the saturation curves on both sides of the snow area together with two phase quality value are needed for all necessary evaluations.
Fig. 4 Illustration of REFPROP data expansion. Cyan: melting line, red: saturated vapour and liquid curves, light green: vapourliquid 2phase quality isolines, dark green: entropy isolines, blue: isotherms, black: snow area bounds, right curve is part of REFPROP data. 
3.3 Mixing process
The mixing of the discharged agent with the gaseous contents of the room is estimated using the ideal gas theory. In the first time step, pure air is mixed with the agent and the successive steps are calculated considering the already existing mixture. The agent is considered to be in a saturated vapour state, usually colder than the rest of the room.
The specific heat (7) and specific gas constant (8) of the mixture is a mass weighted average. The temperature is an average weighted by mass and specific heat of the components (9):$${c}_{pnrmix}=\frac{{c}_{pd}\cdot \mathrm{\Delta}{m}_{n}+{c}_{pnr}\cdot {m}_{\left(n1\right)r}}{\mathrm{\Delta}{m}_{n}+{m}_{\left(n1\right)r}}\text{,}$$(7) $${r}_{nrmix}=\frac{{r}_{d}\cdot \mathrm{\Delta}{m}_{n}+{r}_{nr}\cdot {m}_{\left(n1\right)r}}{\mathrm{\Delta}{m}_{n}+{m}_{\left(n1\right)r}}\text{,}$$(8) $${T}_{nrmix}=\frac{{T}_{d}\cdot {c}_{pd}\cdot \mathrm{\Delta}{m}_{n}+{T}_{nrpremix}\cdot {c}_{pnr}\cdot {m}_{\left(n1\right)r}}{{c}_{pd}\cdot \mathrm{\Delta}{m}_{n}+{c}_{pnr}\cdot {m}_{\left(n1\right)r}}\text{.}$$(9) T_{d}, the agent discharge temperature is already the saturation temperature at room pressure, since the heat exchange during evaporation is considered in the previous step of the calculation described in Section 3.1.
3.4 Heat addition from metal bodies
When the room cools down rapidly, some parts may be able to expel its heat content as well as the gases do. There are a few factors to assess the ability to expel heat rapidly. Bodies for this analysis need to be made of highly conductive materials have thin walls and a large surface area exposed directly to the air. Also, the total amount of heat per Kelvin (mass times specific heat) in such bodies must be high enough to have any impact on the whole process. The eligible parts of a room interior are sheet metal and thin metal profiles.
The temperature distribution in the profile is neglected. A single value of bulk temperature is used. For the initialization of an iterative process, the temperature of the metal part is not expected to change from the previous time step (10).$${T}_{mninit}={T}_{m\left(n1\right)}\text{.}$$(10)
Only two iterations of the process presented in (11), (12) give satisfactory results.$${Q}_{mn}=\alpha \cdot {A}_{m}\cdot \left(\frac{{T}_{m\left(n1\right)}+{T}_{mn}}{2}{T}_{nr}\right)\text{,}$$(11) $${T}_{mn}={T}_{m\left(n1\right)}\frac{{Q}_{mn}}{{c}_{m}\cdot {m}_{m}}\text{.}$$(12)
A qualified guess must be made for the heat transfer coefficient α estimation. It may lie somewhere between 10 and 50 W/(m^{2}K). It can be adjusted concerning the exposure of the metal surfaces to the nozzle wake areas, possible contact with liquid particles of the agent and other factors.
This simple method, without consideration of temperature distribution inside the thinwalls, has been proven to differ less than 1% from a more complex analytical calculation of transient heat conduction in an infinite plate. For the given material parameters of aluminium alloy and a thickness of 1.5 mm, the temperature difference between the exposed surface and opposite adiabatic surface is on the order of 0.01 K.
For an idea of the scale of the results achieved and how notable is the submission of heat from metal parts in the room: 5–10% of heat supplied to the agent usually comes from metal; the rest is donated by air. After the period of agent discharge, the metal parts are cooled down by ∼5 K, whilst mixturetometal temperature difference remains above 20 K. Further heat addition from the metal can be expected.
3.5 Control volume leakage
The room is not manufactured to be completely airtight. Some slots and chinks may be present in the outer bounds of the room structure. Overpressure acting on the inner side of the room causes the gases to leak. This may have the positive effect of reducing the overpressure. Underpressure may be considered too and it actually occurs during the initial part of fire extinguishing process.
The amount of leak per time step of the calculation is based on the current overpressure, the surroundings are considered to stay at a constant pressure equal to the initial room pressure. Due to the low discharge velocities, incompressible gas dynamics theory is adopted.
The first versions of the model have used an idea of a virtual orifice of constant area, where the flow rate was linearly dependent on the pressure difference. The early pressure test on the room was carried out on a single overpressure level and the results of this test were then used to quantify the virtual orifice area. Since the structure of the room is significantly pliable and the slots may stretch and spread, this approach was abandoned. Instead, sets of pressure tests were completed to obtain curves of flow rate dependency on the variable overpressure.
To illustrate the deformations: the room ceiling deflection measured during room prototype pressure test at cca 30 mbar overpressure reached approximately 50 mm.
For simple implementation, the leakage curves were approximated by second order polynomial. This respects the expectations of gaining quadratic rise of slot area with the linear deflection of the room wall panels.
Two sets of tests were worked. The first one with a carefully manufactured room prototype without any structural changes, sealed as well as possible by the standard manufacturing procedure. This room has experienced a nominal value of leak rate at 1.5 mbar overpressure equal to 3.7 L/s. The second one was the same prototype with artificially increased leak rate of 8.4 L/s at the same overpressure level.
Eight overpressure levels ranging from 1 to 30 mbar were then tested for both nominal leak rate cases and the obtained leak rates were used for polynomial regression. Only the 3.7 L/s nominal leak rate case is later used in the calculations, as it is on the safe side considering the higher overpressure is generated more easily with low leak rate values. Secondly, the drilled metal plate, which was used in the case of 8.4 L/s nominal leak rate, may not be an accurate representation of a less tight room, especially considering the dynamic behaviour of the leaks.
A leakage curve polynomial derived from 3.7 L/s measurement is written in (13). The nonrealistic zeroth order coefficient was replaced with zero in calculations. Then the expression does not generate a finite leakage rate at no pressure difference between the room and the environment. The downward shift of the curve, which is applied by this substitution is again an adjustment on the safe side. The overpressure is given in mbar gauge pressure. The values which are generated by this expression are shown in Figures 5 and 6.$${{\displaystyle \dot{V}}}_{l}=f\left({p}_{r}\right)=0.19\cdot {p}_{r}^{2}+0.38\cdot {p}_{r}+3.21\text{.}$$(13)
Underpressure conditions were not tested, so the constant area orifice model was kept for the underpressure cases, where it should be on the safe side of generating intake flow rates possibly above realistic values. The typical slots between flat structural panels should become smaller in this case.
Fig. 5 Results of the CO_{2} case. 
Fig. 6 Results of the R403B case. 
3.6 Overpressure and volumetric concentration
Perfect gas behaviour is considered for final overpressure evaluation. The value of p_{r} is adjusted considering the leakage. No iterative process is implemented to correct the influence of overpressure used for leakage quantification. Equation (14) describes the overpressure as a result of isothermal expansion of the gas mixture due to total mass drop by leaks. Section 3.3 clarifies how the specific gas constant and mixture temperature is obtained. The room volume is an input parameter.
An adiabatic model of expansion is possible to employ and may be more accurate. The room temperature T_{r} would then need to be adjusted. Both the noniterative nature of the calculation and isothermal expansion model call for a short time step to keep the accuracy.$${p}_{r}=\frac{{m}_{r}{r}_{r}{T}_{r}}{{V}_{r}}\text{.}$$(14)
The estimation of the volumetric concentration uses the specific volume of saturated vapour of the agent. Only entropy tables were exported to Excel from REFPROP for the vapour single phase area. It is unnecessary to involve other twoparameter tables in the calculation tool just for the concentration evaluation.
4 Results
Below the examples of results obtained from the above described calculations are shown. The input parameters are presented in Table 1. Two cases are presented, while the room parameters including the initial state of air, metal parts properties and leakage description remain the same. Two different extinguishing agents are used: CO_{2} and R403B with different initial weight and pressure, but with the mass flow rate set to achieve 30 s release time for both cases.
The results are included in Figures 5 and 6.
The agent pressures (p_{e} ) drop almost linearly along with constant flow rate. The initial pressure release is not visible in the graph and curves for both agents start at the saturation pressure for given state. For CO_{2}, the pressure does not reach under 5,1 barA where lays the triple point. The sudden change of phase change enthalpy on the entry to the snow twophase area does not have a significant impact on the other processes.
Room overpressure (p_{r} ) and volumetric flow rate of leaks (${{\displaystyle \dot{V}}}_{l}$) tell that the rapid cooling of the room content has larger effect to the pressure in the room than the addition of agent mass in the initial period of extinguishing. During this period, the pressure is lowered in the room causing outside air to be taken inside the room. The temporal width of this phase is varying based on the agent type used.
The dashed lines suffixed “old” in the figures' legend show the behaviour of the former leakage model. Due to linear response of the flow rate expression, the pressure in the room never stops to grow. The new model shows the room is reaching close to equilibrium state, where at certain overpressure (near 20 mbarG) the volume addition of the agent is almost completely compensated by the leaks of the gas mixture out of the room.
These values of the new model seem to be very close to reality, as the experimental test carried out with the halocarbon agent designated for the application reach very similar values of maximum overpressure. In four consecutive tests, the overpressure reached 18.1, 19.3, 18.4 and 18.5 mbarG.
The temperature of room content (T_{r mix}) drops linearly after passing the initial steeper slope. The metal temperature (T_{m}) follows a similar trend, but as the heat transfer rate from the metal parts is limited, the temperature drops slower. Consequently, this creates a relatively high final temperature difference between metal and gas in the room for both cases, which can cause heatup of the gas after the agent release period resulting in further rise of pressure.
The temperature inside agent pressure vessel (T_{e} ) shows the effect of liquid only release. Despite the pressure curve, the concave shape is significant. The R403B case shows the tables might be too coarse to have stepfree smooth curve, but this is believed to have negligible effects on the important results.
Considering the pressure difference between triple point and critical point for both presented substances (and many others), the change of pressure in the room over time is very low. Therefore, the agent discharge temperature (T_{d} ), being the saturation temperature for a given room pressure, is close to constant. Here the demand for tabulated values density or interpolation algorithm for getting a smooth curve is visible too, but again without a significant effect on the results.
The volumetric concentration estimate reaches far above the needs of the needed levels for reliable function of the system. Environmental analysis results with different inputs may vary.
Input parameters of example calculations.
5 Conclusions
A set of coherent calculations was built for the purpose of designing an integrated fire extinguishing system and analysis of various operating conditions. The overpressure in the extinguished room is checked throughout the timedependent calculation.
There are a few key parts of the work.
At first, a strong constraint of constant mass flow rate is set and carefully proven with experimental results and external References.
The second is the thorough decomposition of the thermodynamic processes and the decisions on how to model them. A summary is given in Figure 3. Many issues have been overcome on the side of using REFPROP tables imported as raw real gas data into MS Excel.
Finally, significant influences were identified. The heat of metal bodies and control volume leakage (Sects. 3.4 and 3.5) are incorporated into the calculation. Two models of leakage were used in the course of development of the tool. The second, more accurate model based on experimental results gives new insight on the process. Since the leakage rate can now reach equilibrium with continual agent addition to the room, the maximum overpressure seems not to ever reach critical values for this particular design.
Results of a case study have been discussed. See Section 4.
The gases may also lose a significant amount of internal energy due to the work being done on the wall's deformation. Analysis of this phenomenon may lead to future enhancements of the design tool.
Nomenclature
c : Specific heat (isobaric) [J/kg/K]
q : Heat per kilogram mass [J/kg]
r : Specific gas constant [J/kg/K]
v : Specific volume [kg/m^{3}]
α : Heat transfer coefficient [W/m^{2}/K]
Indices:
Acknowledgments
This work was supported by the project of specific research UWB SGS2019021. Authors would like to express their special thanks to the company Safran Cabin CZ s.r.o. as well as to professor František Maršík from IT CAS.
References
 L. Hurda, R. Matas, Liquefied extinguishing agent discharge to an overpressuresensitive enclosed volume, MATEC Web Conf. 168 , 07010 (2018) [CrossRef] [Google Scholar]
 A.G. Linde, Retrofit Guide, available at https://www.lindegas.com/internet.global.lindegas.global/en/images/Refrigerant%20Retrofit%20guides17_109809.pdf?v=5.0, 2017 [Google Scholar]
 M. Ahmad, L. Buit, O. Florisson, C. HulsboschDam, M. Bögemannvan Osch, M. Spruijt, F. Davolio, Experimental Investigation of CO2 Outflow from a Highpressure Reservoir, Energy Procedia 37 , 3005–3017 (2013) [Google Scholar]
 National Fire Protection Association, NFPA Code 12A: Standard on Halon 1301 fire extinguishing systems, 2018 [Google Scholar]
 R.M. Stephenson, S. Malanowski, Handbook of the Thermodynamics of Organic Compounds, Elsevier, New York, 1987 [CrossRef] [Google Scholar]
Cite this article as: L. Hurda, R. Matas, M. Redmer, Analysis of room overpressure after launching an integrated extinguishing system, Mechanics & Industry 20, 706 (2019)
All Tables
All Figures
Fig. 1 Discharge process for three experimental tests on a room prototype with CO_{2} used for verification of the constant mass flow rate. Thick grey lines indicate the linear approximation of mass, i.e. approximation of constant mass flow rate. 

In the text 
Fig. 2 Simplified experiment layout. Standard integrated pressure vessel (grey) replaced with external pressure bottles of CO_{2} for parallel discharge controlled initiated by manual release valve (green). 

In the text 
Fig. 3 Schematic pressureenthalpy diagram of the extinguishing agent with annotated examples of phenomena described in (2) and (3). Red − saturation curves; blue − isotherms; magenta − isentropes, green − curves of constant twophase quality. Other curves annotated in the figure by text of corresponding colour. (Own work based on NIST REFPROP fluid thermophysical database program.) 

In the text 
Fig. 4 Illustration of REFPROP data expansion. Cyan: melting line, red: saturated vapour and liquid curves, light green: vapourliquid 2phase quality isolines, dark green: entropy isolines, blue: isotherms, black: snow area bounds, right curve is part of REFPROP data. 

In the text 
Fig. 5 Results of the CO_{2} case. 

In the text 
Fig. 6 Results of the R403B case. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.