Pipeng Toolbox : Critical Flow Modules Login
Short Cuts
GO
Main ±
Beams ±
References ±
Fluid Flow ±
Fluid Properties ±
Maths ±
Materials ±
Pipelines ±
Soils ±
Subsea ±
Data ±
Units ±
Help ±
Demo

Critical Flow Modules

Links : ±
CALCULATOR MODULE : Compressible Flow   ±

Calculate compressible flow ratios and gas properties for isentropic and isothermal flow (critical over stagnation ratios, flowing over stagnation ratios, and flowing over critical flow ratios).

For isentropic flow, critical flow occurs at M = 1. For isothermal flow, critical flow occurs at M = 1 / √k, where k is the specific heat ratio (Cp/Cv). For isothermal flow the isothermal temperature is assumed equal to the stagnation temperature. Phase changes are ignored.

For flow through a throat, the flow upstream from the throat is sub critical (M ≤ Mc). The flow downstream is super critical (M > Mc). The area ratio is inversely proportional to the mass flux ratio. At stagnation conditions, the area ratio is infinite.

Use the Result Plot option to plot flow ratios versus Mach number, or nozzle area ratio and diameter ratio versus Mach number.

Reference : Fluid Mechanics, Frank M White, McGraw Hill

Change Module :

Related Modules :

CALCULATOR MODULE : Compressible Flow Critical Flow   ±
CALCULATOR MODULE : Compressible Flow Pressure Relief Valve   ±

Calculate compressible flow pressure relief valve capacity for isentropic, isothermal, and adiabatic conditions.

For pressure relief valves with no header, the mass flow rate can be calculated for isentropic or isothermal flow. The pressure relief valve is assumed to exit directly to ambient pressure. If the ambient pressure is less than the critical pressure the flow is critical (Mc = 1 for isentropic flow and Mc = √(1/γ) for isothermal flow). If the ambient pressure is greater than the critical nozzle pressure, the flow is sub critical (M < Mc). For isothermal flow the stagnation temperature should be close to or equal to the ambient temperature (for example a gas transmission pipeline). Phase changes are ignored.

For a combined pressure relief valve and pressure relief header, the mass flow rate can be calculated for

  • Isentropic nozzle and adiabatic header
  • Isentropic nozzle and isothermal header
  • Isothermal nozzle and isothermal header

The pressure relief valve is assumed to exit directly into the header. If the header inlet pressure is less than or equal to the nozzle critical pressure the nozzle flow is critical (Mc = 1 for isentropic flow and Mc = √(1/γ) for isothermal flow), and the mass flow rate is restricted by the nozzle. The header inlet pressure is calculated so that the header mass flow rate equals the nozzle mass flow rate. If the header inlet pressure is greater than the critical nozzle pressure, the nozzle flow is sub critical (M < Mc), and the mass flow rate is restricted by the header. The mass flow rate is calculated so that the header inlet pressure is equal to the nozzle pressure. The mass flow rate through the nozzle is always equal to the mass flow rate through the header.

For a pressure relief header, the mass flow rate can be calculated for adiabatic or isothermal flow. If the header is attached directly to the outlet of a pressure relief valve (PRV), the header inlet pressure should be set equal to the PRV nozzle outlet pressure. The header should be sized so that the calculated header mass flowrate is greater than or equal to the PRV mass flowrate. For headers with more than one PRV, the header mass flowrate is divided by the number of PRV's. If the header is oversized, the header inlet pressure will reduce so that the actual header mass flowrate is equal to the nozzle mass flowrate (there is a pressure drop between the PRV outlet and the header inlet).

Note : If the PRV is attached to a small diameter header which feeds into a larger diameter header (possibly with multiple PRVs), the large diameter header should be sized first. The inlet pressure for the large diameter header is used as the ambient pressure for the smaller diameter header (and PRV).

Header pressure losses are calculated from the pressure loss factor (fld = fL/D + K). The Darcy friction factor f is calculated for fully turbulent flow using the rough pipe equation. The header is assumed to be constant diameter. Minor losses can be included by the minor loss K factor, and should include valves and bends etc. The header entry and exit losses should not be included (the fluid dynamic pressure loss is included in the calculation). The discharge coefficient can also be used for minor losses, and as a safety factor. If the ambient pressure is less than the critical header pressure the header exit flow is critical (Mc = 1 for isentropic flow and Mc = √(1/γ) for isothermal flow). If the ambient pressure is greater than the critical header pressure, the header exit flow is sub critical (M < Mc). The header entry flow is assumed to be sub critical for all flow conditions.

The effective PRV valve nozzle orifice diameter and cross section area can be calculated from API letter designation (API 526 type D to T). API effective orifice sizing is used to compensate for the friction pressure losses in the relief valve. The combination of valve and nozzle orifice must be tested with the operating fluid at the design conditions, and certified as having a flow rate greater than or equal to the calculated flow rate for the equivalent size. The API 526 orifice sizing assumes isentropic flow. For certified API 526 valves, the isentropic nozzle calculation option should be used.

Note : The pressure relief header calculation is not suitable for pressure relief vents. Headers are assumed to be part of a PRV system. Vents are constant diameter piping attached to a pipeline or pressure vessel.

Use the Result Plot option to plot pressure, mach number and mass flow rate.

Reference : Fluid Mechanics, Frank M White, McGraw Hill

Change Module :

Related Modules :

CALCULATOR MODULE : Compressible Flow Pressure Relief Vent   ±

Calculate compressible flow pressure relief vent flow rate and pressure drop for either adiabatic or isothermal flow.

The vent is modelled as a frictionless entry, combined with a frictional constant diameter duct. For adiabatic flow the vent entry is assumed to be isentropic. For isothermal flow, the vent entry is assumed to be isothermal. The vent entry is assumed to be subsonic at all conditions. The pipeline is assumed to be at stagnation conditions (M = 0). At high pressure the vent exit flow is critical flow (Mc = 1 for adiabatic low and `Mc = 1 / (√γ)` for isothermal flow : γ = the gas specific heat ratio). At lower pressures the vent exit flow is sub critical (M < Mc).

Vent flow rate is calculated from the vent pressure loss factor (fld).

`fld = fd L/D + K `

where :

fld = vent pressure loss factor
fd = Darcy friction factor
L = vent length
D = vent inside diameter
K = minor loss K factor

The Darcy friction factor is calculated assuming fully turbulent flow. Minor losses should include the vent entry, and valves, bends etc.. The vent exit should not be included (the fluid dynamic pressure is included in the calculation). The discharge coefficient can be used as a safety factor.

Note : The vent calculation is not suitable for pressure relief headers which are part of a pressure relief valve (PRV) system.

Use the Result Plot option to plot inlet and exit pressure versus stagnation pressure, inlet and exit mach number versus stagnation pressure, or mass flow rate versus stagnation pressure and flow type.

Reference : Fluid Mechanics, Frank M White, McGraw Hill

Change Module :

Related Modules :

CALCULATOR MODULE : API 520 Critical Flow Ratio   ±
CALCULATOR MODULE : IAPWS R7-97 Steam Critical Flow   ±

Calculate IAPWS R7-97 isentropic steam critical flow properties from stagnation temperature and pressure.

Flow properties can be calculated for either critical flow, or from a user defined flowing pressure. Flow properties are valid for the vapour phase only. For critical flow the mass flux is a maximum. theoretical critical Mach number equals 1. The Mach number will vary for a user defined flowing pressure. The flowing velocity is calculated from the stagnation enthalpy (ho = h + 1/2 V^2). The anomaly zone is set to region 2 (region 3 does not converge properly). Use the Result Plot option to plot isentropic flowing properties and isentropic temperature versus either flowing pressure or Mach number.

Note : There is an anomaly in the steam calculation for region 3 between the saturated vapour line, the region 2/3 boundary, and the critical pressure. Refer to the region 3 anomaly help page for more details (click the utility button on the data bar). IAPWS R7-97 is intended for industrial use, and is a simplified version of IAPWS R6-95 for scientific use. IAPWS R7-97 was developed as an improvement of the IFC-67 model.

Reference : IAPWS R7-97 Industrial Formulation for thermodynamic Properties of Water and Steam

Change Module :