Hydraulic Subsurface Measurements and Hydrodynamic Modeling as Indicators for Groundwater Flow Systems in the Rotondo Granite, Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man-made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high-relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.

Modeling of subsurface pulsed radar for nondestructive testing of structures

The use of pulsed radar for investigating the integrity of structural elements is gaining popularity and becoming firmly established as a nondestructive test method in civil engineering. Difficulties can often arise in the interpretation of results obtained, particularly where internal details are relatively complex. One approach that can be used to understand and evaluate radar results is through numerical modeling of signal propagation and reflection. By comparing the results of a numerical modeling with those from field measurements, engineers can gain valuable insight into the probable features embedded beneath the surface of a structural element. This paper discusses a series of numerical techniques for modeling subsurface radar and compares the precision of the results with those taken from real field data. It is found that more complex problems require more sophisticated analysis techniques to obtain realistic results, with a consequential increase in the computational resources to carry out the modeling.

A micromechanics approach for damage modeling of polymer matrix composites

A new model for damage evolution in polymer matrix composites is presented. The model is based on a combination of two constituent-level models and an interphase model. This approach reduces the number of empirical parameters since the two constituent- level models are formulated for isotropic materials, namely fiber and matrix. Decomposition of the state variables down to the micro-scale is accomplished by micromechanics. Phenomenological damage evolution models are then postulated for each constituent. Determination of material parameters is made from available experimental data. The required experimental data can be obtained with standard tests. Comparison between model predictions and additional experimental data is presented.

Selective Orthosteric Free Fatty Acid Receptor 2 (FFA2) Agonists IDENTIFICATION OF THE STRUCTURAL AND CHEMICAL REQUIREMENTS FOR SELECTIVE ACTIVATION OF FFA2 VERSUS FFA3

Free fatty acid receptor 2 (FFA2; GPR43) is a G protein-coupled seven-transmembrane receptor for short-chain fatty acids (SCFAs) that is implicated in inflammatory and metabolic disorders. The SCFA propionate has close to optimal ligand efficiency for FFA2 and can hence be considered as highly potent given its size. Propionate, however, does not discriminate between FFA2 and the closely related receptor FFA3 (GPR41). To identify FFA2-selective ligands and understand the molecular basis for FFA2 selectivity, a targeted library of small carboxylic acids was examined using holistic, label-free dynamic mass redistribution technology for primary screening and the receptor-proximal G protein [S-35] guanosine 5'-(3-O-thio) triphosphate activation, inositol phosphate, and cAMP accumulation assays for hit confirmation. Structure-activity relationship analysis allowed formulation of a general rule to predict selectivity for small carboxylic acids at the orthosteric binding site where ligands with substituted sp(3)-hybridized alpha-carbons preferentially activate FFA3, whereas ligands with sp(2)- or sp-hybridized alpha-carbons prefer FFA2. The orthosteric binding mode was verified by site-directed mutagenesis: replacement of orthosteric site arginine residues by alanine in FFA2 prevented ligand binding, and molecular modeling predicted the detailed mode of binding. Based on this, selective mutation of three residues to their non-conserved counterparts in FFA3 was sufficient to transfer FFA3 selectivity to FFA2. Thus, selective activation of FFA2 via the orthosteric site is achievable with rather small ligands, a finding with significant implications for the rational design of therapeutic compounds selectively targeting the SCFA receptors.

Bidirectional, iterative approach to the structural delineation of the functional "Chemoprint" in GPR40 for agonist recognition

GPR40, free fatty acid receptor 1 (FFAR1), is a member of the GPCR superfamily and a possible target for the treatment of type 2 diabetes. In this work, we conducted a bidirectional iterative investigation, including computational modeling and site-directed mutagenesis, aimed at delineating amino acid residues forming the functional "chemoprint" of GPR40 for agonist recognition. The computational and experimental studies revolved around the recognition of the potent synthetic agonist GW9508. Our experimentally supported model suggested that H137(4.56), R183(5.39), N244(6.55), and R258(7.35) are directly involved in interactions with the ligand. We have proposed a polarized NH-pi interaction between H137(4.56) and GW9508 as one of the contributing forces leading to the high potency of GW9508. The modeling approach presented in this work provides a general strategy for the exploration of receptor-ligand interactions in G-protein coupled receptors beginning prior to acquisition of experimental data.

Structural and functional relationships in the virulence-associated cathepsin L proteases of the parasitic liver fluke, Fasciola hepatica

The helminth parasite Fasciola hepatica secretes cysteine proteases to facilitate tissue invasion, migration, and development within the mammalian host. The major proteases cathepsin L1 (FheCL1) and cathepsin L2 (FheCL2) were recombinantly produced and biochemically characterized. By using site-directed mutagenesis, we show that residues at position 67 and 205, which lie within the S2 pocket of the active site, are critical in determining the substrate and inhibitor specificity. FheCL1 exhibits a broader specificity and a higher substrate turnover rate compared with FheCL2. However, FheCL2 can efficiently cleave substrates with a Pro in the P2 position and degrade collagen within the triple helices at physiological pH, an activity that among cysteine proteases has only been reported for human cathepsin K. The 1.4-A three-dimensional structure of the FheCL1 was determined by x-ray crystallography, and the three-dimensional structure of FheCL2 was constructed via homology-based modeling. Analysis and comparison of these structures and our biochemical data with those of human cathepsins L and K provided an interpretation of the substrate-recognition mechanisms of these major parasite proteases. Furthermore, our studies suggest that a configuration involving residue 67 and the "gatekeeper" residues 157 and 158 situated at the entrance of the active site pocket create a topology that endows FheCL2 with its unusual collagenolytic activity. The emergence of a specialized collagenolytic function in Fasciola likely contributes to the success of this tissue-invasive parasite.

Structural Breaks in Mortality Models and their Consequences

In recent years, the issue of life expectancy has become of upmost importance to pension providers, insurance companies and the government bodies in the developed world. Significant and consistent improvements in mortality rates and, hence, life expectancy have led to unprecedented increases in the cost of providing for older ages. This has resulted in an explosion of stochastic mortality models forecasting trends in mortality data in order to anticipate future life expectancy and, hence, quantify the costs of providing for future aging populations. Many stochastic models of mortality rates identify linear trends in mortality rates by time, age and cohort, and forecast these trends into the future using standard statistical methods. The modeling approaches used failed to capture the effects of any structural change in the trend and, thus, potentially produced incorrect forecasts of future mortality rates. In this paper, we look at a range of leading stochastic models of mortality and test for structural breaks in the trend time series.

Advances in GPCR modeling evaluated by the GPCR Dock 2013 assessment:Meeting new challenges

Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems.

Status of GPCR modeling and docking as reflected by community-wide GPCR Dock 2010 assessment

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe.

Modeling study of woody biomass:Interactions of cellulose, hemicellulose, and lignin

Lignocellulosic biomass pretreatment and the subsequent thermal conversion processes to produce solid, liquid, and gas biofuels are attractive solutions for today's energy challenges. The structural study of the main components in biomass and their macromolecular complexes is an active and ongoing research topic worldwide. The interactions among the three main components, cellulose, hemicellulose, and lignin, are studied in this paper using electronic structure methods, and the study includes examining the hydrogen bond network of cellulose-hemicellulose systems and the covalent bond linkages of hemicellulose-lignin systems. Several methods (semiempirical, Hartree-Fock, and density functional theory) using different basis sets were evaluated. It was shown that theoretical calculations can be used to simulate small model structures representing wood components. By comparing calculation results with experimental data, it was concluded that B3LYP/6-31G is the most suitable basis set to describe the hydrogen bond system and B3LYP/6-31G(d,p) is the most suitable basis set to describe the covalent system of woody biomass. The choice of unit model has a much larger effect on hydrogen bonding within cellulose-hemicellulose system, whereas the model choice has a minimal effect on the covalent linkage in the hemicellulose-lignin system. © 2011 American Chemical Society.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

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.

Elastodynamic Modeling and Analysis for an Exechon Parallel Kinematic Machine

As a newly invented parallel kinematic machine (PKM), Exechon has attracted intensive attention from both academic and industrial fields due to its conceptual high performance. Nevertheless, the dynamic behaviors of Exechon PKM have not been thoroughly investigated because of its structural and kinematic complexities. To identify the dynamic characteristics of Exechon PKM, an elastodynamic model is proposed with the substructure synthesis technique in this paper. The Exechon PKM is divided into a moving platform subsystem, a fixed base subsystem and three limb subsystems according to its structural features. Differential equations of motion for the limb subsystem are derived through finite element (FE) formulations by modeling the complex limb structure as a spatial beam with corresponding geometric cross sections. Meanwhile, revolute, universal, and spherical joints are simplified into virtual lumped springs associated with equivalent stiffnesses and mass at their geometric centers. Differential equations of motion for the moving platform are derived with Newton's second law after treating the platform as a rigid body due to its comparatively high rigidity. After introducing the deformation compatibility conditions between the platform and the limbs, governing differential equations of motion for Exechon PKM are derived. The solution to characteristic equations leads to natural frequencies and corresponding modal shapes of the PKM at any typical configuration. In order to predict the dynamic behaviors in a quick manner, an algorithm is proposed to numerically compute the distributions of natural frequencies throughout the workspace. Simulation results reveal that the lower natural frequencies are strongly position-dependent and distributed axial-symmetrically due to the structure symmetry of the limbs. At the last stage, a parametric analysis is carried out to identify the effects of structural, dimensional, and stiffness parameters on the system's dynamic characteristics with the purpose of providing useful information for optimal design and performance improvement of the Exechon PKM. The elastodynamic modeling methodology and dynamic analysis procedure can be well extended to other overconstrained PKMs with minor modifications.

Assessment Of various sensors for structural health monitoring for bridge weigh-in-motion (B-WIM).

In recent years, Structural Health Monitoring (SHM) systems have been developed to monitor bridge deterioration, assess real load levels and hence extend bridge life and safety. A road bridge is only safe if the stresses caused by the passing vehicles are less than the capacity of the bridge to resist them. Conventional SHM systems can be used to improve knowledge of the bridges capacity to resist stresses but generally give no information on the causes of any increase in stresses (based on measuring strain). The concept of in Bridge Weigh-in-Motion (B-WIM) is to establish axle loads, without interruption to traffic flow, by using strain sensors at a bridge soffit and subsequently converting the data to real time axle loads or stresses. Recent studies have shown it would be most beneficial to develop a portable system which can be easily attached to existing and new bridge structures for a specified monitoring period. The sensors could then be left in place while the data acquisition can be moved for various other sites. Therefore it is necessary to find accurate sensors capable of capturing peak strains under dynamic load and suitable methods for attaching these strain sensors to existing and new bridge structures. Additionally, it is important to ensure accurate strain transfer between concrete and steel, the adhesives layer and the strain sensor. This paper describes research investigating the suitably of using various sensors for the monitoring of concrete structures under dynamic vehicle load. Electrical resistance strain (ERS) gauges, vibrating wire (VW) gauges and fibre optic sensors (FOS) are commonly used for SHM. A comparative study will be carried out to select a suitable sensor for a bridge Weigh in Motion System. This study will look at fixing methods, durability, scanning rate and accuracy range. Finite element modeling is used to predict the strains which are then validated in laboratory trials.