Acessibilidade / Reportar erro

Characterization of venturi injector using dimensional analysis

Caracterização de injetor venturi utilizando análise dimensional

ABSTRACT

Venturi injectors are commonly employed for fertigation purposes in agriculture, in which they draw fertilizer from a tank into the irrigation pipeline. The knowledge of the amount of liquid injected by this device is used to ensure an adequate fertigation operation and management. The objectives of this research were (1) to carry out functional tests of Venturi injectors following requirements stated by ISO 15873; and (2) to model the injection rate using dimensional analysis by the Buckingham Pi theorem. Four models of Venturi injectors were submitted to functional tests using clean water as motive and injected fluid. A general model for predicting injection flow rate was proposed and validated. In this model, the injection flow rate depends on the fluid properties, operating hydraulic conditions and geometrical characteristics of the Venturi injector. Another model for estimating motive flow rate as a function of inlet pressure and differential pressure was adjusted and validated for each size of Venturi injector. Finally, an example of an application was presented. The Venturi injector size was selected to fulfill the requirements of the application and the operating conditions were estimated using the proposed models.

Key words:
Buckinghan Pi theorem; irrigation equipment; fertigation

RESUMO

Injetores Venturi são comumente empregados para fins de fertirrigação na agricultura, extraindo fertilizante de um tanque para a linha de irrigação. O conhecimento da quantidade de líquido a ser injetado por este dispositivo é útil para assegurar operação e manejo adequados da fertirrigação. Os objetivos desta pesquisa foram: (1) conduzir testes funcionais de injetores Venturi seguindo as normas estabelecidas pela ISO 15873; e (2) modelar a vazão de injeção utilizando Análise Dimensional pelo teorema de Buckingham Pi. Quatro modelos de injetores Venturi foram submetidos a testes operacionais usando água limpa como fluido motriz e como fluido succionado. Um modelo geral para predição da vazão de injeção foi proposto e validado, sendo este dependente das propriedades do fluido, das condições hidráulicas operacionais e das características geométricas do injetor. Outro modelo para estimar a vazão motriz, em função da pressão de entrada e da pressão diferencial foi ajustado e validado para cada tamanho de Venturi. Finalmente, um exemplo de aplicação foi apresentado. O tamanho do injetor Venturi foi selecionado para atender aos requisitos da aplicação e as condições de operação foram estimadas utilizando os modelos propostos.

Palavras-chave:
teorema de Buckinghan Pi; equipamento de irrigação; fertirrigação

Introduction

Venturi tubes are commonly employed for fertigation purposes in agriculture (Frizzone et al., 2012Frizzone, J. A.; Freitas, P. S. L. de; Rezende, R.; Faria, M. A. de. Microirrigação: Gotejamento e microaspersão. Maringá: EDUEM, 2012. 356p.). Venturi used as injectors rely on the Venturi pressure drop principle to draw chemicals from a stock tank into the irrigation pipeline (Hoffman et al., 2007Hoffman, G. J.; Evans, R. G.; Jensen, M. E.; Martin, D. L.; Elliott, R. L. Design and operation of farm irrigation systems. 2.ed. St. Joseph: American Society of Agricultural and Biological Engineers, 2007. 863p.).

The proper design, operation, and management of fertigation may enable uniform application of fertilizers and increase crop yields or their quality. Proper fertigation may also save fertilizers and reduce water usage by directly delivering water and nutrient near to crop roots, which boost economic and energy profitability to farmers (Sinha et al., 2017Sinha, I.; Buttar, G. S.; Brar, A. S. Drip irrigation and fertigation improve economics, water and energy productivity of spring sunflower (Helianthus annuus L.) in Indian Punjab. Agricultural Water Management, v.185, p.58-64, 2017. https://doi.org/10.1016/j.agwat.2017.02.008
https://doi.org/10.1016/j.agwat.2017.02....
; Tian et al., 2017Tian, D.; Zhang, Y.; Zhou, Y.; Mu, Y.; Liu, J.; Zhang, C.; Liu, P. Effect of nitrification inhibitors on mitigating N2O and NO emissions from an agricultural field under drip fertigation in the North China Plain. Science of the Total Environment, v.598, p.87-96, 2017. https://doi.org/10.1016/j.scitotenv.2017.03.220
https://doi.org/10.1016/j.scitotenv.2017...
). However, inappropriate fertigation may reduce crop growth and yield, as well as may lead to environmental problems and contamination hazards related to soil, groundwater and surface water (Silva et al., 2013Silva, P. F. da; Lima, C. J. G. de S.; Barros, A. C.; Silva, E. M. da; Duarte, S. N. Sais fertilizantes e manejo da fertirrigação na produção de tomateiro cultivado em ambiente protegido. Revista Brasileira de Engenharia Agrícola e Ambiental, v.17, p.1173-1180, 2013. https://doi.org/10.1590/S1415-43662013001100007
https://doi.org/10.1590/S1415-4366201300...
).

Several papers have studied the design, performance and modeling of Venturi injectors for agricultural purposes (Feitosa Filho et al., 1998Feitosa Filho, J. C.; Botrel, T. A.; Pinto, J. M. Desempenho de injetores tipo venturi em condições de descarga livre e pressurizada. Engenharia Agrícola, v.17, p.20-35, 1998., 1999Feitosa Filho, J. C.; Pinto, J. M.; Arruda, N. T. de. Dimensionamento, construção e características hidráulicas de um injetor tipo venturi para uso na quimigação. Irriga, v.4, p.43-53, 1999. https://doi.org/10.15809/irriga.1999v4n2p43-53
https://doi.org/10.15809/irriga.1999v4n2...
; Lima Neto & Porto, 2004Lima Neto, I. E.; Porto, R. M. Performance of low-cost ejectors. Journal of Irrigation and Drainage Engineering, v.130, p.122-128, 2004. https://doi.org/10.1061/(ASCE)0733-9437(2004)130:2(122)
https://doi.org/10.1061/(ASCE)0733-9437(...
). Techniques of computational fluid dynamics (CFD) simulations are also a useful tool to know the flow behavior inside a Venturi injector (Kuldeep & Saharan, 2016Kuldeep, R.; Saharan, V. K. Computational study of different venturi and orifice type hydrodynamic cavitating devices. Journal of Hydrodynamics, v.28, p.293-305, 2016. https://doi.org/10.1016/S1001-6058(16)60631-5
https://doi.org/10.1016/S1001-6058(16)60...
; Manzano et al., 2016Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V.; Vasconcelos, D. V. Geometry and head loss in Venturi injectors though computational fluid dynamics. Engenharia Agrícola, v.36, p.482-491, 2016. https://doi.org/10.1590/1809-4430-Eng.Agric.v36n3p482-491/2016
https://doi.org/10.1590/1809-4430-Eng.Ag...
). However, in these papers, each model of evaluated Venturi injector had a corresponding regression model.

The amount of liquid fertilizer injected by a Venturi tube is influenced by its design, operational parameters, and water properties. A few prediction models for estimating the injection rate are found in the literature, but they are limited for use in gas applications (Liu et al., 2014Liu, K.; Lu, H.; Guo, X.; Sun, X.; Tao, S.; Gong, X. Experimental study on flow characteristics and pressure drop of gas-coal mixture through Venturi. Power Technology, v.268, p.401-411, 2014. https://doi.org/10.1016/j.powtec.2014.07.019
https://doi.org/10.1016/j.powtec.2014.07...
; Xu et al., 2015Xu, Y.; Zhang, Q.; Zhang, T.; Ba, X. An overreading model for nonstandard Venturi meters based on H correction factor. Measurement, v.61, p.100-106, 2015. https://doi.org/10.1016/j.measurement.2014.10.028
https://doi.org/10.1016/j.measurement.20...
; Gupta et al., 2016Gupta, B.; Nayak, A. K.; Kandar, T. K.; Nair, S. Investigation of air-water two phase flow through a Venturi. Experimental Thermal and Fluid Science, v.70, p.148-154, 2016. https://doi.org/10.1016/j.expthermflusci.2015.07.012
https://doi.org/10.1016/j.expthermflusci...
). Moreover, no general models were found in the literature to predict the injection rate of Venturi tubes used for agricultural purposes. ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.) specifies construction and operational requirements, as well as test methods for Venturi tubes employed as liquid injectors in irrigation systems. Functional tests described by that standard do not mention any models for estimation of injection rate.

The objectives of this research were (1) to carry out functional tests of Venturi tubes following requirements stated by ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.); (2) to model the injection rate of Venturi tubes for practical purposes using dimensional analysis by the Buckingham Pi theorem. In practice, the functional tests are required to determine operational characteristics of Venturi injectors, that are necessary to set variables and obtain a target injection flow. Results of those tests basically enable to relate motive flow, inlet pressure, differential pressure and injection flow. In addition, a general model is useful because a single equation enables to relate dimensions and operational characteristics of a range of Venturi sizes operating within a range of conditions.

Material and Methods

Four models of Venturi injectors were submitted to the functional tests stated by ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.) at the Irrigation Testing Laboratory (LEMI) of the Escola Superior de Agricultura Luiz de Queiroz (ESALQ/USP), Piracicaba, São Paulo State, Brazil (22° 42’ S; 47° 30’ W; and altitude of 546 m). Table 1 shows the geometric characteristics of the evaluated devices. Models 1, 2 and 3 were manufactured by Amanco Ltda. while Mondo Irrigação manufactured Venturi Model 4, being both Brazilian companies. The four models represent the commercial sizes of Venturi tubes usually employed for fertigation purposes in Brazil and all of them had a check valve integrated at their injection inlet. The check valve is designed to avoid reverse flow through the injection pipe and is a simple device consisting of a spring, plastic sphere, and a rubber sealing ring. As demonstrated by Santos et al. (2012Santos, L. da C.; Zocoler, J. L.; Justi, A. L.; Silva, A. O.; Correia, J. de S. Estudo comparativo da taxa de injeção em injetor do tipo Venturi com e sem válvula de retenção. Irriga, v.edição especial, p.145-154, 2012. https://doi.org/10.15809/irriga.2012v1n01p145
https://doi.org/10.15809/irriga.2012v1n0...
), when check valves were integrated to Venturi injectors, their influence on hydraulic operational characteristics was not relevant for fertigation purposes. Figure 1 illustrates the geometric characteristics of an ordinary Venturi injector.

Table 1
Geometrical characteristics of evaluated Venturi injectors: nominal outlet diameter (DN), internal inlet diameter (Din), throat diameter (Dt), and injection pipe diameter (Dinj)

Figure 1
Main geometric characteristics of a Venturi injector

Tests were undertaken in a testing bench (Figure 2) that consisted of a pump and an automated system for pressure control and equipped with a variable frequency drive (VFD) and proportional-integral-derivative (PID) controller. The PID controller used a pressure transmitter, installed at the inlet of the Venturi tube, to calculate the error between the pressure set point and the measured pressure, and then applied an output signal to the variable frequency drive.

Figure 2
Instrumentation diagram of the testing bench used for functional tests of Venturi injectors: (TT) Temperature transmitter; (FT) flow transmitter; (PT) pressure transmitter; (PDT) differential pressure transmitter

Each injector was mounted in-line, with the full main flow introduced as motive water at the injector inlet ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.). Clean water was used as motive and injected fluid to comply with ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.). Although experiments evaluating various fluids do not belong to the scope of this research, literature indicates that viscosity and density variation between motive fluid and injected fluid affects operational characteristics of Venturi injectors (Yuan et al., 2000Yuan, Z.; Choi, C. Y. Waller, P. M.; Colaizzi, P. Effects of liquid temperature and viscosity on venturi injectors. Transactions of the American Society Agricultural Engineers, v.43, p.1441-1447, 2000. https://doi.org/10.13031/2013.3042
https://doi.org/10.13031/2013.3042...
; Kumar et al., 2012Kumar, M.; Rajput, T. B. S.; Patel, N. Effect of system pressure and solute concentration on fertilizer injection rate of a venturi for fertigation. Journal of Agricultural Engineering, v.49, p.9-13, 2012.).

ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.) also states that the relative elevations of the suction port and the surface of the suction water must remain constant during the test. During all experiments, the vertical distance between the surface of the suctioned water and the Venturi was set to 1 m (see Δz in Figure 1). Variations of such vertical distances are beyond the scope of this research. Each device was evaluated at five inlet pressures ranging from 100 to 300 kPa in intervals of 50 kPa. The range of inlet pressures was defined based on practical values commonly used in microirrigation and low-pressure sprinkler systems. For each inlet pressure, the motive flow was set in order to provide differential pressures of 20, 40, 60, 80 and 100% in relation to the inlet pressure ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.).

A pressure transmitter (measurement ranging from 0 to 500 kPa, resolution of 0.1 kPa, maximum error of 0.1% of full scale) was installed at the inlet of the Venturi to monitor the inlet pressure and to provide the feedback signal required by the PID controller. A differential pressure transmitter (measurement ranging from 0 to 500 kPa, resolution of 0.01 kPa, and maximum error of 0.1% of full scale) monitored the differential pressure between the inlet and outlet of the Venturi. Two models of electromagnetic flowmeters were used for measuring motive flow rates. The motive flow of Venturi Models 1 and 2 was measured using a flowmeter of measurement ranging from 0 to 4 m3 h-1, while for the other models of Venturi a flowmeter of measurement ranging from 0 to 10 m3 h-1 was employed. Both flowmeters have a resolution of 0.001 m3 h-1 and expanded uncertainty of 0.5% in relation to their full range.

The injection rate was measured by the gravimetric method using a precision scale (measurement ranging from 0 to 100 kg, resolution of 0.01 kg) and a stopwatch. Water temperature was monitored using a temperature transmitter of measurement ranging from 0 to 50 ºC, resolution of 0.1 ºC, and maximum error of 0.1% of full scale. A proportional valve was used for adjusting the testing flow rate (Figure 2).

Models based on dimensional analysis were developed to estimate pressure losses in drip irrigation laterals (Perboni et al., 2015Perboni, A.; Frizzone J. A.; Camargo, A. P. de; Pinto, M. F. Modelling head loss along emitting pipes using dimensional analysis. Engenharia Agrícola, v.35, p.442-457, 2015. https://doi.org/10.1590/1809-4430-Eng.Agric.v35n3p442-457/2015
https://doi.org/10.1590/1809-4430-Eng.Ag...
), pressure losses in filters (Duran-Ros et al., 2010Duran-Ros, M.; Arbat, G.; Barragán, J.; Ramírez de Cartagena, F.; Puig-Bargués, J. Assessment of head loss equations developed with dimensional analysis for micro irrigation filters using effluents. Biosystems Engineering, v.106, p.521-526, 2010. https://doi.org/10.1016/j.biosystemseng.2010.06.001
https://doi.org/10.1016/j.biosystemseng....
; Wu et al., 2014Wu, W.; Chen, W. E. I.; Liu, H; Yin, S.; Niu, Y. A new model for head loss assessment of screen filters developed with dimensional analysis in drip irrigation systems. Irrigation and Drainage, v.63, p.523-531, 2014. https://doi.org/10.1002/ird.1846
https://doi.org/10.1002/ird.1846...
), minor losses in microirrigation connectors (Zitterell et al., 2013Zitterell, D. B.; Frizzone, J. A.; Rettore Neto, O. Dimensional analysis approach to estimate local head losses in microirrigation connectors. Irrigation Science, v.32, p.169-179, 2013. https://doi.org/10.1007/s00271-013-0424-y
https://doi.org/10.1007/s00271-013-0424-...
; Vilaça et al., 2017Vilaça, F. N.; Camargo, A. P. de; Frizzone, J. A.; Mateos, L.; Koech, R. Minor losses in start connectors of microirrigation laterals. Irrigation Science, v.35, p.227-240, 2017. https://doi.org/10.1007/s00271-017-0534-z
https://doi.org/10.1007/s00271-017-0534-...
) and to study hydraulics of microtubes (Vekariya et al., 2010Vekariya, P. B.; Subbaiah, R.; Mashru, H. H. Hydraulics of microtube emitters: A dimensional analysis approach. Irrigation Science, v.29, p.341-350, 2010. https://doi.org/10.1007/s00271-010-0240-6
https://doi.org/10.1007/s00271-010-0240-...
).

In this research, the dimensional analysis was based on the Buckingham Pi theorem (Buckingham, 1914Buckingham, E. On physically similar systems: Illustrations of the use of dimensional equations. Physical Review, v.4, p.345-376, 1914. https://doi.org/10.1103/PhysRev.4.345
https://doi.org/10.1103/PhysRev.4.345...
) and the method of repeating variables, following practical procedures defined by Fox et al. (2011Fox, R. W.; McDonald, A. T.; Pritchard, P. J. Introduction to fluid mechanics. 8.ed. Hoboken: John Wiley & Sons , 2011. 896p.).

The Venturi injection rate is mainly influenced by its geometric characteristics, motive flow rate, differential and inlet pressures, and water properties. Based on such empirical interpretation of this physical process, a mathematical relationship can be defined (Eq. 1):

q inj = f Δ P , P in , q m , ρ , μ , D in , D inj , D t (1)

where:

qinj - injection flow rate, m³ s-1;

ΔP - differential pressure, Pa;

Pin - inlet pressure, Pa;

qm - motive flow rate, m³ s-1;

ρ - water density, kg m-3;

μ - water dynamic viscosity, Pa s;

Din - internal inlet diameter, m;

Dinj - injection pipe diameter, m; and,

Dt- throat diameter, m.

Lemons (2017Lemons, D. S. A student’s guide to dimensional analysis. New York: Cambridge University Press, 2017. 112p.) pointed out that the number of variables necessary to represent a process should be minimized. Including unnecessary predictors in the model complicates the description of the process and may result in poor predictions, while omitting important effects reduce its predictive power (Chatterjee & Simonoff, 2013Chatterjee, S.; Simonoff, J. Handbook of regression analysis. Hoboken: John Wiley & Sons, 2013. 236p. https://doi.org/10.1002/9781118532843
https://doi.org/10.1002/9781118532843...
). A model should be as simple as possible while still accounting for the important relationships in the data (Vilaça et al., 2017Vilaça, F. N.; Camargo, A. P. de; Frizzone, J. A.; Mateos, L.; Koech, R. Minor losses in start connectors of microirrigation laterals. Irrigation Science, v.35, p.227-240, 2017. https://doi.org/10.1007/s00271-017-0534-z
https://doi.org/10.1007/s00271-017-0534-...
). In addition, to facilitate grouping of Pi-terms (Π), injection and motive rates were expressed as flow velocities. Finally, the mathematical relationship assumed to represent the physical process is shown in Eq. 2.

v inj = f Δ P , P in , v m , ρ , μ , D in , D inj , D t (2)

where:

vinj - injection flow velocity, m s-1; and,

vm - motive flow velocity at inlet section, m s-1.

The physical properties ρ and μ varies as a function of temperature and fluid. In this research those properties were assumed to be constant: ρ = 1000 kg m-3 and μ = 1.01 x 10-3 Pa s. The quantification of the effects of solute concentration and liquid viscosity and temperature in injection rate of Venturi injectors are described in Yuan et al. (2000Yuan, Z.; Choi, C. Y. Waller, P. M.; Colaizzi, P. Effects of liquid temperature and viscosity on venturi injectors. Transactions of the American Society Agricultural Engineers, v.43, p.1441-1447, 2000. https://doi.org/10.13031/2013.3042
https://doi.org/10.13031/2013.3042...
) and Kumar et al. (2012Kumar, M.; Rajput, T. B. S.; Patel, N. Effect of system pressure and solute concentration on fertilizer injection rate of a venturi for fertigation. Journal of Agricultural Engineering, v.49, p.9-13, 2012.).

The difference between the number of variables that describe a process (k) and the number of reference dimensions (r) required to define units to the list of variables results in the number of dimensionless groups (i.e., the Pi-terms). In Eq. 2 it can be observed nine variables (k = 9) associated with three reference dimensions (M, L, T; r = 3). So, the number of Pi-terms is equal to six (Eq. 3).

Π 1 = ϕ Π 2 , Π 3 , Π 4 , Π 5 , Π 6 (3)

where:

Πn - Pi-term, dimensionless; and,

ϕ - functional relationship between the variables.

The method of repeating variables was used to form the dimensionless groups. This method requires the selection of the so-called “repeating variables” following a specific set of rules mentioned by Fox et al. (2011Fox, R. W.; McDonald, A. T.; Pritchard, P. J. Introduction to fluid mechanics. 8.ed. Hoboken: John Wiley & Sons , 2011. 896p.). From Eq. 2, the injection flow velocity (vinj) was chosen as the dependent variable and the repeating variables were vm, Din and ρ. Thereafter, the repeating variables were systematically combined with the remainder variables in order to define the Pi-terms (Eq. 4).

v inj v m = ϕ Δ P ρ v m 2 , P in ρ v m 2 , ρ v m D in μ , D inj D in , D t D in (4)

The Eq. 4 represents a fundamental relationship that describes the physical process of injection flow velocity in Venturi injectors. Each dimensionless term shown in Eq. 4 corresponds to a Pi-term indicated by Eq. 3.

A multivariate power-law model (Eq. 5) was fitted to the data using the least square method. Such a model has been employed successfully for modeling processes related to fluid flow in hydraulics and irrigation applications (Vekariya et al., 2010Vekariya, P. B.; Subbaiah, R.; Mashru, H. H. Hydraulics of microtube emitters: A dimensional analysis approach. Irrigation Science, v.29, p.341-350, 2010. https://doi.org/10.1007/s00271-010-0240-6
https://doi.org/10.1007/s00271-010-0240-...
; Zitterel et al., 2013; Perboni et al., 2015Perboni, A.; Frizzone J. A.; Camargo, A. P. de; Pinto, M. F. Modelling head loss along emitting pipes using dimensional analysis. Engenharia Agrícola, v.35, p.442-457, 2015. https://doi.org/10.1590/1809-4430-Eng.Agric.v35n3p442-457/2015
https://doi.org/10.1590/1809-4430-Eng.Ag...
; Vilaça et al., 2017Vilaça, F. N.; Camargo, A. P. de; Frizzone, J. A.; Mateos, L.; Koech, R. Minor losses in start connectors of microirrigation laterals. Irrigation Science, v.35, p.227-240, 2017. https://doi.org/10.1007/s00271-017-0534-z
https://doi.org/10.1007/s00271-017-0534-...
).

Π 1 = β 1 Π 2 β 2 Π 3 β 3 Π n β n (5)

where:

βn - empirical coefficient, dimensionless.

The dataset obtained experimentally was randomly divided into two subsets: the calibrating and the testing dataset. The calibrating dataset accounted for 70% of the whole experimental data and was used to fit the models. The testing dataset consisted of the remaining 30% data and was used to assess the accuracy/performance of the model.

Additionally, for each model of Venturi, the motive flow was estimated as a function of ΔP and Pin (Eq. 6). The model was fitted by the software Sigmaplot® and its coefficients were adjusted using the least square method.

q m = α 1 + α 2 P in + α 3 Δ P + α 4 P in 2 + α 5 Δ P 2 (6)

where:

αn - empirical coefficient, dimensionless.

Models were assessed by the root mean square error (RMSE) and by graphical error analysis. The RMSE is a common index to measure the accuracy of models that quantify differences between observed and estimated values (Duran-Ros et al., 2010Duran-Ros, M.; Arbat, G.; Barragán, J.; Ramírez de Cartagena, F.; Puig-Bargués, J. Assessment of head loss equations developed with dimensional analysis for micro irrigation filters using effluents. Biosystems Engineering, v.106, p.521-526, 2010. https://doi.org/10.1016/j.biosystemseng.2010.06.001
https://doi.org/10.1016/j.biosystemseng....
; Provenzano et al., 2016Provenzano, G. P.; Alagna, V.; Autovino, D.; Juarez, J. M.; Rallo, G. Analysis of geometrical relationships and friction losses in small-diameter lay-flat polyethylene pipes. Journal of Irrigation and Drainage Engineering , v.142, p.1-9, 2016. https://doi.org/10.1061/(ASCE)IR.1943-4774.0000958
https://doi.org/10.1061/(ASCE)IR.1943-47...
). The graphical error analysis is also useful to quantify prediction errors while evaluating the accuracy of models because it provides prediction errors associated with their frequency of occurrence (Vilaça et al., 2017Vilaça, F. N.; Camargo, A. P. de; Frizzone, J. A.; Mateos, L.; Koech, R. Minor losses in start connectors of microirrigation laterals. Irrigation Science, v.35, p.227-240, 2017. https://doi.org/10.1007/s00271-017-0534-z
https://doi.org/10.1007/s00271-017-0534-...
).

Results and Discussion

For each model of Venturi, 25 testing conditions were evaluated. Therefore, the gathered experimental data accounts for 100 records obtained from functional tests of four models of Venturi injectors. Figure 3 shows the experimental results.

Figure 3
Data obtained from the functional tests of four models of Venturi injectors evaluated under 5 conditions of inlet pressure: (A, C, E and G) injection flow as a function of motive flow; (B, D, F and H) injection flow as a function of ΔP/Pin

The injection flow rate was more sensitive to changes in the differential pressure (ΔP) when the Venturi was operated under low values of inlet pressure (Pin - i.e., 100 and 150 kPa). Generally, for Pin higher than 200 kPa and for ΔP/Pin higher than 60%, qinj was relatively constant. According to Lamm et al. (2006Lamm, F. R.; Ayars, J. E.; Nakayama, F. S. Microirrigation for crop production. Amsterdam: Elsevier, 2006. 642p.), the injection rate is relatively insensitive to irrigation pipeline pressure changes when Venturi injectors are operated under inlet pressures higher than about 200 kPa. Moreover, the authors recommend maintaining high values of ΔP/Pin when a constant chemical injection is required.

For inlet pressures ranging from 100 to 300 kPa, qinj varied from 0 to 0.095 m³ h-1, 0 to 0.483 m³ h-1, 0 to 0.652 m³ h-1 and 0.006 to 0.422 m³ h-1 for Venturi Models 1, 2, 3 and 4, respectively (Figure 3). As shown in Table 1, all geometrical characteristics increase for Venturi injector Models 1 to 3. Model 3 presented the highest values of qinj (Figures 3E and F), which is an expected result since that model has the largest dimensions. The injection rate of Venturi Model 4 (Figures 3G and H) represents 64.7 ± 5.2% of Model 3. Models 3 and 4 have the same nominal diameter (DN), but different internal dimensions. The inlet diameter (Din) and the throat diameter (Dt) of Model 4 are 14.5% and 1.9% smaller than Model 3, respectively. Manzano et al. (2018Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Characterization and selection method of Venturi injectors for pressurized irrigation. Revista Ciência Agronômica , v.49, p.201-210, 2018. https://doi.org/10.5935/1806-6690.20180023
https://doi.org/10.5935/1806-6690.201800...
) pointed out that some manufacturers give insufficient information in their catalogs, mainly in relation to the dimensions and diameters. The parameter normally used by manufacturers to differentiate Venturi tubes is a nominal diameter, without specifying the respective internal diameters, the internal throat diameter or the injection diameter (Dinj).

Comparing the same manufacturer Venturi tubes, Model 1 (Figures 3A and B) represents on average 19.6 ± 1.3% and 14.57 ± 1.0% of the injection rate capacity of Models 2 and 3, respectively. Venturi tube Model 2 (Figures 3C and D) represents 74.1 ± 3.3% of the Model 3 injection rate capacity as well. In this case, the motive and injection flow rates increased with the internal and injection diameters. Manzano et al. (2015Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Diseño y alternativas em la instalación de inyectores Venturi en riego localizado. Revista Ciência Agronômica, v.46, p.287-298, 2015.), while evaluating performance and alternative setup of four Venturi tubes, reported a similar trend.

The empirical coefficients of Eq. 5 were fitted using the calibration dataset. A full model was evaluated and shown in Eq. 7.

v inj v m = 0.0046 Δ P ρ v m 2 0.2841 P in ρ v m 2 1.4899 ρ v m D in μ 0.0591 D inj D in 2.7966 D t D in 10.9511 (7)

Eq. 7 is valid for: 12.5258 ≤ ΔP/ρv2m ≤ 59.3003; 22.1238 ≤ Pin/ρv2m ≤ 64.3450; 15563 ≤ ρvmDin/μ ≤ 87544; 0.3827 ≤ Dinj/Din ≤ 0.6581; 0.3109 ≤ Dt/Din ≤ 0.3889. The units of the terms in the equations are: ΔP and Pin = Pa; qinj and qm = m³ s-1; Din, Dinj, and Dt = m; ρ = kg m-3; μ = Pa s.

Expressing the flow velocities as flow rates and rearranging the terms of Eq. 7, the proposed model to estimate the Venturi injection rate is presented in Eq. 8.

q inj = 0.0082 ρ 1.1468 μ 0.0591 q m 3.3527 Δ P 0.2841 P in 1.4899 D in 1.3901 D inj 4.7966 D t 10.9511 (8)

Figure 4A shows the fitting and validation scattering data, and RMSE values as well, of the predicted qinj by the model. The RMSE values were very similar which indicate the accuracy of the model. For this model, relative errors lower than 10% were observed in 79.5% of the predictions, while 95% of predictions presented relative errors of up to 29.1% (Figure 4C).

Figure 4
Comparison between the proposed model on predicting the injection (qinj) and motive (qm) flow rates: (A) observed versus estimated values of qinj; (B) observed versus estimated values of qm; and (C) graphical error analysis presenting relative errors (δ) versus frequency prediction errors

Due to the high dependence among qinj, qm, ΔP, and Pin, and the observed high sensitivity of qinj to variations in qm during the data analysis (Figure 3), it was not possible to fit a general model to predict qm. According to Frizzone et al. (2012Frizzone, J. A.; Freitas, P. S. L. de; Rezende, R.; Faria, M. A. de. Microirrigação: Gotejamento e microaspersão. Maringá: EDUEM, 2012. 356p.) and Manzano et al. (2018Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Characterization and selection method of Venturi injectors for pressurized irrigation. Revista Ciência Agronômica , v.49, p.201-210, 2018. https://doi.org/10.5935/1806-6690.20180023
https://doi.org/10.5935/1806-6690.201800...
), changes in flow rate or operating pressure may affect the other hydraulic parameters. In this way, a multivariate model (Eq. 6) was fitted for each evaluated Venturi injector as a function of ΔP and Pin, and their coefficients are shown in Table 2.

Table 2
Coefficients adjusted for each model of Venturi to estimate qm as a function of ΔP and Pin (Eq. 6)

The RMSE values of the predicted qm for each model are presented in Figure 4B and the frequency of occurrence associated with values of relative errors (Figure 4C) indicates the similar accuracy of both models.

Manzano et al. (2015Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Diseño y alternativas em la instalación de inyectores Venturi en riego localizado. Revista Ciência Agronômica, v.46, p.287-298, 2015.) evaluated four different Venturi models installed in-line. By regressions, the authors obtained two models for predicting the injection rate. However, one of the models only varies by the differential pressure, while the other varies with both the differential pressure and inlet pressure. When these models were obtained, the authors did not take into account the variation of dimensions of the Venturi tubes, the geometric characteristics of the injection pipes, and water properties.

In practice, Venturi injectors are commonly installed in a by-pass setup ISO 15873 (ISO, 2002ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.), which renders controlling the motive flow rate more difficult due to the withdrawing of water from the main irrigation pipeline (Manzano et al., 2018Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Characterization and selection method of Venturi injectors for pressurized irrigation. Revista Ciência Agronômica , v.49, p.201-210, 2018. https://doi.org/10.5935/1806-6690.20180023
https://doi.org/10.5935/1806-6690.201800...
). However, a single manometer would make it possible to control the differential pressure by opening or closing a valve installed across the irrigation pipeline, which generates negative pressure inside the injector and, consequently, the suction of the solution.

To demonstrate the use of the proposed equations, it is necessary to choose the size and to estimate the differential pressure of a Venturi injector that is able to operate in the following application. A microirrigation subunit with an area of 5 ha; the pressure at the inlet of the subunit is 300 kPa; irrigation duration of 5 h; and a need to apply 30 kg ha-1 of nitrogen from a liquid fertilizer with a nutrient concentration of 0.42 kg L-1. The fertigation duration normally corresponds to 80% of the irrigation duration. It is a practical recommendation commonly used for the correct application of nutrients to the crop (Frizzone et al., 2012Frizzone, J. A.; Freitas, P. S. L. de; Rezende, R.; Faria, M. A. de. Microirrigação: Gotejamento e microaspersão. Maringá: EDUEM, 2012. 356p.).

As calculation steps, the amount of liquid fertilizer required is 357.14 L, and the fertigation duration is 4 h; hence the required injection flow is 89.3 L h-1. From this value, the size of the Venturi injector can be defined. Venturi Model 1 fulfills the application since its injection flow ranges from 0 to 0.095 m3 h-1 (Figure 3).

The geometric characteristics of the Venturi Model 1 are known (Table 1). It is assumed that the water dynamic viscosity is 1.003 x 10-3 Pa s and its density is 1000 kg m-3. Therefore, the unknown variables in Eq. 8 are ΔP and qm. From Eq. 6 adjusted for Venturi Model 1 (coefficients given in Table 2), the unknown variables also are ΔP and qm. The differential pressure required to obtain the injection flow rate of 0.0893 m3 h-1 is calculated by iterations. Finally, the differential pressure should be adjusted to 212.8 kPa in order to obtain the target injection flow rate of 0.0893 m3 h-1. Under this operating condition, the motive flow is 1.024 m³ h-1, the ratio ΔP/Pin is 70.9%, and the pressure at the outlet of the Venturi is 87.2 kPa. This operating condition will enable a stable and accurate application of the desired amount of nitrogen. Also, the pressure of 87.2 kPa is enough to properly operate most emitters employed in microirrigation systems.

Conclusions

  1. For all evaluated models, the injection flow rate was positive for ΔP/Pin ratios higher than 20%, regardless of the inlet pressure.

  2. The injection flow rate was more sensitive to changes in the motive flow rate and the differential pressure when the Venturi injector was operated under low values of inlet pressure (i.e., 100 and 150 kPa).

  3. Generally, under inlet pressures higher than 200 kPa and for ΔP/Pin ratios higher than 60%, the injection flow rate was relatively constant.

  4. A general model for predicting injection flow rate was proposed and validated, which depends on the fluid properties, hydraulic operating conditions and geometrical characteristics of the Venturi injector. Another model for estimating motive flow rate as a function of inlet pressure and differential pressure was adjusted and validated for each size of the Venturi injector.

Literature Cited

  • Buckingham, E. On physically similar systems: Illustrations of the use of dimensional equations. Physical Review, v.4, p.345-376, 1914. https://doi.org/10.1103/PhysRev.4.345
    » https://doi.org/10.1103/PhysRev.4.345
  • Chatterjee, S.; Simonoff, J. Handbook of regression analysis. Hoboken: John Wiley & Sons, 2013. 236p. https://doi.org/10.1002/9781118532843
    » https://doi.org/10.1002/9781118532843
  • Duran-Ros, M.; Arbat, G.; Barragán, J.; Ramírez de Cartagena, F.; Puig-Bargués, J. Assessment of head loss equations developed with dimensional analysis for micro irrigation filters using effluents. Biosystems Engineering, v.106, p.521-526, 2010. https://doi.org/10.1016/j.biosystemseng.2010.06.001
    » https://doi.org/10.1016/j.biosystemseng.2010.06.001
  • Feitosa Filho, J. C.; Botrel, T. A.; Pinto, J. M. Desempenho de injetores tipo venturi em condições de descarga livre e pressurizada. Engenharia Agrícola, v.17, p.20-35, 1998.
  • Feitosa Filho, J. C.; Pinto, J. M.; Arruda, N. T. de. Dimensionamento, construção e características hidráulicas de um injetor tipo venturi para uso na quimigação. Irriga, v.4, p.43-53, 1999. https://doi.org/10.15809/irriga.1999v4n2p43-53
    » https://doi.org/10.15809/irriga.1999v4n2p43-53
  • Fox, R. W.; McDonald, A. T.; Pritchard, P. J. Introduction to fluid mechanics. 8.ed. Hoboken: John Wiley & Sons , 2011. 896p.
  • Frizzone, J. A.; Freitas, P. S. L. de; Rezende, R.; Faria, M. A. de. Microirrigação: Gotejamento e microaspersão. Maringá: EDUEM, 2012. 356p.
  • Gupta, B.; Nayak, A. K.; Kandar, T. K.; Nair, S. Investigation of air-water two phase flow through a Venturi. Experimental Thermal and Fluid Science, v.70, p.148-154, 2016. https://doi.org/10.1016/j.expthermflusci.2015.07.012
    » https://doi.org/10.1016/j.expthermflusci.2015.07.012
  • Hoffman, G. J.; Evans, R. G.; Jensen, M. E.; Martin, D. L.; Elliott, R. L. Design and operation of farm irrigation systems. 2.ed. St. Joseph: American Society of Agricultural and Biological Engineers, 2007. 863p.
  • ISO - International Organization for Standardization. ISO 15873: Irrigation equipment - Differential pressure Venturi-type liquid additive injectors. Genebra: International Organization for Standardization, 2002. 12p.
  • Kuldeep, R.; Saharan, V. K. Computational study of different venturi and orifice type hydrodynamic cavitating devices. Journal of Hydrodynamics, v.28, p.293-305, 2016. https://doi.org/10.1016/S1001-6058(16)60631-5
    » https://doi.org/10.1016/S1001-6058(16)60631-5
  • Kumar, M.; Rajput, T. B. S.; Patel, N. Effect of system pressure and solute concentration on fertilizer injection rate of a venturi for fertigation. Journal of Agricultural Engineering, v.49, p.9-13, 2012.
  • Lamm, F. R.; Ayars, J. E.; Nakayama, F. S. Microirrigation for crop production. Amsterdam: Elsevier, 2006. 642p.
  • Lemons, D. S. A student’s guide to dimensional analysis. New York: Cambridge University Press, 2017. 112p.
  • Lima Neto, I. E.; Porto, R. M. Performance of low-cost ejectors. Journal of Irrigation and Drainage Engineering, v.130, p.122-128, 2004. https://doi.org/10.1061/(ASCE)0733-9437(2004)130:2(122)
    » https://doi.org/10.1061/(ASCE)0733-9437(2004)130:2(122)
  • Liu, K.; Lu, H.; Guo, X.; Sun, X.; Tao, S.; Gong, X. Experimental study on flow characteristics and pressure drop of gas-coal mixture through Venturi. Power Technology, v.268, p.401-411, 2014. https://doi.org/10.1016/j.powtec.2014.07.019
    » https://doi.org/10.1016/j.powtec.2014.07.019
  • Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Diseño y alternativas em la instalación de inyectores Venturi en riego localizado. Revista Ciência Agronômica, v.46, p.287-298, 2015.
  • Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V.; Vasconcelos, D. V. Geometry and head loss in Venturi injectors though computational fluid dynamics. Engenharia Agrícola, v.36, p.482-491, 2016. https://doi.org/10.1590/1809-4430-Eng.Agric.v36n3p482-491/2016
    » https://doi.org/10.1590/1809-4430-Eng.Agric.v36n3p482-491/2016
  • Manzano, J.; Palau, C. V.; Azevedo, B. M. de; Bomfim, G. V. do; Vasconcelos, D. V. Characterization and selection method of Venturi injectors for pressurized irrigation. Revista Ciência Agronômica , v.49, p.201-210, 2018. https://doi.org/10.5935/1806-6690.20180023
    » https://doi.org/10.5935/1806-6690.20180023
  • Perboni, A.; Frizzone J. A.; Camargo, A. P. de; Pinto, M. F. Modelling head loss along emitting pipes using dimensional analysis. Engenharia Agrícola, v.35, p.442-457, 2015. https://doi.org/10.1590/1809-4430-Eng.Agric.v35n3p442-457/2015
    » https://doi.org/10.1590/1809-4430-Eng.Agric.v35n3p442-457/2015
  • Provenzano, G. P.; Alagna, V.; Autovino, D.; Juarez, J. M.; Rallo, G. Analysis of geometrical relationships and friction losses in small-diameter lay-flat polyethylene pipes. Journal of Irrigation and Drainage Engineering , v.142, p.1-9, 2016. https://doi.org/10.1061/(ASCE)IR.1943-4774.0000958
    » https://doi.org/10.1061/(ASCE)IR.1943-4774.0000958
  • Santos, L. da C.; Zocoler, J. L.; Justi, A. L.; Silva, A. O.; Correia, J. de S. Estudo comparativo da taxa de injeção em injetor do tipo Venturi com e sem válvula de retenção. Irriga, v.edição especial, p.145-154, 2012. https://doi.org/10.15809/irriga.2012v1n01p145
    » https://doi.org/10.15809/irriga.2012v1n01p145
  • Silva, P. F. da; Lima, C. J. G. de S.; Barros, A. C.; Silva, E. M. da; Duarte, S. N. Sais fertilizantes e manejo da fertirrigação na produção de tomateiro cultivado em ambiente protegido. Revista Brasileira de Engenharia Agrícola e Ambiental, v.17, p.1173-1180, 2013. https://doi.org/10.1590/S1415-43662013001100007
    » https://doi.org/10.1590/S1415-43662013001100007
  • Sinha, I.; Buttar, G. S.; Brar, A. S. Drip irrigation and fertigation improve economics, water and energy productivity of spring sunflower (Helianthus annuus L.) in Indian Punjab. Agricultural Water Management, v.185, p.58-64, 2017. https://doi.org/10.1016/j.agwat.2017.02.008
    » https://doi.org/10.1016/j.agwat.2017.02.008
  • Tian, D.; Zhang, Y.; Zhou, Y.; Mu, Y.; Liu, J.; Zhang, C.; Liu, P. Effect of nitrification inhibitors on mitigating N2O and NO emissions from an agricultural field under drip fertigation in the North China Plain. Science of the Total Environment, v.598, p.87-96, 2017. https://doi.org/10.1016/j.scitotenv.2017.03.220
    » https://doi.org/10.1016/j.scitotenv.2017.03.220
  • Vekariya, P. B.; Subbaiah, R.; Mashru, H. H. Hydraulics of microtube emitters: A dimensional analysis approach. Irrigation Science, v.29, p.341-350, 2010. https://doi.org/10.1007/s00271-010-0240-6
    » https://doi.org/10.1007/s00271-010-0240-6
  • Vilaça, F. N.; Camargo, A. P. de; Frizzone, J. A.; Mateos, L.; Koech, R. Minor losses in start connectors of microirrigation laterals. Irrigation Science, v.35, p.227-240, 2017. https://doi.org/10.1007/s00271-017-0534-z
    » https://doi.org/10.1007/s00271-017-0534-z
  • Wu, W.; Chen, W. E. I.; Liu, H; Yin, S.; Niu, Y. A new model for head loss assessment of screen filters developed with dimensional analysis in drip irrigation systems. Irrigation and Drainage, v.63, p.523-531, 2014. https://doi.org/10.1002/ird.1846
    » https://doi.org/10.1002/ird.1846
  • Xu, Y.; Zhang, Q.; Zhang, T.; Ba, X. An overreading model for nonstandard Venturi meters based on H correction factor. Measurement, v.61, p.100-106, 2015. https://doi.org/10.1016/j.measurement.2014.10.028
    » https://doi.org/10.1016/j.measurement.2014.10.028
  • Yuan, Z.; Choi, C. Y. Waller, P. M.; Colaizzi, P. Effects of liquid temperature and viscosity on venturi injectors. Transactions of the American Society Agricultural Engineers, v.43, p.1441-1447, 2000. https://doi.org/10.13031/2013.3042
    » https://doi.org/10.13031/2013.3042
  • Zitterell, D. B.; Frizzone, J. A.; Rettore Neto, O. Dimensional analysis approach to estimate local head losses in microirrigation connectors. Irrigation Science, v.32, p.169-179, 2013. https://doi.org/10.1007/s00271-013-0424-y
    » https://doi.org/10.1007/s00271-013-0424-y

Publication Dates

  • Publication in this collection
    10 July 2019
  • Date of issue
    July 2019

History

  • Received
    07 May 2018
  • Accepted
    13 May 2019
  • Published
    04 June 2019
Unidade Acadêmica de Engenharia Agrícola Unidade Acadêmica de Engenharia Agrícola, UFCG, Av. Aprígio Veloso 882, Bodocongó, Bloco CM, 1º andar, CEP 58429-140, Campina Grande, PB, Brasil, Tel. +55 83 2101 1056 - Campina Grande - PB - Brazil
E-mail: revistagriambi@gmail.com