Next Article in Journal
Direct and Joint Effects of Genotype, Defoliation and Crop Density on the Yield of Three Inbred Maize Lines
Next Article in Special Issue
Scheduling Regulated Deficit Irrigation with Leaf Water Potential of Cherry Tomato in Greenhouse and its Effect on Fruit Quality
Previous Article in Journal
Determinants of Rice Farmers’ Willingness to Pay for Conservation and Sustainable Use of Swampy Wetlands in Ghana’s Northern and Ashanti Regions
Previous Article in Special Issue
Can Sustained Deficit Irrigation Save Water and Meet the Quality Characteristics of Mango?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Emitter Structure on Its Hydraulic Performance Based on the Vortex

School of Water Resources Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China
*
Author to whom correspondence should be addressed.
Agriculture 2021, 11(6), 508; https://doi.org/10.3390/agriculture11060508
Submission received: 19 April 2021 / Revised: 24 May 2021 / Accepted: 27 May 2021 / Published: 30 May 2021
(This article belongs to the Special Issue Future of Irrigation in Agriculture)

Abstract

:
The rectangular labyrinth emitter is taken as the study object in this article, as we added internal teeth to vortex-free and vortex areas in its lateral channel or lengthened the vertical channel, to change the channel structure. Using the computational fluid dynamics (CFD) method simulates the water flow field, to get the relationship between flow rate and pressure, and the vortexes distribution in channel. The aim of this study is to explore the reasons for the influence of structural change on hydraulic performance of the emitter through the analysis of vortex intensity and its distribution from the perspective of the vortex. The results show that the relative error of simulated results and experimental data was 1.02–2.11%. Adding internal teeth to vortex-free areas in lateral channel can improve hydraulic performance of the emitter; adding them to vortex areas can reduce it. The increase in vortex number and intensity in flow field is the internal reason for the improvement of the emitter’s hydraulic performance. The channel structure changes promote the formation of a larger velocity gradient, and the increase in the velocity gradient in flow field exacerbates vortex formation. Changing channel structure to improve the emitter’s hydraulic performance can promote an increase in the number and intensity of vortexes in the channel.

1. Introduction

Emitters are key components of drip irrigation systems; their hydraulic performance directly determines the working performance of drip irrigation systems [1,2,3]. An emitter with superior hydraulic performance has low flow rate sensitivity to inlet pressure and can better ensure the irrigation uniformity of drip irrigation systems. The hydraulic performance of the emitter is usually expressed as:
q = k · h x
where k is the flow coefficient; x is the flow index; q is the flow rate of the emitter, L/h; and h is the inlet pressure head of the emitter, m [4].
The hydraulic performance of the emitter depends on its channel structure, size, and material [5]. In recent years, studies have shown that emitters of the labyrinth channel have superior hydraulic performance when compared with other types of emitters. For this reason, in recent years, many authors have conducted considerable amounts of research on the relationship between channel structural parameters and the flow index and flow coefficient, especially the flow index, striving to reduce it by changing the channel structure of the emitter in order to improve its hydraulic performance. Yu et al. [6] studied the relationship among five key structural parameters and the flow index and flow coefficient through the orthogonal experimental design method, and gave the regression equations regressed by the linear multivariable regression method regarding the relationship between the flow coefficient and flow index and the structural parameters. Xie et al. [7] analyzed the influence of the tooth angle and pitch on the flow coefficient and flow index, indicating that the tooth angle, tooth pitch, and flow coefficient are positively correlated, and the flow coefficient increased with the increase in the tooth angle and pitch. Tang et al. [8] noted that the flow index of the drip irrigation emitter decreased with the increase in the emitter draft angle, the angle of the tooth tip, and the distance between the opposite tooth tips, which also decreased with the decrease in the fillet radius of the tooth tip, additionally, the flow index is unrelated to the changes in the height of the channel. Tian et al. [9] selected nine major structural parameters, using multivariable regression methods, to establish the regression equations regarding the relationship between the flow index and flow coefficient with the nine structural parameters, in the process of their research.
The sensitivity of flow rate to the inlet pressure not only depends on the flow index but is also affected by the flow coefficient. Therefore, the study of the hydraulic performance of the emitter should analyze the slope of its flow rate-pressure curve. In addition, most of the literature is devoted to studying the relationship between channel structure parameters and the flow index and the significance analysis of the influence of channel structure parameters on the flow index [10,11]. Many scholars have also conducted related research regarding the energy dissipation mechanism in the channel of the labyrinth emitter. Yan et al. [12] revealed the correlation relationship between the dentate structure of the flow path of emitters and energy dissipation. The results showed that the rears of dentations and dentate tips were the primary regions for energy dissipation. Xing et al. [13] analyzed the energy dissipation mechanism of perforated drip irrigation emitters. The results showed that the main energy consumption of PDIE was due to the mixing of two streams formed by the scalariform perforation plate. However, they only revealed the energy dissipation in the flow channel. The research on the internal reason for the flow channel structure causing the change in the hydraulic performance of the emitter, that is, the mechanism study of the influence of channel structure on the emitter’s hydraulic performance, is inadequate.
In this article, we use a larger channel cross section (1.5 mm × 1.5 mm) to ensure the anti-clogging performance of the emitter, and, at the same time, we add internal teeth (0.5 mm high) to the channel to enhance the hydraulic performance. The section size after the internal teeth is bigger than the clogging sensitive size (0.7 mm) [14,15,16,17]. We change the channel structure of the emitter by adding internal teeth to the channel of the commonly used rectangular labyrinth emitter. The aims of this study are as follows: to analyze the internal reasons for the change in the hydraulic performance of the emitter caused by the structural parameter change from the perspective of the vortex; to explore the correlation among the structure of the emitter, the vortex, and the hydraulic performance; and, to provide a certain theoretical and reference basis for the research and development of irrigation emitters.

2. Materials and Methods

2.1. Emitter Structure and Parameters

In this paper, adding different numbers of internal teeth to vortex-free and vortex areas in the lateral channel of different units of the rectangular labyrinth emitter (A emitter) forms the internal-tooth labyrinth emitters (B, C emitters). At the same time, in order to further study the relevance of the emitter’s structure, vortex, and hydraulic performance, the vertical channel of each unit of the A emitter was lengthened to 5 times (D emitter) for comparison. This made the head loss increase caused by the lengthening of the channel equal the head loss increase caused by adding 15 teeth to vortex-free areas in the lateral channel of different units (B1 emitter), at 10 m inlet pressure head. Each model of emitter consisted of 45 units, and the cross section of the flow channel was 1.5 mm × 1.5 mm (width × depth). The dimensions of the internal teeth were as follows: 0.5 mm high and 1 mm wide; the thickness of the internal teeth was consistent with the depth of the channel. The teeth were placed from the beginning of channel units. The parameters of the emitters are shown in Table 1 and Figure 1, and the physical model of A emitter is shown in Figure 2.

2.2. Control Equations and Calculation Model

The water movement in the emitter can be regarded as the movement of incompressible fluid. The basic control equations are composed of the continuity equation and the Navier–Stokes equation [18]. Using FLUENT software for simulation calculation, studies have shown that the standard k-ε model of the labyrinth emitter calculation is consistent with the actual situation [19,20,21]. Thus, the calculation model selects standard k-ε model. The velocity and pressure coupling are solved by the SIMPLE algorithm, and the convergence accuracy is 0.0001.
Continuity equation:
u x + v y + w z = 0
Navier–Stokes equation:
( ( ρ u ) ) t + · ( ρ u U ) = p x + μ 2 u + ρ f x
( ( ρ v ) ) t + · ( ρ v U ) = p y + μ 2 v + ρ f y
( ( ρ w ) ) t + · ( ρ w U ) = p z + μ 2   w + ρ f z
where U is the velocity; u, v, w are the velocity components in the x, y, z coordinate axes; ρ is the density of water, kg/m3; μ is viscosity, Pa·s; p is pressure, Pa; and f x , f y , f z are the components of mass force, regardless of mass force.
k equation:
( ρ k ) t + ( ρ k u i   ) ( x i ) = ( x j )   [ ( μ + μ t σ k ) k ( x j ) ] + G k ρ ε
ε equation:
( ρ ε ) t + ( ρ ε u i   ) ( x i   ) = ( x j   )   [ ( μ + μ t σ ε ) ε ( x j ) ] + ( C 1 ε   ε ) k G k C 2 ε ε 2 k
Among them,
G k = μ t ( u i x j + u j x i ) u i x j
μ t = ρ C μ k 2 ε
where k is turbulent energy; ε is the turbulent dissipation rate; μ t is turbulent viscosity; t is time; μ is viscosity; x i , x j is the coordinate tensor; and G k is the production term of turbulent energy k caused by the average velocity gradient; The empirical constants are C 1 ε = 1.44, C 2 ε = 1.92, C μ   =   0.09 , σ k = 1.0, and σ ε = 1.3 [22].

2.3. Boundary Conditions and Mesh

The boundary layer in the channel of the emitter was set to encrypt the grid. The first boundary layer was 0.01 mm, and the other layers were increased by a ratio of 1.5; there were 6 boundary layers in total with 0.21 mm. In order to ensure high computational efficiency and to reduce the influence of the mesh on the calculation results, the grid sizes of 0.15, 0.12, 0.1, 0.09, 0.08, 0.07, 0.06, 0.05, and 0.04 mm, were selected to calculate the flow rate of the emitter under the inlet pressure head of 5 m, 10 m, and 15 m. Considering the calculation accuracy and efficiency, when the average difference of the calculation results was less than 0.5%, it was assumed that the number of grids does not influence the results. When the grid size was 0.08 mm, the average difference was 0.38%; therefore, the grid size was set as 0.08 mm, as shown in Figure 3.
The initial condition was 0.6–1.1 m/s inlet velocity, and the outlet condition was outflow. The standard wall function method was used to process the wall surface of the channel, and the wall surface roughness was 0.01 mm [23].

3. Results and Discussion

3.1. Analysis of the Accuracy of the Simulated Results

The verification of this numerical simulation calculation results was conducted using the experimental data in the reference document [24]. The relative error of the simulated results and the experimental data of the triangular labyrinth emitter with a flow channel width of 0.7 mm was 1.02–2.11%, therefore, the reliability of the numerical simulated calculation is described. The calculated results are shown in Table 2 and Figure 4.

3.2. Flow Rate-Pressure Curve of the Emitter

The simulation results of A, B, C, and D emitters were analyzed within 5–15 m pressure head loss in this article, obtaining the flow rate-pressure relationship; the flow coefficient, k; and the flow index, x. The data are shown in Table 3. The calculating equation of the slope of the flow rate-pressure curve of each emitter at 6, 8, 10, 12, 14 m pressure head, is q = k · x · h ( x 1 ) ; the calculation results are shown in Table 4.
From Table 3 and Table 4, it can be observed that the flow rate-pressure relationship and the slope of the C emitter curves are approximately the same. The C1 emitter has the biggest flow coefficient (1.995), and the C3 emitter has the smallest flow coefficient (1.935), with a difference of 0.06. The C3 emitter has the biggest flow index (0.568), while the C1 emitter has the smallest flow index (0.556), with a difference of 0.012. The slope of the C3 emitter’s flow rate-pressure curve is the biggest, while that of the C1 emitter is the smallest; the maximum difference is 0.007 (at 14 m pressure head). This shows that after adding internal teeth to the vortex areas in the lateral channel of different units, the number of internal teeth has a minor effect on the hydraulic performance of the emitter. In addition, the smallest flow coefficient of the C emitters (1.935, C3 emitter) is 4.2% smaller than that of the A emitter, which is 2.019. The biggest flow index of C emitters (0.568, C3 emitter) is 3.5% bigger than that of the A emitter, which is 0.549. At the same time, the slope of the C emitter’s flow rate-pressure curve is bigger than that of the A emitter. This indicates that after adding internal teeth to vortex areas in the lateral channel of different units, the emitter’s hydraulic performance does not improve, but decreases.
In summary, the A, B, D, and C1 emitters (which has the best hydraulic performance among C emitters) are selected for comparative analysis. The flow rate-pressure curves of the A, B, C1, and D emitters are shown in Figure 5. The A in the legend represents A emitter, and so on.
From Table 4 and Figure 5, it can be observed that the slopes of the B emitters flow rate-pressure curves are smaller than that of the A emitter, and the maximum curve slopes of the B1, B2, and B3 emitters decrease 16.0%, 24.9%, and 31.5% (at 14 m pressure head), respectively. This shows that the sensitivity of the B emitters flow rate to the inlet pressure is lower than that of the A emitter; the hydraulic performance is improved, and with the increase in the number of internal teeth, there is also an increase in improvement of the hydraulic performance of the emitter. Due to lengthening the vertical channel, the energy dissipation in the channel is increased, and the slope of the D emitters flow rate-pressure curve is smaller than that of the A emitter but bigger than that of the B emitters. From this, we can see that in the case of increasing the same energy loss, the method of lengthening the vertical channel to increase the energy loss has less of an effect on the improvement to emitter hydraulic performance than adding internal teeth to vortex-free areas in the lateral channel of different units does.

3.3. Analysis of Flow Field in Emitter Channels

Figure 6 and Figure 7 show the streamlined diagrams in the channel of the A and D emitter at the velocity of 1.1 m/s. It can be seen that the vortex in the A emitter mainly exists in the lateral channel in Figure 6. From Figure 8, after adding internal teeth to the vortex-free areas in the A emitter’s lateral channel of different units, we can see that although the B emitter’s main flow squeezes the original vortex around the internal teeth; a new vortex appears behind the tooth, and a bigger vortex group appears in the vertical channel behind the tooth. From Figure 9, After adding internal teeth to vortex areas in the A emitter’s lateral channel of different units, not only do the internal teeth destroy the original vortex, but no vortex is formed behind the tooth and in the vertical channel. At the same time, it is known from the previous section that the B emitter’s hydraulic performance is better than that of the A emitter, while the hydraulic performance of the C1 emitter is similar to that of the A emitter, and even less than that of the A emitter. This indicates that after adding internal teeth at different positions in the flow channel of the emitter, the number and distribution of vortexes in the flow channel changes, which, in turn, influences the hydraulic performance of the emitter.
Table 5 shows the sum of the vortex characteristics in a flow channel unit of the A, B, C1, and D emitters. The vortex intensity is defined as ω × A , where A is the cross-sectional area of the vortex; ω is the rotational angular velocity, ω = v r ; r is the equivalent radius of the vortex area; v is the average tangential velocity at the boundary of the vortex area [25]. From Table 5, it can be observed that the number, the sum of area, and the intensity of the vortex in the B emitters are the largest in a flow channel unit. Compared with the A emitter, the sum of intensity of the vortex in the B1, B2, and B3 emitters increases by 1.8 times, 3 times, 3.2 times, respectively. In Figure 7, the number of vortexes in the D emitter is the same as that in the A emitter, but the sum of intensity of the vortex is 7.5% bigger than that of the A emitter from Table 5, therefore, the slope of the D emitters flow rate-pressure curve is slightly smaller than that of the A emitter, and the hydraulic performance is slightly improved. Although the number of vortexes in the C1 emitter is more than that of the A emitter, the sum of intensity of the vortex is 39.6% smaller than that of the A emitter. The number of vortexes in the C1 emitter is the same as that of the B1 emitter, but the sum of intensity of the vortex is only 21.3% of that of the B1 emitter. This indicates that the existence of the vortex in the flow channel and the increase in its intensity are the main reasons for the improvement to the emitter’s hydraulic performance.

3.4. Flow Velocity Distribution in the Emitter

To further analyze the reasons for vortex formation, we extracted the intersection-lines velocity distribution, where P1 and P2’s cross-sections intersect with the center plane of the depth direction of the channel, in the A, B, C1, and D emitters (as shown in Figure 10). The velocity distribution is shown in Figure 11 and Figure 12. The A in the legend represents the A emitter, and so on.
As shown in Figure 11, due to the presence of internal teeth, in the B emitters, the shape of P1’s velocity distribution curve becomes “thin,” and the curve heads’ velocity gradient is bigger. This causes the turning fluid to shrink further and exacerbates the generation of a new vortex in the corner of the channel behind the internal tooth, as shown in Figure 8. The P1’s velocity distribution curve in the C1 and D emitters is similar to that in the A emitter; the curve heads’ velocity gradient is small, so that no new vortex is generated. In addition, as shown in Figure 9, after adding internal teeth to the original vortex area of the C1 emitter, the original vortex is destroyed, which not only fails to make the fluid shrink further; but also causes the velocity distribution at P1’s position to become more uniform and reduces the turbulence degree of the flow field and energy dissipation. This also further explains the reason why the hydraulic performance of the C1 emitter slightly decreases.
As shown in Figure 12, the curves velocity gradient remains large in the B emitter’s vertical channel P2 position behind the internal tooth; the maximum velocity is 30% bigger than that of the A emitter, and new vortexes are generated. This indicates that the internal teeth accelerate the main flow velocity in the vertical channel behind the internal tooth and turn the original low-speed area into vortex areas, which greatly increase the degree of turbulence in the vertical flow channel and intensify energy dissipation. This is also one of the reasons why the hydraulic performance is improved. The vertical channel P2’s velocity distribution of the C1 and the A emitter is basically the same. Due to the lengthening of the vertical channel, the D emitter’s velocity distribution is close to that of the turbulent flow in the direct channel, which tends to be uniform. The maximum velocity is 13% smaller than that of the A emitter.
In this article, we changed the channel structure of the emitter by adding internal teeth to the lateral flow channel and by lengthening the vertical flow channel, and the hydraulic performance of the emitter changed. Based on the above analysis, the change in the flow channel structure of the emitter changes the velocity distribution in the flow field. The addition of internal teeth to the vortex-free area increases the maximum flow velocity and its gradient in the toothed flow channel unit. The increase in the velocity gradient exacerbates the generation of vortexes in the flow channel, and the increase in the number and intensity of the vortexes in the flow channel is the internal reason for the improvement to the hydraulic performance of the emitter. In addition, adding internal teeth at a reasonable position can turn the original low-speed area into a vortex area. The full development of the vortex improves the self-cleaning ability of the water flow to the flow channel wall and improves the anti-clogging ability of the emitter to a certain extent [26].

4. Conclusions

Analyzing the simulation results of the A, B, C, and D labyrinth emitters, the following conclusions can be drawn:
  • Adding internal teeth to vortex-free areas in the lateral channel of different units of rectangular labyrinth emitters can improve hydraulic performance, and doing so to vortex areas not only fails to improve the hydraulic performance but can reduce it.
  • Adding internal teeth is not the reason for the improvement in the emitter’s hydraulic performance. The increase in the number, distribution, and intensity of vortex in the flow field is the internal reason for the improvement in the emitter’s hydraulic performance.
  • In the case of increasing the same energy loss, the method of lengthening the vertical flow channel has less of an effect on the improvement in the emitter’s hydraulic performance than adding internal teeth to increase the vortex characteristics.
  • The change in the flow channel structure promotes the formation of a larger velocity gradient in the flow field, and the increase in the velocity gradient exacerbates the generation of the vortex in the flow field.
  • Changing the flow channel structure to improve the emitter’s hydraulic performance can promote an increase in the number and intensity of vortexes in the emitter and avoid damage to the original vortex structure in the flow field.

Author Contributions

C.W. edited the original draft; Z.L. reviewed and edited the draft; J.M. provided the funding. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (52079085).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the finding of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank the School of Water Resources Science and Engineering for providing the site for experiment and collecting data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, J.; Zhao, W.; Lu, B. Rapid prediction of hydraulic performance for emitters with labyrinth channels. J. Irrig. Drain. Eng. 2013, 139, 414–418. [Google Scholar] [CrossRef]
  2. Yunkai, L.; Ji, F. Progress in measurement of hydrodynamic characteristics in drip irrigation emitters. J. Drain. Irrig. Mach. Eng. 2014, 32, 86–92, (In Chinese with English Abstract). [Google Scholar]
  3. Xiaohong, W. Structural Design and Experimental Study on Hydraulic Performance of Emitters with Labyrinth Channel. Master’s Thesis, Xi’an University of Technology, Xi’an, China, 2006. [Google Scholar]
  4. Keller, J.; Karmeli, D. Trickle irrigation design parameters. Trans. ASAE 1974, 17, 0678–0684. [Google Scholar] [CrossRef]
  5. Mattar, M.A.; Alamoud, A.I.; Al-Othman, A.A.; Elansary, H.O.; Farah, A.-H.H. Hydraulic performance of labyrinth-channel emitters: Experimental study, ann, and gep modeling. Irrig. Sci. 2020, 38, 1–16. [Google Scholar] [CrossRef]
  6. Liming, Y. Relations between the hydraulic performance and the structural parameters of emitter. J. Chang. Univ. Sci. Technol. 2011, 08, 30–35, (In Chinese with English Abstract). [Google Scholar]
  7. Qiaoli, X.; Wenquan, N.; Lianzhong, L. Effect of tooth angle and pitch of labyrinth channel on performance of emitter. J. Drain. Irrig. Mach. Eng. 2013, 31, 449–455, (In Chinese with English Abstract). [Google Scholar]
  8. Ye, T.; Minjie, W.; Ming, L.; Guojun, D. The effect of drip irrigation emitter channel structure parameters of triangular type on its hydraulic performances. Model. Simul. 2018, 7, 55–62, (In Chinese with English Abstract). [Google Scholar]
  9. Jiyang, T.; Dan, B.; Changjiang, R.; Wang, X. Analysis on hydraulic performance of bidirectional flow channel of drip irrigation emitter. Trans. CSAE 2013, 89–94, (In Chinese with English Abstract). [Google Scholar]
  10. Chunjing, L.; Dunbing, T.; Hua, H.; Xingqiang, C. Hydraulic performance robustness analysis for drip irrigation triangle labyrinth channel of emitter. Trans. CSAM 2013, 44, 67–72, (In Chinese with English Abstract). [Google Scholar]
  11. Souza, W.d.J.; Rodrigues Sinobas, L.; Sanchez, R.; Botrel, T.A.; Coelho, R.D. Prototype emitter for use in subsurface drip irrigation: Manufacturing, hydraulic evaluation and experimental analyses. Biosyst. Eng. 2014, 128, 41–51. [Google Scholar] [CrossRef] [Green Version]
  12. Yan, D.; Yang, P.; Ren, S.; Li, Y.; Xu, T. Numerical study on flow property in dentate path of drip emitters. N. Z. J. Agric. Res. 2007, 50, 705–712. [Google Scholar]
  13. Xing, S.; Wang, Z.; Zhang, J.; Liu, N.; Zhou, B. Simulation and verification of hydraulic performance and energy dissipation mechanism of perforated drip irrigation emitters. Water 2021, 13, 171. [Google Scholar] [CrossRef]
  14. Camp, C.R. Subsurface drip irrigation: A review. Trans. ASAE 1998, 41, 1353–1367. [Google Scholar] [CrossRef]
  15. Yunkai, L.; Peiling, Y.; Shumei, R. General review on several fundamental points of design theory about flow path in drip irrigation emitters. Trans. CSAM 2006, 37, 145–149, (In Chinese with English Abstract). [Google Scholar]
  16. Adin, A.; Sacks, M. Dripper-clogging factors in wastewater irrigation. J. Irrig. Drain. Eng. 1991, 117, 813–826. [Google Scholar] [CrossRef]
  17. Nakayama, F.S.; Bucks, D.A. Water quality in drip/trickle irrigation: A review. Irrig. Sci. 1991, 12, 187–192. [Google Scholar] [CrossRef]
  18. Chen, X.; Wu, P.; Fan, X.; Niu, W. Numerical simulation of structural parameters using CFD and anti-clogging design analyse on drip emitter. J. Irrig. Drain. 2008, 27, 35–38, (In Chinese with English Abstract). [Google Scholar]
  19. Chang, Y.; Niu, W.; Wang, W. Numerical simulation and flow analysis of labyrinth path of drip irrigation emitters. J. Northwest A&F Univ. 2009, 37, 203–208, (In Chinese with English Abstract). [Google Scholar]
  20. Xu, T.; Zhang, L. Influence and analysis of structure design and optimization on the performance of a pit drip irrigation emitter*. Irrig. Drain. 2020, 69, 633–645. [Google Scholar] [CrossRef]
  21. Miaoye, K. Influence of Double Internal Tooth Parameters on Hydraulic Performance of Labyrinth Irrigation Device. Master’s Thesis, Taiyuan University of Technology, Taiyuan, China, 2018. [Google Scholar]
  22. Fujun, W. Computational Fluid Dynamics Analysis-CFD Software Principle and Application; Tsinghua University Press: Beijing, China, 2004; Volume 1, pp. 7–12. [Google Scholar]
  23. Yongxin, L.; Guangyong, L.; Xiangyu, Q.; Wang, J. Modeling of hydraulic characteristics through labyrinth emitter in drip irrigation using computational fluid dynamics. Trans. CSAE 2005, 21, 12–16, (In Chinese with English Abstract). [Google Scholar]
  24. Wei, Q.; Shi, Y.; Dong, W.; Gang, L.; Huang, S. Study on hydraulic performance of drip emitters by computational fluid dynamics. Agric. Water Manag. 2006, 84, 130–136. [Google Scholar] [CrossRef]
  25. Zhichang, Z. Hydraulics, 2nd ed.; China Water & Power Press: Beijing, China, 2016; Volume 7, pp. 272–277. [Google Scholar]
  26. Feng, J.; Li, Y.; Wang, W.; Xue, S. Effect of optimization forms of flow path on emitter hydraulic and anti-clogging performance in drip irrigation system. Irrig. Sci. 2017, 36, 37–47. [Google Scholar] [CrossRef]
Figure 1. Dimensions of emitters (mm). (a) Dimensions of A emitter; (b) Dimensions of D emitter.
Figure 1. Dimensions of emitters (mm). (a) Dimensions of A emitter; (b) Dimensions of D emitter.
Agriculture 11 00508 g001
Figure 2. Physical model of A emitter.
Figure 2. Physical model of A emitter.
Agriculture 11 00508 g002
Figure 3. Grid division of the emitter.
Figure 3. Grid division of the emitter.
Agriculture 11 00508 g003
Figure 4. Curves of relationships of pressure versus simulated and experimented flow rate.
Figure 4. Curves of relationships of pressure versus simulated and experimented flow rate.
Agriculture 11 00508 g004
Figure 5. Flow rate-pressure curve of emitters.
Figure 5. Flow rate-pressure curve of emitters.
Agriculture 11 00508 g005
Figure 6. Streamlined diagram of the A emitter.
Figure 6. Streamlined diagram of the A emitter.
Agriculture 11 00508 g006
Figure 7. Streamlined diagram of the D emitter.
Figure 7. Streamlined diagram of the D emitter.
Agriculture 11 00508 g007
Figure 8. Streamlined diagram of the B emitters. (a) Streamlined diagram of the B1 emitter; (b) Scheme 2 emitter; (c) streamlined diagram of the B3 emitter.
Figure 8. Streamlined diagram of the B emitters. (a) Streamlined diagram of the B1 emitter; (b) Scheme 2 emitter; (c) streamlined diagram of the B3 emitter.
Agriculture 11 00508 g008
Figure 9. Streamlined diagram of the C1 emitter.
Figure 9. Streamlined diagram of the C1 emitter.
Agriculture 11 00508 g009
Figure 10. Depth directions center plane schematic diagram.
Figure 10. Depth directions center plane schematic diagram.
Agriculture 11 00508 g010
Figure 11. Velocity distribution in P1 position of A, B, C1 and D emitter.
Figure 11. Velocity distribution in P1 position of A, B, C1 and D emitter.
Agriculture 11 00508 g011
Figure 12. Velocity distribution in P2 position of A, B, C1 and D emitter.
Figure 12. Velocity distribution in P2 position of A, B, C1 and D emitter.
Agriculture 11 00508 g012
Table 1. The structure and parameters of the emitters.
Table 1. The structure and parameters of the emitters.
TypeInternal Teeth
Position
Schematic Tooth
Numbers
A/ Agriculture 11 00508 i0010
B1vortex-free area Agriculture 11 00508 i00215
B2vortex-free area Agriculture 11 00508 i00330
B3vortex-free area Agriculture 11 00508 i00445
C1vortex area Agriculture 11 00508 i00515
C2vortex area Agriculture 11 00508 i00630
C3vortex area Agriculture 11 00508 i00745
D/ Agriculture 11 00508 i0080
Table 2. The relative error of the flow rate.
Table 2. The relative error of the flow rate.
MethodPressure Head (m)
57.51012.51517.52022.525
Simulation
results (L/h)
1.992.422.763.063.323.583.804.084.21
Experimental
results (L/h)
1.972.372.713.003.273.513.733.944.13
Relative error1.02%2.11%1.85%2.00%1.53%1.99%1.88%2.03%1.94%
Table 3. Flow rate-pressure relationship of the emitter.
Table 3. Flow rate-pressure relationship of the emitter.
Typekx
A2.0190.549
B11.9030.523
B21.7700.515
B31.6700.508
C11.9950.556
C21.9580.563
C31.9350.568
D1.5760.603
Table 4. The calculation results of the flow rate-pressure curve slopes.
Table 4. The calculation results of the flow rate-pressure curve slopes.
Pressure Head/mC3C2C1ADB1B2B3
60.5070.5040.5010.4940.4660.4240.3820.351
80.4470.4450.4410.4340.4160.3700.3320.305
100.4060.4040.4000.3920.3810.3320.2980.273
120.3750.3730.3690.3610.3540.3050.2730.250
140.3510.3480.3440.3370.3330.2830.2530.231
Table 5. Vortex characteristics (in a channel unit).
Table 5. Vortex characteristics (in a channel unit).
TypeVortex NumbersTotal Vortex Area/mm2Total Vortex Intensity
A21.4801.689
B154.9704.790
B289.4986.787
B389.5367.125
C150.7201.020
D21.7901.815
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, C.; Li, Z.; Ma, J. Influence of Emitter Structure on Its Hydraulic Performance Based on the Vortex. Agriculture 2021, 11, 508. https://doi.org/10.3390/agriculture11060508

AMA Style

Wang C, Li Z, Ma J. Influence of Emitter Structure on Its Hydraulic Performance Based on the Vortex. Agriculture. 2021; 11(6):508. https://doi.org/10.3390/agriculture11060508

Chicago/Turabian Style

Wang, Cuncai, Zhiqin Li, and Juanjuan Ma. 2021. "Influence of Emitter Structure on Its Hydraulic Performance Based on the Vortex" Agriculture 11, no. 6: 508. https://doi.org/10.3390/agriculture11060508

APA Style

Wang, C., Li, Z., & Ma, J. (2021). Influence of Emitter Structure on Its Hydraulic Performance Based on the Vortex. Agriculture, 11(6), 508. https://doi.org/10.3390/agriculture11060508

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop