Dynamic Energy Exchange Modelling for a Plastic-Covered Multi-Span Greenhouse Utilizing a Thermal Effluent from Power Plant

To utilize the energy in the thermal effluent, many attempts have been made to use the thermal effluent for agricultural facilities such as greenhouses. As the first step, it is important to estimate the energy loads of the greenhouse for deciding a suitable scale for the heating and cooling. Then, it is available to estimate the energy efficiency of the thermal effluent heat pump system installed in the greenhouse. Therefore, the main objectives of this study were to design and validate an energy model of the experimental greenhouse growing Irwin mangoes and to estimate the annual and maximum energy loads using building energy simulation (BES). Field experiments were conducted in a multi-span plastic-covered greenhouse growing Irwin mangoes to measure the internal environments of the greenhouse and crop characteristics. The energy exchange model of the greenhouse considering crop, cladding, heat pump was developed using BES. The BES model was validated using the data measured at field experiments. The designed model was found to be able to provide satisfactory estimates of the changes of the internal air temperature of the greenhouse (R2 = 0.94 and d = 0.97). The hourly energy loads computed by using the validated model were used to analyse the periodic and maximum energy loads according to the growth stage of the cultivated crops. Finally, the energy costs were compared according to the type of energy source based on the calculated annual energy loads. The average energy cost when using the thermal effluent—heat pump system was found to be 68.21% lower than that when a kerosene boiler was used.


Introduction
As South Korea has four distinct seasons, it is difficult to produce crops via field culture throughout the year. Greenhouses can control the growing environment of crops and produce high-quality crops all through the year. Therefore, the greenhouse cultivation in South Korea increased from 23,669 ha in 1990 to 51,226 ha in 2018. To maintain optimum growing environments for the crops, approximately 30% of all greenhouses utilized cooling and heating systems. Among them, approximately 85% of all greenhouse farmers used fossil fuel as an energy source [1]. Additionally, the heating cost of the greenhouse have been a large percentage of the total production cost because of four distinctive seasons in South Korea. Accordingly, farmers have been burdened with production costs such as energy costs.
After the Industrial Revolution, the dependence on fossil fuels increased globally due to the use of petroleum, natural gas, and coal in machines. Moreover, the world population has increased approximately 1.7% per year for the last 50 years [2]. As energy consumption per capita increased from 1.46 ton of oil equivalent (TOE) year −1 in 1981 to 1.86 TOE year −1 in 2017, the energy consumption around the world has also increased from 6582.9 million TOE year −1 to 14,034.9 million TOE year −1 . Along with a dramatic rise in energy consumption, global issues such as global warming, climate change, and energy crunch have arisen, which have attracted interest in new energy sources for replacing fossil fuels. Therefore, attempts have been made to utilize water power, wind power, geothermal heat, and other sources. While the country underwent rapid development, the energy consumption in South Korea increased from 38.9 million TOE year −1 in 1981 to 232.7 million TOE year −1 in 2018. In addition, due to high population density and lack of natural resources, South Korea imported 93.7% of the total energy consumption from abroad in 2018 [3]. Thus, the South Korea government established the Energy Act, Energy Use Rationalization Act, Sustainable Development Act, and other acts with the objective of reducing the total energy consumption and improving the energy efficiency. In addition, the Ministry of Agriculture, Food and Rural Affairs provides financial aid to greenhouse farms that use renewable energy to lower the reliance on imported fossil fuels and improve the efficiency of energy usage [4].
In South Korea, more than 80% of the total energy production was from thermal power generation or nuclear power generation. During energy production, thermal and nuclear power plants waste approximately 56% of the total energy input by heating up the engine. The power plants mostly use seawater to cool down the engine and discharge the heated seawater, which was known as thermal effluent. As the thermal effluent absorbed the waste heat of the power plant, it contains a massive amount of thermal energy. Thermal and nuclear power plants around the country discharge approximately 55.2 billion tons of thermal effluent per year. This includes approximately 388,000 GWh of thermal energy, which can be utilized as an energy source. The total amount of thermal energy in the thermal effluent can easily cover the entire heating energy for greenhouses in South Korea and can reduce the energy cost of greenhouses [5]. To use the thermal energy of the thermal effluent for the greenhouse heating and cooling systems, various types of equipment, such as water storage tanks, heat pumps, heat storage tanks, and fan coil units, must be considered. To use energy more efficiently, each component should be properly designed based on the energy loads of the greenhouses. As the amount of available thermal energy in the thermal effluent changes with time, the energy loads of greenhouses should be estimated as functions of a time variable, such as the season or the crop growth stage.
The methods for estimating the building energy loads can be divided into static methods and dynamic methods according to the time factor. Energy loads are calculated while assuming that the internal and external environments are in a steady state in a static method. In the case of greenhouses, regardless of the type, maximum energy loads can be calculated based on the ambient air temperature and the solar radiation on the warmest and coldest days. Periodic energy loads can be calculated from the accumulated difference between the optimum growth temperature of cultivation crops and the external air temperature [6,7]. The static methods have the disadvantage of not considering the thermal storage of structure by solar radiation and internal heat sources [8]. This can be problematic, especially in greenhouses, as the cladding is much thinner than in general buildings and, hence, is much more sensitive to changes in weather conditions. In addition, the energy exchange between crops in the greenhouse and the ambient air change drastically over time [9]. Therefore, when calculating energy loads of greenhouses via a static method, the calculated energy loads were different from the actual energy loads. In contrast, a dynamic method considers the actual energy flow under an unsteady state. A dynamic method is a method for calculating the energy balance of a structure over time via numerical analysis. With a dynamic method, energy loads are calculated by simulating the energy exchange using real-time weather data. The energy balance equation for a specified time was analysed, and the amounts of energy exchange and thermal storage of each component were passed to the next time step. In architecture, a building energy simulation (BES) technique was typically used as a dynamic method for calculating energy loads with high accuracy and performance. Recently, this technique has also been used in agriculture to simulate energy exchange in greenhouses dynamically and to validate the dynamic energy exchange model of a greenhouse [10][11][12][13][14][15]. Therefore, it is necessary to use a dynamic method to simulate the energy exchange in greenhouses more precisely by considering the internal and external environmental changes of the greenhouse.
The objective of this study was to design and validate an energy model of the experimental greenhouse growing Irwin mangoes and to estimate the annual and maximum energy loads using building energy simulation (BES). The experimental greenhouse was constructed near the power plant to enhance the energy efficiency of the heating and cooling system by utilizing a thermal effluent. The system for using thermal effluent was installed at the experimental greenhouse without quantitative analysis and design. The dynamic energy exchange model of the greenhouse was developed considering the greenhouse structure, the operational conditions of the thermal screens and vent openings, an energy exchange of the crops, and the thermal effluent-heat pump. The periodic and maximum energy loads were calculated using the developed energy exchange model. Then, the impact of the energy exchange by crops was analysed when calculating the energy loads of the greenhouse. Additionally, the energy costs according to the utilization of the thermal effluent were analysed.

Research Flow
A flowchart for the design of a dynamic energy exchange model of the greenhouse and the calculation of the energy loads is presented in Figure 1. The weather data of the target area were analysed and field experiments were conducted to design an energy model of the greenhouse growing Irwin mangoes. Greenhouse structural characteristics, such as the shape, cladding, and heating and cooling system, were investigated for the design of the greenhouse structure model, and crop characteristics, such as the stomatal resistance and leaf area index (LAI), were measured for the design of the crop energy exchange model. The internal and external environments of the greenhouse such as the solar radiation, the air temperature, and relative humidity were also measured to validate the energy model. The experimentally measured transpiration rate of the crops and the ambient air temperature inside the greenhouse were used to validate the energy exchange model of the crops and the entire greenhouse, respectively. Based on the validated greenhouse energy exchange model, the periodic energy loads which refer to the total energy load at certain periods were calculated according to the crop growth stage. Additionally, the peak energy loads referred to as maximum energy at certain periods were also calculated according to the crop growth stage. Energy loads were calculated using the actual weather data of 10 years (1 June 2009-31 May 2019) for accurate simulations [16]. Energy loads were used to analyse the trends according to the weather conditions and to derive the average energy loads. The hourly energy load data were used to calculate periodic energy loads and to analyse the energy cost and maximum energy loads to decide the proper scale of a heat pump in the greenhouse. The effects of crop energy exchange on energy loads of greenhouses were analysed by comparing the cooling and heating loads with versus without consideration of the crops. Finally, a comparative analysis of the energy costs was conducted according to an energy source as an essential prerequisite for the utilization of the thermal energy of the thermal effluent.

Experimental Greenhouse
The experimental greenhouse was located on the west coast of South Korea in Jugyomyeon, Bo-ryeong City, Chungcheongnam-do Province (126°29′ E, 36°23′ N), which was 1 km south of the Bo-ryeong thermal power plant ( Figure 2). The experimental greenhouse was constructed to enhance the energy efficiency of the heating and cooling system by utilizing a thermal effluent. The experimental greenhouse was an eight-span 1-2-W-type greenhouse, as illustrated in Figure 3. Figure 4 shows thermal screens, crops, vent openings, air ducts that were connected with the heat pump, and circulation fans in the greenhouse. The experimental greenhouse had a width of 34.4 m, a length of 30.0 m, an eave height of 4.5 m, and a ridge height of 5.7 m. The covering of the greenhouse was made of a 0.15 mm thick polyolefin film. The greenhouse was divided into a cultivation space (768 m 2 floor area) for cultivating Irwin mangoes and a workspace for controlling the mechanical system (128 m 2 floor area) (Figure 3). The sides of the experimental greenhouse were covered with double-layer polyolefin film for additional insulation. Thermal screens were installed at the sidewall and ceiling height to decrease conduction heat loss that were caused by outdoor versus indoor air temperature differences at the cladding. After sunset, the thermal screen at the sidewall was drawn vertically, and the thermal screen at the ceiling was installed horizontally at a height of 3.8 m from the floor. The greenhouse was occupied with 100 potted Irwin mango fruit trees. In the experimental greenhouse, the mango trees were pruned to a height of approximately 1.5 m, and the branches were pulled into a globular shape to equalize the light-interception. Therefore, the mango trees in the greenhouse were fixed in shape and size during the experimental periods.

Experimental Greenhouse
The experimental greenhouse was located on the west coast of South Korea in Jugyomyeon, Bo-ryeong City, Chungcheongnam-do Province (126 • 29 E, 36 • 23 N), which was 1 km south of the Bo-ryeong thermal power plant ( Figure 2). The experimental greenhouse was constructed to enhance the energy efficiency of the heating and cooling system by utilizing a thermal effluent. The experimental greenhouse was an eight-span 1-2-W-type greenhouse, as illustrated in Figure 3. Figure 4 shows thermal screens, crops, vent openings, air ducts that were connected with the heat pump, and circulation fans in the greenhouse. The experimental greenhouse had a width of 34.4 m, a length of 30.0 m, an eave height of 4.5 m, and a ridge height of 5.7 m. The covering of the greenhouse was made of a 0.15 mm thick polyolefin film. The greenhouse was divided into a cultivation space (768 m 2 floor area) for cultivating Irwin mangoes and a workspace for controlling the mechanical system (128 m 2 floor area) (Figure 3). The sides of the experimental greenhouse were covered with double-layer polyolefin film for additional insulation. Thermal screens were installed at the sidewall and ceiling height to decrease conduction heat loss that were caused by outdoor versus indoor air temperature differences at the cladding. After sunset, the thermal screen at the sidewall was drawn vertically, and the thermal screen at the ceiling was installed horizontally at a height of 3.8 m from the floor. The greenhouse was occupied with 100 potted Irwin mango fruit trees. In the experimental greenhouse, the mango trees were pruned to a height of approximately 1.5 m, and the branches were pulled into a globular shape to equalize the light-interception. Therefore, the mango trees in the greenhouse were fixed in shape and size during the experimental periods.
The thermal effluent from the Bo-ryeong power plant was used as a heat source to operate the heat pump. The facility discharges thermal effluent at a rate of approximately 3 billion tons per year. The thermal effluent was stored initially in a water storage tank and utilized on the farm and in the greenhouse. The internal environment of the greenhouse was controlled using natural ventilation with side and roof vents and was heated or cooled using three heat pumps (ADF-SLX12WHB, A-San Inc., Gimpo, South Korea). These actuators were controlled automatically based on the internal air temperature and the external solar radiation. The performance of the heat pump was 43,276 W in maximum cooling capacity and 36,786 W in maximum heating capacity. To supply heat energy uniformly, air ducts and circulation fans were installed in the greenhouse. The 16 circulation fans (SGA-300P, Shinan Green Tech Co., Sun-cheon, South Korea) have a capacity of 35 m 3 min −1 per unit and were installed at a height of 2.5 m above the ground floor.

Building Energy Simulation (BES)
Estimation of the energy loads is important for designing the greenhouse and enhancing the efficiency of energy usage. The TRNSYS program (Version 18, Solar Energy Laboratory, University of Wisconsin-Madison, Madison, WI, USA), among several BES programs, was used to model the energy exchange of the greenhouse. TRNSYS program has the advantage of a module-based program, which consists of the main program and several sub-modules to analyse energy flow of each component. TRNSYS also has the advantage of availability and compatibility on an enormous energy system because of a number of sub-modules which can compose various systems such as heat pumps, energy exchange of crop, ventilation. The target building component, which was the experimental greenhouse in this study, consists of several zones. The energy balance equation (Equation (1)) was calculated by zone to simulate the thermal behavior for each zone. It was assumed that the energy exchange by the condensation on the cladding of the greenhouse was relatively small in this study because of the double-layer polyolefin film and thermal screen for reducing heat loss.
where, Q i is the total heat gain of zone I (kJ), Q sur f is the convective heat gain or loss from surfaces (kJ h −1 ), Q in f is the heat gain or loss by infiltration (kJ h −1 ), Q vent is the heat gain or loss by ventilation (kJ h −1 ), Q ishcci is the absorbed solar radiation on all internal shading devices of zone and directly transferred as a convective gain to the internal air (kJ h −1 ), Q solar is the fraction of solar radiation entering a zone (kJ h −1 ), Q g,c is the internal convective gains (kJ h −1 ), and Q cplg is the heat gain or loss due to connective air flow from adjacent zone (kJ h −1 ).  The thermal effluent from the Bo-ryeong power plant was used as a heat source to operate the heat pump. The facility discharges thermal effluent at a rate of approximately 3 billion tons per year. The thermal effluent was stored initially in a water storage tank and utilized on the farm and in the greenhouse. The internal environment of the greenhouse was controlled using natural ventilation with side and roof vents and was heated or cooled using three heat pumps (ADF-SLX12WHB, A-San Inc, South Korea). These actuators were controlled automatically based on the internal air temperature and the external solar radiation. The performance of the heat pump was 43,276 W in maximum cooling capacity and 36,786 W in maximum heating capacity. To supply heat energy uniformly, air ducts and circulation fans were installed in the greenhouse. The 16 circulation fans (SGA-300P, Shinan Green Tech Co., Sun-cheon, South Korea) have a capacity of 35 m 3 min −1 per unit and were installed at a height of 2.5 m above the ground floor.

Building Energy Simulation (BES)
Estimation of the energy loads is important for designing the greenhouse and enhancing the efficiency of energy usage. The TRNSYS program (Version 18, Solar Energy Laboratory, University of Wisconsin-Madison, Madison, WI, USA), among several BES programs, was used to model the energy exchange of the greenhouse. TRNSYS program has the advantage of a module-based program, which consists of the main program and several sub-modules to analyse energy flow of each component. TRNSYS also has the advantage of availability and compatibility on an enormous energy system because of a number of sub-modules which can compose various systems such as heat pumps, energy  The internal environments of the experimental greenhouse were monitored to design and validate the energy exchange model of the greenhouse ( Figure 5). A field experiment was conducted from 28 May to 1 June 2016. To measure the radiation transmissivity of the greenhouse cladding, a pyranometer (SP-110; Apogee Inst., Logan, UT, USA) was installed on the center of the top frame, which was 2.5 m above the bottom floor. The internal air temperature of the greenhouse was measured by a thermocouple (T-type; Ondo114 Co., Seoul, South Korea), which was installed 1.0 m above the ground. Fifteen thermocouples were installed for calculating the average air temperature and three thermocouples were installed at the center of the air blower of each heat pump for monitoring the operating time of the heat pump during the experimental period. The installed sensors were connected to a data logger (GL 820; Graphtech Corp., Irvine, CA, USA), and the data were saved at five-second intervals. air temperature of the greenhouse was measured by a thermocouple (T-type; Ondo114 Co., South Korea), which was installed 1.0 m above the ground. Fifteen thermocouples were installed for calculating the average air temperature and three thermocouples were installed at the center of the air blower of each heat pump for monitoring the operating time of the heat pump during the experimental period. The installed sensors were connected to a data logger (GL 820; Graphtech Corp., Irvine, CA, USA), and the data were saved at five-second intervals.

Measurement of the Stomatal Resistance and the Leaf Area Index of Crops
As the Irwin mango is a perennial plant, it was assumed that the stomatal resistance was constant with growth in this study. The stomatal resistances at the ventral and dorsal sides of the top, middle, and bottom leaves of the crops in the experimental greenhouse were measured at the 17-18 May 2016 using a leaf porometer (SC-1; Decagon devices, USA) ( Figure 6). A leaf porometer is a sensor that measures the vapor conductivity per unit area and calculates the stomatal resistance as the reciprocal of the vapor conductivity.

Measurement of the Stomatal Resistance and the Leaf Area Index of Crops
As the Irwin mango is a perennial plant, it was assumed that the stomatal resistance was constant with growth in this study. The stomatal resistances at the ventral and dorsal sides of the top, middle, and bottom leaves of the crops in the experimental greenhouse were measured at the 17-18 May 2016 using a leaf porometer (SC-1; Decagon devices, Pullman, WA, USA) ( Figure 6). A leaf porometer is a sensor that measures the vapor conductivity per unit area and calculates the stomatal resistance as the reciprocal of the vapor conductivity. The radiation above the plants was measured by a pyranometer, and the surface temperature of the leaf was measured by installing a thermocouple on the ventral side of the leaf. A regression equation for estimating the stomatal resistance from the measured radiation and the leaf temperature was derived and applied as a parameter for an energy exchange model of the crops. The radiation above the plants was measured by a pyranometer, and the surface temperature of the leaf was measured by installing a thermocouple on the ventral side of the leaf. A regression equation for estimating the stomatal resistance from the measured radiation and the leaf temperature was derived and applied as a parameter for an energy exchange model of the crops. LAI is the ratio of the whole leaf area to the ground area that is covered by plants, which represents the leaf distribution and density. LAIs of 100 crops were measured at 17 February 2016 and 10 April 2016 using a plant canopy analyser (LAI-2200; LI-COR Inc., Lincoln, NE, USA) ( Figure 7). LAIs were used to calculate the average LAI of the crops in the greenhouse. Based on the Beer-Lambert law, LAI was calculated from the radiation above and under the canopy, which was measured by a plant canopy analyser. The average LAI of 100 crops was applied in the design of the energy exchange model of the crops. LAI is the ratio of the whole leaf area to the ground area that is covered by plants, which represents the leaf distribution and density. LAIs of 100 crops were measured at 17 February 2016 and 10 April 2016 using a plant canopy analyser (LAI-2200; LI-COR Inc., Lincoln, NE, USA) ( Figure 7). LAIs were used to calculate the average LAI of the crops in the greenhouse. Based on the Beer-Lambert law, LAI was calculated from the radiation above and under the canopy, which was measured by a plant canopy analyser. The average LAI of 100 crops was applied in the design of the energy exchange model of the crops. LAI is the ratio of the whole leaf area to the ground area that is covered by plants, which represents the leaf distribution and density. LAIs of 100 crops were measured at 17 February 2016 and 10 April 2016 using a plant canopy analyser (LAI-2200; LI-COR Inc., Lincoln, NE, USA) ( Figure 7). LAIs were used to calculate the average LAI of the crops in the greenhouse. Based on the Beer-Lambert law, LAI was calculated from the radiation above and under the canopy, which was measured by a plant canopy analyser. The average LAI of 100 crops was applied in the design of the energy exchange model of the crops.

Measurement of Transpiration by Crops
Transpiration is directly related to the latent heat and, therefore, was used as data for validation (Equation (2)). To validate the energy exchange model of the crops, the calculated transpiration rates of the crops were compared with the measured transpiration rates. The ambient air temperature, relative humidity, surface temperature of the leaf, stomatal resistance, and LAI were measured from 28 May to 1 June 2016 and used to calculate the transpiration rate. Figure 8 show the diagram and the arrangement of the experimental equipment for measuring the transpiration rate. The crops were three-year-old Irwin mangoes that were cultivated in a 44.5 cm diameter pot. The pot was placed on an electronic scale (CBX32KH; CAS Corp., South Korea) to measure the change in the weight of the pot over time. The weight of the pot was recorded at five-minute intervals, excluding the wa-

Measurement of Transpiration by Crops
Transpiration is directly related to the latent heat and, therefore, was used as data for validation (Equation (2)). To validate the energy exchange model of the crops, the calculated transpiration rates of the crops were compared with the measured transpiration rates. The ambient air temperature, relative humidity, surface temperature of the leaf, stomatal resistance, and LAI were measured from 28 May to 1 June 2016 and used to calculate the transpiration rate. Figure 8 show the diagram and the arrangement of the experimental equipment for measuring the transpiration rate. The crops were three-yearold Irwin mangoes that were cultivated in a 44.5 cm diameter pot. The pot was placed on an electronic scale (CBX32KH; CAS Corp., Seoul, South Korea) to measure the change in the weight of the pot over time. The weight of the pot was recorded at five-minute intervals, excluding the watering periods. To exclude water evaporation from the soil and to measure only transpiration from the plant, the soil part of the pot was sealed with plastic wrap. A pyranometer was installed on the centre of the top frame and was used to measure the radiation, and 3 thermocouples were installed on the dorsal side of the leaf (top, middle, and bottom) and were used to measure the surface temperature of the leaf. Three HOBO data loggers (UX100-003; Onset Corp., Bourne, MA, USA) were also installed near the plant to measure the ambient air temperature and the relative humidity.
where E is the transpiration per area (kg s −1 m −2 ), LE is the latent heat flux (W m −2 ), and λ is the latent heat of vaporization of water (J kg −1 ).

Design of the Energy Exchange Model of the Greenhouse
A dynamic energy exchange model of the greenhouse was designed by adapting four models: A model of the greenhouse structure, a model for generating the operational signals of the thermal screens, heat pumps, and vent openings, an energy exchange model of the crops, and a model of the thermal effluent-heat pump. Several modules in TRNSYS, such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1).
(UX100-003; Onset Corp., Bourne, MA, USA) were also installed near the plant to measure the ambient air temperature and the relative humidity.
where E is the transpiration per area (kg s −1 m −2 ), LE is the latent heat flux (W m −2 ), and λ is the latent heat of vaporization of water (J kg −1 ).

Design of the Energy Exchange Model of the Greenhouse
A dynamic energy exchange model of the greenhouse was designed by adapting four models: a model of the greenhouse structure, a model for generating the operational signals of the thermal screens, heat pumps, and vent openings, an energy exchange model of the crops, and a model of the thermal effluent-heat pump. Several modules in TRNSYS,  Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon
Module Specification PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power Radiation processor Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power As the energy exchange via radiation and conduction at the cladding depends on the Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power As the energy exchange via radiation and conduction at the cladding depends on the Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration) PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power As the energy exchange via radiation and conduction at the cladding depends on the shape and materials of the greenhouse, it was important to model the greenhouse clad-

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power As the energy exchange via radiation and conduction at the cladding depends on the shape and materials of the greenhouse, it was important to model the greenhouse clad-

Switch
Determines on/off signals for each device based on the external/internal weather conditions PEER REVIEW 10 of 28 such as the data reader module for reading weather data as input data and the psychrometric module for calculating the properties of moist air, were used for each part of the model (Table 1). Table 1. Description of the modules in TRNSYS such as data reader, radiation processor, psychrometric calculator, sky temperature calculator, multi-zone, user-defined function, switch, and heat pump for modeling the experimental greenhouse.

Icon Module Specification
Data reader Reads weather and sensor data to send input data to other modules

Radiation processor
Interpolates radiation data, calculates several quantities that are related to the position of the sun and estimates the insolation on surfaces of either fixed or variable orientation Psychrometric calculator Calculates moist air by accepting as input the dry bulb temperature and the relative humidity Sky temperature calculator Determines the effective sky temperature for calculating the longwave radiation exchange between an arbitrary external surface and the atmosphere Multi-zone (greenhouse) Models the thermal behavior inside a greenhouse (heating, cooling, ventilation, and infiltration)

User-defined function
Calculates average greenhouse internal temperature as input data for the switch module, the sensible/latent heat flux by crops, and the energy production and power consumption of the heat pump Switch Determines on/off signals for each device based on the external/internal weather conditions Heat pump Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power As the energy exchange via radiation and conduction at the cladding depends on the shape and materials of the greenhouse, it was important to model the greenhouse cladding precisely. The multi-zone module was used to design the structure of the greenhouse

Heat pump
Models a single-stage liquid-source heat pump based on user-supplied data files that contain catalogue data for the capacity and power As the energy exchange via radiation and conduction at the cladding depends on the shape and materials of the greenhouse, it was important to model the greenhouse cladding precisely. The multi-zone module was used to design the structure of the greenhouse according to geometric information on the experimental greenhouse (Figure 3). The structure data and cladding data for the greenhouse model were applied to the multi-zone module using the TRNSYS plugged-in SketchUp program (ver. 8, Google, Mountain View, CA, USA) and the WINDOW program (ver. 7.4, LBNL, Berkeley, CA, USA), respectively. As the multi-zone module cannot create the internal wall in a single zone, the crop growth space was divided horizontally to set the thermal screen in the greenhouse. The physical properties of the 0.15 mm polyolefin film that was used for the cladding of the greenhouse were specified as follows: thermal conductivity of 0.330 W m −1 K −1 , solar transmittance of 0.797, solar reflectance of 0.106, visible transmittance of 0.935, visible reflectance of 0.106, and infrared emittance of 0.840. The physical properties of the materials that were used for the frame and floor of the greenhouse were as listed in Table 2. To set the thermal screen, each span was divided into three zones; therefore, there were twenty four zones for the crop growth space. A virtual wall was set up between each span to simulate thermal exchange via convection and between the upper sides of the cultivation space to configure the thermal screen on the working schedule. The data reader, radiation processor, psychometric, sky temperature calculator and multi-zone module were used to simulate the energy exchange at the cladding in TRNSYS. The solar radiation on each wall, the dew-point temperature, and the sky temperature were calculated at each module using weather data. These calculated data were linked to the multi-zone module for the simulation of the energy exchange via conduction and radiation (Figure 9a).
Inside the experimental greenhouse, the heat pumps, thermal screens and vent openings were operated according to the internal environments. The operating condition for all equipment was measured by a field experiment. The operating conditions were used as a boundary condition into the greenhouse BES model. Switch modules were used to design the thermal screens, heat pumps and vent openings in the greenhouse (Figure 9b). The internal air temperature, the amount of external solar radiation, and the external air temperature were used as input data to generate the operating signal of each equipment. The operating signal can take a value of 0 or 1, and if the value was 1, the equipment operated. The operating conditions of each equipment were determined via Equations (3)-(5), and the operating temperature was established based on the optimum temperature of Irwin mango crops according to the growth period.
where TSSIGN is the thermal screen signal, HEATSIGN is the heating signal, HEATTEMP is the set temperature of heating, COOLSIGN is the cooling signal, COOLTEMP is the set temperature of cooling, VENTSIGN is the ventilation signal, VENTTEMP is the set temperature of ventilation, T in is the internal air temperature ( • C), and T out is the external air temperature ( • C). In this study, the sensible and latent heat of crops were used to simulate the energy exchange of crops. The sensible and latent heat of crops were calculated via Equations (6) and (7). The BES modules for the energy exchange of crops are illustrated in Figure 9c. The environmental conditions and characteristics of the crops were used to model the energy exchange model of the crops. The external resistance of the crops (r e ) was determined in real time based on Equation (8). The characteristic leaf dimension (l) was measured for three crops to represent each growth stage of the Irwin mangoes in the greenhouse, and the average value was calculated as 9.99 cm. As the average measured air velocity in the greenhouse was less than 0.1 m s −1 , the wind speed (u) was set to 0.1 m s −1 , which is identical to the value that was used by Stanghellini [17]. The regression model that was derived from the data from the field experiments was used to calculate the stomatal resistance (r s ). The measured average LAI of the cultivated crops of 1.78 was used as a constant under the assumption that the crops were pruned during the experimental period. The thermal energy that was stored in the plants (G) was highly variable according to the radiation. In this study, the thermal energy that was stored in the plants was calculated as 15% of the net radiation flux [18,19].
a − e a ) ÷ r e γ * = γ 1 + r s r e (7) where H is the sensible heat of the crop (W m −2 ), LE is the latent heat of the crop (W m −2 ), ρ a is the density of air (kg m −3 ), c p is the specific heat of air at constant pressure (J kg −1 K −1 ), T L is the temperature at the surface of the leaf (K), T a is the ambient air temperature (K), r e is the external resistance of the surface of the leaf to sensible heat transfer (s m −1 ), δ is the slope of the saturation vapor pressure-temperature curve (Pa K −1 ), γ is the thermodynamic psychometric constant (Pa K −1 ), r s is the stomatal resistance (s m −1 ), R G is the solar radiation (W m −2 ), G is the thermal energy that is stored in the crops (W m −2 ), LAI is the leaf area index (dimensionless), e * a is the saturation vapour pressure at the air temperature (Pa), e a is the vapour pressure of the air (Pa), l is the characteristic leaf dimension (m), and u is the wind speed (m s −1 ).
When calculating the energy loads of the greenhouse using historical weather data, a limitation was encountered in simulating the energy exchange of the crop, which was due the absence of real-time environmental data inside the greenhouse. Therefore, the ratio of the sensible and latent heats to the external radiation was used to design the energy exchange model of the greenhouse based on research by [10]. The ratio of the sensible and latent heats from crops to the external radiation was derived from the validated energy exchange model of the crop in this study. Agronomy 2021, 11, x FOR PEER REVIEW 12 of 28 The thermal effluent-heat pump system was modelled based on performance data of the heat pump. Real-time thermal effluent data were used to calculate the rated cooling and heating capacities and the power consumption. The greenhouse internal air temperature and humidity level were considered simultaneously in the calculation of the actual energy production and power consumption levels. The actual calculated energy production and power consumption levels were compared with the energy loads of the experimental greenhouse to simulate the energy consumption of the experimental greenhouse. The BES modules for the thermal effluent-heat pump system are illustrated in Figure 9d. The rated air flow and water flow were set to 2405 and 2.11 l s −1 , respectively, and the cooling capacity, heating capacity, and power consumption were set for several inlet water temperatures. The cooling and heating performance data at inlet water temperatures of 1, 5, 10, 15, 20, 25, and 30 • C that were provided by the company were used ( Table 3). The performance correction data were additionally applied to consider changes in the performance and the power consumption of the heat pump that were due to the changes in the greenhouse internal air temperature and humidity. The dynamic energy exchange model of the experimental greenhouse was designed based on the BES technique, which includes the BES models for the greenhouse structure, generating the operational signals, the energy exchange model of crops, and the thermal effluent heat pump system (Figure 9e). Each part was connected to the greenhouse to simultaneously simulate the energy flow. The simulations of the energy exchange model were conducted at 1 min intervals.

Accuracy Evaluation Method for the BES Model
The internal air temperature computed by the designed BES model was compared with internal air temperature measured at the field experiment to evaluate the accuracy of the greenhouse BES model. The simulated transpiration rate was also compared with the measured transpiration rate to evaluate the energy exchange model of the crop. The measured crop characteristics and environmental conditions, such as the LAI, ambient air temperature, humidity, and radiation above the plants, were used to simulate the transpiration rate.
In this study, the coefficient of determination (R 2 ) was used to analyse the tendency. The index of agreement (d) was used to assess the error between the simulated data and the measured data [20]. R 2 is a coefficient that only indicates the tendency between two groups of data and not the direction or errors. R 2 can take a value that is between '0' and '1'. The more similar the tendencies of two groups of data are, the closer R 2 is to 1 (Equation (9)). The index of agreement can standardize the measure of the degree of the model prediction error, which can take a value between '0' and '1'. The closer the index is to 1, the more closely the simulated data correspond to the measured data (Equation (10)).
where Y obs i is the i-th measured value, Y obs is the average of the measured values, Y si i is the i-th simulation value, and Y si is the average of the simulation values. STDEV obs is the standard deviation of the measured values.

Analysis of the Greenhouse Energy Loads
The experimental greenhouse maintains the growing environment in a suitable range of air temperatures for growing the crops. However, to bear fruit, the crops must undergo several growing stages, which have a diverse range of growing temperatures. The air temperatures for heating, cooling and ventilation were set based on the optimum growing temperature for each growth stage (Table 4). Weather data from 2009 to 2019 collected at the Seo-san weather station of the Korean Meteorological Administration (KMA), which is the nearest weather station from the experimental greenhouse, were used to calculate the energy load of the greenhouse. Weather data such as the air temperature, the relative humidity, the wind environment, and the solar radiation were observed at intervals of 1 h. The periodic and maximum energy loads were calculated by summarizing and comparing the hourly energy loads for each year. Periodic energy loads were used to analyse the energy consumptions of the experimental greenhouse according to the growth stages. Maximum energy loads were used to assess the capacities of heating and cooling systems installed at the experimental greenhouse and suggest the proper capacities of cooling and heating systems. Based on the growth periods, the analysis period was from May 30 to June 1 for each year. As the Bo-ryeong thermal power plant only provided five years of discharge flow data, the temperature of the thermal effluent and the daily oil price for the five-year period (2010 to 2014) were used to compare the energy cost, which depends on the energy source. The internal air temperatures by location and time were analysed according to day (8:00 AM to 8:00 PM) and night (8:00 PM to 8:00 AM) ( Table 5). The average air temperature of the entire period was 29.9 • C, and the standard deviation of each location was 0.57 • C, which was low. The average air temperatures and standard deviations for each location were 33.6 and 0.92 • C during day, and 25.4 and 0.30 • C during night, respectively. The standard deviation of each location at night was lower than that during the day because air ducts and circulation fans were installed to ensure that the thermal environment of the crop growth room remains uniform.  Figure 10 shows a graph of the average air temperature at the centres of three heat pump air blowers and the average internal air temperatures of the greenhouse. According to the graph, the set temperatures of the heat pump on 29-31 May were 25 • C, and the air temperature of an air blast was approximately 42 • C when the heat pump was operated. The air temperatures were changed by the operation of the heat pumps. The heat pump did not operate on 28 May because of a mechanical defect, and the set temperature of the heat pump was increased by approximately 2-3 • C on 1 June. Figure 11 is a graph comparing the internal and external solar radiation of the experimental greenhouse. The trend of the internal solar radiation of the greenhouse usually followed that of the external radiation; however, the internal solar radiation decreased at approximately 2:00 PM every day due to shade from a folded thermal screen. The average radiation transmissivity of the cladding, except for the data at approximately 2:00 PM, was analysed and found to be 83%.

Stomatal Resistance and the Leaf Area Index of the Crops
The field experiment for measuring the stomatal resistance of the crops in the greenhouse was conducted. The stomatal resistances and the leaf temperatures at the ventral and dorsal sides of the top, middle, and bottom leaves were measured. The average stomatal resistances on the two sides of the top, middle, and bottom leaves were analysed, and the average value was found to differ according to the side of the leaf. As the Irwin mango is a perennial plant, the stomata are distributed on the ventral side of the leaf and are scarce on the dorsal leaf. For this reason, there were significantly fewer stomata on the dorsal side of the leaf than on the ventral side. The average stomatal resistance on the dorsal side (top leaf: 2023 s m −1 , middle leaf: 1820 s m −1 , and bottom leaf: 1906 s m −1 ) was analysed and found to be much higher than the average stomatal resistance on the ventral side (top leaf: 850 s m −1 , middle leaf: 168 s m −1 , and bottom leaf: 264 s m −1 ).  Figure 11 is a graph comparing the internal and external solar radiation of the experimental greenhouse. The trend of the internal solar radiation of the greenhouse usually followed that of the external radiation; however, the internal solar radiation decreased at approximately 2:00 P.M. every day due to shade from a folded thermal screen. The average radiation transmissivity of the cladding, except for the data at approximately 2:00 P.M., was analysed and found to be 83%.    Figure 11 is a graph comparing the internal and external solar radiation of the experimental greenhouse. The trend of the internal solar radiation of the greenhouse usually followed that of the external radiation; however, the internal solar radiation decreased at approximately 2:00 P.M. every day due to shade from a folded thermal screen. The average radiation transmissivity of the cladding, except for the data at approximately 2:00 P.M., was analysed and found to be 83%.  A non-linear regression equation for calculating the stomatal resistance from the radiation and the leaf temperature was derived and used as a parameter for the energy exchange model of the crops (Equation (11)). The derived equation was found to simulate the actual stomatal resistance precisely (R 2 = 0.854 and d = 0.917).
To calculate the average LAI of the Irwin mango in the greenhouse, a plant canopy analyser was used to measure the LAIs of 100 Irwin mango pots. As the average LAI of 100 crops was 1.78, the average value of 1.78 was used in the design of the energy exchange model of the crops.

Validation of Transpiration of Crops
To measure the transpiration rate of the crops, the weight of the pot in each case was recorded at five-minute intervals. The results of measuring the transpiration are shown in Figure 12. During the daytime, the crops transpired actively, whereas during the nighttime (8:00 PM to 5:00 AM), the average measured transpiration rate was 0.535 g crop −1 . However, on the night of 28 May, the average transpiration rate was 0.154 g crop −1 , which was 70% lower than those during all night-time periods of the experimental period. This result stemmed from the water content of the plants, which rose after the May 28, as irrigation provided half of a litre of water typically at 9:00 AM, while one litre of water was provided by irrigation every evening at 9:00 PM during the experimental period. For the same reason, the average transpiration rate during the daytime (9:00 AM to 6:00 PM) on 28 May (7.212 g crop −1 ) was found to be 15% lower than all daytime amounts during the experimental period (8.427 g crop −1 ).

Validation of the Air Temperature inside the Greenhouse
To validate the greenhouse dynamic energy exchange model, the measured greenhouse internal air temperature was compared with the simulated data. The hourly internal air temperatures were calculated based on five-second intervals of measured data from fifteen locations and one-minute intervals of simulated data and compared ( Figure  13). The average values of the internal air temperature at all locations were represented by dots in the figure, and the standard deviations were represented by error bars. The simulated data were the average air temperature of the crop growth space overall. To validate the energy exchange model of the crops, the transpiration rate that was measured by the field experiment was compared with the simulated data. Measured data at five-minute intervals and simulated data at one-minute intervals were compared ( Figure 12). The designed energy exchange model of the crops was validated as providing a satisfactory estimate of the change in the transpiration rate throughout the experimental period (R 2 = 0.96, d = 0.98). Considering the error on 28 May, with a dataset of four days that does not include the first day, the energy exchange model of the crops that was designed in this study was found to reflect the actual latent heat from the crops accurately.

Validation of the Air Temperature inside the Greenhouse
To validate the greenhouse dynamic energy exchange model, the measured greenhouse internal air temperature was compared with the simulated data. The hourly internal air temperatures were calculated based on five-second intervals of measured data from fifteen locations and one-minute intervals of simulated data and compared ( Figure 13). The average values of the internal air temperature at all locations were represented by dots in the figure, and the standard deviations were represented by error bars. The simulated data were the average air temperature of the crop growth space overall.

Validation of the Air Temperature inside the Greenhouse
To validate the greenhouse dynamic energy exchange model, the measured greenhouse internal air temperature was compared with the simulated data. The hourly internal air temperatures were calculated based on five-second intervals of measured data from fifteen locations and one-minute intervals of simulated data and compared ( Figure  13). The average values of the internal air temperature at all locations were represented by dots in the figure, and the standard deviations were represented by error bars. The simulated data were the average air temperature of the crop growth space overall. In the daytime, the internal air temperature increased due to an increase in the solar radiation, while at night, the heat pump operated based on a set temperature to maintain the internal air temperature. The simulated data typically follow the measured data; however, on the nights of the 27 and 31 May 2016, the internal air temperature differed substantially from those on other days. According to the data logger, on the night of 27 May, the heat pump did not operate, while on the night of 31 May, the set temperature of the heat pump changed from 25 °C to 27 °C. Therefore, the working schedule of the heat pump was modified, and the designed greenhouse dynamic energy exchange model was shown to provide a satisfactory estimate of the change of the internal air temperature throughout the experimental period (R 2 = 0.94 and d = 0.97) (Figure 14). In the daytime, the internal air temperature increased due to an increase in the solar radiation, while at night, the heat pump operated based on a set temperature to maintain the internal air temperature. The simulated data typically follow the measured data; however, on the nights of the 27 and 31 May 2016, the internal air temperature differed substantially from those on other days. According to the data logger, on the night of 27 May, the heat pump did not operate, while on the night of 31 May, the set temperature of the heat pump changed from 25 • C to 27 • C. Therefore, the working schedule of the heat pump was modified, and the designed greenhouse dynamic energy exchange model was shown to provide a satisfactory estimate of the change of the internal air temperature throughout the experimental period (R 2 = 0.94 and d = 0.97) ( Figure 14).

Validation of the Air Temperature inside the Greenhouse
The simulated sensible and latent heats of crops that were obtained using the validated energy exchange model of the crop are plotted in Figure 15. Positive values of sensible and latent heats of the crops indicate energy flow from the ambient air to the crops, and negative values indicate flow in the opposite direction. As the solar radiation increased in the daytime, the energy exchange by the crops also increased, with energy flowing entirely from the air to the crops. In contrast, in the night-time, the energy exchange by the crops dropped dramatically due to the decreased amount of radiation, with the energy flow mainly due to the vapour pressure deficit between the crops and the surrounding air. Moreover, the surface temperatures of the leaves were higher than the ambient air temperature; therefore, sensible heat flux occurs from the crops to the air. According to the results of simulating the sensible heat and latent heat, the crops absorbed thermal energy at proportions of 67.0% of the latent heat and 5.3% of the sensible heat of the external radiation in the daytime. Other proportions of the external radiation were due to the photosynthesis rate stored in the plants, thermal energy rate stored in the plant, the wide spacing between pots of crops. At night-time, the crops absorbed a thermal en-

Validation of the Air Temperature inside the Greenhouse
The simulated sensible and latent heats of crops that were obtained using the validated energy exchange model of the crop are plotted in Figure 15. Positive values of sensible and latent heats of the crops indicate energy flow from the ambient air to the crops, and negative values indicate flow in the opposite direction. As the solar radiation increased in the daytime, the energy exchange by the crops also increased, with energy flowing entirely from the air to the crops. In contrast, in the night-time, the energy exchange by the crops dropped dramatically due to the decreased amount of radiation, with the energy flow mainly due to the vapour pressure deficit between the crops and the surrounding air. Moreover, the surface temperatures of the leaves were higher than the ambient air temperature; therefore, sensible heat flux occurs from the crops to the air. According to the results of simulating the sensible heat and latent heat, the crops absorbed thermal energy at proportions of 67.0% of the latent heat and 5.3% of the sensible heat of the external radiation in the daytime. Other proportions of the external radiation were due to the photosynthesis rate stored in the plants, thermal energy rate stored in the plant, the wide spacing between pots of crops. At night-time, the crops absorbed a thermal energy of 29.54 W m −2 as latent heat and emitted a thermal energy of 9.49 W m −2 as sensible heat. This ratio for the daytime and the values of the latent heat and the sensible heat were used to analyse the energy load of the experimental greenhouse. creased in the daytime, the energy exchange by the crops also increased, with energy flowing entirely from the air to the crops. In contrast, in the night-time, the energy exchange by the crops dropped dramatically due to the decreased amount of radiation, with the energy flow mainly due to the vapour pressure deficit between the crops and the surrounding air. Moreover, the surface temperatures of the leaves were higher than the ambient air temperature; therefore, sensible heat flux occurs from the crops to the air. According to the results of simulating the sensible heat and latent heat, the crops absorbed thermal energy at proportions of 67.0% of the latent heat and 5.3% of the sensible heat of the external radiation in the daytime. Other proportions of the external radiation were due to the photosynthesis rate stored in the plants, thermal energy rate stored in the plant, the wide spacing between pots of crops. At night-time, the crops absorbed a thermal energy of 29.54 W m −2 as latent heat and emitted a thermal energy of 9.49 W m −2 as sensible heat. This ratio for the daytime and the values of the latent heat and the sensible heat were used to analyse the energy load of the experimental greenhouse.  The periodic energy loads over a period of ten years were calculated. Table 6 presents the calculation results of the periodic energy loads and the external weather conditions according to the growing stage in 1 June 2014-31 May 2015 as a representative case. The cooling and heating loads mainly depended on the external air temperature and the solar radiation. The total energy load per day of the experimental greenhouse was the highest in the fruit-bearing period in winter, whereas the total energy load per day of the experimental greenhouse was the lowest in the generative growth period. The total heating load was estimated to be about 3.1 times larger than the total cooling load over the whole growth period of 1 June 2014-31 May 2015. The cooling loads were relatively high compared to the heating loads in the generative growth period and the floral-initiation period. The heating loads were relatively high compared to the cooling loads after the floral-initiation period. As the generative growth period and the floral- initiation period of the Irwin mango start in summer, heating loads seldom occurred in these two periods. Moreover, the set temperature of the heat pump for cooling was changed from 30 • C to 20 • C as the growth period progressed; therefore, the periodic cooling loads increased from 73,693 MJ to 158,949 MJ. As the season changes from summer to autumn, the cooling loads decrease and the heating loads increase compared to the energy loads after the floral-initiation period. In the flowering period, the heating loads were higher than the cooling loads, while in the fruit-bearing period, no cooling loads were observed. For the periodic heating loads, the set temperature of the heat pump for heating was changed from 18 • C to 22 • C; hence, the periodic heating loads increased from 231,717 MJ in the flowering period to 248,483 MJ in the fruit-bearing period, which represents an increase of approximately 6%. In the fruit-growing period, the internal air temperature of the greenhouse should be maintained at more than 25 • C in the winter and spring seasons. Therefore, higher heating loads were observed relative to the cooling loads. Table 7 lists the maximum cooling and heating loads and the times of their occurrence for each the growth stage in 2014-2015 as a representative case. The maximum cooling loads typically occurred at approximately 2:00 PM to 4:00 PM, when the external air temperature and the solar radiation were high during the day, whereas the maximum heating loads typically occurred at approximately 5:00 AM to 7:00 AM, when the external air temperature and the solar radiation level reached their lowest points of the day. As the internal air temperature of the greenhouse should be maintained in the range of 8-20 • C in spite of the summer season in the floral-initiation period, the cooling load was 456,047 kJ h −1 , which was the highest among the growth stages. Crops in greenhouses mainly absorb thermal energy in the daytime and release thermal energy at night. Many researchers who study energy balances in greenhouses emphasize that a large portion of radiation that penetrates the cladding is used for transpiration by crops [10,[21][22][23][24]. The energy loads of a greenhouse can be defined as the required energy for maintaining the internal air temperature against heat transfer via conduction, radiation, and convection. If crops in the greenhouse were assumed that not to exchange energy with the ambient air, this energy will change the energy loads of the greenhouse. Therefore, to calculate the energy loads of greenhouses more accurately, the energy flow from the crops should be considered precisely. However, as stated in the literature review, most studies that simulate energy balances in greenhouses and calculate the energy loads of greenhouses do not consider the energy exchange by crops or do not consider the characteristics of cultivated crops and instead calculate the energy exchange by crops based on the findings of earlier research. In the present study, the characteristics of Irwin mangoes, such as the LAI, the stomatal resistance, and the shape of the leaves, were considered in the design of the energy exchange model of these crops. Therefore, the periodic and maximum energy loads were compared in consideration of the crops to determine the effect of the crops on the energy loads of the greenhouse. The annual and the maximum energy loads of the experimental greenhouse in consideration of the presence of the crops are plotted in Figures 16 and 17. The annual and the maximum cooling loads were increased by approximately 19% and 12%, respectively, on average when assuming that the crops were not in cultivation. The cultivated crops of the greenhouse BES model that was used in this study were modelled to absorb thermal energy in proportion to the amount of external radiation during the daytime and with a constant value at night. As the cooling loads typically occurred during the daytime, the energy that was absorbed by the crops reduced the cooling loads. For the same reason, the energy that was absorbed by the crops at night increased the heating loads; therefore, the corresponding annual and maximum heating loads were decreased by approximately 5% and 1%, respectively, on average when assuming that the crops were not in cultivation.

Evaluation of the Design Capacities of the Heating and Cooling Systems
The maximum energy loads were calculated and used to evaluate the capacities of the heat pumps that were installed in the experimental greenhouse and to estimate the proper capacities of the cooling and heating systems (Table 8). In the last five years, the maximum cooling load was 590,258 kJ h −1 , which occurred on 16 August 2019 at 3:00 P.M., and the maximum heating load was 471,388 kJ h −1 , which occurred on 18 February 2017 at 8:00 A.M. According to the Rural Development Administration of Korea, the capacity of a heat pump in a greenhouse should be set within 70% of the maximum energy load of the last five years [25]. Therefore, the proper performance of the heat pump in the experimental greenhouse is as follows: approximately 413,181 kJ h −1 for the cooling capacity and 329,972 kJ h −1 for the heating capacity. However, the performance of the heat pump in the experimental greenhouse corresponded to a maximum cooling capacity of 467,381 kJ h −1 and a maximum heating capacity of 397,267 kJ h −1 . The design of the heat pump in the experimental greenhouse was excessive by approximately 54,200 kJ h −1 (13.1% of the proper cooling capacity) in terms of the cooling capacity and by approximately 67,295 kJ h −1 (20.4% of the proper heating capacity) in terms of the heating capacity. This was due

Evaluation of the Design Capacities of the Heating and Cooling Systems
The maximum energy loads were calculated and used to evaluate the capacities of the heat pumps that were installed in the experimental greenhouse and to estimate the proper capacities of the cooling and heating systems (Table 8). In the last five years, the maximum cooling load was 590,258 kJ h −1 , which occurred on 16 August 2019 at 3:00 PM, and the maximum heating load was 471,388 kJ h −1 , which occurred on 18 February 2017 at 8:00 AM According to the Rural Development Administration of Korea, the capacity of a heat pump in a greenhouse should be set within 70% of the maximum energy load of the last five years [25]. Therefore, the proper performance of the heat pump in the experimental greenhouse is as follows: approximately 413,181 kJ h −1 for the cooling capacity and 329,972 kJ h −1 for the heating capacity. However, the performance of the heat pump in the experimental greenhouse corresponded to a maximum cooling capacity of 467,381 kJ h −1 and a maximum heating capacity of 397,267 kJ h −1 . The design of the heat pump in the experimental greenhouse was excessive by approximately 54,200 kJ h −1 (13.1% of the proper cooling capacity) in terms of the cooling capacity and by approximately 67,295 kJ h −1 (20.4% of the proper heating capacity) in terms of the heating capacity. This was due to the method that was used to calculate the energy load of the greenhouse. The calculation method that was provided by the Nigerian Institution of Agricultural Engineers (NIAE) was based on a static design method and does not consider the energy exchange by crops in the greenhouse. Lee [26] determined that the maximum heating loads that were calculated via a static design method were overestimated by 37-55% compared to the maximum heating loads that were calculated via a dynamic design method. In addition, Lee [26] did not consider the crops or the thermal screen in the greenhouse; therefore, the difference between the proper performance of the heat pump as calculated via the greenhouse BES model in this study and the performance of the heat pump in the experimental greenhouse was explainable. However, the difference in the heat pump performance levels likely causes an economic loss. Therefore, it was important to calculate the energy loads of greenhouses using a dynamic design method for the realization of a more feasible heat pump design.

Comparative Analysis of the Energy Cost According to the Energy Source
The energy costs were compared according to the type of energy source, which was kerosene in this case, for the thermal effluent-heat pump system based on the calculated annual energy loads. The analysis period was the period of 2010-2014, when the temperature data of the thermal effluent were logged at the experimental greenhouse. The periodic energy loads, as calculated above, were analysed from 1 June to 31 May based on the growth stage of the Irwin mangoes; however, the energy costs were analysed from 1 January to 31 December based on the thermal effluent data. The boiler was set to use kerosene, which is the most common fuel in the greenhouses of South Korea. The caloric value of the kerosene boiler was set to 9200 kcal L −1 , and the average daily oil price was used to calculate the energy cost. The price of electricity for operating the heat pump was set to 41.9 KRW kWh −1 the same for five years based on the agricultural electricity price. The annual energy usage and cost are presented in Table 9. The average energy cost when using the thermal effluent-heat pump system was found to be 68.21% lower than that when a kerosene boiler was used. Table 10 presents the average of the annual average oil price, the sea water temperature, and the thermal effluent temperature. The effect of the temperature of the sea water on the thermal effluent was minor; however, the oil price has a substantial influence on the energy savings. In 2010, when the oil price was relatively low, the energy consumption reduction ratio when using the thermal effluent-heat pump system was 58.79%, namely, more than 50%. The instability of greenhouses in South Korea, which are fundamentally vulnerable to changes in world oil prices, can be determined with the application of the thermal effluent-heat pump system based on these results.

Conclusions
In this study, an energy exchange model of the experimental greenhouse was designed in four parts: a model of the greenhouse structure, a model for generating the operational signals of the thermal screens, heat pumps, and vent openings, an energy exchange model of the crops, and a model of the thermal effluent-heat pump. To design the energy exchange model, field experiments were conducted to measure the structural characteristics of the greenhouse and the crop characteristics and the working schedule of the thermal screens, heat pumps, and vent openings in the greenhouse. The micro-climate of the experimental greenhouse and the transpiration rate of the cultivated crops were measured to validate the accuracy of the model. The validated model was used to analyse the periodic and the maximum energy loads. The energy loads of the experimental greenhouse were analysed according to the growth stage of the crops in consideration of the crop characteristics. The annual energy loads were used to analyse the energy cost of using the thermal effluent-heat pump system.
The designed energy exchange model of the crops was validated via comparison of the transpiration rate, which provides a satisfactory estimate of the change of the transpiration rate, throughout the experimental period (R 2 = 0.96 and d = 0.98). The energy exchange model of the experimental greenhouse was also validated via comparison with the internal air temperature of the experimental greenhouse that was measured in field experiments. The designed model was found to be able to provide satisfactory estimates of the changes of the internal air temperature of the greenhouse (R 2 = 0.94 and d = 0.97).
The differences in the greenhouse energy loads of the experimental greenhouse according to the presence of the crops were analysed. The maximum cooling loads were overestimated by approximately 19%, and the maximum heating loads were underestimated by approximately 5% when assuming that crops in the greenhouse did not exchange energy with the ambient air. The proper performance of the heat pump was calculated to be 363,092 kJ h −1 in terms of the cooling capacity and 328,910 kJ h −1 in terms of the heating capacity. As the performances of the heat pumps that were installed at the experimental greenhouse were 413,181 kJ h −1 for the maximum cooling capacity and 329,972 kJ h −1 for the maximum heating capacity, the capacities of the heat pumps were excessive by approximately 13.1% and 20.4%, respectively. Finally, the energy costs were compared according to the type of energy source based on the calculated annual energy loads. The