Influence of inter-grain cementation stiffness on the effective elastic properties of porous Bentheim sandstone

Effective elastic properties of porous media are known to be significantly influenced by porosity.In this paper,we investigated the influence of another critical factor,the inter-grain cementation stiffness,on the effective elastic properties of a granular porous rock(Bentheim sandstone)using an advanced numerical workflow with realistic rock microstructure and a theoretical model.First,the disparity between the experimentally tested elastic properties of Bentheim sandstone and the effective elastic properties predicted by empirical equations was analysed.Then,a micro-computed tomography(CT)-scan based approach was implemented with digital imaging software AVIZO to construct the 3D(three-dimensional)realistic microstructure of Bentheim sandstone.The microstructural model was imported to a mechanics solver based on the 3D finite element model with inter-grain boundaries modelled by cohesive elements.Loading simulations were run to test the effective elastic properties for different shear and normal intergrain cementation stiffness.Finally,a relation between the macroscale Young’s modulus and inter-grain cementation stiffness was derived with a theoretical model which can also account for porosity explicitly.Both the numerical and theoretical results indicate the influence of the inter-grain cementation stiffness,on the effective elastic properties is significant for porous sandstone.The calibrated normal and shear stiffnesses at the inter-grain boundaries are 1.2×10^(5) and 4×10^(4) GPa/m,respectively.

Modelling of massive particulates for breakwater engineering using coupled FEMDEM and CFD

The seaward slope of many breakwaters consists of thousands of interlocking units of rock or concrete comprising a massive granular system of large elements each weighing tens of tonnes. The dumped quarry materials in the core are protected by progressively coarser particulates. The outer armour layer of freely placed units is intended to both dissipate wave energy and remain structurally stable as strong flows are drawn in and out of the particulate core. Design guidance on the mass and shape of these units is based on empirical equations derived from scaled physical model tests. The main failure mode for armour layers exposed to severe storms is hydraulic instability where the armour units of concrete or rock are subjected to uplift and drag forces which can in turn lead to rocking, displacement and collisions sufficient to cause breakage of units. Recently invented armour unit designs making up such granular layered system owe much of their success to the desirable emergent properties of interlock and porosity and how these combine with individual unit structural strength and inertial mass. Fundamental understanding of the forces governing such wave-structure interaction remains poor. We use discrete element and combined finite-discrete element methods to model the granular solid skeleton of randomly packed units coupled to a CFD code which resolves the wave dynamics through an interface tracking technique. The CFD code exploits several methods including a compressive advection scheme, node movement, and general mesh optimization. We provide the engineering context and report progress towards the numerical modelling of instability in these massive granular systems.