Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

The analytical description and experiments of the optical bottle generated by an axicon and a lens

Based on the interferential theory, we deduce a new type of analytic expression suitable for describing the evolutions of the optical bottle beam generated from the axicon-lens optical system illuminated by the Gaussian beam for the first time. The theory does not use much approximation in the process of mathematical analysis and can better illustrate the optical bottle beam evolutions at any positions. With the derived expression, the three-dimensional (3D) longitudinal and transverse intensity profiles of the optical bottle beam are simulated numerically. The numerical calculations have been confirmed by the experimental results.

One-dimensional potential for image-potential states on graphene

In the framework of dielectric theory, the static non-local self-energy of an electron near an ultra-thin polarizable layer has been calculated and applied to study binding energies of image-potential states near free-standing graphene. The corresponding series of eigenvalues and eigenfunctions have been obtained by numerically solving the one-dimensional Schrodinger equation. The imagepotential state wave functions accumulate most of their probability outside the slab. We find that the random phase approximation (RPA) for the nonlocal dielectric function yields a superior description for the potential inside the slab, but a simple Fermi-Thomas theory can be used to get a reasonable quasi-analytical approximation to the full RPA result that can be computed very economically. Binding energies of the image-potential states follow a pattern close to the Rydberg series for a perfect metal with the addition of intermediate states due to the added symmetry of the potential. The formalism only requires a minimal set of free parameters: the slab width and the electronic density. The theoretical calculations are compared with experimental results for the work function and image-potential states obtained by two-photon photoemission.

Analytical model for the 3D-RESURF effect

This paper presents an analytical model for the determination of the basic breakdown properties of three-dimensional (3D)-RESURF/CoolMOS/super junction type structures. To account for the two-dimensional (2D) effect of the 3D-RESURF action, 2D models of the electric field distribution are developed. Based on these, expressions are derived for the breakdown voltage as a function of doping concentration and physical dimensions. In addition to cases where the drift regions are fully depleted, the model developed is also applicable to situations involving drift regions which are almost depleted. Accuracy of the analytical approach is verified by comparison with numerical results obtained from the MEDICI device simulator.

A three-dimensional tunnel model for calculation of train-induced ground vibration

The frequency range of interest for ground vibration from underground urban railways is approximately 20 to 100 Hz. For typical soils, the wavelengths of ground vibration in this frequency range are of the order of the spacing of train axles, the tunnel diameter and the distance from the tunnel to nearby building foundations. For accurate modelling, the interactions between these entities therefore have to be taken into account. This paper describes an analytical three-dimensional model for the dynamics of a deep underground railway tunnel of circular cross-section. The tunnel is conceptualised as an infinitely long, thin cylindrical shell surrounded by soil of infinite radial extent. The soil is modelled by means of the wave equations for an elastic continuum. The coupled problem is solved in the frequency domain by Fourier decomposition into ring modes circumferentially and a Fourier transform into the wavenumber domain longitudinally. Numerical results for the tunnel and soil responses due to a normal point load applied to the tunnel invert are presented. The tunnel model is suitable for use in combination with track models to calculate the ground vibration due to excitation by running trains and to evaluate different track configurations. © 2006 Elsevier Ltd. All rights reserved.

Behaviour of unsteady transonic shock / boundary layer interactions with three-dimensional effects

An experimental investigation into the response of transonic SBLIs to periodic down-stream pressure perturbations in a parallel walled duct has been conducted. Tests have been carried out with a shock strength of M ∞ = 1.5 for pressure perturbation frequencies in the range 16-90 Hz. Analysis of the steady interaction at M∞ = 1.5 has also been made. The principle measurement techniques were high speed schlieren photography and laser Doppler anemometry. The structure of the steady SBLI was found to be highly three-dimensional, with large corner flows and sidewall SBLIs. These aspects are thought to influence the upstream transmission of pressure information through the interaction by affecting the post-shock flow field, including the extent of regions of secondary supersonic flow. At low frequency, the dynamics of shock motion can be predicted using an inviscid analytical model. At increased frequencies, viscous effects become significant and the shock exhibits unexpected dynamic behaviour, due to a phase lag between the upstream transmission of pressure information in the core flow and in the viscous boundary layers. Flow control in the form of micro-vane vortex generators was found to have a small impact on shock dynamics, due to the effect it had on the post-shock flow field outside the viscous boundary layer region. The relationship between inviscid and viscous effects is developed and potential destabilising mechanisms for SBLIs in practical applications are suggested. Copyright © 2009 by Paul Bruce and Holger Babinsky.

Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Decay of vibrations along deepwater mooring lines

Model tests for global design verification of deepwater floating structures cannot be made at reasonable scales. An overview of recent research efforts to tackle this challenge is given first, introducing the concept of line truncation techniques. In such a method the upper sections of each line are modelled in detail, capturing the wave action zone and all coupling effects with the vessel. These terminate to an approximate analytical model, that aims to simulate the remainder of the line. The rationale for this is that in deep water the transverse elastic waves of a line are likely to decay before they are reflected at the seabed. The focus of this paper is the verification of this rationale and the ongoing work, which is considering ways to produce a truncation model. Transverse dynamics of a mooring line are modelled using the equations of motion of an inextensible taut string, submerged in still water, one end fixed at the bottom the other assumed to follow the vessel response, which can be harmonic or random. Nonlinear hydrodynamic damping is included; bending and VIV effects are neglected. A dimensional analysis, supported by exact benchmark numerical solutions, has shown that it is possible to produce a universal curve for the decay of transverse vibrations along the line, which is suitable for any kind of line with any top motion. This has a significant engineering benefit, allowing for a rapid assessment of line dynamics - it is very useful in deciding whether a truncated line model is appropriate, and if so, at which point truncation might be applied. Initial efforts in developing a truncated model show that a linearized numerical solution in the frequency domain matches very closely the exact benchmark. Copyright © 2011 by ASME.

High order schemes on three-dimensional general polyhedral meshes - Application to the level set method

In this article, we detail the methodology developed to construct arbitrarily high order schemes - linear and WENO - on 3D mixed-element unstructured meshes made up of general convex polyhedral elements. The approach is tailored specifically for the solution of scalar level set equations for application to incompressible two-phase flow problems. The construction of WENO schemes on 3D unstructured meshes is notoriously difficult, as it involves a much higher level of complexity than 2D approaches. This due to the multiplicity of geometrical considerations introduced by the extra dimension, especially on mixed-element meshes. Therefore, we have specifically developed a number of algorithms to handle mixed-element meshes composed of convex polyhedra with convex polygonal faces. The contribution of this work concerns several areas of interest: the formulation of an improved methodology in 3D, the minimisation of computational runtime in the implementation through the maximum use of pre-processing operations, the generation of novel methods to handle complex 3D mixed-element meshes and finally the application of the method to the transport of a scalar level set. © 2012 Global-Science Press.

Plastic limit solution for the response of plate anchors under combined translational and torsional undrained loading

Plate anchors are increasingly being used to moor large floating offshore structures in deep and ultradeep water. These facilities impart substantial vertical uplift loading to plate anchors. However, extreme operating conditions such as hurricane loading often result in partial system failures, with significant change in the orientation of the remaining intact mooring lines. The purpose of this study is to investigate the undrained pure translational (parallel to plate) and torsional bearing capacity of anchor plates idealized as square and rectangular shaped plates. Moreover, the interaction response of plate anchors under combined translational and torsional loading is studied using a modified plastic limit analysis (PLA) approach. The previous PLA formulation which did not account for shear-normal force interaction on the vertical end faces of the plate provides an exact solution to the idealized problem of an infinitely thin plate but only an approximate solution to the problem of a plate of finite thickness. This is also confirmed by the three-dimensional finite element (FE) results, since the PLA values exceed FE results as the thickness of the plate increases. By incorporating the shear-normal interaction relationship in the modified solution, the torsional bearing capacity factors, as well as the plate interaction responses are enhanced as they show satisfactory agreement with the FE results. The interaction relationship is then obtained for square and rectangular plates of different aspect ratios and thicknesses. The new interaction relationships could also be used as an associated plastic failure locus for combined shear and torsional loading to predict plastic displacements and rotations in translational and torsional loading modes as well. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).

Wave propagation in nonlinear one-dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one-dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross-correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross-correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loop under steady-state condition. Copyright © 2008 John Wiley & Sons, Ltd.

Cabaret finite-difference schemes for the one-dimensional Euler equations

In the present paper we consider second order compact upwind schemes with a space split time derivative (CABARET) applied to one-dimensional compressible gas flows. As opposed to the conventional approach associated with incorporating adjacent space cells we use information from adjacent time layer to improve the solution accuracy. Taking the first order Roe scheme as the basis we develop a few higher (i.e. second within regions of smooth solutions) order accurate difference schemes. One of them (CABARET3) is formulated in a two-time-layer form, which makes it most simple and robust. Supersonic and subsonic shock-tube tests are used to compare the new schemes with several well-known second-order TVD schemes. In particular, it is shown that CABARET3 is notably more accurate than the standard second-order Roe scheme with MUSCL flux splitting.

Binder free three-dimensional sulphur/few-layer graphene foam cathode with enhanced high-rate capability for rechargeable lithium sulphur batteries.

A novel ultra-lightweight three-dimensional (3-D) cathode system for lithium sulphur (Li-S) batteries has been synthesised by loading sulphur on to an interconnected 3-D network of few-layered graphene (FLG) via a sulphur solution infiltration method. A free-standing FLG monolithic network foam was formed as a negative of a Ni metallic foam template by CVD followed by etching away of Ni. The FLG foam offers excellent electrical conductivity, an appropriate hierarchical pore structure for containing the electro-active sulphur and facilitates rapid electron/ion transport. This cathode system does not require any additional binding agents, conductive additives or a separate metallic current collector thus decreasing the weight of the cathode by typically ∼20-30 wt%. A Li-S battery with the sulphur-FLG foam cathode shows good electrochemical stability and high rate discharge capacity retention for up to 400 discharge/charge cycles at a high current density of 3200 mA g(-1). Even after 400 cycles the capacity decay is only ∼0.064% per cycle relative to the early (e.g. the 5th cycle) discharge capacity, while yielding an average columbic efficiency of ∼96.2%. Our results indicate the potential suitability of graphene foam for efficient, ultra-light and high-performance batteries.

Photoluminescence investigation of two-dimensional electron gas in an undoped AlxGa1-xN/GaN heterostructure

Low-temperature photoluminescence measurement is performed on an undoped AlxGa1-xN/GaN heterostructure. Temperature-dependent Hall mobility confirms the formation of two-dimensional electron gas (2DEG) near the heterointerface. A weak photoluminescence (PL) peak with the energy of similar to 79meV lower than the free exciton (FE) emission of bulk GaN is related to the radiative recombination between electrons confined in the triangular well and the holes near the flat-band region of GaN. Its identification is supported by the solution of coupled one-dimensional Poisson and Schrodinger equations. When the temperature increases, the red shift of the 2DEG related emission peak is slower than that of the FE peak. The enhanced screening effect coming from the increasing 2DEG concentration and the varying electron distribution at two lowest subbands as a function of temperature account for such behaviour.

Phylogeny of SARS-CoV based on high-dimensional information geometry

We proposed a novel methodology, which firstly, extracting features from species' complete genome data, using k-tuple, followed by studying the evolutionary relationship between SARS-CoV and other coronavirus species using the method, called "High-dimensional information geometry". We also used the mothod, namely "caculating of Minimum Spanning Tree", to construct the Phyligenetic tree of the coronavirus. From construction of the unrooted phylogenetic tree, we found out that the evolution distance between SARS-CoV and other coronavirus species is comparatively far. The tree accurately rebuilt the three groups of other coronavirus. We also validated the assertion from other literatures that SARS-CoV is similar to the coronavirus species in Group I.