A bond model for DEM simulation of cementitious materials and deformable structures

There is an increasing use of the discrete element method (DEM) to study cemented (e.g. concrete and rocks) and sintered particulate materials. The chief advantage of the DEM over continuum based techniques is that it does not make assumptions about how cracking and fragmentation initiate and propagate, since the DEM system is naturally discontinuous. The ability for the DEM to produce a realistic representation of a cemented granular material depends largely on the implementation of an inter-particle bonded contact model. This paper presents a new bonded contact model based on the Timoshenko beam theory which considers axial, shear and bending behaviour of the bond. The bond model was first verified by simulating both the bending and dynamic response of a simply supported beam. The loading response of a concrete cylinder was then investigated and compared with the Eurocode equation prediction. The results show significant potential for the new model to produce satisfactory predictions for cementitious materials. A unique feature of this model is that it can also be used to accurately represent many deformable structures such as frames and shells, so that both particles and structures or deformable boundaries can be described in the same DEM framework. 

Modeling Frequency-Dependent Boundaries as Digital Impedance Filters in FDTD and K-DWM Room Acoustics Simulations

A new method for modeling-frequency-dependent boundaries in finite-difference time-domain (FDTD) and Kirchhoff variable digital waveguide mesh (K-DWM) room acoustics simulations is presented. The proposed approach allows the direct incorporation of a digital impedance filter (DIF) in the Multidimensional (2D or 3D) FDTD boundary model of a locally reacting surface. An explicit boundary update equation is obtained by carefully constructing a Suitable recursive formulation. The method is analyzed in terms of pressure wave reflectance for different wall impedance filters and angles of incidence. Results obtained from numerical experiments confirm the high accuracy of the proposed digital impedance filter boundary model, the reflectance of which matches locally reacting surface (LRS) theory closely. Furthermore a numerical boundary analysis (NBA) formula is provided as a technique for an analytic evaluation of the numerical reflectance of the proposed digital impedance filter boundary formulation.

Parametric FIR Design of Propagation Loss Filters in Digital Waveguide String Models

One of the attractive features of sound synthesis by physical modeling is the potential to build acoustic-sounding digital instruments that offer more flexibility and different options in its design and control than their real-life counterparts. In order to develop such virtual-acoustic instruments, the models they are based on need to be fully parametric, i.e., all coefficients employed in the model are functions of physical parameters that are controlled either online or at the (offline) design stage. In this letter we show how propagation losses can be parametrically incorporated in digital waveguide string models with the use of zero-phase FIR filters. Starting from the simplest possible design in the form of a three-tap FIR filter, a higher-order FIR strategy is presented and discussed within the perspective of string sound synthesis with digital waveguide models.

A Physically-Constrained Source Model For FDTD Acoustic Simulation

The Finite Difference Time Domain (FDTD) method is becoming increasingly popular for room acoustics simulation. Yet, the literature on grid excitation methods is relatively sparse, and source functions are traditionally implemented in a hard or additive form
using arbitrarily-shaped functions which do not necessarily obey the physical laws of sound generation. In this paper we formulate
a source function based on a small pulsating sphere model. A physically plausible method to inject a source signal into the grid
is derived from first principles, resulting in a source with a near-flat spectrum that does not scatter incoming waves. In the final
discrete-time formulation, the source signal is the result of passing a Gaussian pulse through a digital filter simulating the dynamics of the pulsating sphere, hence facilitating a physically correct means to design source functions that generate a prescribed sound field.

Assessing modern ground survey methods and airborne laser scanning for digital terrain modelling:A case study from the Lake District, England

This paper compares the applicability of three ground survey methods for modelling terrain: one man electronic tachymetry (TPS), real time kinematic GPS (GPS), and terrestrial laser scanning (TLS). Vertical accuracy of digital terrain models (DTMs) derived from GPS, TLS and airborne laser scanning (ALS) data is assessed. Point elevations acquired by the four methods represent two sections of a mountainous area in Cumbria, England. They were chosen so that the presence of non-terrain features is constrained to the smallest amount. The vertical accuracy of the DTMs was addressed by subtracting each DTM from TPS point elevations. The error was assessed using exploratory measures including statistics, histograms, and normal probability plots. The results showed that the internal measurement accuracy of TPS, GPS, and TLS was below a centimetre. TPS and GPS can be considered equally applicable alternatives for sampling the terrain in areas accessible on foot. The highest DTM vertical accuracy was achieved with GPS data, both on sloped terrain (RMSE 0.16. m) and flat terrain (RMSE 0.02. m). TLS surveying was the most efficient overall but veracity of terrain representation was subject to dense vegetation cover. Therefore, the DTM accuracy was the lowest for the sloped area with dense bracken (RMSE 0.52. m) although it was the second highest on the flat unobscured terrain (RMSE 0.07. m). ALS data represented the sloped terrain more realistically (RMSE 0.23. m) than the TLS. However, due to a systematic bias identified on the flat terrain the DTM accuracy was the lowest (RMSE 0.29. m) which was above the level stated by the data provider. Error distribution models were more closely approximated by normal distribution defined using median and normalized median absolute deviation which supports the use of the robust measures in DEM error modelling and its propagation. © 2012 Elsevier Ltd.

A Hydrodynamic Modelling Framework for Strangford Lough Part 1: Tidal Model

Hydrodynamic models are a powerful tool that can be used by a wide range of end users to assist in predicting the effects of both physical and biological processes on local environmental conditions. This paper describes the development of a tidal model for Strangford Lough, Northern Ireland, a body of water renowned for the location of the first grid-connected tidal turbine, SeaGen, as well as the UK’s third Marine Nature Reserve. Using MIKE 21 modelling software, the development, calibration and performance of the modelare described in detail. Strangford Lough has a complex flow pattern with high flows through the Narrows (~3.5 m/s) linking the main body of the Lough to the Irish Sea and intricate flow patterns around the numerous islands. With the aid of good quality tidal and current data obtained throughout the Lough during the model development, the surface elevation and current magnitude between the observed and numerical model were almost identical with model skill >0.98 and >0.84 respectively. The applicability of the model is such that it can be used as an important tool for the prediction of important ecological processes as well as engineering applications within Strangford Lough.

Human body shadowing in cellular device-to-device communications: channel modeling using the shadowed κ−μ fading model

Using device-to-device communications as an underlay for cellular communications will provide an exciting opportunity to increase network capacity as well as improving spectral efficiency. The unique geometry of device-to-device links, where user equipment is often held or carried at low elevation and in close proximity to the human body, will mean that they are particularly susceptible to shadowing events caused not only by the local environment but also by the user's body. In this paper, the shadowed κ - μ fading model is proposed, which is capable of characterizing shadowed fading in wireless communication channels. In this model, the statistics of the received signal are manifested by the clustering of multipath components. Within each of these clusters, a dominant signal component with arbitrary power may exist. The resultant dominant signal component, which is formed by the phasor addition of these leading contributions, is assumed to follow a Nakagami- m distribution. The probability density function, moments, and the moment-generating function are also derived. The new model is then applied to device-to-device links operating at 868 MHz in an outdoor urban environment. It was found that shadowing of the resultant dominant component can vary significantly depending upon the position of the user equipment relative to the body and the link geometry. Overall, the shadowed κ - μ fading model is shown to provide a good fit to the field data as well as providing a useful insight into the characteristics of the received signal.

A DEM study of silo discharge of a cohesive solid

Bulk handling of powders and granular solids is common in many industries and often gives rise to handling difficulties especially when the material exhibits complex cohesive behaviour. For example, high storage stresses in a silo can lead to high cohesive strength of the stored solid, which may in turn cause blockages such as ratholing or arching near the outlet during discharge. This paper presents a Discrete Element Method study of discharge of a granular solid with varying levels of cohesion from a flat-bottomed silo. The DEM simulations were conducted using the commercial EDEM code with a recently developed DEM contact model for cohesive solids implemented through an API. The contact model is based on an elasto-plastic contact with adhesion and uses hysteretic non-linear loading and unloading paths to model the elastic-plastic contact deformation. The adhesion parameter is a function of the maximum contact overlap. The model has been shown to be able to predict the stress history dependent behaviour depicted by a flow function of the material. The effects of cohesion on the discharge rate and flow pattern in the silo are investigated. The predicted discharge rates are compared for the varying levels of cohesion and the effect of adhesion is evaluated. The ability of the contact model to qualitatively predict the phenomena that are present in the discharge of a silo has been shown with the salient feature of mixed flow from a flat bottomed hopper identified in the simulation.

A stacked Late Quaternary fluvio-periglacial sequence from the Axe valley, southern England with implications for landscape evolution and Palaeolithic archaeology

The current model of mid-latitude late Quaternary terrace sequences, is that they are uplift-driven but climatically controlled terrace staircases, relating to both regional-scale crustal and tectonic factors, and palaeohydrological variations forced by quasi-cyclic climatic conditions in the 100 K world (post Mid Pleistocene Transition). This model appears to hold for the majority of the river valleys draining into the English Channel which exhibit 8–15 terrace levels over approximately 60–100 m of altitudinal elevation. However, one valley, the Axe, has only one major morphological terrace and has long-been regarded as anomalous. This paper uses both conventional and novel stratigraphical methods (digital granulometry and terrestrial laser scanning) to show that this terrace is a stacked sedimentary sequence of 20–30 m thickness with a quasi-continuous (i.e. with hiatuses) pulsed, record of fluvial and periglacial sedimentation over at least the last 300–400 K yrs as determined principally by OSL dating of the upper two thirds of the sequence. Since uplift has been regional, there is no evidence of anomalous neotectonics, and climatic history must be comparable to the adjacent catchments (both of which have staircase sequences) a catchment-specific mechanism is required. The Axe is the only valley in North West Europe incised entirely into the near-horizontally bedded chert (crypto-crystalline quartz) and sand-rich Lower Cretaceous rocks creating a buried valley. Mapping of the valley slopes has identified many large landslide scars associated with past and present springs. It is proposed that these are thaw-slump scars and represent large hill-slope failures caused by Vauclausian water pressures and hydraulic fracturing of the chert during rapid permafrost melting. A simple 1D model of this thermokarstic process is used to explore this mechanism, and it is proposed that the resultant anomalously high input of chert and sand into the valley during terminations caused pulsed aggradation until the last termination. It is also proposed that interglacial and interstadial incision may have been prevented by the over-sized and interlocking nature of the sub-angular chert clasts until the Lateglacial when confinement of the river overcame this immobility threshold. One result of this hydrogeologically mediated valley evolution was to provide a sequence of proximal Palaeolithic archaeology over two MIS cycles. This study demonstrates that uplift tectonics and climate alone do not fully determine Quaternary valley evolution and that lithological and hydrogeological conditions are a fundamental cause of variation in terrestrial Quaternary records and landform evolution.

Spatio-Temporal Rich Model For Motion Vector Steganalysis

We propose a spatio-temporal rich model of motion vector planes as a part of a full steganalytic system against motion vector based steganography. Superior detection accuracy of the rich model over the previous methods has been lately demonstrated for digital images in both spatial and DCT domain. It has not been heretofore used for detection of motion vector steganography. We also introduced a transformation so as to extend the feature set with temporal residuals. We carried out the tests along with most recent motion vector steganalysis and steganography methods. Test results show that the proposed model delivers an outstanding performance compared to the previous methods.

Passive Intermodulation of Analog and Digital Signals on Transmission Lines with Distributed Nonlinearities: Modelling and Characterization

Passive intermodulation (PIM) often limits the performance of communication systems with analog and digitally-modulated signals and especially of systems supporting multiple carriers. Since the origins of the apparently multiple physical sources of nonlinearity causing PIM are not fully understood, the behavioral models are frequently used to describe the process of PIM generation. In this paper a polynomial model of memoryless nonlinearity is deduced from PIM measurements of a microstrip line with distributed nonlinearity with two-tone CW signals. The analytical model of nonlinearity is incorporated in Keysight Technology’s ADS simulator to evaluate the metrics of signal fidelity in the receive band for analog and digitally-modulated signals. PIM-induced distortion and cross-band interference with modulated signals are compared to those with two-tone CW signals. It is shown that conventional metrics can be applied to quantify the effect of distributed nonlinearities on signal fidelity. It is found that the two-tone CW test provides a worst-case estimate of cross-band interference for two-carrier modulated signals whereas with a three-carrier signal PIM interference in the receive band is noticeably overestimated. The simulated constellation diagrams for QPSK signals demonstrate that PIM interference exhibits the distinctive signatures of correlated distortion and this indicates that there are opportunities for mitigating PIM interference and that PIM interference cannot be treated as noise. One of the interesting results is that PIM distortion on a transmission line results in asymmetrical regrowth of output PIM interference for modulated signals.