Discrete Element Method Simulation of Wheat Bulk Density as Affected by Grain Drop Height and Kernel Size Distribution

Highlights The predicted bulk density of two wheat varieties varied with drop height, similar to the experiments. The percentage composition of three kernel size fractions in the wheat varieties affected the bulk density. Accurate particle shape representation simulated the heap profile better but required longer computational time. The single-sphere model is more practical to use because of its higher accuracy and lower computational cost. Abstract . Grain bulk density varies widely depending on kernel properties and handling practices. The discrete element method (DEM) can model such behavior at the particle level, including wide-ranging interactions with equipment. The objective of this study was to develop a DEM model to predict wheat bulk density as affected by grain drop height and kernel size distribution. The bulk density of two wheat varieties was measured experimentally for a range of drop heights with a modified test weight per bushel apparatus and was simulated in EDEM v2018.1 using single-sphere and five-sphere particle models that accounted for three kernel size fractions. For both particle models, simulations matched the observed behavior, showing a bulk density increase with increasing drop height and bulk density differences between varieties due to different kernel size fractions. The single-sphere particle model predicted the bulk density with higher accuracy than the five-sphere particle model, while the five-sphere model, which more accurately represented the shape of wheat kernels, allowed better simulation of the heap profile at the cost of longer computation time. These particle models can be used to simulate bulk density of wheat under compaction and to improve prediction models of grain pack factor for wheat. Keywords: Bulk density, DEM, Drop height, Size distribution, Wheat.