Onset And Characterisation Of Chaos In A Two Dimensional Nonlinear Deterministic Model

Nature is full of phenomena which we call "chaotic", the weather being a prime example. What we mean by this is that we cannot predict it to any significant accuracy, either because the system is inherently complex, or because some of the governing factors are not deterministic. However, during recent years it has become clear that random behaviour can occur even in very simple systems with very few number of degrees of freedom, without any need for complexity or indeterminacy. The discovery that chaos can be generated even with the help of systems having completely deterministic rules - often models of natural phenomena - has stimulated a lo; of research interest recently. Not that this chaos has no underlying order, but it is of a subtle kind, that has taken a great deal of ingenuity to unravel. In the present thesis, the author introduce a new nonlinear model, a ‘modulated’ logistic map, and analyse it from the view point of ‘deterministic chaos‘.

Development of Electrochemical and Fluorescent Sensors

Chemical sensors have growing interest in the determination of food additives, which are creating toxicity and may cause serious health concern, drugs and metal ions. A chemical sensor can be defined as a device that transforms chemical information, ranging from the concentration of a specific sample component to total composition analysis, into an analytically useful signal. The chemical information may be generated from a chemical reaction of the analyte or from a physical property of the system investigated. Two main steps involved in the functioning of a chemical sensor are recognition and transduction. Chemical sensors employ specific transduction techniques to yield analyte information. The most widely used techniques employed in chemical sensors are optical absorption, luminescence, redox potential etc. According to the operating principle of the transducer, chemical sensors may be classified as electrochemical sensors, optical sensors, mass sensitive sensors, heat sensitive sensors etc. Electrochemical sensors are devices that transform the effect of the electrochemical interaction between analyte and electrode into a useful signal. They are very widespread as they use simple instrumentation, very good sensitivity with wide linear concentration ranges, rapid analysis time and simultaneous determination of several analytes. These include voltammetric, potentiometric and amperometric sensors. Fluorescence sensing of chemical and biochemical analytes is an active area of research. Any phenomenon that results in a change of fluorescence intensity, anisotropy or lifetime can be used for sensing. The fluorophores are mixed with the analyte solution and excited at its corresponding wavelength. The change in fluorescence intensity (enhancement or quenching) is directly related to the concentration of the analyte. Fluorescence quenching refers to any process that decreases the fluorescence intensity of a sample. A variety of molecular rearrangements, energy transfer, ground-state complex formation and collisional quenching. Generally, fluorescence quenching can occur by two different mechanisms, dynamic quenching and static quenching. The thesis presents the development of voltammetric and fluorescent sensors for the analysis of pharmaceuticals, food additives metal ions. The developed sensors were successfully applied for the determination of analytes in real samples. Chemical sensors have multidisciplinary applications. The development and application of voltammetric and optical sensors continue to be an exciting and expanding area of research in analytical chemistry. The synthesis of biocompatible fluorophores and their use in clinical analysis, and the development of disposable sensors for clinical analysis is still a challenging task. The ability to make sensitive and selective measurements and the requirement of less expensive equipment make electrochemical and fluorescence based sensors attractive.

On Parallel Communicating Grammar Systems and Correctness Preserving Restarting Automata

This paper contributes to the study of Freely Rewriting Restarting Automata (FRR-automata) and Parallel Communicating Grammar Systems (PCGS), which both are useful models in computational linguistics. For PCGSs we study two complexity measures called 'generation complexity' and 'distribution complexity', and we prove that a PCGS Pi, for which the generation complexity and the distribution complexity are both bounded by constants, can be transformed into a freely rewriting restarting automaton of a very restricted form. From this characterization it follows that the language L(Pi) generated by Pi is semi-linear, that its characteristic analysis is of polynomial size, and that this analysis can be computed in polynomial time.

MA230exam00

Exam questions and solutions in LaTex

Radar studies of the vertical distribution of insects migrating over southern Britain: the influence of temperature inversions on nocturnal layer concentrations

Insects migrating over two sites in southern UK (Malvern in Worcestershire, and Harpenden in Hertfordshire) have been monitored continuously with nutating vertical-looking radars (VLRs) equipped with powerful control and analysis software. These observations make possible, for the first time, a systematic investigation of the vertical distribution of insect aerial density in the atmosphere, over temporal scales ranging from the short (instantaneous vertical profiles updated every 15 min) to the very long (profiles aggregated over whole seasons or even years). In the present paper, an outline is given of some general features of insect stratification as revealed by the radars, followed by a description of occasions during warm nights in the summer months when intense insect layers developed. Some of these nocturnal layers were due to the insects flying preferentially at the top of strong surface temperature inversions, and in other cases, layering was associated with higher-altitude temperature maxima, such as those due to subsidence inversions. The layers were formed from insects of a great variety of sizes, but peaks in the mass distributions pointed to a preponderance of medium-sized noctuid moths on certain occasions.

Stability of stationary solutions of the forced Navier-Stokes equations on the two-torus

We study the linear and nonlinear stability of stationary solutions of the forced two-dimensional Navier-Stokes equations on the domain [0,2π]x[0,2π/α], where α ϵ(0,1], with doubly periodic boundary conditions. For the linear problem we employ the classical energy{enstrophy argument to derive some fundamental properties of unstable eigenmodes. From this it is shown that forces of pure χ2-modes having wavelengths greater than 2π do not give rise to linear instability of the corresponding primary stationary solutions. For the nonlinear problem, we prove the equivalence of nonlinear stability with respect to the energy and enstrophy norms. This equivalence is then applied to derive optimal conditions for nonlinear stability, including both the high-and low-Reynolds-number limits.

Interface dynamics in planar neural field models

Neural field models describe the coarse-grained activity of populations of interacting neurons. Because of the laminar structure of real cortical tissue they are often studied in two spatial dimensions, where they are well known to generate rich patterns of spatiotemporal activity. Such patterns have been interpreted in a variety of contexts ranging from the understanding of visual hallucinations to the generation of electroencephalographic signals. Typical patterns include localized solutions in the form of traveling spots, as well as intricate labyrinthine structures. These patterns are naturally defined by the interface between low and high states of neural activity. Here we derive the equations of motion for such interfaces and show, for a Heaviside firing rate, that the normal velocity of an interface is given in terms of a non-local Biot-Savart type interaction over the boundaries of the high activity regions. This exact, but dimensionally reduced, system of equations is solved numerically and shown to be in excellent agreement with the full nonlinear integral equation defining the neural field. We develop a linear stability analysis for the interface dynamics that allows us to understand the mechanisms of pattern formation that arise from instabilities of spots, rings, stripes and fronts. We further show how to analyze neural field models with linear adaptation currents, and determine the conditions for the dynamic instability of spots that can give rise to breathers and traveling waves.

The stability of a two-dimensional vorticity filament under uniform strain

The quantitative effects of uniform strain and background rotation on the stability of a strip of constant vorticity (a simple shear layer) are examined. The thickness of the strip decreases in time under the strain, so it is necessary to formulate the linear stability analysis for a time-dependent basic flow. The results show that even a strain rate γ (scaled with the vorticity of the strip) as small as 0.25 suppresses the conventional Rayleigh shear instability mechanism, in the sense that the r.m.s. wave steepness cannot amplify by more than a certain factor, and must eventually decay. For γ < 0.25 the amplification factor increases as γ decreases; however, it is only 3 when γ e 0.065. Numerical simulations confirm the predictions of linear theory at small steepness and predict a threshold value necessary for the formation of coherent vortices. The results help to explain the impression from numerous simulations of two-dimensional turbulence reported in the literature that filaments of vorticity infrequently roll up into vortices. The stabilization effect may be expected to extend to two- and three-dimensional quasi-geostrophic flows.

Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

Phytoplankton–zooplankton dynamics in the ‘presence’ or ‘absence’ of toxic phytoplankton

The inhibitory effects of toxin-producing phytoplankton (TPP) on zooplankton modulate the dynamics of marine plankton. In this article, we employ simple mathematical models to compare theoretically the dynamics of phytoplankton–zooplankton interaction in situations where the TPP are present with those where TPP are absent. We consider two sets of three-component interaction models: one that does not include the effect of TPP and the other that does. The negative effects of TPP on zooplankton is described by a non-linear interaction term. Extensive theoretical analyses of the models have been performed to understand the qualitative behaviour of the model systems around every possible equilibria. The results of local-stability analysis and numerical simulations demonstrate that the two model-systems differ qualitatively with regard to oscillations and stability. The model system that does not include TPP is asymptotically stable around the coexisting equilibria, whereas, the system that includes TPP oscillates for a range of parametric values associated with toxin-inhibition rate and competition coefficients. Our analysis suggests that the qualitative dynamics of the plankton–zooplankton interactions are very likely to alter due to the presence of TPP species, and therefore the effects of TPP should be considered carefully while modelling plankton dynamics.

Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Tyrosine hydroxylase immunoreactivity as indicator of sympathetic activity: simultaneous evaluation in different tissues of hypertensive rats

Burgi K, Cavalleri MT, Alves AS, Britto LRG, Antunes VR, Michelini LC. Tyrosine hydroxylase immunoreactivity as indicator of sympathetic activity: simultaneous evaluation in different tissues of hypertensive rats. Am J Physiol Regul Integr Comp Physiol 300: R264-R271, 2011. First published December 9, 2010; doi: 10.1152/ajpregu.00687.2009.-Vasomotor control by the sympathetic nervous system presents substantial heterogeneity within different tissues, providing appropriate homeostatic responses to maintain basal/stimulated cardiovascular function both at normal and pathological conditions. The availability of a reproducible technique for simultaneous measurement of sympathetic drive to different tissues is of great interest to uncover regional patterns of sympathetic nerve activity (SNA). We propose the association of tyrosine hydroxylase immunoreactivity (THir) with image analysis to quantify norepinephrine (NE) content within nerve terminals in arteries/arterioles as a good index for regional sympathetic outflow. THir was measured in fixed arterioles of kidney, heart, and skeletal muscle of WistarKyoto rats (WKY) and spontaneously hypertensive rats (SHR) (123 +/- 2 and 181 +/- 4 mmHg, 300 +/- 8 and 352 +/- 8 beats/min, respectively). There was a differential THir distribution in both groups: higher THir was observed in the kidney and skeletal muscle (similar to 3-4-fold vs. heart arterioles) of WKY; in SHR, THir was increased in the kidney and heart (2.4- and 5.3-fold vs. WKY, respectively) with no change in the skeletal muscle arterioles. Observed THir changes were confirmed by either: 1) determination of NE content (high-performance liquid chromatography) in fresh tissues (SHR vs. WKY): +34% and +17% in kidney and heart, respectively, with no change in the skeletal muscle; 2) direct recording of renal (RSNA) and lumbar SNA (LSNA) in anesthetized rats, showing increased RSNA but unchanged LSNA in SHR vs. WKY. THir in skeletal muscle arterioles, NE content in femoral artery, and LSNA were simultaneously reduced by exercise training in the WKY group. Results indicate that THir is a valuable technique to simultaneously evaluate regional patterns of sympathetic activity.

Stability of numerical schemes on staggered grids

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthemore, we consider the cases when both explicit and implicit approximations of the boundary conditions arc employed. Why we choose to do this is clearly motivated and arises front solving fluid flow equations with free surfaces when the Reynolds number can be very small. in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t =n delta t rather than t = (n + 1)delta t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar. thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally, stable. Copyright (C) 2008 John Wiley & Sons, Ltd.

Local influence in estimating equations

Local influence diagnostics based on estimating equations as the role of a gradient vector derived from any fit function are developed for repeated measures regression analysis. Our proposal generalizes tools used in other studies (Cook, 1986: Cadigan and Farrell, 2002), considering herein local influence diagnostics for a statistical model where estimation involves an estimating equation in which all observations are not necessarily independent of each other. Moreover, the measures of local influence are illustrated with some simulated data sets to assess influential observations. Applications using real data are presented. (C) 2010 Elsevier B.V. All rights reserved.

Comparative electrochemical study of 5832-9 stainless steel at different media of biological interest

The electrochemical behavior of ISO 5832-9 stainless steel at 37 degrees C in 0.9% NaCl, Ringer Lactate and minimum essential medium (MEM) has been studied, using linear voltammetry, and surface analysis by SEM and EDS. Mechanical and toxicity tests were made. ISO 5832-9 is passivated at corrosion potential (E) and it does not present pitting corrosion on the media studied from to 50 in V above the transpassivation potential (Ei). SEM and EDS analysis have shown that the sample previously immersed in MEM presents a diffirent behavior at 50 in V above El: the manganese oxide inclusions are absent in the surface. E. values and passivation current density values j(pass) changed according to the following. E(corr, RL) < E(corr,NaCl) < E(corr, MEM) and J (MEM) << j(RL) congruent to j(NaCl) The stainless steel was characterized as non toxic in the cytotoxicity assay