Compositional data analysis in archaeometry

We shall call an n × p data matrix fully-compositional if the rows sum to a constant, and sub-compositional if the variables are a subset of a fully-compositional data set1. Such data occur widely in archaeometry, where it is common to determine the chemical composition of ceramic, glass, metal or other artefacts using techniques such as neutron activation analysis (NAA), inductively coupled plasma spectroscopy (ICPS), X-ray fluorescence analysis (XRF) etc. Interest often centres on whether there are distinct chemical groups within the data and whether, for example, these can be associated with different origins or manufacturing technologiesGeologische Vereinigung; Universitat de Barcelona, Equip de Recerca Arqueomètrica; Institut d’Estadística de Catalunya; International Association for Mathematical Geology; Patronat de l’Escola Politècnica Superior de la Universitat de Girona; Fundació privada: Girona, Universitat i Futur

Compositional data analysis of archaeological glass: problems and possible solutions

Presentation in CODAWORK'03, session 4: Applications to archeometryGeologische Vereinigung; Universitat de Barcelona, Equip de Recerca Arqueomètrica; Institut d’Estadística de Catalunya; International Association for Mathematical Geology; Patronat de l’Escola Politècnica Superior de la Universitat de Girona; Fundació privada: Girona, Universitat i Futur

Adaptive mesh optimization for simulation of immiscible viscous fingering

Viscous fingering can be a major concern when waterflooding heavy-oil reservoirs. Most commercial reservoir simulators use low-order finite-volume/-difference methods on structured grids to resolve this phenomenon. However, this approach suffers from a significant numerical-dispersion error because of insufficient mesh resolution, which smears out some important features of the flow. We simulate immiscible incompressible two-phase displacements and propose the use of unstructured control-volume finite-element (CVFE) methods for capturing viscous fingering in porous media. Our approach uses anisotropic mesh adaptation where the mesh resolution is optimized on the basis of the evolving features of flow. The adaptive algorithm uses a metric tensor field dependent on solution-interpolation-error estimates to locally control the size and shape of elements in the metric. The mesh optimization generates an unstructured finer mesh in areas of the domain where flow properties change more quickly and a coarser mesh in other regions where properties do not vary so rapidly. We analyze the computational cost of mesh adaptivity on unstructured mesh and compare its results with those obtained by a commercial reservoir simulator on the basis of the finite-volume methods

A Discontinuous Control Volume Finite Element Method for Multi-Phase Flow in Heterogeneous Porous Media

We present a new, high-order, control-volume-finite-element (CVFE) method for multiphase porous media flow with discontinuous 1st-order representation for pressure and discontinuous 2nd-order representation for velocity. The method has been implemented using unstructured tetrahedral meshes to discretize space. The method locally and globally conserves mass. However, unlike conventional CVFE formulations, the method presented here does not require the use of control volumes (CVs) that span the boundaries between domains with differing material properties. We demonstrate that the approach accurately preserves discontinuous saturation changes caused by permeability variations across such boundaries, allowing efficient simulation of flow in highly heterogeneous models. Moreover, accurate solutions are obtained at significantly lower computational cost than using conventional CVFE methods. We resolve a long-standing problem associated with the use of classical CVFE methods to model flow in highly heterogeneous porous media

A force-balanced control volume finite element method for multi-phase porous media flow modelling

A novel method for simulating multi-phase flow in porous media is presented. The approach is based on a control volume finite element mixed formulation and new force-balanced finite element pairs. The novelty of the method lies in: (a) permitting both continuous and discontinuous description of pressure and saturation between elements; (b) the use of arbitrarily high-order polynomial representation for pressure and velocity and (c) the use of high-order flux-limited methods in space and to time avoid introducing non-physical oscillations while achieving high-order accuracy where and when possible. The model is initially validated for two-phase flow. Results are in good agreement with analytically obtained solutions and experimental results. The potential of this method is demonstrated by simulating flow in a realistic geometry composed of highly permeable meandering channels

Predicting live birth, preterm and low birth weight infant after in-vitro fertilisation: a prospective study of 144018 treatment cycles

Background The extent to which baseline couple characteristics affect the probability of live birth and adverse perinatal outcomes after assisted conception is unknown. Methods and Findings We utilised the Human Fertilisation and Embryology Authority database to examine the predictors of live birth in all in vitro fertilisation (IVF) cycles undertaken in the UK between 2003 and 2007 (n = 144,018). We examined the potential clinical utility of a validated model that pre-dated the introduction of intracytoplasmic sperm injection (ICSI) as compared to a novel model. For those treatment cycles that resulted in a live singleton birth (n = 24,226), we determined the associates of potential risk factors with preterm birth, low birth weight, and macrosomia. The overall rate of at least one live birth was 23.4 per 100 cycles (95% confidence interval [CI] 23.2–23.7). In multivariable models the odds of at least one live birth decreased with increasing maternal age, increasing duration of infertility, a greater number of previously unsuccessful IVF treatments, use of own oocytes, necessity for a second or third treatment cycle, or if it was not unexplained infertility. The association of own versus donor oocyte with reduced odds of live birth strengthened with increasing age of the mother. A previous IVF live birth increased the odds of future success (OR 1.58, 95% CI 1.46–1.71) more than that of a previous spontaneous live birth (OR 1.19, 95% CI 0.99–1.24); p-value for difference in estimate <0.001. Use of ICSI increased the odds of live birth, and male causes of infertility were associated with reduced odds of live birth only in couples who had not received ICSI. Prediction of live birth was feasible with moderate discrimination and excellent calibration; calibration was markedly improved in the novel compared to the established model. Preterm birth and low birth weight were increased if oocyte donation was required and ICSI was not used. Risk of macrosomia increased with advancing maternal age and a history of previous live births. Infertility due to cervical problems was associated with increased odds of all three outcomes—preterm birth, low birth weight, and macrosomia. Conclusions Pending external validation, our results show that couple- and treatment-specific factors can be used to provide infertile couples with an accurate assessment of whether they have low or high risk of a successful outcome following IVF

Cone-beam CT reconstruction with gravity-induced motion.

Fixed-gantry cone-beam computed tomography (CBCT), where the imaging hardware is fixed while the subject is continuously rotated 360° in the horizontal position, has implications for building compact and affordable fixed-gantry linear accelerators (linacs). Fixed-gantry imaging with a rotating subject presents a challenging image reconstruction problem where the gravity-induced motion is coupled to the subject's rotation angle. This study is the first to investigate the feasibility of fixed-gantry CBCT using imaging data of three live rabbits in an ethics-approved study. A novel data-driven motion correction method that combines partial-view reconstruction and motion compensation was developed to overcome this challenge. Fixed-gantry CBCT scans of three live rabbits were acquired on a standard radiotherapy system with the imaging beam fixed and the rabbits continuously rotated using an in-house programmable rotation cradle. The reconstructed images of the thoracic region were validated against conventional CBCT scans acquired at different cradle rotation angles. Results showed that gravity-induced motion caused severe motion blur in all of the cases if unaccounted for. The proposed motion correction method yielded clinically usable image quality with  <1 mm gravity-induced motion blur for rabbits that were securely immobilized on the rotation cradle. Shapes of the anatomic structures were correctly reconstructed with  <0.5 mm accuracy. Translational motion accounted for the majority of gravity-induced motion. The motion-corrected reconstruction represented the time-averaged location of the thoracic region over a 360° rotation. The feasibility of fixed-gantry CBCT has been demonstrated. Future work involves the validation of imaging accuracy for human subjects, which will be useful for emerging compact fixed-gantry radiotherapy systems

Isolation of a wide range of minerals from a thermally treated plant: Equisetum arvense, a Mare’s tale

Silica is the second most abundant biomineral being exceeded in nature only by biogenic CaCO3. Many land plants (such as rice, cereals, cucumber, etc.) deposit silica in significant amounts to reinforce their tissues and as a systematic response to pathogen attack. One of the most ancient species of living vascular plants, Equisetum arvense is also able to take up and accumulate silica in all parts of the plant. Numerous methods have been developed for elimination of the organic material and/or metal ions present in plant material to isolate biogenic silica. However, depending on the chemical and/or physical treatment applied to branch or stem from Equisetum arvense; other mineral forms such glass-type materials (i.e. CaSiO3), salts (i.e. KCl) or luminescent materials can also be isolated from the plant material. In the current contribution, we show the chemical and/or thermal routes that lead to the formation of a number of different mineral types in addition to biogenic silica