Numerical prediction of the effect of vent arrangements on the ventilation and energy transfer in a multi-span glasshouse using a bi-band radiation model

The distributed climate inside a greenhouse mainly results from mass, momentum and energy transfers occurring through the openings or through the cover. These transfers can be assessed using computational fluid dynamics (CFD) techniques. However, to date very few studies considered properly the radiative exchanges between the surfaces inside the calculation domain by solving simultaneously the radiative transfer equation (RTE) and the convective ones. The object of this work is to implement this equation in order to get an overview of the resulting climate inside the greenhouse. The analysis focuses on the heterogeneity of the climatic parameters in the canopy vicinity as well as on the ventilation rate. For the purpose of the study, a two-dimensional steady-state CFD model was developed. It is based on the resolution of the Navier–Stokes equations with the Boussinesq assumption and a k–e closure. Solar and atmospheric radiations were included by solving the RTE and distinguishing short [0–3 μm] and long wavelength [3–100 μm] contributions. The virtual tracer gas technique was implemented to assess the ventilation rate. The model was first successfully validated against experimental data considering a glasshouse divided into two compartments separated by a plastic partition and comparing simulated and recorded temperatures inside the shelter and along the walls (roof and sides). Simulations were then carried out for a 2500 m 2 four-span sloped-roof glasshouse covered with a 4-mm-thick horticulture glass and equipped with continuous roof vents and one side-opening. Results show that roof-opening configurations combined with side-vent location strongly affect inside ventilation and microclimate parameters. Thus, for hot summer conditions with relatively low wind, computed ventilation rates varied from 2 to 40 air changes per hour whereas temperature differences varied from 3 to 13 °C. This study also showed that other characteristics such as climate heterogeneity must be investigated in order to define the best ventilation configuration and that strong velocity and temperature heterogeneity can occur at the plant level. Simulations, however, suggest that the symmetric roof-vent configuration seems to be a good compromise between cooling and homogeneity performances.

[1]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[2]  Wilfried Brutsaert,et al.  Evaporation into the atmosphere : theory, history, and applications , 1982 .

[3]  Nicola Pearsall,et al.  Heat transfer modelling in thermophotovoltaic cavities using glass media , 2005 .

[4]  Diego L. Valera,et al.  Measurement and simulation of climate inside Almerı́a-type greenhouses using computational fluid dynamics , 2004 .

[5]  P. Richards,et al.  Appropriate boundary conditions for computational wind engineering models using the k-ε turbulence model , 1993 .

[6]  Dafydd Gibbon,et al.  1 User’s guide , 1998 .

[7]  S. Patankar Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[8]  Pietro Picuno,et al.  Computational analysis of ventilation in greenhouses at zero- and low-wind speeds. , 1997 .

[9]  Thierry Boulard,et al.  Effect of Vent Arrangement on Windward Ventilation of a Tunnel Greenhouse , 2004 .

[10]  Pedro J. Coelho Fundamentals of a new method for the solution of the radiative transfer equation ? ? A preliminary v , 2005 .

[11]  Juan Ignacio Montero,et al.  COMPUTATIONAL FLUID DYNAMIC MODELLING OF NIGHT-TIME ENERGY FLUXES IN UNHEATED GREENHOUSES , 2005 .

[12]  Pietro Picuno,et al.  Analysis of the efficiency of greenhouse ventilation using computational fluid dynamics , 1997 .

[13]  Pierre-Emmanuel Bournet,et al.  MEASUREMENT AND NUMERICAL SIMULATION OF MICROCLIMATE IN A NATURALLY-VENTILATED LARGE MULTI-SPAN GLASSHOUSE , 2005 .

[14]  Thierry Boulard,et al.  Analysis of Greenhouse Ventilation Efficiency based on Computational Fluid Dynamics , 2006 .

[15]  J. B Campen,et al.  Determination of greenhouse-specific aspects of ventilation using three-dimensional Computational Fluid Dynamics , 2003 .

[16]  T. H. Short,et al.  TWO-DIMENSIONAL NUMERICAL SIMULATION OF NATURAL VENTILATION IN A MULTI-SPAN GREENHOUSE , 2000 .

[17]  Thierry Boulard,et al.  Effect of Roof and Side Opening Combinations on the Ventilation of a Greenhouse Using Computer Simulation , 2007 .

[18]  A. Maldonado,et al.  Physical properties of ZnO:F obtained from a fresh and aged solution of zinc acetate and zinc acetylacetonate , 2006 .

[19]  T. H. Short,et al.  Verification of Computational Fluid Dynamic Temperature Simulations in a Full-Scale Naturally Ventilated Greenhouse , 2001 .