Development and application of a novel algorithm for steady-state optimising control using dynamic information

A novel optimising controller is designed that leads a slow process from a sub-optimal operational condition to the steady-state optimum in a continuous way based on dynamic information. Using standard results from optimisation theory and discrete optimal control, the solution of a steady-state optimisation problem is achieved by solving a receding-horizon optimal control problem which uses derivative and state information from the plant via a shadow model and a state-space identifier. The paper analyzes the steady-state optimality of the procedure, develops algorithms with and without control rate constraints and applies the procedure to a high fidelity simulation study of a distillation column optimisation.

Ensemble prediction for nowcasting with a convection-permitting model - II: forecast error statistics

A 24-member ensemble of 1-h high-resolution forecasts over the Southern United Kingdom is used to study short-range forecast error statistics. The initial conditions are found from perturbations from an ensemble transform Kalman filter. Forecasts from this system are assumed to lie within the bounds of forecast error of an operational forecast system. Although noisy, this system is capable of producing physically reasonable statistics which are analysed and compared to statistics implied from a variational assimilation system. The variances for temperature errors for instance show structures that reflect convective activity. Some variables, notably potential temperature and specific humidity perturbations, have autocorrelation functions that deviate from 3-D isotropy at the convective-scale (horizontal scales less than 10 km). Other variables, notably the velocity potential for horizontal divergence perturbations, maintain 3-D isotropy at all scales. Geostrophic and hydrostatic balances are studied by examining correlations between terms in the divergence and vertical momentum equations respectively. Both balances are found to decay as the horizontal scale decreases. It is estimated that geostrophic balance becomes less important at scales smaller than 75 km, and hydrostatic balance becomes less important at scales smaller than 35 km, although more work is required to validate these findings. The implications of these results for high-resolution data assimilation are discussed.

A method for merging flow-dependent forecast error statistics from an ensemble with static statistics for use in high resolution variational data assimilation

The background error covariance matrix, B, is often used in variational data assimilation for numerical weather prediction as a static and hence poor approximation to the fully dynamic forecast error covariance matrix, Pf. In this paper the concept of an Ensemble Reduced Rank Kalman Filter (EnRRKF) is outlined. In the EnRRKF the forecast error statistics in a subspace defined by an ensemble of states forecast by the dynamic model are found. These statistics are merged in a formal way with the static statistics, which apply in the remainder of the space. The combined statistics may then be used in a variational data assimilation setting. It is hoped that the nonlinear error growth of small-scale weather systems will be accurately captured by the EnRRKF, to produce accurate analyses and ultimately improved forecasts of extreme events.

High frequency variability of the Atlantic meridional overturning circulation

We compare the variability of the Atlantic meridional overturning circulation (AMOC) as simulated by the coupled climate models of the RAPID project, which cover a wide range of resolution and complexity, and observed by the RAPID/MOCHA array at about 26N. We analyse variability on a range of timescales, from five-daily to interannual. In models of all resolutions there is substantial variability on timescales of a few days; in most AOGCMs the amplitude of the variability is of somewhat larger magnitude than that observed by the RAPID array, while the time-mean is within about 10% of the observational estimate. The amplitude of the simulated annual cycle is similar to observations, but the shape of the annual cycle shows a spread among the models. A dynamical decomposition shows that in the models, as in observations, the AMOC is predominantly geostrophic (driven by pressure and sea-level gradients), with both geostrophic and Ekman contributions to variability, the latter being exaggerated and the former underrepresented in models. Other ageostrophic terms, neglected in the observational estimate, are small but not negligible. The time-mean of the western boundary current near the latitude of the RAPID/MOCHA array has a much wider model spread than the AMOC does, indicating large differences among models in the simulation of the wind-driven gyre circulation, and its variability is unrealistically small in the models. In many RAPID models and in models of the Coupled Model Intercomparison Project Phase 3 (CMIP3), interannual variability of the maximum of the AMOC wherever it lies, which is a commonly used model index, is similar to interannual variability in the AMOC at 26N. Annual volume and heat transport timeseries at the same latitude are well-correlated within 15--45N, indicating the climatic importance of the AMOC. In the RAPID and CMIP3 models, we show that the AMOC is correlated over considerable distances in latitude, but not the whole extent of the north Atlantic; consequently interannual variability of the AMOC at 50N, where it is particularly relevant to European climate, is not well-correlated with that of the AMOC at 26N, where it is monitored by the RAPID/MOCHA array.

The synthesis and properties of [Fe3(CO)11(RCN)] derivatives, and a high-yield general route to mono- and di-substituted derivatives of [Fe3(CO)12]; the crystal and molecular structure of [Fe3(CO)11(NCC6H4-2-Me)]

The clusters [Fe3(CO)11(RCN)] (1: R = Me, C3H5, C6H5, or C6H4-2-Me) have been prepared at low temperature from [Fe3(CO)12] and RCN in the presence of Me3NO. Compounds 1 react essentially quantitatively with a wide range of two-electron donors, L, (viz.: CO, PPh3, P(OMe)3, PPh2H, PPh2Me, PF3, CyNC (Cy = cyclohexyl), P(OEt)3, SbPh3, PBu3, AsPh3, or SnR2 (R = CH(SiMe3)2)) to give [Fe3(CO)11L] (2). In some cases (2), on treatment with Me3NO and then L′ (L′ = a second two-electron donor) yields [Fe3(CO)10LL′] in high yield. The crystal and molecular structures of 1 (L = NCC6H4Me-2) have been determined by a full single crystal structure analysis, and shown to have an axial nitrile coordinated at the unique iron atom, with two CO groups bridging the other two metal atoms.

Generic unified modelling process for developing semantically rich, dynamic and temporal models

Models play a vital role in supporting a range of activities in numerous domains. We rely on models to support the design, visualisation, analysis and representation of parts of the world around us, and as such significant research effort has been invested into numerous areas of modelling; including support for model semantics, dynamic states and behaviour, temporal data storage and visualisation. Whilst these efforts have increased our capabilities and allowed us to create increasingly powerful software-based models, the process of developing models, supporting tools and /or data structures remains difficult, expensive and error-prone. In this paper we define from literature the key factors in assessing a model’s quality and usefulness: semantic richness, support for dynamic states and object behaviour, temporal data storage and visualisation. We also identify a number of shortcomings in both existing modelling standards and model development processes and propose a unified generic process to guide users through the development of semantically rich, dynamic and temporal models.

Thermal and high hydrostatic pressure inactivation of myrosinase from green cabbage: a kinetic study

Myrosinase, a family of enzymes which coexist with glucosinolates in all Brassica vegetables, catalyses the hydrolysis of glucosinolates to yield compounds that can have beneficial effects on human health. In this study, the thermal and pressure inactivation of myrosinase from green cabbage was kinetically investigated. Thermal inactivation started at 35 C and inactivation kinetics was studied in the temperature range 35–55 C. Thermal inactivation of green cabbage myrosinase followed the well known consecutive step model. Pressure inactivation started at 300 MPa, even at 10 C, and the consecutive step model effectively described pressure inactivation in the range 300–450 MPa at 10 C. The combined effects of applying various pressures and temperatures on myrosinase inactivation kinetics were studied in the ranges 35–50 C and, 100–400 MPa. The inactivation followed first-order kinetics at all of the applied combinations. This study demonstrates that myrosinase from green cabbage is highly susceptible to both thermal and high pressure processing. Furthermore, it is also noted that myrosinase stability during processing appears to vary widely between different Brassica species.

Electro-optical properties of liquid crystal copolymers and their relationship to structural order

A range of side chain liquid crystal copolymers have been prepared using mesogenic and non-mesogenic units. It is found that high levels of the non-mesogenic moieties may be introduced without completely disrupting the organization of the liquid crystal phase. Incorporation of this comonomer causes a marked reduction in the glass transition temperature (Tg), presumably as a result of enhanced backbone mobility and a corresponding lowering of the nematic transition temperature, thereby restricting the temperature range for stability of the liquid crystal phase. The effect of the interactions between the various components of these side-chain polymers on their electro-optic responses is described. Infrared (i.r.) dichroism measurements have been made to determine the order parameters of the liquid crystalline side-chain polymers. By identifying a certain band (CN stretching) in the i.r. absorption spectrum, the order parameter of the mesogenic groups can be obtained. The temperature and composition dependence of the observed order parameter are related to the liquid crystal phase transitions and to the electro-optic response. It is found that the introduction of the non-mesogenic units into the polymer chain lowers the threshold voltage of the electro-optic response over and above that due to the reduction in the order parameter. The dynamic electro-optic responses are dominated by the temperature-dependent viscosity and evidence is presented for relaxation processes involving the polymer backbone which are on a time scale greater than that for the mesogenic side-chain units.

Chaos in high performance digital robot controllers

Active robot force control requires some form of dynamic inner loop control for stability. The author considers the implementation of position-based inner loop control on an industrial robot fitted with encoders only. It is shown that high gain velocity feedback for such a robot, which is effectively stationary when in contact with a stiff environment, involves problems beyond the usual caveats on the effects of unknown environment stiffness. It is shown that it is possible for the controlled joint to become chaotic at very low velocities if encoder edge timing data are used for velocity measurement. The results obtained indicate that there is a lower limit on controlled velocity when encoders are the only means of joint measurement. This lower limit to speed is determined by the desired amount of loop gain, which is itself determined by the severity of the nonlinearities present in the drive system.

Effect of high-hydrostatic pressure and pH on the rheological properties of gum arabic

The effect on the viscoelastic behaviour, of pressure-treating hydrated gumarabic samples (800 MPa) at different pH values (2.8, 4.2 and 8.0) was investigated, using controlled stress rheometry. The treated samples were analysed for their complex (G∗), storage (G′) and loss (G″) moduli as a function of frequency, using dynamic oscillatory testing. Significant changes in the rheologicalproperties were observed in both the pressurised gum solutions and in those previously buffered at pH 2.8. The gum, at its natural pH (4.25) and at alkaline pH (8.0), was enhanced by pressure treatment, but only for the already “good” quality gum samples. High-pressure treatment had substantial effects on the frequency-dependence of the moduli of both the pressurised and the pressurised/pH-treated solutions, with the latter being more pronounced, suggesting differing structures or changes in the overall degree of interaction of the gum systems after pressure treatment.

High prevalence of multidrug-tolerant bacteria and associated antimicrobial resistance genes isolated from ornamental fish and their carriage water

Background: Antimicrobials are used to directly control bacterial infections in pet (ornamental) fish and are routinely added to the water these fish are shipped in to suppress the growth of potential pathogens during transport. Methodology/Principal Findings: To assess the potential effects of this sustained selection pressure, 127 Aeromonas spp. isolated from warm and cold water ornamental fish species were screened for tolerance to 34 antimicrobials. Representative isolates were also examined for the presence of 54 resistance genes by a combination of miniaturized microarray and conventional PCR. Forty-seven of 94 Aeromonas spp. isolates recovered from tropical ornamental fish and their carriage water were tolerant to >= 15 antibiotics, representing seven or more different classes of antimicrobial. The quinolone and fluoroquinolone resistance gene, qnrS2, was detected at high frequency (37% tested recent isolates were positive by PCR). Class 1 integrons, IncA/C broad host range plasmids and a range of other antibiotic resistance genes, including floR, blaTEM21, tet(A), tet(D), tet(E), qacE2, sul1, and a number of different dihydrofolate reductase and aminoglycoside transferase coding genes were also detected in carriage water samples and bacterial isolates. Conclusions: These data suggest that ornamental fish and their carriage water act as a reservoir for both multi-resistant bacteria and resistance genes.

A high prevalence of antimicrobial resistant Escherichia coli isolated from pigs and a low prevalence of antimicrobial resistant E-coli from cattle and sheep in Great Britain at slaughter

The incidence of antimicrobial resistance and expressed and unexpressed resistance genes among commensal Escherichia coli isolated from healthy farm animals at slaughter in Great Britain was investigated. The prevalence of antimicrobial resistance among the isolates varied according to the animal species; of 836 isolates from cattle tested only 5.7% were resistant to one or more antimicrobials, while only 3.0% of 836 isolates from sheep were resistant to one or more agents. However, 92.1% of 2480 isolates from pigs were resistant to at least one antimicrobial. Among isolates from pigs, resistance to some antimicrobials such as tetracycline (78.7%), sulphonamide (66.9%) and streptomycin (37.5%) was found to be common, but relatively rare to other agents such as amikacin (0.1%), ceftazidime ( 0.1%) and coamoxiclav (0.2%). The isolates had a diverse range of resistance gene profiles, with tet(B), sul2 and strAB identified most frequently. Seven out of 615 isolates investigated carried unexpressed resistance genes. One trimethoprim-susceptible isolate carried a complete dfrA17 gene but lacked a promoter for it. However, in the remaining six streptomycin-susceptible isolates, one of which carried strAB while the others carried aadA, no mutations or deletions in gene or promoter sequences were identified to account for susceptibility. The data indicate that antimicrobial resistance in E. coli of animal origin is due to a broad range of acquired genes.

Properties of singular vectors using convective available potential energy as final time norm

We study the feasibility of using the singular vector technique to create initial condition perturbations for short-range ensemble prediction systems (SREPS) focussing on predictability of severe local storms and in particular deep convection. For this a new final time semi-norm based on the convective available potential energy (CAPE) is introduced. We compare singular vectors using the CAPE-norm with SVs using the more common total energy (TE) norm for a 2-week summer period in 2007, which includes a case of mesoscale extreme rainfall in the south west of Finland. The CAPE singular vectors perturb the CAPE field by increasing the specific humidity and temperature of the parcel and increase the lapse rate above the parcel in the lower troposphere consistent with physical considerations. The CAPE-SVs are situated in the lower troposphere. This in contrast to TE-SVs with short optimization times which predominantly remain in the high troposphere. By examining the time evolution of the CAPE singular values we observe that the convective event in the south west of Finland is clearly associated with high CAPE singular values.

Simulated anthropogenic changes in the Brewer–Dobson circulation, including its extension to high latitudes

Recent studies using comprehensive middle atmosphere models predict a strengthening of the Brewer-Dobson circulation in response to climate change. To gain confidence in the realism of this result it is important to quantify and understand the contributions from the different components of stratospheric wave drag that cause this increase. Such an analysis is performed here using three 150-yr transient simulations from the Canadian Middle Atmosphere Model (CMAM), a Chemistry-Climate Model that simulates climate change and ozone depletion and recovery. Resolved wave drag and parameterized orographic gravity wave drag account for 60% and 40%, respectively, of the long-term trend in annual mean net upward mass flux at 70 hPa, with planetary waves accounting for 60% of the resolved wave drag trend. Synoptic wave drag has the strongest impact in northern winter, where it accounts for nearly as much of the upward mass flux trend as planetary wave drag. Owing to differences in the latitudinal structure of the wave drag changes, the relative contribution of resolved and parameterized wave drag to the tropical upward mass flux trend over any particular latitude range is highly sensitive to the range of latitudes considered. An examination of the spatial structure of the climate change response reveals no straightforward connection between the low-latitude and high-latitude changes: while the model results show an increase in Arctic downwelling in winter, they also show a decrease in Antarctic downwelling in spring. Both changes are attributed to changes in the flux of stationary planetary wave activity into the stratosphere.