A Comparison between Separately Calibrated P-{\alpha} and Mesoscale Models for Weak Shock Compaction of Granular Sugar

TL;DR

This paper compares two modeling approaches—continuum P-alpha and mesoscale simulations—for predicting particle velocity in granular sugar under weak shock, highlighting their different calibration mechanisms and sensitivities.

Contribution

It demonstrates that continuum and mesoscale models require different calibration strategies, with continuum models sensitive to parameters like crush-out pressure and mesoscale models focusing on porosity.

Findings

Both models match experimental particle velocities.

Continuum model sensitivity to calibration parameters.

Mesoscale model primarily depends on porosity.

Abstract

This study compares calibration strategies for predicting particle velocity in granular sugar subjected to weak shock loading, using measurements from flyer-plate impact experiments as a benchmark. Two computational approaches are evaluated: a continuum-based P-alpha Menikoff model requiring calibration of effective constitutive parameters, and mesoscale simulations that explicitly resolve grain geometry and porosity. Both models can match the measured particle-velocity histories, but only through fundamentally different calibration mechanisms. In the P-alpha model, a pressure-dependent yield strength is essential and the response remains highly sensitive to parameter choices such as the crush-out pressure. In contrast, mesoscale simulations are far less sensitive to parameter tuning and instead depend primarily on the physical state variable of porosity, represented in 2D through an equivalent mapping of the 3D specimen. These results show that continuum parameters act as effective surrogates for underlying grain-scale processes, whereas mesoscale modeling identifies porosity as the dominant control on macroscopic wave onset, highlighting distinct calibration pathways and interpretive implications for each modeling approach.