Fluid flow through porous media using distinct element based numerical method

Many analytical and numerical methods have been developed to describe and analyse fluid flow through the reservoir’s porous media. The medium considered by most of these models is continuum based homogeneous media. But if the formation is not homogenous or if there is some discontinuity in the formation, most of these models become very complex and their solutions lose their accuracy, especially when the shape or reservoir geometry and boundary conditions are complex. In this paper, distinct element method (DEM) is used to simulate fluid flow in porous media. The DEM method is independent of the initial and boundary conditions, as well as reservoir geometry and discontinuity. The DEM based model proposed in this study is appeared to be unique in nature with capability to be used for any reservoir with higher degrees of complexity associated with the shape and geometry of its porous media, conditions of fluid flow, as well as initial and boundary conditions. This model has first been developed by Itasca Consulting Company and is further improved in this paper. Since the release of the model by Itasca, it has not been validated for fluid flow application in porous media, especially in case of petroleum reservoir. In this paper, two scenarios of linear and radial fluid flow in a finite reservoir are considered. Analytical models for these two cases are developed to set a benchmark for the comparison of simulation data. It is demonstrated that the simulation results are in good agreement with analytical results. Another major improvement in the model is using the servo controlled walls instead of particles to introduce tectonic stresses on the formation to simulate more realistic situations. The proposed model is then used to analyse fluid flow and pressure behaviour for hydraulically induced fractured and naturally fractured reservoir to justify the potential application of the model

Discrete element models of soil-geogrid interaction

Geogrids are the geosynthetics of choice for soil reinforcement applications. To evaluate the efficiency of geogrid reinforcement, several methods are used including field tests, laboratory tests and numerical modeling. Field studies consume long period of time and conducting these investigations may become highly expensive because of the need for real-size structures. Laboratory studies present also significant difficulties: large-size testing machines are required to accommodate realistic geogrid designs. The discrete element method (DEM) may be used as a complementary tool to extend physical testing databases at lower cost. Discrete element models do not require complex constitutive formulations and may be fed with particle scale data (size, strength, shape) thus reducing the number offree calibration parameters. Discrete element models also are well suited to problems in which large displacements are present, such as geogrid pullout. This paper reviews the different approaches followed to model soil-geogrid interaction in DEM and presents preliminary results from pull-out conditions.Peer ReviewedPostprint (author's final draft

Numerical modeling of the tension stiffening in reinforced concrete members via discontinuum models

[prova tipográfica]This study presents a numerical investigation on the fracture mechanism of tension stiffening phenomenon in reinforced concrete members. A novel approach using the discrete element method (DEM) is proposed, where three-dimensional randomly generated distinct polyhedral blocks are used, representing concrete and one-dimensional truss elements are utilized, representing steel reinforcements. Thus, an explicit representation of reinforced concrete members is achieved, and the mechanical behavior of the system is solved by integrating the equations of motion for each block using the central difference algorithm. The inter-block interactions are taken into consideration at each contact point with springs and cohesive frictional elements. Once the applied modeling strategy is validated, based on previously published experimental findings, a sensitivity analysis is performed for bond stiffness, cohesion strength, and the number of truss elements. Hence, valuable inferences are made regarding discontinuum analysis of reinforced concrete members, including concrete-steel interaction and their macro behavior. The results demonstrate that the proposed phenomenological modeling strategy successfully captures the concrete-steel interaction and provides an accurate estimation of the macro behavior

Friction phenomena and their impact on the shear behaviour of granular material

In the discrete element simulation of granular materials, the modelling of contacts is crucial for the prediction of the macroscopic material behaviour. From the tribological point of view, friction at contacts needs to be modelled carefully, as it depends on several factors, e.g. contact normal load or temperature to name only two. In discrete element method (DEM) simulations the usage of Coulomb’s law of friction is state of the art in modelling particle–particle contacts. Usually in Coulomb’s law, for all contacts only one constant coefficient of friction is used, which needs to reflect all tribological effects. Thus, whenever one of the influence factors of friction varies over a wide range, it can be expected that the usage of only one constant coefficient of friction in Coulomb’s law is an oversimplification of reality. For certain materials, e.g. steel, it is known that a dependency of the coefficient of friction on the contact normal load exists. A more tribological tangential contact law is implemented in DEM, where the interparticle friction coefficient depends on the averaged normal stress in the contact. Simulations of direct shear tests are conducted, using steel spheres of different size distributions. The strong influence of interparticle friction on the bulk friction is shown via a variation of the constant interparticle friction coefficient. Simulations with constant and stress-dependent interparticle friction are compared. For the stress-dependent interparticle friction, a normal stress dependency of the bulk friction is seen. In the literature, measurements of different granular materials and small normal loads also show a stress dependency of the bulk friction coefficient. With increasing applied normal stress, the bulk friction coefficient reduces both in the experiments and in the simulations

Discrete element modelling of scaled railway ballast under triaxial conditions

The aim of this study is to demonstrate the use of tetrahedral clumps to model scaled railway ballast using the discrete element method (DEM). In experimental triaxial tests, the peak friction angles for scaled ballast are less sensitive to the confining pressure when compared to full-sized ballast. This is presumed to be due to the size effect on particle strength, whereby smaller particles are statistically stronger and exhibit less abrasion. To investigate this in DEM, the ballast is modelled using clumps with breakable asperities to produce the correct volumetric deformation. The effects of the quantity and properties of these asperities are investigated, and it is shown that the strength affects the macroscopic shear strength at both high and low confining pressures, while the effects of the number of asperities diminishes with increasing confining pressure due to asperity breakage. It is also shown that changing the number of asperities only affects the peak friction angle but not the ultimate friction angle by comparing the angles of repose of samples with different numbers of asperities

Evaluation of Mode I Fracture Toughness Assisted by the Numerical Determination of K-Resistance

The fracture toughness of a rock often varies depending on the specimen shape and the loading type used to measure it. To investigate the mode I fracture toughness using semi-circular bend (SCB) specimens, we experimentally studied the fracture toughness using SCB and chevron bend (CB) specimens, the latter being one of the specimens used extensively as an International Society for Rock Mechanics (ISRM) suggested method, for comparison. The mode I fracture toughness measured using SCB specimens is lower than both the level I and level II fracture toughness values measured using CB specimens. A numerical study based on discontinuum mechanics was conducted using a two-dimensional distinct element method (DEM) for evaluating crack propagation in the SCB specimen during loading. The numerical results indicate subcritical crack growth as well as sudden crack propagation when the load reaches the maximum. A K-resistance curve is drawn using the crack extension and the load at the point of evaluation. The fracture toughness evaluated by the K-resistance curve is in agreement with the level II fracture toughness measured using CB specimens. Therefore, the SCB specimen yields an improved value for fracture toughness when the increase of K-resistance with stable crack propagation is considered

Numerical Modeling of a Granular Collapse Immersed in a Viscous Fluid

The three-dimensional unsteady collapse of the granular column in a viscous fluid has been investigated with an IBM/DEM approach. Present numerical simulations allow one to confirm quantitatively several experimental observations of Rondon et al. regarding morphology, characteristic sizes of granular deposits and the basal pressure below the column. In the presented simulations, the collapse dynamics is controlled by the viscous time Tv. To our knowledge, a numerical approach, e.g. the IBM/DEM method, is able for the first time to capture the pore pressure feedback phenomenon in flowing fluid-grains mixture. The effect of the initial packing fraction has a great influence of the dynamics of granular collapse in the simulation results as in the experiments of Rondon et al. Furthermore, the IBM/DEM permits to investigate the inner state of the granular column during the collapse, in particular, the evolution of the pressure field inside the granular column can be analyzed which is difficult to do in experiments. Simulations of the collapse of a granular column immersed in a fluid can be performed in the inertial and free-fall regimes as well varying the nature of the fluid and/or the particles

Macro deformation and micro structure of 3D granular assemblies subjected to rotation of principal stress axes

This paper presents a numerical investigation on the behavior of three dimensional granular materials during continuous rotation of principal stress axes using the discrete element method. A dense specimen has been prepared as a representative element using the deposition method and subjected to stress rotation at different deviatoric stress levels. Significant plastic deformation has been observed despite that the principal stresses are kept constant. This contradicts the classical plasticity theory, but is in agreement with previous laboratory observations on sand and glass beads. Typical deformation characteristics, including volume contraction, deformation non-coaxiality, have been successfully reproduced. After a larger number of rotational cycles, the sample approaches the ultimate state with constant void ratio and follows a periodic strain path. The internal structure anisotropy has been quantified in terms of the contact-based fabric tensor. Rotation of principal stress axes densifies the packing, and leads to the increase in coordination numbers. A cyclic rotation in material anisotropy has been observed. The larger the stress ratio, the structure becomes more anisotropic. A larger fabric trajectory suggests more significant structure re-organization when rotating and explains the occurrence of more significant strain rate. The trajectory of the contact-normal based fabric is not centered in the origin, due to the anisotropy in particle orientation generated during sample generation which is persistent throughout the shearing process. The sample sheared at a lower intermediate principal stress ratio (b=0.0) (b=0.0) has been observed to approach a smaller strain trajectory as compared to the case b=0.5 b=0.5 , consistent with a smaller fabric trajectory and less significant structural re-organisation. It also experiences less volume contraction with the out-of plane strain component being dilative