Water Tee
Tee junction for incompressible flow
Description
Equations
Variables
Connections
Parameters
See Also
The Water Tee component models a tee junction for incompressible flow, and has multiple methods to calculate a local resistance for each path, incl. Idelchik[1], Rennels[2], and Blevins[3]. In addition, the calculation of these local resistances can be specified as use of simulated (Time-variant) flow rates or based on specified flow rate at design point (Time-invariant). This component calculates mainly pressure difference and mass flow rate.
Important note : This component supports only "Dynamics of mass = Dynamic".
Basic concept
The following diagram shows the relationship between parameters and the geometry. Noted that port_a and port_b are the same diameter.
This component consists of three flow calculation components and an internal control volume at the junction. The local resistances are calculated based on the direction of flows, that are determined with pressure difference between the junction (control volume) and ports.
Regarding the equations of Control volume, you can find more information in the Water Volume help page. The each flow component calculate the mass flow rate with the following equation:
mflow=α⋅A⋅{inStream`port_a.rho`dp≥0inStream`port_b.rho`others⋅dp
A is cross section area, α is total flow coefficient. The total flow coefficient is calculated with:
α=α__friction+α__local=2⋅D__hλ⋅L+2ζ
D__h is diameter of the section, L is the length of the section, and λ is Friction coefficient. ζ is the local resistance coefficient.
Friction coefficient calculation
The following table shows the coverage of flow directions for each method.
Note : dp[1] means the pressure difference between port_a and junction. dp[2] is for between port_b and junction. dp[3] is also for port_c and junction.
Type of flow
Condition
(pressure difference)
Availability
Comment
Converging 1 : port_a and port_c to port_b
port_a - junction node
dp1>0
port_b - junction node
dp2<0
port_c - junction node
dp3>0
Idelchik : available
Rennels : available
Blevins : available
Converging 2: port_b and port_c to port_a
dp1<0
dp2>0
Converging 2: port_a and port_b to port_c
dp3<0
Idelchik : ---
For Idelchik, parameter zeta__con_in is used.
Diverging 1: port_a to port_b and port_c
Diverging 2: port_b to port_a and port_c
Diverging 3: port_c to port_a and port_b
For Idelchik, parameter zeta__div_ex is used.
Constant zeta
This is applied when the parameter setting is Tee calculation = Constant zeta.
The local resistance ζ is defined by a parameter ζ__constant.
ζ1=ζ__constant1
ζ2=ζ__constant2
ζ3=ζ__constant3
Note : [1] port_a, [2] port_b, [3] port_c (branch)
Idelchik's method
This is applied when the parameter setting is Tee calculation = Idelchik with design point, or Idelchik with actual flow ratio.
The following equations are defined in Idelchik's method for each case. The calculation uses Volume flow rate vflow.
If Tee calculation = Idelchik with design point, these volume flow rates are specified by parameter values,
vflow1=vflow__dp1
vflow2=vflow__dp2
vflow3=vflow__dp3
And, if Idelchik with actual flow ratio, the actual calculated volume flow rates are used for vflow.
ζ1=1.55⋅vflow3vflow2−vflow3vflow221−vflow3vflow22
ζ2=0
ζ3=({1.0A__sA__m≤0.35{0.9⋅1−vflow3vflow2vflow3vflow2≤0.40.55otherwiseotherwise)⋅1+vflow3vflow2⋅A__mA__s−2⋅1−vflow3vflow22vflow3vflow2⋅A__mA__s2
Converging 1 : port_b and port_c to port_a
ζ1=0
ζ2=1.55⋅vflow3vflow1−vflow3vflow121−vflow3vflow12
ζ3=({1.0A__sA__m≤0.35{0.9⋅1−vflow3vflow1vflow3vflow1≤0.40.55otherwiseotherwise)⋅1+vflow3vflow1⋅A__mA__s−2⋅1−vflow3vflow12vflow3vflow1⋅A__mA__s2
Converging 3 : port_a and port_b to port_c
ζ1=zeta__con_in
ζ2=zeta__con_in
ζ3=0
Diverging 1 : port_a to port_b to port_c
ζ2=vflow3vflow12⋅({0.4A__sA__m≤0.4{2⋅2⋅vflow3vflow1−1vflow3vflow1≤0.50.3⋅2⋅vflow3vflow1−1otherwiseotherwise)1−vflow3vflow12
ζ3=1+0.3⋅vflow3vflow1⋅A__mA__s2vflow3vflow1⋅A__mA__s2
Diverging 2 : port_b to port_a to port_c
ζ1=vflow3vflow22⋅({0.4A__sA__m≤0.4{2⋅2⋅vflow3vflow2−1vflow3vflow2≤0.50.3⋅2⋅vflow3vflow2−1otherwiseotherwise)1−vflow3vflow22
ζ3=1+0.3⋅vflow3vflow2⋅A__mA__s2vflow3vflow2⋅A__mA__s2
Diverging 3 : port_c to port_a to port_b
ζ1=ζ__div_ex
ζ2=ζ__div_ex
Other : pressure difference ports and junction are close to zero [Pa] (dp≤dp__small2)
ζ1=ζ__0
ζ2=ζ__0
ζ3=ζ__0
Rennels' method
This is applied when the parameter setting is Tee calculation = Rennels with design point, or Rennels with actual flow ratio.
The following equations are defined in Idelchik's method for each case. The calculation uses Mass flow rate mflow.
If Tee calculation = Rennels with design point, these mass flow rates are specified by parameter values
mflow1=mflow__dp1
mflow2=mflow__dp2
mflow3=mflow__dp3
And, if Rennels with actual flow ratio, the actual calculated mass flow rates are used for mflow.
ζ1=0.54−1.12⋅rD__side+0.28⋅rD__side⋅mflow1mflow2−2+0.38+0.42⋅rD__side+0.56⋅rD__side⋅mflow1mflow2−1−0.88+0.7⋅rD__side−0.84⋅rD__side
ζ3=−1.92+1.4⋅rD__side−0.84⋅rD__side+3.46−2.7⋅rD__side+1.12⋅rD__side⋅mflow3mflow2−1+−0.92+0.2⋅rD__side+0.07⋅rD__side⋅mflow3mflow2−2⋅D__sideD__main4+1.0−0.5⋅D__sideD__main1+D__side/D__main
ζ2=0.54−1.12⋅rD__side+0.28⋅rD__side⋅mflow2mflow1−2+0.38+0.42⋅rD__side+0.56⋅rD__side⋅mflow2mflow1−1−0.88+0.7⋅rD__side−0.84⋅rD__side
ζ3=−1.92+1.4⋅rD__side−0.84⋅rD__side+3.46−2.7⋅rD__side+1.12⋅rD__side⋅mflow3mflow1−1+−0.92+0.2⋅rD__side+0.07⋅rD__side⋅mflow3mflow1−2⋅D__sideD__main4+1.0−0.5⋅D__sideD__main1+D__side/D__main
ζ1=0.81−1.16⋅rD__side+0.5⋅rD__side⋅mflow1mflow3−2−0.95−1.65⋅rD__side⋅mflow1mflow3−1+1.34−1.69⋅rD__side
ζ2=0.81−1.16⋅rD__side+0.5⋅rD__side⋅mflow2mflow3−2−0.95−1.65⋅rD__side⋅mflow2mflow3−1+1.34−1.69⋅rD__side
ζ2=0.62−0.98⋅mflow2mflow1−1+0.36⋅mflow2mflow1−2+0.04⋅mflow2mflow16
ζ3=0.81−1.13−0.16⋅rD__side⋅mflow3mflow1−1+1.00−0.24⋅rD__side⋅mflow3mflow1−2⋅D__sideD__main4+1.08⋅D__sideD__main−1.06⋅D__sideD__main3+zeta__entrance
zeta__entrance=0.57−1.07⋅rD__side−2.13⋅rD__side+8.24⋅rD__side3/2−8.48⋅rD__side2+2.90⋅rD__side5/2
ζ1=0.62−0.98⋅mflow1mflow2−1+0.36⋅mflow1mflow2−2+0.04⋅mflow1mflow26
ζ3=0.81−1.13−0.16⋅rD__side⋅mflow3mflow2−1+1.00−0.24⋅rD__side⋅mflow3mflow2−2⋅D__sideD__main4+1.08⋅D__sideD__main−1.06⋅D__sideD__main3+ζ__entrance
ζ__entrance=0.57−1.07⋅rD__side−2.13⋅rD__side+8.24⋅rD__side3/2−8.48⋅rD__side2+2.90⋅rD__side5/2
ζ1=0.59⋅mflow1mflow3−2+1.18−1.84⋅rD__side+1.16⋅rD__side⋅mflow1mflow3−1−0.68+1.04⋅rD__side−1.16⋅rD__side
ζ2=0.59⋅mflow2mflow3−2+1.18−1.84⋅rD__side+1.16⋅rD__side⋅mflow2mflow3−1−0.68+1.04⋅rD__side−1.16⋅rD__side
Blevins' method
This is applied when the parameter setting is Tee calculation = Blevins with design point, or Blevins with actual flow ratio.
The following equations are defined in Idelchik's method for each case. The calculation uses flow velocity v.
If Tee calculation = Blevins with design point, these flow velocities are specified by parameter values
v1=v__dp1
v2=v__dp2
v3=v__dp3
And, if Blevins with actual flow ratio, the actual calculated flow velocity are used for v.
ζ1=0.045+1.38−1.94⋅rD__side+1.34⋅rD__side⋅v3v2−0.90−0.95⋅rD__side+1.23⋅rD__side⋅v3v221−v3v2⋅A__sideA__main2
ζ3=1.09−0.8⋅rD__side−0.53+1.27⋅rD__side−1.86⋅rD__side⋅v1v2−1.48−2.28⋅rD__side+1.80⋅rD__side⋅v1v22v3v2
ζ2=0.045+1.38−1.94⋅rD__side+1.34⋅rD__side⋅v3v1−0.90−0.95⋅rD__side+1.23⋅rD__side⋅v3v121−v3v1⋅A__sideA__main2
ζ3=1.09−0.8⋅rD__side−0.53+1.27⋅rD__side−1.86⋅rD__side⋅v2v1−1.48−2.28⋅rD__side+1.80⋅rD__side⋅v2v12v3v1
ζ1=1.19−1.16⋅rD__side+0.46⋅rD__side−1.73⋅1.0−rD__side⋅v1v3+1.34−1.69⋅rD__side⋅v1v32v1v32
ζ2=1.19−1.16⋅rD__side+0.46⋅rD__side−1.73⋅1.0−rD__side⋅v2v3+1.34−1.69⋅rD__side⋅v2v32v2v32
ζ2=({1.55⋅0.22−v3v12−0.03v3v1<0.220.65⋅v3v1−0.222−0.03otherwise)⋅11−v3v12
ζ3=0.99−0.23⋅rD__side−0.82+0.29⋅rD__side+0.3⋅rD__side⋅v3v1+1.02−0.64⋅rD__side+0.76⋅rD__side⋅v3v12v3v12
ζ1=({1.55⋅0.22−v3v22−0.03v3v2<0.220.65⋅v3v2−0.222−0.03otherwise)⋅11−v3v22
ζ3=0.99−0.23⋅rD__side−0.82+0.29⋅rD__side+0.3⋅rD__side⋅v3v2+1.02−0.64⋅rD__side+0.76⋅rD__side⋅v3v22v3v22
ζ1=0.59+1.18−1.84⋅rD__side+1.16⋅rD__side⋅v1v3−0.68−1.04⋅rD__side+1.16⋅rD__side⋅v1v32v1v32
ζ2=0.59+1.18−1.84⋅rD__side+1.16⋅rD__side⋅v2v3−0.68−1.04⋅rD__side+1.16⋅rD__side⋅v2v32v2v32
Friction calculation
The calculation method can be specified with 3 options.
Friction calculation = Constant
If this option is selected, the friction coefficient λ is specified by a parameter λ__constant.
λ=λ__constant
Friction calculation = Darcy-Weisbach with Constant velocity (Design point)
If this option is selected, the friction coefficient λ is calculated with a parameter of velocity at design point v__design_f.
Re=max{ρ__adp≥0ρ__bothers⋅v__design_f⋅D__h_act{μ__adp≥0μ__bothers,0.1
λ=`HeatTransfer.Functions.lambda_Re`Re,roughness,D__h_act,Re__CoT,IF__speed,Geo__act
ρ__a, ρ__b are density at each port of flow components. μ__a, μ__b are dynamic viscosity at each port of flow components. This means, one is at port of Water Tee junction, and the another is at the junction node. D__h_act is the diameter of flow components. Re is Reynolds number. roughness is pipe roughness of flow components.
Friction calculation = Darcy-Weisbach
If this option is selected, the friction coefficient λ is calculated with the following equations:
Re__target=max{ρ__adp≥0ρ__bothers⋅v⋅D__h_act{μ__adp≥0μ__bothers,0.1
ⅆReⅆt=Re__target−ReT__const
ρ__a, ρ__b are density at each port of flow components. μ__a, μ__b are dynamic viscosity at each port of flow components. D__h_act is the diameter of flow components. Re is Reynolds number. roughness is pipe roughness of flow components.
(*) The above function `HeatTransfer.Functions.lambda_Re` is to calculated friction factor for Laminar and Turbulent flow. The fundamental implementation is based on the following equations. Especially, the equation of Turbulent flow is Swamee and Jain's approximation[1] .
You can find more information in the Water Detailed Flow help page.
References
[1] : I.E. Idelchik, "Handbook of Hydraulic Resistance", 4th edition, begell house, inc.
[2] : D.C. Rennels (2022), "Pipe Flow: A Practical and Comprehensive Guide", 2nd Edition, Wiley
[3] : R.D. Blevins (1984), "Applied Fluid Dynamics Handbook", Van Nostrand Reinhold Co.
Symbol
Units
Modelica ID
dp3
Pa
Pressure difference
[1] port_a, [2] port_b, [3] port_c
dp
mflow3
kgs
Mass flow rate
mflow
v3
ms
Velocity of flow
v
vflow3
m3s
Volume flow rate
vflow
ρ__act3
kgm3
Density
rho_act
μ__act3
Pa⋅s
Dynamic viscosity
vis_act
Re3
−
Reynolds number
Re
ζ3
Local flow resistance
zeta
α3
Flow coefficient (sqrt(2*Dh/lambda/L) + sqrt(2/zeta))
alpha
λ3
Friction coefficient for Darcy-Weisbach equation
lambda
α__local3
Flow coefficient for local resistance
alpha_local
α__friction3
Flow coefficient for Friction
alpha_friction
vflow3vflow1
Ratio of volume flow rate port_c/port_a
vflow_ratio_a
vflow3vflow2
Ratio of volume flow rate port_c/port_b
vflow_ratio_b
mflow1mflow2
Ratio of volume flow rate port_a/port_b
mflow_ratio_a_b
mflow2mflow1
Ratio of volume flow rate port_b/port_a
mflow_ratio_b_a
mflow2mflow3
Ratio of volume flow rate port_b/port_c
mflow_ratio_b_c
mflow1mflow3
Ratio of volume flow rate port_a/port_c
mflow_ratio_a_c
mflow3mflow1
mflow_ratio_c_a
mflow3mflow2
mflow_ratio_c_b
v1v2
Ratio of flow velocity rate port_a/port_b
v_ratio_a_b
v2v1
Ratio of flow velocity port_b/port_a
v_ratio_b_a
v2v3
Ratio of flow velocity port_b/port_c
v_ratio_b_c
v1v3
Ratio of flow velocity port_a/port_c
v_ratio_a_c
v3v1
Ratio of flow velocity port_c/port_a
v_ratio_c_a
v3v2
Ratio of flow velocity port_c/port_b
v_ratio_c_b
Name
port__a
Water Port
port_a
port__b
port_b
port__c
port_c
Default
Watersimulationsettings
WaterSettings1
Specify a component of Water simulation settings
Settings
Teecalculation
Select calculation mode of Tee junction
- Constant zeta
- Idelchik with design point
- Idelchik with actual flow ratio
- Rennels with design point
- Rennels with actual flow ratio
- Blevins with design point
- Blevins with actual flow ratio
modelTee
Frictioncalculation
Constant
Select calculation mode of friction loss
- Constant
- Darcy-Weisbach with Constant velocity (Design point)
- Darcy-Weisbach
modelFriction
D__main
0.1
m
Inner diameter of pipe, Main channel
Dh_m
A__main
π⋅D__main22
m2
Cross section area, Main channel
A_m
D__side
Inner diameter of pipe, Side channel
Dh_s
A__side
π⋅D__side22
Cross section area, Side channel
A_s
r
0.001
rad
Radius of curvature of Tee junction, up to r/Dh_s = 0.5
L3
0.1,0.1,0.1
Length of seach section, 1:port_a-junction 2:port_b-junction 3:port_c-junction
L
roughness
2.5ⅇ−5
Absolute roughness of pipe, with a default for a smooth steel pipe
zeta__constant3
1,1,1
Coeffiecient of Local resistance, 1:port_a 2:port_b 3:port_c
(When Tee calculation = Constant zeta)
zeta_constant
vflow__dp3
1.0,0.5,0.5
Design points in Volume Flow Rate, 1:port_a 2:port_b 3:port_c
(When Tee calculation = Idelchik with design point or Idelchik with actual flow ratio)
vflow_dp
mflow__dp3
Design points in Mass Flow Rate, 1:port_a 2:port_b 3:port_c
(When Tee calculation = Rennels with design point or Rennels with actual flow ratio)
mflow_dp
v__dp3
Design points in Flow velocity, 1:port_a 2:port_b 3:port_c
(When Tee calculation = Blevins with design point or Blevins with actual flow ratio)
v_dp
zeta__con_in
0.7
Loss coefficient for case of converging to branch, intake side, Default value is obtained from Rannels
zeta_branch_converging_in
ζ__div_ex
30
Loss coefficient for case of diverging to branch, exhaust side, Default value is obtained from Rannels
zeta_branch_diverging_ex
ζ__0
1
Loss coefficient, around zero pressure difference
zeta_zero
zeta__min
Loss coefficient, minimum value
zeta_min
λ__const3
0.15ⅇ−4, 0.15ⅇ−4, 0.15ⅇ−4
(When Friction calculation = constant)
lambda_const
v__dp_fric3
(When Friction calculation = Darcy-Weisbach with Constant velocity (Design point))
v_dp_fric
V
m3
Volume of control volume
p__start
101325
Initial condition of pressure, junction node
p_start
T__start
293.15
K
Initial condition of temperature, junction node
T_start
T__const
s
Time constant for Reynolds number changing
T_const
sharpness
1.0
Sharpness of approximation for sqrt(dp) and sqrt(rho * dp)
dp__small
Regularization of zero flow if |dp| < dp_small (dummy if use_dp_small = false)
dp_small
Re__CoT
3500
Reynolds number of the center of Transition zone
Re_CoT
IF__spread
0.007
Spread of Intermittency factor
IF_spread
Heat Transfer Library Overview
Water Overview
Water Shapes Overview
Download Help Document