When fluid flows through a pipe bend or elbow, the fluid changes direction and hence there is a change in in momentum. As per Newtons Second Law: "The rate of change of momentum of a body is equal to the net force acting on the body". Slug forces are generated when there is 'two-phase' (gaseous and liquid phase) fluid flowing through pipes and slugs of liquid that are formed intermittently travel at relatively high velocity.

The liquid slug can cause large reaction forces at changes in direction such as elbows or tee connections. High slug forces have the potential to cause large pipe deflections or damage to the supports if they are not designed to withstand them. When the velocity of slug flow is high, the slug force is considered to act at one location at any instant of time for the piping system under consideration. The slug load lasts until the slug has traversed the elbow and then drops to near zero. The time duration of the load is calcuated as slug length divided by the fluid velocity.

Force on Pipe Bend due to Fluid Flow

Let us assume a horizontal pipe bend turning through an angle θ as shown in figure below.

The following details apply:

  • V = Velocity of flow
  • Q = Flow Rate
  • ρ = Fluid Density
  • A = Cross-section area of pipe ID
  • θ = Change in direction at Bend

The forces exerted on the bend by the flowing fluid are Fx and Fy respectively.

Force in X-direction

Net force acting on the bend in x-direction (Fx)  = Rate of change of momentum in x-direction

= Mass Flow * (Final velocity in x-direction - Initial velocity in x-direction)

= Density(ρ) * Flow Rate(Q) * (Vcosθ - V) = (ρ) * (Q) * (Vcosθ - V)

Net force acting on the bend in x-direction (Fx) =  ρQV(cosθ - 1)...........(1)

Force in Y-direction

Net force acting on the bend in x-direction (Fy)  = Rate of change of momentum in y-direction

= Mass Flow * (Final velocity in y-direction - Initial velocity in y-direction)

= Density(ρ) * Flow Rate(Q) * (Vsinθ - 0) = (ρ) * (Q) * (Vsinθ)

Net force acting on the bend in y-direction (Fy) =  ρQ(Vsinθ)...........(2)

Resultant Force acting on the Bend

The resultant force acting on the bend (FR) = \(\sqrt[]{(F_x)^2 + (F_y)^2 }\)..........(3)

Solving equation (3) gives

The resultant force acting on the bend (FR) = ρQV\(\sqrt[]{2(1-cosθ)}\)..........(4)

Substituting Q = A*V

The resultant force acting on the bend (FR) = ρAV2\(\sqrt[]{2(1-cosθ)}\)..........(5)

For 90 elbow or bend, cosθ = 0, hence 

The resultant force acting on the bend (FR) = ρAV2\(\sqrt[]{2}\) .........(6)

For equivalent static analysis in Caesar-II, the slug load calculated by above formula shall be multiplied by a Dynamic Load Factor (DLF) of two, unless more accurate data is available.

Example Problem:

A slug of length 3.05m passes a 12" sch 40, 90° elbow with a velocity of 15m/s. The density of the fluid is 800 kg/m3. Calculate slug load acting on the elbow. Consider a DLF of 2.

Solution:

The slug load is calculated using equation (6)

Density (ρ) = 800kg/m3

ID of pipe = 323.8 - 2*10.31 = 303.18mm = 0.303m

Cross-section Area (A) = π * (0.303)2 / 4 = 0.07211 m2

Velocity (V) = 15 m/s

Length of Slug (L) = 3.05m

Resultant Force acting on the bend = (DLF) * ρAV2\(\sqrt[]{2}\) = 2 * 800 * 0.07211 * 152 * \(\sqrt[]{2}\) = 36707 N.

If R is the radius of elbow, then Rθ(in radian) is the length of bend. 

Elbow Radius (Long Radius) of 90° (R) = 457mm = 0.457m.

90° = π/2 = 1.57 radians, hence the length of bend (Rθ) = 0.457 * 1.57 = 0.7175m.

It takes the slug \(\frac{Rθ}{V}\) seconds to exit the bend and attain the maximum slug force of 36707 N.

Slug Force Time History

The slug force - time history is as indicated in graph above.

The slug force increases to maximum value in \(\frac{Rθ}{V}\) seconds = \(\frac{0.7175}{15}\) = 0.04783 sec.

The maximum slug force acts for \(\frac{L-Rθ}{V}\) seconds = \(\frac{(3.05-0.7175)}{15}\) = 0.1555 sec.

The slug force drops to zero in the next \(\frac{Rθ}{V}\) seconds = \(\frac{0.7175}{15}\) = 0.04783 sec.

Thus t1 = 0.04783s, t2 = 0.04783 + 0.1555 = 0.2033s and t3 = 0.2033 + 0.04783 = 0.2512s.

Other articles you may like to read...

Trunnion (Dummy) Support - Stress On...
Trunnion (Dummy) Support - Stress On... The dummy support calculation is based on M.W.Kellogg Design Manual. As per the Kellogg manual, all attachments to the p... Stress Analysis Related Calculations 8848 views
Allowable Sress Values for Piping Ma...
Allowable Sress Values for Piping Ma... Allowable stress values for piping materials are provided in Table A-1 of ASME B31.3. This page provides a convenient an... Stress Analysis Related Calculations 6194 views