Combining multi-scale geophysical techniques for robust hydro-structural characterisation in catchments underlain by hard rock in post-glacial regions

Accurate conceptual models of groundwater systems are essential for correct interpretation of monitoring data in catchment studies. In surface-water dominated hard rock regions, modern ground and surface water monitoring programmes often have very high resolution chemical, meteorological and hydrological observations but lack an equivalent emphasis on the subsurface environment, the properties of which exert a strong control on flow pathways and interactions with surface waters. The reasons for this disparity are the complexity of the system and the difficulty in accurately characterising the subsurface, except locally at outcrops or in boreholes. This is particularly the case in maritime north-western Europe, where a legacy of glacial activity, combined with large areas underlain by heterogeneous igneous and metamorphic bedrock, make the structure and weathering of bedrock difficult to map or model. Traditional approaches which seek to extrapolate information from borehole to field-scale are of limited application in these environments due to the high degree of spatial heterogeneity. Here we apply an integrative and multi-scale approach, optimising and combining standard geophysical techniques to generate a three-dimensional geological conceptual model of the subsurface in a catchment in NE Ireland. Available airborne LiDAR, electromagnetic and magnetic data sets were analysed for the region. At field-scale surface geophysical methods, including electrical resistivity tomography, seismic refraction, ground penetrating radar and magnetic surveys, were used and combined with field mapping of outcrops and borehole testing. The study demonstrates how combined interpretation of multiple methods at a range of scales produces robust three-dimensional conceptual models and a stronger basis for interpreting groundwater and surface water monitoring data.

Uranium halide complexes in ionic liquids: an electrochemical and structural study

The electrochemistry of the salts, [emim](2)[UBr6] and [emim](2)[UO2Br4] ([emim] = 1-ethyl-3-methylimidazolium), has been investigated in both a basic and an acidic bromoaluminate(III) ionic liquid. In the basic ionic liquid, the hexabromo salt undergoes a one-electron reversible reduction process at a stationary glassy carbon disc electrode, while the tetrabromodioxo salt was reduced to a uranium(IV) species by an irreversible two-electron process with the simultaneous transfer of oxide to the ionic liquid. On the other hand, dissolution of either of the salts in an acidic bromoaluminate( III) ionic liquid resulted in the formation of the same electroactive species. The solid state structures of the uranium chloride salts, [emim](2)[UCl6] and [emim](2)[UO2Cl4], have previously been reported, but have now been re-evaluated using a new statistical model developed in our group, to determine the presence or absence of weak hydrogen bonding interactions in the crystalline state.

Verification of a novel material model to predict damage in composite structures

An intralaminar damage model (IDM), based on continuum damage mechanics, was developed for the simulation of composite structures subjected to damaging loads. This model can capture the complex intralaminar damage mechanisms, accounting for mode interactions, and delaminations. Its development is driven by a requirement for reliable crush simulations to design composite structures with a high specific energy absorption. This IDM was implemented as a user subroutine within the commercial finite element package, Abaqus/Explicit[1]. In this paper, the validation of the IDM is presented using two test cases. Firstly, the IDM is benchmarked against published data for a blunt notched specimen under uniaxial tensile loading, comparing the failure strength as well as showing the damage. Secondly, the crush response of a set of tulip-triggered composite cylinders was obtained experimentally. The crush loading and the associated energy of the specimen is compared with the FE model prediction. These test cases show that the developed IDM is able to capture the structural response with satisfactory accuracy

Kinetostatic-model-based Stiffness Analysis of Exechon PKM

As a comparative newly-invented PKM with over-constraints in kinematic chains, the Exechon has attracted extensive attention from the research society. Different from the well-recognized kinematics analysis, the research on the stiffness characteristics of the Exechon still remains as a challenge due to the structural complexity. In order to achieve a thorough understanding of the stiffness characteristics of the Exechon PKM, this paper proposed an analytical kinetostatic model by using the substructure synthesis technique. The whole PKM system is decomposed into a moving platform subsystem, three limb subsystems and a fixed base subsystem, which are connected to each other sequentially through corresponding joints. Each limb body is modeled as a spatial beam with a uniform cross-section constrained by two sets of lumped springs. The equilibrium equation of each individual limb assemblage is derived through finite element formulation and combined with that of the moving platform derived with Newtonian method to construct the governing kinetostatic equations of the system after introducing the deformation compatibility conditions between the moving platform and the limbs. By extracting the 6 x 6 block matrix from the inversion of the governing compliance matrix, the stiffness of the moving platform is formulated. The computation for the stiffness of the Exechon PKM at a typical configuration as well as throughout the workspace is carried out in a quick manner with a piece-by-piece partition algorithm. The numerical simulations reveal a strong position-dependency of the PKM's stiffness in that it is symmetric relative to a work plane due to structural features. At the last stage, the effects of some design variables such as structural, dimensional and stiffness parameters on system rigidity are investigated with the purpose of providing useful information for the structural optimization and performance enhancement of the Exechon PKM. It is worthy mentioning that the proposed methodology of stiffness modeling in this paper can also be applied to other overconstrained PKMs and can evaluate the global rigidity over workplace efficiently with minor revisions.

The structural response of an oscillating wave surge converter to wave slamming

This paper describes the problems in experimentally obtaining hydrodynamic loads on an oscillating wave surge converter during slamming events, with the aim of furthering understanding of full scale hydrodynamic loads that flap type devices must be designed to withstand. Including how hydro-elastic effects and structural response are linked and why they are essential to the measurement of impulsive hydrodynamic loads. A combined experimental and numerical structural response study carried out on a 40th scale Oyster model drew conclusions on the structural vibration observed in the strain gauge load cell measurement. A further structural response study on a piezo electric load measurement device gave an insight into the advantages it could bring to reducing hydro-elastic effects.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Efficient Marginal Likelihood Computation for Gaussian Process Regression

In a Bayesian learning setting, the posterior distribution of a predictive model arises from a trade-off between its prior distribution and the conditional likelihood of observed data. Such distribution functions usually rely on additional hyperparameters which need to be tuned in order to achieve optimum predictive performance; this operation can be efficiently performed in an Empirical Bayes fashion by maximizing the posterior marginal likelihood of the observed data. Since the score function of this optimization problem is in general characterized by the presence of local optima, it is necessary to resort to global optimization strategies, which require a large number of function evaluations. Given that the evaluation is usually computationally intensive and badly scaled with respect to the dataset size, the maximum number of observations that can be treated simultaneously is quite limited. In this paper, we consider the case of hyperparameter tuning in Gaussian process regression. A straightforward implementation of the posterior log-likelihood for this model requires O(N^3) operations for every iteration of the optimization procedure, where N is the number of examples in the input dataset. We derive a novel set of identities that allow, after an initial overhead of O(N^3), the evaluation of the score function, as well as the Jacobian and Hessian matrices, in O(N) operations. We prove how the proposed identities, that follow from the eigendecomposition of the kernel matrix, yield a reduction of several orders of magnitude in the computation time for the hyperparameter optimization problem. Notably, the proposed solution provides computational advantages even with respect to state of the art approximations that rely on sparse kernel matrices.

Testing of a Market Fraction Model and Power-law Behaviour in the DAX 30

This paper tests a simple market fraction asset pricing model with heterogeneous
agents. By selecting a set of structural parameters of the model through a systematic procedure, we show that the autocorrelations (of returns, absolute returns and squared returns) of the market fraction model share the same pattern as those of the DAX 30. By conducting econometric analysis via Monte Carlo simulations, we characterize these power-law behaviours and find that estimates of the power-law decay indices, the (FI)GARCH parameters, and the tail index of the selected market fraction model closely match those of the DAX 30. The results strongly support the explanatory power of the heterogeneous agent models.

Modelling Static and Dynamic FRP-Concrete Bond Behaviour Using a Local Concrete Damage Model

This paper presents a study on the bond behaviour of FRP-concrete bonded joints under static and dynamic loadings, by developing a meso-scale finite element model using the K&C concrete damage model in LS-DYNA. A significant number of single shear experiments under static pull-off loading were modelled with an extensive parametric study covering key factors in the K&C model, including the crack band width, the compressive fracture energy and the shear dilatation factor. It is demonstrated that the developed model can satisfactorily simulate the static debonding behaviour, in terms of mesh objectivity, the load-carrying capacity and the local bond-slip behaviour, provided that proper consideration is given to the selection of crack band width and shear dilatation factor. A preliminary study of the effect of the dynamic loading rate on the debonding behaviour was also conducted by considering a dynamic increase factor (DIF) for the concrete strength as a function of strain rate. It is shown that a higher loading rate leads to a higher load-carrying capacity, a longer effective bond length, and a larger damaged area of concrete in the single shear loading scenario.

FE Simulation of Viscoplastic Consistency Model

Finite element (FE) simulations provides an inexpensive alternative for material testingof new metal alloys. Carrying out experimental testing is expensive. Nanoindentation is particularly costly due to the equipment needed to work on such a scale. FE simulations provide an inexpensive means of material testing if accurately carried out. This paper will demonstrate the applicability and accuracy of using FE modelling for basic material tests and will propose that the viscoplastic model may be used for nanoindentation testing. The simulations will test the Young’s modulus of materials during analysis when an Abaqus VUMAT is used. The viscoplastic model is incorporated into a subroutine and is tested at the macroscopic scale against previous published results.

Oxygen vacancy formation in CeO2 and Ce1-xZrxO2 solid solutions:electron localization, electrostatic potential and structural relaxation

Ceria (CeO2) and ceria-based composite materials, especially Ce1-xZrxO2 solid solutions, possess a wide range of applications in many important catalytic processes, such as three-way catalysts, owing to their excellent oxygen storage capacity (OSC) through the oxygen vacancy formation and refilling. Much of this activity has focused on the understanding of the electronic and structural properties of defective CeO2 with and without doping, and comprehending the determining factor for oxygen vacancy formation and the rule to tune the formation energy by doping has constituted a central issue in material chemistry related to ceria. However, the calculation on electronic structures and the corresponding relaxation patterns in defective CeO2-x oxides remains at present a challenge in the DFT framework. A pragmatic approach based on density functional theory with the inclusion of on-site Coulomb correction, i.e. the so-called DFT + U technique, has been extensively applied in the majority of recent theoretical investigations. Firstly, we review briefly the latest electronic structure calculations of defective CeO2(111), focusing on the phenomenon of multiple configurations of the localized 4f electrons, as well as the discussions of its formation mechanism and the catalytic role in activating the O-2 molecule. Secondly, aiming at shedding light on the doping effect on tuning the oxygen vacancy formation in ceria-based solid solutions, we summarize the recent theoretical results of Ce1-xZrxO2 solid solutions in terms of the effect of dopant concentrations and crystal phases. A general model on O vacancy formation is also discussed; it consists of electrostatic and structural relaxation terms, and the vital role of the later is emphasized. Particularly, we discuss the crucial role of the localized structural relaxation patterns in determining the superb oxygen storage capacity in kappa-phase Ce1-xZr1-xO2. Thirdly, we briefly discuss some interesting findings for the oxygen vacancy formation in pure ceria nanoparticles (NPs) uncovered by DFT calculations and compare those with the bulk or extended surfaces of ceria as well as different particle sizes, emphasizing the role of the electrostatic field in determining the O vacancy formation.

A LEED STRUCTURAL STUDY OF THE CO INDUCED RECONSTRUCTION OF PD(110) - EVIDENCE FOR A MISSING ROW STRUCTURE

Molecularly adsorbed CO on Pd{110} has been shown (R. Raval et al., Chem. Phys. Lett. 167 (1990) 391, ref. [1]) to induce a substantial reconstruction of the surface in the coverage range 0.3 <theta less-than-or-equal-to 0.75. Throughout this coverage range, the adsorbate-covered reconstructed surface exhibits a (4 x 2) LEED pattern. However, the exact nature of the reconstruction remains uncertain. We have conducted a LEED I(E) "fingerprinting" analysis of the CO/Pd{110}-(4 x 2) structure in order to establish the type of reconstruction induced in the metal surface. This study shows that the LEED I(E) profiles of the integral order and appropriate half-order beams of the CO/Pd{110}-(4 x 2) pattern closely resemble the I(E) profiles theoretically calculated for a Pd{110}-(1 x 2) missing-row structure. Additionally, there is a strong resemblance to the experimental LEED I(E) profiles for the Cs/Pd{110}-(1 x 2) structure which has also been shown to exhibit the missing-row structure. On the basis of this evidence we conclude that the CO/Pd{110}-(4 x 2) LEED pattern arises from a missing-row reconstruction of the Pd{110} surface which gives rise to a strong underlying (1 x 2) pattern plus a poorly ordered CO overlayer which produces weak, diffuse fourth-order spots in the LEED pattern.

Concrete damage plasticity model for modeling FRP-to-concrete bond behavior

The technique of externally bonding fiber-reinforced polymer (FRP) composites has become very popular worldwide for retrofitting existing reinforced concrete (RC) structures. Debonding of FRP from the concrete substrate is a typical failure mode in such strengthened structures. The bond behavior between FRP and concrete thus plays a crucial role in these structures. The FRP-to-concrete bond behavior has been extensively investigated experimentally, commonly using a single or double shear test of the FRP-to-concrete bonded joint. Comparatively, much less research has been concerned with numerical simulation, chiefly due to difficulties in the accurate modeling of the complex behavior of concrete. This paper presents a simple but robust finite-element (FE) model for simulating the bond behavior in the entire debonding process for the single shear test. A concrete damage plasticity model is proposed to capture the concrete-to-FRP bond behavior. Numerical results are in close agreement with test data, validating the model. In addition to accuracy, the model has two further advantages: it only requires the basic material parameters (i.e., no arbitrary user-defined parameter such as the shear retention factor is required) and it can be directly implemented in the FE software ABAQUS.