Analysis of Energy Exchange with the Ground in a Two-Chamber Vegetable Cold Store, Assuming Di ﬀ erent Lengths of Technological Break, with the Use of a Numerical Calculation Method—A Case Study

: The paper deals with the impact of the technological break duration during the cold storage cycle on the energy demand of the cold store for vegetables and fruit and the temperature distribution in the ground under the cold store. The studied facility was a two-chamber vegetable cold store located in southern Poland used to store carrots ( Daucus carota ) for nine months a year. The experiments were conducted for 12 months (1 May 2017–30 April 2018). The technological break during this period lasted three months (from 1 July 2018 to 30 September 2018). Continuous measurements (with 1-h frequency) were made in order to determine the boundary conditions for numerical analysis. The measured parameters included indoor air temperature, outdoor air temperature, ground temperature under the building and in its vicinity. There were 22 measuring points andPT100 sensors were used. The numerical analysis was based on the elementary balances method. WUFIplus ® software was used as a calculation supporting tool. The numerical analysis was conducted for 14 calculation variants, with di ﬀ erent duration of technological break. The calculation model validation was performed and the results showed a good correlation with the experimental data. The results of experimental studies and of calculations showed a signiﬁcant impact of the technological break duration on the soil distribution in the ground and the building energy demand. A technological break of less than 4 weeks is the most optimal in the summer. The technological break longer than 4 weeks signiﬁcantly a ﬀ ects the cooling energy demand in the ﬁrst days of the cooling cycle and signiﬁcantly extends the time necessary for the ground and the ﬂoor to reach the optimum temperature. The analysis of the ﬂoor temperature results (points A1–C1) showed that the technological break longer than four weeks causes the average ﬂoor temperature to exceed 4.0 ◦ C. Therefore, the optimum solution is technological break lasting 7–35 days. Absence of technological break results in a decrease of energy gains from the ground by 20% relative to a three-month technological break. The impact of technological break duration was clearly seen in terms of energy losses from the cold store to the ground. In case of a 91-day technological break, the energy losses to the ground were 1289.5 kWh / a, while in case of absence of technological break this value was ninefold lower (147.5 kWh / a).


Introduction
In order to preserve the suitable commercial nutritional and biological quality, the vegetables and fruit should be stored at a specific temperature and humidity, with appropriate gas composition and

Materials and Methods
The experimental studies were conducted in a vegetable cold store in southern Poland. The building has two cold chambers totalling 170 m 2 (chamber K1-78 m 2 , chamber K2-96 m 2 ) and a 58 m 2 sorting and packing area. The walls are made of steel H-sections, with sandwich panels between them. The insulating core in the panels is 15 cm XPS. The roof is also steel structure and sandwich panels as in case of the walls. The floor on the ground was made in traditional technology (10 cm cement screed, 4 cm XPS for thermal insulation). The dimensions of cold chambers were adapted to pallet boxes in which carrots (Daucus carota) were stored. During the cooling cycle (from 15 October to 30 June), the indoor air temperature in the chambers was kept at 0.0-4.0 • C, in accordance with Energies 2020, 13, 4970 3 of 15 recommendations for the stored crops. The cooling system was based on a 20 kW refrigerating unit. An organoleptic soil examination showed that the cold store surrounding includes a 40-cm thick humus (fertile soil) layer, with silty clay below it down to the depth of 2.0 m below ground level. The examination was made using a 2.50 m long geotechnical auger.
The scope of experimental studies included continuous measurements of indoor air temperature, outdoor air temperature, floor temperature, ground temperature under the floor down to the depth of 1.50 m below the floor level, ground temperature in the cold store surrounding down to the depth of 1.50 m below ground level. Twenty-two PT100 sensors with 0.1 • C measurement accuracy were used. These parameters were measured for the chamber K2. The preliminary tests by the authors who analyzed the temperature and humidity conditions in both chambers did not show any statistically significant differences, so it was decided to continue with more accurate studies for the larger chamber (K2). The measurements were automatically recorded at 60-minintervals on a computer hard drive connected to the measuring instruments. The location of measurement points and the cold store layout are presented in Figure 1. The results of experimental studies allowed determining the thermal conditions inside the examined chamber and in the ground below it. The unit heat losses to the ground through the floor were estimated based on the relationship (1), where: Θ i -indoor air temperature, ( • C), Θ g -ground temperature interpolated linearly between points (A1-C1), ( • C), λ s -thermal conductivity coefficient of soil, the used value was 1.9 W/mK, (EN ISO 6946), R si -heat penetration resistance, the used value was 0.14 m 2 K/W. The calculations were made based on the floor temperature measured 0.05 m below its surface and the actual curve of indoor air temperature over time was accounted for. Then, the unit values (W/m 2 ) were added and presented in kWh.
Energies 2020, 13, x FOR PEER REVIEW 3 of 15 pallet boxes in which carrots (Daucus carota) were stored. During the cooling cycle (from 15 October to 30 June), the indoor air temperature in the chambers was kept at 0.0-4.0 °C, in accordance with recommendations for the stored crops. The cooling system was based on a 20 kW refrigerating unit. An organoleptic soil examination showed that the cold store surrounding includes a 40-cm thick humus (fertile soil) layer, with silty clay below it down to the depth of 2.0 m below ground level. The examination was made using a 2.50 m long geotechnical auger. The scope of experimental studies included continuous measurements of indoor air temperature, outdoor air temperature, floor temperature, ground temperature under the floor down to the depth of 1.50 m below the floor level, ground temperature in the cold store surrounding down to the depth of 1.50 m below ground level. Twenty-two PT100 sensors with 0.1 °C measurement accuracy were used. These parameters were measured for the chamber K2. The preliminary tests by the authors who analyzed the temperature and humidity conditions in both chambers did not show any statistically significant differences, so it was decided to continue with more accurate studies for the larger chamber (K2). The measurements were automatically recorded at 60-minintervals on a computer hard drive connected to the measuring instruments. The location of measurement points and the cold store layout are presented in Figure 1. The results of experimental studies allowed determining the thermal conditions inside the examined chamber and in the ground below it. The unit heat losses to the ground through the floor were estimated based on the relationship (1), where: Θi-indoor air temperature, (°C ), Θg-ground temperature interpolated linearly between points (A1-C1), (°C ), λs-thermal conductivity coefficient of soil, the used value was 1.9 W/mK, (EN ISO 6946), Rsi-heat penetration resistance, the used value was 0.14 m 2 K/W. The calculations were made based on the floor temperature measured 0.05 m below its surface and the actual curve of indoor air temperature over time was accounted for. Then, the unit values (W/m 2 ) were added and presented in kWh.
(1) (a) The next step involved the validation of the calculation model based on the elementary balances numerical method. The complex area of the cold store and the ground was divided into cuboidal balance-differential elements (Figure 2), using the model discretization. For each element, it is possible to calculate the increase of its internal energy, which depends on the temperature change (2), where: ∆Q-internal energy increase of the element, (J), C-total heat capacity of the element, (J·K −1 ), ∆θ-element temperature rise, (K).

∆ = • ∆ [J]
(2) The next step involved the validation of the calculation model based on the elementary balances numerical method. The complex area of the cold store and the ground was divided into cuboidal balance-differential elements (Figure 2), using the model discretization. For each element, it is possible to calculate the increase of its internal energy, which depends on the temperature change (2), where: ∆Q-internal energy increase of the element, (J), C-total heat capacity of the element, (J·K −1 ), ∆θ-element temperature rise, (K).
The heat flux flowing through the analyzed area is calculated in non-stationary terms, assuming the time step ∆τ (3), where: i-elements number, Φ i x i y i z , . . . -heat flux (heat flows) flowing between element i x i y i z and adjacent elements (W).
The computational analysis also takes into account the heat flux density, which is calculated on the basis of Fourier's law (4), where: → .
In the case of stationary heat flow through the balance-differential element, heat balance ∆Q = 0. For non-stationary flow, the temperature increase with the heat increase phenomenon should be taken Energies 2020, 13, 4970 5 of 15 into account. The formula for the heat balance can therefore be written as (5), where: ΣΦ-sum of heat flows between individual elements, (J), θ k+1 -temperature in the time step k + 1, (K), θ k -temperature in the time step k, (K), ∆τ-time step, [s].
In the case of stationary heat flow through the balance-differential element, heat balance ΔQ = 0. For non-stationary flow, the temperature increase with the heat increase phenomenon should be taken into account. The formula for the heat balance can therefore be written as (5), where: ΣΦ-sum of heat flows between individual elements, (J), θ k+1 -temperature in the time step k+1, (K), θ ktemperature in the time step k, (K), Δτ-time step, [s].
The next step of activities based on the balance-difference method is to write down the system of balance-difference equations. At this stage, the occurrence of a start condition that describes the temperature distribution in the previous time step should also be taken into account. In addition, the system of balance and difference equations should take into account: the current values of air temperature inside and outside the analyzed object, the temperature of the soil medium, as well as adiabatic planes. The above-mentioned parameters are called boundary conditions. As a result of the discretization of the area in the method of elementary balances, a discrete geometric model with a division grid is obtained. For each of these elements, the energy balance equation can be written using explicit difference quotient, implicit differential quotient and Crank-Nicolson differential quotient. Energy and temperature balance calculations were made for each such element, using a 60-min time step. The calculation model validation was made based on the temperature measured inside and The next step of activities based on the balance-difference method is to write down the system of balance-difference equations. At this stage, the occurrence of a start condition that describes the temperature distribution in the previous time step should also be taken into account. In addition, the system of balance and difference equations should take into account: the current values of air temperature inside and outside the analyzed object, the temperature of the soil medium, as well as adiabatic planes. The above-mentioned parameters are called boundary conditions. As a result of the discretization of the area in the method of elementary balances, a discrete geometric model with a division grid is obtained. For each of these elements, the energy balance equation can be written using explicit difference quotient, implicit differential quotient and Crank-Nicolson differential quotient. Energy and temperature balance calculations were made for each such element, using a 60-min time step. The calculation model validation was made based on the temperature measured inside and outside the cold store. In relation to the ground, the temperature was the third-type boundary condition. The initial temperature used was the average annual temperature for Krakow which is 8.8 • C. The calculations were made using the WUFIplus ® ver. 3.2.0 software (Fraunhofer Institut fur Bauphysik, Holzkirchen Branch, Germany).The validation results were subjected to a statistical correlation analysis, after checking the distribution normality with the Kolmogorov-Smirnov test, and then using the Spearman's Rank Test. After the completion of the validation and validation calculations, the variant analysis was started.  Table 1. In addition, in all variants the assumed minimum indoor air temperature was Θ i,min = 0 • C, maximum indoor air temperature Θ i,max = 4 • C, minimum indoor air relative humidity Rh i,min = 80%,and maximum indoor air relative humidity Rh i,max = 99%. Assumed was also natural air change 0.3 h −1 and infiltration 0.2 h −1 . The climatic data for a typical meteorological year for Krakow (TRY) available in WUFIplus ® were used in the calculation variant analysis. Based on the results obtained by Czapp and Charun [26] and by Jedynak [27], the following breathing heat values for the stored carrots were used: for 0 • C-18.9 W/t, for 4 • C-30.5 W/t. The intermediate values were linearly interpolated.
The obtained model validation results were analyzed for mean arithmetic difference in comparison to the measured data. The mean error was determined based on the Formula (2), where: ∆θ-mean error, θp-measured temperature, θo-calculated temperature, N-number of measurements.
Precise curves for indoor air temperature, floor temperature and ground temperature were obtained for each variant. The detailed analysis was performed of the cold store energy demand and of the heat exchange with the ground, with particular emphasis on the initial period of the cooling cycle and the whole yearly operation.

Experimental Studies
The average outdoor air temperature during the experimental studies conducted from 1 May 2017 to 30 April 2018 was 9.4 • C. The lowest temperature recorded in winter was −15.3 • C. The temperature Energies 2020, 13, 4970 7 of 15 inside the cold store was kept at 0.4-5.6 • C. During the technological break from the beginning of July to the end of September, the indoor microclimate followed the outdoor microclimate changes. No vegetables were stored in the chamber during that period and the refrigerating units were off, causing the air temperature in the chamber to increase to 15.3 • C in August. The maximum outdoor air temperature in that period was 22.3 • C. In the winter, the average temperature inside the cold store was from 1.0 • C to 1.6 • C, but in the remaining months it was near 4.0 • C (Figure 3a). The floor temperature analysis in point A1 showed that on the last day of the cooling cycle it was 4.2 • C. On the first day of the next cooling cycle, the floor temperature was 11.6 • C and gradually decreased. It was determined that the period necessary to lower the floor temperature to the value from the end of the cooling cycle in the analyzed variant is 18 days (Figure 3b-d). The annual ground temperature amplitude in point E1 was 29.4 • C, but at the depth of 1.5 m in the same measuring line (E4) it was 9.2 • C, which is very close the average annual outdoor air temperature(9.4 • C). The calculation results showed the highest ground temperature in point E1, which was 18.8 • C and was reached in June. The lowest ground temperature (−10.6 • C) was recorded in January in the same measurement point (Figure 3e,f).

Experimental Studies
The average outdoor air temperature during the experimental studies conducted from 1 May 2017 to 30 April 2018 was 9.4 °C. The lowest temperature recorded in winter was −15.3 °C. The temperature inside the cold store was kept at 0.4-5.6 °C. During the technological break from the beginning of July to the end of September, the indoor microclimate followed the outdoor microclimate changes. No vegetables were stored in the chamber during that period and the refrigerating units were off, causing the air temperature in the chamber to increase to15.3 °C in August. The maximum outdoor air temperature in that period was 22.3 °C. In the winter, the average temperature inside the cold store was from 1.0 °C to 1.6 °C, but in the remaining months it was near 4.0 °C (Figure 3a). The floor temperature analysis in point A1 showed that on the last day of the cooling cycle it was 4.2°C. On the first day of the next cooling cycle, the floor temperature was 11.6 °C and gradually decreased. It was determined that the period necessary to lower the floor temperature to the value from the end of the cooling cycle in the analyzed variant is 18 days (Figure  3b-d). The annual ground temperature amplitude in point E1 was 29.4 °C, but at the depth of 1.5 m in the same measuring line (E4) it was 9.2 °C, which is very close the average annual outdoor air temperature(9.4 °C). The calculation results showed the highest ground temperature in point E1, which was 18.8 °C and was reached in June. The lowest ground temperature (−10.6 °C) was recorded in January in the same measurement point (Figure 3e,f).   The analysis of heat exchange between the indoor and outdoor air and the ground air showed that the greatest heat gains of 27.5 kWh occurred in December. Due to technological break in summer, the heat flow to the ground in July was intensive, reaching the maximum of 37.2 kWh (Figure 4). In total, the ground delivered 1807.9 kWh of energy to the cold store, and the heat losses from the cold store to the ground were 1289.5 kWh in the whole period of study, mainly during the technological break. The beginning of technological break also marks the beginning of the heating process of the floor and the ground below it. In a cooling facility, this is undesirable because the accumulated energy surplus in the ground directly affects the increased cooling energy demand on the first days of the next cycle. After the cooling cycle ended, the floor in point A1 reached the temperature close to the indoor air temperature after 8 days and 14 h (Figure 5a). The temperature equalization in points B1 and C1 with the indoor air temperature was much quicker. Twenty-two hours had passed since the refrigeration unit was stopped until the temperature in these points reached the comparable value ( Figure 5b). Such large differences in the heating time allows a clear separation of the wall zone (up to 1.0 m from the outside wall) and the central zone. The wall zone is very susceptible to the outdoor weather, also to the action of the ground surrounding the cold store building. The results from measurement points located at other depths were also analyzed. The presence of wall zone was equally noticeable there. The analysis of heat exchange between the indoor and outdoor air and the ground air showed that the greatest heat gains of 27.5 kWh occurred in December. Due to technological break in summer, the heat flow to the ground in July was intensive, reaching the maximum of 37.2 kWh (Figure 4). In total, the ground delivered 1807.9 kWh of energy to the cold store, and the heat losses from the cold store to the ground were 1289.5 kWh in the whole period of study, mainly during the technological break.  The analysis of heat exchange between the indoor and outdoor air and the ground air showed that the greatest heat gains of 27.5 kWh occurred in December. Due to technological break in summer, the heat flow to the ground in July was intensive, reaching the maximum of 37.2 kWh (Figure 4). In total, the ground delivered 1807.9 kWh of energy to the cold store, and the heat losses from the cold store to the ground were 1289.5 kWh in the whole period of study, mainly during the technological break. The beginning of technological break also marks the beginning of the heating process of the floor and the ground below it. In a cooling facility, this is undesirable because the accumulated energy surplus in the ground directly affects the increased cooling energy demand on the first days of the next cycle. After the cooling cycle ended, the floor in point A1 reached the temperature close to the indoor air temperature after 8 days and 14 h (Figure 5a). The temperature equalization in points B1 and C1 with the indoor air temperature was much quicker. Twenty-two hours had passed since the refrigeration unit was stopped until the temperature in these points reached the comparable value (Figure 5b). Such large differences in the heating time allows a clear separation of the wall zone (up to 1.0 m from the outside wall) and the central zone. The wall zone is very susceptible to the outdoor weather, also to the action of the ground surrounding the cold store building. The results from measurement points located at other depths were also analyzed. The presence of wall zone was equally noticeable there. The beginning of technological break also marks the beginning of the heating process of the floor and the ground below it. In a cooling facility, this is undesirable because the accumulated energy surplus in the ground directly affects the increased cooling energy demand on the first days of the next cycle. After the cooling cycle ended, the floor in point A1 reached the temperature close to the indoor air temperature after 8 days and 14 h (Figure 5a). The temperature equalization in points B1 and C1 with the indoor air temperature was much quicker. Twenty-two hours had passed since the refrigeration unit was stopped until the temperature in these points reached the comparable value (Figure 5b). Such large differences in the heating time allows a clear separation of the wall zone (up to 1.0 m from the outside wall) and the central zone. The wall zone is very susceptible to the outdoor weather, also to the action of the ground surrounding the cold store building. The results from measurement points located at other depths were also analyzed. The presence of wall zone was equally noticeable there.

Calculation Model Validation
The calculations for the analyzed variants were performed using the cold store geometrical model built in the WUFI ® plus, which corresponded to actual dimensions. The first step included the calculation model validation (Figure 6).

Calculation Model Validation
The calculations for the analyzed variants were performed using the cold store geometrical model built in the WUFI ® plus, which corresponded to actual dimensions. The first step included the calculation model validation ( Figure 6).

Calculation Model Validation
The calculations for the analyzed variants were performed using the cold store geometrical model built in the WUFI ® plus, which corresponded to actual dimensions. The first step included the calculation model validation (Figure 6).  Correlation analysis was performed for results obtained from all measurement points based on the Spearman's Rank Test (the data distribution was not normal). The data distribution normality was tested previously using the Shapiro-Wilk Test. The highest correlation coefficient (0.95-0.99) was determined for measurement points E2, E3, E4, D2, D3, D4, C2, C3, C4, B2, B3, B4and A2, A3 and A4. Consequently, the conformity of calculation and measurement data at the depth of 0.50 m, 1.00 m and 1.50 m can be considered full. In case of points E1 and D1, located 0.05 m below the ground level, and A1, B1 and C1 located under the floor, the correlation coefficient was 0.85-0.88, so the conformity of measurement and calculation data should be considered very high. It should be emphasized that in winter the actual soil temperature at the depth of 0.10 m differed from the calculated values by 11.2 °C in line E and 11.0 °C in line D. The reason was the impossibility of accounting for the snow cover. The differences occurred only for a short time however, and were not significant in relation to the entire period of study. The statistical significance of differences was tested using the Kruskal-Wallis test (p < 0.05). The mean error for all measurement points was in the 1.0-3.0 °C range. The lowest error values in the 1.0-1.8 °C were obtained for the lines located inside the cold store(A, B and C), and the highest in the 1.4-3.0 °C for the outside lines (E and D). The greatest mean error was found in point E1 (3.0°C).

Calculation Variant Analysis
The calculation analysis was performed on 14 variants differing in terms of the technological break duration (Figure 7). Due to a large amount of data, the variant analysis results are given as average values and for the whole year, for better illustration of global differences between individual variants. Correlation analysis was performed for results obtained from all measurement points based on the Spearman's Rank Test (the data distribution was not normal). The data distribution normality was tested previously using the Shapiro-Wilk Test. The highest correlation coefficient (0.95-0.99) was determined for measurement points E2, E3, E4, D2, D3, D4, C2, C3, C4, B2, B3, B4and A2, A3 and A4. Consequently, the conformity of calculation and measurement data at the depth of 0.50 m, 1.00 m and 1.50 m can be considered full. In case of points E1 and D1, located 0.05 m below the ground level, and A1, B1 and C1 located under the floor, the correlation coefficient was 0.85-0.88, so the conformity of measurement and calculation data should be considered very high. It should be emphasized that in winter the actual soil temperature at the depth of 0.10 m differed from the calculated values by 11.2 • C in line E and 11.0 • C in line D. The reason was the impossibility of accounting for the snow cover. The differences occurred only for a short time however, and were not significant in relation to the entire period of study. The statistical significance of differences was tested using the Kruskal-Wallis test (p < 0.05). The mean error for all measurement points was in the 1.0-3.0 • C range. The lowest error values in the 1.0-1.8 • C were obtained for the lines located inside the cold store(A, B and C), and the highest in the 1.4-3.0 • C for the outside lines (E and D). The greatest mean error was found in point E1 (3.0 • C).

Calculation Variant Analysis
The calculation analysis was performed on 14 variants differing in terms of the technological break duration (Figure 7). Due to a large amount of data, the variant analysis results are given as average values and for the whole year, for better illustration of global differences between individual variants.
In all variants the action of low temperature inside the cold store was clearly noticeable down to the depth of 0. The simulation results have showed that the technological break duration has a significant impact on the cold store energy demand and the heat exchange with the ground, but also on the cooling time of the ground under the floor. Absence of technological break (variant 1) results in a decrease of energy gains from the ground by 20% relative to a three-month technological break (variant 14). The impact of technological break duration was clearly seen in terms of energy losses from the cold store to the ground (Figure 9). In case of a 91-day technological break (variant 14), the energy losses to the ground were 1289.5 kWh/a, while in case of absence of technological break (variant 1) this value was ninefold lower (147.5 kWh/a). On the cooling cycle end , the ground temperature in point A2 was 5.2 °C. During 91 days of technological break, the floor and the ground under it was gradually heated. The free flow of outdoor air into the cold store was assumed in order to reflect the actual conditions. The assumed end of technological break was 30 September. On this day, the soil temperature in point A2 was 11.0 °C. After the start of the cooling cycle, an intensive flow of the heat to the ground took place and the time necessary to cool down the ground in point A2 to the same temperature as on 30 June was 143 days. For comparison, in case of a seven-day technological break (1-7 July), the time necessary to cool down the ground to the  The simulation results have showed that the technological break duration has a significant impact on the cold store energy demand and the heat exchange with the ground, but also on the cooling time of the ground under the floor. Absence of technological break (variant 1) results in a decrease of energy gains from the ground by 20% relative to a three-month technological break (variant 14). The impact of technological break duration was clearly seen in terms of energy losses from the cold store to the ground ( Figure 9). In case of a 91-day technological break (variant 14), the energy losses to the ground were 1289.5 kWh/a, while in case of absence of technological break (variant 1) this value was ninefold lower (147.5 kWh/a). On the cooling cycle end , the ground temperature in point A2 was 5.2 °C. During 91 days of technological break, the floor and the ground under it was gradually heated. The free flow of outdoor air into the cold store was assumed in order to reflect the actual conditions. The assumed end of technological break was 30 September. On this day, the soil temperature in point A2 was 11.0 °C. After the start of the cooling cycle, an intensive flow of the heat to the ground took place and the time necessary to cool down the ground in point A2 to the same temperature as on 30 June was 143 days. For comparison, in case of a seven-day technological break (1-7 July), the time necessary to cool down the ground to the The simulation results have showed that the technological break duration has a significant impact on the cold store energy demand and the heat exchange with the ground, but also on the cooling time of the ground under the floor. Absence of technological break (variant 1) results in a decrease of energy gains from the ground by 20% relative to a three-month technological break (variant 14). The impact of technological break duration was clearly seen in terms of energy losses from the cold store to the ground ( Figure 9). In case of a 91-day technological break (variant 14), the energy losses to the ground were 1289.5 kWh/a, while in case of absence of technological break (variant 1) this value was ninefold lower (147.5 kWh/a). On the cooling cycle end , the ground temperature in point A2 was 5.2 • C. During 91 days of technological break, the floor and the ground under it was gradually heated. The free flow of outdoor air into the cold store was assumed in order to reflect the actual conditions. The assumed end of technological break was 30 September. On this day, the soil temperature in point A2 was 11.0 • C. After the start of the cooling cycle, an intensive flow of the heat to the ground took place and the time necessary to cool down the ground in point A2 to the same temperature as on 30 June was 143 days. For comparison, in case of a seven-day technological break (1-7 July), the time necessary to cool down the ground to the temperature as before the technological break was 1 day and 14 h (variant 2). An interesting relationship was obtained for a 70-day technological break (variant 11). In this case, the time necessary to cool down the ground was 141 days and was similar to results in variants 12-14. The reason is mainly the impact of the outdoor microclimate on the microclimate inside the cold store. In variant 11, the technological break ends in early September. In this time of the year, the outdoor air temperature falls, so even without the cooling system operation, the ground under the floor and the ground surrounding the building cools down spontaneously. relationship was obtained for a 70-day technological break (variant 11). In this case, the time necessary to cool down the ground was 141 days and was similar to results in variants 12-14. The reason is mainly the impact of the outdoor microclimate on the microclimate inside the cold store. In variant 11, the technological break ends in early September. In this time of the year, the outdoor air temperature falls, so even without the cooling system operation, the ground under the floor and the ground surrounding the building cools down spontaneously. The analysis of indoor air temperature simulation results showed that a six-week technological break (variant 6) allows maintaining the average indoor air temperature at the 3.95 °C level ( Table 2).  [7][8][9][10][11][12][13][14] contributes to the average annual indoor air temperature exceeding 4.0 °C, which significantly increases the cooling energy demand, particularly on the first days of the cooling cycle. The analysis of the floor temperature results (points A1-C1) showed that the technological break longer than four weeks causes the average floor temperature to The analysis of indoor air temperature simulation results showed that a six-week technological break (variant 6) allows maintaining the average indoor air temperature at the 3.95 • C level ( Table 2).

Discussion
Agricultural buildings used periodically allow numerous technological and operational modifications. They can contribute to rational management of energy necessary to maintain a given building type [28,29]. Current state of knowledge shows that the share of the ground in the total energy balance of residential and industrial buildings is about 5-10% [17]. As a result of continuous development of digital technologies, numerical methods are increasingly used for modelling and forecasting of physical phenomena [30][31][32]. However, computer simulations have some limitations and do not account for all factors which in reality can affect the studied parameters [23]. Validation of calculation models requires the results of field measurements encompassing the period of at least one year. Without such results, the simulation results could be burdened with a large error, particularly in case of untypical buildings as agricultural ones, and the obtained data should be considered estimates. There are studies that indicate the need for verification of the impact of the cold store thermal insulation in different floor location variants on the intensity of heat exchange with the ground and the share of the soil in the total building energy balance [30,33,34]. There has been no scientific paper to date which discusses the impact of the technological break in the cooling cycle on the heat exchange with the ground and the building energy management. This paper should be considered a case study and a foundation for further, more elaborated analyses. The studies in this regard on agricultural buildings are difficult because there are no relevant standards which could regulate some issues as it is the case in residential and industrial buildings. The agricultural buildings have different indoor condition, particularly air temperature and relative humidity. Depending on the building structure, area of its envelopes and most importantly its intended use, the percent share of envelopes in contact with the ground can significantly affect the energy gains and losses and the building energy management. In the future, the authors plan to use the model for further energy analyses related to the influence of other factors, such as. the type of stored products and their distribution on the shaping of selected parameters of the microclimate of the cold store. For this it is planned to use the CFD model.

Conclusions
The paper includes a multivariant analysis of the duration of technological break between the cooling cycles in a vegetable cold store located in southern Poland on the heat exchange with the ground. Due to absence of available standards which would allow a precise determination of the heat flux through the ground in such facilities, a numerical elementary balances method was used. The advantage of this approach was obtaining results for 14 variants in non stationary conditions, taking into account actual climatic conditions. The results of experimental studies showed that the conditions inside the cold store were stable and appropriate for stored crops. The technological break took place from the beginning of July to the end of September and the microclimate inside the cold store was shaped passively depending on the changes of the outdoor microclimate. The results of field measurements of the indoor and outdoor air temperatures and the ground temperature were used to validate the calculation model in the WUFI ® plus software, and then to verify the conformity of the calculated results with the real results. After numerous validation calculations, a satisfactory level of conformity of calculated and measured results was obtained. The validated model was used in a simulation of the indoor microclimate, ground temperature and the heat exchange with the ground for all 14 variants. The calculations were made in a non stationary approach. The calculation model validation was performed based on the results of field studies of the indoor and outdoor air temperature and the ground temperature in five measurement lines at the depth of 0.05, 0.50, 1.00 and 1.50 m. The obtained validation results showed a very high correlation between the measured and calculated data and absence of significant differences (p < 0.05). The annual floor temperature amplitude for variant 1 (cooling all year long) in point A1 was 3.9 • C. In case of three-month technological break (variant 14) the temperatures of the floor and the ground under it were in a much wider range. The annual temperature amplitude in point A1 was 13.9 • C, in point A2-8.3 • C, A3-7.2 • C and A4-6.3 • C. The technological break duration has therefore a significant impact on the thermal conditions of the floor and ground under it. The analysis of simulation results showed also that absence of technological break (variant1) results in a decrease of energy gains from the ground by 20% relative to a three-month technological break. The analysis of the floor temperature calculations (points A1-C1) showed that the technological break longer than four weeks causes the average floor temperature to exceed 4.0 • C. Therefore, the optimum solution is technological break lasting 7-35 days. At a shorter technological break, the facility can be used for storing young vegetables in the summer period.