Modeling and Simulation of Temperature and Relative Humidity Inside a Growth Chamber

: Modeling and simulation of internal variables such as temperature and relative humidity are relevant for designing future climate control systems. In this paper, a mathematical model is proposed to predict the internal variables temperature and relative humidity (RH) of a growth chamber (GCH). Both variables are incorporated in a set of ﬁrst-order di ﬀ erential equations, considering an energy-mass balance. The results of the model are compared and assessed in terms of the coe ﬃ cients of determination (R 2 ) and the root mean squared error (RMSE). The R 2 and RMSE computed were R 2 = 0.96, R 2 = 0.94, RMSE = 0.98 ◦ C, and RMSE = 1.08 ◦ C, respectively, for the temperature during two consecutive weeks; and R 2 = 0.83, R 2 = 0.81, RMSE = 5.45%RH, and RMSE = 5.48%RH, respectively, for the relative humidity during the same period. Thanks to the passive systems used to control internal conditions, the growth chamber gives average di ﬀ erences between inside and outside of + 0.34 ◦ C for temperature, and + 15.7%RH for humidity without any climate control system. Operating, the GCH proposed in this paper produces 3.5 kg of wet hydroponic green forage (HGF) for each kilogram of seed (corn or barley) harvested on average. This study centers on the prediction and measurement of the internal air temperature and relative humidity in a GCH built on a farm in the open air, in Quer é taro, M é xico.


Introduction
The development of mathematical models and their subsequent simulations are used to predict the behavior of certain variables in a determinate system that changes over time. In our case, such a system is the growth chamber (GCH) presented in Figure 1. Due to the complex interactions carried out among the variables, such as temperature, humidity, radiation, carbon dioxide, water, nutrients, pests, diseases, and weeds, among others, it is necessary to identify and analyze them in subsystems. Several authors [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20] have characterized the subsystems from different perspectives-for example, climate, nutrition, and irrigation. Models are used to achieve detailed knowledge of all interactions and processes. In some research, through modeling and simulation, the effects of environmental control On the other hand, one would think that a GCH is like a greenhouse; however, there are several differences. For example, a greenhouse is a closed, static, and accessible on foot, and the light crosses the translucent glass or the plastic, letting use of the effect of solar radiation. This structure covers the crop and has at least one device that can be manipulated to control the variables of the environment, nutritional or biotic or both, that influence the growth and development of the crop. Another marked difference between a greenhouse and the GCH is in the commercial aspect because a greenhouse seeks large-scale production that generates economic value, whereas a GCH is used to ensure the production of animal food.

Instrumentation and Environmental Data
The GCH counts with a temperature sensor and a data logger factory certified and calibrated In our proposal, physical principles are used to obtain a mathematical model of the GCH. Due to this, its internal environment is described employing differential equations [26]. This methodology permits researchers to obtain an environmental model of the GCH, relating outside variables (air temperature, relative humidity, global radiation, and wind speed), inside variables (for example, evapotranspiration), and the optical and thermal properties of building materials, along with the volume and surface of the enclosure. In our research, the resolution of equations and the implementation of the simulations were carried out with Matlab®(R2010A, MathWorks, Natick, MA, USA) software.
In this study, the authors propose a novel GCH that integrates the concept of renewable energy taking advantage of natural resources for its operation because it uses a roof-mounted solar panel. Due to its design, the GCH uses a passive system to control temperature and humidity. The GCH does not contain artificial heating, cooling, or lighting systems, and might be relevant to the production of hydroponic green forage. A mathematical model has been developed to estimate the inside air temperature and relative humidity of the GCH. The confidence in the accuracy and reliability of the model is established comparing the simulated values and measured values through the coefficient of determination (R 2 ), and the root mean squared error (RMSE). Results of the comparison reveal good agreement between experimental measurement data and those obtained by the mathematical model. It needs clarifying that it is not the main objective of the current research to develop a complex model with a high level of detail and accuracy, but a simple model that gives satisfactory predictions for realistic environmental responses of the air temperature and humidity inside the GCH. The application of the simulation results of this study may be regarded as a reference to optimize the prediction model of temperature and humidity in other growth chambers under study.

Materials and Methods
This research was carried out on a farm in Querétaro, México (20 • 32'57.4"N, 100 • 30'11.7"W), where the GCH is in the open air. To clarify the explanation, this section consists of the following four parts: (1) Growth chamber, in this part describes the camera and all its components; (2) Instrumentation and environmental data, which deals with the sensors used in the study, along with the equation to filter the information; (3) Environmental model, here the mathematical expressions corresponding to the equations of temperature and relative humidity are introduced to analyze the interior of the GCH; and (4) Initial conditions, in which the initial conditions are established to solve the differential equations.

Growth Chamber
This study centers on the prediction and measurement of the internal air temperature and relative humidity in a GCH built on a farm in the open air, in Querétaro, México.
Its purpose is the growth and development of temporary crops (TCs) of hydroponic green forage (HGF). It is worth mentioning that the crop term used in our work is employed according to the FAO definition (available from http://www.fao.org/es/faodef/fdef11e.htm), i.e., forage used mainly as part of the animal diet, and that can become temporary. Figure 1 presents the scheme of the proposed chamber, and a description is given as follows: Due to its design, the GCH uses passive systems to control the temperature, humidity, and wind speed based on filtering outside conditions without operating any climate control system (heating, cooling, or evaporative).
The GCH integrates the concept of renewable energy because it counts with technology that takes advantage of the natural resources for its operation. Furthermore, it is built with an anticorrosive metallic structure with a tubular rectangular profile in stainless steel of 40 × 40 mm, caliber 12, resulting in rigidity. Additionally, covered with a sandwich-type panel of galvanized steel metal faces and a rigid insulating core of high-density polyurethane PUR (40 kg·m −3 ) with a thickness of 40 mm, model PF600-40L of the company PAINEL2000, which is useful for isolating the interior from the extreme conditions of the elements. On the roof are placed a solar collector, a polycrystalline solar panel of the 125 W brand EPCOM model EPL12512, and a wind extractor brand ATCFANS model BK2000-14 of anodized aluminum. The photovoltaic system ensures the implementation and operation of the GCH anywhere, without relying on conventional electricity, complemented by the following three elements: (1) A solar charge controller with pulse width modulation (PMW), manufactured by PHOCOS model CX40; (2) a RALLY 1200 W inverter model 7466 to convert 12 V DC to 120 V AC; and (3) two solar batteries of a deep cycle of absorbent glass mat (AGM), brand CONER model 31H, 12 V, and 115 Ah. The extractor only rotates because of the external wind present in the GCH installation site (the movement of the wind extractor is ignored, due to the temperature differential). Therefore, ventilation is said to be natural. The lower part of each wall of the chamber (right and left) has two openings, with a stainless-steel mesh 40, SORFORD brand, that prevents the pests or objects entry since the GCH is outdoors. Due to the wind extractor, air enters through these openings and flows upwards. A high-density geomembrane (HDPE) of 1 mm thickness manufactured by GSE covers the wood floor.
The interior of the GCH is divided into two zones: The first part offers semi-dark conditions and must be kept like this during the first eight days, when the seed begins to germinate, and a second part, where the forage shoots are exposed to sunlight to complete their development from the ninth day to the fifteenth day. The lighting is natural thanks to the sun rays that cross the cellular polycarbonate panel (SKYFLEX U-profile 6 mm). In the interior of the GCH, there are 12 rails used to displace the trays horizontally, distributed in two sections, with six rails per section, and rails are separated vertically by 31 cm.
The growth of the HGF takes place in trays of 60 × 40 cm and 7.5 cm high, manufactured by HYDROCULTURA. Each tray has a 35.5 • tilt angle with holes in the bottom, allowing water to flow and drain; unlike conventional production systems that use the practice of run and waste irrigation, i.e., the water on the floor is not reused and wasted. Besides, there is a system for the recirculation of water, which lets considerable water savings. This is one of the advantages of the proposed production system.
Additionally, the GCH has a system of irrigation with the following two components: (1) A storage tank of 1000 L, and (2) a submersible water pump of the 250 W brand ORANGEPUMPS model CSP250C-4. The distribution of the water is carried out with PVC tubes that are 13 mm in diameter, with perforations every 64 cm, where there are installed sprinklers manufactured by the EVERGREEN model MS 8009. The irrigation is conducted using pure water, without nutrients and fertilizers.
The GCH is defined as follows: The GCH is a prototype aimed at mitigating or preventing production losses for small and medium livestock producers, or producers of small animals. It has the following characteristics: It is a small chamber, people cannot walk inside it, and only a part of it contains translucent material. Thanks to its shape and size, the chamber can be relocated without disarming it in its entirety, and for its installation and startup, little time is required. The GCH makes efficient use of space due to the structure composed of several layers in a vertical direction. Due to the optimal isolation of the camera, it exhibits an efficient energy consumption.
On the other hand, one would think that a GCH is like a greenhouse; however, there are several differences. For example, a greenhouse is a closed, static, and accessible on foot, and the light crosses the translucent glass or the plastic, letting use of the effect of solar radiation. This structure covers the crop and has at least one device that can be manipulated to control the variables of the environment, nutritional or biotic or both, that influence the growth and development of the crop. Another marked difference between a greenhouse and the GCH is in the commercial aspect because a greenhouse seeks large-scale production that generates economic value, whereas a GCH is used to ensure the production of animal food.

Instrumentation and Environmental Data
The GCH counts with a temperature sensor and a data logger factory certified and calibrated brand EXTECH model RHT10.The sensor location is 2.08 m behind from the doors, 1.25 m above the floor, and 0.7 m in front of the walls. This device is useful to register and store temperature and relative humidity each day. The resolution and accuracy of the temperature sensor are 0.1 • C and ±1.0 • C, respectively. The resolution of the relative humidity sensor is 0.1%RH. The accuracy changes according to the range of measurement; for the range of 60 to 80%RH, it is of ±3.5%RH, and for 80 to 10%RH, it is of ±5.0%RH. The external weather conditions are obtained from a station Davis Vantage Pro weather, model 6162. This is located outside of the GCH, 4 m to the left from the wall and 3 m from the ground, and permits to capture the following data: Temperature, relative humidity, solar radiation, and wind speed.
In this study, data collected on 11 November 2015 and 2 February 2016 were used to carry out the characterization of the model for simulation purposes. These data were selected because when compared with the information obtained on the other days, they presented major disturbances. The temperature and humidity were acquired every 5 min (measurement interval), obtaining a total of 288 measurements per variable in a day. Other data collected during two consecutive weeks, from 3 July 2016 to 16 July 2016 (15 days were selected because this period corresponds to a cycle of production of HGF), were used to validate the model, i.e., 4320 measurements in 15 days per variable. This measurement interval for data collection was established based on the work reported in [2,8].
It is also worth mentioning that the recorded data were filtered to avoid unwanted peaks during the simulation process. The filter used is called median-mean, and it is useful in statistical analysis. This filter computes the mean filter first, followed by the median operator. Formally, this is expressed by: where, ξ µB represents the median-mean filter size µ; XM µB is the median filter size µ; X µB denotes the mean filter size µ; f represents the input vector; x is the evaluation point; B symbolizes the elemental structuring element; and µ is a size criterion, e.g., if µ = 1, the number of elements is 2µ + 1 = 3, and if µ = 2, then the number of elements is 2µ + 1 = 5, etc. The filter expressed in Equation (1) was implemented in Matlab® R2010A.

Environmental Model
The internal environment of the GCH can be described using a system of two differential equations of the first order, which characterize the behavior of the air temperature and relative humidity [21,27,28] following the mass and energy balances considering the interactions illustrated in Figure 2. External climatic conditions influence this model, as well as the location, orientation, and size of the GCH. The simplified environmental model described in this section takes into consideration next assumptions:  (1) In the temperature model of the air, the heat sources are the heat absorbed through the translucent surfaces, the evapotranspiration of the TC, and the heat transfer through walls and ceiling, while the main source of heat loss is the ventilation. The condensation of the walls (as a result of latent heat) is considered as a parameter neither the water nebulization nor evaporation on the surface of the soil.
(2) The model of humidity in the air considers the evapotranspiration of the crop and ventilation as the main cause of humidity gain and loss, respectively.
(3) The crop evapotranspiration variable utilized in the environmental model comes from some relationships proposed by other authors. Such variables satisfy to be a linear function of the global radiation.
(4) The mathematical relationships of physical processes such as ventilation are simplified. The wind speed within the GCH is considered homogeneous, without considering the multiple crop layers in the GCH. Natural ventilation temperature is despised, due to the small effect of it, with respect to that produced for the wind. Empirical relationships are used to model this physical phenomenon.
(5) The internal temperature and humidity models consider a uniform homogeneous distribution of variables inside the GCH. Figure 2 illustrates that the GCH gains heat due to solar radiation absorbed through the translucent surfaces (Q SR ), the heat provided by crop evapotranspiration (Q CE ), and the heat transfer through walls and roof (Q WR ), while the heat losses are due to infiltrations and the effect of wind (Q IW ). The equation proposed for modeling the behavior of the internal temperature is the following:

The Internal Air Temperature Model
where c p is the specific heat at constant pressure (J·kg −1 dry air ·K −1 ), ρ a is the air density (kg dry air ·m −3 ), H represents the average height of the GCH (m).
The heat absorbed by the translucent surface of the GCH due to solar radiation is estimated using the next equation [29,30]: where Q GR is the external global radiation (W·m −2 ), and C 1 is a dimensionless coefficient that relates to the energy transmittance through the translucent material and it is a property of the materials used during its construction. According to the literature, there are several mathematical models of evapotranspiration for crops [31][32][33], some of which are more complex than others. However, the factor that has demonstrated the highest correlation with evapotranspiration is radiation. Based on the work carried out by Jolliet [34], we assume that crop evapotranspiration is a linear function of solar radiation. For this reason, Equation (4) expresses the term that associates the heat contributed by evapotranspiration (Q CE ) [29]: The heat gain through the walls and the ceiling to the interior of the GCH (Q WR ) consider (1) the conduction through the wall and roof, and (2) the convection and radiation associated with its surfaces. Equation (5) is defined as [35]: where T in − T out , defines the temperature difference between inside and outside of the GCH, R total is the total thermal resistance in m 2 ·°C·W −1 . Q WR can be handled using the thermal resistance concept. Figure 2 (Detail A) illustrates the analogy of the thermal resistance networks.
The Equations (6)-(11) represent the concept of the thermal resistance parameter. This concept is employed, e.g., by Sarafraz et al. [36], for estimating the overall heat transfer coefficient (1/R). The equations defined in [35] are presented as follows: where h conv and h rad are the convection and radiation heat transfer coefficients at the surfaces; ε is the emissivity coefficient; σ corresponds to the Stefan-Boltzmann constant; T cov, in T cov, out are the temperatures in the interior and the outside surfaces, respectively; R 1 represents the thermal resistance for the inner and outer sheets of the galvanized steel metal faces with a thickness of 0.5 mm, respectively; b is the thickness of the insulating core; K mat is the declared thermal conductivity with a value of 0.0211 W·m −1 ·°C −1 . The values are in the technical data sheet of the materials.
Since the convection and radiation resistances are parallel to each other, this may cause some complications in the thermal resistance network. By replacing h conv in Equations (6) and (10) by h conv, rad, in and h conv, rad, out , which include the effects of radiation on the inner and outer sheets of the galvanized steel metal surfaces of walls and roofs (called surface conductance, which represents the combined convection and radiation heat transfer coefficients), avoids all the complications associated with radiation effect. The common values for h conv, rad, in and h conv, rad, out , considering a maximum load, are computed as follows [37]: h conv, rad, in = 8.29 W·m −2 ·°C −1 (summer or winter) (12) h conv, rad, out = 34.0 W·m −2 ·°C −1 (winter) 22.7 W·m −2 ·°C −1 (summer) (13) According to the consideration made for the radiation effect, the resistance configuration is in series, hence the total thermal resistance is obtained below: The energy flow caused by infiltrations and the effect of ventilation is computed below [30]: where ϕ wind represents the airflow through the apertures of the GCH (m 3 ·s −1 ), and A LC is the longitudinal area of the camera (m 2 ), which is the result of multiplying the width by the height of the GCH.
Subsequently, when replacing Equations (3)- (5) and (15) in Equation (2), the mathematical model to predict the internal air temperature of the GCH is:

The Internal Air Relative Humidity Model
For the relative humidity, the source of water vapor was evapotranspiration from TC (E), and the only water loss from the system was due to ventilation. Then, Equation (17) represents the changes in water vapor produced inside the GCH, i.e., Equation (4) establishes the crop evapotranspiration E (kg water ·m −2 ·s −1 ), ω in − ω out represents the difference between the inside and outside absolute humidity of the GCH, and it is calculated utilizing the next equations [38]: where RH in is the internal relative humidity of the chamber, and T in sim is the value simulated by the temperature obtained through Equation (16).

Initial Conditions
The initial conditions (time t = 0) for each one of the state variables of the differential Equations (12) and (14) were assumed to be equal to the measured outside temperature (T out ) and humidity (ω out ). For instance, when recording measurement data at the time 00:00, with the outside temperature sensor reading 15 • C, then T in = T out = 15 • C; similarly, it is done for ω in = ω out . Other initial input values used in these equations were the outside measured variables, such as global radiation and wind speed. Table 1 presents the parameters used in the model.

Results
To evaluate the performance of the mathematical model, we used two criteria: (1) The coefficient of determination (R 2 ) obtained through a regression analysis (for the humidity, a fifth-degree polynomial function is used, due to the nonlinear behavior), and (2) the root means squared error (RMSE). These criteria assessed the correlation between simulated and measured values. They are defined as follows: where T in sim and RH in sim correspond to the simulated values of the temperature and the relative humidity inside the GCH, respectively, and T in mea and RH in mea represent the measured values of the temperature and the relative humidity inside the GCH, respectively. The uncertainties of the experimental results are obtained using the Kline-McClintock equation [39] employed in [36]. Based on the uncertainty analysis, the uncertainties of the T in sim and RH in sim are 12.8% and 8.1%.
The results section consists of the following four parts: (1) Two days are analyzed in different seasons. The parameter R 2 is computed analyzing the differences between the measured and simulated data; (2) parameters RMSE and R 2 are obtained using the data of two consecutive weeks corresponding to the summer season; (3) internal conditions of the GCH, where a percentage is computed to know how much the temperature and humidity maintain ideal conditions without operating any heating or cooling system; and (4) GCH operation example. Here is illustrated the practical functionality of the proposed GCH.

Two Days in Different Seasons
The results of the simulation corresponding to the temperature and relative humidity of two days, 11 November 2015 and 2 February 2016, can be seen in Figure 3a,c,e,g, respectively. For each day are taken into account 288 samples of temperature and relative humidity, 1152 samples are considered in total. From the data graphed in Figure 3, the average difference between the simulated value and the measured value is computed, giving a result of T in sim − T in mea = 0.06 • C for the temperature and RH in sim − RH in mea = −6.1%RH for the relative humidity on a day in autumn, and T in sim − T in mea = 1.2 • C and RH in sim − RH in mea = −7.9%RH for a winter day. Negative signs reveal an underestimation of the mathematical model, i.e., the simulated value was less than the measured value. Figure 3b,d,f,h shows the graphs corresponding to regression analysis for the two days considered in the simulation. The coefficients of determination deduced from the data are the following: (1) R 2 = 0.94 for the temperatures on an autumn day and R 2 = 0.99 for temperatures on a winter day, and (2) R 2 = 0.97 for the relative humidity on an autumn day and R 2 = 0.87 for the relative humidity on a winter day. Comparing the coefficients obtained in this research concerning those reported in [40], whose values are R 2 = 0.92 and R 2 = 0.87 for the temperature and relative humidity in a greenhouse, a little improvement was observed in respect to the temperature, and no difference to the relative humidity.

Two Consecutive Weeks
Similar simulations to those presented in Figure 3 are provided in this section, but now considering 15 days. In Figure 4, the corresponding graphs for the temperature (n = 2016) and relative humidity (n = 2016) for the two weeks (from 3 July 2016 to 16 July 2016) are displayed. In Figure 4, the horizontal lines indicate the days. Figure 4a,c shows a visual check of the measured internal temperature and the simulated (or predicted) internal temperature of the week from July 3 to 9 (first week) and the week from July 10 to 16 (second week), respectively. Figure 4b,d shows the internal temperature regression analysis, with the regression lines that best adjust to the scatter diagram and the linear regression equations with their coefficient of determination values. Similarly, Figure 4e,g shows a visual check of the measured internal relative humidity, and the simulated (or predicted) internal relative humidity of the first and the second week, respectively. Figure 4f,h shows the internal relative humidity regression analysis using a polynomial approximation (to get a best adjust of the data), with the regression lines that best adjust to the scatter diagram and the regression equations with their coefficient of determination values. The RMSE values for the first week are 0.98 • C and 5.5%RH, whereas those for the second week are 1.08 • C and 5.5%RH. In the same way, the RMSE value of the simulated temperature is less than the value reported in [15], of 1.9 • C, and a bit higher than the RMSE value of the relative humidity reported in the same paper, with a value of 3.2%RH for three consecutive days in a greenhouse. The computed RMSE value obtained in our research is slightly higher than that reported in [26], of 0.9 • C, and lower when compared to the relative humidity of 6.6%RH for a single day of simulation inside a greenhouse. Table 2 presents the coefficient of determination and the mean squared errors associated with the obtained simulations through the proposed humidity model, and those obtained from a humidity model that includes (i) a different method to obtain evapotranspiration [32], and (ii) an equation for condensation [41]. In Appendix A, the tests under different scenarios using the mathematical model are presented. In several studies [6,40], the regression analysis applied to the model has as a linear correspondence, but in our case, due to the fluctuations presented by the data, a polynomial approximation is computed since the behavior does not correspond with a straight line. Table 3 shows the R 2 values calculated using linear regression and a fifth-degree polynomial.   In Appendix B, the units of Equations (2) and (17) are verified.

Internal Conditions of the GCH
The variables that most influence the development of the TC of the HGF are temperature and humidity. In the literature, several studies have been carried out under different conditions of controlled temperature and relative humidity, depending on the application-for example, 21-25 • C and 65-75%RH [42]; 22-26 • C and 50-73%RH [43]; 23-25 • C and 6-70%RH [44]; 18-21 • C and 70%RH [45]. In such studies are reported yields of 7-12 kg of wet HGF per kilogram of harvested seed. The temperature and humidity must keep as stable as possible within a specific range to obtain a good crop inside the GCH. If we establish the ideal range for internal conditions (temperature and humidity), based on the previously cited authors, which is between 20 • C and 25 • C for temperature and 70% and 80% for relative humidity, we find that only 31% of the temperature data measured inside the chamber remain in that ideal internal temperature range, although the average internal temperature measured is 20 • C. In the case of relative humidity, only 14% of the humidity data measured inside the chamber are in the ideal humidity range, while the average of internal relative humidity measured is 82%.

GCH Operation Example
The TC examples provided in this subsection serve uniquely to illustrate the functionality of the proposed GCH. This paper just reports the results obtained in the field of forage production using the GCH. Our approach is the modeling and simulation of the internal air temperature and internal relative humidity. In this way, any economic, practical, or nutritional analysis is outside of the scope of the paper. The TCs correspond to the Zea mays and Hordeum vulgare in production cycles of 15 days and 12 days, respectively. The total water used inside the chamber has been generally defined in agronomy [44] as the water added during the irrigation minus the water drained out of the trays. The amount of water used to produce a daily average of forage is equal to the kilograms of wet forage produced on average daily by the GCH among the total water used inside the chamber. Due to the use of recirculation system water, the GCH can produce HGF while saving water. On average, it produces 1.4 kg of corn HGF and 1.7 kg of barley HGF using a liter of water per day in cycles of 15 and 12 days, respectively. Figure 5 illustrates a comparison of the HGF production from Zea mays and Hordeum vulgare. Additionally, the mass of the seed at the time of harvesting (wet seed after being immersed in water for 48 hours) was registered. During the seeding procedure, the seed was weighed dry and recorded. Another measure considered was the mass of wet forage at the time of being recollected. In Appendix C, there are further details on the obtention of HGF.    Table 4 provides information about the kilograms of wet forage obtained for a TC of corn in a cycle of 15 days. Equally, in Table 5, the results for a TC of barley in a cycle of 12 days are presented. The conversion factor relates dry seed to fresh forage of the same sample produced in the GCH, and differs from the conversion ratio (CR), which considers the concentration of dry matter. A conversion factor of 3.5 is obtained on average, i.e., 3.5 kg of fresh HGF is produced for each kilogram of dry seed planted. The GCH presented in this paper only fulfills the specified objectives, and there is not any interest to carry out a commercial or economic analysis. The performance is low when compared with the results reported in [45], with respect to the proposed hydroponic system, which presented values between 4.93 kg and 7.21 kg of fresh forage for each kilogram of barley grain. In [42], 1 kg of forage seed produces a quantity of wet forage between 7 kg and 9 kg. This occurs because the seed utilized in our project is not certified. Therefore, the percentage of germination is unknown; besides, irrigation was carried out using pure water.

Discussion
The results of the model give satisfactory predictions for the air temperature and the humidity in the GCH. Overestimation during the day and underestimation during the night may be due to the unsuitable consideration of a heat transfer coefficient for walls-air or crop-air. Due to the fact that the value of relative humidity depends on the internal airflow through the apertures of the GCH (ϕ wind ) and the evapotranspiration inside of the GCH (E), which is not measured directly in this research, they are only estimated using equations that were developed with experimental data for different conditions than those used in this study.
The results obtained with the proposed humidity model were compared with another system that includes (i) a different method of calculating evapotranspiration [32], and (ii) an equation for condensation [41]. This comparison shows that there are no significant differences in the R 2 parameters obtained by both methods. Our results are similar when applying the mathematical model proposed by Hargreaves [32] where, the authors do not consider analyzing samples in short lapses of time; in our case, every 5 min, radiation was measured-i.e., variations of temperature occur in a longer time, at the beginning and the end of the day, instead of a sampling period. While the equation proposed here is not general, it can be used to compute determined parameters using the captured information of solar radiation. On the other hand, the model introduced by Hargreaves [32] employs tabulated values derived from other experimental procedures.
Moreover, Cengel and Boles [46] have affirmed that, contrary to the relative humidity that varies as the temperature varies, the absolute humidity remains constant in the presence of these variations.
Since the estimation of it is performed using absolute humidity, the results may not reflect the variations that relative humidity would display, and this could explain the low coefficient of determination values for the humidity model.
Similarly, from the point cloud shown in the scatter diagrams for the humidity model in Figure 3f,h and Figure 4f,h, there is a better correspondence between the simulated values and the measured values with a nonlinear function. According to that expressed in [25] and [15], relative humidity is a nonlinear phenomenon because nonlinear equations involved in its modeling govern the transfer of heat and mass. Therefore, linear models like ours suffer from many of the shortcomings by not considering nonlinear effects.
Besides, two factors we have considered for explaining the deviation remaining between the values simulated by the model and the measurements are the following: (1) The complexity of choosing the correct calibration parameters that allow all the state variables that are involved in the model to be adjusted simultaneously, and (2) the assumption that the GCH is a perfect mixing tank, when in fact, the air movement is somewhat limited as a result of the multiple trays inside the GCH. This impacts on the temperature and moisture distribution inside the chamber, since we assume that the air is like a perfectly mixed fluid when there is a spatial distribution.
The simulated data and the values measured experimentally have a good fit. However, the model can be improved. Nevertheless, the use of more complex equations to model humidity within cultivated enclosures does not guarantee better results. For example, in [6], the authors obtained lower coefficients of determination than ours, or the research reported in [40], who obtained slightly higher coefficients than ours. It is so complex to model humidity conditions, that there are even authors who, in their work, only model and simulate temperature, or perhaps they do not show interest in this [11,15,28,30,38].
For 11 November, which corresponds to the autumn season, the measured maximum difference between the internal air temperature of the chamber T in and the external temperature T out was 2.27 • C. However, on average, the internal air temperature was only 0.01 • C lower than the external temperature. In the case of relative humidity, the measured maximum difference between the internal relative humidity (φ in ) and external relative humidity (φ out ) was 20.8%. The internal relative humidity was 11.5% higher than the external relative humidity on average. On 2 February, corresponding to the winter season, the measured maximum difference between the internal air temperature of the chamber and the external temperature was 0.69 • C. However, the internal temperature remains 1.1 • C lower than the external temperature on average. On the other hand, the maximum difference between the internal relative humidity of the chamber and the external relative humidity was 33.2%, and the internal humidity was 11.5% higher than the external humidity on average. In the case of the summer season, two consecutive weeks from 3 July to 16 July were analyzed, the measured maximum differences between the interior and exterior of the chamber were 3.82 • C and 38.8% for temperature and relative humidity, respectively. Similarly, the internal air temperature was only higher than the external temperature by an average of 0.34 • C, and the internal relative humidity was maintained at an average of 15.7% higher than the external humidity.
On the other hand, the R 2 data computed with values 0.96 and 0.94 for the temperature model and of 0.81 and 0.80 for the humidity model for two consecutive weeks are not necessarily bad when compared, concerning other index values published in other papers. For example, in the research presented in [6], the R 2 values computed for the internal air temperature were 0.89, 0.80, and 0.80; similarly, R 2 values of 0.52, 0.73, and 0.77 were obtained for simulated humidity in several runs. As can be seen, these R 2 values are lower than those obtained by the model proposed in this paper. Another example is presented in [40], who obtained an R 2 maximum of 0.92 for the internal temperature and 0.87 for the internal relative humidity. These results are slightly higher than those reported in our research. Notice that the R 2 value of 0.80 is quite common among the published works, which derives from the complexity to predict humidity in closed enclosures for crops.

Conclusions
The simple and efficient mathematical model proposed in this paper to predict the internal air temperature and relative humidity is a useful tool to know the behavior of these variables within the GCH since, from the obtained results, the predicted values present a good correlation with the measured values and with those reported in the literature when comparing the R 2 and RMSE parameters. The simulation results of this study are useful for the design and optimization of environmental control systems applied to a growth chamber. Next, some conclusions from the results are presented: (1) During two consecutive weeks of the summer, the coefficients of determination values were R 2 = 0.96 and R 2 = 0.94 for the model of temperature, and R 2 = 0.83 and R 2 = 0.81 for the model of the relative humidity. For a winter day, these coefficients were R 2 = 0.99 and R 2 = 0.87 for air temperature and relative humidity, respectively. Similarly, the results for an autumn day were R 2 = 0.94 and R 2 = 0.97.
(2) The RMSE values for the temperature model were 0.98 ( • C) and 1.08 ( • C) during two consecutive weeks of summer, and 5.45 (%RH) and 5.48 (%RH) for the model of the humidity in the same period of measure. For an autumn day, the RMSE values computed were 0.78 ( • C) and 6.5 (%RH) for temperature and humidity, respectively. For a winter day, the RMSE values were 1.47 ( • C) and 11.6 (%RH).
(3) The average temperature differences between the simulated and measured air temperatures and simulated and measured relative humidity were no less than 0.1 • C and 1.2%RH for the period from the 3 to 9 July, and for the period from the 10 to 16 July, the differences were no less than 0.2 • C and less than 0.7%RH. Equally, the average differences were more than 0.06 • C (below the resolution of the sensor) and 6.1%RH, corresponding to an autumn day, and for a winter day, the results were no more than 1.2 • C and less than 7.9%RH.
(4) With the use of GCH, the growth of forage in different places and times of the year can be guaranteed, regardless of climatic conditions. Although the interior of the chamber, on average, was only a little colder in autumn ((T in − T out ) = −0.01 • C, below the resolution of the sensor) and winter ((T in − T out ) = −1.1 • C) and only slightly hotter in summer ((T in − T out ) = +0.3 • C), this did not prevent it from fulfilling its objective of producing HGF. It is worth emphasizing that the interior was always wet compared to the external conditions that occurred in the place of experimentation, so materials selected for the construction of the GCH are adequate. (5) The results obtained in the production of HGF are interesting because, despite not having an environmental control, not using certified seed, and using water without nutrients, the average production rate of 3.5 kg of wet HGF was achieved for each kilogram of seed (corn or barley) harvested.  The first scenario, the temperature action does not consider ventilation, irrigation, or forage inside the GCH, i.e., only the heat flowing through the walls, ceilings, and translucent surfaces intervene. Figure A1a present the result, and Equation (2) is written as follows: In the second scenario, natural ventilation without applying irrigation and forage is considered. Here, the variable that relates to the losses due to infiltrations of ventilation intervenes. In Figure A1b, is presented the result, and Equation (2) is expressed as: In the first and second scenarios, an empty GCH is analyzed, so the heat provided by crop evapotranspiration does not intervene.
Third and fourth scenarios: These include all the variables and the contribution from crop evapotranspiration. The results are displayed in Figure A1c,d, respectively. Notice that better simulations occur when increasing the variable number and this indicates that at least four inputs are necessary for the development of a general model. Table A1 shows the coefficients of determination and the variables that intervene in each one of the scenarios. Third and fourth scenarios: These include all the variables and the contribution from crop evapotranspiration. The results are displayed in Figure A1c,d, respectively. Notice that better simulations occur when increasing the variable number and this indicates that at least four inputs are necessary for the development of a general model. Table A1 shows the coefficients of determination and the variables that intervene in each one of the scenarios.

Appendix B
In this appendix the units of Equations (2) and (17) of the mathematical model are explained in detail.

Appendix C
Appendix C describes the seeding procedure. Consider maize seed (Zea mays) in a production cycle of 15 days. The process of the production of HGF was carried out through the following procedure: (1) Screening the seed: The seed was sieved to eliminate foreign agents; (2) cleaning of the seed: Impurities and fragments of seed not removed with the sieve were manually eliminated; (3) disinfecting and washing: Seeds were washed using tap water, and they were immersed posteriorly, in a solution of sodium hypochlorite for 10 min; (4) seed immersion: Corn seed was soaked in a container with water for 48 h to speed up germination, and posteriorly, water was changed after 24 h; (5) sowing in trays: The seed was moved to the cleaned and disinfected plastic trays and distributed without exceeding 1 cm in height; (6) harvesting: Once the trays contained the seeds, they were transported immediately to the interior of the GCH. The trays were situated in a semi-dark place on the rails (see Figure 1), advancing horizontally every time when placed inside; (7) irrigation: The micro-sprinklers were activated for 2 min each, with a total of eight irrigations per day; (8) growth: From the third day, the buds of the seedlings started growing, until completing their cycle on the fifteenth day.
During the production cycle, the variables of interest, humidity, and temperature in the GCH were collected.