Thermal and CFD Analyses of Sustainable Heat Storage-Based Passive Greenhouse Dryer Operating in No-Load Condition

<jats:p>This article presents a comprehensive study on thermal and computational fluid dynamics (CFD) analysis of an innovative greenhouse dryer designed for passive operation under a no-load condition. The dryer incorporates hybrid thermal storage at the floor and a reflective mirror with thermocoal as the north wall, transforming a classical even-span greenhouse dryer into an efficient and effective system. The experimentation was conducted under clear sky conditions, with variations in global solar radiation (GSR) ranging from 166.6 to 1209 W/m2, resulting in an average value of 875.9 W/m2. The variations in GSR influenced other ambient parameters, including ambient temperature (28.7 °C to 35.6 °C), ambient relative humidity (33.2% to 45.7%), and ambient wind speed (0.1 to 1.02 m/s). Indoor parameters of the proposed dryer, such as inside temperature (31 °C to 47.35 °C), inside relative humidity (31.1% to 39.1%), ground temperature (44.2 °C to 70.6 °C), and outlet temperature (29 °C to 45.35 °C), were measured hourly. The average values of these parameters were 41.25 °C, 35.31%, 61.65 °C, and 39.25 °C, respectively. Quantitative parameters, including heat loss, overall heat transfer coefficient, coefficient of diffusion, and instantaneous efficiency, were calculated to evaluate the dryer’s performance. The proposed dryer exhibited an improved range of overall heat transfer coefficients (3.87 to 5.03 W/m2 K) compared to the modified greenhouse dryer under passive mode and the conventional greenhouse under passive mode. CFD analysis provided temperature distribution plots showing a progressively increasing range of temperatures near the trays, ranging from 310 K to 335 K, suitable for natural convection drying. The findings highlight the superior performance of the innovative dryer compared to contemporary systems. This research contributes to the advancement of drying technology and holds potential for applications in the agriculture and food processing industries.</jats:p>

section_content: Drying is a process to reduce moisture from a product, which is one of the ancient techniques used for food or agricultural products for safe storage [1, 2]. Food products, particularly fruits and vegetables require hot climatic conditions in the temperature range of 45 to 65 • C for safe drying to keep their edible and nutritional stuffs unchanged [3, 4]. Food preservation is the need to reduce food losses and cost-effective transportation, which results from economic final products [2, 3]. It has been widely accepted as an important tool for a long time [5]. Therefore, well-designed energy-efficient drying technologies are required [6, 7]. Developing countries are facing an energy crisis; hence, solar drying is the most attractive drying technology since sun-radiant energy is available in abundant amounts [8]. However, solar radiation is limited to use during sunshine hours due to its intermittent nature and non-availability at night [9]. Hence, heat storage is required to store or absorb solar energy during sunshine hours and that can be utilized after sunshine hours [10]. The hybrid single tank thermocline designed by incorporating a heat concrete block in between D mannitol and adipic acid is an efficient thermal storage system [11]. It can be three types (a) sensible heat storage, (b) latent heat storage, and (c) thermochemical (combination of sensible and latent heat storage) storage [12]. In sensible heat storage, the temperature increases of any liquid or solid without any phase changes but in latent heat storage phase change occurs in heat-storing material [13]. While in a thermochemical heat storage system the heat energy is stored in the form of chemical bonds, and it releases this stored energy by reversible chemical reaction [14]. Solar drying systems can be classified as direct drying and indirect drying systems. A direct drying system is where crops are directly exposed to sun radiation inside the dryer, whereas in an indirect system, an air heating system is required and crops are kept in another cabinet (Figure 1a, b ) [15]. Rock beds were used as sensible thermal storage. Simultaneously, experiments were performed without a rock bed to compare its result. The utilization of heat storage material improved the drying efficiency of the dryer with rock bed by 2-3%. A natural convection greenhouse dryer was developed using different types of heat storage materials like concrete, rock-bed, and sand to evaluate their performance [18]. The drying system reduced the moisture content of coconuts from 52% to 7% wet basis using concrete heat storage material for a total of 78 h and successively reduced drying time by 55% as compared to open sun drying. For the same heat storage material, sand took 66 h and reduced 62% drying time; however, for heat storage material, the rock-bed acquired only 53 h and reduced drying time by 69%. The drying efficiency of the dryer with thermal Generally, greenhouse dryers and tunnel dryers come under the category of direct drying systems. Various researchers have applied heat storage material in greenhouse drying systems. A greenhouse dryer was developed and integrated with PCM (CaCl 2 •6H 2 O) on the north wall [16]. The rise of inside air temperature was found to be in the range of 6 to 12 • C and the cover temperature was 4 to 5 • C due to the use of 4 cm thick PCM on the north wall as a thermal storage. A natural convection solar tunnel drier was developed using a rock bed to examine its performance for copra drying [17]. 

Rock beds were used as sensible thermal storage. Simultaneously, experiments were performed without a rock bed to compare its result. The utilization of heat storage material improved the drying efficiency of the dryer with rock bed by 2-3%. A natural convection greenhouse dryer was developed using different types of heat storage materials like concrete, rock-bed, and sand to evaluate their performance [18]. The drying system reduced the moisture content of coconuts from 52% to 7% wet basis using concrete heat storage material for a total of 78 h and successively reduced drying time by 55% as compared to open sun drying. For the same heat storage material, sand took 66 h and reduced 62% drying time; however, for heat storage material, the rock-bed acquired only 53 h and reduced drying time by 69%. The drying efficiency of the dryer with thermal storage material was found to always be higher at 11. 6%, 11%, and 9. 5% using rock-bed, sand, and concrete, respectively. A north-wall insulated greenhouse dryer was developed with a heating collector for gooseberry and bitter gourd drying for complete dehydration of moisture [15]. The study yielded significant results by effectively reducing the drying time. In this research, a greenhouse dryer was developed and integrated with various types of heat storage beds, including gravel beds, black-painted gravel beds, ground beds, and concrete floor beds [19]. One of the notable findings was that the black-painted gravel bed achieved a maximum room air temperature of 64. 4 • C, which corresponds to a heat gain percentage of 53%. This demonstrates the effectiveness of the integration in terms of heat transfer and performance improvement. 

Another aspect of the study focused on investigating the impact of different corrugation interruptions on thermohydrodynamic characteristics and heat transfer performance using 3D corrugated tubes. The researchers explored the flow field and heat transfer enhancement by combining corrugated tubes with twisted tape within a 3D circular tube, considering variations in dimple configurations. The objective was to develop correlations for improved thermo-hydraulic flow and heat transfer performance. Furthermore, the researchers conducted an analysis to examine the flow structure and heat transfer improvement in a 3D circular tube utilizing various axial groove turbulator configurations. They carried out a numerical study to investigate the effect of turbulators on thermal flow and heat performance in a 3D pipe [20]. 

Despite the existing studies on greenhouse dryers and heat storage materials, there is a research gap in understanding the specific effects of different heat storage materials on the drying efficiency of greenhouse dryers. While previous research has explored the use of concrete, rock-bed, and sand as heat storage materials, there is a need for further investigation to compare their performance and determine the optimal heat storage material for improved drying efficiency. Additionally, there is a lack of studies focusing on the integration of reflective mirrors, insulation, and thermal storage to minimize heat loss and maximize solar radiation utilization in greenhouse dryers. This research gap highlights the need for a comprehensive study that combines these elements to enhance the energy efficiency and overall performance of greenhouse dryers. 

The present experimental work focuses on the design and modernization of a greenhouse dryer with the objective of reducing energy consumption under passive mode. The novelty of this study lies in the application of a reflective mirror on the north side and the use of a 10 mm thick polystyrene sheet as insulation to maximize solar radiation utilization and minimize heat loss from the north wall. Additionally, thermal storage material is incorporated to mitigate conductive heat loss from the ground. These modifications aim to enhance the effectiveness and performance of the dryer. 

The determination of the overall heat transfer coefficient, calculation of heat loss, evaluation of the instantaneous efficiency factor, and the use of computational fluid dynamics (CFD) analysis are key objectives of this research. These parameters provide essential insights into the performance and effectiveness of the dryer, which are crucial for optimizing its operation and energy efficiency. 

The applications of this research extend to the field of sustainable agriculture, specifically in the area of crop drying. Efficient and effective greenhouse dryers can contribute to preserving the quality and nutritional value of agricultural produce [21, 22]. The findings of this study can benefit the agricultural community, farmers, researchers, and policymakers working towards sustainable food production and post-harvest management [23] [24] [25]. The knowledge gained from this research can help optimize drying processes, reduce energy consumption, and improve the overall quality of dried crops, leading to economic and environmental benefits.
section_name: Introduction

section_name: Material and Method

section_content: The proposed system is the modified version of the classical even-span greenhouse dryer. There are two major modifications incorporated in the classical system, namely in the north wall and ground of the dryer. In the north wall of the greenhouse dryer where maximum heat loss happens in the classical dryer, the transparent wall is replaced by a reflecting mirror with thermocoal. A detailed description of the system is presented in Table 1. A schematic diagram and a real view of the proposed system are presented in Figure 2a, b. Reflecting mirrors are used to restrict outgoing radiation and heat. Thermocoal is kept outside of the glass should no heat loss happen through the glass. By this arrangement, a substantial amount of heat can be preserved. This enhances the drying efficiency of the proposed system [26] [27] [28]. The ground floor of the dryer is modified. Heat storage material is applied on the floor of the dryer so that system thermal performance can be enhanced. Dual-layer thermal storage material is used where the bottom layer is the phase change material (latent heat storage material) and the upper layer is gravel (sensible heat storage material). 

Firstly, paraffin wax is mixed with black painted gravel and put above the black PVC sheet. This makes a 50 mm thicker layer of this storage material at the bottom of the greenhouse dryer. Thermo-physical properties of paraffin wax are presented in Table 2.
section_name: Experimental Setup

section_content: The proposed system is examined critically on an hourly basis with the help of technically effective instruments. Details of the instruments are presented in Table 3. Both indoor and outdoor/ambient parameters are observed to clearly understand the thermal behavior of the proposed system. The following ambient parameters are observed: global solar radiation, ambient temperature, ambient relative humidity, and wind speed. However, the following parameters are observed at the outlet point: outlet temperature, outlet relative humidity, and air speed. The following parameters are observed inside the dryer namely ground temperature, inside temperature, and inside relative humidity.
section_name: Instrumentation

section_content: The proposed system is observed on an hourly basis from 9 a. m. to 5 p. m. to identify its thermal behavior. An experiment was conducted on 15 February 2022 at the University campus of Birla Institute of Technology Mesra, Ranchi-India.
section_name: Experimentation

section_name: Performance and CFD Analysis

section_content: The following mathematical parameters need to be calculated to properly evaluate the thermal performance of the proposed system. Such parameters are the overall heat transfer coefficient, dimensionless numbers, coefficient of diffusion, heat loss, and instantaneous thermal loss efficiency factor. Overall heat transfer coefficient factor further depends upon other factors namely radiative heat transfer coefficient, and convective heat transfer coefficients [29, 30]. These parameters are evaluated based on the following Equation (1). 

Connellan [8] proposed a different relation for 'U' as in Equation ( 2 ):

It can be analyzed by using Equations ( 3 )-( 6 ): T grd -T rm

(5)

Here, h evp can be neglected as it is very small. The dimensionless number can be evaluated by using Equations ( 7 )- (10). 

The parameter W can be evaluated by using Equation ( 11 ) as the floor of the dryer is in rectangular shape:

Instantaneous heat loss can be evaluated by Equation ( 12 ):

Coefficient of diffusion can be evaluated by Equation ( 13 ):

Heat loss can be calculated by Equation ( 14 ):
section_name: Data Reduction

section_content: The governing equations considered in the present CFD analysis are based on the following Navier-Stokes equations [20, 31, 32] :

Momentum Equation

Energy equation

where v is the fluid flow velocity, ρ is the density, F is the applied external force, S m is the added mass source, ρg is the gravitational force, S h represents the sensible enthalpy source, E represents the total energy contained by the fluid, and p is the pressure.
section_name: Mathematical Expressions for CFD Modelling

section_content: The transport equations for turbulence kinetic energy (k) and rate of dissipation (ε) are given as:

G k and G b denote the development of turbulence kinetic energy due to the mean velocity gradients and buoyancy, respectively. The contribution of the fluctuating dilatation is represented by Y M. The turbulence Prandtl numbers for k and ε are represented, respectively, by σ k and σ ε. S k and S ε denote the source terms for k and ε, respectively, and the constants are represented by C 1ε, C 2ε, and C 3ε. 

The following combination of k and ε yields the eddy viscosity,

where C µ = 0. 09, is the model constant. 

The term G k is modeled identically for all the k and ε models (Standard, RNG, and Realizable). 

G k =µ t S 2 (21)

where S stands for the modulus of the mean rate-of-strain tensor and is given as,

The turbulence developed by buoyancy is given by

where Pr t stands for the turbulent-Prandtl number, g i stands for the ith direction gravitational vector component. Pr t is 0. 85 for the standard and realizable k-ε models, whereas it is 1/α for the RNG k-ε model. 

where α 0 = 1. 0. The contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate

The turbulent Mach Number M t is defined as

where α is the speed of the sound (a = √ γRT).
section_name: The Standard k-ε Model

section_content: where → r is the position vector, → s stands for the direction vector, → s represents the vector along scattering direction, s denotes the path length, a is the absorption coefficient, n is the refractive index, σ s stands for the scattering coefficient, σ is the Stefan-Boltzmann constant (5. 669 × 10 -8 W/m 2 K 4 ), I represents the radiation intensity, which relies on → r and → s, T stands for the local temperature, ϕ denotes the phase function, and Ω denotes the solid angle.
section_name: The radiative transfer equation (RTE) for an absorbing, emitting, and scattering medium at position

section_content: The boundary conditions considered have been evaluated from the average of the experimental measurements of the present study. In Ansys 19. 2, a solar ray tracing model was employed to evaluate solar heat flux. Figure 3a,b shows the inside and outside heat flux variation. Table 4 shows the boundary conditions on different surfaces of the greenhouse dryer. The turbulent Mach Number Mt is defined as

where α is the speed of the sound ( = �)

The radiative transfer equation (RTE) for an absorbing, emitting, and scattering medium at position ⃗ in the direction ⃗ is,

where ⃗ is the position vector, ⃗ stands for the direction vector, ′ ��⃗ represents the vector along scattering direction, s denotes the path length, a is the absorption coefficient, n is the refractive index, σs stands for the scattering coefficient, σ is the Stefan-Boltzmann constant (5. 669 × 10 -8 W/m 2 -K 4 ), I represents the radiation intensity, which relies on ⃗ and ⃗, T stands for the local temperature, Փ denotes the phase function, and Ω denotes the solid angle.
section_name: Boundary Conditions

section_content: The boundary conditions considered have been evaluated from the average of the experimental measurements of the present study. In Ansys 19. 2, a solar ray tracing model was employed to evaluate solar heat flux. Figure 3a,b shows the inside and outside heat flux variation. Table 4 shows the boundary conditions on different surfaces of the greenhouse dryer.
section_name: Boundary Conditions

section_content: The geometric model for the studied configurations has been developed using the 3D geometry creator "Design Modular" of the Ansys Workbench (V19. 2). The working fluid is considered air, which has a density and viscosity of 1. 225 kg/m 3 and 1. 885 × 10 -5 kg/ms, respectively. At all solid surfaces, the no-slip boundary condition is prescribed. The simulations are carried out using the Standard K-ε turbulence model, which is also utilized to model the entire computational domain and estimate the values of performance parameters. The SIMPLEC (semi-implicit method for pressure-linked equation consistent) algorithm is used to assure the pressure and velocity coupling. The discretization of pressure was achieved with a standard scheme and for momentum discretization, the QUICK scheme was considered. The second-order upwind scheme discretizes the energy equations and the dissipation rate. For all the conservation equations, the solver's tolerance is set to 10 -5. Although the convergence of residuals is reached at about 1000 iterations, the simulations are run for additional 300 iterations to ensure higher accuracy.
section_name: Numerical Method

section_content: The computational domain for each of the considered models has been meshed using the tetrahedral meshing elements. Three levels of meshes have been considered to ensure grid independence of the solutions, viz: 262,856, 522,887, and 758,982 meshing elements. However, the solutions in terms of outlet air velocity remain unchanged when the mesh size is increased from 522,887 to 758,982 (Table 5 ). This suggests the attainment of grid independence of the results. Thus, for the grid generation of the numerical models, the element size corresponding to 522,887 is considered for the present analysis. The generated grid for the entire computational domain and the cut section model has been
section_name: Grid Independency

section_content: The computational domain for each of the considered models has been meshed using the tetrahedral meshing elements. Three levels of meshes have been considered to ensure grid independence of the solutions, viz: 262,856, 522,887, and 758,982 meshing elements. However, the solutions in terms of outlet air velocity remain unchanged when the mesh size is increased from 522,887 to 758,982 (Table 5 ). This suggests the attainment of grid independence of the results. Thus, for the grid generation of the numerical models, the element size corresponding to 522,887 is considered for the present analysis. The generated grid for the entire computational domain and the cut section model has been shown in Figure 4. Table 6 has exhibited the variation in maximum air outlet velocity with mesh size. then compared with the CFD results. The maximum average deviation is below 5%, which confirms the reliability of the considered computational model for the performance estimation of the greenhouse dryer. Table 7 shows the comparison of the experimental model and the numerical model. The temperature at the outlet of the dryer is measured with a digital thermometer for different hours on 21 June 2021. This experimental result for the outlet air temperature is then compared with the CFD results. The maximum average deviation is below 5%, which confirms the reliability of the considered computational model for the performance estimation of the greenhouse dryer.
section_name: Grid Independency

section_name: Result and Discussion

section_content: The experiment was conducted in a clear sky condition without any abrupt weather change. Global solar radiation (GSR) varies from 166. 6 to 1209 W/m 2 with an average value of 875. 9 W/m 2 (as shown in Figure 5a ). GSR is the core parameter from which all other parameters depend, which may be ambient temperature, ambient relative humidity, or ambient wind speed. Since the variation of GSR is in a regular pattern so variation of other dependent parameters is also in a regular pattern. Variations of ambient temperature, ambient relative humidity, and ambient wind speed are 28. 7-35. 6 •, 33. 2-45. 7%, 0. 1-1. 02 m/s, and average values of these parameters are 31. 98 • C, 40. 08%, 0. 32 m/s, respectively (as shown in Figure 5b ).
section_name: Variation of Ambient Parameters during Experimentation

section_content: The experiment was conducted in a clear sky condition without any abrupt weather change. Global solar radiation (GSR) varies from 166. 6 to 1209 W/m 2 with an average value of 875. 9 W/m 2 (as shown in Figure 5a ). GSR is the core parameter from which all other parameters depend, which may be ambient temperature, ambient relative humidity, or ambient wind speed. Since the variation of GSR is in a regular pattern so variation of other dependent parameters is also in a regular pattern. Variations of ambient temperature, ambient relative humidity, and ambient wind speed are 28. 7-35. 6°, 33. 2-45. 7%, 0. 1-1. 02 m/s, and average values of these parameters are 31. 98 °C, 40. 08%, 0. 32 m/s, respectively (as shown in Figure 5b ).
section_name: Variation of Ambient Parameters during Experimentation

section_content: Inside temperature, inside relative humidity, ground temperature, and outlet temperature are measured on an hourly basis. The variation of these parameters is 31-47. 35 • C, 31. 1-39. 1%, 44. 2-70. 6 • C, and 29-45. 35 • C, respectively. The average value of these parameters is 41. 25 • C, 35. 31%, 61. 65 • C, and 39. 25 • C, respectively. Variations of these parameters are regular because these parameters are dependent upon GSR. Due to the storage concept and reflective mirror, inside parameters are quite high which is highly favorable for drying and other thermal applications.
section_name: Variation of Indoor Parameters of Proposed Dryer during Experimentation

section_content: The proposed system properly examined their thermal performance with the help of standard parameters like convective heat transfer coefficient from the ground (h c,gr ), convective heat transfer coefficient from the canopy (h c,ca ), radiative heat transfer coefficient from the ground (h r,gr ), overall heat transfer coefficient (U), heat loss, coefficient of diffusion, dimensionless number and instantaneous thermal loss efficiency factor. Overall heat transfer coefficient depends upon convective and radiative heat transfer coefficients. Hourly variations of these parameters are shown in Figure 6. The variation of h c,gr, h c,ca, h r,gr, and U are 2. 47-3. 55 W/m 2 K, 7. 58-11. 08 W/m 2 K, 6. 05-7. 39 W/m 2 K, and 3. 87-5. 03 W/m 2 K, respectively, also the average value of these parameters are 3. 14 W/m 2 K, 8. 43 W/m 2 K, 6. 91 W/m 2 K, and 4. 28 W/m 2 K, respectively.
section_name: Variation of Mathematical Parameters

section_content: Inside temperature, inside relative humidity, ground temperature, and outlet temperature are measured on an hourly basis. The variation of these parameters is 31-47. 35 °C, 31. 1-39. 1%, 44. 2-70. 6 °C, and 29-45. 35 °C, respectively. The average value of these parameters is 41. 25 °C, 35. 31%, 61. 65 °C, and 39. 25 °C, respectively. Variations of these parameters are regular because these parameters are dependent upon GSR. Due to the storage concept and reflective mirror, inside parameters are quite high which is highly favorable for drying and other thermal applications.
section_name: Variation of Indoor Parameters of Proposed Dryer during Experimentation

section_content: The proposed system properly examined their thermal performance with the help of standard parameters like convective heat transfer coefficient from the ground (hc,gr), convective heat transfer coefficient from the canopy (hc,ca), radiative heat transfer coefficient from the ground (hr,gr), overall heat transfer coefficient (U), heat loss, coefficient of diffusion, dimensionless number and instantaneous thermal loss efficiency factor. Overall heat transfer coefficient depends upon convective and radiative heat transfer coefficients. Hourly variations of these parameters are shown in Figure 6. The variation of hc,gr, hc,ca, hr,gr, and U are 2. 47-3. 55 W/m 2 K, 7. 58-11. 08 W/m 2 K, 6. 05-7. 39 W/m 2 K, and 3. 87-5. 03 W/m 2 K, respectively, also the average value of these parameters are 3. 14 W/m 2 K, 8. 43 W/m 2 K, 6. 91 W/m 2 K, and 4. 28 W/m 2 K, respectively. As per Prakash and Kumar, the characteristic graph for an efficient greenhouse dryer should intercept at zero [33]. A characteristic graph is drawn between instantaneous thermal loss efficiency factor verses (Tr -Ta/Ir). A characteristic graph for the proposed system is presented in Figure 7 and intercepts at zero, which validated the modification in the greenhouse dryer. As per Prakash and Kumar, the characteristic graph for an efficient greenhouse dryer should intercept at zero [33]. A characteristic graph is drawn between instantaneous thermal loss efficiency factor verses (T r -T a /I r ). A characteristic graph for the proposed system is presented in Figure 7 and intercepts at zero, which validated the modification in the greenhouse dryer. The heat loss and coefficient of diffusivity are calculated for the proposed system based on observed experimental data. Hourly variations of these parameters are shown in Figure 8. Heat loss varies from 105. 58 to 549. 42 W and the average value is 428. 33 W. However, the variation of the coefficient of diffusivity is 0. 1531 × 10 -3 to 0. 1547 × 10 -3 and the average value is 0. 1541 × 10 -3. The heat loss and coefficient of diffusivity are calculated for the proposed system based on observed experimental data. Hourly variations of these parameters are shown in Figure 8. Heat loss varies from 105. 58 to 549. 42 W and the average value is 428. 33 W. However, the variation of the coefficient of diffusivity is 0. 1531 × 10 -3 to 0. 1547 × 10 -3 and the average value is 0. 1541 × 10 -3. The heat loss and coefficient of diffusivity are calculated for the proposed system based on observed experimental data. Hourly variations of these parameters are shown in Figure 8. Heat loss varies from 105. 58 to 549. 42 W and the average value is 428. 33 W. However, the variation of the coefficient of diffusivity is 0. 1531 × 10 -3 to 0. 1547 × 10 -3 and the average value is 0. 1541 × 10 -3. The temperature distribution for the GHD for the selected day (21 June 2021) is shown in Figure 10a,b. For the simulations in Ansys Fluent, the direct and diffuse irradiation data obtained from the solar ray tracing is coupled with the discrete ordinate model and the energy model. From the inside view of the greenhouse dryer, it is seen that the maximum temperature is at the base of the GHD, i. e., of the phase change material. The temperature near the trays varies in the range of 335 K to 310 K, which is appropriate for product drying with natural convection. The maximum velocity at the outlet is 0. 837 m/s whereas the inlet velocity is 0. 32 m/s. obtained from the solar ray tracing is coupled with the discrete ordinate model and the energy model. From the inside view of the greenhouse dryer, it is seen that the maximum temperature is at the base of the GHD, i. e., of the phase change material. The temperature near the trays varies in the range of 335 K to 310 K, which is appropriate for product drying with natural convection. The maximum velocity at the outlet is 0. 837 m/s whereas the inlet velocity is 0. 32 m/s.
section_name: Variation of Mathematical Parameters

section_content: The chemical and physical features of the products to be dried are directly influenced by relative humidity. The relative humidity is high where the temperature is low and viceversa. The spatial distribution of relative humidity is shown in Figure 11. The maximum relative humidity in GHD is 59. 8% whereas the minimum relative humidity is 9. 9%. It is observed that the higher relative humidity regions overlap with the regions of wet air entering from the outside. Thus, the focus should be made on managing the ventilation
section_name: Relative Humidity inside the Solar GHD

section_content: The chemical and physical features of the products to be dried are directly influenced by relative humidity. The relative humidity is high where the temperature is low and viceversa. The spatial distribution of relative humidity is shown in Figure 11. The maximum relative humidity in GHD is 59. 8% whereas the minimum relative humidity is 9. 9%. It is observed that the higher relative humidity regions overlap with the regions of wet air entering from the outside. Thus, the focus should be made on managing the ventilation areas during higher humidity periods. Good ventilation management ensures higher efficiency of the greenhouse dryer.
section_name: Relative Humidity inside the Solar GHD

section_content: The comparative analysis presented in Table 8 validates the superiority of the proposed system compared to the existing systems. The modified greenhouse dryer under passive mode, implemented with the innovative modifications discussed in this study, demonstrates superior performance compared to the greenhouse under passive mode. The proposed system showcases an improved range of overall heat transfer coefficients (U) compared to both the modified greenhouse dryer under passive mode and the conventional greenhouse under passive mode. This enhanced range (3. 87-5. 03 W/m 2 K) in the proposed system reflects its superior heat transfer capabilities and improved efficiency. Additionally, the characteristics graph of the proposed system intersects at zero, similar to the modified greenhouse dryer under the passive mode, indicating the validity and effectiveness of the modifications implemented in the proposed system. This comparative analysis further strengthens the conclusions drawn in this study, highlighting the superior performance of the proposed system compared to existing systems. The innovative modifications in the greenhouse dryer design contribute to its enhanced heat transfer efficiency and overall effectiveness in solar drying processes. These findings emphasize the significance and novelty of the proposed system in advancing the field of solar drying technologies.
section_name: Validation

section_content: The study presented herein focused on the design, fabrication, experimentation, and analysis of an innovative greenhouse dryer operating under natural convection mode in no-load conditions, with a particular emphasis on clear sky conditions. The primary objective was to enhance the classical even-span greenhouse dryer by incorporating hybrid thermal storage and a reflective mirror with thermocoal, thereby improving its
section_name: Conclusions

section_content: The study presented herein focused on the design, fabrication, experimentation, and analysis of an innovative greenhouse dryer operating under natural convection mode in no-load conditions, with a particular emphasis on clear sky conditions. The primary objective was to enhance the classical even-span greenhouse dryer by incorporating hybrid thermal storage and a reflective mirror with thermocoal, thereby improving its drying efficiency. Through a comprehensive analysis of the experimental data and calculations of various heat transfer parameters, the following conclusions were drawn: i. 

The variation of key parameters such as global solar radiation (GSR), ambient temperature, ambient relative humidity, and ambient wind speed followed regular patterns, with average values of 875. 9 W/m 2, 31. 98 335 K to 310 K, providing an appropriate environment for natural convection-based product drying. The maximum relative humidity recorded in the greenhouse dryer was 59. 8%, while the minimum relative humidity was 9. 9%. Proper management of ventilation areas is crucial during periods of higher humidity to ensure optimal drying performance. vii. Comparison with existing systems validated the superiority of the proposed system, exhibiting an improved range of overall heat transfer coefficients compared to the modified greenhouse dryer under passive mode and the conventional greenhouse under passive mode. 

These findings highlight the enhanced heat transfer efficiency and overall effectiveness of the proposed innovative greenhouse dryer. The modifications implemented in the dryer design contribute to the advancement of solar drying technologies. The quantitative data and analysis presented in this study offer valuable insights for researchers, scientists, and entrepreneurs in the field of solar drying, paving the way for the development of sustainable and efficient drying systems.
section_name: Conclusions

section_content: The present work is focused on no-load conditions under natural convection mode and clear sky conditions. While these conditions are important to evaluate the performance of the innovative greenhouse dryer, they may not fully capture its behavior under different operational scenarios and real-world conditions. Future studies should consider expanding the scope to include varying loads, weather conditions, and geographical locations to more comprehensively assess the dryer's performance. Also, the present study does not reflect the needs and requirements of large-scale applications. Investigating the system's performance with different crop types and varying drying capacities would provide a more comprehensive understanding of its versatility and applicability.
section_name: Limitation of the Present Work

section_content: The study on the innovative greenhouse dryer presented in this work paves the way for future research and development in the field of solar drying. Several potential areas for further exploration arise from this study. First, future investigations can focus on evaluating the performance of the greenhouse dryer under extended operational conditions, including varying loads, diverse weather patterns, and different geographic locations. This will provide a more comprehensive understanding of its versatility and applicability. Second, researchers can explore the use of advanced materials, such as high-capacity phase change materials (PCMs), and innovative designs to further enhance the energy storage capabilities and overall efficiency of the system. Third, the integration of renewable energy sources, such as solar photovoltaic systems or wind turbines, presents an opportunity to augment the energy supply of the greenhouse dryer and reduce dependence on conventional energy sources. Fourth, the implementation of intelligent control systems and automation technologies can optimize the drying process by continuously monitoring and adjusting key parameters, resulting in improved efficiency, reduced energy consumption, and enhanced product quality. Finally, conducting comprehensive economic analyses and feasibility studies will be crucial to assess the cost-effectiveness and scalability of the proposed system, considering factors such as initial investment, operational costs, and potential savings. By addressing these areas of future research, solar drying technologies can be further advanced, leading to more sustainable, efficient, and economically viable drying solutions for various agricultural and industrial applications.
section_name: Future Scope

section_name: Nomenclature

section_name: Data Availability Statement:

section_content: No data were used to support this study.
section_name: Data Availability Statement:

section_content: The authors would like to thank King Saud University, Riyadh, Saudi Arabia, with researchers supporting project number RSP2023R145.
section_name: Acknowledgments:

section_name: Institutional Review Board Statement: Not applicable.

section_content: Author Contributions: Conceptualization, A. A., O. P., S. K. S., P. S. C., R. C. and S. S. methodology, A. A., O. P., S. K. S., P. S. C., R. C. and S. S. ; formal analysis, A. A., O. P., S. K. S., P. S. C., R. C. and S. S. ; investigation, A. A., O. P., S. K. S., P. S. C., R. C. and S. S. ; writing-original draft preparation, A. A., O. P., S. K. S., P. S. C., R. C. and S. S. ; writing-review and editing, S. S., R. K., S. M. T., A. K., B. S. and S. S. U. ; supervision, S. S., R. K., S. M. T., A. K., B. S. and S. S. U. ; project administration, S. S., B. S. and S. S. U. ; funding acquisition, S. S., B. S. and S. S. U. All authors have read and agreed to the published version of the manuscript. 

Funding: This study received funding from King Saud University, Saudi Arabia through researchers supporting project number (RSP2023R145). Additionally, the APCs were funded by King Saud University, Saudi Arabia through researchers supporting project number (RSP2023R145).

section_content: Informed Consent Statement: Not applicable.
section_name: Institutional Review Board Statement: Not applicable.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: