## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Scientia Agricola

##
*On-line version* ISSN 1678-992X

### Sci. agric. (Piracicaba, Braz.) vol.69 no.1 Piracicaba Jan./Feb. 2012

#### http://dx.doi.org/10.1590/S0103-90162012000100002

**BIOMETRY,MODELING AND STATISTICS**

**Fractal analysis of agricultural nozzles spray**

**Francisco Agüera ^{I, *}; David Nuyttens^{II}; Fernando Carvajal^{I}; Julián Sánchez-Hermosilla^{I}**

^{I}Universidad de Almería/Escuela de Ingeniería - Depto. de Ingeniería Rural - Ctra. Sacramento s/n - 04120 - La Cañada de San Urbano Almería - España

^{II}Institute for Agricultural and Fisheries Research (ILVO)/Technology & Food Science Unit - Agricultural Engineering, Burg. Van Gansberghelaan, 115 - 9820 - Merelbeke - Belgium

**ABSTRACT**

Fractal scaling of the exponential type is used to establish the cumulative volume (V) distribution applied through agricultural spray nozzles in size x droplets, smaller than the characteristic size X. From exponent d, we deduced the fractal dimension (D_{f}) which measures the degree of irregularity of the medium. This property is known as 'self-similarity'. Assuming that the droplet set from a spray nozzle is self-similar, the objectives of this study were to develop a methodology for calculating a D_{f} factor associated with a given nozzle and to determine regression coefficients in order to predict droplet spectra factors from a nozzle, taking into account its own D_{f} and pressure operating. Based on the iterated function system, we developed an algorithm to relate nozzle types to a particular value of D_{f}. Four nozzles and five operating pressure droplet size characteristics were measured using a Phase Doppler Particle Analyser (PDPA). The data input consisted of droplet size spectra factors derived from these measurements. Estimated D_{f} values showed dependence on nozzle type and independence of operating pressure. We developed an exponential model based on the D_{f} to enable us to predict droplet size spectra factors. Significant coefficients of determination were found for the fitted model. This model could prove useful as a means of comparing the behavior of nozzles which only differ in not measurable geometric parameters and it can predict droplet spectra factors of a nozzle operating under different pressures from data measured only in extreme work pressures.

**Keywords:** pesticides application, modeling

**Introduction**

Pesticides play an important role in modern agriculture, and spray nozzles are considered a key element in the application of phytosanitary products given that their task is to generate droplets through which those products will be distributed. Droplet size may influence the biological efficacy of the pesticide applied. As pressure determines the droplet characteristics of a given nozzle, is necessary to know the ideal combination of pressure-nozzle in order to optimize spray efficiency (Zhu et al., 1994; Baetens et al., 2007, 2009; Nuyttens et al., 2007, 2010, 2011)

Knowing the D_{v0.1}, D_{v0.5}, and D_{v0.9} droplet size spectra factors corresponding with 10%, 50% and 90%, respectively, of the cumulative volume spray liquid volume contained in droplets up to the indicated diameter (ASAE Standards, 1997), of nozzles operating under different pressure, orifice characteristics or liquid flow rate could be useful in helping us to compare their behavior (Womac et al., 1999) and might also be used as inputs for computer models that would enable us to predict, for example, the dispersion and deposition of aerially released spray material (Teske et al., 2000) or the amount of spray drift (Baetens et al., 2007, 2009).

Mandelbrot (1982) derived the term 'fractal' from the Latin *fractus*, which refers to the appearance of shattered rock. Traditional Euclidean geometry describes categories of objects such as points, curves, surfaces and cubes using dimensions 0, 1, 2, and 3, respectively. Nevertheless, fractal geometry is an improvement of and a development upon Euclidean geometry because it provides new ideas and concepts for the mathematical description of highly irregular and heterogeneous media. Fractal theory postulates an intrinsic symmetry law underlying the apparent disorder of these media, which consists of the repetition of the disorder itself over a certain range of scales. This property is called 'self-similarity' and objects or sets which exhibit this property are referred to as 'fractals'.

In the field of agricultural nozzles,fractal concepts and ideas have been used by Agüera et al. (2006), to predict droplet size spectra factors.

Scaling laws of the type shown below, have been applied to the cumulative number of soil aggregates (Perfect et al., 1992)

where: *N _{x>X}* is the cumulative number of size

*x*aggregates greater than a characteristic size

*X*. Exponent D

*, called 'fractal dimension', measures the degree of irregularity of the medium. These laws are based on the assumption that the behavior of certain soil properties is invariable regardless of the scale used for their study. Thus, Eghball et al. (1993) quantified changes in soil structure by detecting changes in the D*

_{f}*values associated with them.*

_{f}Taking expression (1) into account, the number-size distribution can be inferred from the volume-size distribution function *V(x < X)* of the cumulative volume of aggregates with a characteristic size lower than *X*. Thus, if spray droplets are grouped into different classes according to their sizes, droplet size spectra function shape (droplet diameter vs. accumulated volume fraction) at this scale is defined by the volume in each of these classes. If the scale in one of these classes is changed, new and different classes will appear and it might be expected that that the initial structure will be repeated, in other words: the proportion of volume in each class will agree with the above scale. Fractal ideas could be applied to the description of such distributions based on the scale: 'invariable behavior of spray properties'. If self-similarity is accepted for droplet sets, this process could be repeated at different scales.

If we accept the assumptions of constant density and spherical shape for droplet spray, the law scale expressed in (1) may be related to

where: the relation between the exponents is *d*= 3-D* _{f}* (Tyler and Wheatcraft, 1992).

Our study, therefore, had two objectives: i) to develop a methodology for calculating a D* _{f}* factor associated with a given nozzle and independent of operating conditions; and ii) to determine regression coefficients in order to predict droplet spectra factors (D

_{v0.1}, D

_{v0.5}, and D

_{v0.9}) from a nozzle, taking into account its own D

*and pressure operating conditions.*

_{f}

**Materials and Methods**

Four nozzle types widely used in agriculture were tested: i) Teejet DG 11002 (Spraying Systems Co., Wheaton, Il), a drift guard even flat fan nozzle, which we will refer to here as 'DG-110'; ii)Teejet XR-11001 (Spraying Systems Co., Wheaton, Il), an extended range flat fan nozzle, hereafter referred to as 'XR-110'; iii) Teejet TXA 8001 (Spraying Systems Co., Wheaton, Il), a hollow cone nozzle, to which we will refer in the text as 'TXA-80'; and iv) Teejet TP-9501 E (Spraying Systems Co., Wheaton, Il), an even flat fan nozzle, hereafter referred to simply as 'TP-9501'. These nozzles are widely used in the greenhouse crop-production system in the South-East of Spain, but the methodology described in this paper could be applied to any nozzle at any operating pressure. All four nozzles were tested at five pressure levels (0.2, 0.4, 0.6, 0.8 and 1 MPa ) with three replications for each pressure. This makes a total of 5 × 3 × 4 = 60 measurements. These data were obtained using a Phase Doppler Particle Analyser (PDPA).

The PDPA laser was an Aerometrics PDPA one-dimensional system. All measurements were carried out by spraying water with a temperature of 20 ºC. Environmental conditions were kept constant at 20 ºC and a relative humidity of 60-70 %. The nozzle was positioned 0.50 m above the measuring point of the PDPA. To enable the whole of the spray cloud to be sampled, the nozzle was mounted on a transporter which allowed to move it in a transverse range of 1.5 × 1.0 m, and ∆x and ∆y of 0.1 m. Details of this instrument and measurement protocol can be found in Nuyttens et al. (2007).

The data sets derived from these measurements and used in the present study were 20 values of droplet size spectra factors (D_{v0.05}, D_{v0.1}, D_{v0.15}, ..., D_{v0.9}, D_{v0.95}, D_{v1}) for each nozzle, pressure and repetition.

D* _{f}* was calculated in two ways. Firstly, as data corresponding to droplet diameter were available from each experiment, fractal dimension for each experiment was calculated, taking into account equation (1), using the diameter droplets measured with the PDPA instrument for each nozzle, repetition and operating pressure. To do so, we arranged every data set in ascending order from minimum to maximum diameter. We then grouped each set into sub-sets of 500 droplets and their average diameter was calculated, thereby obtaining a set of average diameters. D

*of every set was computed, taking into account equation (1), which can be particularized with Mandelbrot's (1982) expression for self-similar sets*

_{f}where *k*_{1} is a constant. This can be rewritten as

Therefore, D* _{f}* will be derived from the regression coefficient of the linear regression fitted between log(

*N*) and log(

_{x>X}*X*), being

*N*the cumulative number of diameter

_{x>X}*x*droplets greater than a characteristic average diameter

*X*of a sub-set containing 500 droplets. This calculation was carried out considering all operating pressure (D

*) and considering only the extreme operating pressure (D*

_{fr}*).*

_{fr1}Secondly, D* _{f}* was estimated as follows: let I

_{1}= [D

_{min}, D

_{v1}], I

_{2}= (D

_{v1}, D

_{v2}], ..., I

_{j}= (D

_{vj-1}, D

_{vj}], ..., I

_{n}= (D

_{vn-1}, D

_{max}], which are the sub-intervals of the sizes corresponding to the n droplet size spectra factor, where D

_{min}and D

_{max}are the minimum and maximum droplet diameter, respectively. Thus, the whole range of droplet diameter is I= [D

_{min}, D

_{max}]. In addition,

*q*

_{1},

*q*

_{2}, ...,

*q*

_{j}, ...,

*q*

_{n}are the relative volume proportions or probabilities (

*q*

_{1}+

*q*

_{2}+ ...+

*q*

_{j}+ ...+

*q*

_{n}= 1) of volume intervals I

_{1}, I

_{2}, ..., I

_{j}, ..., I

_{n}, respectively.

The following functions may be considered:

where: *x*∈[D_{min}, D_{max}]. In the equations 5 to 8, Ψ_{φ}, j= {1, 2, ..., n} are the linear functions that transform the [D_{min}, D_{max}] points into I_{j}, j= {1, 2, ...,n} points. The set {Ψ_{φ}; *q*_{j}; j= {1, 2, ...,n}} is called an 'iterated function system' (Barnsley and Demko, 1985). By means of the similarities Ψ_{j} and the volume proportions *q*_{j}, an iterated function system determines how a fractal distribution reproduces its structure at different scales. Figure 1 shows a graphical representation of the performance of the iterated functions system: each function Ψ_{φ}, with j= {1, 2, ...,n} transforms values of the whole range of drop diameters [D_{min}, D_{max}] into points of the ranges I_{j}, with j= {1, 2, ...,n}. For example, function Ψ_{1} transform the values of [D_{min}, D_{max}] into values of I_{1}, taking into account equation (1), and so for the rest of Ψ_{φ}s. Furthermore, the probability of use the function Ψ_{j} is *q _{j}*, which is the volume fraction accumulated in the range I

_{j}.

This second method to estimate the fractal dimension has been developed because usually, data supplied by the nozzle manufacturer is the volumetric droplet size distribution. So, available data are D_{v}'s instead of cloud droplet diameters. Using fractal geometry, an algorithm can be defined to obtain D* _{f}* of a nozzle (Elton, 1987; Turcotte, 1992): i) take any value

*x*

_{0}of [D

_{min}, D

_{max}]; ii)take a random number j, j= {1, 2, ..., n}, being

*q*

_{j}= {

*q*

_{1},

*q*

_{2}, ...,

*q*

_{n}} their respective probabilities, and compute

*x*

_{1}= Ψ

_{j}(

*x*

_{0}); and iii) repeat step 2 taking

*x*

_{1}to compute the next value. In this way, a set S= {

*x*

_{1},

*x*

_{2}, ...,

*x*

_{m}} is obtained.

This is not a set of droplets, instead these elements allow us to obtain the droplet diameter versus the cumulative volume fraction function and from this function we can then obtain the droplet size distribution, as follows: if g is the number of *x*'s values that belong to any interval I_{j}, the ratio g/m approaches the relative volume of the interval I_{j} as the number of iterations m goes to infinite. For example, if droplet size spectra factors measurements were D_{v0.1}= 147, D_{v0.5}= 464, and D_{v0.9}= 1000 µm and the elements generated by the algorithm corresponding to this test {*x*_{1}, *x*_{2},..., *x*_{m}} are sorted by 'size', there will be 0.1 × m elements with a size equal or less than 147 (D_{v0.1}); 0.5 × m elements with a size equal or less than 464 (D_{v0.5}); and 0.9 × m elements with a size equal or less than 1000 (D_{v0.9}).

D* _{f}* of this set will be computed taking into account equation (2), which can be particularized as follows

where k is a constant. It can be rewritten as

Therefore, D* _{f }*= 3-d will be derived from the regression coefficient of the linear regression fitted between log(

*N*) and log(

_{x<X}*X*).

The number m of elements generated in each application of this algorithm were 10000, and the initial value for *x* (*x*_{0}) was D_{v0.9}. Similar results could be expected if m __>__ 3000 and *x*_{0} takes any value of [D_{min}, D_{max}] (Taguas et al., 1999).

A total of ten iterated function systems differing in number and combination of droplet size spectra factors (Table 1) were used to estimate D* _{f}* (D

*) in order to study the accuracy of the algorithm described when applied to each of these combinations. Fractal dimension was also estimated in two ways: 1) by taking into account all operating pressures (D*

_{fe}*); and 2) by taking into account only the extreme operating pressure of 2 and 10 bar (D*

_{fe}*).*

_{fe1}The size frequency distribution of every set S_{q}= {x_{1}, x_{2}, ..., x_{m} } was calculated by dividing the size range into 20 classes (Steel and Torrie, 1980). The number of elements in each class was counted and the linear regression between log(N_{ > x}) and log(x) was carried out using the droplet size spectra factors corresponding to extreme operating pressures (2 and 10 bar). This task was carried out by a computer program developed by the authors using Visual Basic 6.0.

Droplet size spectra factors were predicted using an exponential function as follows:

where: D_{vh} can be D_{v0.1}, D_{v0.5 }or D_{v0.9 }(µm) of a given nozzle, p is the operating pressure (bar), and D* _{fe1}* is the fractal dimension, estimated as described above ('A fractal approach to droplet size distribution' section). The data used to carry out this regression were those from operating pressures of 0.4, 0.6 and 0.8 Mpa, which were not used to estimate D

*.*

_{fe1}

**Results and Discussion**

For a given nozzle, differences between functions representing each operating pressure presented no effect (Figure 2). Nevertheless, the DG-110 (Figure 2a) behavior at 0.8 and 1 MPa was very similar, and was quite different from its behavior at other pressures. Furthermore, the TXA-80 nozzle (Figure 2c), shows differences with different operating pressures. In general, a higher spray pressure corresponded to a smaller droplet size spectrum. These results are similar to those found by Womac (2001) working with nozzles under a range of pressures. Furthermore, when the nozzles were compared, the cumulative volumetric droplet size distributions were quite different. All these affirmations about the cumulative droplet size distribution are based on an ANOVA carried out comparing D_{v10}, D_{v25}, D_{v50}, D_{v75} and D_{v90} of each nozzle working at 0.2, 0.4, 0.6, 0.8 and 1 MPa (data not presented). Similarity cumulative volumetric droplet size distribution behavior was related with nozzle type: TP-9501 and XR-110 are flat fan nozzles whereas DG-110 is a drift guard flat fan and TXA-80 is a hollow cone nozzle. Droplet sizes vary from a few micrometers up to some hundreds of micrometers depending on the nozzle type. Cone nozzle (TXA-80), produced the finest droplet size spectrum, followed by flat fan nozzles (TP-9501 and XR-110) and drift guard nozzle (DG-110). Nuyttens et al. (2007) observed the same behavior for these kinds of nozzles under a range of pressures.

Fractal dimension, calculated taking into account equation (1), can be considered the 'real D* _{f}*' given that it is derived from the droplet diameter measurements from each experiment. Figure 3 shows an example of relationship between log(D) and log(

*N*

_{d>D}), and the fitted line for each tested nozzle. To clarify these graphics, only 10 % approximately of data of each nozzle have been plotted. Anyway, the regressions have taken into account the whole data set. Although data had not a perfect straight tendency, which means no perfect fractal behavior, all studied cases (nozzle, operation pressure and replication) shown similar behavior: the relationship between variables was correlated (

*p <*0.01). Furthermore, all regression coefficients indicate that the model as fitted explained more than 90 % of the variability in dependent variable (log(

*N*

_{d>D})).

Table 2 shows the average 'real fractal dimension' for each nozzle, operating pressure and repetitions studied, which have been calculated taking into account measured droplet diameters. The D* _{fr}* column shows the D

*average calculated based on the values corresponding to the five operating pressures, while the D*

_{f}*column shows the D*

_{fr}*average calculated from the extreme pressures (0.2 and 1 MPa). Since D*

_{f}*is not a measurable physical magnitude, we have decided to report these values to an accuracy of four decimal places, although we have not looked at the effect of decimal places on our results.*

_{f}

Values shown in each row are different (*p < *0.05), whereas the values shown in each column are not different (*p *> 0.05) (Table 2). Thus, the D* _{f}* calculated in this way and associated with nozzles, results in values which are independent of operating pressure. Similar results were found by Agüera et al. (2006) when working with hydraulic nozzles with variable geometry.

Table 3 shows the D* _{fe}* for every nozzle and D

_{v}combination (Table 1) when the algorithm based on iterated function systems described as the second method in 'Calculating the fractal dimension' section was applied. The column labeled 'sd' indicates the standard deviation of D

*derived from the five operating pressures and the three repetitions, whereas r*

_{fe}^{2}

_{max}and r

^{2}

_{min}indicate the maximum and minimum regression coefficients of linear regression between log(

*N*) and log(

_{x<X}*X*). Column D

*represents the fractal dimension for every nozzle, considering data from operating pressures of 0.2 and 1 MPa, and D*

_{fe1}_{v}combinations when the algorithm based on iterated function was applied. For a given nozzle and column, D

*or D*

_{fe}*values followed by the same letter are not different (*

_{fe1}*p <*0.05). D

*and D*

_{fe}*estimated taking into account D*

_{fe1 }_{v}combination 1 can be regarded as being the nearest to the 'real D

*' shown in Table 2, given that the 20 D*

_{f}_{v}values were used with the iterate function system. In fact, the values shown in columns D

*or D*

_{fr}*of Table 2 are not different (*

_{fr1}*p <*0.05) from those of Table 3 corresponding to combination 1. Relationship between log(

*N*) and log(

_{x<X}*X*) were correlated (

*p <*0.01). Furthermore, all regression coefficients indicate that the model as fitted explained the variability in dependent variable (log(

*N*

_{d>D})) between 0.826 (observed in the TXA-80 nozzle) and 0.979 (observed in the DG-110 nozzle).

D_{v} combination numbers 2, 4 and 5 yielded D* _{fe}* values which were different from those of the remaining D

_{v}combinations. Combination 2 includes only three D

_{v}values: D

_{v0.1}, D

_{v0.5}and D

_{v0.9}, which are probably not enough to obtain precise D

*values. Combination 4 adds two values of D*

_{f}_{v}: D

_{v0.15}and D

_{v0.85}, with respect to combination 2, while combination 5 adds four values with respect to combination 2: D

_{v0.15}, D

_{v0.25}, D

_{v0.75}and D

_{v0.85}. All these values are between D

_{v0.1}and D

_{v0.9}, situated in the constant-slope zone of the curve diameter vs. cumulative volume. Similar results were obtained for combination 3 which includes D

_{v}values in these areas, or combinations which include values between D

_{v0.25}and D

_{v0.75}, except combination 1. Nevertheless, a combination including D

_{v0.1}, D

_{v0.5}, D

_{v0.9}and other values around D

_{v0.1}and D

_{v0.9}(combinations 7 and 10) yielded the D

*values nearest to that yielded by combination 1. In Table 3, combinations 1, 7 and 10 are those which are included only in group 'a'.*

_{fe}D_{v} combinations are not as influential in the estimated fractal dimension as in the D* _{fe}* column: the D

*value for the XR-110 and TP-9501 nozzles was independent of D*

_{fe1}_{v}combinations, while the D

*value for the DG-110 and TXA-80 nozzles does not show differences with the D*

_{fe1}*values (Table 3). Using iterated function systems and several combinations of soil particles size spectra factors to simulate and to test of self-similar structures for soil-size distributions, Taguas et al. (1999) found relationship between the used combination and the results accuracy.*

_{fe}As the D* _{fe}* values from D

_{v}combination 10 are the nearest to those from D

_{v}combination 1, except for the TXA-80 nozzle, D

*values from D*

_{fe1}_{v}combination 10 were used to fit data in equation (11). Coefficients of determination and regression coefficients as determined by regression modeling of droplet size factors were:

where: D* _{fe1}* is the own fractal dimension value of each nozzle shown in Table 3 for combination 10, D

_{v0.1}, D

_{v0.5 }or D

_{v0.9 }are µm) and p is the operating pressure (0.1×MPa).

Coefficients of determination which were significant (*p* < 0.05) tended to have decreasing r^{2} values with increasing droplet size factors, which is a similar behavior to that observed by Agüera et al. (2006) using polynomial and multilinear fitted regressions and also to that observed by Womac (2001) using a multilinear fitted regression, although in both these studies the authors found some non-significant coefficients of determination.

**Conclusions**

The algorithm resulted in D* _{f}* values similar to those calculated from measured droplet diameter ('real D

*'), regardless of the operating pressure and related to the nozzle type. D*

_{f}*for a given nozzle can be calculated from two droplet size spectra measurements performed at extreme operating pressure (0.2 and 1 MPa), thus eliminating the need for droplet size measurements at intermediate operating pressures. D*

_{f}_{v0.1}, D

_{v0.5}and D

_{v0.9}for the nozzles tested at 0.2 and 1 MPa pressure can be estimated at any operating pressure from this exponential model.

**Acknowledgements**

This work was supported by grant P07-AGR-02995 from CICE - Junta de Andalucia (Spain), co-financed with FEDER funds of the European Union.

**References**

Agüera, F.; Aguilar, F.J.; Aguilar, M.A.; Carvajal, F. 2006. Atomization characteristics of hydraulic nozzles using fractal geometry. Transaction of the ASABE 49: 581-587. [ Links ]

ASAE Standards. 1997. S327.2: Terminology and definitions for agricultural chemical applications. ASAE Standards. 44th ed. St. Joseph, Mich.: ASAE. [ Links ]

Baetens, K.; Nuyttens, D.; Verboven, P.; De Schampheleire, M.; Nicolaï, B.; Ramon, H. 2007. Predicting drift from field spraying by means of a 3D computational fluid dynamics model. Computers and Electronics in Agriculture. 56: 161-173. [ Links ]

Baetens, K.; Ho, Q. T.; Nuyttens, D.; De Schampheleire, M.; Melese Endalew, M.; Hertog, M.; Nicolaï, B.; Ramon, H; Verboven, P. 2009. A Validated 2-D diffusion-advection model for prediction of drift from ground boom sprayers. Atmospheric Environment 43: 1674-1682. [ Links ]

Eghball, B.; Mielke, N.M.; Calvo, G.A.; Wilhelm, W.W. 1993. Fractal description of soil fragmentation for various tillage methods and crop sequences. Soil Science Society of America Journal 57: 1337-1341. [ Links ]

Elton, J. 1987. An ergodic theorem for iterated maps. Journal of Ergodic Theory and Dynamical Systems 7: 481-488. [ Links ]

Mandelbrot, B.B. 1982. The fractal geometry of nature. W.H. Freeman, New York, NY, USA. [ Links ]

Nuyttens, D.; Baetens, K.; De Schampheleire, M.; Sonck, B. 2007. Effect of nozzle type, size and pressure on spray droplet characteristics. Biosystems Engineering 97: 333-345. [ Links ]

Nuyttens, D.; De Schampheleire, M.; Verboven, P.; Sonck, B. 2010. Comparison between indirect and direct spray drift assessment methods. Biosystems Engineering. 105: 2-12. [ Links ]

Nuyttens, D.; De Schampheleire, M.; Baetens, K.; Brusselman, E.; Dekeyser, D.; Verboven, P. 2011. Drift from field crop sprayers using an integrated approach: results of a 5 year study. Transactions of the ASABE 54: 403-408 [ Links ]

Perfect, E.; Bird, N.; Rieu, M. 1992. Generalizing the fractal model of soil structure: the pore-solid fractal approach. Geoderma 88: 137-164. [ Links ]

Steel, R.G.D; Torrie, J.H. 1980. Principles and Procedures of Statistics: biometrial approach. McGraw-Hill, New York, NY, USA. [ Links ]

Taguas, F.J.; Martín, M.A.; Perfect, E. 1999. Simulation and testing of self-similar structures, for soil particle-size distributions using iterated function systems. Geoderma 88: 191-203. [ Links ]

Turcotte, D.L. 1992. Fractals and Chaos in Geology and Geophysics. Cambridge University Press, Cambridge, UK. [ Links ]

Tyler, S.W.; Wheatcraft, S.W. 1992. Fractal scaling for soil particle-size distributions: analysis and limitations. Soil Science Society of America Journal 56: 362-369. [ Links ]

Teske, M.E.; Thistle, H.W.; Mickle, R.E. 2000. Modeling finer droplet aerial spray drift and deposition. Applied Engineering in Agriculture 16: 351-357. [ Links ]

Womac, A.R. 1999. Measurements variations in reference sprays for nozzle classification. Transaction of ASAE 42: 609-616. [ Links ]

Womac, A.R. 2001. Atomization characteristics of high-flow variable-orifice flooding nozzles. Transaction of ASAE 44: 463-471. [ Links ]

Zhu, H.; Reichard, D.L.; Fox, R.D.; Brazee, R.D.; Ozkam, H.E. 1994. Simulation of drift of discrete size of water droplets from field sprayers. Transaction of ASAE 37: 1401-1407. [ Links ]

*** Corresponding author**

faguera@ual.es

**Edited by:** José Euclides Stipp Paterniani

Received June 09, 2010 Accepted September 29, 2011