The implementation of the function is based on "Handbook of Hydraulic Resistance" in its first translated Version from 1960! The book has been republished in several updated versions since then!
Function calculating the pressure loss of a Y-shaped splitter of type I (Fig. 1) as f(w_c,w_b,w_s, alpha, rho). (Currently not yet available)
where:
Calculation according to Idelchik (1960). The pressure loss is calculated as
dp_b = rho/2 * zeta_cb * w_c^2
dp_s = rho/2 * zeta_cs * w_c^2
As you can see above both pressure loss calculation are with respect to the velocity at the common branch. The pressure loss coefficient of the branch zeta_cb is calculated as: zeta_cb = A' * zeta'_cb
According to Idel'chik while w_relbc = w_b/w_c <= 0.8, A'= 1
and if w_relbc > 0.8
then A' = 0.8
.
Since the velocity ratio may change during simulations, the discontinuity in A' is circumvented by the stepsmoother function used in region 0.8 < w_relbc < 1.05
.
The uncorrected pressure loss coefficient zeta'_cb
is determined by the formula:
zeta'_cb = 1 + w_relbc^2 - 2*cos(alpha)
The pressure loss coefficient at the straight channel for all brnaching angles is described as:
zeta_cs = 0.4*(1-w_relsc)^2
The following figure Fig.2, pressure loss coefficients of the branching channel are shown. (Currently not yet available)
[P. Jordan; HTWG Konstanz; 01/24]
function dp_SplitterWyeType1_DP extends Modelica.Icons.Function; import SMOOTH = ThermofluidStream.Processes.Pipes.Internal.Utilities.Stepsmoother; input SI.Velocity w_c "Velocity of common inlet"; input SI.Velocity w_b "Velocity of branching outlet"; input SI.Velocity w_s "Velocity of straight outlet"; input SI.Angle alpha "Branching angle"; input SI.Density rho "Medium density"; input SI.Velocity eps "To avoid division by zero"; output SI.Pressure dp_b "Pressure loss of branching outlet"; output SI.Pressure dp_s "Pressure loss of straight outlet"; output Internal.Types.PressureLossCoefficient zeta_cb "Pressure loss coefficient of branching pipe w.r.t. common inlet velocity"; output Internal.Types.PressureLossCoefficient zeta_cs "Pressure loss coefficient of straight pipe w.r.t. common inlet velocity"; end dp_SplitterWyeType1_DP;