This paper presents pore unit assembly-discrete element model (PUA-DEM), a pore-scale hydromechanical framework that resolves interactions between mobile granular particles and multiphase fluids in unsaturated granular media. The framework uniquely integrates DEM with pore-scale hydrodynamic models to capture unsaturated flow dynamics, while leveraging a two-way coupling mechanism to ensure bidirectional fluid-grain feedback through stabilized domain partitioning. Further innovations include a dynamic pore-merging and retriangulation algorithm that enhances computational efficiency for large-scale systems. Validated against experimental data for glass beads and Ottawa sand, PUA-DEM accurately reproduces critical hydromechanical phenomena-including capillary/viscous fingering, wetting-induced granular swelling/collapse, and quasi-static deformation-under diverse saturation and loading regimes. Numerical case studies reveal how capillary forces and wetting fluid saturation collectively govern granular response, from pore-scale meniscus evolution to macroscale flow instabilities. By bridging pore-and particle-scale physics, PUA-DEM advances predictive modeling of partially saturated granular systems, offering transformative insights for geohazard mitigation, sustainable agriculture, pharmaceutical manufacturing, and energy-related engineering applications.

The dynamic behaviour of granular media can be observed widely in nature and in many industrial processes. Yet, the modelling of such media remains challenging as they may act with solid-like and fluid-like properties depending on the rate of the flow and can display a varying apparent friction, cohesion and compressibility. Over the last two decades, the mu(I)-rheology has become well established for modelling granular liquids in a fluid mechanics framework where the apparent friction mu depends on the inertial number I. In the geo-mechanics community, modelling the deformation of granular solids typically relies on concepts from critical state soil mechanics. Along the lines of recent attempts to combine critical state and the mu(I)-rheology, we develop a continuum model based on modified cam-clay in an elastoplastic framework which recovers the mu(I)-rheology under flow. This model permits a treatment of plastic compressibility in systems with or without cohesion, where the cohesion is assumed to be the result of persistent inter-granular attractive forces. Implemented in a two- and three-dimensional material point method, it allows for the trivial treatment of the free surface. The proposed model approximately reproduces analytical solutions of steady-state cohesionless flow and is further compared with previous cohesive and cohesionless experiments. In particular, satisfactory agreements with several experiments of granular collapse are demonstrated, albeit with shear bands which can affect the smoothness of the surface. Finally, the model is able to qualitatively reproduce the multiple steady-state solutions of granular flow recently observed in experiments of flow over obstacles.