Mathematical analysis of activation thresholds in enzyme-catalyzed positive feedbacks: application to the feedbacks of blood coagulation.

A hierarchy of enzyme-catalyzed positive feedback loops is examined by mathematical and numerical analysis. Four systems are described, from the simplest, in which an enzyme catalyzes its own formation from an inactive precursor, to the most complex, in which two sequential feedback loops act in a cascade. In the latter we also examine the function of a long-range feedback, in which the final enzyme produced in the second loop activates the initial step in the first loop. When the enzymes generated are subject to inhibition or inactivation, all four systems exhibit threshold properties akin to excitable systems like neuron firing. For those that are amenable to mathematical analysis, expressions are derived that relate the excitation threshold to the kinetics of enzyme generation and inhibition and the initial conditions. For the most complex system, it was expedient to employ numerical simulation to demonstrate threshold behavior, and in this case long-range feedback was seen to have two distinct effects. At sufficiently high catalytic rates, this feedback is capable of exciting an otherwise subthreshold system. At lower catalytic rates, where the long-range feedback does not significantly affect the threshold, it nonetheless has a major effect in potentiating the response above the threshold. In particular, oscillatory behavior observed in simulations of sequential feedback loops is abolished when a long-range feedback is present.

Calibração de redes de distribuição de água pelo método inverso aplicado a escoamento transiente

Desenvolve-se um método para estimar os parâmetros de uma rede hidráulica a partir de dados observados de cargas hidráulicas transientes. Os parâmetros físicos da rede como fatores de atrito, rugosidades absolutas, diâmetros e a identificação e quantificação de vazamentos são as grandezas desconhecidas. O problema transiente inverso é resolvido utilizando uma abordagem indireta que compara os dados disponíveis de carga hidráulica transiente observados com os calculados através de um método matemático. O Método Transiente Inverso (MTI) com um Algoritmo Genético (AG) emprega o Método das Características (MOC) na solução das equações do movimento para escoamento transiente em redes de tubos. As condições de regime permanente são desconhecidas. Para avaliar a confiabilidade do MTI-AG desenvolvido aqui, uma rede-exemplo é usada para os vários problemas de calibração propostos. O comportamento transiente é imposto por duas manobras distintas de uma válvula de controle localizada em um dos nós da rede. Analisam-se, ainda, o desempenho do método proposto mediante a variabilidade do tamanho do registro transiente e de possíveis erros de leitura nas cargas hidráulicas. Ensaios numéricos realizados mostram que o método é viável e aplicável à solução de problema inverso em redes hidráulicas, sobretudo recorrendo-se a poucos dados observados e ao desconhecimento das condições iniciais de estado permanente. Nos diversos problemas de identificação, as informações transientes obtidas da manobra mais brusca produziu estimações mais eficientes.

Aplicação e avaliação de desempenho de um sistema de otimização em tempo real em uma unidade de produção de propeno.

Com o objetivo de aumentar o lucro de plantas químicas, a Otimização em Tempo Real (RTO) é uma ferramenta que busca determinar as condições ótimas operacionais do processo em estado estacionário, respeitando as restrições operacionais estabelecidas. Neste trabalho foi realizada a implementação prática de um ciclo RTO em um processo de destilação por recompressão de vapor (VRD), propileno-propano, da Refinaria de Paulínia (Petrobras S.A.), a partir de dados históricos da planta. Foram consideradas as principais etapas de um ciclo clássico de RTO: identificação de estado estacionário, reconciliação de dados, estimação de parâmetros e otimização econômica. Essa unidade foi modelada, simulada e otimizada em EMSO (Environment for Modeling, Simulation and Optimization), um simulador de processos orientado a equações desenvolvido no Brasil. Foram analisados e comparados dois métodos de identificação de estado estacionário, um baseado no teste estatístico F e outro baseado em wavelets. Ambos os métodos tiveram resultados semelhantes e mostraram-se capazes de identificar os estados estacionários de forma satisfatória, embora seja necessário o ajuste de parâmetros na sua implementação. Foram identificados alguns pontos estacionários para serem submetidos ao ciclo RTO e foi possível verificar a importância de partir de um estado estacionário para a continuidade do ciclo, já que essa é uma premissa do método. A partir dos pontos analisados, os resultados deste estudo mostram que o RTO é capaz de aumentar o ganho econômico entre 2,5-24%, dependendo das condições iniciais consideradas, o que pode representar ganhos de até 18 milhões de dólares por ano. Além disso, para essa unidade, verificou-se que o compressor é um equipamento limitante no aumento de ganho econômico do processo.

Poland and Greece - two contrasting EU enlargement experiences. ZEI Discussion Paper C213, 2012

The aim of this contribution is a comparative analysis of the challenges Poland and Greece (and more broadly – CEE-10 and GIPS countries) had to face in the past as latecomers to the European Union and are facing now, in the aftermath of the world financial and economic crisis of 2008-09. The main underlying message conveyed in this text is two-fold. Firstly, the author is going to argue that the breadth and complexity of the challenges Poland and other CEE-10 countries had to face while entering the road of systemic transformation was by far greater compared to past and in particular – current problems of Greece (and the remaining GIPS countries) in the aftermath of the global financial and economic crisis of 2008-09. Secondly, a resilience of Poland and other CEE-10 economies, relative to Greece and other GIPS, to the recent crisis was due to a comparatively higher level of institutional development of the former group at the time of their EU accession and at present. The ensuing discussion is organized as follows. Section 2 below provides comparative background information on the two reference groups. In Section 3 we discuss the most salient features of the design of the command economy and its legacy, as a key determinant of the initial conditions of systemic transformation. Next, in Section 4 we overview the basic indicators of growth performance and institutional reforms in CEE-10 countries between 1990 and 2011. Section 5 offers a picture of economic growth and real economic convergence in Greece and the remaining GIPS countries. In Section 6 we embark on comparative analysis of the institutional quality of Greece and Poland against a broader background of GIPS, CEE-10 and the remaining EU member countries. Section 7 concludes with a summary of major findings.

Constraining storm-scale forecasts of deep convective initiation with surface weather observations

Thesis (Ph.D.)--University of Washington, 2016-06

Multiquantum eigenstates of a linear chain of coupled qubits

We present a technique to identify exact analytic expressions for the multiquantum eigenstates of a linear chain of coupled qubits. A choice of Hilbert subspaces is described that allows an exact solution of the stationary Schrodinger equation without imposing periodic boundary conditions and without neglecting end effects, fully including the dipole-dipole nearest-neighbor interaction between the atoms. The treatment is valid for an arbitrary coherent excitation in the atomic system, any number of atoms, any size of the chain relative to the resonant wavelength and arbitrary initial conditions of the atomic system. The procedure we develop is general enough to be adopted for the study of excitation in an arbitrary array of atoms including spin chains and one-dimensional Bose-Einstein condensates.

Finite-size effects in on-line learning of multilayer neural networks

We complement recent advances in thermodynamic limit analyses of mean on-line gradient descent learning dynamics in multi-layer networks by calculating fluctuations possessed by finite dimensional systems. Fluctuations from the mean dynamics are largest at the onset of specialisation as student hidden unit weight vectors begin to imitate specific teacher vectors, increasing with the degree of symmetry of the initial conditions. In light of this, we include a term to stimulate asymmetry in the learning process, which typically also leads to a significant decrease in training time.

Belief propagation vs. TAP for decoding corrupted messages

We employ two different methods, based on belief propagation and TAP,for decoding corrupted messages encoded by employing Sourlas's method, where the code word comprises products of K bits selected randomly from the original message. We show that the equations obtained by the two approaches are similar and provide the same solution as the one obtained by the replica approach in some cases K=2. However, we also show that for K>=3 and unbiased messages the iterative solution is sensitive to the initial conditions and is likely to provide erroneous solutions; and that it is generally beneficial to use Nishimori's temperature, especially in the case of biased messages.

The belief in TAP

We show the similarity between belief propagation and TAP, for decoding corrupted messages encoded by Sourlas's method. The latter is a special case of the Gallager error- correcting code, where the code word comprises products of K bits selected randomly from the original message. We examine the efficacy of solutions obtained by the two methods for various values of K and show that solutions for K>=3 may be sensitive to the choice of initial conditions in the case of unbiased patterns. Good approximations are obtained generally for K=2 and for biased patterns in the case of K>=3, especially when Nishimori's temperature is being used.

Ultra high-speed transmission using dispersion managed solitons

This thesis presents the results of numerical modelling of ultra high-speed transmission using DM solitons. The theory of propagation in optical fibres is presented with specific reference to optical communication systems. This theory is then expanded to. incorporate dispersion-managed transmission and the dispersion managed soliton. The first part of this work focuses on ultra high-speed dispersion managed soliton propagation in short period dispersion maps. Initially, the cbaracteristics .of dispersion managed soliton propagation in short period dispersion maps are contrasted to those of the more conventional dispersion managed regime. These properties are then utilised to investigate transmission at single channel data rates of 80 Gbit/s, 160 Gbit/s and 320 Gbit/s. For all three data rates, the tolerable limits for transmission over 1000 km, 3000 km and·transoceanic distances are defined. A major limitation of these higher bjt rate systems arises from the problem of noise-induced interactions, which is where the.accumulation of timing jitter causes neighbouring dispersion-managed solitons to interact. In addition, the systems become more sensitive to initial conditions as the data rate increases, .. The second part of the work focuses on contrasting the performance of a range of propagation regimes, from quasi-linear through to soliton-like propagation at 40 Gbit/s for both single channel and WDM dispersion managed transmission. The results indicated that whilst the optimal single channel performance was achieved for soliton-like propagation, the optimal WDM performance was achieved for propagation regime that lay between quasi-linear and soliton-like.

Studies on the deactivation and reactivation of spent hydroprocessing catalyst

A study is reported on the deactivation of hydroprocessing catalysts and their reactivation by the removal of coke and metal foulants. The literature on hydrotreating catalyst deactivation by coke and metals deposition, the environmental problems associated with spent catalyst disposal, and its reactivation/rejuvenation process were reviewed. Experimental studies on catalyst deactivation involved problem analysis in industrial hydroprocessing operations, through characterization of the spent catalyst, and laboratory coking studies. A comparison was made between the characteristics of spent catalysts from fixed bed and ebullating bed residue hydroprocessing reactor units and the catalyst deactivation pattern in both types of reactor systems was examined. In the laboratory the nature of initial coke deposited on the catalyst surface and its role on catalyst deactivation were studied. The influence of initial coke on catalyst surface area and porosity was significant. Both catalyst acidity and feedstock quality had a remarkable influence on the amount and the nature of the initial coke. The hydroenitrogenation function (HDN) of the catalyst was found to be deactivated more rapidly by the initial coke than the hydrodesulphurization function (HDS). In decoking experiments, special attention was paid to the initial conditions of coke combustion, since the early stages of contact between the coke on the spent catalyst surface and the oxygen are crucial in the decoking process. An increase in initial combustion temperature above 440oC and the oxygen content of the regeneration gas above 5% vanadium led to considerable sintering of the catalyst. At temperatures above 700oC there was a substantial loss of molybdenum from the catalyst, and phase transformations in the alumina support. The preferred leaching route (coked vs decoked form of spent catalyst) and a comparison of different reagents (i.e., oxalic acid and tartaric acid) and promoters (i.e., Hydrogen Peroxide and Ferric Nitrate) for better selectivity in removing the major foulant (vanadium), characterization and performance evaluation of the treated catalysts and modelling of the leaching process were addressed in spent catalyst rejuvenation studies. The surface area and pore volume increased substantially with increasing vanadium extraction from the spent catalyst; the HDS activity showed a parallel increase. The selectivity for leaching of vanadium deposits was better, and activity recovery was higher, for catalyst rejuvenated by metal leaching prior to decoking.

Investigation of chaotic mixing in laminar fluid systems by the use of computational fluid dynamics

This work presents significant development into chaotic mixing induced through periodic boundaries and twisting flows. Three-dimensional closed and throughput domains are shown to exhibit chaotic motion under both time periodic and time independent boundary motions, A property is developed originating from a signature of chaos, sensitive dependence to initial conditions, which successfully quantifies the degree of disorder withjn the mixing systems presented and enables comparisons of the disorder throughout ranges of operating parameters, This work omits physical experimental results but presents significant computational investigation into chaotic systems using commercial computational fluid dynamics techniques. Physical experiments with chaotic mixing systems are, by their very nature, difficult to extract information beyond the recognition that disorder does, does not of partially occurs. The initial aim of this work is to observe whether it is possible to accurately simulate previously published physical experimental results through using commercial CFD techniques. This is shown to be possible for simple two-dimensional systems with time periodic wall movements. From this, and subsequent macro and microscopic observations of flow regimes, a simple explanation is developed for how boundary operating parameters affect the system disorder. Consider the classic two-dimensional rectangular cavity with time periodic velocity of the upper and lower walls, causing two opposing streamline motions. The degree of disorder within the system is related to the magnitude of displacement of individual particles within these opposing streamlines. The rationale is then employed in this work to develop and investigate more complex three-dimensional mixing systems that exhibit throughputs and time independence and are therefore more realistic and a significant advance towards designing chaotic mixers for process industries. Domains inducing chaotic motion through twisting flows are also briefly considered. This work concludes by offering possible advancements to the property developed to quantify disorder and suggestions of domains and associated boundary conditions that are expected to produce chaotic mixing.

Solids mixing and segregation in a gas spout-fluid bed-effect of the distributor design

This work is concerned with the assessment of a newer version of the spout-fluid bed where the gas is supplied from a common plenum and the distributor controls the operational phenomenon. Thus the main body of the work deals with the effect of the distributor design on the mixing and segregation of solids in a spout-filled bed. The effect of distributor design in the conventional fluidised bed and of variation of the gas inlet diameter in a spouted bed were also briefly investigated for purpose of comparison. Large particles were selected for study because they are becoming increasingly important in industrial fluidised beds but have not been thoroughly investigated. The mean particle diameters of the fraction ranged from 550 to 2400 mm, and their specific gravity from 0.97 to 2.45. Only work carried out with binary systems is reported here. The effect of air velocity, particle properties, bed height, the relative amount of jetsam and flotsam and initial conditions on the steady-state concentration profiles were assessed with selected distributors. The work is divided into three sections. Sections I and II deal with the fluidised bed and spouted bed systems. Section III covers the development of the spout-filled bed and its behaviour with reference to distributor design and it is shown how benefits of both spouting and fluidising phenomena can be exploited. In the fluidisation zone, better mixing is achieved by distributors which produce a large initial bubble diameter. Some common features exist between the behaviour of unidensity jetsam-rich systems and different density flotsam-rich systems. The shape factor does not seem to have an affect as long as it is only restricted to the minor component. However, in the case of the major component, particle shape significantly affects the final results. Studies of aspect ratio showed that there is a maximum (1.5) above which slugging occurs and the effect of the distributor design is nullified. A mixing number was developed for unidensity spherical rich systems, which proved to be extremely useful in quantifying the variation in mixing and segregation with changes in distributor design.

Nonlinear spectral shaping and optical rogue events in fiber-based systems

We provide an overview of our recent work on the shaping and stability of optical continua in the long pulse regime. Fibers with normal group-velocity dispersion at all-wavelengths are shown to allow for highly coherent continua that can be nonlinearly shaped using appropriate initial conditions. In contrast, supercontinua generated in the anomalous dispersion regime are shown to exhibit large fluctuations in the temporal and spectral domains that can be controlled using a carefully chosen seed. A particular example of this is the first experimental observation of the Peregrine soliton which constitutes a prototype of optical rogue-waves.

Parabolic optical pulses under the action of the third-order dispersion

Recent developments in nonlinear optics reveal an interesting class of pulses with a parabolic intensity profile in the energy-containing core and a linear frequency chirp that can propagate in a fiber with normal group-velocity dispersion. Parabolic pulses propagate in a stable selfsimilar manner, holding certain relations (scaling) between pulse power, width, and chirp parameter. In the additional presence of linear amplification, they enjoy the remarkable property of representing a common asymptotic state (or attractor) for arbitrary initial conditions. Analytically, self-similar (SS) parabolic pulses can be found as asymptotic, approximate solutions of the nonlinear Schr¨odinger equation (NLSE) with gain in the semi-classical (largeamplitude/small-dispersion) limit. By analogy with the well-known stable dynamics of solitary waves - solitons, these SS parabolic pulses have come to be known as similaritons. In practical fiber systems, inherent third-order dispersion (TOD) in the fiber always introduces a certain degree of asymmetry in the structure of the propagating pulse, eventually leading to pulse break-up. To date, there is no analytic theory of parabolic pulses under the action of TOD. Here, we develop aWKB perturbation analysis that describes the effect of weak TOD on the parabolic pulse solution of the NLSE in a fiber gain medium. The induced perturbation in phase and amplitude can be found to any order. The theoretical model predicts with sufficient accuracy the pulse structural changes induced by TOD, which are observed through direct numerical NLSE simulations.