Next Article in Journal
Crop Organ Segmentation and Disease Identification Based on Weakly Supervised Deep Neural Network
Next Article in Special Issue
An Integrated Yield Prediction Model for Greenhouse Tomato
Previous Article in Journal
Nondestructive Estimation of the Chlorophyll b of Apple Fruit by Color and Spectral Features Using Different Methods of Hybrid Artificial Neural Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Semi-Empirical Ventilation Models in A Mediterranean Greenhouse with Opposing Thermal and Wind Effects. Use of Non-Constant Cd (Pressure Drop Coefficient Through the Vents) and Cw (Wind Effect Coefficient)

by
Alejandro López-Martínez
1,*,
Francisco D. Molina-Aiz
1,
Diego L. Valera-Martínez
1,
Javier López-Martínez
1,
Araceli Peña-Fernández
1 and
Karlos E. Espinoza-Ramos
2
1
CIAIMBITAL Research Centre, University of Almería, Ctra. de Sacramento s/n, 04120 Almería, Spain
2
Department of Engineerings, University Center of the South Coast, University of Guadalajara. Av. Independencia Nacional 151, Autlán de Navarro, 48900 Jalisco, Mexico
*
Author to whom correspondence should be addressed.
Agronomy 2019, 9(11), 736; https://doi.org/10.3390/agronomy9110736
Submission received: 14 October 2019 / Revised: 5 November 2019 / Accepted: 7 November 2019 / Published: 10 November 2019
(This article belongs to the Special Issue Greenhouse Technology)

Abstract

:
The present work analyses the natural ventilation of a multi-span greenhouse with one roof vent and two side vents by means of sonic anemometry. Opening the roof vent to windward, one side vent to leeward, and the other side vents to windward (this last vent obstructed by another greenhouse), causes opposing thermal GT (m3 s−1) and wind effects Gw (m3 s−1), as outside air entering the greenhouse through the roof vent circulates downward, contrary to natural convection due to the thermal effect. In our case, the ventilation rate RM (h−1) in a naturally ventilated greenhouse fits a second order polynomial with wind velocity uo (RM = 0.37 uo2 + 0.03 uo + 0.75; R2 = 0.99). The opposing wind and thermal effects mean that ventilation models based on Bernoulli’s equation must be modified in order to add or subtract their effects accordingly—Model 1, in which the flow is driven by the sum of two independent pressure fields G M 1 = | G T 2 ± G w 2 | , or Model 2, in which the flow is driven by the sum of two independent fluxes G M 2 = | G T ± G w | . A linear relationship has been obtained, which allows us to estimate the discharge coefficient of the side vents (CdVS) and roof vent (CdWR) as a function of uo [CdVS = 0.028 uo + 0.028 (R2 = 0.92); CdWR = 0.036 uo + 0.040 (R2 = 0.96)]. The wind effect coefficient Cw was determined by applying models M1 and M2 proved not to remain constant for the different experiments, but varied according to the ratio uo/Tio0.5 or δ [CwM1 = exp(−2.693 + 1.160/δ) (R2 = 0.94); CwM2 = exp(−2.128 + 1.264/δ) (R2 = 0.98)].

1. Introduction

Natural ventilation is perhaps the main means of climate control in greenhouses [1,2], particularly in regions such as the Mediterranean, where new technologies have not been widely incorporated [3]. During most of the year, good management of natural ventilation may prove sufficient to maintain suitable levels of temperature, humidity, and CO2 concentration inside the greenhouse [1]. It is important to have both quantitative and qualitative knowledge of how natural ventilation functions in greenhouses in order to correctly design and use the vents [4].
The earliest studies on the circulation of air in greenhouses date back to the mid-20th century [5,6]. Since then, many authors have shown interest in studying and understanding natural ventilation in greenhouses. Numerous methods have been used: scale models, tracer gas method, CFD (Computational Fluid Dynamics) simulations, and direct measurements in the greenhouse using many types of sensors.
In a greenhouse with side and roof openings, the main driving forces of natural ventilation are caused by a combination of pressure differences [4,7] motivated by: (i) the static wind effect due to the mean component of the wind velocity, which generates pressure differences between the side and roof openings [8] and between the windward and leeward parts of the greenhouse [1]; (ii) the buoyancy forces (stack or chimney effect) generating a vertical distribution of pressures between the side and roof vents [9]; (iii) and the turbulent effect of the wind, which generates pressure fluctuations along and across the greenhouse openings [1,7].
In summary, natural ventilation of the greenhouse is the result of combining the airflows generated by the wind (wind/eolic effect) and the buoyancy of the air (thermal/chimney effect). These airflows can be combined, producing different models derived from Bernoulli’s equation [4,7,10]. Applying these semi-empirical models requires using two fundamental parameters: the pressure drop coefficient through the vents (vent coefficient, Cd) and the coefficient of the wind effect (eolic coefficient, Cw), which comprises of both the stationary and turbulent aspects of the wind’s action. These parameters may be obtained by field experiments, which determine the ventilation rate.
Though many researchers have obtained the coefficients Cd and Cw by applying semi-empirical models [2,4,10,11,12,13,14,15,16,17,18,19], there is considerable uncertainty regarding the relationship between the wind and thermal effects. Papadakis et al. [11] found that the thermal effect was fundamental for winds <1.8 m s−1, while Boulard and Baille [7] speak of 1.5 m s−1 and Sase et al. [20] speak of 1 m s−1. According to Boulard and Baille [7] and Fatnassi et al. [13], in greenhouses with roof vents, the thermal effect is thought to be negligible for air velocity > 2 m s−1. Also, in greenhouses with side and roof vents, Kittas et al. [4] established that the thermal effect is important when the ratio uo/∆Tio0.5 is less than 1, while Bot [21] put this limit at 0.3.
The relationship between wind and thermal effects must be known in order to correctly design greenhouse ventilation. Few authors mention this fact, and most works apply models that take these effects as being complementary to one another. This is not always the case, however, as depending on the arrangement of ventilation openings, wind can assist or oppose the thermal force in natural ventilation. It is most important to bear this in mind when applying semi-empirical models of ventilation [22].
To study the natural ventilation in greenhouses, the gas tracer method can allow us to quantify ventilation rate correctly, but not characterise the real airflow inside the greenhouse [1]. One disadvantage of the tracer gas method is that the position of the sensors can affects the results obtained, as Van Buggenhout et al. [23] observed, with errors of up to 86% according to the position of the sensors.
With sonic anemometers, we can measure the different components of the air velocity vector at the greenhouse vents, determining where the air enters or exits the greenhouse through each vent. This methodology allows us to study when the wind and thermal effects are complementary or opposing. The objective of the present work is to study natural ventilation in a greenhouse with two side vents and one roof vent with the use of sonic anemometry. The roof vent was windward, the windward side vent was protected from the wind by another greenhouse located in close proximity, and the leeward side vent was free of obstacles. Under such conditions, wind and thermal effects maybe fail to contribute to natural ventilation in the same manner. Rather, they can have opposing effects.
In this work, different ventilation models have been applied in a greenhouse with opposition of the wind and thermal effects. The models have been modified to take into account the opposition of both effects. Identifying which model best fits the experimental values was done by statistical analysis.

2. Materials and Methods

2.1. Experimental Setup

The experimental work took place in a three-span Mediterranean greenhouse (1080 m2) located at the agricultural research farm belonging to the University of Almería (36°51′ N, 2°16′ W, and 87 MASL). The greenhouse was divided into two independent halves (Figure 1). In this work, we studied the natural ventilation in the western half (24 × 20 m2), which has less obstruction to the wind.
In the western sector of the greenhouse, the vent surface area was 1.05 × 17.5 m2 for each of the two-side vents (2.1 m in height from the ground to the mid-point of the vent) and 0.97 × 17.5 m2 for the roof vent (6.2 m in height from the ground to the mid-point of the vent), and 4.1 m of vertical distance between the mid-point of the side and the roof vents. The ventilation surface SV was 11.2% of the greenhouse base area (SV/SA = 0.112). To prevent insects entering, all the greenhouse vents were fitted with insect-proof screens of 13 × 30 threads cm−2 (0.39 porosity; 164.6 µm pore width; 593.3 µm pore height; 165.5 µm thread diameter). Measurement tests were carried out under prevailing winds from the southwest (Figure 1), popularly known as “Poniente wind” (in the province of Almería). The outside climatic conditions remained relatively stable over the ten measurement tests (Table 1). The greenhouse contained a tomato crop (Solanum lycopersicum L. var. cerasiforme Hort., cv. Salomee) with an average height of approximately 0.85 m for the first test and 1.88 m for the final one, and a leaf area index (m2 leaf per m2 ground.) of about 0.37 and 2.75, respectively. The measurements were performed from 7 April 2009 to 2 July 2009.
Airflow was measured at 21 points (tests 1 to 4) and 12 points (tests 5 to 10), evenly distributed at the side vents (Figure 2a,b). To measure air velocity at the side vents, two 3D sonic anemometers were used (Figure 2c), recording data over three minutes at each point [19].
Given the difficulty of placing the sonic anemometers in the roof vent, it was divided into three equal surfaces and air velocity measurements were taken continuously at the centre of each one, with two (measurement tests 1 to 4) or one (measurement tests 5 to 10) 2D sonic anemometers (Figure 2a,b). At the roof vent, these measurements were taken continuously using six (tests 1 to 4) or three (tests 5 to 10) 2D sonic anemometers (Figure 2c).
The distribution of measurement points used in the present study was very similar to the one made in Espinoza et al. [24], in the same experimental greenhouse. For lateral vents, the surface corresponding to each point was 0.9 m2 (tests with 21 points) and 1.5 m2 (tests with 12 points). These are similar values to those used by other authors. Boulard et al. [25] used 2.6 m2 per point in a tunnel greenhouse with one roof vent, Teitel et al. [18] used 1.1 m2 per point in a mono-span greenhouse with two side vent openings, and Molina-Aiz et al. [19] used 2.1 m2 per point in a five-span Almería-type greenhouse. At the roof vent, the mean surface corresponding to each section was 5.8 m2. Teitel et al. [26] used a larger surface area for each measurement point, with 8.5 m2 per point in a four-span greenhouse with three roof vents.

2.2. Equipment and Instrumentation

The three components of air velocity and air temperature were measured with two 3D sonic anemometers. The horizontal components of air velocity (x and y) were measured with six 2D sonic anemometers (Table 2). The data measured by the sonic anemometers were recorded by two CR3000 Microloggers (Campbell Scientific Spain S.L.). The data registration frequency was 10 Hz [27] for 3D sonic anemometers and 1 Hz for 2D sonic anemometers, respectively.
Outside climatic conditions were recorded by a meteorological station at a frequency of 0.5 Hz (Figure 1). It included a measurement box with a Pt1000 temperature sensor and a capacitive humidity sensor (BUTRON II). Wind speed and direction were recorded by a cup anemometer and a vane. Solar radiation was measured using a Kipp Solari sensor. Air temperature and humidity inside the greenhouse were measured using six autonomous dataloggers HOBO Pro Temp-RH U23-001 (Table 2). The dataloggers were placed in a vertical profile under the ridge of the three greenhouse spans at heights of 1 and 2 m (Figure 1). These autonomous dataloggers were protected against direct solar radiation with a passive solar radiation open shield.

2.3. Ventilation Models

In general, we can consider that the natural ventilation in a greenhouse is the result of combining two different fluxes: one generated by wind and the other generated by buoyancy forces. We can combine these fluxes and obtain different semi-empirical models based on Bernoulli’s equation [4,7,19,28]. These models (used by Molina-Aiz et al. [19]) can be modified by applied using a different discharge coefficient for the side vents (CdVS) and for the roof vent (CdVR):
(1)
Model 1. In this model, the flow is driven by the sum of two independent pressure fields. In a greenhouse with side and roof openings, G is calculated as the vector sum of the free component of the flux induced by buoyancy forces GT and the flux induced by wind forces Gw:
G M 1 = ( ( C d V R × S V R ) ( C d V S × S V S ) ( C d V R × S V R ) 2 + ( C d V S × S V S ) 2 ) 2 ( 2 g Δ T i o T o h S R ) + ( C d V R × S V R + C d V S × S V S 2 ) 2 C w × u o 2
where SVR and SVS are the roof and side total vent areas, respectively [m2], g is the gravitational constant [m s−2], ∆Tio is the inside-outside temperature difference [K], To is the outside temperature [K], uo is the outside wind speed [m s−1], Cd is the discharge coefficient (CdVR, roof vent; CdVS, side vents), Cw is the wind effect coefficient, and hSR is the vertical distance between the mid-point of the side and the roof vents [m].
(2)
Model 2. In this model, the flow is driven by the sum of two independent fluxes. The greenhouse volumetric flow rate is calculated as the algebraic sum of the free component of the flux induced by buoyancy forces GT and the flux induced by wind forces Gw:
G M 2 = [ ( C d V R × S V R ) ( C d V S × S V S ) ( C d V R × S V R ) 2 + ( C d V S × S V S ) 2 2 g Δ T i o T o h S R + C d V R × S V R + C d V S × S V S 2 C w u o ]
(3)
Model 3. In this model, only the wind effect is considered, neglecting the thermal effect:
G M 3 = C d V R × S V R + C d V S × S V S 2 C w u o
The height of the meteorological station and its distance from the greenhouse condition the results obtained when applying the models. This aspect is important to keep in mind when applying the models described in the literature.

2.4. Anemometric Measurement of Volumetric Flow Rate

To describe the air circulation through the greenhouse vents, the mean and turbulent volumetric flow rates can be calculated considering only the component of air velocity perpendicular to the vent, ux [25]:
G ¯ j = j = 1 n ( S V , j × u x , j )
G j = j = 1 n ( S V , j × | u x , j | )
where ux,j are the time average value of air velocity perpendicular to the vent, |u’x,j| are the time average value of the absolute values of the turbulent component, and SV,j are elementary surface corresponding to each measurement point. With only two 3D sonic anemometers measuring at the lateral vents, at different intervals of 3 min (Figure 2), the external conditions can change or fluctuate during the test. This problem can be overcome by correcting (scaling with the wind speed) the air velocities measured at each position j at the lateral vents [19]. The corrected air velocity u*x,j(t) has been calculated by multiplying measured values of ux,j(t), at time t and at each point j, by the ratio between the average wind speed uo for the overall test period, and the measured values of uo(t) at time t [19]:
u x , j * ( t ) = u x , j ( t ) u o u o ( t )

2.5. Estimation of the Pressure Drop/discharge Coefficient of the Greenhouse Openings

The individual pressure drop coefficient of each vent Cd was determined following the methodology outlined by Molina-Aiz et al. [19]. Calculation of Cd requires knowing the specific permeability (Kp = 1.851·10−9) and inertial factor (Y = 0.155) of the insect-proof screens installed in the vents; these values were determined in wind-tunnel tests. A detailed description of the wind-tunnel experiment was provided by Valera et al. [29] or López et al. [30]. To determine Kp and Y, the thickness of the screen is required (e = 391.7 µm), which was obtained by measuring a transversal section of the net with an optical measurement unit equipped with a video system (TESA-VISIO 300, TESA SA, Switzerland; with a resolution of 0.05 μm and precision of ±7 μm).

2.6. Statistical Analysis

We carried out regression analyses to compare the different variables for statistically significant relationships (p-value < 0.05) using Statgraphics®Centurion 18 v18.1 (Statgraphics Technologies, Inc., The Plains, VA, USA). In order to analyse the fit of the average volumetric flow rate for the greenhouse GM,sim (simulated by models M1, M2 and M3) with the experimentally observed values of GM,obs, different statistics were used, as well as the determining coefficient R2. Two of the most commonly used statistics based on the deviation of GM,sim from GM,obs are RMSD (root mean squared deviation) and MD (or bias) [31]:
R M S D = 1 n i = 1 n ( G M , o b s G M , s i m ) 2
M D = 1 n i = 1 n ( G M , o b s G M , s i m )
RMSD represents the average distance between the simulated and observed values. MD corresponds to the mean value of the differences between the simulated and observed values, though in this case negative differences are compensated by positive ones, which may lead to erroneous interpretation. Both statistics represent different aspects of the deviation from simulated values, but the relationship between them has not been well defined [31]. The Nash-Sutcliffe efficiency (NSE) is a normalised statistic indicating the relative magnitude of the residual variance of the model (“noise”) with respect to the variance of the observed values (“information”) [32]. NSE indicates how well the plot of the observed versus simulated values fits the 1:1 line [33]:
N S E = 1 [ i = 1 n ( G M , o b s G M , s i m ) 2 i = 1 n ( G M , o b s G ¯ M , o b s ) 2 ]
where G ¯ M , o b s is the average value of all the GM,obs values. NSE can take values of between −∞ and 1.0, the latter being the optimal value. Values between 0.0 and 1.0 are generally viewed as acceptable levels of performance [33]. Percent bias PBIAS represents the mean trend of simulated values to greater or lower than their respective observed values [34]. The optimal value of PBIAS is 0.0, and values close to 0 indicate high precision of the model. Positive values indicate that the model provides values that are lower than those observed (model underestimation bias), while negative ones indicate the contrary (model overestimation bias) [34]. PBIAS can be expressed as a percentage:
P B I A S = [ i = 1 n ( G M , o b s G M , s i m ) × 100 i = 1 n ( G M , o b s ) ]
RMSD-observations standard deviation ratio (RSR): this statistic [33] was developed based on the recommendations of [35]. RSR standardises the values of RMSD with standard deviation from the observed values. It combines both an error index and the additional information recommended by Legates and McCabe [36]:
R S R = R M S D S T D E V o b s = [ i = 1 n ( G M , o b s G M , s i m ) 2 i = 1 n ( G M , o b s G ¯ M , o b s ) 2 ]
RSR takes values from 0 to a large positive value, and when RMSD is equal to zero, the model is considered perfect. Small values of RSR indicate a better performance of the simulation model [33].

3. Results and Discussion

3.1. Airflow Inside the Greenhouse

Given the particular situation of the greenhouse and the location of the vents, in conditions of natural ventilation with prevailing southwest wind (SW) or “Poniente winds”, the eolic and thermal effects oppose each other. The former causes air to enter through the windward roof vent and to leave through the leeward and windward side vents (the windward side vent was protected from the wind by another greenhouse located in close proximity). The thermal effect, on the other hand, causes warm air to rise and leave through the roof vent, which favours the entrance of air through the side vents. Figure 3a,b shows the entrance and exit of air through the roof vent (see the polar histograms), demonstrating the negative interaction of the wind and thermal effects. This flow pattern with opposition of the wind effect and thermal effect in the roof vents, in the same naturally-ventilated three-span Mediterranean greenhouse, was described in more detail, including measurements among the crop lines, in López et al. [37]. And, similar flow patterns, in the same greenhouse and with the same methodology, was described in Espinoza et al. [24], opening the lateral vents combined with two and three roof vents.
The greenhouse natural ventilation rate is affected by the thermal effect due to the inside-outside temperature difference, ∆Tio. For greenhouses with side and roof vents, Kittas et al. [4] established that the thermal effect is important when the ratio uo/Tio0.5 < 1, while Bot [21] put this limit at 0.3. In tests 1, 2, 3, 7, and 10, with uo/Tio0.5 ≥ 1.5, hardly any air left the greenhouse through the roof vent, indicating that the wind effect clearly prevailed over the thermal effect, and the air left in quite a uniform fashion through the side vents (Figure 4).
In the other tests, with uo/Tio0.5 < 1.5, the opposing wind and thermal effects gave rise to less uniformity of air flow at the ventilation surfaces (Figure 3a,b), with alternating positive (entrance) and negative (exit) airflow at the roof vent (see polar histograms). In some of the tests, there was a certain degree of discrepancy between the wind direction and the airflow direction entering the greenhouse through the roof vent. This is likely due to the location of the meteorological station, to the north of the greenhouse, which therefore recorded the wind characteristics once it had passed through the experimental greenhouse.

3.2. Evaluation of the Mean and Turbulent Ventilation Flows

The longitudinal component ux at the side vents is corrected according to Equation 6 to compensate for the change in wind speed during the tests. This did not prove necessary for the roof vent, as the measurements were taken continuously (Table 3).
By applying Equations (4) and (5), we determined the volumetric flow rates at the three greenhouse vents. Knowing the volumetric flow rate of the greenhouse, GM allows the ventilation rate RM to be calculated (Table 4 and Figure 5a).
The precision of the mean values of GM can be assessed by comparing them with the inflow and outflow volumes measured at the different ventilation surfaces. To verify the degree of satisfaction of the Mass Conservation Law in the greenhouse [1,28], the error in the calculation of the ventilation volumes has been estimated as follows:
E G = Δ G G M = G L S ( u x * ) + G W S ( u x * ) + G W R ( u x ) G M × 100
If we calculate the airflow without correcting the inside air velocity with the outside wind speed, we obtain a mean error for all the assays of EG = 17.2%, but on correcting it, the mean error of the assays is reduced to EG = 15.1%. Following a similar methodology in an Almería greenhouse with two side and two roof vents, Molina-Aiz et al. [19] obtained errors of EG between 3.0% and 37.0% (calculating airflows with ux*) and between 8.4% and 70.5% (with ux). In a multi-tunnel greenhouse with a roof vent, other authors obtained errors of 2.2% and 2.6% [1], and of 31.6% [28].
The accuracy of the method of calculating the ventilation volumetric flow rates will depend mainly on the stability of the wind conditions (intensity and direction) and on the influence of the thermal effect on greenhouse ventilation. When wind conditions are not stable or when the dominance of the thermal effect or the wind effect is not clear, a greater degree of error should be expected. Maximum error was recorded in the assays for which the ratio uo/Tio0.5 was close to 1.5 (Figure 5b).
The greenhouse ventilation rate RM under conditions of natural ventilation (Table 4) presents a better fit to a second order polynomial with the air velocity (RM = 0.37 uo2 + 0.03 uo +0.75; R2 = 0.99 and p-value = 0.0000) (Figure 5a). Many authors present linear fits between wind velocity and the ventilation rate [4,7,18]. If the fit is carried out in this way, the R2 coefficient obtained is somewhat lower (RM = 3.20 uo – 5.04; R2 = 0.95 and p-value = 0.0000) and the resulting straight line would cross the abscissa at uo = 1.61 m s−1 (Figure 5a). According to this fit, for air velocity below this value, there would be no exchange of air with the outside, and so this cannot be considered a valid fit.

3.3. Application of the Semi-empirical Ventilation Models Based on Bernoulli’s Equation

3.3.1. Determining the Discharge Coefficient of the Vents Cd

The discharge coefficients (Table 5) have been obtained for each of the vents, following the methodology outlined by Molina-Aiz et al. [19]. The windward roof vent had the highest average discharge coefficient (CdWR = 0.161 ± 0.054), followed by the windward side vent (CdWS = 0.143 ± 0.023), and finally the leeward side vent (CdLS = 0.100 ± 0.062); the average discharge coefficient for both side vents was CdVS = 0.121 ± 0.042. The discharge coefficients of side vents CdVS and of roof vent CdWS can be expressed as a function of air velocity [CdVS = 0.028·uo + 0.028 (R2 = 0.92 and p-value = 0.0000); CdWR = 0.036·uo + 0.040 (R2 = 0.96 and p-value = 0.0000)] (Figure 6).

3.3.2. Determining the Wind Effect Coefficient Cw

Once the discharge coefficient of the greenhouse vents, CdVS and CdWR (Table 5), and the average volumetric flow rate, GM, are known (Table 4), we can calculate the wind effect coefficients Cw for each of the experiments (Table 6).
This dimensionless wind effect coefficient expresses the relationship of the field of velocities measured on a reference level (the meteorological station) and on another level close to the greenhouse vents [2].
Equations (1)–(3) have been used to obtain Cw in accordance with the three models of ventilation, M1, M2, and M3, respectively, derived from Bernoulli’s equation [4,7,10]. Ventilation models M1 and M2 consider that the wind and thermal effects complement each other. However, this does not always occur, and Li and Delsante [22] considered two situations: “fully assisting and fully opposing”. According to these authors, the latter can occur when there are only two ventilation openings. Although the greenhouse in the present study has three ventilation openings, the ventilation models present a better fit considering the fully opposing scenario (subtracting the wind and thermal effect) than considering the fully assisting scenario (adding the wind and thermal effect). It is necessary to modify models M1 and M2 so as to indicate the relationship between the two effects:
G M 1 = | G T 2 ± G w 2 | = | ( ( C d V R × S V R ) ( C d V S × S V S ) ( C d V R × S V R ) 2 + ( C d V S × S V S ) 2 ) 2 ( 2 g Δ T i o T o h S R ) ± ( C d V R × S V R + C d V S × S V S 2 ) 2 C w u o 2 |
G M 2 = | G T ± G w | = | ( C d V R × S V R ) ( C d V S × S V S ) ( C d V R × S V R ) 2 + ( C d V S × S V S ) 2 2 g Δ T i o T o h S R ± C d V R × S V R + C d V S × S V S 2 C w u o |
Equation (13) is similar to that proposed by Li and Delsante [22]. Table 6 shows the wind coefficients obtained with the original models M1 and M2 (Equations (1) and (2)) and with the modified models (Equations (13) and (14)). When we do not consider the appropriate relationship between the wind and thermal effects, very low wind effect coefficients are obtained, and even negative coefficients are found in some of the experiments for model M1 (Table 6).
To optimise the proposed ventilation models, it would be necessary to evaluate the real relationship between the two effects in each specific situation. The two simplest scenarios are those studied by Li and Delsante [22], “fully assisting and fully opposing”, and we understand that the latter occurs in our greenhouse. In commercial greenhouses with more than three ventilation surfaces, an intermediate situation can occur, making it difficult to determine an appropriate model.
In Espinoza et al. [24], the flow patterns were determined in the east sector of the experimental greenhouse (Figure 1), opening 2 lateral vents combined with 2 or 3 roof vents (making a total of 4 or 5 vents). In the windward roof vent, opposition of the wind effect and thermal effect was observed, but in the leeward roof vent, these two effects were added. By combining windward and leeward roof vents in the greenhouse, an intermittent situation between “fully assisting and fully opposing” would take place.
Table 6 and Figure 7 illustrate how the coefficient Cw obtained by model M3 (which ignores the thermal effect) maintains relatively similar values for the different tests (between 0.134 and 0.218). This means that a mean value of Cw can be considered for model M3 to predict the ventilation rate according to wind velocity, as has been done in previous research works [10,12,13,15,17,18,19]. However, the same does not apply to models M1 and M2, in which the thermal effect is taken into account in the equations (and with a wind effect opposition). For model M1, the coefficient Cw, obtained by applying Equation (13), takes values of between 0.111 and 0.433; while for model M2, applying Equation (14), it takes values of between 0.184 and 0.824 (Figure 7). This variation in Cw when applying models M1 and M2 was also observed in an Almería-type greenhouse with two side vents and two roll-up roof vents [38]. The values of Cw for models M1 and M2 have been seen to vary according to the ratio uo/Tio0.5 (Figure 7).
Low values of the ratio uo/Tio0.5 would indicate that the thermal effect has a major bearing on natural ventilation of the greenhouse. In this case, since the thermal and wind effects are opposing, on applying models M1 and M2, higher values of Cw are obtained, possibly to compensate for the strong opposition of the two effects.
High values of the ratio uo/Tio0.5 would indicate that the thermal effect has less bearing on natural ventilation and that the wind effect predominates. In this case, by applying models M1 and M2, the Cw values obtained are close to those obtained by model M3 (Figure 7), which ignores the thermal effect. It has been observed that the Cw coefficient that corresponds to each test can be calculated according to the ratio uo/Tio0.5 or δ [Cw,M1 = exp(−2.693 + 1.160/δ) (R2 = 0.94 and p-value = 0.0000); Cw,M2 = exp(−2.128+1.264/δ) (R2 = 0.98 and p-value = 0.0000)] (Figure 7). Kittas et al. [2] observed how this coefficient increased at low air velocities in a tunnel greenhouse with continuous side openings.
Although models M1, M2, and M3 can be used in their current form to predict the greenhouse ventilation rate, as is explained below, more work is required to ascertain the dependence of Cw on the ratio uo/Tio0.5 for models M1 and M2.
Using the original models M1 and M2, the mean values of Cw are 0.013 ± 0.079 (M1) and 0.014 ± 0.018 (M2) (Table 6), which are much lower than those obtained by other authors (Table 7). On applying the modified ventilation models, these values are 0.186 ± 0.098 (M1) and 0.358 ± 0.192 (M2), which are closer to the values found by other researchers in multi-span type greenhouses and are considerably higher than those found by Molina-Aiz et al. [19] in an Almería-type greenhouse (Table 7).
The use of sonic anemometry has allowed us to observe the negative interaction between wind and thermal effects in the experimental greenhouse. Other techniques, such as the tracer gas method, allow correct quantification of the ventilation rates in the greenhouse, but they provide no information on the characteristics of the airflow [1], nor do they allow us to determine the way in which the two effects interact. Had the tracer gas method been used, applying the original ventilation methods would give much lower wind effect coefficients than those obtained taking into account the correct relationship between the wind and thermal effects. The wind effect coefficients obtained by other authors using the tracer gas method (Table 7) may be affected by a negative interaction between the two effects.
Since M3 does not consider the thermal effect, in general, lower wind effect coefficients are obtained (Table 6) than with models M1 and M2. Model M3 includes the thermal effect in the wind effect coefficients, and as the relationship between both effects is negative, lower wind effect coefficients are obtained than when the two effects are treated separately. Model M3 assumes that the volumetric caudal of air in a greenhouse is proportional to: the surface area of the vents SV, the air velocity uo, and the coefficient of effectiveness of the openings Ev = Cd·Cw0.5. The latter is frequently used in the literature to characterise the wind effect in the ventilation of greenhouses [19]. The present study has used a different expression in model M3 to that used by Molina-Aiz et al. [19]; rather than the coefficient Cd and the total vent surface area SV, this study opts for the coefficients corresponding to the side vents CdVS, to the roof vent CdWR, and their corresponding surfaces. In our case, the coefficient Ev of the vents has been obtained by calculating Cd as the weighted average value between CdVS and CdWR, according to the surface of each vent.
The coefficients of effectiveness obtained in the present work with M3 (Table 6), in some tests, are lower than those obtained by other authors in greenhouses with side and roof vents, and with insect-proof screens (Table 7), possibly due to the negative interaction between the wind and thermal effects in our greenhouse.

3.3.3. Fitting the Semi-Empirical Models to Experimental Data

To assess the fit of the three models with the experimental values, the flow rates were obtained for each experiment applying the three models. These ventilation flow rates have been estimated by applying the ventilation models in four ways (Figure 8): (i) first, using the mean wind coefficient Cw (0.186 for M1, 0.358 for M2 and 0.177 for M3) and the particular discharge coefficients, CdVS and CdWR, for each test (Table 5); (ii) second, using the mean wind coefficient Cw (idem) and the mean discharge coefficients CdVS and CdWR (equal to 0.121 and 0.161, respectively); (iii) third, using the mean wind coefficient Cw (idem) and the fits CdVS = 0.028·uo + 0.028 and CdWR = 0.036·uo + 0.040; (iv) fourth, also using the fits Cw,M1 = exp(−2.693 + 1.160/δ) and Cw,M2 = exp(−2.128 + 1.264/δ).
Based solely on the determining coefficient R2, in the former case, the three models show a very good fit with the experimental data (with R2 = 0.992 for M1 and M3, R2 = 0.991 for M2).
In the second case, model M3 best fits the experimental data with R2 = 0.952, followed by M1 with R2 = 0.925, and finally M2 with R2 = 0.919. Applying the models with the fits CdVS(uo) and CdWR(uo) improves the fit of models M1 and M2 (R2 = 0.986) and of M3 (R2 = 0.987). Applying the models with the fits Cw,M1(δ), Cw,M2(δ), CdVS(uo), and CdWR(uo), the values of R2 go down slightly for models M1 and M2 (R2 = 0.973 for M1 and R2 = 0.980 for M2). In Molina-Aiz et al. [19], the models that obtain a better fit to the experimental data are Model M1 (R2 = 0.984), considering the resulting pressure distribution as the sum of the pressure fields due to stack and wind effects, and the most simplified M3 (R2 = 0.985). In the case of Molina-Aiz et al. [19], the relationship between wind and thermal effects was “fully assisting”, with two side vents (one to windward and one to leeward) and two roll-up roof vents; the air entered through the side vents to exit through the roof vents.
Analysis based solely on the coefficient R2 may give the impression that the behaviour of the three models is very good, given the high R2 values obtained in all cases. On comparing the simulated and observed GM values one by one (Figure 8), the simulated values are seen to vary considerably from the observed ones in certain cases. Table 8 presents the different statistics obtained for the three models, M1, M2, and M3, and the four different applications of the coefficients Cw, CdVS, and CdWR (Applications i, ii, iii and iv in Table 8). With model M3, higher values of GM are obtained than those observed experimentally, irrespective of how the coefficients are applied (Figure 8), always obtaining values of PBIAS<0 (model overestimation bias). Models M1 and M2 present values of PBIAS<0 (model overestimation bias) for the first three applications of the coefficients (i, ii, and iii in Table 8) and values of PBIAS>0 (model underestimation bias) on applying the models using the fits Cw,M1(δ), Cw,M2(δ), CdVS(uo), and CdWR(uo). For the remaining statistics obtained, no great differences are observed between the different applications of the coefficients Cw, CdVS, and CdWR (Applications i, ii, iii and iv in Table 8). We consider the best application of the models uses the fits obtained for CdVS and CdWR as a function of uo (for models M1, M2, and M3), and the fits obtained for Cw as a function of the ratio uo/Tio0.5 (for models M1 and M2). In this way, the values of the coefficients applied are closest to the values obtained experimentally.
In summary, model M3 considers a linear relationship between the wind velocity and the ventilation volumetric flow rates GM, a statistical device which improves the quality of the fit, and which therefore must be used and interpreted with great care [39]. It has also been observed that the ventilation rate fits better a second order polynomial than a linear relationship with air velocity. We therefore recommend using the modified model M1, since it includes both the wind and thermal effects and it is the one with the greatest physical foundation, as it combines the pressure fields rather than the flow rates [4,7]. For model M1, the lowest values of RMSD and RSR are obtained when applying the model using the fits Cw,M1(δ), Cw,M2(δ), CdVS(uo), and CdWR(uo) (Application iv in Table 8), rather than Applications i, ii, and iii. Moreover, before applying the model, it is essential to determine whether the interaction between the wind and thermal effects is positive or negative.
In view of the results obtained in this work, we recommend to experimentally determine the coefficients Cd and Cw coefficients before applying the semi-empirical models M1, M2, or M3 in a greenhouse. Special care must be taken when they are not determined experimentally and when the coefficients Cd and Cw coefficients of the bibliography have to be selected. The coefficient Cd depends on the geometry of the vent and on the characteristics of the insect-proof screen installed in the vent [19]; a coefficient obtained for a vent with similar characteristics must be selected. But as noted, the coefficient Cd varies with air velocity through the vents, and this should also be taken into account. In the case of the coefficient Cw, it has been observed that in the current form of the M1 and M2 models, this coefficient varies depending on the interaction of the thermal effect and thermal effect (ratio uo/Tio0.5), which complicates the task of selecting coefficient Cw of the bibliography that adapts well to the greenhouse where you want to apply these semi-empirical models.

3.3.4. Estimating the Contribution of the Thermal Effect in the Natural Ventilation of the Greenhouse

By applying models M1 and M2, calculating the free component of the flux induced by buoyancy forces GT and the flux induced by wind forces Gw, we estimated the contribution of the thermal effect to the natural ventilation of the greenhouse (using fits Cw,M1(δ), Cw,M2(δ), CdVS(uo), and CdWR(uo)).
The contribution of the thermal effect in all the experiments (Table 9) is sufficient for it to be considered in any study of natural ventilation of greenhouses in regions with a similar climate to the Mediterranean one, in which the temperature gradients between inside and outside the greenhouse are considerable. Indeed, in test 1, with a wind velocity of 6.86 m s−1 and ratio uo/Tio0.5 = 2.98, the thermal effect would be equivalent to 36% or 27% of the wind effect applying models M1 and M2, respectively (Table 9). The limits established by other authors should be revised. For instance, the value of 2 m s−1, above which the thermal effect may be ignored, as suggested by Boulard and Baille [7], Kittas et al. [39] and Fatnassi et al. [13], does not fit the results of the present study.

3.3.5. Estimating the Reduction in the Ventilation Rate Caused by the Insect-Proof Screens

The reduction in greenhouse ventilation due to the insect-proof screens may be considered proportional to the reduction in air velocity through the vents [40]. The reduction in the ventilation rate caused by the insect-proof screens (φ = 0.39) can be estimated by applying the ventilation models using the discharge coefficients of the openings calculated without the screens (CdVS = 0.657, CdWR = 0.712); determined following the methodology outlined by Molina-Aiz et al. [19]. Applying the mean value of Cw (equal at 0.177) for model M3, and the fits Cw,M1(δ) and Cw,M2(δ) for models M1 and M2, the insect-proof screens (φ = 0.39) have been estimated to cause a mean reduction of 70% (M1), 61% (M2), and 85% (M3).
Natural ventilation is extremely important for optimal plant growth during the summer in Mediterranean countries. For most of the year, a good system of natural ventilation will allow growers to maintain suitable microclimate conditions inside the greenhouse for the crops [41]. However, Valera et al. [3] point out that the 14.4% average ventilation surface (SV/SA) in Almería’s greenhouses is well below the minimum recommended value of 30% [42,43]. Indeed, it is also considerably lower than the 25% recommended in the Andalusian Regulations on Integrated Production of Protected Horticultural Produce [44]. Von Zabeltitz [45] recommends between 18% and 25% for the Mediterranean Basin, while Kittas et al. [41] quote a value of between 15% and 30%. It should be noted that all of the above recommendations refer to greenhouses that are not equipped with insect-proof screens over the vent openings. The value of ventilation rate RM required for a temperature rise of 5 °C between outside-inside air temperatures vary from 0.02 to 0.09 m3 s–1 m–2 (15 to 65 h–1 for a greenhouse with an average height of 5 m) depending on the solar radiation and the crop transpiration [46], with an optimum value of 45-60 h–1 [47].
In this work, the experiments were carried out for a ventilation surface, SV/SA, of 0.112. By applying the ventilation models as follows: with the mean value of Cw = 0.177 for model M3; the fits Cw,M1(δ) and Cw,M2(δ) for models M1 and M2; and the fits CdVS(uo) and CdWR(uo) for the three models; it has been estimated that for the vent arrangement in the experimental greenhouse, a ventilation surface SV/SA of 0.62 (M1), 0.54 (M2), and 0.36 (M3) would be necessary to reach a value of 45 h−1 (under conditions of southwest wind and uo = 5 m s−1).

4. Conclusions

Sonic anemometry has allowed us to identify the entrance and exit vents of greenhouse ventilation air, thus allowing us to establish natural ventilation flow patterns for greenhouses. According to most of the works reviewed in the literature, the ventilation rate of a greenhouse with natural ventilation fits a linear relationship with air velocity. In the present case, it was observed that this relationship fits better a second order polynomial, RM = 0.37 uo2 + 0.03 uo + 0.75 (R2 = 0.99 and p-value = 0.0000). Opening the roof vent to windward, one side vent to leeward, and the other side vents to windward (this last vent obstructed by another greenhouse) brings about the opposition of the thermal and wind effects in natural ventilation of the greenhouse, as it drives the air entering through it downwards, opposing natural convection due to the thermal effect. The modification of the ventilation models obtained from Bernouilli’s equation in order to add or subtract the airflows due to wind and thermal effects improves their precision.
G ( M 1 ) = | G T 2 ± G w 2 |   and   G ( M 2 ) = | G T ± G w |
A linear relationship has been obtained, which allows us to estimate the discharge coefficient of the side vents (CdVS) and roof vent (CdWR) as a function of wind velocity [CdVS = 0.028·uo + 0.028 (R2 = 0.92); CdWR = 0.036·uo + 0.040 (R2 = 0.96)]. The wind effect coefficient determined by applying models M1 and M2 has not proved constant for the different tests carried out, but rather it varies depending on the ratio uo/Tio0.5 or δ [Cw,M1 = exp(−2.693 + 1.160/δ) (R2 = 0.94); Cw,M2 = exp(−2.128+1.264/δ) (R2 = 0.98]. The fits for Cd and Cw allow us to estimate the ventilation flow rate from the outside and inside temperatures and from the wind velocity using the ventilation models.
The contribution of the thermal effect to the natural ventilation of the greenhouse has been quantified based on models M1 and M2. The opposition of the thermal effect accounted for at least 36% or 27% of the wind effect (applying models M1 and M2, respectively). The thermal effect must therefore always be taken into account in studies on the natural ventilation of greenhouses. Applying models M1, M2, and M3, it has been estimated that the presence of insect-proof screens on the vents (φ = 39.0%) can reduce the ventilation flow rate of the greenhouse by 70%, 61%, and 85%, respectively. Likewise, by applying these models, it has been ascertained that the surface area of ventilation with respect to the surface area of the soil, SV/SA, should be 0.36 (M3), 0.54 (M2), and 0.62 (M1) to reach a ventilation rate of 45 h−1 for an air velocity of 5 m s−1.

Author Contributions

Conceived and research design: A.L.-M., F.D.M.-A. and D.L.V.-M.; acquisition data and statistical analysis: A.L.-M., K.E.E.-R. and J.L.-M.; analysis and interpretation of data: A.L.-M. and F.D.M.-A.; drafting the manuscript: A.L.-M.; critical revision of the manuscript for important intellectual content: F.D.M.-A., D.L.V.-M. and A.P.-F.; coordinating and supervising the research work: D.L.V.-M. and A.P.-F.; project administration: D.L.V.-M.; funding acquisition: D.L.V.-M.

Funding

This research was funded by the Spanish Ministry of Economy and Competitiveness and the European Regional Development Fund (ERDF) by means of research grant AGL2015-68050-R.

Acknowledgments

The authors wish to express their gratitude to the Research Centre CIAIMBITAL of the University of Almería (Spain) for their support throughout the development of this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MDbias
NSENash-Sutcliffe efficiency
PBIASpercent bias
RMSDroot mean squared deviation
RSRRMSD-observations standard deviation ratio
Nomenclature
Cdtotal discharge coefficient of the opening
Cwwind effect coefficient for the mean airflow
ethickness of the screen [m]
EGerror in the calculation of volumetric flow rates [%]
Evcoefficient of effectiveness of the openings
Ggravitational constant [m s−2]
Gvolumetric flow rate [m3 s−1]
GTfree component of the ventilation flux induced by buoyancy forces [m3 s−1]
Gwforced component of the ventilation flux induced by wind forces [m3 s−1]
hSRdifference in height between side and roof openings [m]
RHrelative humidity [%]
Kpscreen permeability [m2]
nnumber of measurement points in vents
Rgsolar radiation [W m−2]
RMventilation rate for greenhouse [h−1]
SAgreenhouse area [m2]
SVsurface area of the vent openings [m2]
ttime [s]
Ttemperature [°C]
uair velocity [m s−1]
Yinertial factor
Greek Letters
Δdifference
δratio for wind/thermal effect
θwind direction [°]
φporosity [%]
Subscripts
iinside
jmeasurement point
Lleeward
Maverage value for the greenhouse
M1ventilation model 1
M2ventilation model 2
M3ventilation model 3
ooutside
obsobserved
Rroof vent
Sside vent
simsimulated
Vvent
Wwindward
xlongitudinal component
ytransversal component
zvertical component
Superscripts
*corrected
fluctuating component

References

  1. Boulard, T.; Meneses, J.F.; Mermier, M.; Papadakis, G. The mechanisms involved in the natural ventilation of greenhouses. Agric. Forest. Meteorol. 1996, 79, 61–77. [Google Scholar] [CrossRef]
  2. Kittas, C.; Boulard, T.; Mermier, M.; Papadakis, G. Wind induced air exchange rates in a greenhouse tunnel with continuous side openings. J. Agric. Eng. Res. 1996, 65, 37–49. [Google Scholar] [CrossRef]
  3. Valera, D.L.; Belmonte, L.J.; Molina-Aiz, F.D.; López, A. Greenhouse Agriculture in Almería. A Comprehensive Techno-Economic Analysis, 1st ed.; Cajamar Caja Rural: Almería, Spain, 2016; p. 408. [Google Scholar]
  4. Kittas, C.; Boulard, T.; Papadakis, G. Natural ventilation of a greenhouse with ridge and side openings: Sensitivity to temperature and wind effects. Trans. ASAE 1997, 40, 415–425. [Google Scholar] [CrossRef]
  5. Morris, L.G.; Neale, F.E. The Infra-Red Carbon Dioxide Gas Analyser and Its use in Glasshouse Research; Tech. Memo. No. 99; National Institute of Agricultural Engineering: Silsoe, UK, 1954; p. 13. [Google Scholar]
  6. Businger, J.A. The glasshouse (greenhouse) climate. In Physics of Plant Environmen, 1st ed.; van Wijk, W.R., Ed.; North-Holland Publishing Co.: Amsterdam, The Netherlands, 1963; pp. 277–318. [Google Scholar]
  7. Boulard, T.; Baille, A. Modelling of air exchange rate in a greenhouse equipped with continuous roof vents. J. Agric. Eng. Res. 1995, 61, 37–48. [Google Scholar] [CrossRef]
  8. Bruce, J.M. Natural convection through openings and its applications to cattle building ventilation. J. Agric. Eng. Res. 1978, 23, 151–167. [Google Scholar] [CrossRef]
  9. Bruce, J.M. Ventilation of a model livestock building by thermal buoyancy. Trans. ASAE 1982, 25, 1724–1726. [Google Scholar] [CrossRef]
  10. Boulard, T.; Feuilloley, P.; Kittas, C. Natural ventilation performance of six greenhouse and tunnel types. J. Agric. Eng. Res. 1997, 67, 249–266. [Google Scholar] [CrossRef]
  11. Papadakis, G.; Mermier, M.; Meneses, J.F.; Boulard, T. Measurement and analysis of air exchange rates in a greenhouse with continuous roof and side openings. J. Agric. Eng. Res. 1996, 63, 219–228. [Google Scholar] [CrossRef]
  12. Demrati, H.; Boulard, T.; Bekkaoui, A.; Bouirden, L. Natural ventilation and microclimatic performance of a large-scale banana greenhouse. J. Agric. Eng. Res. 2001, 80, 261–271. [Google Scholar] [CrossRef]
  13. Fatnassi, H.; Boulard, T.; Demrati, H.; Bouirden, L.; Sappe, G. SE-Structures and Environment: Ventilation Performance of a Large Canarian-Type Greenhouse equipped with Insect-proof Nets. Biosyst. Eng. 2002, 82, 97–105. [Google Scholar] [CrossRef]
  14. Sase, S.; Reiss, E.; Both, A.J.; Roberts, W.J. A natural ventilation model for open-roof greenhouses; Paper No. 024010. In Proceedings of the ASAE International Meeting/CIGR XVth World Congress, ASAE, Chicago, IL, USA, 28–31 July 2002. [Google Scholar]
  15. Fatnassi, H.; Boulard, T.; Lagier, J. Simple indirect estimation of ventilation and crop transpiration rates in a greenhouse. Biosyst. Eng. 2004, 88, 467–478. [Google Scholar] [CrossRef]
  16. Abreu, P.E.; Boulard, T.; Mermier, M.; Meneses, J.F. Parameter estimation and selection of a greenhouse natural ventilation model and its use on an energy balance model to estimate the greenhouse air temperature. Acta Hortic. 2006, 691, 611–618. [Google Scholar] [CrossRef]
  17. Katsoulas, N.; Bartzanas, T.; Boulard, T.; Mermier, M.; Kittas, C. Effect of vent openings and insect screens on greenhouse ventilation. Biosyst. Eng. 2006, 93, 427–436. [Google Scholar] [CrossRef]
  18. Teitel, M.; Liran, O.; Tanny, J.; Barak, M. Wind driven ventilation of a mono-span greenhouse with a rose crop and continuous screened side vents and its effect on flow patterns and microclimate. Biosyst. Eng. 2008, 101, 111–122. [Google Scholar] [CrossRef]
  19. Molina-Aiz, F.D.; Valera, D.L.; Peña, A.A.; Gil, J.A.; López, A. A Study of Natural Ventilation in an Almería-Type Greenhouse with Insect Screens by Means of Tri-sonic Anemometry. Biosyst. Eng. 2009, 104, 224–242. [Google Scholar] [CrossRef]
  20. Sase, S.; Takakura, T.; Nara, M. Wind tunnel testing on airflow and temperature distribution of a naturally ventilated greenhouse. Acta Hortic. 1984, 148, 329–336. [Google Scholar] [CrossRef]
  21. Bot, G.P.A. Greenhouse Climate: From Physical Processes to a Dynamic Model. Ph.D. Thesis, Agricultural University, Wageningen, The Netherlands, 1983. [Google Scholar]
  22. Li, Y.; Delsante, A. Natural ventilation induced by combined wind and thermal forces. Build. Environ. 2001, 36, 59–71. [Google Scholar] [CrossRef]
  23. Van Buggenhout, S.; Van Brecht, A.; Eren Özcan, S.; Vranken, E.; Van Malcot, W.; Berckmans, D. Influence of sampling positions on accuracy of tracer gas measurements in ventilated spaces. Biosyst. Eng. 2009, 104, 216–223. [Google Scholar] [CrossRef]
  24. Espinoza, K.; López, A.; Valera, D.L.; Molina-Aiz, F.D.; Torres, J.A.; Peña, A. Effects of ventilator configuration on the flow pattern of a naturally-ventilated three-span Mediterranean greenhouse. Biosyst. Eng. 2017, 164, 13–30. [Google Scholar] [CrossRef]
  25. Boulard, T.; Kittas, C.; Papadakis, G.; Mermier, M. Pressure Field and Airflow at the Opening of a Naturally Ventilated Greenhouse. J. Agric. Eng. Res. 1998, 71, 93–102. [Google Scholar] [CrossRef]
  26. Teitel, T.; Tanny, J.; Ben-Yakir, D.; Barak, M. Airflow patterns through roof openings of a naturally ventilated greenhouse and their effect on insect penetration. Biosyst. Eng. 2005, 92, 341–353. [Google Scholar] [CrossRef]
  27. Shilo, E.; Teitel, M.; Mahrer, Y.; Boulard, T. Air-flow patterns and heat fluxes in roof-ventilated multi-span greenhouse with insect-proof screens. Agric. Forest. Meteorol. 2004, 122, 3–20. [Google Scholar] [CrossRef]
  28. Boulard, T.; Papadakis, G.; Kittas, C.; Mermier, M. Air flow and associated sensible heat exchanges in a naturally ventilated greenhouse. Agric. Forest. Meteorol. 1997, 88, 111–119. [Google Scholar] [CrossRef]
  29. Valera, D.L.; Álvarez, A.J.; Molina-Aiz, F.D. Aerodynamic analysis of several insect-proof screens used in greenhouses. Span. J. Agris. Res. 2006, 4, 273–279. [Google Scholar] [CrossRef]
  30. López, A.; Molina-Aiz, F.D.; Valera, D.L.; Peña, A. Wind tunnel analysis of the airflow through insect-proof screens and comparison of their effect when installed in a Mediterranean greenhouse. Sensors 2016, 16, 690. [Google Scholar] [CrossRef] [PubMed]
  31. Kobayashi, K.; Salam, M.U. Comparing simulated and measured values using mean squared deviation and its components. Agron. J. 2000, 92, 345–352. [Google Scholar] [CrossRef]
  32. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models: Part 1. A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  33. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  34. Gupta, H.V.; Sorooshian, S.; Yapo, P.O. Status of automatic calibration for hydrologic models: Comparison with multilevel expert calibration. J. Hydrol. Eng. 1999, 4, 135–143. [Google Scholar] [CrossRef]
  35. Singh, J.; Knapp, H.V.; Demissie, M. Hydrologic Modeling of the Iroquois River watershed using HSPF and SWAT. In ISWS CR 2004-08; Illinois State Water Survey: Champaign, IL, USA, 2004. [Google Scholar]
  36. Legates, D.R.; McCabe, G.J. Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  37. López, A.; Valera, D.L.; Molina-Aiz, F.D.; Peña, A. Sonic anemometry measurements to determine airflow patterns in multi-tunnel greenhouse. Span. J. Agris. Res. 2012, 10, 631–642. [Google Scholar] [CrossRef]
  38. Molina-Aiz, F.D.; Valera, D.L.; López, A. Numerical and experimental study of heat and mass transfer in Almería-type greenhouse. Acta Hortic. 2015, 1170, 209–218. [Google Scholar] [CrossRef]
  39. Kittas, C.; Draoui, B.; Boulard, T. Quantifcation du taux d’aération d’une serre à ouvrant continu en toiture. Agric. Forest. Meteorol. 1995, 77, 95–111. [Google Scholar] [CrossRef]
  40. Teitel, M.; Barak, M.; Berlinger, M.J.; Lebiush-Mordechai, S. Insect-proof screens in greenhouses: Their effect on roof ventilation and insect penetration. Acta Hortic. 1999, 507, 25–34. [Google Scholar] [CrossRef]
  41. Kittas, C.; Katsoulas, N.; Bartzanas, T.; Bakker, S. Greenhouse climate control and energy use. In Good Agricultural Practices for Greenhouse Vegetable Crops. Principles for Mediterranean Climate Areas; Baudoin, W., Nono-Womdim, R., Lutaladio, N., Hodder, A., Castilla, N., Leonardi, C., De Pascale, S., Qaryouti, M., Duffy, R., Eds.; Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 2013; pp. 63–96. [Google Scholar]
  42. Verheye, P.; Verlodt, H. Comparison of different systems for static ventilation of hemispheric plastic greenhouses. Acta Hortic. 1990, 281, 183–197. [Google Scholar] [CrossRef]
  43. Castilla, N.P. Invernaderos de Plástico. Tecnología y Manejo, 1st ed.; Mundi-Prensa: Madrid, Spain, 2007; p. 462. (In Spanish) [Google Scholar]
  44. Junta, D.A. Reglamento Específico de Producción Integrada de Cultivos Hortícolas Protegidos; Boletín Oficial de la Junta de Andalucía (BOJA) nº 211 de 25 de Octubre de 2007: Sevilla, Spain, 2007. (In Sanish) [Google Scholar]
  45. Von Zabeltitz, C. L’efficacité énergétique dans la conception des serres méditerranéennes. Plasticulture 1993, 96, 6–16. [Google Scholar]
  46. ASABE Standards. EP406.4: Heating, Ventilating and Cooling Greenhouses; ASAE: St. Joseph, MI, USA, 2008. [Google Scholar]
  47. Hellickson, M.A.; Walker, J.N. Ventilation of Agricultural Structures, ASAE Monograph No. 6, 1st ed.; ASABE: St. Joseph, MI, USA, 1983. [Google Scholar]
Figure 1. Location of the experimental greenhouse at the farm.
Figure 1. Location of the experimental greenhouse at the farm.
Agronomy 09 00736 g001
Figure 2. Measurement points at the lateral vents and at the roof vent (view for south side). Measurement tests 1 to 4 (a) and 5 to 10 (b). Separation of the anemometers from the vents (c).
Figure 2. Measurement points at the lateral vents and at the roof vent (view for south side). Measurement tests 1 to 4 (a) and 5 to 10 (b). Separation of the anemometers from the vents (c).
Agronomy 09 00736 g002
Figure 3. Polar histograms of the horizontal angle and projection of the air velocity (plane XY) at the side vents (middle of the vent) and at the roof vent in the western sector of the experimental greenhouse in conditions of light southwest wind (SW). Tests 4 on 07/05/2009 (a) and 6 on 23/04/2009 (b).
Figure 3. Polar histograms of the horizontal angle and projection of the air velocity (plane XY) at the side vents (middle of the vent) and at the roof vent in the western sector of the experimental greenhouse in conditions of light southwest wind (SW). Tests 4 on 07/05/2009 (a) and 6 on 23/04/2009 (b).
Agronomy 09 00736 g003
Figure 4. Polar histograms of the horizontal angle and projection of the air velocity (plane XY) at the side vents (middle of the vent) and at the roof vent of the western sector of the experimental greenhouse in conditions of moderate southwest wind (SW). Test 2 on 08/04/2009.
Figure 4. Polar histograms of the horizontal angle and projection of the air velocity (plane XY) at the side vents (middle of the vent) and at the roof vent of the western sector of the experimental greenhouse in conditions of moderate southwest wind (SW). Test 2 on 08/04/2009.
Agronomy 09 00736 g004
Figure 5. Ventilation rate for greenhouse as a function of air velocity (a). Error of closure of the flows EG in function of the ratio uo/Tio0.5 (b). Subfigure (a): second order polynomial fit (); linear fit (---).
Figure 5. Ventilation rate for greenhouse as a function of air velocity (a). Error of closure of the flows EG in function of the ratio uo/Tio0.5 (b). Subfigure (a): second order polynomial fit (); linear fit (---).
Agronomy 09 00736 g005
Figure 6. Discharge coefficients of side vents CdVS (◊) and roof vent CdWR (□), obtained for each test, expressed as a function of air velocity uo.
Figure 6. Discharge coefficients of side vents CdVS (◊) and roof vent CdWR (□), obtained for each test, expressed as a function of air velocity uo.
Agronomy 09 00736 g006
Figure 7. Greenhouse wind effect coefficient Cw, obtained for each test, expressed as a function of the ratio uo/Tio0.5 or δ (b): , Cw-M1; □, Cw-M2; ×, Cw-M3.
Figure 7. Greenhouse wind effect coefficient Cw, obtained for each test, expressed as a function of the ratio uo/Tio0.5 or δ (b): , Cw-M1; □, Cw-M2; ×, Cw-M3.
Agronomy 09 00736 g007
Figure 8. Values of simulated vs. observed GM. Applying the models with mean values of Cw and the particular discharge coefficients CdVS and CdWR for each test (a), with mean values of Cw, CdVS and CdWR (b); with mean values of Cw and the fits CdVS(uo) and CdWR(uo) (c); with fits Cw,M1(δ), Cw,M2(δ), CdVS(uo) and CdWR(uo) (d). , model M1; □, model M2; ×, model M3.
Figure 8. Values of simulated vs. observed GM. Applying the models with mean values of Cw and the particular discharge coefficients CdVS and CdWR for each test (a), with mean values of Cw, CdVS and CdWR (b); with mean values of Cw and the fits CdVS(uo) and CdWR(uo) (c); with fits Cw,M1(δ), Cw,M2(δ), CdVS(uo) and CdWR(uo) (d). , model M1; □, model M2; ×, model M3.
Agronomy 09 00736 g008
Table 1. Outside climatic conditions for the measurement tests. Average wind speed uo [m s−1], wind direction θ [°], outside and inside temperature To and Ti [°C], outside and inside humidity RHo and RHi [%], outside radiation Rg [W m−2] and ratio determining the relative importance of the wind and thermal buoyancy forces uo/Tio0.5.
Table 1. Outside climatic conditions for the measurement tests. Average wind speed uo [m s−1], wind direction θ [°], outside and inside temperature To and Ti [°C], outside and inside humidity RHo and RHi [%], outside radiation Rg [W m−2] and ratio determining the relative importance of the wind and thermal buoyancy forces uo/Tio0.5.
TestDateTimeuoθaRHoRHiToTiRguo/Tio0.5
107/04/200911:52–14:476.86 ± 1.41300 ± 767 ± 252 ± 717.5 ± 0.422.8 ± 1.9527 ± 2592.98
208/04/200910:49–13:314.16 ± 1.06295 ± 1529 ± 838 ± 218.3 ± 0.425.7 ± 1.4562 ± 2371.53
314/04/200911:21–14:014.01 ± 1.13294 ± 1072 ± 257 ± 216.3 ± 0.623.7 ± 0.7692 ± 1081.47
407/05/200910:56–12:362.03 ± 0.84264 ± 1436 ± 360 ± 522.8 ± 0.727.6 ± 0.6784 ± 600.93
517/04/200911:06–13:071.94 ± 0.70226 ± 2559 ± 553 ± 416.9 ± 0.425.7 ± 1.0584 ± 1750.65
623/04/200911:23–13:162.34 ± 0.98267 ± 1438 ± 255 ± 622.1 ± 0.628.3 ± 0.1809 ± 560.94
722/06/200911:14–12:583.42 ± 0.50258 ± 856 ± 470 ± 225.5 ± 0.427.2 ± 0.5617 ± 712.57
826/06/200911:17–13:042.67 ± 0.72227 ± 2064 ± 273 ± 324.0 ± 0.528.2 ± 0.2751 ± 1021.29
902/07/200911:00–12:452.62 ± 0.63239 ± 1765 ± 376 ± 327.0 ± 0.930.8 ± 0.4725 ± 681.34
1002/07/200914:41–16:283.22 ± 0.48242 ± 1460 ± 269 ± 228.0 ± 0.832.1 ± 0.5868 ± 371.58
a Direction perpendicular to the roof window at windward is 208° for a southwest wind (SW).
Table 2. Sensors used and their specifications.
Table 2. Sensors used and their specifications.
TypeModelManufacturerMeasurementsMeasurement RangeResolutionAccuracy
3D sonic AnemometerCSAT3Campbell Scientific Spain S.LAir velocity (ux, uy and uz) and air temperature0–30 m s−1
−30–50 °C
0.001 m s−1
0.002 °C
±0.001 m s−1
±0.002 °C
2D sonic AnemometerWindsonicGill Instrument LTDAir velocity (ux and uy)0–60 m s−1
0–359°
0.01 m s−1
2%
Temperature Pt1000 and capacitive humidity sensorBUTRON IIHortimax S.L.Air temperature and humidity−25–75 °C
0–100%
0.01 °C
1%
±0.01 °C
±3%
Cup anemomer and vaneMeteostation IIHortimax S.L.Wind velocity
Wind direction
0–40 m s−1
0–359°
0.1 m s−1
±5%
±5%
Radiation sensorKipp SolariHortimax S.L.Solar radiation0–2000 W m−20.1 W m−2 ±20 W m−2
Termistor and capacitive humidity sensorHOBO Pro Temp-RH U23-001Onset Computer Corp.Air temperature and humidity−40–75 °C
0–100%
0.02 °C
0.03%
±0.18 °C
±3%
Table 3. Average values of the ux component of the air velocity perpendicular to the greenhouse vents and the corrected ux*. LS, leeward side vent; WS, windward side vent; WR, windward roof vent.
Table 3. Average values of the ux component of the air velocity perpendicular to the greenhouse vents and the corrected ux*. LS, leeward side vent; WS, windward side vent; WR, windward roof vent.
Test Numberux, [m s−1]ux*, [m s−1]
LSWSWRLSWS
1−0.42 ± 0.13−0.29 ± 0.110.83 ± 0.29−0.42−0.29
2−0.14 ± 0.06−0.14 ± 0.060.30 ± 0.09−0.14−0.15
3−0.17 ± 0.06−0.16 ± 0.070.27 ± 0.04−0.18−0.17
40.03 ± 0.04−0.11 ± 0.050.07 ± 0.020.03−0.11
5−0.02 ± 0.06−0.08 ± 0.060.09 ± 0.03−0.01−0.08
6−0.01 ± 0.06−0.12 ± 0.080.11 ± 0.020.02−0.11
7−0.09 ± 0.07−0.17 ± 0.050.25 ± 0.06−0.08−0.17
80.01 ± 0.08−0.15 ± 0.050.12 ± 0.050.03−0.14
90.00 ± 0.06−0.11 ± 0.040.11 ± 0.05−0.01−0.11
10−0.06 ± 0.13−0.14 ± 0.030.17 ± 0.05−0.06−0.14
Positive values for an inflow and negative values for an outflow.
Table 4. Average values G (Equation (8)) and turbulent component G’ (Equation (9)) of the ventilation volumetric flow rates through each vent opening. Error in the calculation of ventilation volumetric flow rates EG and average values of ventilation rate R. LS, leeward side vent; WS, windward side vent; WR, windward roof vent; M, average values for the greenhouse.
Table 4. Average values G (Equation (8)) and turbulent component G’ (Equation (9)) of the ventilation volumetric flow rates through each vent opening. Error in the calculation of ventilation volumetric flow rates EG and average values of ventilation rate R. LS, leeward side vent; WS, windward side vent; WR, windward roof vent; M, average values for the greenhouse.
TestGLS
[m3 s−1]
GWS
[m3 s−1]
GWR
[m3 s−1]
GM
[m3 s−1]
EG, [%]G’LS
[m3 s−1]
G’WS
[m3 s−1]
G’WR
[m3 s−1]
RM
[h−1]
1−7.8−5.314.113.67.54.14.08.218.2
2−2.6−2.74.04.6−27.82.21.94.06.2
3−3.3−3.24.65.5−33.62.62.34.07.4
40.6−2.01.21.9−11.41.71.22.62.6
5−0.1−1.51.591.61.82.11.83.02.1
60.3−2.01.942.112.31.61.32.92.8
7−1.5−3.14.304.4−6.21.61.33.36.0
80.5−2.62.032.6−5.31.81.23.33.5
9−0.2−2.01.802.0−20.81.41.03.02.7
10−1.2−2.62.933.3−24.31.91.33.44.5
Positive values for an inflow and negative values for an outflow.
Table 5. Discharge coefficients obtained at the vents of the experimental greenhouse. CdLS, leeward vent; CdWS, windward vent; CdVS, average coefficient for side vents; CdWR, windward roof vent.
Table 5. Discharge coefficients obtained at the vents of the experimental greenhouse. CdLS, leeward vent; CdWS, windward vent; CdVS, average coefficient for side vents; CdWR, windward roof vent.
TestCdLSCdWSCdVSCdWR
10.2270.1960.2110.288
20.1440.1470.1450.181
30.1610.1580.1600.193
40.0730.1270.1000.103
50.0300.1110.0710.120
60.0500.1260.0880.130
70.1100.1530.1310.183
80.0630.1430.1030.132
90.0400.1260.0840.124
100.0970.1400.1190.154
Average0.100
± 0.062
0.143
± 0.023
0.121
± 0.042
0.161
± 0.054
Table 6. Wind effect coefficients Cw obtained for the experimental greenhouse according to the original ventilation models M1, M2, and M3 (Equations (1), (2) and (3)) and the modified models M1 and M2 (Equations (13) and (14)). Ev, coefficient of the effectiveness of the opening Ev = CdM Cw0.5 based on model M3.
Table 6. Wind effect coefficients Cw obtained for the experimental greenhouse according to the original ventilation models M1, M2, and M3 (Equations (1), (2) and (3)) and the modified models M1 and M2 (Equations (13) and (14)). Ev, coefficient of the effectiveness of the opening Ev = CdM Cw0.5 based on model M3.
TestCwEv
OriginalModifiedM3M3
M1M2M1M2
10.0850.0390.1110.1840.2180.110
20.0230.0020.1170.2330.1340.087
30.0410.0060.1410.2770.1710.052
40.0150.0000.2280.4550.1730.064
5−0.1880.0430.4330.8240.1980.038
6−0.0300.0020.2480.4940.1840.044
70.0890.0370.1240.2110.2040.101
80.0350.0040.1680.3330.1730.035
90.0210.0010.1550.3090.1580.058
100.0440.0070.1330.2580.1610.079
Average0.013
± 0.079
0.014
± 0.018
0.186
± 0.098
0.358
± 0.192
0.177
± 0.025
0.060
± 0.023
Table 7. Values of the wind effect coefficient Cw, discharge coefficient of the openings Cd, coefficient of effectiveness of the openings Ev = Cd Cw0.5 obtained by different authors in different types of greenhouses with side and roof vents, with and without insect-proof screens (all values are taken from Molina-Aiz et al. [19]). SA, greenhouse base surface area [m2]; uo, air velocity [m s−1]; M, model used.
Table 7. Values of the wind effect coefficient Cw, discharge coefficient of the openings Cd, coefficient of effectiveness of the openings Ev = Cd Cw0.5 obtained by different authors in different types of greenhouses with side and roof vents, with and without insect-proof screens (all values are taken from Molina-Aiz et al. [19]). SA, greenhouse base surface area [m2]; uo, air velocity [m s−1]; M, model used.
MCwEvCdGreenhouse TypeMethod *SAuoSource
Without Insect-Proof Screens
M10.0790.2100.754Multi-spanGas4160–8[11]
M10.090.2040.68Multi-spanGas4160–8[4]
M10.0980.1750.56Multi-spanGas2421–3[14]
M10.103 a0.210-Multi-spanGas1790.5–2.7[16]
M10.1160.3030.89Multi-spanGas2421–3[14]
M10.120.2080.6TunnelGas4162[2]
M10.130.2520.7Multi-spanGas4160–8[4]
M20.10.1780.42Multi-spanGas4160–8[4]
M30.012 a0.07–0.18-TunnelGas2401–2[10]
M30.04 a0.131-TunnelOm3685[15]
M30.085 a0.190-Multi-tunnelGas4165.3[10]
M30.173 a0.270-Multi-spanGas100002–4[12]
MφCwEvCdGreenhouse TypeMethod *SAuoSource
With Insect-Proof Screens
M10.340.0480.0430.194Almería3D16945.71[19]
M20.340.0210.0280.194Almería16945.71[19]
M30.340.0660.0500.194Almería16945.71[19]
M30.50.022 a0.096-Single tunnelGas1602.2[17]
M30.350.0710.0690.253Single tunnelGas/H2O/3D74.44.5[18]
M30.690.1100.1400.42CanaryGas56001-2[13]
φ, porosity of the insect-proof screens; SA, greenhouse base surface area; ue, air velocity. a, wind effect coefficient Cw calculated using Cd = 0.65 (the discharge coefficient of the openings was not indicated by these authors). * Method: Gas (N2O tracer gas), H2O (Mass balance on water vapour), 3D (trisonic anemometer), Om (omnidirectional hot-ball anemometer).
Table 8. Statistical parameters obtained for the average volumetric flow rate for the greenhouse GM,sim (simulated with the different models M1, M2 and M3) with the experimentally observed values GM,obs. R2, determining coefficient; RMSD, root mean squared error; MD, bias; NSE, the Nash-Sutcliffe efficiency; PBIAS, percent bias; RSR, RMSD-observations standard deviation ratio. Applying the models with mean values of Cw and the particular discharge coefficients CdVS and CdWR for each test (Application i), with mean values of Cw, CdVS and CdWR (Application ii); with mean values of Cw and the fits CdVS(uo) and CdWR(uo) (Application iii); with fits Cw,M1(δ), Cw,M2(δ), CdVS(uo) and CdWR(uo) (Application iv).
Table 8. Statistical parameters obtained for the average volumetric flow rate for the greenhouse GM,sim (simulated with the different models M1, M2 and M3) with the experimentally observed values GM,obs. R2, determining coefficient; RMSD, root mean squared error; MD, bias; NSE, the Nash-Sutcliffe efficiency; PBIAS, percent bias; RSR, RMSD-observations standard deviation ratio. Applying the models with mean values of Cw and the particular discharge coefficients CdVS and CdWR for each test (Application i), with mean values of Cw, CdVS and CdWR (Application ii); with mean values of Cw and the fits CdVS(uo) and CdWR(uo) (Application iii); with fits Cw,M1(δ), Cw,M2(δ), CdVS(uo) and CdWR(uo) (Application iv).
Application iApplication ii
M1M2M3M1M2M3
R20.9920.9910.9920.9250.9190.952
RMSD1.72.62.01.31.01.7
RMSD (%)406248312440
MD−0.9−1.2−1.5−0.1−0.2−0.9
NSE0.80.40.70.90.90.8
PBIAS−22.0−28.0−36.9−2.7−5.1−21.0
RSR1.62.51.91.21.01.6
Application iiiApplication iv
M1M2M3M1M2M3
R20.9860.9860.9870.9730.9800.987
RMSD1.82.72.10.91.32.1
RMSD (%)426550233050
MD−0.9−1.2−1.50.81.1−1.5
NSE0.70.40.60.90.90.6
PBIAS−20.7−27.7−36.918.425.9−36.9
RSR1.72.62.00.91.22.0
Table 9. Ventilation flow rates generated by the wind effect (Gw) and the thermal effect (GT) estimated by applying ventilation models M1 and M2 using the fits CwM2(δ), CdVS(uo), and CdWR(uo).
Table 9. Ventilation flow rates generated by the wind effect (Gw) and the thermal effect (GT) estimated by applying ventilation models M1 and M2 using the fits CwM2(δ), CdVS(uo), and CdWR(uo).
Model M1Model M2
TestuoGT *GwGT/GwGT *GwGT/Gw
16.86−5.114.10.36−5.119.00.27
24.16−4.05.60.70−4.07.60.52
34.01−3.95.30.73−3.97.10.54
42.03−1.91.61.15−1.92.90.85
51.94−2.51.51.66−2.52.01.23
62.34−2.42.11.14−2.42.80.85
73.42−1.64.00.41−1.65.40.31
82.67−2.12.60.83−2.13.50.61
92.62−2.02.50.79−2.03.40.59
103.22−2.43.60.67−2.44.80.50
Average0.85 Average0.63
* negative values indicate that the flow generated by the thermal effect opposes that which is generated by the wind effect.

Share and Cite

MDPI and ACS Style

López-Martínez, A.; Molina-Aiz, F.D.; Valera-Martínez, D.L.; López-Martínez, J.; Peña-Fernández, A.; Espinoza-Ramos, K.E. Application of Semi-Empirical Ventilation Models in A Mediterranean Greenhouse with Opposing Thermal and Wind Effects. Use of Non-Constant Cd (Pressure Drop Coefficient Through the Vents) and Cw (Wind Effect Coefficient). Agronomy 2019, 9, 736. https://doi.org/10.3390/agronomy9110736

AMA Style

López-Martínez A, Molina-Aiz FD, Valera-Martínez DL, López-Martínez J, Peña-Fernández A, Espinoza-Ramos KE. Application of Semi-Empirical Ventilation Models in A Mediterranean Greenhouse with Opposing Thermal and Wind Effects. Use of Non-Constant Cd (Pressure Drop Coefficient Through the Vents) and Cw (Wind Effect Coefficient). Agronomy. 2019; 9(11):736. https://doi.org/10.3390/agronomy9110736

Chicago/Turabian Style

López-Martínez, Alejandro, Francisco D. Molina-Aiz, Diego L. Valera-Martínez, Javier López-Martínez, Araceli Peña-Fernández, and Karlos E. Espinoza-Ramos. 2019. "Application of Semi-Empirical Ventilation Models in A Mediterranean Greenhouse with Opposing Thermal and Wind Effects. Use of Non-Constant Cd (Pressure Drop Coefficient Through the Vents) and Cw (Wind Effect Coefficient)" Agronomy 9, no. 11: 736. https://doi.org/10.3390/agronomy9110736

APA Style

López-Martínez, A., Molina-Aiz, F. D., Valera-Martínez, D. L., López-Martínez, J., Peña-Fernández, A., & Espinoza-Ramos, K. E. (2019). Application of Semi-Empirical Ventilation Models in A Mediterranean Greenhouse with Opposing Thermal and Wind Effects. Use of Non-Constant Cd (Pressure Drop Coefficient Through the Vents) and Cw (Wind Effect Coefficient). Agronomy, 9(11), 736. https://doi.org/10.3390/agronomy9110736

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

Article Metrics

Back to TopTop