Effects of Particle Segregation and Grain Pressure on Ventilation Airflow and Temperature–Humidity Distribution in Maize Pilot Silo

The distribution of grain particles within a silo influences heat and moisture transfer during stored grain ventilation, leading to grain quality losses. A study on porosity distribution analysis and ventilation tests was conducted in a pilot silo with a height of 3 m, a diameter of 1.5 m, and a conical dome height of 0.85 m. The E-B constitutive model was incorporated into the secondary development of FLAC3D 5.0 to analyze the vertical pressure distribution in the grain bulk. An anisotropic porosity distribution model for the maize bulk was developed, accounting for both vertical pressure and segregation mechanisms. The differences in airflow and heat transfer during ventilation between isotropic and anisotropic porosity distributions were quantified. A nonlinear model was innovatively proposed to predict the temperature front curve (TFC) during ventilation as affected by porosity variation. The results indicate that friction between the maize kernel and the silo wall led to vertical pressure at the center of the bottom that was 10.7% higher than that near the wall. The average surface porosity of the maize bulk was 2.8% higher than at the bottom. This led to a minimum porosity of 0.409 at the center of the silo bottom, due to the combined effect of impact during the loading process and vertical pressure. The numerical simulation demonstrated excellent consistency with the experimental data. At a supply vent air velocity of 0.126 m/s, an increase in the maize bulk height from 0.725 m to 2.9 m resulted in reductions in airflow rate and average relative humidity of 20.3% and 9.67%. The airflow velocity near the wall was 13.4% higher than that in the center, leading to a faster cooling rate in the peripheral region compared to the center of the maize bulk. The airflow velocity based on the isotropic porosity model was higher at the center than that predicted by the anisotropic model, whereas the opposite trend was observed near the wall. The temperature front during ventilation based on the anisotropic porosity model exhibited a concave curve. A nonlinear model was developed to predict this temperature front, showing strong agreement with computational data.