Three-Dimensional Modeling of Natural Convection During Postharvest Storage of Corn and Wheat in Metal Silos in the Bajío Region of Mexico

This study presents a three-dimensional numerical analysis of natural convection during the postharvest storage of corn and wheat in a galvanized steel silo with a conical roof and floor, measuring 3 m in radius and 18.7 m in height, located in the Bajío region of Mexico. Simulations were carried out specifically for December, a period characterized by cold ambient temperatures (10–20 °C) and comparatively lower solar radiation than in warmer months, yet still sufficient to induce significant heating of the silo’s metallic surfaces. The governing conservation equations of mass, momentum, energy, and species were solved using the finite volume method under the Boussinesq approximation. The model included grain–air sorption equilibrium via sorption isotherms, as well as metabolic heat generation: for wheat, a constant respiration rate was assumed due to limited biochemical data, whereas for corn, respiration heat was modeled as a function of grain temperature and moisture, thereby more realistically representing metabolic activity. The results, obtained for December storage conditions, reveal distinct thermal and hygroscopic responses between the two grains. Corn, with higher thermal diffusivity, developed a central thermal core reaching 32 °C, whereas wheat, with lower diffusivity, retained heat in the upper region, peaking at 29 °C. Radial temperature profiles showed progressive transitions: the silo core exhibited a delayed response relative to ambient temperature fluctuations, reflecting the insulating effect of grain. In contrast, grain at 1 m from the wall displayed intermediate amplitudes. In contrast, zones adjacent to the wall reached 40–41 °C during solar exposure. In comparison, shaded regions exhibited minimum temperatures close to 15 °C, confirming that wall heating is governed primarily by solar radiation and metal conductivity. Axial gradients further emphasized critical zones, as roof-adjacent grain heated rapidly to 38–40 °C during midday before cooling sharply at night. Relative humidity levels exceeded 70% along roof and wall surfaces, leading to condensation risks, while core moisture remained stable (~14.0% for corn and ~13.9% for wheat). Despite the cold ambient temperatures typical of December, neither temperature nor relative humidity remained within recommended safe storage ranges (10–15 °C; 65–75%). These findings demonstrate that external climatic conditions and solar radiation, even at reduced levels in December, dominate the thermal and hygroscopic behavior of the silo, independent of grain type. The identification of unstable zones near the roof and walls underscores the need for passive conservation strategies, such as grain redistribution and selective ventilation, to mitigate fungal proliferation and storage losses under non-aerated conditions.