Friction Correlations

The friction in the Flow Simulator element can be calculated by the given equation, where:

stands for the calculated friction based on the user-selected friction mode (Abauf, Swamee), friction type (Fanning or Darcy), and Re (to determine whether it is in a turbulent region, or in a laminar region).

If Re_Dh < Re_Turb, laminar friction calculations take place, otherwise, the turbulent friction calculation routine is used.


Nomenclature	Subscripts
f: friction	F: fanning
Re: Reynolds number	D: darcy
FMULT: friction multiplier	turb: turbulent flow
ε: sand grain roughness	lam: laminar flow
A: Cross sectional area	Abuaf: Abuaf friction relation
L₊: Inlet station + 1/9 of 2^nd station	Smooth: smooth surface
XMU: dynamic viscosity	Rough: rough surface
W: mass flow rate	SJ: Swamee-Jain approximation
X: station length	Dh: hydraulic diameter
	L: equivalent diameter

Laminar Friction

Calculates the friction coefficient for laminar flow in shaped ducts based on the references Yunus A. Cengel, 2006 and Bruce Munson, 2005.

Friction coefficient for hydrodynamically fully developed flow can be calculated as:

For a Tube Element, Laminar Friction Inlets effects can be accounted. Friction coefficient for hydrodynamically developing flow with “Muzychka Yovanovich Laminar Inlet Effects” can be calculated as:

The friction coefficient for combining developing flow and fully developed flow can be calculated as:

Darcy type friction is calculated as

Fanning type friction is calculated as:

Turbulent Friction

Calculates the turbulent friction for smooth or rough walls.

Abuaf Friction Relation

The Abuaf friction relation should generally be used for smooth walled tubes.

In Flow Simulator, you have the option to use the Abuaf friction relation together with wall roughness. The following adjustment equation is used:

Swamee-Jain Approximation of the Colebrook-White Equation (Moody Diagram)

The Darcy and Fanning type frictions are calculated as:

User-specified Friction Factor

Roughness

Surface roughness values can be entered in four different measurement types. The roughness values are converted to sand grain roughness equivalents using the following equations from table 1 of reference 63.


ε=5.863∗R_a,	R_a=Average Absolute Roughness
ε=3.100∗R_rms,	R_rms=Root Mean Square Roughness
ε=0.978∗R_zd,	R_zd=Peak to Valley Roughness

Non-Circular Shapes in Flow Simulator Tubes

The friction factor and heat transfer coefficient (HTC) correlations were developed for circular pipes. The traditional method to use these correlations on non-circular shapes is to calculate a hydraulic diameter based on the shape area and perimeter.

D h = \frac{4 * A r e a}{P e r i m e t e r}

The errors associated with this method can be +/-40% for laminar flow but less for turbulent flow, +/-15% (see White, ref 3).

A more accurate option is to adjust the hydraulic diameter with a friction factor ratio (see White, ref 3). The effective hydraulic diameter can then be used in the friction factor and HTC correlations.

The following table summarizes the relationship between the Dh based on 4*A/P and the effective hydraulic diameter.


Shape	Effective Dh Equation	Aspect Ratio (AR)
Circle	$\frac{{D h}_{e f f}}{D h} = 1$	AR=1
Rectangle	$\frac{{D h}_{e f f}}{D h} = \frac{2}{3} + \frac{11}{24} * A R * (2 - A R)$	AR=b/a
Ellipse	$\frac{{D h}_{e f f}}{D h} = 1 - . 2109 * {(1 - A R)}^{2}$	AR based on area and perimeter of the ellipse.
Isosceles Triangle	$\frac{{D h}_{e f f}}{D h} = \frac{64}{48 + 11.442 * A R - 6.0026 * {A R}^{2}}$	AR based on area and perimeter of the triangle.
Annulus	$\frac{{D h}_{e f f}}{D h} = \frac{1 + {A R}^{2} + \frac{(1 - {A R}^{2})}{l n (A R)}}{{(1 - A R)}^{2}}$	AR=b/a
Freeform (Arbitrary- Shape)	$\frac{{D h}_{e f f}}{D h} = \frac{\sqrt{A R} * (1 + A R) * (1 - (\frac{192}{π^{5}}) * A R * \tanh (\frac{π}{2 * A R}))}{. 75 * \frac{D h}{\sqrt{A r e a}}}$	AR based on area and perimeter of freeform shape using a rectangle equation. $A R = \frac{(P - \sqrt{(P^{2} - 16 * A)})}{(P + \sqrt{(P^{2} - 16 * A)})}$

See Blevins (ref 15) and Muzychka et al. (ref 50) for additional information.

If the compressible tube, advanced orifice, and incompressible tube have a cross sectional shape that is not circular, the equations in this table are used for the effective hydraulic diameter equation.

Two-phase Flow Friction

The friction factor for two-phase flow (liquid and gas) in an incompressible tube can be calculated using two options. The homogenous approach uses the laminar and turbulent equations shown above with fluid properties based on the liquid/gas mixture. The second approach uses friction equations developed by Friedel (ref 64).

\emptyset = \frac{{(\frac{Δ P}{Δ L})}_{2 - p h a s e}}{{(\frac{Δ P}{Δ L})}_{l i q u i d o n l y}}

\emptyset = E + \frac{3.24 F H}{F r^{0.045} W e^{0.035}}

E = {(1 - x)}^{2} + x^{2} (\frac{ρ_{l} f_{g}}{ρ_{g} f_{l}})

F = x^{0.78} {(1 - x)}^{0.224}

H = {(\frac{ρ_{l}}{ρ_{g}})}^{0.91} {(\frac{μ_{g}}{μ_{l}})}^{0.19} {(1 - \frac{μ_{g}}{μ_{l}})}^{0.7}

F r = F r o u d e N u m b e r = \frac{U^{2}}{g L}

W e = W e b e r N u m b e r = \frac{ρ U^{2} D_{h}}{σ} (ρ i s d e n s i t y o f t h e m i x t u r e)

x = \frac{\dot{m_{g}}}{\dot{m_{g}} + \dot{m_{l}}}

(quality)

Recommended to use if $\frac{μ_{l}}{μ_{g}} < 1000$