High pressure jet
A high pressure jet is a stream of pressurized fluid that is released from an environment at a significantly higher pressure than ambient pressure from a nozzle or orifice, due to operational or accidental release. In the field of safety engineering, the release of toxic and flammable gases has been the subject of many R&D studies because of the major risk that they pose to the health and safety of workers, equipment and environment. Intentional or accidental release may occur in an industrial settings like natural gas processing plants, oil refineries and hydrogen storage facilities.
A main focus during a risk assessment process is the estimation of the gas cloud extension and dissipation, important parameters that allow to evaluate and establish safety limits that must be respected in order to minimize the possible damage after a high pressure release.
Mechanism and structure of a gaseous jet
Subsonic and sonic flow
When a pressurized gas is released, the velocity of the flow will heavily depend on the pressure difference between stagnant pressure and downstream pressure. By assuming an isentropic expansion of an ideal gas from its stagnant conditions to downstream conditions, the subsonic flow rate of the source term is given by Ramskill's formulation:File:Gas flow rate variation as the pressure ratio decreases.png|thumb| The flow rate of an ideal gas can be represented by the graphed line. As the pressure ratio decreases and the critical value is approached, the flow changes from non-choked to choked, setting an upper limit to the velocity of the gas to the speed of sound of the medium.
As the ratio between downstream condition pressure and stagnant condition pressure decreases, the flow rate of the ideal gas will increase. This behavior will continue until a critical value is reached, changing the condition of the jet from a non-choked flow to a choked flow. This will lead to the a newly defined expression for the aforementioned pressure ratio and, sub-sequentially, the flow rate equation.
The critical value for the pressure ratio is defined as:
This newly defined ratio can then be used to determine the flow rate for a sonic choked flow:
The flow rate equation for a choked flow will have a fixed velocity, which is the speed of sound of the medium, where the Mach number is equals to 1:
It is important to note that if P1 keeps on decreasing, no flow rate change will occur if the ratio is already below the critical value, unless P0 also changes.
Under-expanded jet structure
An under-expanded jet is one that manifests when the pressure at downstream conditions is greater that the pressure of the environment where the gas is being released in. It is said to be under-expanded since the gas will expand, trying to reach the same pressure of its surroundings. When under-expanded, the jet will have characteristics of a compressible flow, a condition in which pressure variations are significant enough to have a strong effect on the velocity, density and temperature. It is important to note that as the jet expands and incorporates gases from the surrounding medium, the jet will behave more and more like an incompressible fluid, allowing for a general definition of the structure of a jet to be the following:- Nearfield zone: this zone is composed of a core layer that is isolated from the surrounding medium, with its behavior being mostly dominated by compressible effects, and an outer layer that is in contact with the surrounding medium fluid. Due to turbulent effects, the outer layer, nominated as mixing layer, permits gas entrainment as it is facilitated, diluting the jet. In this shearing zone, a subsonic and supersonic section may be distinguished, where temperature, density and pressure vary wildly in a few centimeters of distance from the source. This zone has the characteristics of a compressible fluid.
- Transition zone: the beginning of this zone represents the ending of the nearfield zone, where variations are small compared to the previous one. Density and temperature variations are mostly because of mixing with the surrounding fluid.
- Farfield zone: this final zone is one of a fully expanded and incompressible jet. Longitudinal velocity and temperature are now inversely proportional to the distance from the source and radial evolution can be described by a gaussian dispersion model. It is important to note that this zone can be further split into inertial, buoyant and turbulent zones.
Under-expanded jet classification
This process will repeat until the pressure difference between ambient pressure and jet pressure is null.
Compression and expansion are accomplished through a series of shock waves, formed as a result of Prandlt-Meyer expansion and compression waves.
Development of the aforementioned shock waves will be related to the difference in pressure between the stagnant conditions or downstream conditions and the ambient conditions, as well as the mach number. With varying pressure ratios, under-expanded jets can be classified as:
- Moderately under-expanded jet: Nearfield with diamond shaped structures. A Prandlt-Meyer expansion generates oblique expansion waves that expand the fluid downstream from the exit orifice. As these waves attain constant pressure from the surrounding fluid, they are deflected back as compression waves, converging in oblique shock waves. When they meet on the axis of the jet, reflected shock waves move outwardly until they attain constant pressure from the surrounding fluid, repeating the process, and in turn, recreating the cell structure.
- Highly under-expanded jet: Nearfield with barrel shaped structures. As the pressure ratio increases, the intercepting shock waves can't meet on the axis of the jet anymore, which forces the generation of a normal shock wave when the intercepting shock waves go beyond a certain critical angle. From the interception point of the mach disk and the intercepting shock, a residual slipstream will reflect outwardly, until it reaches constant pressure from the surrounding fluid, repeating the process, recreating the barrel shaped cell structure.
- Extremely under-expanded jet: Nearfield with a single cell structure. When the pressure ratio goes beyond a critical value, cell numbers within the nearfield of the jet decrease, until they all coalesce into a single cell with a single mach disk. Due to the increase in velocity and lower pressure zones around the jet, ambient fluid entrainment will increase.
Natural gas release
Three parameters that must be considered when dealing with flammable gasses are their flash point, upper flammability limit and lower flammability limit, as they are set values for any compound at a specific pressure and temperature. If we consider the fire triangle model, to induce a combustion reaction three components are needed: a fuel, an oxidizing agent and heat.
When release happens in an ambient filled with air, the oxidizing agent will be oxygen. At an almost pure concentration, a few centimeters from the exit plane, the concentration of natural gas is too high and oxygen too low to generate any kind of combustion reaction, but as the high pressure jet develops, the concentration of its components will dilute as air entrainment increases, allowing an enrichment of oxygen within the jet. Assuming a constant concentration for oxygen, the jet must dilute enough to enter within its flammability range; below its UFL. Within this range, a flammable mixture can be made and any source of heat can jump-start the reaction.
To properly judge the damage and potential risk that the jet fire can generate, several studies regarding the maximum distance that the cloud generated by the jet can reach have been made. As dilution of the jet continues due to air entrainment in the farfield, going below its UFL, the maximum distance that the flammable mixture can reach is at the point in which the concentration of the cloud is equals to the LFL of the gas, as it is the lowest concentration allowable that permits the formation of a flammable mixture between air and natural gas at standard conditions.
Considering a free jet at sub-critical pressure, its mean volume fraction axial concentration decay of any gas released in air can be defined as follows:
Computational Fluid Dynamics">Computational fluid dynamics">Computational Fluid Dynamics
of high pressure jets have to be limited in terms of size and complexity of the scenario due to the inherit dangers and expenses correlated to the experiment itself. Alternative methods to gather data, such as representative models, can be used in order to predict what the maximum extend of the gas cloud at its LFL concentration can reach. Simpler models like a gaussian gas dispersion model or integral model can be useful to have a quick and qualitative overview on how the jet may extend. Unfortunately, their inability to properly simulate jet-obstacle interactions make them impossible to use beyond preliminary calculations. This is the reason why Computational Fluid Dynamic simulations are generally preferred for more complex scenarios.Although there exists several approaches for CFD simulations, a common approach is the use of a finite volume method that discretizes the volume into smaller cells of varying shapes. Every single cell will represent a fluid-filled volume where the scenarios parameters will be applied. Every cell that was modeled solves a set of conservation equations of mass, momentum and energy, along with the continuity equation. Fluid-obstacle interaction is then modeled with varying algorithms based on the closure turbulent model used.
Depending on the number of total cells within the volume, the better the quality of the simulation, the longer the simulation time. Convergence problems can arise within the simulation as large momentum, mass and energy gradients appear in the volume. The points where these problems are expected to appear need to have a higher number of cells to achieve gradual changes between one cell and another. Ideally, through CFD simulations, a simpler model can be derived which, for a specific set of scenarios, allows to have results with an accuracy and precision level similar to the CFD simulation itself.