Micromechanical analysis of cohesive granular materials using the discrete element method with an adhesive elasto-plastic contact model

An adhesive elasto-plastic contact model for the discrete element method with three dimensional non-spherical particles is proposed and investigated to achieve quantitative prediction of cohesive powder flowability. Simulations have been performed for uniaxial consolidation followed by unconfined compression to failure using this model. The model has been shown to be capable of predicting the experimental flow function (unconfined compressive strength vs. the prior consolidation stress) for a limestone powder which has been selected as a reference solid in the Europe wide PARDEM research network. Contact plasticity in the model is shown to affect the flowability significantly and is thus essential for producing satisfactory computations of the behaviour of a cohesive granular material. The model predicts a linear relationship between a normalized unconfined compressive strength and the product of coordination number and solid fraction. This linear relationship is in line with the Rumpf model for the tensile strength of particulate agglomerate. Even when the contact adhesion is forced to remain constant, the increasing unconfined strength arising from stress consolidation is still predicted, which has its origin in the contact plasticity leading to microstructural evolution of the coordination number. The filled porosity is predicted to increase as the contact adhesion increases. Under confined compression, the porosity reduces more gradually for the load-dependent adhesion compared to constant adhesion. It was found that the contribution of adhesive force to the limiting friction has a significant effect on the bulk unconfined strength. The results provide new insights and propose a micromechanical based measure for characterising the strength and flowability of cohesive granular materials. 

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.

Discrete Conditional Phase-type Model Utilising a Multiclass Support Vector Machine for the Prediction of Retinopathy of Prematurity

Retinopathy of prematurity (ROP) is a rare disease in which retinal blood vessels of premature infants fail to develop normally, and is one of the major causes of childhood blindness throughout the world. The Discrete Conditional Phase-type (DC-Ph) model consists of two components, the conditional component measuring the inter-relationships between covariates and the survival component which models the survival distribution using a Coxian phase-type distribution. This paper expands the DC-Ph models by introducing a support vector machine (SVM), in the role of the conditional component. The SVM is capable of classifying multiple outcomes and is used to identify the infant's risk of developing ROP. Class imbalance makes predicting rare events difficult. A new class decomposition technique, which deals with the problem of multiclass imbalance, is introduced. Based on the SVM classification, the length of stay in the neonatal ward is modelled using a 5, 8 or 9 phase Coxian distribution.

Discrete-Time Conserved Quantities for Damped Oscillators

Numerical sound synthesis is often carried out using the finite difference time domain method. In order to analyse the stability of the derived models, energy methods can be used for both linear and nonlinear settings. For Hamiltonian systems the existence of a conserved numerical energy-like quantity can be used to guarantee the stability of the simulations. In this paper it is shown how to derive similar discrete conservation laws in cases where energy is dissipated due to friction or in the presence of an energy source due to an external force. A damped harmonic oscillator (for which an analytic solution is available) is used to present the proposed methodology. After showing how to arrive at a conserved quantity, the simulation of a nonlinear single reed shows an example of an application in the context of musical acoustics.

Geotechnical monitoring using a combination of LIDAR, real-time dGPS, and electrical tomography to assess geotechnical risk of a coastal landslide.

The Antrim Coast Road stretching from the seaport of Larne in the East of Northern Ireland has a well-deserved reputation for being one of the most spectacular roads in Europe (Day, 2006). However the problematic geology; Jurassic Lias Clay and Triassic Mudstone overlain by Cretaceous Limestone and Tertiary Basalt, and environmental variables result in frequent instances of slope instability manifested in both shallow debris flows and occasional massive rotational movements, creating a geotechnical risk to this highway. This paper describes how a variety of techniques are being used to both assess instability and monitor movement of these active slopes near one site at Straidkilly Point, Glenarm. An in-depth understanding of the geology was obtained via boreholes, resistivity surveys and laboratory testing. Environmental variables recorded by an on-site weather station were correlated with measured pore water pressure and soil moisture infiltration data. Terrestrial LiDAR (TLS), with surveys carried out on a bi-monthly basis allowed for the generation of Digital Elevation Models (DEMs) of difference, highlighting areas of recent movement, accumulation and depletion. Morphology parameters were generated from the DEMs and include slope, curvature and multiple measures of roughness. Changes in the structure of the slope coupled with morphological parameters were characterised and linked to progressive failures from the temporal monitoring. In addition to TLS monitoring, Aerial LiDAR datasets were used for the spatio-morphological characterisation of the slope on a macro scale. A Differential Global Positioning System (dGPS) was also deployed on site to provide a real-time warning system for gross movements, which were also correlated with environmental conditions. Frequent electrical resistivity tomography (ERT) surveys were also implemented to provide a better understanding of long-term changes in soil moisture and help to define the complex geology. The paper describes how the data obtained via a diverse range of methods has been combined to facilitate a more informed management regime of geotechnical risk by the Northern Ireland Roads Service.

Using a Discrete Hidden Markov Model Kernel for lip-based biometric identification

In this paper, a novel and effective lip-based biometric identification approach with the Discrete Hidden Markov Model Kernel (DHMMK) is developed. Lips are described by shape features (both geometrical and sequential) on two different grid layouts: rectangular and polar. These features are then specifically modeled by a DHMMK, and learnt by a support vector machine classifier. Our experiments are carried out in a ten-fold cross validation fashion on three different datasets, GPDS-ULPGC Face Dataset, PIE Face Dataset and RaFD Face Dataset. Results show that our approach has achieved an average classification accuracy of 99.8%, 97.13%, and 98.10%, using only two training images per class, on these three datasets, respectively. Our comparative studies further show that the DHMMK achieved a 53% improvement against the baseline HMM approach. The comparative ROC curves also confirm the efficacy of the proposed lip contour based biometrics learned by DHMMK. We also show that the performance of linear and RBF SVM is comparable under the frame work of DHMMK.

Passive Intermodulation in Distributed Circuits with Cascaded Discrete Nonlinearities

The principle aspects of passive intermodulation (PIM) characterisation in distributed printed circuits with cascaded lumped nonlinearities are presented. Mechanisms of PIM generations have been investigated experimentally and modelled using the formalism of X-parameters. The devised equivalent circuit models are applied to the analysis of microstrip lines with distributed and cascaded lumped sources of nonlinearity. The dynamic measurements have revealed that PIM generation rates in straight and meandered microstrip lines differ and significantly deviate from those expected for the respective discrete sources of nonlinearity. The obtained results indicate that multiple physical sources of nonlinearity contribute to PIM generation in printed circuits. Finally, it is demonstrated that the electrical discontinuities can have significant effect on the overall PIM response of the distributed passive circuits and cause PIM product leakage and parasitic coupling between isolated circuit elements.

Comparisons of anterior segment biometry between Chinese and Caucasians using anterior segment optical coherence tomography.

PURPOSE: To compare anterior segment parameters between eyes of Chinese and Caucasians using anterior segment optical coherence tomography and to evaluate the association between these parameters and anterior chamber angle width between the two ethnic groups. METHODS: 60 Chinese and 60 Caucasians, 30 with open angles and 30 with narrow angles (defined as Shaffer grade < or =2 in > or =3 quadrants during dark room gonioscopy) in each group, were consecutively enrolled. One eye of each subject was randomly selected for imaging in a completely darkened room. Measurements, including anterior chamber depth (ACD), scleral spur-to-scleral spur distance (anterior chamber width (ACW)), anterior chamber angle width, iris convexity and iris thickness, were compared between the groups. The associations between angle opening distance and biometric measurements were evaluated with univariate and multivariate regression analyses. RESULTS: There were no differences in age, axial length, anterior chamber angle measurements, pupil diameter and iris convexity between Chinese and Caucasians in both open-angle and narrow-angle groups. However, the ACD and ACW were smaller and the iris was thicker in Chinese. In the multivariate analysis, the ACD was the most influential biometric parameter for angle opening distance in both Chinese and Caucasians. After adjusting the effects of axial length, age and sex, ACD and ACW were significantly smaller in Chinese. CONCLUSIONS: Chinese eyes had smaller ACD, smaller ACW and greater iris thickness than Caucasians. ACD was the most influential parameter in determining the angle width in both ethnic groups.

Comparison of anterior ocular biometry between optical low-coherence reflectometry and anterior segment optical coherence tomography in an adult Chinese population.

PURPOSE: To evaluate the agreement between optical low-coherence reflectometry (OLCR) and anterior segment optical coherence tomography (AS-OCT) for biometry of the anterior segment. SETTING: State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China. DESIGN: Evaluation of diagnostic technology. METHODS: A series of OLCR (Lenstar LS 900) and AS-OCT measurements of the anterior segment were taken for consecutive subjects aged 35 years and older in a population-based study. The differences and correlations between the 2 methods of ocular biometry were assessed. Agreement was calculated as the 95% limits of agreement (LoA). RESULTS: The mean age of the 776 subjects was 55.2 years ± 12.0 (SD); 54.6% were women. The mean central corneal thickness (CCT) was smaller with OLCR than with AS-OCT (537.84 ± 31.46 μm versus 559.39 ± 32.02 μm) as was anterior chamber depth (ACD) (2.60 ± 0.37 mm versus 2.72 ± 0.37 mm) and anterior chamber width (ACW) (11.76 ± 0.47 mm versus 12.04 ± 0.55 mm) (all P<.001). The 95% LoA between the 2 instruments were -44.80 to 1.71 μm for CCT, -0.17 to -0.06 mm for ACD, and -1.28 to 0.72 mm for ACW. CONCLUSION: Optical low-coherence reflectometry and AS-OCT yielded potentially interchangeable ACD measurements, while the CCT and ACW measurements acquired by the 2 devices showed clinically significant differences.

Calculus of variations on time scales and discrete fractional calculus

Estudamos problemas do cálculo das variações e controlo óptimo no contexto das escalas temporais. Especificamente, obtemos condições necessárias de optimalidade do tipo de Euler–Lagrange tanto para lagrangianos dependendo de derivadas delta de ordem superior como para problemas isoperimétricos. Desenvolvemos também alguns métodos directos que permitem resolver determinadas classes de problemas variacionais através de desigualdades em escalas temporais. No último capítulo apresentamos operadores de diferença fraccionários e propomos um novo cálculo das variações fraccionário em tempo discreto. Obtemos as correspondentes condições necessárias de Euler– Lagrange e Legendre, ilustrando depois a teoria com alguns exemplos.

Development of a single photon counting computed tomography system using MPGDs

The development of computed tomography systems with energy resolving detectors is a current challenge in medical physics and biomedical engineering. A computed tomography system of this kind allows getting complementary informations relatively to conventional systems, that can help the medical diagnosis, being of great interest in medicine. The work described in this thesis is related to the development of a computed tomography system using micropattern gaseous detectors, which allow storing, simultaneously, information about the interaction position and the energy of each single photon that interacts with the detector. This kind of detectors has other advantages concerning the cost and characteristics of operation when compared with solid state detectors. Tomographic acquisitions were performed using a MicroHole & Strip Plate based detector, which allowed reconstructing cross-sectional images using energy windows, applying the energy weighting technique and performing multi-slice and tri-dimensional reconstructions. The contrast-to-noise ratio was improved by 31% by applying the energy weighting technique, comparing with the corresponding image obtained with the current medical systems. A prototype of a computed tomography with flexibility to change the detector was developed, making it possible to apply different detectors based on Thick-COBRA. Several images acquired with these detectors are presented and demonstrate their applicability in X-ray imaging. When operating in NeCH4, the detector allowed a charge gain of 8 104, an energy resolution of 20% (full width at half maximum at 8 keV), a count rate of 1 106 Hz/mm2, a very stable operation (gain fluctuations below 5%) and a spacial resolution of 1.2 mm for an energy photon of 3.6 keV. Operating the detector in pure Kr allowed increasing the detection efficiency and achieving a charge gain of 2 104, an energy resolution of 32% (full width at half maximum at 22 keV), a count rate of 1 105 Hz/mm2, very stable operation and a spatial resolution of 500 m. The software already existing in the group was improved and tools to correct geometric misalignments of the system were also developed. The reconstructions obtained after geometrical correction are free of artefacts due to the referred misalignments.

An autonomous system for ocean acoustic tomography

Vertical line arrays (VLA) are a widely used apparatus in underwater acoustics with applications in sonar prediction, underwater communications and acoustic tomography, among others. Recent developments in digital electronics and communications allow for off-the-shelf development of VLA systems, with a large number of embedded acoustic and non-acoustic sensors able to fulfill application requirements, as opposed to single or few receiver configurations available until only a few years ago. Very often, the flexibility in water column sampling is achieved by splitting the VLA into modules that can be assembled according to the application. Such systems can be deployed and recovered from small vessels with a shorthanded crew, and make it possible for research labs with reduced budgets and operational means (ships and manpower) to gain control over the whole development process, from data acquisition to post-processing.