G. Marketos

Bio: G. Marketos is an academic researcher from Utrecht University. The author has contributed to research in topics: Discrete element method & Granular material. The author has an hindex of 12, co-authored 25 publications receiving 451 citations. Previous affiliations of G. Marketos include National Oceanography Centre, Southampton & University of Cambridge.

Quantifying the extent of crushing in granular materials: A probability-based predictive method

Abstract: Grain crushing is one of the micromechanisms that governs the stress–strain behaviour of a granular material, and also its permeability by altering the grain size distribution. It is therefore advantageous to be able to predict the point of onset of crushing and to quantify the subsequent evolution of crushing. This paper uses the data of Discrete Element Method (DEM) simulations to inform a statistical model of granular crushing. Distributions of normalised contact forces are first obtained. If the statistical distribution of the crushing strength of the grains is then known, the onset of crushing within an assembly of grains should be predictable. Two different cases, one in which grain strength was statistically independent of grain size and one showing an arbitrary trend, were used to compare with DEM results and so confirm the validity of the statistical method.

Flat boundaries and their effect on sand testing

Abstract: A study of the effect of the use of flat boundaries on the stressing of a sample of an idealized granular material with no applied shear is presented. Discrete element method (DEM) data of 1D compression were analysed and the local strain field inside the sample was investigated as the sample was stressed. A best-fit strain was seen to best describe the material behaviour free from boundary effects. The individual particle displacements were probed, providing insight into the behaviour of particles adjacent to the boundaries. In addition, the porosity and force distribution inside the sample were observed, allowing for estimates of the width of a boundary region to be made. This region, non-representative of far-field material behaviour, will affect the behaviour of a granular sample in DEM or laboratory tests, with local porosity differences leading to a change in the transport properties of the sample, and force distribution changes leading to a bias in the location of grain cracking or crushing events for sufficiently high stress levels. Nevertheless, the largest effect of the boundary region was a severe underestimation of the stiffness of a granular material.

Micromechanics of seismic wave propagation in granular materials

Abstract: In this study experimental data on a model soil in a cubical cell are compared with both discrete element (DEM) simulations and continuum analyses. The experiments and simulations used point source transmitters and receivers to evaluate the shear and compression wave velocities of the samples, from which some of the elastic moduli can be deduced. Complex responses to perturbations generated by the bender/extender piezoceramic elements in the experiments were compared to those found by the controlled movement of the particles in the DEM simulations. The generally satisfactory agreement between experimental observations and DEM simulations can be seen as a validation and support the use of DEM to investigate the influence of grain interaction on wave propagation. Frequency domain analyses that considered filtering of the higher frequency components of the inserted signal, the ratio of the input and received signals in the frequency domain and sample resonance provided useful insight into the system response. Frequency domain analysis and analytical continuum solutions for cube vibration show that the testing configuration excited some, but not all, of the system’s resonant frequencies. The particle scale data available from DEM enabled analysis of the energy dissipation during propagation of the wave. Frequency domain analysis at the particle scale revealed that the higher frequency content reduces with increasing distance from the point of excitation.

Two-dimensional discrete element modelling of bender element tests on an idealised granular material

Abstract: The small-strain (elastic) shear stiffness of soil is an important parameter in geotechnics. It is required as an input parameter to predict deformations and to carry out site response analysis to predict levels of shaking during earthquakes. Bender element testing is often used in experimental soil mechanics to determine the shear (S-) wave velocity in a given soil and hence the shear stiffness. In a bender element test a small perturbation is input at a point source and the propagation of the perturbation through the system is measured at a single measurement point. The mechanics and dynamics of the system response are non-trivial, complicating interpretation of the measured signal. This paper presents the results of a series of discrete element method (DEM) simulations of bender element tests on a simple, idealised granular material. DEM simulations provide the opportunity to study the mechanics of this testing approach in detail. The DEM model is shown to be capable of capturing features of the system response that had previously been identified in continuum-type analyses of the system. The propagation of the wave through the sample can be monitored at the particle-scale in the DEM simulation. In particular, the particle velocity data indicated the migration of a central S-wave accompanied by P-waves moving along the sides of the sample. The elastic stiffness of the system was compared with the stiffness calculated using different approaches to interpreting bender element test data. An approach based upon direct decomposition of the signal using a fast-Fourier transform yielded the most accurate results.

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

The Quest for the Africa-Eurasia plate boundary West of the Strait of Gibraltar

Abstract: The missing link in the plate boundary between Eurasia and Africa in the central Atlantic is presented and discussed. A set of almost linear and sub parallel dextral strike–slip faults, the SWIM1 Faults, that form a narrow band of deformation over a length of 600 km coincident with a small circle centred on the pole of rotation of Africa with respect to Eurasia, was mapped using a new swath bathymetry compilation available in the area offshore SW Portugal. These faults connect the Gloria Fault to the Rif–Tell Fault Zone, two segments of the plate boundary between Africa and Eurasia. The SWIM faults cut across the Gulf of Cadiz, in the Atlantic Ocean, where the 1755 Great Lisbon earthquake, M ~ 8.5–8.7, and tsunami were generated, providing a new insight on its source location.

A DEM model for soft and hard rocks: Role of grain interlocking on strength

Abstract: The Discrete Element Method (DEM) is increasingly used to simulate the behavior of rock. Despite their intrinsic capability to model fracture initiation and propagation starting from simple interaction laws, classical DEM formulations using spherical discrete elements suffer from an intrinsic limitation to properly simulate brittle rock behavior characterized by high values of UCS/TS ratio associated with non-linear failure envelopes, as observed for hard rock like granite. The present paper shows that the increase of the interaction range between the spherical discrete elements, which increases locally the density of interaction forces (or interparticle bonds), can overcome this limitation. It is argued that this solution represents a way to implicitly take into account the degree of interlocking associated to the microstructural complexity of rock. It is thus shown that increasing the degree of interlocking between the discrete elements which represent the rock medium, in addition to enhancing the UCS/TS ratio, results in a non-linear failure envelop characteristic of low porous rocks. This approach improves significantly the potential and predictive capabilities of the DEM for rock modeling purpose. A special emphasis is put on the model ability to capture the fundamental characteristics of brittle rocks in terms of fracture initiation and propagation. The model can reproduce an essential component of brittle rock failure, that is, cohesion weakening and frictional strengthening as a function of rock damage or plastic strain. Based on model predictions, it is finally discussed that frictional strengthening may be at the origin of the brittle ductile transition occurring at high confining pressures.

Fundamentals of soil stabilization

Abstract: Clayey soils are usually stiff when they are dry and give up their stiffness as they become saturated. Soft clays are associated with low compressive strength and excessive settlement. This reduction in strength due to moisture leads to severe damages to buildings and foundations. The soil behavior can be a challenge to the designer build infrastructure plans to on clay deposits. The damage due to the expansive soils every year is expected to be $1 billion in the USA, £150 million in the UK, and many billions of pounds worldwide. The damages associated with expansive soils are not because of the lack of inadequate engineering solutions but to the failure to identify the existence and magnitude of expansion of these soils in the early stage of project planning. One of the methods for soil improvement is that the problematic soil is replaced by suitable soil. The high cost involved in this method has led researchers to identify alternative methods, and soil stabilization with different additives is one of those methods. Recently, modern scientific techniques of soil stabilization are on offer for this purpose. Stabilized soil is a composite material that is obtained from the combination and optimization of properties of constituent materials. Adding cementing agents such as lime, cement and industrial byproducts like fly ash and slag, with soil results in improved geotechnical properties. However, during the past few decades, a number of cases have been reported where sulfate-rich soils stabilized by cement or lime underwent a significant amount of heave leading to pavement failure. This research paper addressed the some fundamental and success soil improvement that used in civil engineering field.