1.

PRESSURE LOSS AT FLUID FLOW AND ITS CALCULATION

Jiří Škorpík, skorpik.jiri@email.cz
1.3 . . . . . . . . . . . . . .
1.4 . . . . . . . . . . . . . .
1.5 . . . . . . . . . . . . . .
1.7 . . . . . . . . . . . . . .
1.8 . . . . . . . . . . . . . .
1.11 . . . . . . . . . . . . . .
1.13 . . . . . . . . . . . . . .
1.14 . . . . . . . . . . . . . .
1.14 . . . . . . . . . . . . . .
1.17 . . . . . . . . . . . . . .
1.19 . . . . . . . . . . . . . .
1.21 . . . . . . . . . . . . . .
1.22 . . . . . . . . . . . . . .
1.2
Author:
ŠKORPÍK, Jiří, ORCID: 0000-0002-3034-1696
Issue date:           
April 2010, May 2021, June 2023, April 2024 (4th ed.)
Title:
Pressure loss at fluid flow and its calculation
Journal:
Transformační technologie (transformacni-technolgie.cz; fluid-dynamics.education; turbomachinery.education; engineering-sciences.education; stirling-engine.education)
ISSN:
1804-8293

Copyright©Jiří Škorpík, 2024
All rights reserved.

 Chapter: Description of pressure loss development and basic concepts
1.3

Description of pressure loss development and basic concepts

Pressure loss

Internal friction

During the fluid flow, friction on the surface of the flow channel and the flowing bodies as well as friction inside the fluid (so-called internal friction). Through friction, the fluid loses kinetic energy and in order to flow through the channel at the required velocity (flow rate), it must gain kinetic energy at the expense of pressure energy - a pressure loss Lp is created, or at the expense of other energy, such as potential energy, etc.

Pressure loss in pipe

Incompressible fluid

Figure 1 shows the simplest case of pressure drop occurring when an incompressible fluid flows in a constant cross-sectional area pipe. Since the inlet and outlet of the duct must have the same flow, and therefore velocity, with no change in potential energy, the pressure drop Lp is equal to the pressure drop between the inlet and outlet, see Equation 1a.

Fluid friction in pipe and its consequences
1: Fluid friction in pipe and its consequences
A [m2] flow area; F [N] frictional force acting between channel wall and fluid; l [m] investigated length of channel; Lp [Pa] pressure loss on investigated length of pipe; Lq [J·kg-1] heat loss due to internal fluid friction; p [Pa] pressure; ρ [kg·m-3] density of working fluid. The index i indicates the inlet, the index e the outlet. The derivation of the equations is shown in Appendix 4.

Loss heat

The loss heat Lq, which is generated during friction, heats the working fluid. The loss heat for the case of Figure 1 corresponds to the pressure energy of the pressure loss, see Equation 1(b).

Frictional force

The fluid acts with a frictional force F on the channel in the direction of flow. The friction force for the case of Figure 1 corresponds to the product of the pressure difference between the inlet and outlet of the channel (pressure loss) and the flow area of the channel, see Equation 1(c).

Pressure loss in pipeline

Cryogenics

Heart

For normal technical practice, addressing pressure loss in piping networks with valves is essential.Determining the pressure loss helps to calculate the work of the pump or fan - part of this work is consumed by the pressure loss. Heat loss calculations are crucial in cryogenics for transporting liquefied gases via pipelines, preventing property loss or evaporation. In the body, blood flow loss is compensated by the heart's activity, with output linked to pressure differentials (diastolic and systolic pressures).

 Chapter: Description of pressure loss development and basic concepts
1.4

Calculation of pressure loss

During fluid transport, incompressible fluid theories suffice for calculation pressure loss due to stable density. Yet, for gas transport in lengthy pipelines, density variations occur. Here, pressure loss calculations rely on mean gas density [Mikula et al., 1974, p. 71] or more precisely on friction-induced pressure loss equations in later paper chapters.

Laminar flow

Turbulent flow

Mean flow velocity

Reynolds number

characteristic length

Kinematic viscosity

Critical Reynolds number

Upper critical Reynolds number

The procedure for calculating the pressure drop in the channel under investigation depends on whether the channel is in laminar or turbulent flow. This can be determined by the value of the Reynolds number for the case in question, for the calculation of which it is necessary to know the mean flow velocity of the fluid, the characteristic lenght of the channel (in the case of pipes this is the diameter) and the value of the kinematic viscosity. If the value of the Reynolds number is less than the value of the critical Reynolds number, then the flow is likely to be laminar, if the value of the Reynolds number is less than the value of the upper critical Reynolds number, then the flow is likely to be turbulent.

Nozzles

Diffusers

Profile cascades

We also identify the pressure loss in channels designed for the transformation of pressure and kinetic energy of the fluid such as nozzles and diffusers or profile cascades, but in these cases the pressure loss is defined indirectly - the problem of losses in these channels is described in the articles Flow of gases and steam through nozzles4., Flow of gases and vapours through diffusers5., Aerodynamika profilových mříží [Škorpík, 2022a].

Equations for calculating pressure loss in pipe

Navier-Stokes equation

Darcy-Weisbach equation

Henry Darcy

Julius Weisbach

Pipe fittings

Valves

The relationship for calculating the pressure loss for the case of laminar steady flow as function of dynamic pressure can be derived from the Navier-Stokes equations. This equation is called the Darcy-Weisbach equation, which was developed by the French engineer Henry Darcy (1803-1858) for pipelines, see Equation 2. Later, on the basis of long term experiments and deduction, the German engineer Julius Weisbach (1806-1871) confirmed the validity of this relationship for transient and turbulent flows and even for losses in pipe fittings and valves.

2: Darcy-Weisbach equation for calculating pressure loss
ζ [1] loss coefficient of section related to kinetic energy of mean velocity (defined by Weisbach [Maštovský, 1964, p. 82]); V‾ [m·s-1] mean velocity of mass flow (mean flow velocity).
 Chapter: Equations for calculating pressure loss in pipe
1.5

Loss coefficient

From the Darcy-Weisbach equation it follows that the pressure loss is a certain fraction of the dynamic pressure, this fraction is given by the value of the loss coefficient. For channels of constant flow area, the loss coefficient can be calculated using the equations given in the chapter Calculation of pipe loss coefficient. For other types of channels, such as elbows, valves, etc., see chapter Pressure loss coefficient of local losses.

Calculation of pipe loss coefficient

Friction coefficient

Characteristic length

The loss coefficient for a constant cross-section pipe is computed using Equation 3. It is therefore a function of the length and diameter of the pipe (d is taken to be the characteristic length if the pipe is of non-circular cross-section [Mikula et al., 1974, p. 91]) and a quantity called the friction coefficient.

3: Equation of pipe loss coefficient
d [m] internal diameter of pipe; l [m] length of pipe; λ [1] friction coefficient in pipe on pipe section under investigation.

Friction coefficient

Laminar flow

Johann Nikuradse

Nikuradse chart

Pipe roughness

The equation for the pipe friction coefficient in laminar flow λLF can be easily derived from the Navier-Stokes equations, see Equation 4. The determination of the value of the friction coefficient at turbulent flow is based on the conclusion of measurements made by Johann Nikuradse on a series of glass pipes with artificial roughness using a sand film. Nikuradse measured the pressure drop of several pipes with different relative surface roughnesses for selected Reynolds numbers and from there calculated the values of the friction coefficient λ according to the Darcy-Weisbach equation (Equation 2). From these values he produced a chart of the dependence of the friction coefficient on the Reynolds number and confirmed the existence of four regions with different dependences of the friction coefficient on the Reynolds number, see Figure 4.

Nikuradse chart
4: Nikuradse chart
 Chapter: Calculation of pipe loss coefficient
1.6
left-practical division of Nikuradse chart into basic areas; right-view of original Nikuradse chart [Nikuradse, 1933]. (A) the friction coefficient is linear function of only Reynolds numbers without influence of pipe roughness - laminar flow region; (B) transition region of flow from laminar to turbulent - both laminar and turbulent flow can occur; (C) turbulent flow region, in which friction coefficient is function of both Reynolds numbers and relative roughness of pipe; (D) turbulent region, in which friction coefficient is function of relative roughness of pipe – the higher relative roughness, the greater coeicient of friction. C [1] relative pipe roughness, see also Nomogram 8; Re [1] Reynolds number; ReC [1] critical Reynolds number; λLF [1] friction coefficient for laminar flow, see Appendix 5 for derivation of equation; λSP [1] friction coefficient for turbulent flow in hydraulically smooth pipes (C→0) [Bauer et al., 1950, p.148] ; λRP [1] limit from which friction coefficient does not change with increasing Reynolds number, the so-called flow in hydraulically rough pipe [Bašta, 2003, p. 23]; ε [m] absolute roughness of inner walls of pipe (for values see, for example, [Mikula et al., 1974], Table 7).

Friction coefficient

Turbulent flow

Colebrook equation

Cyril Colebrook

Lewis Moody

Moody chart

To calculate the friction coefficient in the regions (C-D) in Figure 4 (C-D) regions, semi-empirical relationships approximate Nikuradse chart values, supplemented by other measurements. An overview of these equations is given, for example, in [Štefan, 2009]. A universal and accurate equation for engineering, by Cyril Colebrook (1910-1997) [Míka, 1977, p. 150] (Equation 5), is central. The American engineer Lewis Moody (1880-1954) then created the widely used Moody chart, derived from the Colebrook equation [Cihelka et al., 1975, p. 684], [Roček, 2002, p. 230].

5: Colebrook equation for calculation of pipe friction coefficient

Smooth pipe

Rough pipe

In region (C), turbulent velocity profile develops. In region (D), the evolution is already complete and even with increasing Reynolds number, the ratio of the kinetic energy of the fluid in the boundary layer to the kinetic energy in the flow core does not change.

Marginal Reynolds number

The values of the marginal Reynolds numbers ReRP, i.e. the approximate boundary between regions (C) and (D), can be calculated by substituting the equation for λRP into the Colebrook equation. Selected values of the marginal Reynolds numbers calculated in this way are given in Tabulce 6.

C 1·10-6 1·10-5 1·10-4 0,001 0,01 0,01 0,04 0,05
ReRP 2,62·109 2,22·108 1,82·107 1,42·106 2,28·105 1,02·105 1,95·104 1,48·104
6: Approximate values of marginal Reynolds number
C [1]; ReRP [1] marginal Reynolds number at which the friction coefficient ceases to be sensitive to the change in Re
 Chapter: Calculation of pipe loss coefficient
1.7
     ε      ε
Drawn pipes (new) from: Copper, Brass, glass 0,001...0,002 Cast iron 0,2..0,6
Plastic and ruber 0,0015...0,007 Galvanized steel 0,07...0,1
Steel 0,04...0,1 Corroded steel pipes cleaned 0,15...0,2
Weld steel 0,04...0,1
7: Approximate values of absolute pipe roughness
[mm]. Selection from [Mikula et al., 1974].
8: Nomogram for calculating relative roughness
d [mm], ε [mm], C [1].

Gas Pipeline

The Darcy-Weisbach equation implies that for minimum pressure loss it is advantageous to transport gas at higher pressures and densities than at low pressures and high velocities. Therefore, the pressures in transit pipelines are around 7 MPa and the gas pressure is reduced before the appliances (see Table 9), which are designed for lower pressures for safety reasons.

p p
Transit pipeline 7,5 Medium pressure gas pipeline 0,1...0,3
High pressure gas pipeline 4 Low pressure (household) 0,002
9: Overpressures in natural gas pipelines
p [MPa] overpressure in gas pipeline.

Pressure loss per unit lengt of pipe

For basic pipe route designs, designers use the quantity specific pressure loss in the pipe corresponding to the pressure loss in a 1 m long pipe, see alsoNomogram 10.

 Chapter: Pressure loss per unit lengt of pipe
1.8
10: Nomogram for calculation of specific pressure loss, dynamic pressure and specific kinetic energy of fluid in pipe
pd [Pa] medium dynamic flow pressure; d [mm], Q [m3·s-1], V‾ [m·s-1], ρ [kg·m-3], λ [1], πL specific pressure loss [Pa·m-1].

Pressure loss coefficient of local losses

Pipe network

The pipe route (pipe network) is not usually straight and may consist of other pipe elements (branch pipes of various shapes, bends, constrictions), fittings, filters, meters and other flow parts, see Figure 11. These elements are local resistances and local pressure loss occurs in them.

 Chapter: Pressure loss coefficient of local losses
1.9
Example of a pipeline route with local resistences
11: Example of a pipeline route with local resistences
a-gate valve; b-closing valve (generally has higher pressure loss than gate valve); c-standard tee; d-narrowing of pipe; e-elbow.

Local resistance

Throttling

In local resistances, the pressure losses are similar to those in a straight pipe. These pressure losses tend to be much more intense than in a straight section of pipe due to the fact that the flow through these sections also changes the shape of the flow channel, the direction of the flow and often the fluid throttling. Inlets and outlets of the pipe can also be considered as a special case of local resistance. At the edges, the flow is usually unsteady and influenced by the shape of the beginning or end of the pipe.

Pressure loss in local resistance

Mean flow velocity in local resistance

The local resistance pressure loss can also be calculated using Equation 2. The calculation of the pressure loss occurring in a given local resistance is based on the mean flow velocity before the element and the loss coefficient of the local resistance.

Loss coefficients of local resistances

Reynolds number

The loss coefficient ζ of some types of local resistances can be calculated [Maštovský, 1964, p. 85], but more often it is based on measurements of the local resistance for different Reynolds numbers. However, for some types of local resistances the influence of the Reynolds number is not significant and tabulated values can be used, especially for valves and pipe fittings, e.g. in [Cihelka et al., 1975, p. 672], [Miller et al., 1972, p. 252], [Řasa and Švercl, 2004, p. 737]. The corresponding loss factor is provided by the manufacturer of the local resistor in question. Loss coefficients for different types of pipe edges are given in [Ibler et al., 2002, p. 268].

 Chapter: Pressure loss coefficient of local losses
1.10

Valve loss coefficient

Valve flow coefficient

In the case of valves, the manufacturer usually also directly supplies charts of the dependence of their pressure loss on the flow rate (depending on the type of flowing medium). If the nominal flow coefficient of a KVS valve is known, the loss versus flow can be calculated from the definition of the flow coefficient, see Equation 12. The nominal flow coefficient is measured on the 2·d pipe section upstream of the valve and the 8·d pipe section downstream of the valve, so the loss coefficient calculated in this way includes this length of pipe - so the actual loss coefficient of the valve is lower by the loss coefficient corresponding to a 10·d smooth pipe. Approximate values of loss coefficients of some valves are given in [Roček, 2002, p. 231, 232]. However, there are other types of coefficients, usually derived from the pressure loss of the valve. It depends on the manufacturer what methodology he uses to compare valves. The relevant relationships are then given in the valve catalogue.

12: Calculation of valve loss coefficient
d [mm] internal diameter of inlet of valve; KVS [m3·h-1] nominal flow coefficient of valve. The relation is derived for the water flow rate in [Roček, 2002, p. 236].

Closing valve

When selecting the most suitable closing valve, the allowable pressure loss Lp is first determined at the volume flow rate Q and the density of the flowing medium at the inlet ρ. The nominal flow coefficient KVS is calculated. Next, the valve with the next higher KVS is selected from the valve catalogue of the relevant manufacturer.

Equivalent pipe length

Smooth pipe

For approximate calculation of the local resistance pressure loss, a quantity called the equivalent pipe length can also be used. This quantity gives the length of smooth pipe (expressed as the number of diameters of smooth pipe) of the same diameter as the input diameter of the local resistance under investigation with the same pressure loss. Equivalent pipe lengths of some valves and pipe fittings are given in [Izard, 1961], [Fraas, 1989], and a selection is given in Table 13. The advantage is that in the calculation it is sufficient to add the individual equivalent lengths and calculate their total pressure loss as if they were the same length of hydraulically smooth pipe, see Problem 1.

l·d-1 l·d-1
GLOBE VALVES
with no obstruction 340   Y-pattern with stem 60° from run of pipeline 175  
 Chapter: Pressure loss coefficient of local losses
1.11
l·d-1 l·d-1
with guided in flow area (under seat) 450 Y-pattern with stem 45° from run of pipeline 145
ANGLE VALVES
with no obstruction 145 with guided in flow area (under seat) 200
GATE VALVES
convetional wedge 13 conduit pipeline 3
pulp stock 17
CHECK VALVES
conventional swing 35 in-line ball 150
clearway swing 50 foot valves with strainer with poppet lift-type disc 420
globe 340 foot valves with strainer with leather-hinged disc 75
angle 145 butterfly valves 20
COCKS
rectangular plug port area equal to 100% of pipe area 18 three-way 140
FITTINGS
90° standard elbow 30 square corner elbow 57
45° standard elbow 16 180° close pattern return bend 50
90° long radius elbow 20 standard tee with flow through run 20
90° street elbow 50 standard tee with flow through branch 60
45° street elbow 26
FLOW MEASUREMENT
turbine flow meter 150   orifice plates 200  
piston meter 400
13: Equivalent pipe length l·d-1 some valves and pipe fittings
l·d-1 [1] equivalent pipe length. Choice of [Fraas, 1989], [Izard, 1961].

Economic velocity in pipe

Pipeline costs vs. pumping work

The Darcy-Weisbach equation links higher mean flow velocity to increased pressure loss. This, in turn, impacts the cost of acquiring and operating machinery (e.g., pump, fan). Larger pipeline diameters, reducing mean flow velocity, raise costs for pipeline routes and fittings. Economic velocity, balancing pressure loss and acquisition expenses, is precisely calculated in [Krbek et al., 1999, p. 187]. Usual economic velocities, balancing costs, are derived from this compromise [Mikula et al., 1974, p. 141], as shown in Table 14. However, factors like layout considerations may influence velocities beyond economic reasons.

 Chapter: Economic velocity in pipe
1.12
V V
oil 1...2 steam superheated to 4 MPa 20...40
water 1...4 steam superheated at high pressure 30...60, 80
low pressure heating steam 10...15 exhaust steam (after expansion in machine) 15...30
sat. steam up to 1 MPa 15...20 air (compressed) 2...15
14: Economic velocity values in pipes of different working substances
V‾ [m·s-1]

Calculation of pipe diameter

The pipe diameter d is calculated from the design economic velocity, density and required specific flow rate, see Nomogram 15. The calculated pipe diameter must be rounded off according to the manufactured pipe diameters corresponding to the pressure and temperature at which the pipe will be operated.

15: Nomogram for pipe diameter calculation
 Chapter: Economic velocity in pipe
1.13
V‾ [m·s-1], ρ [kg·m-3], m [kg·s-1] mass flow; mm [kg·min-1], mh [kg·h-1], Q [m3·s-1] volume flow; Qm [m3·min-1], Qh [m3·h-1] volumetric flow rate through pipe, d [mm] pipe diameter.

Characteristics pipeline

Definition of piping system characteristics

The characteristic pipeline is the dependence of the pressure loss of the pipeline route on the volumetric flow rate. From the equation for calculating the pressure loss it is clear that at ρ=const. the pressure loss will be a quadratic function with a parameter CS called the pipeline system constant, see Equation 16.

16: Characteristics pipeline
n [-] number of individual pipeline sections (each section has constant diameter along entire length); k [-] number of local resistances; Lpipe [Pa] pressure loss of pipeline section; Lcomponent [Pa] pressure loss of local resistance; CS [kg·m-7] pipeline system constant; Q [m3·s-1] volumetric flow. Lp,n [Pa] pressure loss at nominal flow Qn through system. The equation is also valid for pipelines of non-circular flow area.

Pipeline system constant

The piping system constant CS is usually considered as a constant for a given opening of individual valves, but since the friction coefficient λ is a function of the Reynolds number, CS must also change with the flow rate. However, this change is not very large if we are interested in the pressure loss in the nominal flow region. For calculations in a larger range of flow rates, a correction can be applied by not multiplying the volumetric flow rate by 2, but by another exponent, see more in [Bašta, 2003, p. 25].

The pipeline system constant can be calculated according to Equation 16 from the individual pressure losses of the pipeline system for a known (nominal) flow rate (see Problem 1) or it can be calculated from the measured pressure loss at a particular volumetric flow rate, see Problem 2.

 Chapter: Determination of characteristics pipeline from measurements
1.14

Determination of characteristics pipeline from measurements

Pressure loss measurement

Engineering mathematics

The characteristic pipeline or fuction Lp=f(Q) can be established by measuring different cases. This data can then be processed using computer software or determined by plotting on logarithmic paper, creating a straight line whose slope corresponds to the flow rate exponent, as seen in Problem 2 and the article Engineering mathematics [Škorpík, 2023].

Change in pressure loss due to pipe fouling or corrosion

Pressure loss vs. fouling thickness

A fouling can form in the pipe if the liquid is not clean. Fouling in a pipeline system causes a reduction in the flow area of the pipeline and therefore a change in the characteristics pipeline and an increase in the pressure loss. Figure 17 shows the change in pressure loss in a pipeline when there is a uniform fouling in the pipeline - approximately the same percentage increase in pressure loss will increase the pumping work. The relationship in this figure was developed by substituting the Darcy-Weisbach equation into the pressure loss ratio Lp after reducing the flow area and the pressure loss Lp,n. From here it can be seen the reduction in diameter per pressure loss increases with the fifth power. On the other hand, even when absolute roughness is maintained, the effect of the change in friction coefficient is several orders of magnitude smaller.

Change of pressure loss of pipe due to fouling
17: Change of pressure loss of pipe due to fouling
Created for dn=100 mm; Vn=3 m·s-1; εn=0,05 mm; υn=553,2 nm2·s-1 (water at 50 °C); Q=const. F-fouling. The index n indicates the parameters before fouling.
 Chapter: Change in pressure loss due to pipe fouling or corrosion
1.15

Crystallization of minerals

Biofouling

Particles fouling

Fouling of the pipeline may be caused by chemical or biological action or by solid particles in the liquid. In the case of a chemical or electrochemical process, minerals precipitate and crystallise on the internal surfaces of the pipe. The biofouling on the pipe can be of plant or animal origin - usually some kind of algae or crustacean - and are highly dependent on water temperature, nutrient content of the water and, in the case of algae, light conditions. A typical sign of fouling by solid particles in the liquid is that it is not evenly distributed along the length of the pipe. The solid particles are deposited in areas of low flow velocity, at the lowest points of the pipeline route where the fluid flow is unable to displace them and upstream of constrictions.

Flow velocity vs. pipe fouling

Scale deposition on the pipe walls does not occur at velocities of approximately 1,5 to 2,5 m·s-1 [Vosmík, 2023]. However, at certain combinations of pH and temperature, this velocity may not be sufficient. Deposition of solid particles can be prevented from velocities as low as around 1,5 m·s-1, but also depends on the orientation of the pipe and the size and mass of the individual particles, according to [Pugh et al., 2009]. Biofouling of pipes can be prevented at velocities above 2 m·s-1.

Tangential tension

Viscosity

These velocities are for water. For other liquids, the limiting velocity may vary because a certain tangential tension, which is a function of viscosity, is required to prevent fouling at lower velocities and vice versa. Details on fouling of pipes and exchangers are given in the referenced paper [Pugh et al., 2009].

Prevence

Constant flow velocity during irregular pipeline operation can be maintained by creating loops on the exposed parts of the pipeline in which the fluid will flow at a constant velocity regardless of the flow rate between the inlet and outlet of the pipeline. Alternatively, a significantly higher nominal flow rate must be provided at switch-on to clear the pipe after partial operation (e.g. after a night of light operation). Prevention of pipe fouling can also be done by modifying the working fluid, filters or changing temperatures, but this is beyond the scope of this article.

 Chapter: Change in pressure loss due to pipe fouling or corrosion
1.16

Calculation of pipe fouling

Calculation of heat exchanger fouling

Cleaning

Shutdown

Rayleigh distribution

Fouling of pipes and heat exchangers usually gradually causes such problems that they need to be cleaned. The period when the pipes will need to be cleaned, that is, the shutdown time, can be predicted using statistics. This statistical method is based on the assumption that the increase in pressure drop follows a Rayleigh distribution, see Figure 18. To predict the increase in pressure loss due to pipeline fouling, it is sufficient to know an estimate of the operating time after which the pressure loss begins to increase, the expected modus of the rate at which the pressure loss increases most rapidly, and also the rate of increase in pressure loss at the start of fouling, see Problem 3. These estimates can be refined in real operation by measuring the pressure loss and thus refining the prediction of the increase in pressure loss over time.

Rayleigh distribution applied to pressure loss change
18: Rayleigh distribution applied to pressure loss change
n [s-1] change in pressure loss over time; t [s] time. The horizontal axis denotes the difference (V-Vn) because the Rayleigh distribution starts at zero and deposits form only after some time when the flow velocity is nominal Vn. The index n indicates the parameters before fouling.

Corrosion

Roughness

Pipe corrosion increases the absolute roughness of the pipe and causes a loss of pipe wall thickness. If the material loss does not cause a significant change in the flow surface area of the pipe, then, given the other parameters in the Darcy-Weisbach equation, the ratio of the pressure loss Lp to the pressure loss at nominal (initial) Lp,n can be expressed as a ratio of the coefficients of friction. The data in Table 7 shows that corrosion can increase the pressure loss by tens of percent. Therefore, when calculating the pipe that will not be cleaned of corrosion, the pressure loss must be calculated as if the pipe were corroded.

 Chapter: Pressure loss at significant density change
1.17

Pressure loss at significant density change

General equation

Critical gas velocity

In addition to fluid transport, we encounter dynamic gas flow in which the density of the gas can change significantly. If it is an adiabatic flow, then the pressure loss can be determined by assuming that the stagnation enthalpy of the gas remains constant and equal to the stagnation enthalpy at the inlet, but the entropy will increase due to internal friction. Based on this assumption, general Equations 19, which describes gas flow in the presence of friction in all types of channels, can be derived from the continuity, energy balance and momentum conservation equations for the assumption of constant specific heat capacity of the gas. However, in engineering practice, these equations are only used in calculations of flows with large density variations in narrow channel seals.

19: General equations of adiabatic gas flow in presence of friction
V*i [m·s-1] critical velocity for case of isentropic flow; κ [1] heat capacity ratio; A [m2] flow area of the channel; V [m·s-1] velocity of gas in investigated point of the channel (this velocity corresponds to velocity during isentropic expansion from stagnation pressure ps to static pressure p). If channel is not circular, characteristic leght L is used instead of d as in incompressible flow. Derivation in [Dejč, 1967, p. 209].

Friction coefficient

The friction coefficient λ in Equation 19 is a constant along the length of the channel, but in actual fact is more or less dependent on Re and the Mach number at the channel location under investigation. Thus, it depends on how much the channel flow area and Mach number vary. Experimental verification of the changes in the friction coefficient during compressible flow and the validity of Equation 19 is carried out in [Dejč, 1967, p. 217].

Special equation

In the case of compressible adiabatic flow in a channel of constant flow area, the pressure loss can be calculated using Equation 20, which is based on a modification of the general Equation 19 for the condition dA=0.

 Chapter: Pressure loss at significant density change
1.18
20: Equations for calculating the pressure loss when gas flows through channel with constant flow area
(a) velocity equation; (b) pressure loss equation; (c) continuity equation. Equations (a) and (b) are derived from Equation 19 for dA=0, the other assumptions of derivation are identical. Equation (c) is based on the continuity equation, where G=const.

Acceleration of flow

Fanno lines

In adiabatic gas flow, friction heats the gas, increasing its specific volume and velocity in a constant flow area channel. This leads to a gradual decrease in gas pressure and specific enthalpy. The Fanno line on the h-s diagram plots gas states along the channel axis. Figure 21 depicts three Fanno lines for a channel of length l with varying friction coefficients λ, influencing pressure changes as the channel lengthens (the same effect as changes in the friction coefficient has on the pressure change as the channel lengthens). At the maximum λ1, outlet flow doesn't reach critical velocity; λ2 just reaches critical velocity, and λ3, less than λ2, also reaches critical velocity at the outlet.

21: Fanno lines
h [J·kg-1] enthalpy; s [J·kg-1·K-1] entropy; hs [J·kg-1] stagnation gas enthalpy; h* [J·kg-1] critical enthalpy; psur [Pa] surrounding pressure at outlet of channel. The subscript i denotes the initial gas state, the subscript e the final gas state (at the end of the section/process under study). The subscript s denotes the stagnation gas state.
 Chapter: Pressure loss at significant density change
1.19

Seals

Seal friction coefficient

In engineering practice, the theory is particularly applicable to the investigation of flow in non-contact seals. The principle of dry-running gas seals is also based on the high pressure loss associated with gas flow in a very small gap. However, even labyrinth seals can be likened to a smooth seal with a constant flow area and a particular coefficient of friction.

Problems

Problem 1:

Pipeline system constant

Pipeline system constant

Pressure loss

Find the characteristics pipeline at the discharge of a condensate pump (see attached figure) in which condensate is pumped from the auxiliary condensate tank CT1 to the feed tank through the condensate heater H1. A parallel pipeline system with a redundant pump (blue) is connected to the route. The water temperature at the outlet of the pump is 60 °C and 105 °C after the H1 heater. The flow rate through the pump is 2,4 m3·h-1. The flow coefficient of ball valve 001 is 48,5 m3·h-1. The check valve has a pressure loss of 5 kPa. The minimum pressure loss of the balancing valve is 750 Pa. The pressure loss of the water meter is 18 kPa. The pressure loss of heater H1 is 12 kPa. The piping is standard one-inch water main. The solution to the problem is shown in Appendix 1.
CT1-auxiliary condensate tank No. 1; H1-heater No. 1; WM1-water meter No. 1. Marking according to [Krbek et al., 1999, p. 178]. The lengths of the individual sections of the piping system are given in metres.
§1   entry:   ti; te; Qn; KVS,001; Lp,002; Lp,003; Lp,WM1; Lp,H1; l §6   calculation:   ζpipe; Lpipe
§2   read off:   ν; ρ §7   calculation:   Lp,001
§3   read off:   d; ε §8   calculation:   Lp,elbow
§4   calculation:   ; Re §9   calculation:   Lp, n; CS
§5   calculation:   λ    
The procedure for solving Problem 1. Symbol descriptions are in Appendix 1.
Problem 2:

Pipeline system constant

Find the approximate value of the constant of the heating piping system. Hot water flows through the pipe. There are the measured flows through the system and the corresponding pressure loss given in the table below. Measured values adapted from [Pleskot, 1947, p. 17]. The solution to the problem is shown in Appendix 2.
 Chapter: Problems
1.20
Lp 10 25,1 62 140 320 700 1400
Q 19,64 29,64 50,07 74,61 113,9 161 233,7
Table of measured values
Hodnoty naměřené tlakové ztráty v potrubí zakreslené na logaritmickém papíře
Lp [Pa], Q [m3·s-1].
§1   entry:   Lp; Q §4   read off:   x; P
§2   marking:   Lp; Q do grafu §5   calculation:   CS
§3   approximation:   hodnot Lp; Q přímkou    
The procedure for solving Problem 2. Symbol descriptions are in Appendix 2.
Problem 3:

Exchanger fouling

Calculate the expected increase in pressure loss of the plate water/water exchanger using the statistical method. Scale crystallizes in the exchanger. The nominal flow velocity in the exchanger is 1 m·s-1 and the nominal pressure loss is 0.185 bar. Based on experience with the operation of previous exchangers, the pressure loss starts to increase after 500 minutes with an initial rate of 0.2703·10-3 min-1, and the parameters of the (V-Vn)-n curve in Figure 18 are: nmax=2.1622·10-3 min-1; (V-Vn)mod=1.1911 m·s-1. During operation, the flow rate remains constant. The solution to the problem is shown in Appendix 3.
Increased pressure loss in heat exchanger due to clogging
Lp [bar]; t [min] time.
 Chapter: Problems
1.21
§1   entry:   Vn; Lp,n; t0; n0; nmax;
(V-Vn)mod
  calculation:   Δt
§2   calculation:   C §4   calculation:   (Lp/Lp,n)k=1; Vk=1; Lp, k=1
§3   proposal:   tmax; kmax §5   calculation:   Vk; nk; (Lp/Lp,n)k; Lp, k
The procedure for solving Problem 3. Symbol descriptions are in Appendix 3.

References

ŠKORPÍK, Jiří, 2022a, Aerodynamika profilových mříží, Transformační technologie, Brno, ISSN 1804-8293, https://turbomachinery.education/aerodynamika-profilovych-mrizi.html.
ŠKORPÍK, Jiří, 2022b, Essential equations of turbomachines, Transformační technologie, Brno, ISSN 1804-8293, https://turbomachinery.education/essential-equations-of-turbomachines.html.
ŠKORPÍK, Jiří, 2023, Technická matematika, Transformační technologie, Brno, [online], ISSN 1804-8293. Dostupné z https://engineering-sciences.education/technicka-matematika.html.
ŠKORPÍK, Jiří, 2024, Introduction to turbomachinery, Transformační technologie, Brno, ISSN 1804-8293, https://turbomachinery.education/introduction-to-turbomachinery.html.
BAŠTA, Jiří, 2003, Hydraulika otopných soustav, Vydavatelství ČVUT, Praha, ISBN 80-01-02808-9.
BAUER, František, Oldřich BRŮHA a Zbyněk JAŇOUR, PEŠEK, Rudolf, ed., 1950, Základy proudění, Vědecko-technické nakladatelství, Praha.
CIHELKA, Jaromír, BRANDA, Jaroslav, CIKHART, Jiří, ČERMÁK, Jan, CHYSKÝ, Jaroslav, PITTER, Jaroslav, VALÁŠEK, Jiří, 1975, Vytápění a větrání, SNTL, Praha.
DEJČ, Michail, 1967, Technická dynamika plynů, SNTL, Praha.
FRAAS, Arthur, 1989, Heat exchanger design, John Wiley&Sons, Inc., ISBN 0-471-62868-9.
IBLER, Zbyněk, KARTÁK, Jan, MERTLOVÁ, Jiřina, IBLER, Zbyněk ml., 2002, Technický průvodce energetika, BEN-technická literatura, Praha, ISBN 80-7300-026-1.
IZARD, Julien, 1961, Příručka technické fyziky, Státní nakladatelství technické literatury, Praha.
JÍCHA, Miroslav, 2001, Přenos tepla a látky, Vysoké učení technické v Brně, Brno, ISBN 80-214-2029-4.
KRBEK, Jaroslav, POLESNÝ, Bohumil, FIEDLER, Jan, 1999, Strojní zařízení tepelných centrál-Návrh a výpočet, PC-DIR Real, s.r.o., Brno, ISBN 80-214-1334-4.
MAŠTOVSKÝ, Otakar, 1964, Hydromechanika, Statní nakladatelství technické literatury, Praha.
MÍKA, Vladimír, 1977, Základy chemického inženýrství, Statní nakladatelství technické literatury, Praha.
MIKULA, Julius, KOČKA, Jaroslav, ŠKRAMLÍK, Emanuel, ŠTAUBER, Zdeněk, VESELÝ Adolf, OBR, Jan, 1974, Potrubí a armatury, Státní nakladatelství technické literatury, Praha.
MILLER, Rudolf, HOCHRAINER, A., LÖHNER, K., PETERMANN, H., 1972, Energietechnik und Kraftmaschinen, Rowohlt taschenbuch verlag GmbH, Hamburg, ISBN 3-499-19042-7.
NIKURADSE, Johann, 1933, Strömungsgesetze in rauen Rohren, V. D. I. Forschungsheft, 361: 1–22, Berlin.
PLESKOT, Václav, 1947, Nomografie v technické praxi, Praha, SPASEI.
POLESNÝ, Bohumil a kol., 1990, Termodynamická data pro výpočet tepelných a jaderných energetických zařízení, Vysoké učení technické v Brně, Brno, ISBN 80-214-0160-5.
PUGH, Simon, HEWITT, Geoffrey, MÜLLER-STEINHAGEN, Hans, 2009, Fouling During the Use of “Fresh” Water as Coolant—The Development of a “User Guide”, Heat Transfer Engineering, 30:10-11, 851-858, DOI: 10.1080/01457630902753706.
ROČEK, Jaroslav, 2002, Průmyslové armatury, INFORMATORIUM, Praha, ISBN 80-7333-000-8.
ŘASA, Jaroslav, ŠVERCL, Josef, 2004, Strojnické tabulky, jednotky, matematika, mechanika, technické kreslení, strojní součásti, Scientia, spol. s.r.o. Praha, ISBN 80-7183-312-6.
ŠTEFAN, David, 2009, Hydraulické ztráty v potrubí, Vysoké učení technické v Brně, Brno, Bakalářská práce.
VOSMÍK, Zdenek, 2023, Nerezové výměníky tepla – koroze a zanášení úsadami 2. část, https://www.vosmik-vymeniky.cz/, https://www.vosmik-vymeniky.cz/novinky.html.
©Jiří Škorpík, LICENCE