ABSTRACT:
DIMSUB is a computer program to complement a decision support tool (DST) to effectively study different hydraulic design alternatives in microirrigation systems. We developed environments in Visual Basic for applications for Microsoft Excel^{®} that allow specific stepbystep functions to be created for the design of irrigation subunits. Different alternatives can be considered, such as types of emitter, lateral and submain pipe sizes, different feeding points, irregular subunit shapes and topography slopes. Furthermore, specific uniformity criteria need to be considered to achieve efficient water applications and proper design systems. Lengths of run lateral and submain pipes, position of the hydrant connection, pressure head and head loss in pipes or pressurecompensating emitters can be assigned to evaluate the results and choose the best design alternative. This userfriendly tool to study hydraulic variables is expected to be a valuable aid for the decisionmaking process in designing irrigation systems. Some examples of practical cases under specific crop conditions to design drip irrigation subunits are given using DIMSUB.
Keywords:
irrigation; emission uniformity; dripper; manifold
Introduction
Modernization of irrigation systems in the Mediterranean Basin is necessary in rural areas, due to water scarcity. Changing to pressure irrigation systems, in which water can be effectively controlled and supplied, requires systems to be appropriately designed and managed.
In microirrigation, a subunit is the set of lateral pipes connected to a manifold and controlled by a manual or automatic pressureregulating valve (Figure 1). The correct design of this system is essential to distribute water and fertilizers uniformly in the field (Wu and Barragan, 2000Wu, I.P.; Barragan, J. 2000. Design criteria for microirrigation systems. Transactions of the ASAE 43: 11451154., Holzapfel et al., 2009, Carrión et al., 2014Carrión, F.; Montero, J.; Tarjuelo, J.M.; Moreno, M. 2014. Design of sprinkler irrigation subunit of minimum cost with proper operation: application at corn crop in Spain. Water Resources Management 28: 50735089.). Attempts have been made by several researchers to design the hydraulics of subunits. Myers and Bucks (1972)Myers, L.E.; Bucks D.A. 1972. Uniform irrigation with low pressure trickle system. Journal of Irrigation and Drainage 98: 341346., Howell and Hiler (1974), Wu and Gitlin (1975Wu, I.P.; Gitlin, H.M. 1975. Energy gradient line for drip irrigation laterals. Journal of Irrigation and Drainage 101: 323–326, ^{1983}Wu, I.P.; Gitlin, H.M. 1983. Drip irrigation application efficiency and schedules. Transactions of the ASAE 26: 9299.) studied the behavior of flow conditions along lateral pipes and the manifold. Wu (1997)Wu, I.P. 1997. An assessment of hydraulic design of microirrigation systems. Agricultural Water Management 32: 275284. studied the analysis of single diameter laterals based on the energy line analytical approach. Bralts and Segerlind (1985)Bralts, V.F.; Segerlind, L.J. 1985. Finite elements analysis of drip irrigation submain unit. Transactions of ASAE 28: 809–814. and Bralts et al. (1993)Bralts, V.F.; Driscoll, M.A.; Shayya, W.H.; Cao, L. 1993. An expert system for the hydraulic analysis of microirrigation systems. Computers and Electronics in Agriculture 9: 275287. used the finite element numerical approach to analyze microirrigation systems. Pitts et al. (1986)Pitts, D.J.; Zajueta, F.S.; Smajstrla, A.G. 1986. Micro Irrigation System Submain Design Evaluation. University of Florida, Gainesville, FL, USA. used the stepbystep numerical method for the trickle irrigation design. Hoffman et al. (2007)Hoffman, G.J.; Evans, R.G.; Marvin, E.J.; Derrel, L.M.; Elliott, R.L. 2007. Design and Operation of Farm Irrigation Systems. 2ed. ASABE, St. Joseph, MI, USA. widely reviewed the design of microirrigation systems. Recently, researchers have worked on this issue to improve the design using computeraided design techniques. Ravindra et al. (2008)Ravindra, V.K.; Singh, R.P.; Mahar, P.S. 2008. Optimal design of pressurized irrigation subunit. Journal of Irrigation and Drainage Engineering 134:137146. presented a linear programming model for optimal design of pressurized irrigation subunits. Carrión et al. (2013)Carrión, F.; Tarjuelo, J.; Hernandez, D.; Moreno, M. 2013. Design of microirrigation subunit of minimum cost with proper operation. Irrigation Science 31: 11991211. developed Presud to identify the optimum microirrigation subunit design using the total annual water application cost per unit area (investment, maintenance, energy and water consumption). Pedras et al. (2009)Pedras, C.M.G.; Pereira, L.S.; Gonçalves, J.M. 2009. MIRRIG: a decision support system for design and evaluation of microirrigation systems. Agricultural Water Management 96: 691701. proposed a complex multicriteria decision support tool (DST) to generate different design scenarios.
The main objective of the subunit design is therefore to reach maximum application efficiency by satisfying the plant needs, for which highly uniform water distribution is required.
This paper describes a decision support tool, DIMSUB, aimed to design irrigation subunits. Case studies are evaluated with this computer application entering all the necessary variables required for the design of irrigation subunits in different terrains. Main hydraulic parameters to achieve uniform irrigation are obtained and designers can decide the best solution in terms of water distribution, pressure head requirements and material cost.
Materials and Methods
The design of irrigation subunits involves the correct selection and size of all components in the system, especially the type of emitter, lateral and manifold pipes, ensuring their technical characteristics.
Field uniformity is affected by a combination of a number of design variables that influence pressure and flow rate at the emitter as well as the pressure variation permitted in the system (Wu, 1997Wu, I.P. 1997. An assessment of hydraulic design of microirrigation systems. Agricultural Water Management 32: 275284.; Pereira et al., 2002Pereira, L.S.; Oweis, T.; Zairi, A. 2002. Irrigation management under water scarcity. Agricultural Water Management 57: 175206.). Even quality of emitters, expressed with the manufacturer's coefficient of variation (CV), affects the uniformity of field emission (Carrión et al., 2013Carrión, F.; Tarjuelo, J.; Hernandez, D.; Moreno, M. 2013. Design of microirrigation subunit of minimum cost with proper operation. Irrigation Science 31: 11991211.; Zhang et al., 2013Zhang, L.; Wu, P.; Zhu, D. 2013. Hydraulic design procedure for drip irrigation submain unit based on relative flow difference. Irrigation Science 31: 10651073.).
The relationship between the emitter flow rate and pressure head is established by emitter hydraulic characteristic:
where: q = emitter flow rate (L h^{−1}); H = pressure head (m); k_{e} = emitter discharge coefficient; x = exponent of discharge, close to 0 for pressurecompensating (PC) emitters and 0.5 for nonpressurecompensating (NPC) emitters.
Variation in the emitter flow rate of an irrigation subunit is defined by the expression:
where: q_{max} = maximum emitter flow rate (L h^{−1}); q_{min} = minimum emitter flow rate (L h^{−1});
The hydraulic principles involved in designing subunits are based on the fundamental conservation of energy and continuity equation. Pressure evolution in manifolds and laterals is mainly due to the slope and head losses in pipelines.
The total energy in a pipe, H, between sections 1 and 2 per unit weight of water (m), can be expressed by the steady flow energy equation:
where: z = elevation as potential energy (m);
Friction head losses can be estimated by the DarcyWeisbach equation, in which the modified Blasius equation (1913)Blasius, H. 1913. The law of similarity for friction in fluids = Das aehnlichkeitsgesetz bei reibungsvorgangen in flussigkeiten. Forshhungsarbeiten auf dem Gebiete des Ingenieurwesen 131(1) (in German).) can be assumed to estimate the friction factor f for practically turbulent flow in smooth pipes:
where: f = friction factor and Re = Reynolds number within the range 3,000<Re<10^{5} (Von Bernuth, 1990Von Bernuth, R.D. 1990. Simple and accurate friction loss equation for plastic pipe. Journal of Irrigation and Drainage Engineering 116: 294–298.; Bagarello et al., 1995Bagarello, V.; Ferro, V.; Provenzano, G.; Pumo, D. 1995. Experimental study on flow resistance law for smalldiameter plastic pipes. Journal of Irrigation and Drainage Engineering 121: 313–316.; Provenzano et al., 2005Provenzano, G.; Pumo, D.; Di Dio, P. 2005. Simplified procedure to evaluate head losses in drip irrigation laterals. Journal of Irrigation and Drainage Engineering 131: 525–532.).
For the specific case in which head losses are estimated for a discrete flow distribution throughout a multioutlet lateral or manifold, a reduction coefficient F or Christiansen (1942)Christiansen, J.E. 1942. Irrigation by Sprinkling. UC Davis, Davis, CA, USA. (California Agricultural Experimental Station Bulletin, 670). factor is introduced into head losses by the formula:
where: C = head loss coefficient varying with water temperature, which is assumed ranging from 0.516 to 0.381 for 5 °C < T < 60 °C (Montalvo, 2007Montalvo, T. 2007. Trickle irrigation: Installations Design = Riego Localizado: Diseño de Instalaciones. Ediciones VJ, Valencia, Spain (in Spanish).); D = inside diameter (mm); K_{m} = increasing coefficient accounting for minor head losses; F = Christiansen reduction factor (Christiansen, 1942Christiansen, J.E. 1942. Irrigation by Sprinkling. UC Davis, Davis, CA, USA. (California Agricultural Experimental Station Bulletin, 670).); L= lateral or manifold length; Q = initial flow rate.
In this sense, the insertion of each emitter to the lateral and the connection to manifold introduce some additional losses that should be considered. Several researchers carried out different approaches to estimate this local loss in microirrigation laterals and proposed a constant equivalent length (l_{e}) for emitter connection to estimate local loss. In the last decade, researchers proposed new approaches to determine local losses (Provenzano et al., 2007Provenzano, G.; Di Dio, P.; PalauSalvador, G. 2007. New computational fluid dynamic procedure to estimate friction and local losses in coextruded drip laterals. Journal of Irrigation and Drainage Engineering 133: 520527, ^{2014}Provenzano, G.; Di Dio, P.; Leone, R. 2014. Assessing a local losses evaluation procedure for lowpressure layflat drip laterals. Journal of Irrigation and Drainage Engineering 140: 04014017.; Wang et al., 2016Wang, Y.; Zhu, D.; Chen, J.; Zheng, C.; Liu, Y. 2016. Numerical study on minor head loss of online emitters in drip irrigation laterals. Journal of Hydroelectric Engineering 35: 2231.; Ding et al., 2016Ding, F.; Wang. W.; Hu, X. 2016. Prediction model for head loss in drip irrigation laterals based on dimensional analysis. Journal of Hydroelectric Engineering 35: 6774.). DIMSUB allows establishing the value of l_{e} and calculating an increasing coefficient to estimate local losses:
where: s = emitter spacing; l_{e}= equivalent length due to local head losses.
Consequently, iterative calculation methods help to determine the minimum pressure point along the lateral length and solve the energy gradient line produced in a lateral or manifold (Arviza and Palau, 2016Arviza, J.; Palau, C.V. 2016. Determination of the feeding point of an irrigation subunit with slope = Determinación del punto de alimentación de una subunidad de riego con desnivel. Available at: http://hdl.handle.net/10251/68355 [Accessed July 28, 2016] (in Spanish).
http://hdl.handle.net/10251/68355...
). Design and sizing of subunits are based on emitter uniformity criteria by means of controlling their allowable pressure head variation.
Correct and precise design requires complex calculation procedures (Dandy and Hassanli, 1996Dandy, G.C.; Hassanly, A.M. 1996. Optimum design and operation of multiple subunit drip irrigation systems. Journal of Irrigation and Drainage Engineering 122: 165175.; Lamm et al., 2007Lamm, F.R.; Ayars, J.E.; Nakayama, F.S. 2007. Microirrigation for Crop Production: Design, Operation, and Management. Elsevier, Amsterdam, The Netherlands.; Ravindra et al., 2008Ravindra, V.K.; Singh, R.P.; Mahar, P.S. 2008. Optimal design of pressurized irrigation subunit. Journal of Irrigation and Drainage Engineering 134:137146.). Given the complexity of the design of the microirrigation subunit, in many cases sizing has been oversimplified to the detriment of precision and oversizing subunits. In many cases, even nowadays, many engineers use tables supplied by manufacturers to determine the maximum lateral lengths and neglect calculating the size of the manifold or the required initial pressure of the subunit.
DIMSUB was developed to provide freeware to engineers and technicians to design and size irrigation subunits appropriately and rapidly.
Microirrigation optimum design using DIMSUB
DIMSUB is a decision tool for designing microirrigation subunits on Microsoft Excel® by developing a series of user forms on Visual Basic for Applications.
Optimum design and sizing procedures start with specifying the possible maximum pressure head variation, ΔH_{sub}, in the irrigation system to ensure adequate uniformity (Valiantzas, 2003Valiantzas, J.D. 2003. Explicit hydraulic design of microirrigation submains units with tapered manifold and laterals. Journal of Irrigation and Drainage Engineering 129: 227236.; Yildirim, 2008Yildirim, G. 2008. Determining operating inlet pressure head incorporating uniformity parameters for multioutlet plastic pipelines. Journal of Irrigation and Drainage Engineering 134: 341348.). The slope, terrain geometry, emitter type (nominal flow rate and hydraulic characteristic curve), lateral diameter and emitter spacing allow the analysis of the optimal microirrigation system design. It can also be used to estimate the initial pressure and flow rate of a subunit to keep the energy requirements as low as possible. Figure 2 shows DIMSUB synthesized by the stepbystep method for the hydraulic design of subunits.
Stepbystep method for irrigation subunits design using NPC emitters (A) and PC emitters (B). F is the Christiansen's reduction factor; L is the length of manifold (m); Q is the initial flow rate (L h^{−1}); Km is the increasing coefficient for local losses; h_{adm} are admissible head losses in manifold; α and β are adjustment coefficients to obtain average inlet pressure.
Maximum pressure head variation in the subunit for the nonpressurecompensating emitter, ΔH_{subunit}, can be estimated by x, discharge exponent in the emitter hydraulic characteristic, and
The size of the manifold pipe depends on the remaining maximum head loss after designing the laterals, as shown in the following equations (Keller and Bliesner, 1990Keller, J.; Bliesner, R.D. 1990. Sprinkle and Trickle Irrigation. Van Nostrand Reinhold, New York, NY, USA.):
where: ΔH= maximum pressure head variation in the subunit, lateral or manifold, respectively (m); h_{max manifold} =maximum pressure head loss in manifold (m); ΔZ = elevation difference (m).
At this point, introducing all data related to the subunit, DIMSUB is able to calculate the remaining maximum head loss for the manifold (eq. 8) after designing the laterals and the minimum inside diameter (D_{min}) to keep the pressure head variation considered.
Additionally, this tool allows estimating the pressure head required at the subunit entrance, P/γ, considering head losses and subunit slopes to ensure flow uniformity. Keller and Karmeli (1975)Keller, J.; Karmeli, D. 1975 Trickle Irrigation Design Parameters. Rain Bird Sprinkler Manufacturing, Glendora, CA, USA., in a study to achieve a better flow distribution, introduced α and β coefficients to adjust initial pressure head in multiple outlet pipes for nonpressurecompensating emitters. Keller and Bliesner (1990)Keller, J.; Bliesner, R.D. 1990. Sprinkle and Trickle Irrigation. Van Nostrand Reinhold, New York, NY, USA. proposed the formula to obtain the average head along the pipe. Montalvo (2007)Montalvo, T. 2007. Trickle irrigation: Installations Design = Riego Localizado: Diseño de Instalaciones. Ediciones VJ, Valencia, Spain (in Spanish). justified adequate values for α and β, correcting difference in elevation and head losses along laterals and manifold. Values for each type of pipe configuration are calculated by DIMSUB considering that pressurecompensating emitters do not require this adjustment (Figure 2).
DIMSUB computer application estimates the emission uniformity coefficient (EU) in the microirrigation system with the following expression proposed by Keller and Karmeli (1975)Keller, J.; Karmeli, D. 1975 Trickle Irrigation Design Parameters. Rain Bird Sprinkler Manufacturing, Glendora, CA, USA., and used later by ASAE (2005)American Society of Agricultural and Engineers [ASAE]. 2005. EP405.1: Design and Installation of Microirrigation Systems. ASAE, St. Joseph, MI, USA. with designing purposes:
where: n=number of emitters per plant; q_{min} and
In order to demonstrate versatility and easiness of the DIMSUB program, two case studies are described as carried out in a citrus orchard and an almond orchard in Valencia (Spain) (Figures 3A, 3B).
Case studies. A) Regular citrus orchard in Picassent; B) Irregular geometry of almond orchard in Chulilla (Valencia, Spain).
Citrus case study
The citrus orchard in Picassent (Lat. 39°20’59.5″ N, Long. 0°28’16.67″ W, altitude 58.5 m), Valencia (Spain), covers approximately 0.96 ha with a rectangular geometry of 128 × 75 m (Figure 3A).
For this example, a commercial dripperline with outside diameter (OD) of 16 mm (Internal diameter, ID = 14.2 mm) and 3.5 L h^{−1} emitters are used (coefficient of variation, CV <7 %, ISO, 2004). Two laterals per plant row with emitter spacing of 1.0 m were installed. The laterals are on a downhill slope of −1 % while the manifold is on the level (0 % slope). In this example, an equivalent length of 0.25 m due to lateral local losses is considered according to manufacturer's recommendation.
Different alternatives were considered with different emitter types and feeding points to analyze the operating mode of the DIMSUB application (Table 1).
The first three alternatives use PC emitters to achieve wider pressure variations, since the emitter flow rate is always close to the nominal flow rate. In this case, engineers establish the allowable pressure head variation for the subunit considering the terrain slope, lengths of laterals and manifolds and distance between emitters. The range of pressure variation can influence in the initial pressure head and consequently energy consumption by the irrigation system.
The last three alternatives include NPC emitters, whose flow rate depends on their operating pressure, for pressure variation, which depends on the x exponent of the emitter, to be limited to achieve acceptably uniform irrigation. For instance, a nonpressurecompensating emitter with a discharge exponent around x = 0.5, with working pressure of 10.0 m (CV < 7 %), admits pressure head variations from 2.0 to 3.0 m to control the maximum flow rate variation within 10 to 15 %, respectively (Wu and Gitlin, 1983Wu, I.P.; Gitlin, H.M. 1983. Drip irrigation application efficiency and schedules. Transactions of the ASAE 26: 9299.).
Computer application allows hydraulic variables to be selected and calculated from different feeding points in either the manifold or lateral pipes (Figure 4). The manifold and laterals can be fed from the extreme point, from a midpoint in flat terrain with no slope, or from an intermediate point in pipes with certain slope. Recently, some approaches have been proposed to find the optimal feeding point. Baiamonte et al. (2015)Baiamonte, G.; Provenzano, G.; Rallo, G. 2015. Analytical approach determining the optimal length of paired drip laterals in uniformly sloped fields. Journal of Irrigation and Drainage Engineering 141: 04014042., Ju et al. (2017)– Ju, K.L; Wu, P.T.; Weckler, P.R.; Zhu, D. 2017. Simplified approach for designing length of microirrigation laterals. Applied Engineering in Agriculture 33: 7582. and Monserrat et al. (2018)Monserrat, J.; Barragan, J.; Cots, L. 2018. Design of paired laterals on uniformly sloping fields. Journal of Irrigation and Drainage 144: 04018008. presented analytical approaches determining the optimal length of paired drip laterals. DIMSUB with a stepbystep procedure calculates the required uphill and downhill pipe lengths to balance pressures in all emitters or laterals (Arviza and Palau, 2016Arviza, J.; Palau, C.V. 2016. Determination of the feeding point of an irrigation subunit with slope = Determinación del punto de alimentación de una subunidad de riego con desnivel. Available at: http://hdl.handle.net/10251/68355 [Accessed July 28, 2016] (in Spanish).
http://hdl.handle.net/10251/68355...
).
Subunit feeding points: A) Laterals and manifold fed from an extreme point; B) Laterals fed from extreme point and manifold fed from midpoint; C) Laterals fed from intermediate point and manifold fed from extreme point.
The pressurecompensating interval for PC emitters varies for each manufacturer, usually ranging from 5.0 to 40.0 m. This interval allows a maximum pressure head variation of 35.0 m to ensure flow uniformity. However, a 5.0 m pressure head variation was established as design criteria, considering the plot characteristics, in alternatives PC1, PC2 and PC3, in order to get an adequate initial pressure head at the beginning of the subunit. Nevertheless, this value can be adjusted in DIMSUB establishing a minimum and maximum inlet operating pressure, as shown in DIMSUB (Figure 5).
DIMSUB hydraulic data display. PC (Pressurecompensating) emitter, 3.5 L h^{−1} (pressure range 5.0 – 40.0 m).
Figure 5 shows the data entered for the first alternative, for which the program provides the hydraulic results for the irrigation lateral in box 1, including: initial flow rate, Christiansen coefficient, lateral head losses, pressure variation, α and β coefficients, required initial pressure, final extreme pressure, point of minimum pressure head and emission uniformity (EU) coefficient (eq. 9). Box 2 provides the lateral main calculation variables, such as number of emitters, lateral slope, maximum head losses admitted for manifold, local loss coefficient of the drippers.
The next step is to input the parameters required for designing and sizing the manifold.
The maximum pressure head variation established for the subunit (ΔH_{subunit}) is then shared between the lateral and manifold head losses and slopes (eq. 7). DIMSUB begins by calculating lateral parameters. According to lateral head losses and differences in elevation, the remaining energy is available for manifold sizing.
In order to estimate the manifold length, manifold data are inserted, including slope of the terrain, number of plant rows, distance between laterals in the rows or distance from laterals of adjacent rows, as indicated in box 3 (Figure 5).
Based on the admissible pressure head variation in the subunit and laterals, DIMSUB calculates the parameters required to design the manifold pipe (box 4) (Figure 5): maximum pressure variation and maximum head loss remaining for the manifold (eq. 8), flow rate per lateral, total length and slope. It also estimates the minimum pipe diameter required to limit head losses within the admissible value.
Finally, box 5 (Figure 5) shows the results for the manifold pipe: nominal pipe diameter, internal diameter, real head loss, required initial pressure head, final pressure and pressure head variation in the manifold.
Almond orchard case study
This orchard is in Chulilla, Valencia (Spain) (Lat. 39° 56’ 56.73″N, Longitude 0° 50’ 3.25″W, altitude 387.8 m), and covers an area of approximately 0.75 ha. Figure 3B shows that the subunit geometry is irregular and the lateral length is variable as it follows the plot shape.
Designing and sizing irregular landscapes requires a special treatment, which can be carried out by outlining terrain geometry from an image or orthophoto in DIMSUB. The manifold network is graphically defined thus as water inlet points to the laterals and final points of laterals of different lengths (Figure 6).
The study considered a commercial lateral having outside diameter (OD) of 16 mm (Internal diameter, ID = 14.2 mm) with an integral pressurecompensating emitters discharging 3.5 L h^{−1}. The emitter spacing is 1.0 m with a double lateral per plant row. It is assumed that the insertion of each emitter produces the head losses corresponding to an additional pipe equivalent length of 0.25 m. This orchard has a 2.2 % uphill slope along the line of almond trees and −1.4 % crosswise from the water hydrant. The feeding point was from the extreme in the manifold and laterals (Figure 3B). Admissible pressure head variation was assumed equal to 5.0 m in all the alternatives considered.
This case of study evaluates the design results adjusting to the actual irregular shape (IS) of the subunit and those obtained with a simplified regular plot shape (RS). IS1 and IS2 use DIMSUB irregular subunit scheme to solve the case (Figure 6). RS1 and RS2 consider a simplified procedure designing a rectangular subunit of 104 × 107 m. Alternatives IS2 and RS2 are improved by a twodiameter tapered manifold.
Results and Discussion
The proposed computer application is capable of efficiently analyzing different solutions for a given irrigation system and allowing the choice for the best pricewise and technical alternative.
DIMSUB can handle data files from the user form, as well as visualize a scheme of any type of lateral and manifold feed. An Excel sheet with the results can be printed and exported for further analysis of the irrigation distribution network. Moreover, the computer program measures the total lateral and manifold length and number of emitters for each subunit design analyzed and calculates the corresponding total cost in pipe material.
Results of the hydraulic analysis and cost of material are shown in Table 2 and Table 3, respectively, for all alternatives proposed in the case of the citrus orchard subunit.
This case study shows that, based on the subunit conditions, the use of noncompensating emitters offers better hydraulic and economic conditions. Lower initial subunit pressure is required (between 10.7 and 12.6 m) and the cost of laterals and emitters is considerably lower (44 % of cost reduction) (Tables 2 and 3).
In the second case, irregular and regular geometries are studied in an almond orchard. The first cases analyzed (Table 4), IS1 and IS2, show the results obtained from the DIMSUB tool for designing irregular plots. In the regular geometry study, sizing is done with a single pipe size, while in the regular study, design is improved by adjusting the diameters with a twodiameter tapered manifold and achieving a more economical solution with lower energy requirements for the same established pressure variation.
This study shows that the application can deal with irregularly shaped plots by means of an image file that must be scaled to enable designing the subunit.
The results proved that adjusting the design to the terrain geometry appreciably improves the solution and reduces costs, due to the adjustment of manifold pipe sizes. Furthermore, multioutlet tapered manifolds in long pipe runs, in both IR2 and RS2, achieve lower energy and cost requirements at the subunit entrance.
Conclusions
DIMSUB is a computer application that can be used as a decision support tool. This program can aid in designing and sizing precisely and rapidly individual irrigation subunits. The lateral and manifold pipes in a subunit can be graphically drawn for any type of geometry.
Irregular geometries, sloped terrains and different feeding points of laterals and manifold can be easily studied to achieve accurate results and evaluate design alternatives.
Case studies have been evaluated with DIMSUB, showing the main performance variables of irrigation systems that allow technicians to make decisions. Results obtained from this application have been supported by real cases where engineers designed and executed irrigation subunits successfully.
References
 American Society of Agricultural and Engineers [ASAE]. 2005. EP405.1: Design and Installation of Microirrigation Systems. ASAE, St. Joseph, MI, USA.
 Arviza, J.; Palau, C.V. 2016. Determination of the feeding point of an irrigation subunit with slope = Determinación del punto de alimentación de una subunidad de riego con desnivel. Available at: http://hdl.handle.net/10251/68355 [Accessed July 28, 2016] (in Spanish).
» http://hdl.handle.net/10251/68355  Bagarello, V.; Ferro, V.; Provenzano, G.; Pumo, D. 1995. Experimental study on flow resistance law for smalldiameter plastic pipes. Journal of Irrigation and Drainage Engineering 121: 313–316.
 Baiamonte, G.; Provenzano, G.; Rallo, G. 2015. Analytical approach determining the optimal length of paired drip laterals in uniformly sloped fields. Journal of Irrigation and Drainage Engineering 141: 04014042.
 Blasius, H. 1913. The law of similarity for friction in fluids = Das aehnlichkeitsgesetz bei reibungsvorgangen in flussigkeiten. Forshhungsarbeiten auf dem Gebiete des Ingenieurwesen 131(1) (in German).
 Bralts, V.F.; Segerlind, L.J. 1985. Finite elements analysis of drip irrigation submain unit. Transactions of ASAE 28: 809–814.
 Bralts, V.F.; Driscoll, M.A.; Shayya, W.H.; Cao, L. 1993. An expert system for the hydraulic analysis of microirrigation systems. Computers and Electronics in Agriculture 9: 275287.
 Carrión, F.; Tarjuelo, J.; Hernandez, D.; Moreno, M. 2013. Design of microirrigation subunit of minimum cost with proper operation. Irrigation Science 31: 11991211.
 Carrión, F.; Montero, J.; Tarjuelo, J.M.; Moreno, M. 2014. Design of sprinkler irrigation subunit of minimum cost with proper operation: application at corn crop in Spain. Water Resources Management 28: 50735089.
 Christiansen, J.E. 1942. Irrigation by Sprinkling. UC Davis, Davis, CA, USA. (California Agricultural Experimental Station Bulletin, 670).
 Dandy, G.C.; Hassanly, A.M. 1996. Optimum design and operation of multiple subunit drip irrigation systems. Journal of Irrigation and Drainage Engineering 122: 165175.
 Ding, F.; Wang. W.; Hu, X. 2016. Prediction model for head loss in drip irrigation laterals based on dimensional analysis. Journal of Hydroelectric Engineering 35: 6774.
 Hoffman, G.J.; Evans, R.G.; Marvin, E.J.; Derrel, L.M.; Elliott, R.L. 2007. Design and Operation of Farm Irrigation Systems. 2ed. ASABE, St. Joseph, MI, USA.
 – Ju, K.L; Wu, P.T.; Weckler, P.R.; Zhu, D. 2017. Simplified approach for designing length of microirrigation laterals. Applied Engineering in Agriculture 33: 7582.
 Keller, J.; Bliesner, R.D. 1990. Sprinkle and Trickle Irrigation. Van Nostrand Reinhold, New York, NY, USA.
 Keller, J.; Karmeli, D. 1975 Trickle Irrigation Design Parameters. Rain Bird Sprinkler Manufacturing, Glendora, CA, USA.
 Lamm, F.R.; Ayars, J.E.; Nakayama, F.S. 2007. Microirrigation for Crop Production: Design, Operation, and Management. Elsevier, Amsterdam, The Netherlands.
 Myers, L.E.; Bucks D.A. 1972. Uniform irrigation with low pressure trickle system. Journal of Irrigation and Drainage 98: 341346.
 Monserrat, J.; Barragan, J.; Cots, L. 2018. Design of paired laterals on uniformly sloping fields. Journal of Irrigation and Drainage 144: 04018008.
 Montalvo, T. 2007. Trickle irrigation: Installations Design = Riego Localizado: Diseño de Instalaciones. Ediciones VJ, Valencia, Spain (in Spanish).
 Pedras, C.M.G.; Pereira, L.S.; Gonçalves, J.M. 2009. MIRRIG: a decision support system for design and evaluation of microirrigation systems. Agricultural Water Management 96: 691701.
 Pereira, L.S.; Oweis, T.; Zairi, A. 2002. Irrigation management under water scarcity. Agricultural Water Management 57: 175206.
 Pitts, D.J.; Zajueta, F.S.; Smajstrla, A.G. 1986. Micro Irrigation System Submain Design Evaluation. University of Florida, Gainesville, FL, USA.
 Provenzano, G.; Pumo, D.; Di Dio, P. 2005. Simplified procedure to evaluate head losses in drip irrigation laterals. Journal of Irrigation and Drainage Engineering 131: 525–532.
 Provenzano, G.; Di Dio, P.; PalauSalvador, G. 2007. New computational fluid dynamic procedure to estimate friction and local losses in coextruded drip laterals. Journal of Irrigation and Drainage Engineering 133: 520527
 Provenzano, G.; Di Dio, P.; Leone, R. 2014. Assessing a local losses evaluation procedure for lowpressure layflat drip laterals. Journal of Irrigation and Drainage Engineering 140: 04014017.
 Ravindra, V.K.; Singh, R.P.; Mahar, P.S. 2008. Optimal design of pressurized irrigation subunit. Journal of Irrigation and Drainage Engineering 134:137146.
 Valiantzas, J.D. 2003. Explicit hydraulic design of microirrigation submains units with tapered manifold and laterals. Journal of Irrigation and Drainage Engineering 129: 227236.
 Von Bernuth, R.D. 1990. Simple and accurate friction loss equation for plastic pipe. Journal of Irrigation and Drainage Engineering 116: 294–298.
 Wang, Y.; Zhu, D.; Chen, J.; Zheng, C.; Liu, Y. 2016. Numerical study on minor head loss of online emitters in drip irrigation laterals. Journal of Hydroelectric Engineering 35: 2231.
 Wu, I.P.; Gitlin, H.M. 1974. Drip irrigation design based on uniformity. Transactions of the ASAE 17: 429432.
 Wu, I.P.; Gitlin, H.M. 1975. Energy gradient line for drip irrigation laterals. Journal of Irrigation and Drainage 101: 323–326
 Wu, I.P.; Gitlin, H.M. 1983. Drip irrigation application efficiency and schedules. Transactions of the ASAE 26: 9299.
 Wu, I.P. 1997. An assessment of hydraulic design of microirrigation systems. Agricultural Water Management 32: 275284.
 Wu, I.P.; Barragan, J. 2000. Design criteria for microirrigation systems. Transactions of the ASAE 43: 11451154.
 Yildirim, G. 2008. Determining operating inlet pressure head incorporating uniformity parameters for multioutlet plastic pipelines. Journal of Irrigation and Drainage Engineering 134: 341348.
 Zhang, L.; Wu, P.; Zhu, D. 2013. Hydraulic design procedure for drip irrigation submain unit based on relative flow difference. Irrigation Science 31: 10651073.
Appendix
DIMSUB is available by contacting jarviza@agf.upv.es.
Publication Dates

Publication in this collection
05 Sept 2019 
Date of issue
2020
History

Received
28 June 2018 
Accepted
22 Oct 2018