Many widely used numerical models of rock fracture based on mesoscale laboratory test characterisation of effective ‘intact’ strength parameters neglect microstructure effects. They therefore cannot explain grain boundary and pore effects on crack propagation and consequently are inadequate for models of rock destruction that exploit point and indentation stresses. Understanding deep drilling processes involving drill-bit buttons and/or water-jetting where rock loading is concentrated in domains with fewer mineral grains will therefore require models with microstructure. To investigate microscale failure mechanisms of granular rocks in diverse scenarios, we target a porous sandstone and introduce a novel workflow consisting of a computerized tomography (CT) based microstructure construction approach and a complementary mechanical numerical approach. The construction approach extracts the realistic rock microstructure and transforms the large voxel number CT-scan data into significantly fewer triangular elements. The finite-discrete element method (FDEM) with grain-based model (GBM) is adopted to solve the mechanics. The microscale failure mechanism of sandstone during the Brazilian test was thoroughly analysed using the numerical results together with the post failure CT-scan test data. The build-up of compressive and tensile stress chains, micro-crack nucleation, local relaxation, chain switching and final crack-path development exploiting pores was illustrated, revealing the micro-to-macro failure mechanism in time and space. Fracture paths in the specimens during Brazilian tensile test were dominated by the pores and the inter-grain boundaries. The tensile strength of the inter-grain joints was estimated to be at least 3.67 times the mesoscale specimen's intact tensile strength, while the pores account for 72.76% of the fracture path. The influence of the cementation distribution and microscale discontinuities was investigated with numerical cases.@en

Radial water jet drilling (RJD) is a method of enhancing heat recovery by accessing and connecting to high permeable zones within geothermal reservoirs. The wall rock geometry behind an advancing water jet borehole under in-situ conditions is largely unknown. Water jet drilling tests were performed on 300 mm cubical blocks of weak porous sandstone under true-triaxial boundary stress conditions at the Delft Technical University (DTU) rock mechanics laboratory. Some of these tests showed distinct breakout features depending on the applied stress field. Geometries of resulting boreholes are recovered using X-Ray CT scans, and are analysed using segmentation software (Avizo). The code Solidity, using a combined finite-discrete element method with a cohesive zone fracture model, simulates stress take-up and wall shearing giving breakouts comparable to the experiments. The results lead to the suggestion that criteria based on Kirsch solutions would be suitable to provide general guidance on in-situ stress and rock strength conditions free of breakouts. FEMDEM models appear well-suited to examine geometries and dimensions that can be sustained by given strengths under deeper in-situ conditions. Wall-rock failure and a process of jet-hole enlargement together with the potential benefits of greater heat recovery arising from larger holes is also briefly discussed.

@en

Numerical modelling has the potential to probe the complexity of the interacting physics of rubble mound armour systems. Through forward modelling of armour unit packs, stochastic variables such as unit displacement and maximum contact force per unit during an external oscillatory disturbance can be predicted. The combined finite-discrete element method (FEMDEM) is a multi-body method ideally suited to model the behaviour of the armour layer system and the stresses generated within complex shape units. In this paper we highlight the latest developments made with the application of FEMDEM technology to breakwater modelling including realistic rock underlayer and concrete unit layer topologies, maximum contact force distributions, internal unit stresses, fracture and unit breakages. Finally, fully coupled wave and multi-body armour unit motion with internal dynamic stress generation is illustrated.

@en

The influence of in-situ stresses on flow processes in fractured rock is investigated using a novel modelling approach. The combined finite-discrete element method (FEMDEM) is used to model the deformation of a fractured rock mass. The fracture wall displacements and aperture changes are modelled in response to uniaxial and biaxial stress states. The resultant changes in flow properties of the rock mass are investigated using the Complex Systems Modelling Platform (CSMP++). CSMP++ is used to model single-phase flow through fractures with variable aperture and a permeable rock matrix. The study is based on a geological outcrop mapping of a low density fracture pattern that includes the realism of intersections, bends and segmented features. By applying far-field (boundary) stresses to a square region, geologically important phenomena are modelled including fracture-dependent stress heterogeneity, the re-activation of pre-existing fractures (i.e. opening, closing and shearing), the propagation of new fractures and the development of fault zones. Flow anisotropy is investigated under various applied stresses and matrix permeabilities. In-situ stress conditions that encourage a closing of fractures together with a more pervasive matrix-dominated flow are identified. These are compared with conditions supporting more localised flow where fractures are prone to dilatational shearing and can be more easily exploited by fluids. The natural fracture geometries modelled in this work are not perfectly straight, promoting fracture segments that dilate as they shear. We have demonstrated the introduction of several realistic processes that have an influence on natural systems: fractures can propagate with wing cracks; there is the potential for new fractures to connect with existing fractures, thus increasing the connectivity and flow; blocks can rotate when bounded by fractures, bent fractures lead to locally different aperture development; highly heterogeneous stress distributions emerge naturally. Results presented in this work provide a mechanically rigorous demonstration that a change in the stress state can cause reactivation of pre-existing fractures and channelling of flow in critically stressed fractures.

@en

FEMDEM modelling which combines the multi-body particle interaction and motion modelling (i.e. Discrete Element Model, DEM) with the ability to model internal deformation of arbitrary shape (Finite Element Model, FEM) has been applied to breakwater models. There are two versions of a FEMDEM solver developed; Y3D-D is for deformable materials and is required for dynamic and static stress analysis and Y3D-R is the rigid version often used to numerically construct the armour unit packs. This paper also reports the placement protocols: POSITIT. FEMDEM modelling deals with solids interactions and is one modelling component that is to be coupled to other modelling technologies e.g. CFD, interface tracking, wave models, porous media etc. so that the key fluid-solid interactions can be modelled in a full scale virtual breakwater alongside work on scaled hydraulic laboratory models and prototype structures. The latest developments of two-way coupled interactions of waves with coastal structures are also described in this paper.

@en