Acessibilidade / Reportar erro

NUMERICAL SIMULATION AND EXPERIMENTAL ANALYSIS OF SPRAY FIELD OF IMPINGEMENT NOZZLE

ABSTRACT

The shape of a spray droplet from a flat-fan nozzle in the plant protection sector is similar to a ‘U shape’, the middle particle size being small but the diameter of both sides being too large. So, flat-fan nozzles have the defects of wide droplet size and poor uniformity of droplet size. This paper simulates and tests the impingement nozzle, based on the coupling effect of the jet and impinging stream, so as to improve the droplet size uniformity and droplet spectrum. After jet impingement, it was found that the droplet size of the spray fan surface of the impingement nozzle showed the characteristics of large, uniform and smaller sides. This was beneficial in solving the uneven distribution of droplet size distribution in the field operation of the sprayer. The droplet distribution of the impingement nozzle is more concentrated, forming a narrow droplet spectrum. The experimental results are in good agreement with the results of numerical simulations and theoretical analysis. The impingement nozzle solves the shortcomings of wider fan droplets and smaller particle size uniformity, and provides feasible technical support for developing intelligent precision spray equipment.

impingement nozzle; numerical simulation; spray field; particle size uniformity; droplet spectrum

INTRODUCTION

As an important tool for protecting agricultural, forestry, fruit and livestock production, plant protection machinery plays an important role in protecting stable and high crop yields ( Dong et al., 2019Dong F, Zhou H, Shi M, Ru Y, Yi K (2019) Design and Atomisation Test of Impinging nozzle. Scientia Silvae Sinicae 55(1): 83-90. ). As one of the important components of the spraying system of plant protection machinery, the performance of the nozzle is related to the reliability and economy of the whole plant protection machinery and the operation of the whole system; it is also a key element in the improvement of spraying technology ( Wang, 2019Wang G (2019) Spraying Quality Evaluation on Plant Protection UAV. Guangzhou: South China Agricultural University. 2019.6. ; Dong & Zhou, 2018Dong F, Zhou H (2018) Development of foreign plant protection nozzles. Acta Agriculturae Universitatis Jiangxiensis 40(4):866-874. ). Nozzle jet spray atomisation is a multi-phase flow stochastic process with gas-liquid phase interaction, and there is no universally applicable relevant theory for the atomisation mechanism. Researchers at home and abroad have conducted a lot of experimental and theoretical research in this area (He, 2019; Song et al., 2017Song S, Chen J, Hong T, Xue X, Xia H, Song Y (2017) Variation of droplet diameter in wind field for long-range air-assisted sprayer. Transactions of the Chinese Society of Agricultural Engineering 33(6):59-66. ; Salcedo et al., 2020)Salcedo R, Ramon L, Jordi C, Javier C, Michael G, Montserrat O, Paula G, Emi L (2020) Evaluation of leaf deposit quality between electrostatic and conventional multi-row sprayers in a trellised vineyard. Crop Protection 127:104964. . Lü et al. (2007)Lü X, He X, Song J, Zeng A, Andreas H (2007) Analysis of spray process produced by agriculture flat-fan nozzles. Transactions of the Chinese Society of Agricultural Engineering 23(9):95-100. studied the spatial distribution of droplet size and motion of a standard fan nozzle and found that the droplet particle size distribution showed a concave ellipsoidal distribution, with a small middle edge

and a large edge, when gradually moving away from the nozzle terminal outlet section. The droplet velocity in the axial direction was a hill-shaped distribution with a large middle edge. Zheng et al. (2015)Zheng G, Nie W, He B, Xue C (2015) Study on the influence of impact angle on the atomisation characteristics of impinging nozzle. Journal of Propulsion Technology 36(4):608-613. used the CLSVOF method to study the impact angle of the collision nozzle on the atomisation. It was found that the fragmentation length of the liquid film decreases as the impact angle increases, while the surface wave amplitude and fragmentation of the liquid film also increases. The velocity difference at the impact position plays a key role in the fragmentation and surface wave amplitude of the liquid film. Fei et al. (2015)Fei J, Sun F, Yang W, Fu Y, Wang Y (2015) Experimental analysis on movement and size distribution of atomised droplets from impinging liquid jet. Journal of Rocket Propulsion 41(1):10-14. studied the droplet motion process and particle size distribution characteristics of jet impingement nozzle atomisation under atmospheric conditions. The analysis of the process of impingement jet atomisation and liquid film fragmentation into droplets found that the velocity of atomised droplets in the radial velocity component (relative to the origin as the centre of symmetry) away from the origin increases. The radial velocity distribution and the initial atomisation stage droplet size distribution was the opposite.

The hedge nozzle is a fan-impact nozzle based on the coupling effect of impinging the flow and the jet, which has the characteristics of a narrow droplet spectrum and better particle size uniformity, relative to the fan nozzle ( Dong et al., 2018Dong F, Zhou H, Ru Y Shi M, Chen Q, Yi K (2018) Experiment on droplet distribution characteristics in spray field of impinging nozzle. Transactions of the Chinese Society for Agricultural Machinery 49(12):116-121+128. ). The analysis of the fluid flow field in the impinging region outside the nozzle, using Fluent software, is of great significance to the spray mechanism and nozzle structure optimisation.

Hedge nozzle working principle and geometric modelling

Working principle

Many researchers ( Zhang et al., 2022Zhang J, Zhang Z, Feng Y, Dong X (2022) Advances in multiscale analysis of impinging streams flow field and mixing characteristics. Journal of Engineering Thermophysics 43(11):2944-2956. ; Cai et al., 2022Cai Y, Duan Z, Niu J, Li D, Zhang J, Qiao G (2022) Numerical simulation of flow field characteristics in impinging stream continuous crystallizer. Chemical Engineering (China) 50(11):43-49. ) found that the impact of the two jets can make the droplet particle size distribution narrower, homogenising the droplet group particle size; the more intense the impact, the more obvious the particle size homogenisation effect. As shown in Figure 1 , the hedge nozzle has two jets from the terminal outlet at a certain angle, which hit each other at the impact surface, and the impact makes the droplet particle size distribution narrower and homogenises the droplet group. This improves the current fan spray as the droplet spectrum is wider but the particle size uniformity is poor.

FIGURE 1
Schematic diagram of impinging jet.

Geometric modelling

Figure 2(a) shows part of the experimental test set-up: the outlet orifice diameter was 1 mm, the distance between the two outlets was 7 mm, the notch angle of the hedge nozzle was 30° (two jets impacting before the spray surface) and the vertical impact surface angle was 15°. The internal fluid region of the model nozzle is shown in Figure 2(b) and the black area in Figure 2 (c) shows the fluid domain inside the hedge nozzle.

FIGURE 2
Geometric model.

According to the model of the hedge nozzle geometry, the model was meshed using the ANSYS mesh pre-processing software ‘ICEM CFD’ ( Gu et al., 2015Gu Z, Zhao W, Zhang Y, Wei H (2015) A Research on the effects of mesh quality evaluation criteria on the results of vehicle external flow field simulation[J]. Automotive Engineering 37(11):1265-1269. ), as shown in Figure 3 . In order to obtain more accurate calculation results, the spray field mesh was divided into three parts in the back section of the spray nozzle terminal orifice outlet. According to the actual situation, and taking into account the dense nozzle jet impact region, jet impact

FIGURE 3
Meshing of the model.

and droplet fragmentation occurs in this region. The second part of the mesh for local grid encryption processing (mm), the first part of the grid second dense (mm), and the third part of the grid is the coarsest (mm). In order to facilitate the calculation and to obtain the results, a hexahedral-structured grid topology was used, with a mesh count of 695557.

Model and parameter settings

Mathematical model

VOF model

A traditional fan nozzle only has one jet, and the ‘DPM’ discrete phase processing method is usually used for simulation; it is believed that the spray droplets will be formed directly at the outlet, allowing visual simulation of the spray nozzle atomisation process ( Wang et al., 2022Wang X, Liu J, Zhang Q (2022) Water–Pesticide Integrated Micro-Sprinkler Design and Influence of Key Structural Parameters on Performance. Agriculture 12(10):1532. ). The hedge nozzle originated from the jet impact, and the outlet cannot be used in this spray mode model. The water flow just comes out of the nozzle’s terminal orifice outlet and can be considered to be jet-like. This then collides with the impact surface or impact centre near the outlet, to form many fine droplets. The use of VOF two-phase flow model processing can better fit the actual situation. However, the VOF model can only show the velocity distribution and shape of the spray field, it cannot obtain the spray nozzle atomisation process condition. The model uses volume fraction to characterise the volume percentage of the second phase fluid in the computational grid cell. In the calculation, when F=0, it is considered that the grid is all gas; when F=1, it is considered that the grid is all liquid and, when 0<F<1, it is considered that the grid is a gas-liquid mixture.

F t + v F = 0 (1)

In the equation, time, Speed vector, Fluid Density, Fluid viscosity. In the gas-liquid mixing zone, fluid density and viscosity can be expressed as ρ=(1F)ρg+Fρl and μ=(1F)μg+Fμl , respectively.

Gas-liquid phase control equations

In the study of the flow field process outside the hedge nozzle, the gas and liquid are considered to be incompressible fluids, and the controlling equations have a mass and momentum conservation equation. The mass conservation equation is:

ρ ¯ t + ( ρ ¯ v i ) x i = 0 (2)

The conservation of momentum equation is given by:

( ρ ¯ u ¯ i ) t + ( ρ ¯ u ¯ j u ¯ i ) x j = p ¯ x i + μ ¯ 2 u ¯ i + μ ¯ 3 2 u ¯ j x j x i + S ¯ (3)

In the equation, ρ̄ is density, μ̄ is power viscosity, ūi is the speed vector, and is a generalised source term of the momentum equation.

Turbulence model

A model which achieved(Realizable) k-ε was chosen for the study.

The turbulent kinetic energy( К )equation is:

( ρ k ) t + ( ρ k U ) x i = x j [ ( μ + μ t σ k ) k x j ] + G k + G b ρ ε Y M (4)

The turbulence dissipation rate(ε)equation is:

( ρ ε ) t + ( ρ v i ) x i = x j [ ( μ + μ 4 σ ε ) ε x j ] + ρ C 1 E ε C 2 ε ρ ε 2 k + v ε + C 1 ε ε k C 3 ε C b (5)

In the equation, C1=max(0.43,η/(η+5))η=(2EijEij)1/2k/ε , and the turbulent viscosity coefficient μt=ρCμk2/ε . Gk is the generation term of turbulent kinetic energy К due to the mean velocity gradient, Gb is the generation term of turbulent kinetic energy К due to buoyancy, YM is the contribution of the pulsation expansion to the total dissipation rate in compressible turbulent flows, ( Gμ,G,G ) and G indicates the empirical constant. σ k is the Prandtl number corresponding to turbulent kinetic energy, and σ ε is the Prandtl number corresponding to the dissipation rate. In FLUENT, Gε1= 1.44, G= 1.9, G= 0.09, σ k = 1.0 , and σ ε = 1.2 .

Boundary condition setting

The nozzle’s external flow field was modelled using the following basic assumptions: nozzle external flow is a two-stream transient flow problem, the fluid and air are an incompressible Newtonian fluid, the physical properties are constant and there are no chemical reactions, and the effect of gravity on atomisation is ignored without considering convective heat transfer.

The experiment investigated the effect of three different spray pressures on the spray condition of the nozzle, setting the nozzle inlet boundary in the boundary definition as the pressure-inlet during the simulation; the pressure-outlet was connected to the external atmospheric pressure. Therefore, it is necessary to know the external atmospheric pressure. The velocity at the inlet of the external flow field (i.e. the pressure-outlet of the internal flow field) was obtained from the calculation of the internal flow field and the data are the same as in the internal flow field simulation. The nozzle inlet in the boundary condition was the nozzle inlet during the internal flow field simulation; the inlet pressure parameters were set to 0.3, 0.5 and 0.7 MPa, respectively. Air pressure was the relative static pressure at 0.0 MPa (i.e. atmospheric pressure). Air density was set at 1.225 kg/m3 and viscosity was 1.789 kg/(ms). Using no-slip boundary conditions, the jet medium was water with a density of 998 kg/m3 and a viscosity of 0.001 kg/(ms); a time step of 1e-5 was used in the calculations.

Simulation results and analysis

Cloud analysis of impact models

The impact model clouds of the hedge nozzles at spraying pressures of 0.3, 0.5 and 0.7 MPa are shown in Figure 4 . It is found that the droplet particle size, after the impact of the two jets, becomes smaller and more uniform as the spray application pressure increases, which is consistent with the conclusion that the stronger the impact is, the more favourable the jet atomisation. As reflected in Figure 4 , it can be seen that the splitting and uniformity of the liquid stream cluster has more obvious changes with increasing pressure; Figure 4 (c) shows the best uniformity of the hedge nozzle atomisation particle size at 0.7 MPa.

FIGURE 4
Impact model cloud picture.

Distribution analysis of droplet velocity

Lu Xiaolan, Liu Qiusheng and many other scholars found that droplet velocity plays a significant role in the spray atomisation process. The commonly used fan nozzle droplet particle size in the spray sector is not the same, in the same height horizontal direction presents the middle of the droplet diameter is small to gradually increase on both sides, while the measured droplet velocity in the axial direction shows a trend of gradually decreasing in the middle and at both sides. This indicates that there is a certain relationship between droplet velocity and droplet diameter and analysis and speculation suggests that

droplets with similar velocities in the axial direction, at the same height and horizontal direction, will tend to have the same droplet size.

For spraying pressures of 0.3, 0.5 and 0.7 MPa, the test points were defined using FLUENT software. The coordinates of the points were entered and calculated, and the resulting data files were processed to derive the velocities in the x, y and z directions for each test point, i.e. the velocity components of the droplets in each direction, as shown in Figure 5 .

FIGURE 5
Droplet velocity output.

The specific test point location of the hedging nozzle is shown in Figure 6 . Figure 6 (c) shows the geometric model of the hedging nozzle body at the bottom end of the centre of the circle: the vertical direction represents the z-axis, the horizontal direction of the spray sector is the x-axis, perpendicular to the spray sector direction for the y-axis, and x, y and z are the three axes of direction in line with the right-hand spiral relationship. By testing the velocity components of the droplet velocity in the x, y and z directions in the XOZ plane, and at vertical heights of 18, 23 and 28 mm from the origin, the test point locations are shown in Figures 6(a) and 6(b). The measured velocity components of the droplet velocity in the x, y and z directions are shown in Tables 1, 2 and 3, for each test point symmetrical to the left and right, i.e. the -2.0, -1.5, -1.0, -0.5, 0.0, 0.5, 1.0, 1.5 and 2.0 mm positions measured from the centreline position (i.e. z-axis) of the spray sector.

FIGURE 6
Test location diagram.

From Table 1 , it can be seen that, in the vertical direction (i.e. the height or z-direction), at a spraying pressure of 0.3 MPa, the droplet velocity component at 18 mm is smaller than that at 23 mm (at the point -0.5, 0.0, 0.5). Because the droplets are closer to the impact centre, the droplets may be more dense and not yet fully atomised. The velocity components of droplets at different pressures in the z-direction, i.e. vertical height, are close to those in the horizontal direction (e.g. at 0.5 Mpa pressure). At a horizontal height of 23 mm, the velocity components of droplets are between 70.33 and 73.70 m/s. As the spraying pressure increases, the droplet velocity increases accordingly, and at the outlet away from the nozzle, the velocity decreases. For example, at 0 mm on the z-axis, and with a spray pressure of 0.5 Mpa (at 18, 23 and 28 mm), the speeds are 74.821, 73.701, and 72.960 m/s, respectively. At a spray pressure of 0.7 Mpa, at 18, 23 and 28 mm, the speeds are 94.108, 93.452, and 90.698 m/s, respectively. As the droplet moves away from the nozzle’s terminal outlet, the droplet velocity will be reduced, mainly because of the role of air disturbance/resistance. As the hedge nozzle after impact droplet velocity in the z-axis direction of the velocity component is similar, then the same height in the spray field in the horizontal direction of the droplet diameter should tend to the same trend (i.e. uniform particle size).

TABLE 1
Fog droplet velocity component in vertical direction.

As can be seen from Table 2 , in the left and right direction (i.e. the x-direction), and for the horizontal direction of the same height at different pressures, the droplet velocity in the x-axis component shows a trend of gradually increasing from the middle towards both sides. With an increase in pressure, the droplet velocity component also increases but in a position closer to the central axis (the z-axis), especially near the z-axis line. The droplet velocity component changes irregularly and the value is small; the gradual increase in the droplet velocity component near both sides of the edge is conducive to the formation of the spray fan surface.

TABLE 2
Fog droplet velocity component in horizontal direction.

From Table 3 , it can be seen that in the y direction, the velocity components of droplets at the same height in the horizontal direction, under different pressures, are close to 0, with small values and irregular state motion. It can be considered that the droplets are basically unaffected in the front-to-back direction.

TABLE 3
Fog droplet velocity component in front and rear directions.

In summary, after the two jets of the hedging nozzle impacted, the velocity component of the droplet velocity is similar in the x-axis direction, characterised by the same height direction in the spray field, and the droplet diameter should tend to the same trend (i.e. uniform particle size). The velocity of the droplets at the same height in the horizontal direction on the x-axis, shows a trend of being small in the middle to gradually increasing on both sides; the gradual increase of the velocity component of the droplets near both sides of the edge facilitates the formation of a spray fan. For the traditional fan nozzle measured in the direction of the axis, the droplet velocity shows a trend of gradually decreasing from the middle to the sides. The droplet size in the horizontal direction at the same height shows a gradual increase from a small droplet diameter in the middle to the sides, there are differences between the numerical simulation results of hedge nozzles and conventional fan nozzles.

Test verification

Testing was conducted with clean water under room temperature and windless indoor conditions, using a test system to perform nozzle atomisation ( Figure 7 ) at pressures of 0.3, 0.5 and 0.7 MPa, the pressure droplet particle size of the hedge nozzle spray field was tested, see Figure 5 .

FIGURE 7
Test system for atomisation performance of the nozzle.

Figure 8 shows that the spray pressure varies in the range 0.3-0.7 MPa and the droplet size gradually becomes smaller as the spray pressure increases in the positive direction along the axis (vertically downward), i.e. the test surface at 100, 200, 300, and 400 mm from the nozzle outlet. If the spraying pressure is 0.7 MPa, the droplet size will change from 252 μm to 148 μm. At the same level, the droplet size is larger and uniform in the middle but smaller on the two outermost sides. With the increase in spraying pressure, the droplet size at the test points on both sides will be close to the droplet size in the middle area. If the spraying pressure is 0.5 MPa on the 300 mm test surface, the droplet size between -70 mm and 70 mm is 227-231 μm. The droplet diameter at the edge of the spray fan, 105 mm away from the axis, is 210-211 μm. When the spraying pressure increases to 0.7 MPa, the droplet size between -140 mm and 140 mm is 156-165 μm, the variation range being only 9 μm. With the increase in spraying pressure, the spray amplitude of the fan-shaped surface will slightly increase but the droplet size will generally tend to be consistent at the same level. With the increase in spraying pressure, the droplets in the spray field tend to be homogenised. For example, when the spraying pressure is 0.7 MPa, at 200 mm and 400 mm horizontal heights, the droplet size tends to be close to the droplet size at 300 mm.

FIGURE 8
Particle size distribution of impingement nozzle.

The droplet distribution for the sector nozzle and impingement nozzle, measured under the same spraying pressure (0.5 MPa), is shown in Figure 9 . Curves 1 and 2 represent the frequency distribution curve and cumulative distribution curve, respectively. It is obvious from the figure that the frequency distribution of the spray droplets of the fan nozzle is wide and presents a normal distribution, while the frequency distribution of the spray droplets of the impingement nozzle is narrow, which indicates that the distribution of the spray droplets is relatively concentrated and the spectrum of the spray droplets is narrow. This is because the two jets produce intense momentum transfer on the impact surface, leading to drastic changes in energy. Collision, extrusion and shearing may occur between flow clusters, causing them to deform and split. The impact makes the large fog droplets break and the smaller fog droplets combine to form larger fog droplets. In general, the droplet size tends to be uniform, reaching homogenisation. The greater the spraying pressure and the more intense the impact, the more obvious the droplet size homogenisation.

FIGURE 9
Droplet distribution.

CONCLUSIONS

  • After the impact of the two jets, with the increase of pressure, the splitting and uniformity of the liquid flow cluster improves. This is conducive to a smaller and more uniform droplet size at the impingement nozzle, a more concentrated droplet distribution, and the formation of a narrow droplet spectrum.

  • The velocity component of the droplet in the vertical downward direction (z-axis) of the opposing jet nozzle, after the impact of the jet, is close to the same horizontal height (for example, the velocity component of the droplet is 70.33-73.70 m/s on the horizontal test surface 23 mm from the bottom of the nozzle body under a pressure of 0.5 MPa, and decreases along the positive axis direction. The velocity component of the droplet in the horizontal direction (x-axis) is small in the middle and large on both sides, and the velocity near the central axis is small and moves in an irregular direction. The velocity component of droplets in the front and rear directions (y-axis) is very small. Due to the close velocity component of the droplet velocity in the positive direction of the z-axis and the jet impact, the droplet group size will be homogenised, which is conducive to generating droplets with uniform size at the same level in the spray field. The velocity component of the droplet velocity in the x-axis direction has an effect on the spray amplitude of the nozzle.

  • The droplet size on the fan shaped surface of the opposing jet spray nozzle is larger in the middle and uniform, but smaller on both sides. For example, the droplet size in the middle is 227-231 μm when the spraying pressure is 0.5 MPa on the test surface, at a horizontal height of 300 mm. The droplet size near the edge is 210-211 μm. This characteristic is beneficial to solve the problem of uneven droplet size distribution when the spray bar operates in the field. The experimental results are in agreement with the numerical simulation and theoretical analysis results.

REFERENCES

  • Cai Y, Duan Z, Niu J, Li D, Zhang J, Qiao G (2022) Numerical simulation of flow field characteristics in impinging stream continuous crystallizer. Chemical Engineering (China) 50(11):43-49.
  • Dong F, Zhou H, Shi M, Ru Y, Yi K (2019) Design and Atomisation Test of Impinging nozzle. Scientia Silvae Sinicae 55(1): 83-90.
  • Dong F, Zhou H (2018) Development of foreign plant protection nozzles. Acta Agriculturae Universitatis Jiangxiensis 40(4):866-874.
  • Dong F, Zhou H, Ru Y Shi M, Chen Q, Yi K (2018) Experiment on droplet distribution characteristics in spray field of impinging nozzle. Transactions of the Chinese Society for Agricultural Machinery 49(12):116-121+128.
  • Fei J, Sun F, Yang W, Fu Y, Wang Y (2015) Experimental analysis on movement and size distribution of atomised droplets from impinging liquid jet. Journal of Rocket Propulsion 41(1):10-14.
  • Gu Z, Zhao W, Zhang Y, Wei H (2015) A Research on the effects of mesh quality evaluation criteria on the results of vehicle external flow field simulation[J]. Automotive Engineering 37(11):1265-1269.
  • He X (2019) Research and development of crop protection machinery and chemical application technology in china.Chinese Journal of Pesticide Science. 21(5):921-930.
  • Lü X, He X, Song J, Zeng A, Andreas H (2007) Analysis of spray process produced by agriculture flat-fan nozzles. Transactions of the Chinese Society of Agricultural Engineering 23(9):95-100.
  • Salcedo R, Ramon L, Jordi C, Javier C, Michael G, Montserrat O, Paula G, Emi L (2020) Evaluation of leaf deposit quality between electrostatic and conventional multi-row sprayers in a trellised vineyard. Crop Protection 127:104964.
  • Song S, Chen J, Hong T, Xue X, Xia H, Song Y (2017) Variation of droplet diameter in wind field for long-range air-assisted sprayer. Transactions of the Chinese Society of Agricultural Engineering 33(6):59-66.
  • Wang G (2019) Spraying Quality Evaluation on Plant Protection UAV. Guangzhou: South China Agricultural University. 2019.6.
  • Wang X, Liu J, Zhang Q (2022) Water–Pesticide Integrated Micro-Sprinkler Design and Influence of Key Structural Parameters on Performance. Agriculture 12(10):1532.
  • Zhang J, Zhang Z, Feng Y, Dong X (2022) Advances in multiscale analysis of impinging streams flow field and mixing characteristics. Journal of Engineering Thermophysics 43(11):2944-2956.
  • Zheng G, Nie W, He B, Xue C (2015) Study on the influence of impact angle on the atomisation characteristics of impinging nozzle. Journal of Propulsion Technology 36(4):608-613.

Edited by

Area Editor: João Paulo Arantes Rodrigues da Cunha

Publication Dates

  • Publication in this collection
    04 Dec 2023
  • Date of issue
    Nov 2023

History

  • Received
    13 Nov 2022
  • Accepted
    2 Oct 2023
Associação Brasileira de Engenharia Agrícola SBEA - Associação Brasileira de Engenharia Agrícola, Departamento de Engenharia e Ciências Exatas FCAV/UNESP, Prof. Paulo Donato Castellane, km 5, 14884.900 | Jaboticabal - SP, Tel./Fax: +55 16 3209 7619 - Jaboticabal - SP - Brazil
E-mail: revistasbea@sbea.org.br