Bifurcation Analysis of Reaction Diffusion Systems on Arbitrary Surfaces

In this article we present a computational framework for isolating spatial patterns arising in the steady states of reaction-diffusion systems. Such systems have been used to model many different phenomena in areas such as developmental and cancer biology, cell motility and material science. Often one is interested in identifying parameters which will lead to a particular pattern. To attempt to answer this, we compute eigenpairs of the Laplacian on a variety of domains and use linear stability analysis to determine parameter values for the system that will lead to spatially inhomogeneous steady states whose patterns correspond to particular eigenfunctions. This method has previously been used on domains and surfaces where the eigenvalues and eigenfunctions are found analytically in closed form. Our contribution to this methodology is that we numerically compute eigenpairs on arbitrary domains and surfaces. Here we present various examples and demonstrate that mode isolation is straightforward especially for low eigenvalues. Additionally we see that if two or more eigenvalues are in a permissible range then the inhomogeneous steady state can be a linear combination of the respective eigenfunctions. Finally we show an example which suggests that pattern formation is robust on similar surfaces in cases that the surface either has or does not have a boundary.

Multipartite entanglement and quantum state exchange

We investigate multipartite entanglement in relation to the process of quantum state exchange. In particular, we consider such entanglement for a certain pure state involving two groups of N trapped atoms. The state, which can be produced via quantum state exchange, is analogous to the steady-state intracavity state of the subthreshold optical nondegenerate parametric amplifier. We show that, first, it possesses some 2N-way entanglement. Second, we place a lower bound on the amount of such entanglement in the state using a measure called the entanglement of minimum bipartite entropy.

Coherent state construction of representations of osp(2|2) and primary fields of osp(2|2) conformal field theory

Representations of the superalgebra osp(2/2)(k)((1)) and current superalgebra. osp(2/2)k in the standard basis are investigated. All finite-dimensional typical and atypical representations of osp(2/2) are constructed by the vector coherent state method. Primary fields of the non-unitary conformal field theory associated with osp(2/2)(k)((1)) in the standard basis are obtained for arbitrary level k. (C) 2004 Elsevier B.V. All rights reserved.

Doxorubicin pharmacokinetics following a single-dose infusion to sulphur-crested cockatoos (Cacatua galerita)

Objective To determine the pharmacokinetics of doxorubicin in sulphur-crested cockatoos, so that its use in clinical studies in birds can be considered. Design A pharmacokinetic study of doxorubicin, following a single intravenous (IV) infusion over 20 min, was performed in four healthy sulphur-crested cockatoos (Cacatua galerita). Procedure Birds were anaesthetised and both jugular veins were cannulated, one for doxorubicin infusion and the other for blood collection. Doxorubicin hydrochloride (2 mg/kg) in normal saline was infused IV over 20 min at a constant rate. Serial blood samples were collected for 96 h after initiation of the infusion. Plasma doxorubicin concentrations were assayed using an HPLC method involving ethyl acetate extraction, reverse-phase chromatography and fluorescence detection. The limit of quantification was 20 ng/mL. Established non-parametric methods were used for the analysis of plasma doxorubicin data. Results During the infusion the mean +/- SD for the C-max of doxorubicin was 4037 +/- 2577 ng/mL. Plasma concentrations declined biexponentially immediately after the infusion was ceased. There was considerable intersubject variability in all pharmacokinetic variables. The terminal (beta-phase) half-life was 41.4 +/- 18.5 min, the systemic clearance (Cl) was 45.7 +/- 18.0 mL/min/kg, the mean residence time (MRT) was 4.8 +/- 1.4 min, and the volume of distribution at steady state (V-SS) was 238 131 mL/kg. The extrapolated area under the curve (AUC(0-infinity)) was 950 +/- 677 ng/mL.h. The reduced metabolite, doxorubicinol, was detected in the plasma of all four parrots but could be quantified in only one bird with the profile suggesting formation rate-limited pharmacokinetics of doxorubicinol. Conclusions and clinical relevance Doxorubicin infusion in sulphur-crested cockatoos produced mild, transient inappetence. The volume of distribution per kilogram and terminal half-life were considerably smaller, but the clearance per kilogram was similar to or larger than reported in the dog, rat and humans. Traces of doxorubicinol, a metabolite of doxorubicin, were detected in the plasma.

A numerical approach to modeling the catalytic voltammetry of surface-confined redox enzymes

A finite difference method for simulating voltammograms of electrochemically driven enzyme catalysis is presented. The method enables any enzyme mechanism to be simulated. The finite difference equations can be represented as a matrix equation containing a nonlinear sparse matrix. This equation has been solved using the software package Mathematica. Our focus is on the use of cyclic voltammetry since this is the most commonly employed electrochemical method used to elucidate mechanisms. The use of cyclic voltammetry to obtain data from systems obeying Michaelis-Menten kinetics is discussed, and we then verify our observations on the Michaelis-Menten system using the finite difference simulation. Finally, we demonstrate how the method can be used to obtain mechanistic information on a real redox enzyme system, the complex bacterial molybdoenzyme xanthine dehydrogenase.

Overexpressed Ca-v beta 3 inhibits N-type (Ca(v)2.2) calcium channel currents through a hyperpolarizing shift of "ultra-slow" and "closed-state" inactivation

It has been shown that P auxiliary subunits increase current amplitude in voltage-dependent calcium channels. In this study, however, we found a hovel inhibitory effect of beta3 Subunit on macroscopic Ba2+ currents through recombinant N- and R-type calcium channels expressed in Xenopus oocytes. Overexpressed beta3 (12.5 ng/ cell cRNA) significantly suppressed N- and R-type, but not L-type, calcium channel currents at physiological holding potentials (HPs) of -60 and -80 mV At a HP of -80 mV, coinjection of various concentrations (0-12.5 ng) of the beta3 with Ca,.2.2alpha(1) and alpha(2)delta enhanced the maximum conductance of expressed channels at lower beta3 concentrations but at higher concentrations (>2.5 ng/cell) caused a marked inhibition. The beta3-induced Current suppression was reversed at a HP of - 120 mV, suggesting that the inhibition was voltage dependent. A high concentration of Ba-2divided by (40 mM) as a charge carrier also largely diminished the effect of P3 at -80 mV Therefore, experimental conditions (HP, divalent cation concentration, and P3 subunit concentration) approaching normal physiological conditions were critical to elucidate the full extent of this novel P3 effect. Steady-state inactivation curves revealed that N-type channels exhibited closed-state inactivation without P3, and that P3 caused an similar to40 mV negative shift of the inactivation, producing a second component with an inactivation midpoint of approximately -85 mV The inactivation of N-type channels in the presence of a high concentration (12.5 ng/cell) of P3 developed slowly and the time-dependent inactivation curve was best fit by the sum of two exponential functions with time constants of 14 s and 8.8 min at -80 mV Similar ultra-slow inactivation was observed for N-type channels Without P3. Thus, P3 can have a profound negative regulatory effect on N-type (and also R-type) calcium channels by Causing a hyperpolarizing shift of the inactivation without affecting ultra-slow and closed-state inactivation properties.

Prediction of slurry transport in SAG mills using SPH fluid flow in a dynamic DEM based porous media

DEM modelling of the motion of coarse fractions of the charge inside SAG mills has now been well established for more than a decade. In these models the effect of slurry has broadly been ignored due to its complexity. Smoothed particle hydrodynamics (SPH) provides a particle based method for modelling complex free surface fluid flows and is well suited to modelling fluid flow in mills. Previous modelling has demonstrated the powerful ability of SPH to capture dynamic fluid flow effects such as lifters crashing into slurry pools, fluid draining from lifters, flow through grates and pulp lifter discharge. However, all these examples were limited by the ability to model only the slurry in the mill without the charge. In this paper, we represent the charge as a dynamic porous media through which the SPH fluid is then able to flow. The porous media properties (specifically the spatial distribution of porosity and velocity) are predicted by time averaging the mill charge predicted using a large scale DEM model. This allows prediction of transient and steady state slurry distributions in the mill and allows its variation with operating parameters, slurry viscosity and slurry volume, to be explored. (C) 2006 Published by Elsevier Ltd.

Entangled-state cycles from conditional quantum evolution

A system of cascaded qubits interacting via the one-way exchange of photons is studied. While for general operating conditions the system evolves to a superposition of Bell states (a dark state) in the long-time limit, under a particular resonance condition no steady state is reached within a finite time. We analyze the conditional quantum evolution (quantum trajectories) to characterize the asymptotic behavior under this resonance condition. A distinct bimodality is observed: for perfect qubit coupling, the system either evolves to a maximally entangled Bell state without emitting photons (the dark state) or executes a sustained entangled-state cycle-random switching between a pair of Bell states while emitting a continuous photon stream; for imperfect coupling, two entangled-state cycles coexist, between which a random selection is made from one quantum trajectory to another.

Instantaneous RRMI of iron rods in reduced gravity

Standard upward-burning promoted ignition tests (“Standard Test Method for Determining the Combustion Behavior of Metallic Materials in Oxygen-Enriched Atmospheres,” ASTM G4-124 [1] or “Flammability, Odor, Offgassing, and Compatibility Requirements and Test Procedures for Materials in Environments that Support Combustion,” NASA-STD-6001, NASA Test 17 [2]) were performed on cylindrical iron (99.95% pure) rods in various oxygen purities (95.0–99.98%) in reduced gravity onboard NASA JSC's KC-135 to investigate the effect of gravity on the regression rate of the melting interface. Visual analysis of experiments agrees with previous published observations showing distinct motions of the molten mass attached to the solid rod during testing. Using an ultrasonic technique to record the real-time rod length, comparison of the instantaneous regression rate of the melting interface and visual recording shows a non-steady-state regression rate of the melting interface for the duration of a test. Precessional motion is associated with a higher regression rate of the melting interface than for test periods in which the molten mass does not show lateral motion. The transition between the two types of molten mass motion during a test was accompanied by a reduced regression rate of the melting interface, approximately 15–50% of the average regression rate of the melting interface for the entire test.

Indagine modellistica sulla tenuta idraulica delle difese arginali deputate alla protezione idraulica del territorio

La prima parte di questo lavoro di tesi tratta dell’interazione tra un bacino di laminazione e il sottostante acquifero: è in fase di progetto, infatti, la costruzione di una cassa di espansione sul torrente Baganza, a monte della città di Parma. L’obiettivo di tale intervento è di ridurre il rischio di esondazione immagazzinando temporaneamente, in un serbatoio artificiale, la parte più pericolosa del volume di piena che verrebbe rilasciata successivamente con portate che possono essere agevolmente contenute nel tratto cittadino del torrente. L’acquifero è stato preliminarmente indagato e monitorato permettendone la caratterizzazione litostratigrafica. La stratigrafia si può riassumere in una sequenza di strati ghiaioso-sabbiosi con successione di lenti d’argilla più o meno spesse e continue, distinguendo due acquiferi differenti (uno freatico ed uno confinato). Nel presente studio si fa riferimento al solo acquifero superficiale che è stato modellato numericamente, alle differenze finite, per mezzo del software MODFLOW_2005. L'obiettivo del presente lavoro è di rappresentare il sistema acquifero nelle condizioni attuali (in assenza di alcuna opera) e di progetto. La calibrazione è stata condotta in condizioni stazionarie utilizzando i livelli piezometrici raccolti nei punti d’osservazione durante la primavera del 2013. I valori di conducibilità idraulica sono stati stimati per mezzo di un approccio geostatistico Bayesiano. Il codice utilizzato per la stima è il bgaPEST, un software gratuito per la soluzione di problemi inversi fortemente parametrizzati, sviluppato sulla base dei protocolli del software PEST. La metodologia inversa stima il campo di conducibilità idraulica combinando osservazioni sullo stato del sistema (livelli piezometrici nel caso in esame) e informazioni a-priori sulla struttura dei parametri incogniti. La procedura inversa richiede il calcolo della sensitività di ciascuna osservazione a ciascuno dei parametri stimati; questa è stata valutata in maniera efficiente facendo ricorso ad una formulazione agli stati aggiunti del codice in avanti MODFLOW_2005_Adjoint. I risultati della metodologia sono coerenti con la natura alluvionale dell'acquifero indagato e con le informazioni raccolte nei punti di osservazione. Il modello calibrato può quindi essere utilizzato come supporto alla progettazione e gestione dell’opera di laminazione. La seconda parte di questa tesi tratta l'analisi delle sollecitazioni indotte dai percorsi di flusso preferenziali causati da fenomeni di piping all’interno dei rilevati arginali. Tali percorsi preferenziali possono essere dovuti alla presenza di gallerie scavate da animali selvatici. Questo studio è stato ispirato dal crollo del rilevato arginale del Fiume Secchia (Modena), che si è verificato in gennaio 2014 a seguito di un evento alluvionale, durante il quale il livello dell'acqua non ha mai raggiunto la sommità arginale. La commissione scientifica, la cui relazione finale fornisce i dati utilizzati per questo studio, ha attribuito, con molta probabilità, il crollo del rilevato alla presenza di tane di animali. Con lo scopo di analizzare il comportamento del rilevato in condizioni integre e in condizioni modificate dall'esistenza di un tunnel che attraversa il manufatto arginale, è stato realizzato un modello numerico 3D dell’argine mediante i noti software Femwater e Feflow. I modelli descrivono le infiltrazioni all'interno del rilevato considerando il terreno in entrambe le porzioni sature ed insature, adottando la tecnica agli elementi finiti. La tana è stata rappresentata da elementi con elevata permeabilità e porosità, i cui valori sono stati modificati al fine di valutare le diverse influenze sui flussi e sui contenuti idrici. Per valutare se le situazioni analizzate presentino o meno il verificarsi del fenomeno di erosione, sono stati calcolati i valori del fattore di sicurezza. Questo è stato valutato in differenti modi, tra cui quello recentemente proposto da Richards e Reddy (2014), che si riferisce al criterio di energia cinetica critica. In ultima analisi è stato utilizzato il modello di Bonelli (2007) per calcolare il tempo di erosione ed il tempo rimanente al collasso del rilevato.

Jobs, Unemployment, and Macroeconomic Transmission

Measuring Job Openings: Evidence from Swedish Plant Level Data. In modern macroeconomic models “job openings'' are a key component. Thus, when taking these models to the data we need an empirical counterpart to the theoretical concept of job openings. To achieve this, the literature relies on job vacancies measured either in survey or register data. Insofar as this concept captures the concept of job openings well we should see a tight relationship between vacancies and subsequent hires on the micro level. To investigate this, I analyze a new data set of Swedish hires and job vacancies on the plant level covering the period 2001-2012. I find that vacancies contain little power in predicting hires over and above (i) whether the number of vacancies is positive and (ii) plant size. Building on this, I propose an alternative measure of job openings in the economy. This measure (i) better predicts hiring at the plant level and (ii) provides a better fitting aggregate matching function vis-à-vis the traditional vacancy measure. Firm Level Evidence from Two Vacancy Measures. Using firm level survey and register data for both Sweden and Denmark we show systematic mis-measurement in both vacancy measures. While the register-based measure on the aggregate constitutes a quarter of the survey-based measure, the latter is not a super-set of the former. To obtain the full set of unique vacancies in these two databases, the number of survey vacancies should be multiplied by approximately 1.2. Importantly, this adjustment factor varies over time and across firm characteristics. Our findings have implications for both the search-matching literature and policy analysis based on vacancy measures: observed changes in vacancies can be an outcome of changes in mis-measurement, and are not necessarily changes in the actual number of vacancies. Swedish Unemployment Dynamics. We study the contribution of different labor market flows to business cycle variations in unemployment in the context of a dual labor market. To this end, we develop a decomposition method that allows for a distinction between permanent and temporary employment. We also allow for slow convergence to steady state which is characteristic of European labor markets. We apply the method to a new Swedish data set covering the period 1987-2012 and show that the relative contributions of inflows and outflows to/from unemployment are roughly 60/30. The remaining 10\% are due to flows not involving unemployment. Even though temporary contracts only cover 9-11\% of the working age population, variations in flows involving temporary contracts account for 44\% of the variation in unemployment. We also show that the importance of flows involving temporary contracts is likely to be understated if one does not account for non-steady state dynamics. The New Keynesian Transmission Mechanism: A Heterogeneous-Agent Perspective. We argue that a 2-agent version of the standard New Keynesian model---where a ``worker'' receives only labor income and a “capitalist'' only profit income---offers insights about how income inequality affects the monetary transmission mechanism. Under rigid prices, monetary policy affects the distribution of consumption, but it has no effect on output as workers choose not to change their hours worked in response to wage movements. In the corresponding representative-agent model, in contrast, hours do rise after a monetary policy loosening due to a wealth effect on labor supply: profits fall, thus reducing the representative worker's income. If wages are rigid too, however, the monetary transmission mechanism is active and resembles that in the corresponding representative-agent model. Here, workers are not on their labor supply curve and hence respond passively to demand, and profits are procyclical.

Amplitude response of a unilaterally constrained nonlinear micromechanical resonator

Dynamical systems that involve impacts frequently arise in engineering. This Letter reports a study of such a system at microscale that consists of a nonlinear resonator operating with an unilateral impact. The microresonators were fabricated on silicon-on-insulator wafers by using a one-mask process and then characterised by using the capacitively driving and sensing method. Numerical results concerning the dynamics of this vibro-impact system were verified by the experiments. Bifurcation analysis was used to provide a qualitative scenario of the system steady-state solutions as a function of both the amplitude and the frequency of the external driving sinusoidal voltage. The results show that the amplitude of resonant peak is levelled off owing to the impact effect and that the bandwidth of impacting is dependent upon the nonlinearity and the operating conditions.

Computer control of chemical process plants with special reference to distillation

A sieve plate distillation column has been constructed and interfaced to a minicomputer with the necessary instrumentation for dynamic, estimation and control studies with special bearing on low-cost and noise-free instrumentation. A dynamic simulation of the column with a binary liquid system has been compiled using deterministic models that include fluid dynamics via Brambilla's equation for tray liquid holdup calculations. The simulation predictions have been tested experimentally under steady-state and transient conditions. The simulator's predictions of the tray temperatures have shown reasonably close agreement with the measured values under steady-state conditions and in the face of a step change in the feed rate. A method of extending linear filtering theory to highly nonlinear systems with very nonlinear measurement functional relationships has been proposed and tested by simulation on binary distillation. The simulation results have proved that the proposed methodology can overcome the typical instability problems associated with the Kalman filters. Three extended Kalman filters have been formulated and tested by simulation. The filters have been used to refine a much simplified model sequentially and to estimate parameters such as the unmeasured feed composition using information from the column simulation. It is first assumed that corrupted tray composition measurements are made available to the filter and then corrupted tray temperature measurements are accessed instead. The simulation results have demonstrated the powerful capability of the Kalman filters to overcome the typical hardware problems associated with the operation of on-line analyzers in relation to distillation dynamics and control by, in effect, replacirig them. A method of implementing estimator-aided feedforward (EAFF) control schemes has been proposed and tested by simulation on binary distillation. The results have shown that the EAFF scheme provides much better control and energy conservation than the conventional feedback temperature control in the face of a sustained step change in the feed rate or multiple changes in the feed rate, composition and temperature. Further extensions of this work are recommended as regards simulation, estimation and EAFF control.

Simultaneous chemical reaction and distillation of formaldehyde

The thesis is concerned with the development and testing of a mathematical model of a distillation process in which the components react chemically. The formaldehyde-methanol-water system was selected and only the reversible reactions between formaldehyde and water giving methylene glycol and between formaldehyde and methanol producing hemiformal were assumed to occur under the distillation conditions. Accordingly the system has been treated as a five component system. The vapour-liquid equilibrium calculations were performed by solving iteratively the thermodynamic relationships expressing the phase equilibria with the stoichiometric equations expressing the chemical equilibria. Using optimisation techniques, the Wilson single parameters and Henry's constants were calculated for binary systems containing formaldehyde which was assumed to be a supercritical component whilst Wilson binary parameters were calculated for the remaining binary systems. Thus the phase equilibria for the formaldehyde system could be calculated using these parameters and good accuracy was obtained when calculated values were compared with experimental values. The distillation process was modelled using the mass and energy balance equations together with the phase equilibria calculations. The plate efficiencies were obtained from a modified A.I.Ch.E. Bubble Tray method. The resulting equations were solved by an iterative plate to plate calculation based on the Newton Raphson method. Experiments were carried out in a 76mm I.D., eight sieve plate distillation column and the results were compared with the mathematical model calculations. Overall, good agreement was obtained but some discrepancies were observed in the concentration profiles and these may have been caused by the effect of limited physical property data and a limited understanding of the reactions mechanism. The model equations were solved in the form of modular computer programs. Although they were written to describe the steady state distillation with simultaneous chemical reaction of the formaldehyde system, the approach used may be of wider application.

Mathematical methods in power system stability studies

In this thesis various mathematical methods of studying the transient and dynamic stabiIity of practical power systems are presented. Certain long established methods are reviewed and refinements of some proposed. New methods are presented which remove some of the difficulties encountered in applying the powerful stability theories based on the concepts of Liapunov. Chapter 1 is concerned with numerical solution of the transient stability problem. Following a review and comparison of synchronous machine models the superiority of a particular model from the point of view of combined computing time and accuracy is demonstrated. A digital computer program incorporating all the synchronous machine models discussed, and an induction machine model, is described and results of a practical multi-machine transient stability study are presented. Chapter 2 reviews certain concepts and theorems due to Liapunov. In Chapter 3 transient stability regions of single, two and multi~machine systems are investigated through the use of energy type Liapunov functions. The treatment removes several mathematical difficulties encountered in earlier applications of the method. In Chapter 4 a simple criterion for the steady state stability of a multi-machine system is developed and compared with established criteria and a state space approach. In Chapters 5, 6 and 7 dynamic stability and small signal dynamic response are studied through a state space representation of the system. In Chapter 5 the state space equations are derived for single machine systems. An example is provided in which the dynamic stability limit curves are plotted for various synchronous machine representations. In Chapter 6 the state space approach is extended to multi~machine systems. To draw conclusions concerning dynamic stability or dynamic response the system eigenvalues must be properly interpreted, and a discussion concerning correct interpretation is included. Chapter 7 presents a discussion of the optimisation of power system small sjgnal performance through the use of Liapunov functions.