Experimentally Determined Solute Mixing under Laminar and Transitional Flows at Junctions in Water Distribution Systems

Associate Professor, Department of Civil Engineering, Zhejiang University, Hangzhou 310058, China M.S. Student, Department of Civil Engineering, Zhejiang University, Hangzhou 310058, China Senior Scientist and Advisor, U.S.EPA, Office of Research and Development, ORISE Fellowship, 26 W. Martin Luther King Dr., Cincinnati, OH 45268, USA Professor, Department of Civil Engineering, Zhejiang University, Hangzhou 310058, China Associate Professor, Department of Civil Engineering, Zhejiang University, Hangzhou 310058, China


Introduction
Drinking water distribution network modeling using EPANET program [1] or EPANET Multi-Species Extension is widely used to predict the fate and transport of disinfectant residues and contaminants [2], the optimization of water quality monitoring sensors, optimal pump operation, and valve control for water quality purpose [3], and the determination of pollutant source.EPANET and those adopting the core simulation models assume perfect mixing of solutes in water parcels at pipe junctions [4].Later, using numerical and experimental studies, Austin et al. [5] and Choi et al. [6] found that the complete mixing assumption may underestimate substance concentration leading to a large simulation error.Similarly, Ho et al. [7], Choi [8], and Mckenna et al. [9] found that the degree of mixing was much lower than that of complete mixing at pipe cross junctions.When studying the effects of Reynolds numbers on the amount of mixing, Austin et al. [10] did not find any effect of Reynolds numbers on mixing at Re > 10000.At the same time, Romero-Gomez et al. [11] used numerical simulations to find that the main factor influencing the degree of solute mixing at high Reynolds numbers is the Reynolds number ratio rather than the Reynolds number itself.
Notably, these previous studies were focused on turbulent flow of the large Reynolds number (e.g., Re > 10000).Yu et al. [12] investigated the solute mixing at cross junctions and found that the inlet Reynolds number ratio was the most important factor affecting the mixing, followed by the pipe diameter ratio and the outlet Reynolds number ratio.Yu et al. [13] studied the solute mixing at double-Tee junctions and found that the mixing at double-Tee junctions is better than that at cross junctions.ey also found that the joint with a large connecting pipe length (L) has a more complete mixing than the one with a small connecting pipe length.Shao et al. [14] further used computational fluid modeling and solute mixing to analyze the mixing at cross junction and double-Tee junction, from which they introduced dimensionless parameter of Reynolds ratios in pipe flows to quantify the degree of solute mixing.Yu et al. [15] studied the mixing characteristics of the double-Tee junction with different diameters and found that the larger the pipe diameter ratio is, the better the mixing degree is.Song et al. [16] conducted laboratory experiments to characterize solute mixing patterns at double-Tee and wye junctions.
Comparatively, the junction mixing of the laminar and transitional flow has not been thoroughly investigated, partially because of difficulties in controlling experimental conditions.Among several studies in literature, Austin et al. [10] found that different Reynolds number has the different mixing degree at cross junction when the Reynolds number is less than 10000.Choi et al. [6] also found that the Reynolds number had a certain effect on the solute mixing at double-Tee junctions at low Reynolds numbers (Re < 10000).Although most pipe flows in WDSs are turbulent, a small part of pipe flows are in the laminar and transitional flow region in water distribution systems [17].
e importance of studying solute mixing of the laminar and transitional flow can be underestimated.e purpose of this work is to explore the mixing behavior at cross junctions and double-Tee junctions at Re < 10000, especially the laminar and transitional flow.e paper is presented in detail as follows: first, the definition of the mixing parameters and experimental approach, followed by experimental observation and results; finally, some discussions and conclusions are presented.

Mixing Parameters.
A typical cross joint and double-Tee joint used in experiments are shown in Figure 1. e pipes with geographic notations were labeled as W (west inlet, tracer water with high concentration C W ), S (south inlet, clean water with low concentration C S ), E (east outlet, concentration C E ), and N (north outlet, concentration C N ).D and L/D were also defined as pipe diameter and dimensionless connecting pipe length.Romero-Gomez et al. [18] defined the dimensionless concentration (C * ) to describe the degree of mixing: Where (3) Based on the numerical research, Shao et al. [14] found that the flow field in the mixing zone is complicated (Figure 2(a)).e flow rate is unevenly distributed at the mixing interface.e mixing interface line keeps 45 °as shown in Figure 2(b).erefore, the interaction time can be defined as where L * is a collection of line segments of the mixing interface line.l i and v i are the length of line segment and the velocity at this segment, respectively.Since the velocity at the interface is difficult to be measured and calculated, a simplified indicator mixing time is used to replace T: where D and v are the diameter and the average velocity of pipes connected to the junction.

Experimental Approach.
To better understand the mixing behavior, a set of experiments were conducted to research the influence of average Reynolds number, Reynolds number ratios of outflows (Re N/E ), and pipe diameter at cross junctions.e Reynolds numbers of the four pipes (Re S , Re W , Re N , and Re E ) in the experiments were controlled in the range from 500 to 7500 (laminar, transition or uncompleted turbulence, and turbulence flow), and the pipe diameters are DN16, DN25, DN32, and DN50.e experimental cases to investigate the influence of pipe diameters for cross junctions are presented in Table 1.e Reynolds numbers of four pipes are kept equal.To investigate the influence of Reynolds number ratios of outflows, the cross junction with diameter of DN25 is used, and the experimental cases are shown in Table 2. Experiments on double-Tee junctions were also carried out, and the cases are presented in Tables 3 and 4. Note that the experimental setup and approach were based on the work of [13].Every experiment was repeated thrice to obtain the average results.e percent mass fraction error (PMFE) is used to check the experimental accuracy [9], as shown below: e closer the absolute value of PMFE is to 0, the more reliable the data are.Each experimental case was repeated thrice independently by turning off the power of experimental system.Most of the PMFE values are −1%∼3%, which indicate that the experimental errors are well controlled.

Influence of Reynolds Number on Cross Junction.
e dimensionless concentrations of north outflow versus average Reynolds number are presented in Figure 3. e experiments were done at cross junctions of equal diameter for every 4 pipe legs.e incoming and outgoing flow rates are control to be equal, and therefore Re � Re ave for this case.e pipe diameters of the cross junctions are DN16, DN25, DN32, and DN50.e experimental results of 4 junction sizes appear to have a similar trend.e data can be divided into 3 regions, Re < 1500, 1500 < Re < 3000, and Re > 3000, respectively.When Re < 1500, the degree of mixing decreases as the Reynolds number increases.
e possible reason may be that the mixing time of the junction in this region is the primary factor to govern the mixing instead of The mixing interface  Re ave Re N/E 500 0.25, 0.5, 0.67, 1.0, 1.5, 2.0, 4.0 1500 0.25, 0.5, 0.67, 1.0, 1.5, 2.0, 4.0 2500 0.25, 0.5, 0.67, 1.0, 1.5, 2.0, 4.0 4500 0.25, 0.5, 0.67, 1.0, 1.5, 2.0, 4.0 Advances in Civil Engineering disturbance diffusion.e flow rate increases and the mixing time decreases with the increase of the Reynolds number, the mixing becomes worse, and C * N reaches the peak near Re � 1,500.Figure 4 shows that mixing time dominates junction mixing under laminar flow especially in the junction with larger pipe diameter.When 1500 < Re < 3000, the mixing increases with the increase of the Reynolds number because the turbulence intensity is starting to be a determining factor on mixing, and the effect of mixing time weakens.It is logical that with the increase of the Reynolds number, the increased turbulence intensity enhances mixing.When Re > 3000, the effects of the mixing time and turbulence intensity on mixing are reduced further, until the mixing is not affected by them.Because the mixing time and turbulence intensity are two opposed factors, increasing the flow rate will enhance turbulence intensity but reduce mixing time.With the increase of Re, the above two factors will reach a balance, and then the change of the Reynolds number will not affect the junction mixing.
Laminar flow occurs when the calculated Reynolds number is less than 2300, and turbulent flow occurs when the Reynolds number calculation exceeds 4000 for pipe flows.In the actual mixing experiment, the flows are not well controlled to be very stable and the junction configuration breaks the laminar flow condition.
e flow cannot be considered as laminar flow, and it is more like a jet stream.However, for the pipe diameters of DN16-DN50, the C * N reaches the peak at Re � 1500 and reaches stable value at Re > 3000.When Re � 500, the experimental data show a best mixing.Due to the limitation of flow measurement, the experiments with the Reynolds number less than 500 have not been done.eoretically, the mixing time can be infinitely long if the flow velocity is close to zero, and then the molecular diffusion at the junction can also produce complete mixing.
To visually analyze the mixing phenomenon, a Plexiglass cross junction was used.As is shown in Figure 5, when Re ave � 500, the mixing interface (the diagonal line of 45 °) is blurry and the mixing happened in the entire junction space (Figure 5(a)). is observed phenomenon supports the good mixing result at the low Reynolds number.When Reave � 1500, the 45 °mixing interface line is very clean (shown in Figure 5(c)), and the mixing only happened near the 45 °line.Obviously, the mixing is less sufficient than that of Re � 500.When Re � 3,000, the mixing gets better than Re � 1500 because of the effect of turbulent disturbance (Figure 5(d)), although, it is still worse than that of Re � 500.
e experimental concentration data (Figure 3) can fully support the mixing phenomenon seen in Figure 5.

Influence of Pipe Diameter on Cross Junction.
e normalized north outlet concentration C * N versus pipe diameters D under Re ave � 500∼5000 is shown in Figure 6.Due to the limitation of the flow measurement, the experiment with a Reynolds number greater than 3000 at DN50 is not carried out.From Figure 6, it appears that at the same average Reynolds number, the larger the diameter, the better mixing happens at the cross junction.is is because at the same average Reynolds number, the larger the diameter is, the larger the mixing zone (the junction geometry space) is, and the lower flow rate leads to more mixing time.erefore, mixing at cross junction of a large diameter is better than that of a small diameter.However, the effect of the diameter size is weakened with the increase of the pipe diameter.Especially, when Re > 3000, the size of the pipe diameter almost has no effect on the mixing.

Influence of Reynolds Number Ratios at Outlets (Re N/E ) for
Cross Junction.Romero-Gomez et al. [11] had studied the effect of Reynolds numbers ratios of outflows on cross junction at high Reynolds numbers.Our experimental results focused on Re < 10000 are shown in Figure 7. e experiment was done at cross junction with equal diameter Advances in Civil Engineering of DN25.e flow rates of incoming south and west flows were set to be equal.e north and east outlet flow rates were changed to follow 0.25 < Re N/E < 4.0; however, the average Reynolds numbers were kept to be 500, 1500, 2500, and 4500.With the increase of the Reynolds number ratio at outlets (Re N/E ), the mixing becomes better at the north outlet, but worse at the east outlet.is is because, the inlet flow rates are equal, which flow together towards the north outlet and largely blocked the water flow from the west inlet to the east outlet, when the north outlet Reynolds number is greater than the east outlet Reynolds number (i.e., Re N/E > 1).e average Reynolds number has almost no effect on the mixing degree when Re N/E increases above 1.5.Because when Re N/E increases over 1.5, the two outflows have a large flow rate difference, and the Reynolds number ratio of outflows (Re N/E ) is the primary factor to govern the mixing behavior instead of the average Reynolds number.When Re N/E is about 1, the average Reynolds number has an obvious impact on mixing, which is also presented in Figure 2.

Influence of Dimensionless Connecting Pipe Lengths (L/D) on Double-Tee Junctions.
e influence of the average Reynolds number and dimensionless connecting pipe length on mixing at the double-Tee junction with equal pipe size (DN25) is shown in Figure 8.It appears that the dimensionless connecting pipe length (L/D) can significantly affect the mixing ratio; mixing at double-Tee junction becomes better with the increase of L/D. is is because the long connecting pipe length provides more mixing space and time for the water flow mixing.Yet, at large L/D values, the influence of the average Reynolds number on mixing is little evident.e average Reynolds number is not the determined factor on mixing for double-Tee junctions.

Influence of Reynolds Number Ratio at Outlets (Re N/E )
on Double-Tee Junction.Yu et al. [13] found that the Reynolds number ratio of outflows (Re N/E ) will influence the mixing ratio on double-Tee junctions at high Reynolds numbers (Re > 50000).To understand the effect of Reynolds number ratio of outflows (Re N/E ) on double-Tee junction at low Reynolds numbers (Re < 10000), a series of experiments (L/D � 2.5, 5.0, 7.0, and 10.0) have been carried out, and the results are shown in Figure 9.
e experimental result of L/D � 2.5 is shown in   Advances in Civil Engineering 3.6.Comparison with Previous Studies.Austin et al. [10] only showed the mixing result on a cross junction with 16 mm inner diameter; hence, a comparison with our experimental result with 16 mm inner diameter is shown in Figure 10.When Re > 3,000, the overall trend of the two curves is close.When 1500 < Re < 3000, these two curves show the same trend that the mixing enhances with the increase of the Reynolds number.e biggest difference appears when Re < 1500.For the curve of 16 mm inner diameter, because of the experimental condition limit of the smallest Re � 500, the descent stage of the curve when Re < 500 is lost in Figure 10.However, we can infer that there must be a descent stage from Re � 500 to Re � 0. is tendency can be found for the DN50, DN32, and DN25 junction, as shown in Figure 3.Because the smallest Re in Austin et al. [10] is 700, it is not strange that the descent stage is not presented in their mixing curve.
Moreover, the geometries of the two junctions used in these two experiments are different.
is study uses an integrated manufacturing Plexiglass junction shown in Figure 11.e inner surface of the junction is smoother than that of Austin et al. [10].Figure 12 shows the detail of the experimental junction used by Austin et al. [10].e junction is made by combining a four-way connection and four pipe legs.ere is a small obstacle in the connecting location, which will change the flow regime locally.It has a significant influence in mixing, especially when the flow regime is laminar flow or transitional flow.at is to say, the junction of Austin et al. [10] has better mixing than this study for the same Reynolds number because of the enhanced disturbance by the junction of Austin et al. [10].Figure 10 supports this inference when 1000 < Re < 4000.
e diameter of the pipe has a great influence on the mixing ratio when the Re < 3000.It can be seen from Figure 3 that when the Reynolds number is between 500 and 1500, the larger the pipe diameter, the better the mixing at the same Reynolds number.
e experimental results of the junction with 16 mm inner diameter are similar to those of Austin et al. [10].e pipe diameter determines the size of the mixing zone.When the flow rate is low and the flow state is laminar, the mixing time dominates the junction mixing instead of turbulent diffusion.When the Reynolds number increases, the flow rate increases, the mixing time decreases, and the degree of mixing is deteriorated.Moreover, when the Reynolds number increases and the flow state approaches the transitional flow, the influence of the turbulent diffusion becomes larger and exceeds the effect of the mixing time.
erefore, the experimental results with a diameter bigger than 25 mm in this paper appear as a convex curve in the interval of 500 to 1500.On the other hand, when Re ⟶ 0, molecular diffusion dominates the mixing, and sufficient mixing time causes the mixing to be complete.

Conclusions
Solute mixing phenomena under laminar and transitional flow conditions were investigated at cross and double-Tee junctions.ese physical configurations and flow conditions prevail at small water demand areas, suburbs, mostly in perimeters of a water distribution system (WDS).us, study of solute mixing of the laminar and transitional flow can be important to the junction mixing model, particularly for water quality simulations.
For the cross junction, C * N versus Re ave relationship can be characterized in three hydraulic regions.Under Re < 1500 or laminar flow, mixing decreases as the Reynolds number increases.In transitional flow (500 < Re < 3000), mixing increases as the Reynolds number increases; and lastly in turbulent flow of Re > 3000, the degree of mixing is constant invariant of flow variations.Mixing at cross junction of a large diameter is better than that of a small diameter when Re < 3000.Under this condition, mixing becomes better at the north outlet, but worse at the east outlet, with the increase of the Reynolds number ratio at outlets (Re N/E ).
For the double-Tee junction, the influence of the average Reynolds number on junction mixing is less prevalent than for the cross junction under laminar and transitional flow conditions.Mixing at double-Tee junction becomes better with the increase of L/D.At higher values > 5.0, L/D becomes a determined factor on mixing, while Reynolds number ratios of outflows (Re N/E ) and the average Reynolds number become less important.e longer connecting pipe length provides more mixing space and time for the water flows mixing.
In this study, experimental research is used for understanding the junction mixing under laminar and transitional flow.Different pipe diameters and Reynolds number ratios of outflows are considered in this paper.e geometry characteristic difference of the junctions used in this study and Austin et al. [10] causes the slight difference of these results.e diameter of the pipe has a great influence on the mixing ratio when Re < 3000.e lack of experimental data of Re below 700 and only the 16 mm diameter junctions used in Austin et al. [10] caused the mixing curve difference.Advances in Civil Engineering

Figure 2 :
Figure 2: Flow velocity vectors and contours map for cross junction.Contours and arrows indicate the normalized flow velocity and direction.

Figure 1 :
Figure 1: Hydraulic and water quality configurations of the junctions: (a) cross junction; (b) double-Tee junction.

Figure 3 :
Figure 3: e dimensionless concentration of north outflow C * N

Figure 8 :Figure 9 :Figure 10 :
Figure 8: Normalized north outlet concentration C * N with dif-L/D for varied average Reynolds number at double-Tee junction.

8
Advances in Civil EngineeringFuture work should be focused on the mechanism of the junction mixing, especially under laminar and transitional flow.

Figure 11 :
Figure 11: e geometry of joint used in this study.

Figure 12 :
Figure 12: e geometry of joint used by Austin et al. [10].

Table 1 :
Summary of the experimental cases at cross junctions with Re W/S � 1 and Re N/E � 1.

Table 2 :
Summary of the experimental cases at cross junctions with different Re N/E and Re W/S � 1.

Table 3 :
Summary of the experimental cases at double-Tee junctions with Re W/S � 1 and Re N/E � 11.

Table 4 :
Summary of the experimental cases at double-Tee junctions with varied Reynolds number ratio at outlets (Re N/E ) and Re W/E � 1.