Predictive Control Applied to a Solar Desalination Plant Connected to a Greenhouse with Daily Variation of Irrigation Water Demand

The water deficit in the Mediterranean area is a known matter severely affecting agriculture. One way to avoid the aquifers’ exploitation is to supply water to crops by using thermal desalination processes. Moreover, in order to guarantee long-term sustainability, the required thermal energy for the desalination process can be provided by solar energy. This paper shows simulations for a case study in which a solar multi-effect distillation plant produces water for irrigation purposes. Detailed models of the involved systems are the base of a predictive controller to operate the desalination plant and fulfil the water demanded by the crops.


Introduction
Modern agricultural systems are characterized by the intensive and optimal use of land and water, turning agricultural exploitation into a semi-industrial concept.Greenhouses are systems suitable both for zones with unfavourable climatic conditions, allowing crop growth regardless of the ambient temperature, and for regions with less restrictive weather, with the aim in this case of increasing crop productivity and improving fruit quality.The greenhouse environment is ideal for farming, because there are variables that can be manipulated to achieve optimal growth and plant development.Control of greenhouse crop growth requires energy consumption, depending on the crop's physiological requirements and, additionally, depending on the production patterns adopted for yield quantity and timing.This crop growth control is achieved through adequate manipulation of climatic variables of the environment and the amount of water and fertilizers applied through irrigation.
Productivity optimization through efficient and adequate irrigation is a basic objective in those countries with water limitations.This resource scarcity is increasing due to the recent rapid expansion of the surface area occupied by greenhouses in the Mediterranean Basin.Consequently, this has also led water to become an important actor in the sustainability of the greenhouse-based system in southeastern Spain.
This water consumption has been progressively depleting the aquifers in the area [1].Eighty percent of the irrigation water used in Almería (Spain) comes from underground sources, leading to localized over exploitation of aquifers [2].Over the last few years, as in other arid and semi-arid areas of the world, the use of alternative water sources has been promoted, such as purified water, rain and condensed water collection as a secondary source, the reuse of drainage water, the development of new technologies related to water-use efficiency, such as advanced irrigation controllers, and sea water desalination.
In this line, the idea of integrating solar desalination systems in the agricultural environment has been significantly considered with the aim of dealing with water limitations in some regions of the planet.One of the most simple and inexpensive techniques, solar stills [3], can be easily combined with greenhouses [4].Nevertheless, as explained in [4], water produced by a solar still is not enough for growing a crop.This paper deals with the combination of a greenhouse and solar multi-effect distillation (MED) unit.It consists of taking advantage of the water produced in the MED unit to feed a greenhouse, both systems being located in the southeast of Spain.The challenge is to properly operate the desalination plant to produce the daily water demanded by the crop.
Most of the optimization algorithms applied to desalination processes deal with the design of cogeneration plants [5][6][7][8].The objective function includes thermodynamic and environmental issues.However, few optimization problems are focused on satisfying a variable water demand.In the present paper, a two-layer controller that combines basic controllers with model predictive control (MPC) is applied to find and follow proper temperature setpoints in the MED unit and in the solar field.The idea is to analyse if it is possible to maintain a desired volume of distillate obtained by means of a solar desalination process despite the variable consumed water by irrigation in the greenhouse.

Case Study
The case study explored in this paper is a micro-grid framework in which two interconnected plants must be managed; a greenhouse and a solar desalination plant (see Figure 1).On the one hand, the greenhouse demands fresh water daily for irrigation purposes, and on the other hand, a solar desalination plant produces distillate water in a multi-effect distillation unit.An intermediate storage tank is assumed to be located between the production process and the consumer system.An explanation of the facilities involved is included in the following sections.

Solar Desalination Facility
The desalination plant used in this study is the AQUASOL system [9] (see Figure 2) situated at Plataforma Solar de Almería (PSA), a dependency of the Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT) in the south of Spain.This pilot plant includes a forward-feed vertically-stacked MED unit with 14 cells coupled with a solar collector field (although the system was designed to operate also with fossil energy, this work only deals with the solar operation mode).As shown in Figure 1, seawater is pumped and preheated towards the first cell evaporator (or heater) of the MED.There, water is sprayed on the external surface of the horizontal tube bundle, and part of the seawater is evaporated by the release of the sensible heat of the hot water.A small part of the steam generated in this evaporator preheats the seawater, whereas the main part reaches the evaporator of the next stage, providing the required thermal power to continue the partial evaporation of the feedwater.Similar condensation/evaporation processes are repeated from the second to the 14th evaporators.The steam generated in the last effect is condensed in the final condenser (cooled by seawater).A more detailed description of the MED unit can be found in [10].
The required hot water for the heater is provided by the solar field that supplies hot water to the storage system (two 12-m 3 water storage tanks).A three-way regulation valve, γ M , is used to reach the nominal temperature at the inlet of the heater, by mixing water from the hot tank with that returned from the heater.
When the solar field temperature exceeds the one in the hot tank, the on-off valve position, γ t , is changed to connect these components.Otherwise, the solar field should be connected to the bottom part of the cold water tank to avoid cooling down the hot water tank.

Greenhouse Environment
The research greenhouse used in this work (see Figure 3) is located at the Experimental Station of Cajamar Foundation (El Ejido) in the province of Almería in SE Spain (2 • 43'W, 36 • 48'N, and 151-m elevation).The crops were grown in two multi-tunnel greenhouses.Each greenhouse has an area of 800 m 2 , with a crop area of 616 m 2 with polyethylene cover; each has automated ventilation with lateral windows in the northern and southern walls, a flap roof window in each span and mesh-protected anti-trip bio-nets, 20 × 10 in thickness.The orientation of the greenhouses is E-W; with crop rows aligned N-S.Cropping conditions and crop management were very similar to those in commercial greenhouses.
A meteorological station is available outside the greenhouses, where air temperature and relative humidity, solar and photosynthetic active radiation, rain detection, wind direction and velocity measurements are taken.During the experiments, some indoor climate variables were also measured: air temperature, relative humidity, solar radiation, photosynthetic active radiation (PAR), CO 2 monitoring, soil temperature at 3 cm and 40 cm, leaf and substrate temperature, leaf wetness, as well as electrical conductivity and pH monitoring in irrigation and drainage water.All sensors were located in the centre of the greenhouse; the psychrometer and CO 2 sensor at a height of 2.5 m are just above the mature crop; and the PAR and solar radiation sensors at a height of 3.5 m.Table 1 shows the different sensors installed in the greenhouse, its model and brand, range and accuracy.The crop grows in coconut coir bags with six plants and three droppers each.The irrigation is automated by a demand tray, which is formed by two crop bags.Drainage water is set at a 20% volume.All of the actuators are driven by relays designed for this task.All data are recorded every 30 s with a personal computer.

Dynamic Model of the Case Study
The proposed case study is used to simulate a control system, which manages the solar desalination plant operation to fulfil the water demanded by the crops.Therefore, it has been necessary to develop a model able to represent the behaviour of both processes properly.A description of the variables and subscripts involved are included in Nomenclature

Greenhouse Model
In general terms, without taking into account the specific subsystems for each greenhouse, the air temperature T a,int can be modelled using the following balance [11]: where Q cnv,cv−a is the convective flux with the cover, Q cnv,ss−a is the convective flux with the soil surface, Q ven is the heat lost by natural ventilation and the heat lost by infiltration losses, Q trp,cr is the latent heat effect from crop transpiration and (c p,a ρ a V g )/A ss is the product of specific air heat, the air density and the effective height of the greenhouse defined as the greenhouse volume and soil surface area.
The convective heat transfer per unit area from the inside air to the cover is calculated based on the difference between the cover temperature, T cv , and the greenhouse air temperature, T a,int [11]: where V cv is the cover volume and h cnv,cv−a is the cover's inside convective heat transfer coefficient based on the difference between the cover temperature and the internal air temperature and the mean greenhouse air speed.In order to simplify the model, the approach proposed in [12] has been used: if the internal air temperature is higher than the cover temperature, the heat transfer is turbulent; otherwise, the heat flux is laminar.
The convective heat transfer per unit area from inside air to soil surface, T ss , is calculated in the same way as the cover convective fluxes using the following equation [11]: where h cnv,ss−a is the inside soil surface convective heat transfer coefficient based on the difference between the soil surface temperature and the internal air temperature and the mean greenhouse air speed over the soil surface.Using the studies in [13], the conductive flux between the soil surface and the first soil layer is calculated based on the assumption that the heat flux is one-dimensional (Z axis).

The heat lost by natural ventilation and infiltration (Q ven ).
The heat lost per unit area by the natural ventilation term is modelled according to the American Society of Agricultural Engineers (ASAE) Standard Heating, Ventilating, and Cooling Greenhouses (EP406.3)[14] as a function of volumetric flow rate, q ven , greenhouse temperature and external air temperature (T a,ext ) [15]: where q ven is described in Equation ( 5) and T a,exh is the exhaust air temperature, calculated as a linear combination of external and internal air temperatures [16].This term includes the heat lost by infiltration losses, as shown in the volumetric flow rate equation, which is based on the thermal buoyancy and wind forces and is described as: where q loss is the leakage when the vent is closed, n ven is the number of vents, l ven is the length of the vents, d ven is the discharge coefficient, g is the gravity constant, k w is the wind-effect coefficient, l ven,cord is the cord joining the two extremities of the vent based on the position of the vent, U ven , using the following equation: where w ven is the width of the vent.The leakage when the vent is closed, q loss , is based on the external wind speed (v w,ext ): where v w,ext,lim is the wind speed, considered as the limit between high and low wind speeds.
The crop itself affects the greenhouse air temperature.As no leaf area measurements are available, it is not possible to use a convective factor in the heat balance equation functioning as a boundary variable.One way to model the effect of the crop on the air temperature is based on the latent heat caused by plant transpiration, described by the following equation [17]: where Ṁtrp,cr is the crop transpiration and λ is the latent heat of evaporation.Most transpiration estimators are based on the Penman-Monteith equation [17,18].This equation essentially combines the heat transfer equation between the crop and the surrounding air mass.A simplified pseudophysical transpiration model can be used, based on two main variables: the solar radiation arriving at a particular depth in the plant canopy and the vapour pressure deficit [17]: where k l is the light extinction coefficient for the crop (this is related to the leaf inclination angle and the leaf arrangement with regard to the Leaf Area Index, LAI, and provides an indication of the plant's efficiency at intercepting solar radiation), DPV is the vapour pressure deficit and I crop is the global radiation reaching the crop.The model of absolute humidity, H a,int (water vapour content of the greenhouse air), is based on a water vapour mass balance equation.The vapour outflow takes place through condensation on the internal side of the cover, crop transpiration, ventilation and vapour lost by infiltration losses.As artificial water influxes (cooling, fogging, etc.) are not installed in the greenhouses in which the experiments were carried out, H a,int is modelled using the water mass balance equation given by Equation (10).Based on the mass flux processes, the greenhouse air humidity can be modelled using the equation developed in [11]: where Ṁtrp,cr is the crop transpiration flux described in Equation ( 9), Ṁevp,ss is the soil surface evaporation, Ṁcd,cv is the condensation flux from the cover, Ṁven,int−ext is the outflow by natural ventilation, Ṁloss are the humidity losses by infiltration, and it can be obtained from Equation (7).Ṁven,int−ext is described by the following equation [11]: where H a,ext is the external air absolute humidity, and the volumetric flow rate, q ven , is described in Equation (5).

Solar Desalination Plant Model
The model of the solar desalination plant is divided into three main components; the solar field, the storage system and the MED unit.Figure 4 shows a diagram with the connections between these models.Solar field outlet temperature, T oF , is obtained with a lumped-parameter model based on an energy balance [19], and the storage tanks have been characterized with energy and mass balances [20].For the MED unit, the dynamic model developed in Modelica [21] and validated with data from the PSA facility has been used as the reference.With the aim of reducing the computational effort, this model has been simulated to obtain two first-order models that give the distillate mass flow rate, ṁd , and the outlet MED temperature, T oM , as a function of the inlet MED temperature, T iM , assuming nominal conditions in the heater flow rate, ṁM =12 L/s, seawater temperature of 22.3 °C, 4.4 kg/s as the condenser mass flow rate and 2.2 kg/s as the seawater mass flow rate: To obtain the above models, the step response method has been applied [22]; the input variable, T iM , was changed with different amplitudes, and the output signals, ṁd and T oM were recorded to identify the relation between the input and the outputs.The dynamic responses have been represented with the Laplace transform, because its algebraic form is simple [23]; the transfer functions are useful to extract systems information; and it is easy to include them in the simulation tool (MATLAB).
A comparison between the results obtained with both models is shown in Figure 5.

Global Irradiance
The study presented in this paper combines two facilities (the greenhouse and the solar distillate plant), which are located in two different research centres and whose irradiance measurements cannot be compared.Since the greenhouse model requires more meteorological variables as inputs, the solar desalination model uses solar irradiance measurements from the greenhouse as input.Nevertheless, this irradiance is measured on a horizontal surface, and it has been transformed to solar irradiance on a tilted surface [24] to take into account the solar field inclination (36 • ).

Control Scheme
The idea behind the controller proposed in this section is to maintain a desired volume of distillate, V d,re f , taking into account the variable quantity of water demanded by the greenhouse.The proposed scheme (see Figure 6) includes a state machine, a reference layer and a regulatory layer.The models presented in Section 3 are used to estimate the future water consumption and desalination system temperatures that are used by the reference layer to evaluate optimal setpoints for the temperatures at the inlet of the MED heater and at the outlet of the solar field.These two variables can be regulated, in the regulatory layer, by means of the valve γ t and water mass flow rate, ṁF , respectively.Both layers are activated when the MED unit is running (CM=1), where CM is defined by the state machine on the basis of the rules explicitly imposed (see Section 4.1).

State Machine
The state machine chooses one of these four operating modes: • MED in solar mode: the MED unit is running, and the heated water by the solar field flows to the hot tank.CM = 1, γ t = 0. • MED in solar recirculation mode: the MED unit is running, and the valve γ t is switched, so that the tanks and the solar field are decoupled.CM = 1, γ t = 1.This mode is used to avoid cooling situations in the storage system due to a downfall in the solar field temperature.• Solar-field mode: the MED unit is stopped because the temperature is not high enough in the hot tank, and the heated water by the solar field flows to the hot tank.CM = 0, γ t = 0. • Solar-field recirculation mode: the MED unit is stopped, and the solar field water recirculates to increase its temperature.CM = 0, γ t = 1.
The inputs required by this state machine are the signal PlantReady, temperature in the hot tank, T h , and temperature at the outlet of the solar field, T oF .
This state machine includes the following conditions: where PlantReady is a manual input provided by the plant operator who must take into account the irradiance conditions and check that the vacuum system (to remove the air and the non-condensable gases generated during the desalination process) is ready at the MED unit.

Reference Layer
To reach the desired volume of distillate water, two previous steps are recommended: to calculate the required MED inlet temperature and to choose the appropriate solar field temperature setpoint with the aim of reaching a high enough temperature in the hot tank.
As described in Section 3.2, the distillate production and MED inlet temperature are linearly dependent, so a linear model predictive controller (LMPC) can be used to calculate the setpoint, T iM,re f , taking into account the temperature constraints: there is a maximum temperature, to avoid scale formation in the heater, 72 • C, and a minimum temperature is also defined to assure a minimum distillate production: 55 On the other hand, a nonlinear model predictive controller (NMPC) is used to solve an optimization problem, to calculate a reference for the outlet solar-field temperature and to assure the required distillate production.This controller must be applied over the temperature instead of the water flow rate to include the following constraints: The solar-field outlet-inlet temperature difference must be lower than 20 • C to avoid stress in the absorber tubes; a minimum temperature difference must be guaranteed to avoid cooling down the water; T oF must be under 95 • C to avoid evaporation, and it should be higher than T h not to cool down the stored water.
Although most of the MPCs are applied to linear models, some techniques have been developed to obtain the future control actions using a nonlinear model of the process [25].In this work, the Nonlinear Extended Prediction Self-Adaptive Control (NEPSAC) approach [26] has been chosen, because it directly uses the nonlinear prediction model without local linearization, and it solves the optimization problem in a low computational time.The idea behind this technique is to use the nonlinear model to predict the base and the step responses.

Regulatory Layer
The outlet solar-field temperature, T oF , can be controlled by varying the solar-field water mass flow rate, ṁF , in its operating range [1, 3.2] kg/s.Several algorithms have been tested in this plant, obtaining successful results [19,[27][28][29][30][31].In this case, the feedback linearization control (FLC) is used [19].A PI controller, with the valve γ M as the control variable, is implemented to reach the desired T iM .

Simulation Results
This section shows the results of the proposed control scheme under different meteorological conditions and with different curves of water volume demanded by the greenhouse.A comparison with the results obtained when no reference layer is included is also presented at the end of this section.With the aim of scaling the greenhouse water consumption to the solar desalination production, 10 greenhouses have been considered in the simulation.Therefore, the real consumption values obtained from the greenhouse have been multiplied by 10.The setpoint in the distillate tank has been established to 36 m 3 in order to store water for three days in case of cloudy days.

Sunny Day
In this example, the solar irradiance follows the Sun curve almost perfectly (Figure 7a).At the beginning, the stored water is equal to the setpoint (Figure 7d), and the temperatures in the cold and hot tanks are 72 °C and 76 °C, respectively (see Figure 7c).When irradiance reaches 500 W/m 2 , the PlantReady signal is switched to one, and the MED unit is turned on (at 10:30 in Figure 7b; CM = 1), as well as the reference layer is activated.Since the estimated consumed water curve is under the real consumption (Figure 7g, blue and red lines, respectively), the inlet MED temperature setpoints, T iM,re f , have lower values most of the time (between 55 °C and 60 °C, as can be observed in Figure 7h).The highest setpoint (68 °C at 11:00) is a consequence of the highest difference (3.7 m 3 ) between the distillate volume and the setpoint (V d and V d,re f in Figure 7d).In the regulatory layer, the inlet MED temperatures are properly controlled with the valve aperture, which regulates the mass flow rate (Figure 7i).
To reach the desired temperatures at the inlet of the MED unit (Figure 7h, dashed line), the hot tank must also reach these temperatures, T h,re f = T iM,re f .However, in this case, T h (Figure 7c, dashed line) is quite higher than T iM,re f (Figure 7e, dashed line), so the NMPC calculates that, with those irradiance conditions, the outlet temperature setpoint, T oF,re f , must be the maximum one, so the water mass flow rate is close to the minimum one (Figure 7e), reducing the thermal energy delivered to the hot tank.As a consequence, the electricity consumption from the solar field pump is also reduced.
It is important to mention that, when the solar field is in recirculation mode (γ t = 1 in Figure 7b), the NMPC output is different.In those situations, the calculated setpoints (T oF,re f in Figure 7e) are close to the minimum allowed to extend the recirculation mode period and reduce the temperature in the tanks.Although these situations happen in small periods of time (usually at startup and shut-down), this is a non-desirable consequence that should be avoided.

Partly Cloudy Day
In this example, the irradiance curve shows low fluctuations due to light clouds (Figure 8a), and the water consumption slightly overestimates the real demand (Figure 8g).
At the beginning, when the MED unit is started at 10:30 (Figure 8b, red line), the volume of distillate water is lower than the desired one (Figure 8g), so the inlet MED temperature setpoint is 69 °C (Figure 8h, dashed line) to compensate the estimated consumption.Due to the low temperature in the hot tank (Figure 8b, red line), this MED temperature is not feasible.The mass flow rate from the hot tank to the MED unit is the maximum one producing a quick stored thermal energy reduction (T h falls down 9 °C in 1 h, as can be observed in Figure 8c).At 11:30, the stored distillate volume, V d , is higher than the one desired, V d,re f (Figure 8d), and the inlet MED temperature setpoint (Figure 8h, dashed line) falls down to the minimum value (55 °C).As can be observed, this temperature, which assures low water production, is maintained during six hours.Nevertheless, the water production exceeds by 4 m 3 the one demanded, highlighting the challenge of coupling the production-demand variabilities.In order to maintain the temperature in the hot tank close to T iM,re f , the reference layer must change the solar field temperature setpoint (Figure 8e, dashed line), depending on the irradiance level, the solar field mode and the tanks temperatures.In this case, the solar field water mass flow rate (Figure 8e) shows more variation due to the irradiance disturbances.This example also shows a second situation that must be avoided.When the temperature in the hot tank is too much higher than the one desired in the MED unit, the mass flow rate from the hot tank to the MED unit can be higher than the solar field water mass flow rate, and part of the water flows from the hot to the cold tank (inverted flow situation).Under these conditions, the NMPC output gives lower solar field temperature setpoints to increase the solar field water mass flow rate, increase the quantity of water that flows from the hot to the cold tank and reduce the hot tank temperature more quickly.Although this fact happens in very reduced moments, obviously, this situation is undesirable, because stored thermal energy is being lost, and the electrical energy consumption from the solar field pump is being increased.

Ten Days
The proposed control strategy has been simulated during 10 days using greenhouse irradiance data measured in April 2015 (see Figure 9a).As Figure 9a shows, during the first five days, the distillate volume is maintained around the setpoint, with the higher deviation (Day 3) caused, in part, by a water consumption overestimation.During the following four days, the low solar irradiance level causes a severe fall in the distillate volume, and two weeks with good solar irradiance will be required to recover the desired level.This situation reveals the necessity of considering an auxiliary system to feed the MED unit or to have a backup of water from the district network if available.It must be pointed out that on the seventh day, the MED unit is not operated due to the low irradiance level.On the 10th day, the clouds disappear; the solar distillate plant can be operated again during almost six hours, and the volume of distillate starts its recovery.
Figure 9 shows the water and the solar field electrical energy consumptions.The mean error between the real and simulated water consumption is 0.07 m 3 , and the maximum deviation is 0.16 m 3 on the third day.The electrical energy consumption depends on the irradiance level, volume of distillate and the MED unit operating hours; on sunny days in which the volume of distillate is close to the setpoint, such as the first day, the MPC provides low inlet MED temperature setpoints to reduce the water production.Therefore, lower tank temperatures are required, and the NMPC gives high solar field temperature setpoints, because it requires lower mass flow rates, which means a thermal energy reduction.For this reason, the electrical energy related to the pump is low.Nevertheless, if the distillate level is far from the setpoint and the water consumption estimation is high, the controller tries to increase the thermal energy from the solar field, and the electrical energy increases (see the second day).On the third day, the situation could be similar to the first day, but the electrical energy increases due to the undesirable situation (inverted flow) explained in Section 5.2.When the irradiance level is low, the electrical energy is also reduced, because the operating time is considerably lower, and the mass flow rate must be around the minimum value.In days, such as the 10th one (good irradiance level, but poor quantity of stored distillate), the system requires high solar field mass flow rates to increase the thermal energy and the water production, causing also an increment in the electrical energy consumption.

Results without the Reference Layer
In order to observe the benefits of using the reference layer, the previous results have been compared to a case without a reference layer (see Table 2).In these simulations, it is considered T iM,re f = T oF,re f =72 °C.The advantages are: lower electrical energy consumed in the solar field, P F , higher energy storage in the tanks, E a , and a lower integral of the absolute error (IAE), which is a performance index defined as |V d,re f − V d |dt.The disadvantage is that the MED unit must be operated more hours, increasing then the electrical energy consumed in the plant.

Conclusions
The use of an appropriate control scheme in a solar desalination plant for greenhouse irrigation purposes could assure the water demand, reduce electricity costs in the solar field pump and maintain more thermal energy in the storage system.Trying to reduce the cost of the water produced in a solar desalination plant, a hierarchical controller has been proposed to maintain a distillate volume for greenhouse irrigation.Promising results have been obtained from the fact that dynamic models of both demand and production systems can be the key for reaching more competitive distillate prices when model predictive controllers are applied.Nevertheless, some remarks have to be highlighted to improve the results: (i) the control objective should be changed to minimize the distillate cost taking into account distillate volume constraints; (ii) an external source of energy could be included to assure a continuous production of water (when a water backup is not available); (iii) simulation in different seasons should be included to consider the variability in the crop state and in the hydric requirements.

Figure 1 .
Figure 1.Solar desalination plant coupled with a greenhouse.

Figure 2 .
Figure 2. Solar desalination facility at PSA-CIEMAT.From left to right and from top to bottom: cold and hot tanks, MED unit and solar field.

Figure 3 .
Figure 3. Greenhouse facilities used for the experiments performed in this work.From left to right and from top to bottom: greenhouse, dropper, solar and PAR radiation device at the outside, irrigation system, solar and PAR radiation device inside the greenhouse and the tomato crop lines.

Figure 4 .
Figure 4. Connection between the submodels of the solar desalination plant.

Figure 5 .
Figure 5. Results obtained with the MED unit first-order model.

Figure 6 .
Figure 6.Proposed control scheme to obtain a desired volume of distillate

Figure 7 .
Figure 7. Simulation results obtained for a sunny day.(a) Irradiance; (b) Valve γ t , MED unit state and PlantReady; (c) Temperature in the tanks; (d) Volume of distillate; (e) Solar field temperatures; (f) Solar field mass flow rate; (g) 800-m 2 greenhouse water consumption; (h) Heater temperatures; (i) Heater mass flow rate;

Figure 8 .
Figure 8. Simulation results obtained for a partly cloudy day.(a) Irradiance; (b) Valve γ t , MED unit state and PlantReady; (c) Temperature in the tanks; (d) Volume of distillate; (e) Solar field temperatures; (f) Solar field mass flow rate; (g) 800-m 2 greenhouse water consumption; (h) Heater temperatures; (i) Heater mass flow rate;

Figure 9 .
Figure 9. Simulation results obtained during ten days.(a) Irradiance and distillate volume; (a) 800-m 2 greenhouse water consumption and solar field electrical energy consumption

Table 1 .
Different variables sampled with the main characteristics of the sensors installed.