Stochastic protein folding simulation in the three-dimensional HP-model

We present results from three-dimensional protein folding simulations in the HP-model on ten benchmark problems. The simulations are executed by a simulated annealing-based algorithm with a time-dependent cooling schedule. The neighbourhood relation is determined by the pull-move set. The results provide experimental evidence that the maximum depth D of local minima of the underlying energy landscape can be upper bounded by D < n(2/3). The local search procedure employs the stopping criterion (In/delta)(D/gamma) where m is an estimation of the average number of neighbouring conformations, gamma relates to the mean of non-zero differences of the objective function for neighbouring conformations, and 1-delta is the confidence that a minimum conformation has been found. The bound complies with the results obtained for the ten benchmark problems. (c) 2008 Elsevier Ltd. All rights reserved.

Increasing entanglement through engineered disorder in the random Ising chain

The ground-state entanglement entropy between block of sites in the random Ising chain is studied by means of the Von Neumann entropy. We show that in presence of strong correlations between the disordered couplings and local magnetic fields the entanglement increases and becomes larger than in the ordered case. The different behavior with respect to the uncorrelated disordered model is due to the drastic change of the ground state properties. The same result holds also for the random three-state quantum Potts model.

Plant-wide predictive control for a thermal power plant based on a physical plant model

A constrained non-linear, physical model-based, predictive control (NPMPC) strategy is developed for improved plant-wide control of a thermal power plant. The strategy makes use of successive linearisation and recursive state estimation using extended Kalman filtering to obtain a linear state-space model. The linear model and a quadratic programming routine are used to design a constrained long-range predictive controller One special feature is the careful selection of a specific set of plant model parameters for online estimation, to account for time-varying system characteristics resulting from major system disturbances and ageing. These parameters act as nonstationary stochastic states and help to provide sufficient degrees-of-freedom to obtain unbiased estimates of controlled outputs. A 14th order non-linear plant model, simulating the dominant characteristics of a 200 MW oil-fired pou er plant has been used to test the NPMPC algorithm. The control strategy gives impressive simulation results, during large system disturbances and extremely high rate of load changes, right across the operating range. These results compare favourably to those obtained with the state-space GPC method designed under similar conditions.

The Effects of Anti-Hypertensive Drugs Evaluated Using Markov Modelling for Northern Ireland Chronic Kidney Disease Patients

The aim of this paper is to use Markov modelling to
investigate survival for particular types of kidney patients
in relation to their exposure to anti-hypertensive treatment
drugs. In order to monitor kidney function an intuitive three
point assessment is proposed through the collection of blood
samples in relation to Chronic Kidney Disease for Northern
Ireland patients. A five state Markov Model was devised
using specific transition probabilities for males and
females over all age groups. These transition probabilities
were then adjusted appropriately using relative risk scores
for the event death for different subgroups of patients. The
model was built using TreeAge software package in order to
explore the effects of anti-hypertensive drugs on patients.

The 2p interloper resonances in collisions

Ab initio electron scattering calculations using the R -matrix approach have been performed for within a three-state valence configuration-interaction model (VCI). The lowest three electronic target states ( , and the ) of this molecular nitrogen cation are included in the close-coupling method, with each state being represented by a valence CI approximation. From a detailed analysis of the resonance structure found in our work for the symmetries we find four prominent Rydberg series of the type , , , and a interloper resonance. This interloper molecular resonance associated with the B state of is seen to cause distortions of the resulting resonance spectra. A comparison of our total cross sections for the X - B transition shows excellent agreement with the available experimental data.

Dynamic Multivariate Statistical Process Control using Subspace Identification

This paper points out a serious flaw in dynamic multivariate statistical process control (MSPC). The principal component analysis of a linear time series model that is employed to capture auto- and cross-correlation in recorded data may produce a considerable number of variables to be analysed. To give a dynamic representation of the data (based on variable correlation) and circumvent the production of a large time-series structure, a linear state space model is used here instead. The paper demonstrates that incorporating a state space model, the number of variables to be analysed dynamically can be considerably reduced, compared to conventional dynamic MSPC techniques.

Kinetic Monte Carlo study of electrochemical growth in a heteroepitaxial system

Structural and kinetic aspects of 2-D irreversible metal deposition under potentiostatic conditions are analyzed by means of dynamic Monte Carlo simulations employing embedded atom potentials for a model system. Three limiting models, all considering adatom diffusion, were employed to describe adatom deposition. The first model (A) considers adatom deposition on any free substrate site on the surface at the same rate. The second model (B) considers adatom deposition only on substrate sites which exhibit no neighboring sites occupied by adatoms. The third model (C) allows deposition at higher rates on sites presenting neighboring sites occupied by adatoms. Under the proper conditions, the coverage (theta) versus time (t) relationship for the three cases can be heuristically fitted to the functional form theta = 1 - exp(-betat(alpha)), where alpha and beta are parameters. We suggest that the value of the parameter alpha can be employed to distinguish experimentally between the three cases. While model A trivially delivers a = 1, models B and C are characterized by alpha 1, respectively.

Electron-scale electrostatic solitary waves and shocks: the role of superthermal electrons

The propagation of electron-acoustic solitary waves and shock structures is investigated in a plasma characterized by a superthermal electron population. A three-component plasma model configuration is employed, consisting of inertial (“cold”) electrons, inertialess ? (kappa) distributed superthermal (“hot”) electrons and stationary ions. A multiscale method is employed, leading to a Korteweg-de Vries (KdV) equation for the electrostatic potential (in the absence of dissipation). Taking into account dissipation, a hybrid Korteweg-de Vries-Burgers (KdVB) equation is derived. Exact negative-potential pulse- and kink-shaped solutions (shocks) are obtained. The relative strength among dispersion, nonlinearity and damping coefficients is discussed. Excitations formed in superthermal plasma (finite ?) are narrower and steeper, compared to the Maxwellian case (infinite ?). A series of numerical simulations confirms that energy initially stored in a solitary pulse which propagates in a stable manner for large ? (Maxwellian plasma) may break down to smaller structures or/and to random oscillations, when it encounters a small-? (nonthermal) region. On the other hand, shock structures used as initial conditions for numerical simulations were shown to be robust, essentially responding to changed in the environment by a simple profile change (in width).

Towards a better understanding of the nature and contribution of the interflow to catchment dynamics.

A detailed understanding of flow and contaminant transfer along each of the key hydrological pathways within a catchment is critical for designing and implementing cost effective Programmes of Measures under the Water
Framework Directive.
The Contaminant Movement along Pathways Project (’The Pathways Project’) is an Irish, EPA STRIVE funded, large multi-disciplinary project which is focussed on understanding and modelling flow and attenuation along each of these pathways for the purposes of developing a catchment management tool. The tool will be used by EPA and RBD catchment managers to assess and manage the impacts of diffuse contamination on stream aquatic ecology. Four main contaminants of interest — nitrogen, phosphorus, sediment and pathogens — are being
investigated in four instrumented test catchments. In addition to the usual hydrological and water chemistry/quality parameters typically captured in catchment studies, field measurements at the test catchments include ecological
sampling, sediment dynamics, soil moisture dynamics, and groundwater levels and chemistry/quality, both during and between significant rainfall events. Spatial and temporal sampling of waters directly from the pathways of
interest is also being carried out.
Sixty-five percent of Ireland is underlain by poorly productive aquifers. In these hydrogeological settings, the main pathways delivering flow to streams are overland flow, interflow and shallow bedrock flow. Little is
known about the interflow pathway and its relative importance in delivery of flow and contaminants to the streams. Interflow can occur in both the topsoil and subsoil, and may include unsaturated matrix flow, bypass or macropore
flow, saturated flow in locally perched water tables and artificial field drainage.
Results to date from the test catchment experiments show that artificial field drains play an important role in the delivery of interflow to these streams, during and between rainfall events when antecedent conditions are
favourable. Hydrochemical mixing models, using silica and SAC254 (the absorbance of UV light at a wavelength of 254 nm which is a proxy for dissolved organic matter) as tracers, show that drain flow is an important end
member contributing to the stream and that proportionally, its contribution is relatively high.
Results from the study also demonstrate that waters originating from one pathway often mix with the waters from another, and are subsequently delivered to the stream at rates, and with chemical/quality characteristics,
that are not typical of either pathway. For example, pre-event shallow groundwater not far from the catchment divide comes up to the surface as rejected recharge during rainfall events and is rapidly delivered to the stream
via overland flow and/or artificial land drainage, bringing with it higher nitrate than would often be expected from a quickflow pathway contribution. This is contrary to the assumption often made in catchment studies that the
deeper hydrological pathways have slower response times in stream hydrographs during a rainfall event, and it emphasizes that it is critical to have a strong three-dimensional conceptual model as the basis for the interpretation
of catchment data.

Numerically Investigating the Effects of Cross-Links in Scaled Microchannel Heat Sinks

Thermal management as a method of heightening performance in miniaturized electronic devices using microchannel heat sinks has recently become of interest to researchers and the industry. One of the current challenges is to design heat sinks with uniform flow distribution. A number of experimental studies have been conducted to seek appropriate designs for microchannel heat sinks. However, pursuing this goal experimentally can be an expensive endeavor. The present work investigates the effect of cross-links on adiabatic two-phase flow in an array of parallel channels. It is carried out using the three dimensional mixture model from the computational fluid dynamics software, FLUENT 6.3. A straight channel and two cross-linked channel models were simulated. The cross-links were located at 1/3 and 2/3 of the channel length, and their widths were one and two times larger than the channel width. All test models had 45 parallel rectangular channels, with a hydraulic diameter of 1.59 mm. The results showed that the trend of flow distribution agrees with experimental results. A new design, with cross-links incorporated, was proposed and the results showed a significant improvement of up to 55% on flow distribution compared with the standard straight channel configuration without a penalty in the pressure drop. Further discussion about the effect of cross-links on flow distribution, flow structure, and pressure drop was also documented.

Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Numerical prediction of the low-velocity impact damage and compression after impact strength of composite laminates

Low-velocity impact damage can drastically reduce the residual mechanical properties of the composite structure even when there is barely visible impact damage. The ability to computationally predict the extent of damage and compression after impact (CAI) strength of a composite structure can potentially lead to the exploration of a larger design space without incurring significant development time and cost penalties. A three-dimensional damage model, to predict both low-velocity impact damage and compression after impact CAI strength of composite laminates, has been developed and implemented as a user material subroutine in the commercial finite element package, ABAQUS/Explicit. The virtual tests were executed in two steps, one to capture the impact damage and the other to predict the CAI strength. The observed intra-laminar damage features, delamination damage area as well as residual strength are discussed. It is shown that the predicted results for impact damage and CAI strength correlated well with experimental testing.

Modelling of Organic Rankine Cycle for Waste Heat Recovery Process in Supercritical Condition

Organic Rankine Cycle (ORC) is the most commonly used method for recovering energy from small sources of heat. The investigation of the ORC in supercritical condition is a new research area as it has a potential to generate high power and thermal efficiency in a waste heat recovery system. This paper presents a steady state ORC model in supercritical condition and its simulations with a real engine’s exhaust data. The key component of ORC, evaporator, is modelled using finite volume method, modelling of all other components of the waste heat recovery system such as pump, expander and condenser are also presented. The aim of this paper is to investigate the effects of mass flow rate and evaporator outlet temperature on the efficiency of the waste heat recovery process. Additionally, the necessity of maintaining an optimum evaporator outlet temperature is also investigated. Simulation results show that modification of mass flow rate is the key to changing the operating temperature at the evaporator outlet.

Predicting the crushing behaviour of composite material using high-fidelity finite element modelling

The capability to numerically model the crushing behaviour of composite structures will enable the efficient design of structures with high specific energy absorption capacity. This is particularly relevant to the aerospace and automotive industries where cabin structures need to be shown to be crashworthy. In this paper, a three-dimensional damage model is presented, which accurately represents the behaviour of composite laminates under crush loading. Both intralaminar and interlaminar failure mechanisms are taken into account. The crush damage model was implemented in ABAQUS/Explicit as a VUMAT subroutine. Numerical predictions are shown to agree well with experimental results, accurately capturing the intralaminar and interlaminar damage for a range of stacking sequences, triggers and composite materials. The use of measured material parameters required by the numerical models, without the need to ‘calibrate’ this input data, demonstrates this computational tool's predictive capabilities