Optimization of chemical plant simulation using double collocation

A method has been constructed for the solution of a wide range of chemical plant simulation models including differential equations and optimization. Double orthogonal collocation on finite elements is applied to convert the model into an NLP problem that is solved either by the VF 13AD package based on successive quadratic programming, or by the GRG2 package, based on the generalized reduced gradient method. This approach is termed simultaneous optimization and solution strategy. The objective functional can contain integral terms. The state and control variables can have time delays. Equalities and inequalities containing state and control variables can be included into the model as well as algebraic equations and inequalities. The maximum number of independent variables is 2. Problems containing 3 independent variables can be transformed into problems having 2 independent variables using finite differencing. The maximum number of NLP variables and constraints is 1500. The method is also suitable for solving ordinary and partial differential equations. The state functions are approximated by a linear combination of Lagrange interpolation polynomials. The control function can either be approximated by a linear combination of Lagrange interpolation polynomials or by a piecewise constant function over finite elements. The number of internal collocation points can vary by finite elements. The residual error is evaluated at arbitrarily chosen equidistant grid-points, thus enabling the user to check the accuracy of the solution between collocation points, where the solution is exact. The solution functions can be tabulated. There is an option to use control vector parameterization to solve optimization problems containing initial value ordinary differential equations. When there are many differential equations or the upper integration limit should be selected optimally then this approach should be used. The portability of the package has been addressed converting the package from V AX FORTRAN 77 into IBM PC FORTRAN 77 and into SUN SPARC 2000 FORTRAN 77. Computer runs have shown that the method can reproduce optimization problems published in the literature. The GRG2 and the VF I 3AD packages, integrated into the optimization package, proved to be robust and reliable. The package contains an executive module, a module performing control vector parameterization and 2 nonlinear problem solver modules, GRG2 and VF I 3AD. There is a stand-alone module that converts the differential-algebraic optimization problem into a nonlinear programming problem.

Resonance studies of artificial earth satellites

Orbit determination from artificial satellite observations is a key process in obtaining information about the Earth and its environment. A study of the perturbations experienced by these satellites enables knowledge to be gained of the upper atmosphere, the gravity field, ocean tides, solid-Earth tides and solar radiation. The gravity field is expressed as a double infinite series of associated Legendre functions (tesseral harmonics). In contemporary global gravity field models the overall geoid is well determined. An independent check on these gravity field harmonics of a particular order may be made by analysis of satellites that pass through resonance of that order. For such satellites the perturbations of the orbital elements close to resonance are analysed to derive lumped harmonic coefficients. The orbital parameters of 1984-106A have been determined at 43 epochs, during which time the satellite was close to 14th order resonance. Analysis of the inclination and eccentricity yielded 6 lumped harmonic coefficients of order 14 whilst analysis of the mean motion yielded additional pairs of lumped harmonics of orders 14, 28 and 42, with the 14th order harmonics superseding those obtained from analysis of the inclination. This thesis concentrates in detail on the theoretical changes of a near-circular satellite orbit perturbed by the Earth's gravity field under the influence of minimal air-drag whilst in resonance with the Earth. The satellite 1984-106A experienced the interesting property of being temporarily trapped with respect to a secondary resonance parameter due to the low air-drag in 1987. This prompted the theoretical investigation of such a phenomenon. Expressions obtained for the resonance parameter led to the determination of 8 lumped harmonic coefficients, coincidental to those already obtained. All the derived lumped harmonic values arc used to test the accuracy of contemporary gravity field models and the underlying theory in this thesis.

The DHR1 domain of DOCK180 binds to SNX5 and regulates cation-independent mannose 6-phosphate receptor transport

DOCK180 is the archetype of the DOCK180-family guanine nucleotide exchange factor for small GTPases Rac1 and Cdc42. DOCK180-family proteins share two conserved domains, called DOCK homology region (DHR)-1 and -2. Although the function of DHR2 is to activate Rac1, DHR1 is required for binding to phosphoinositides. To better understand the function of DHR1, we searched for its binding partners by direct nanoflow liquid chromatography/tandem mass spectrometry, and we identified sorting nexins (SNX) 1, 2, 5, and 6, which make up a multimeric protein complex mediating endosome-to-trans-Golgi-network (TGN) retrograde transport of the cation-independent mannose 6-phosphate receptor (CI-MPR). Among these SNX proteins, SNX5 was coimmunoprecipitated with DOCK180 most efficiently. In agreement with this observation, DOCK180 colocalized with SNX5 at endosomes. The RNA interference-mediated knockdowns of SNX5 and DOCK180, but not Rac1, resulted in the redistribution of CI-MPR from TGN to endosomes. Furthermore, expression of the DOCK180 DHR1 domain was sufficient to restore the perturbed CI-MPR distribution in DOCK180 knockdown cells. These data suggest that DOCK180 regulates CI-MPR trafficking via SNX5 and that this function is independent of its guanine nucleotide exchange factor activity toward Rac1.

System identification of linear vibrating structures

Methods of dynamic modelling and analysis of structures, for example the finite element method, are well developed. However, it is generally agreed that accurate modelling of complex structures is difficult and for critical applications it is necessary to validate or update the theoretical models using data measured from actual structures. The techniques of identifying the parameters of linear dynamic models using Vibration test data have attracted considerable interest recently. However, no method has received a general acceptance due to a number of difficulties. These difficulties are mainly due to (i) Incomplete number of Vibration modes that can be excited and measured, (ii) Incomplete number of coordinates that can be measured, (iii) Inaccuracy in the experimental data (iv) Inaccuracy in the model structure. This thesis reports on a new approach to update the parameters of a finite element model as well as a lumped parameter model with a diagonal mass matrix. The structure and its theoretical model are equally perturbed by adding mass or stiffness and the incomplete number of eigen-data is measured. The parameters are then identified by an iterative updating of the initial estimates, by sensitivity analysis, using eigenvalues or both eigenvalues and eigenvectors of the structure before and after perturbation. It is shown that with a suitable choice of the perturbing coordinates exact parameters can be identified if the data and the model structure are exact. The theoretical basis of the technique is presented. To cope with measurement errors and possible inaccuracies in the model structure, a well known Bayesian approach is used to minimize the least squares difference between the updated and the initial parameters. The eigen-data of the structure with added mass or stiffness is also determined using the frequency response data of the unmodified structure by a structural modification technique. Thus, mass or stiffness do not have to be added physically. The mass-stiffness addition technique is demonstrated by simulation examples and Laboratory experiments on beams and an H-frame.

Molecular regulation of iron uptake in pseudomonas aeruginosa

The microbial demand for iron is often met by the elaboration of siderophores into the surrounding medium and expression of cognate outer membrane receptors for the ferric siderophore complexes. Conditions of iron limitation, such as those encountered in vivo, cause Pseudomonas aeruginosa to express two high-affinity iron-uptake systems based on pyoverdin and pyochelin. These systems will operate both in the organism's natural habitat, soil and water, where the solubility of iron at neutral pH is extremely low, and in the human host where the availability of free iron is too low to sustain bacterial growth due to the iron-binding glycoproteins transferrin and lactoferrin. Cross-feeding and radiolabelled iron uptake experiments demonstrated that pyoverdin biosynthesis and uptake were highly heterogeneous amongst P.aeruginosa strains, that growth either in the presence of pyoverdin or pyochelin resulted in induction of specific IROMPs, and that induction of iron uptake is siderophore-specific. The P.aeruginosa Tn5 mutant PH1 is deficient in ferripyoverdin uptake and resistant to pyocin Sa, suggesting that the site of interaction of pyocin Sa is a ferripyoverdin receptor. Additional Tn5 mutants appeared to exploit different strategies to achieve pyocin Sa-resistance, involving modifications in expression of pyoverdin-mediated iron uptake, indicating that complex regulatory systems exist to enable these organisms to compete effectively for iron. Modulation of expression of IROMPs prompted a study of the mechanism of uptake of a semi-synthetic C(7) α-formamido substituted cephalosporin BRL 41897A. Sensitivity to this agent correlated with expression of the 75 kDa ferri-pyochelin receptor and demonstrated the potential of high-affinity iron uptake systems for targeting of novel antibiotics. Studies with ferri-pyoverdin uptake-deficient mutant PH1 indicated that expression of outer membrane protein G (OprG), which is usually expressed under iron-rich conditions and repressed under iron-deficient conditions, was perturbed. Attempts were made to clone the oprG gene using a degenerate probe based on the N-terminal amino acid sequence. A strongly hybridising HindIll restriction fragment was cloned and sequenced, but failed to reveal an open reading frame correspondmg to OprG. However, there appears to be good evidence that a part of the gene codmg for the hydrophilic membrane-associated ATP-binding component of a hitherto uncharacterised periplasmic- binding-protein-dependent transport system has been isolated. The full organisation and sequence of the operon, and substrate for this putative transport system, are yet: to be elucidated,

Influence of iron on bacterial infections in leukaemia

The influence of iron metabolism, both on the invading bacterial pathogen and in the host is widespread and often appears to be crucial in determining the outcome of an infection. This study involved the investigation of leukaemia, a clinical disease where abnormal availability of iron may play a part in predisposing patients to bacterial infection. The iron status throughout a Gram-negative septicaemia and in 20 random, newly diagnosed leukaemic patients was assessed. The results revealed that the majority of the patients exhibited high serum iron levels and serum transferrin saturation often at 100%, with an inability to reduce the latter to within normal values during an infection episode. The antibody response to P.aeruginosa, E.coli and K.pneumoniae outer membrane protein (OMP) antigens were investigated by immunoblotting with sequential serum samples during infection in the leukaemic host. Antibodies to all the major OMPs, were observed, although recognition of iron-regulated membrane proteins (IRMPs) was in many cases weak. Results from the enzyme-linked immunosorbent assay indicated that in all patients antibody titre in response to infection was poor. Sub-MICs of mitomycin C significantly altered the surface characteristics of P.aeruginosa. The silver-stained SDS-PAGE gels of proteinase K digested whole cell lysates of strains PAO1, 6750, M7 and PAJ indicated that core LPS was affected in the presence of mitomycin C. In contrast, the rough strain AK1012 showed no observable differences. Results obtained using quantitative gas-liquid chromatographic analysis showed the amount of LPS fatty acids to be unaffected, however, the KDO and carbohydrate content in strains PAO1, 6750 and M7 under Fe+ and Fe- growth conditions were decreased by up to 4-fold in the presence of mitomycin C, indicating perturbed expression of LPS. The cell surface became significantly more hydrophobic in the P.aeruginosa strains, except AK1012 which was comparatively unaffected. The induction of protein G (OprG) in P.aeruginosa was found to be a sensitive indicator of media iron. The data indicated that expression of OprG can be modulated by growth rate/phase, availability of iron and by the presence of ciprofloxacin in the growth medium.

Application of SOA-NOLM in all-optical processing

This thesis presents experimental investigations of the use of semiconductor optical amplifiers in a nonlinear loop mirror (SOA-NOLM) and its application in all-optical processing. The techniques used are mainly experimental and are divided into three major applications. Initially the semiconductor optical amplifier, SOA, is experimentally characterised and the optimum operating condition is identified. An interferometric switch based on a Sagnac loop with the SOA as the nonlinear element is employed to realise all-optical switching. All-optical switching is a very attractive alternative to optoelectronic conversion because it avoids the conversion from the optical to the electronic domain and back again. The first major investigation involves a carrier suppressed return to zero, CSRZ, format conversion and transmission. This study is divided into single channel and four channel WDM respectively. The optical bandwidth which limits the conversion is investigated. The improvement of the nonlinear tolerance in the CSRZ transmission is shown which shows the suitability of this format for enhancing system performance. Second, a symmetrical switching window is studied in the SOA-NOLM where two similar control pulses are injected into the SOA from opposite directions. The switching window is symmetric when these two control pulses have the same power and arrive at the same time in the SOA. Finally, I study an all-optical circulating shift register with an inverter. The detailed behaviour of the blocks of zeros and ones has been analysed in terms of their transient measurement. Good agreement with a simple model of the shift register is obtained. The transient can be reduced but it will affect the extinction ratio of the pulses.

A perturbative analysis of dispersion-managed solitons

We study soliton solutions of the path-averaged propagation equation governing the transmission of dispersion-managed (DM) optical pulses in the (practical) limit when residual dispersion and nonlinearity only slightly affect the pulse dynamics over one compensation period. In the case of small dispersion map strengths, the averaged pulse dynamics is governed by a perturbed form of the nonlinear Schrödinger equation; applying a perturbation theory – elsewhere developed – based on inverse scattering theory, we derive an analytic expression for the envelope of the DM soliton. This expression correctly predicts the power enhancement arising from the dispersion management. Theoretical results are verified by direct numerical simulations.

Phenotypic effect of mutations in evolving populations of RNA molecules

Background. The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. Results. We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). Conclusions. The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate µ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of µ, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks. © 2010 Stich et al; licensee BioMed Central Ltd.

A preliminary investigation into the effects of ocular lubricants on higher order aberrations in normal and dry eye subjects

Purpose: To study the effects of ocular lubricants on higher order aberrations in normal and self-diagnosed dry eyes. Methods: Unpreserved hypromellose drops, Tears Again™ liposome spray and a combination of both were administered to the right eye of 24 normal and 24 dry eye subjects following classification according to a 5 point questionnaire. Total ocular higher order aberrations, coma, spherical aberration and Strehl ratios for higher order aberrations were measured using the Nidek OPD-Scan III (Nidek Technologies, Gamagori, Japan) at baseline, immediately after application and after 60. min. The aberration data were analyzed over a 5. mm natural pupil using Zernike polynomials. Each intervention was assessed on a separate day and comfort levels were recorded before and after application. Corneal staining was assessed and product preference recorded after the final measurement for each intervention. Results: Hypromellose drops caused an increase in total higher order aberrations (p= <0.01 in normal and dry eyes) and a reduction in Strehl ratio (normal eyes: p= <0.01, dry eyes p= 0.01) immediately after instillation. There were no significant differences between normal and self-diagnosed dry eyes for response to intervention and no improvement in visual quality or reduction in higher order aberrations after 60. min. Differences in comfort levels failed to reach statistical significance. Conclusion: Combining treatments does not offer any benefit over individual treatments in self-diagnosed dry eyes and no individual intervention reached statistical significance. Symptomatic subjects with dry eye and no corneal staining reported an improvement in comfort after using lubricants. © 2013 British Contact Lens Association.

A perturbative analysis of dispersion-managed solitons

We study soliton solutions of the path-averaged propagation equation governing the transmission of dispersion-managed (DM) optical pulses in the (practical) limit when residual dispersion and nonlinearity only slightly affect the pulse dynamics over one compensation period. In the case of small dispersion map strengths, the averaged pulse dynamics is governed by a perturbed form of the nonlinear Schrödinger equation; applying a perturbation theory – elsewhere developed – based on inverse scattering theory, we derive an analytic expression for the envelope of the DM soliton. This expression correctly predicts the power enhancement arising from the dispersion management. Theoretical results are verified by direct numerical simulations.

Asymmetric reactions of organosilicon compounds and properties of novel optically active organosilanes

Partial reduction of racemic methoxysilanes by 1:1 complexes of lithium aluminium hydride with optically active cinchona and ephedra alkaloids give optically active silanes and methoxysilanes. Optical yields depend on the groups attached to silicon and the alkaloid used but in some cases approach 50%, The method has been used to prepare novel optically active organosilanes, possessing an asymmetric silicon centre, which are either inaccessible by any of the other available routes or would require a time consuming preparation. Such compounds are of use in the study of the mechanism of substitutions at silicon. Attempts have been made to rationalize the results of the asymmetric reductions in terms of differences in sterio and electronic interactions in diastereoisomeric transition states. Circular dichroism and optical rotatory dispersion spectra have been obtained for the optically active products in an attempt to elucidate the absolute configurations of the novel asymmetric organosilanes. The results from these studies provide a useful addition to the data so far accumulated for asymmetrically perturbed aromatic chromophores. Nuclear magnetic resonanoe studies of diastereoisomaric (-)-menthoxysilanes show that these compounds possess resonances extremely useful in the determination of optical purities for asymmetric organosilanes which possess an aromatic group. The effect of variable temperature on the spectra has revealed evidence for the conformational preferences in these compounds. Other diastereoisomeric alkoxysilanes have been prepared and their n.m.r.spectra studied in the hope of establishing trends. Exploratory studies for other asymmetric reactions proceeding at silicon have proved unfruitful.

Determination of a spacewise dependent heat source

This paper investigates the inverse problem of determining a spacewise dependent heat source in the parabolic heat equation using the usual conditions of the direct problem and information from a supplementary temperature measurement at a given single instant of time. The spacewise dependent temperature measurement ensures that the inverse problem has a unique solution, but this solution is unstable, hence the problem is ill-posed. For this inverse problem, we propose an iterative algorithm based on a sequence of well-posed direct problems which are solved at each iteration step using the boundary element method (BEM). The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for various typical benchmark test examples which have the input measured data perturbed by increasing amounts of random noise.

Measuring the spatial extent of texture pooling using reverse correlation

The local image representation produced by early stages of visual analysis is uninformative regarding spatially extensive textures and surfaces. We know little about the cortical algorithm used to combine local information over space, and still less about the area over which it can operate. But such operations are vital to support perception of real-world objects and scenes. Here, we deploy a novel reverse-correlation technique to measure the extent of spatial pooling for target regions of different areas placed either in the central visual field, or more peripherally. Stimuli were large arrays of micropatterns, with their contrasts perturbed individually on an interval-by-interval basis. By comparing trial-by-trial observer responses with the predictions of computational models, we show that substantial regions (up to 13 carrier cycles) of a stimulus can be monitored in parallel by summing contrast over area. This summing strategy is very different from the more widely assumed signal selection strategy (a MAX operation), and suggests that neural mechanisms representing extensive visual textures can be recruited by attention. We also demonstrate that template resolution is much less precise in the parafovea than in the fovea, consistent with recent accounts of crowding. © 2014 The Authors.

Communication between the nucleotide binding domains of P-glycoprotein occurs via conformational changes that involve residue 508

Our aim is to provide molecular understanding of the mechanisms underlying the (i) interaction between the two nucleotide binding domains (NBDs) and (ii) coupling between NBDs and transmembrane domains within P-glycoprotein (Pgp) during a transport cycle. To facilitate this, we have introduced a number of unique cysteine residues at surface exposed positions (E393C, S452C, I500C, N508C, and K578C) in the N-terminal NBD of Pgp, which had previously been engineered to remove endogenous cysteines. Positions of the mutations were designed using a model based on crystallographic features of prokaryotic NBDs. The single cysteine mutants were expressed in insect cells using recombinant baculovirus and the proteins purified by metal affinity chromatography by virtue of a polyhistidine tag. None of the introduced cysteine residues perturbed the function of Pgp as judged by the characteristics of drug stimulated ATP hydrolysis. The role of residues at each of the introduced sites in the catalytic cycle of Pgp was investigated by the effect of covalent conjugation with N-ethyl-maleimide (NEM). All but one mutation (K578C) was accessible to labeling with [3H]-NEM. However, perturbation of ATPase activity was only observed for the derivitized N508C isoform. The principle functional manifestation was a marked inhibition of the "basal" rate of ATP hydrolysis. Neither the extent nor potency to which a range of drugs could affect the ATPase activity were altered in the NEM conjugated N508C isoform. The results imply that the accessibility of residue 508, located in the alpha-helical subdomain of NBD1 in Pgp, is altered by the conformational changes that occur during ATP hydrolysis.