Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Modelling of induction motors for system faults and transient stability studies

The research carried out in this thesis was mainly concerned with the effects of large induction motors and their transient performance in power systems. Computer packages using the three phase co-ordinate frame of reference were developed to simulate the induction motor transient performance. A technique using matrix algebra was developed to allow extension of the three phase co-ordinate method to analyse asymmetrical and symmetrical faults on both sides of the three phase delta-star transformer which is usually required when connecting large induction motors to the supply system. System simulation, applying these two techniques, was used to study the transient stability of a power system. The response of a typical system, loaded with a group of large induction motors, two three-phase delta-star transformers, a synchronous generator and an infinite system was analysed. The computer software developed to study this system has the advantage that different types of fault at different locations can be studied by simple changes in input data. The research also involved investigating the possibility of using different integrating routines such as Runge-Kutta-Gill, RungeKutta-Fehlberg and the Predictor-Corrector methods. The investigation enables the reduction of computation time, which is necessary when solving the induction motor equations expressed in terms of the three phase variables. The outcome of this investigation was utilised in analysing an introductory model (containing only minimal control action) of an isolated system having a significant induction motor load compared to the size of the generator energising the system.

Thin film bearings with phase change of lubricant

The possible evaporation of lubricant in fluid film bearings has been investigated theoretically and by experiment using a radial flow hydrostatic bearing supplied with liquid refrigerant R114. Good correlation between measured and theoretical values was obtained using a bespoke computational fluid dynamic model in which the flow was assumed to be laminar and adiabatic. The effects of viscous dissipation and vapour generation within the fluid film are fully accounted for by applying a fourth order Runge-Kutta routine to satisfy the radial and filmwise transverse constraints of momentum, energy and mass conservation. The results indicate that the radial velocity profile remains parabolic while the flow remains in the liquid phase and that the radial rate of enthalpy generation is then constant across the film at a given radius. The results also show that evaporation will commence at a radial location determined by geometry and flow conditions and in fluid layers adjacent to the solid boundaries. Evaporation is shown to progress in the radial direction and the load carrying capacity of such a bearing is reduced significantly. Expressions for the viscosity of the liquid/vapour mixture found in the literature survey have not been tested against experimental data. A new formulation is proposed in which the suitable choice of a characteristic constant yields close representation to any of these expressions. Operating constraints imposed by the design of the experimental apparatus limited the extent of the surface over which evaporation could be obtained, and prevented clear identification of the most suitable relationship for the viscosity of the liquid/vapour mixture. The theoretical model was extended to examine the development of two phase flow in a rotating shaft face seal of uniform thickness. Previous theoretical analyses have been based on the assumption that the radial velocity profile of the flow is always parabolic, and that the tangential component of velocity varies linearly from the value at the rotating surface, to zero at the stationary surface. The computational fluid dynamic analysis shows that viscous shear and dissipation in the fluid adjacent to the rotating surface leads to developing evaporation with a consequent reduction in tangential shear forces. The tangential velocity profile is predicted to decay rapidly through the film, exhibiting a profile entirely different to that assumed by previous investigators. Progressive evaporation takes place close to the moving wall and does not occur completely at a single radial location, as has been claimed in earlier work.

Optimal Control of a Second Order Parabolic Heat Equation

In this paper, we are concerned with the optimal control boundary control of a second order parabolic heat equation. Using the results in [Evtushenko, 1997] and spatial central finite difference with diagonally implicit Runge-Kutta method (DIRK) is applied to solve the parabolic heat equation. The conjugate gradient method (CGM) is applied to solve the distributed control problem. Numerical results are reported.

Degradação fotocatalítica oxidativa do fenol utilizando carvão obtido da pirólise de diferentes biomassas

The modern industrial progress has been contaminating water with phenolic compounds. These are toxic and carcinogenic substances and it is essential to reduce its concentration in water to a tolerable one, determined by CONAMA, in order to protect the living organisms. In this context, this work focuses on the treatment and characterization of catalysts derived from the bio-coal, by-product of biomass pyrolysis (avelós and wood dust) as well as its evaluation in the phenol photocatalytic degradation reaction. Assays were carried out in a slurry bed reactor, which enables instantaneous measurements of temperature, pH and dissolved oxygen. The experiments were performed in the following operating conditions: temperature of 50 °C, oxygen flow equals to 410 mL min-1 , volume of reagent solution equals to 3.2 L, 400 W UV lamp, at 1 atm pressure, with a 2 hours run. The parameters evaluated were the pH (3.0, 6.9 and 10.7), initial concentration of commercial phenol (250, 500 and 1000 ppm), catalyst concentration (0, 1, 2, and 3 g L-1 ), nature of the catalyst (activated avelós carbon washed with dichloromethane, CAADCM, and CMADCM, activated dust wood carbon washed with dichloromethane). The results of XRF, XRD and BET confirmed the presence of iron and potassium in satisfactory amounts to the CAADCM catalyst and on a reduced amount to CMADCM catalyst, and also the surface area increase of the materials after a chemical and physical activation. The phenol degradation curves indicate that pH has a significant effect on the phenol conversion, showing better results for lowers pH. The optimum concentration of catalyst is observed equals to 1 g L-1 , and the increase of the initial phenol concentration exerts a negative influence in the reaction execution. It was also observed positive effect of the presence of iron and potassium in the catalyst structure: betters conversions were observed for tests conducted with the catalyst CAADCM compared to CMADCM catalyst under the same conditions. The higher conversion was achieved for the test carried out at acid pH (3.0) with an initial concentration of phenol at 250 ppm catalyst in the presence of CAADCM at 1 g L-1 . The liquid samples taken every 15 minutes were analyzed by liquid chromatography identifying and quantifying hydroquinone, p-benzoquinone, catechol and maleic acid. Finally, a reaction mechanism is proposed, cogitating the phenol is transformed into the homogeneous phase and the others react on the catalyst surface. Applying the model of Langmuir-Hinshelwood along with a mass balance it was obtained a system of differential equations that were solved using the Runge-Kutta 4th order method associated with a optimization routine called SWARM (particle swarm) aiming to minimize the least square objective function for obtaining the kinetic and adsorption parameters. Related to the kinetic rate constant, it was obtained a magnitude of 10-3 for the phenol degradation, 10-4 to 10-2 for forming the acids, 10-6 to 10-9 for the mineralization of quinones (hydroquinone, p-benzoquinone and catechol), 10-3 to 10-2 for the mineralization of acids.

Análise da biorremediação de compostos monoaromáticos em água através da pseudomonas aeruginosa

The monoaromatic compounds are toxic substances present in petroleum derivades and used broadly in the chemical and petrochemical industries. Those compounds are continuously released into the environment, contaminating the soil and water sources, leading to the possible unfeasibility of those hydrous resources due to their highly carcinogenic and mutagenic potentiality, since even in low concentrations, the BTEX may cause serious health issues. Therefore, it is extremely important to develop and search for new methodologies that assist and enable the treatment of BTEX-contaminated matrix. The bioremediation consists on the utilization of microbial groups capable of degrading hydrocarbons, promoting mineralization, or in other words, the permanent destruction of residues, eliminating the risks of future contaminations. This work investigated the biodegradation kinetics of water-soluble monoaromatic compounds (benzene, toluene and ethylbenzene), based on the evaluation of its consummation by the Pseudomonas aeruginosa bacteria, for concentrations varying from 40 to 200 mg/L. To do so, the performances of Monod kinetic model for microbial growth were evaluated and the material balance equations for a batch operation were discretized and numerically solved by the fourth order Runge-Kutta method. The kinetic parameters obtained using the method of least squares as statistical criteria were coherent when compared to those obtained from the literature. They also showed that, the microorganism has greater affinity for ethylbenzene. That way, it was possible to observe that Monod model can predict the experimental data for the individual biodegradation of the BTEX substrates and it can be applied to the optimization of the biodegradation processes of toxic compounds for different types of bioreactors and for different operational conditions.

Análise da biorremediação de compostos monoaromáticos em água através da pseudomonas aeruginosa

The monoaromatic compounds are toxic substances present in petroleum derivades and used broadly in the chemical and petrochemical industries. Those compounds are continuously released into the environment, contaminating the soil and water sources, leading to the possible unfeasibility of those hydrous resources due to their highly carcinogenic and mutagenic potentiality, since even in low concentrations, the BTEX may cause serious health issues. Therefore, it is extremely important to develop and search for new methodologies that assist and enable the treatment of BTEX-contaminated matrix. The bioremediation consists on the utilization of microbial groups capable of degrading hydrocarbons, promoting mineralization, or in other words, the permanent destruction of residues, eliminating the risks of future contaminations. This work investigated the biodegradation kinetics of water-soluble monoaromatic compounds (benzene, toluene and ethylbenzene), based on the evaluation of its consummation by the Pseudomonas aeruginosa bacteria, for concentrations varying from 40 to 200 mg/L. To do so, the performances of Monod kinetic model for microbial growth were evaluated and the material balance equations for a batch operation were discretized and numerically solved by the fourth order Runge-Kutta method. The kinetic parameters obtained using the method of least squares as statistical criteria were coherent when compared to those obtained from the literature. They also showed that, the microorganism has greater affinity for ethylbenzene. That way, it was possible to observe that Monod model can predict the experimental data for the individual biodegradation of the BTEX substrates and it can be applied to the optimization of the biodegradation processes of toxic compounds for different types of bioreactors and for different operational conditions.

Techniques en appui des formats de modulation avancés pour les futurs réseaux optiques

Les systèmes de communication optique avec des formats de modulation avancés sont actuellement l’un des sujets de recherche les plus importants dans le domaine de communication optique. Cette recherche est stimulée par les exigences pour des débits de transmission de donnée plus élevés. Dans cette thèse, on examinera les techniques efficaces pour la modulation avancée avec une détection cohérente, et multiplexage par répartition en fréquence orthogonale (OFDM) et multiples tonalités discrètes (DMT) pour la détection directe et la détection cohérente afin d’améliorer la performance de réseaux optiques. Dans la première partie, nous examinons la rétropropagation avec filtre numérique (DFBP) comme une simple technique d’atténuation de nonlinéarité d’amplificateur optique semiconducteur (SOA) dans le système de détection cohérente. Pour la première fois, nous démontrons expérimentalement l’efficacité de DFBP pour compenser les nonlinéarités générées par SOA dans un système de détection cohérente porteur unique 16-QAM. Nous comparons la performance de DFBP avec la méthode de Runge-Kutta quatrième ordre. Nous examinons la sensibilité de performance de DFBP par rapport à ses paramètres. Par la suite, nous proposons une nouvelle méthode d’estimation de paramètre pour DFBP. Finalement, nous démontrons la transmission de signaux de 16-QAM aux taux de 22 Gbaud sur 80km de fibre optique avec la technique d’estimation de paramètre proposée pour DFBP. Dans la deuxième partie, nous nous concentrons sur les techniques afin d’améliorer la performance des systèmes OFDM optiques en examinent OFDM optiques cohérente (CO-OFDM) ainsi que OFDM optiques détection directe (DDO-OFDM). Premièrement, nous proposons une combinaison de coupure et prédistorsion pour compenser les distorsions nonlinéaires d’émetteur de CO-OFDM. Nous utilisons une interpolation linéaire par morceaux (PLI) pour charactériser la nonlinéarité d’émetteur. Dans l’émetteur nous utilisons l’inverse de l’estimation de PLI pour compenser les nonlinéarités induites à l’émetteur de CO-OFDM. Deuxièmement, nous concevons des constellations irrégulières optimisées pour les systèmes DDO-OFDM courte distance en considérant deux modèles de bruit de canal. Nous démontrons expérimentalement 100Gb/s+ OFDM/DMT avec la détection directe en utilisant les constellations QAM optimisées. Dans la troisième partie, nous proposons une architecture réseaux optiques passifs (PON) avec DDO-OFDM pour la liaison descendante et CO-OFDM pour la liaison montante. Nous examinons deux scénarios pour l’allocations de fréquence et le format de modulation des signaux. Nous identifions la détérioration limitante principale du PON bidirectionnelle et offrons des solutions pour minimiser ses effets.

Aplicação da metodologia numérica “fast bounds crack” para uma estimativa eficiente da evolução do tamanho de trinca

A significant part of the life of a mechanical component occurs, the crack propagation stage in fatigue. Currently, it is had several mathematical models to describe the crack growth behavior. These models are classified into two categories in terms of stress range amplitude: constant and variable. In general, these propagation models are formulated as an initial value problem, and from this, the evolution curve of the crack is obtained by applying a numerical method. This dissertation presented the application of the methodology "Fast Bounds Crack" for the establishment of upper and lower bounds functions for model evolution of crack size. The performance of this methodology was evaluated by the relative deviation and computational times, in relation to approximate numerical solutions obtained by the Runge-Kutta method of 4th explicit order (RK4). Has been reached a maximum relative deviation of 5.92% and the computational time was, for examples solved, 130,000 times more higher than achieved by the method RK4. Was performed yet an Engineering application in order to obtain an approximate numerical solution, from the arithmetic mean of the upper and lower bounds obtained in the methodology applied in this work, when you don’t know the law of evolution. The maximum relative error found in this application was 2.08% which proves the efficiency of the methodology "Fast Bounds Crack".

Método Gráfico para establecer el campo de pendientes de una Ecuación Diferencial

Trabajo orientado en el área de ecuaciones diferenciales enfocándose en el método gráfico para establecer el campo de pendiente de una ecuación diferencial y el método de aproximaciones numéricas para aproximar la solución de una ecuación diferencial. Presenta los métodos de Euler, Runge-Kutta de cuarto orden y el método multipasos de Adams-Bashforth-Moulton. Asimismo, se explica las ecuaciones mediante el uso del software para los métodos gráficos tales como el Maple y Geogebra.

Fast formulas for computing cryptographic pairings

The most powerful known primitive in public-key cryptography is undoubtedly elliptic curve pairings. Upon their introduction just over ten years ago the computation of pairings was far too slow for them to be considered a practical option. This resulted in a vast amount of research from many mathematicians and computer scientists around the globe aiming to improve this computation speed. From the use of modern results in algebraic and arithmetic geometry to the application of foundational number theory that dates back to the days of Gauss and Euler, cryptographic pairings have since experienced a great deal of improvement. As a result, what was an extremely expensive computation that took several minutes is now a high-speed operation that takes less than a millisecond. This thesis presents a range of optimisations to the state-of-the-art in cryptographic pairing computation. Both through extending prior techniques, and introducing several novel ideas of our own, our work has contributed to recordbreaking pairing implementations.

Formulas for algorithmic computation of impedance matrix using general circuit constants

A set of formulas is derived from general circuit constants which facilitates formation of the impedance matrix of a power system by the bus-impedance method. The errors associated with the lumpedparameter representation of a transmission line are thereby eliminated. The formulas are valid for short lines also, if the relevant general circuit constants are employed. The mutual impedance between the added line and the existing system is not considered, but the approach suggested can well be extended to it.

Approximate synthesis formulas for microstrip line with aperture in ground plane

Approximate closed-form expressions for the propagation characteristics of a microstrip line with a symmetrical aperture in its ground plane are reported in this article. Well-known expressions for the characteristic impedance of a regular microstrip line have been modified to incorporate the effect of this aperture. The accuracy of these expressions for various values of substrate thickness, permittivity and line width has been studied in detail by fullwave simulations. This has been further verified by measurements. These expressions are easier to compute and find immense use in the design of broadband filters, tight couplers, power dividers, transformers, delay lines, and matching circuits. A broadband filter with aperture in ground plane is demonstrated in this article. (c) 2011 Wiley Periodicals, Inc. Int J RF and Microwave CAE, 2012.

An Exponential Lower Bound for Homogeneous Depth Four Arithmetic Formulas

We show here a 2(Omega(root d.log N)) size lower bound for homogeneous depth four arithmetic formulas. That is, we give an explicit family of polynomials of degree d on N variables (with N = d(3) in our case) with 0, 1-coefficients such that for any representation of a polynomial f in this family of the form f = Sigma(i) Pi(j) Q(ij), where the Q(ij)'s are homogeneous polynomials (recall that a polynomial is said to be homogeneous if all its monomials have the same degree), it must hold that Sigma(i,j) (Number of monomials of Q(ij)) >= 2(Omega(root d.log N)). The above mentioned family, which we refer to as the Nisan-Wigderson design-based family of polynomials, is in the complexity class VNP. Our work builds on the recent lower bound results 1], 2], 3], 4], 5] and yields an improved quantitative bound as compared to the quasi-polynomial lower bound of 6] and the N-Omega(log log (N)) lower bound in the independent work of 7].