This study analyzes the dynamics of temperature and moisture in a cylindrical silo with a conical roof and floor used for storing corn in the Bajío region of Mexico, considering conditions both with and without aeration. The model incorporates external temperature fluctuations, solar radiation, grain moisture equilibrium with air humidity through the sorption isotherm (water activity), and grain respiration to simulate real storage conditions. The model is based on continuity, momentum, energy, and moisture conservation equations in porous media. This model was solved using the finite element method (FEM) to evaluate temperature and interstitial humidity variations during January and May, representing cold and warm environmental conditions, respectively. The simulations show that, without aeration, grain temperature progressively accumulates in the center and bottom region of the silo, reaching critical values for safe storage. In January, the low ambient temperature favors the natural dissipation of heat. In contrast, in May, the combination of high ambient temperatures and solar radiation intensifies thermal accumulation, increasing the risk of grain deterioration. However, implementing aeration periods allowed for a reduction in the silo’s internal temperature, achieving more homogeneous cooling and reducing the threats of mold and insect proliferation. For January, an airflow rate of 0.15 m³/(min·ton) was optimal for maintaining the temperature within the safe storage range (≤17 °C). In contrast, in May, neither this airflow rate nor the accumulation of 120 h of aeration was sufficient to achieve optimal storage temperatures. This indicates that, under warm conditions, the aeration strategy needs to be reconsidered, assessing whether a higher airflow rate, longer periods, or a combination of both could improve heat dissipation. The results also show that interstitial relative humidity remains stable with nocturnal aeration, minimizing moisture absorption in January and preventing excessive drying in May. However, it was identified that aeration period management must be adaptive, taking environmental conditions into account to avoid issues such as re-wetting or excessive grain drying.

This study presents a numerical comparison of dimensional and non-dimensional implementations of the classical benchmark problem of steady natural convection in a two-dimensional differentially heated square cavity over the Rayleigh-number range 10³ ≤ Ra ≤ 10¹². The novelty of this work lies in the systematic comparison of both formulations under identical numerical conditions, providing an implementation-oriented assessment of their benchmark accuracy, mesh sensitivity, continuation strategy, and computational efficiency. The governing equations of mass, momentum, and energy conservation were solved under the Boussinesq approximation using primitive variables and the finite-element method. Since both formulations are theoretically equivalent descriptions of the same physical problem, the purpose of this work is not to reassess their physical validity, but to examine their numerical behavior under identical benchmark conditions in terms of mesh sensitivity, continuation strategy, benchmark accuracy, and computational efficiency. In the dimensional implementation, temperature differences of 1, 10, 25, and 50 K were considered, and the corresponding cavity lengths were determined from the Rayleigh-number definition. The main calculations were performed with ΔT = 25 K, while ΔT = 50 K was retained only as an exploratory sensitivity case. Boundary-layer refinement was applied along the vertical walls over the range 10³ ≤ Ra ≤ 10¹². The results show that, once the near-wall gradients are properly resolved, both implementations predict essentially identical average Nusselt numbers, with a maximum relative difference of 0.022%. Temperature contours, stream-function distributions, and centerline profiles also exhibited the same structural behavior in both cases. For the reference mesh, the dimensional implementation exhibited a lower computational cost than the non-dimensional implementation, resulting in an approximately 31% reduction in computation time. These results demonstrate that solving the governing equations directly in dimensional variables provides a numerically efficient and physically interpretable alternative while preserving the benchmark accuracy of the classical non-dimensional formulation.