Modelagem e simulação da propagação de ondas em barras não homogêneas envolvendo materiais elásticos não lineares.

O objetivo deste trabalho é tratar da simulação do fenômeno de propagação de ondas em uma haste heterogênea elástico, composta por dois materiais distintos (um linear e um não-linear), cada um deles com a sua própria velocidade de propagação da onda. Na interface entre estes materiais existe uma descontinuidade, um choque estacionário, devido ao salto das propriedades físicas. Empregando uma abordagem na configuração de referência, um sistema não-linear hiperbólico de equações diferenciais parciais, cujas incógnitas são a velocidade e a deformação, descrevendo a resposta dinâmica da haste heterogénea. A solução analítica completa do problema de Riemann associado são apresentados e discutidos.

Bessel functions and equations of mathematical physics

The aim of this dissertation is to introduce Bessel functions to the reader, as well as studying some of their properties. Moreover, the final goal of this document is to present the most well- known applications of Bessel functions in physics.

Distribution theory and fundamental solutions of differential operators

The objective of this dissertation is to study the theory of distributions and some of its applications. Certain concepts which we would include in the theory of distributions nowadays have been widely used in several fields of mathematics and physics. It was Dirac who first introduced the delta function as we know it, in an attempt to keep a convenient notation in his works in quantum mechanics. Their work contributed to open a new path in mathematics, as new objects, similar to functions but not of their same nature, were being used systematically. Distributions are believed to have been first formally introduced by the Soviet mathematician Sergei Sobolev and by Laurent Schwartz. The aim of this project is to show how distribution theory can be used to obtain what we call fundamental solutions of partial differential equations.

Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

Building flexible user interfaces for solving PDEs

FEniCS is a collection of software tools for the automated solution of differential equations by finite element methods. In this note, we describe how FEniCS can be used to solve a simple nonlinear model problem with varying levels of automation. At one extreme, FEniCS provides tools for the fully automated and adaptive solution of nonlinear partial differential equations. At the other extreme, FEniCS provides a range of tools that allow the computational scientist to experiment with novel solution algorithms. © 2010 American Institute of Physics.

Reconstruction of arbitrary biochemical reaction networks: A compressive sensing approach

Reconstruction of biochemical reaction networks (BRN) and genetic regulatory networks (GRN) in particular is a central topic in systems biology which raises crucial theoretical challenges in system identification. Nonlinear Ordinary Differential Equations (ODEs) that involve polynomial and rational functions are typically used to model biochemical reaction networks. Such nonlinear models make the problem of determining the connectivity of biochemical networks from time-series experimental data quite difficult. In this paper, we present a network reconstruction algorithm that can deal with ODE model descriptions containing polynomial and rational functions. Rather than identifying the parameters of linear or nonlinear ODEs characterised by pre-defined equation structures, our methodology allows us to determine the nonlinear ODEs structure together with their associated parameters. To solve the network reconstruction problem, we cast it as a compressive sensing (CS) problem and use sparse Bayesian learning (SBL) algorithms as a computationally efficient and robust way to obtain its solution. © 2012 IEEE.

Delayed control of a Moore-Greitzer axial compressor model

Several feedback control laws have appeared in the literature concerning the stabilization of the nonlinear Moore-Greitzer axial compression model. Motivated by magnitude and rate limitations imposed by the physical implementation of the control law, Larsen et al. studied a dynamic implementation of the S-controller suggested by Sepulchre and Kokotović. They showed the potential benefit of implementing the S-controller through a first-order lag: while the location of the closed-loop equilibrium achieved with the static control law was sensitive to poorly known parameters, the dynamic implementation resulted in a small limit cycle at a very desirable location, insensitive to parameter variations. In this paper, we investigate the more general case when the control is applied with a time delay. This can be seen as an extension of the model with a first-order lag. The delay can either be a result of system constraints or be deliberately implemented to achieve better system behavior. The resulting closed-loop system is a set of parameter-dependent delay differential equations. Numerical bifurcation analysis is used to study this model and investigate whether the positive results obtained for the first-order model persist, even for larger values of the delay.

Spontaneous emission from radiative chiral nematic liquid crystals at the photonic band-gap edge: An investigation into the role of the density of photon states near resonance

In this article, we investigate the spontaneous emission properties of radiating molecules embedded in a chiral nematic liquid crystal, under the assumption that the electronic transition frequency is close to the photonic edge mode of the structure, i.e., at resonance. We take into account the transition broadening and the decay of electromagnetic field modes supported by the so-called "mirrorless"cavity. We employ the Jaynes-Cummings Hamiltonian to describe the electron interaction with the electromagnetic field, focusing on the mode with the diffracting polarization in the chiral nematic layer. As known in these structures, the density of photon states, calculated via the Wigner method, has distinct peaks on either side of the photonic band gap, which manifests itself as a considerable modification of the emission spectrum. We demonstrate that, near resonance, there are notable differences between the behavior of the density of states and the spontaneous emission profile of these structures. In addition, we examine in some detail the case of the logarithmic peak exhibited in the density of states in two-dimensional photonic structures and obtain analytic relations for the Lamb shift and the broadening of the atomic transition in the emission spectrum. The dynamical behavior of the atom-field system is described by a system of two first-order differential equations, solved using the Green's-function method and the Fourier transform. The emission spectra are then calculated and compared with experimental data. © 2013 American Physical Society.

On the application of differential geometry to MDO

Multidisciplinary Design Optimization (MDO) is a methodology for optimizing large coupled systems. Over the years, a number of different MDO decomposition strategies, known as architectures, have been developed, and various pieces of analytical work have been done on MDO and its architectures. However, MDO lacks an overarching paradigm which would unify the field and promote cumulative research. In this paper, we propose a differential geometry framework as such a paradigm: Differential geometry comes with its own set of analysis tools and a long history of use in theoretical physics. We begin by outlining some of the mathematics behind differential geometry and then translate MDO into that framework. This initial work gives new tools and techniques for studying MDO and its architectures while producing a naturally arising measure of design coupling. The framework also suggests several new areas for exploration into and analysis of MDO systems. At this point, analogies with particle dynamics and systems of differential equations look particularly promising for both the wealth of extant background theory that they have and the potential predictive and evaluative power that they hold. © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

A differential lyapunov framework for contraction analysis

Lyapunov's second theorem is an essential tool for stability analysis of differential equations. The paper provides an analog theorem for incremental stability analysis by lifting the Lyapunov function to the tangent bundle. The Lyapunov function endows the state-space with a Finsler structure. Incremental stability is inferred from infinitesimal contraction of the Finsler metrics through integration along solutions curves. © 2013 IEEE.

Solitary magnon excitations in a one-dimensional antiferromagnet with Dzyaloshinsky-Moriya interactions

We investigate solitary excitations in a model of a one-dimensional antiferromagnet including a single-ion anisotropy and a Dzyaloshinsky-Moriya antisymmetric exchange interaction term. We employ the Holstein-Primakoff transformation, the coherent state ansatz and the time variational principle. We obtain two partial differential equations of motion by using the method of multiple scales and applying perturbation theory. By so doing, we show that the motion of the coherent amplitude must satisfy the nonlinear Schrodinger equation. We give the single-soliton solution.

A semi-analytical theory of coastal upwelling and jet

The general forms of the conservation of momentum, temperature and potential vorticity of coastal ocean are obtained in the x-z plane for the nonlinear ocean circulation of Boussinesq fluid, and a elliptic type partial differential equations of second order are derived. Solution of the partial differential equations are obtained under the conditions that the fluid moves along the topography. The numerical results show that there exist both upwelling and downwelling along coastline that mainly depends on the large scale ocean condition. Numerically results of the upwelling (downwelling), coastal jet and temperature front zone are favorable to the observations.

Algorithmic derivation of isochronicity conditions

Hill, Joe M., Lloyd, Noel G., Pearson, Jane M., 'Algorithmic derivation of isochronicity conditions', Nonlinear Analysis (2007) 67, 52-69.

SOVEREIGN: A Self-Organizing, Vision, Expectation, Recognition, Emotion, Intelligent, Goal-Oriented Navigation System

Both animals and mobile robots, or animats, need adaptive control systems to guide their movements through a novel environment. Such control systems need reactive mechanisms for exploration, and learned plans to efficiently reach goal objects once the environment is familiar. How reactive and planned behaviors interact together in real time, and arc released at the appropriate times, during autonomous navigation remains a major unsolved problern. This work presents an end-to-end model to address this problem, named SOVEREIGN: A Self-Organizing, Vision, Expectation, Recognition, Emotion, Intelligent, Goal-oriented Navigation system. The model comprises several interacting subsystems, governed by systems of nonlinear differential equations. As the animat explores the environment, a vision module processes visual inputs using networks that arc sensitive to visual form and motion. Targets processed within the visual form system arc categorized by real-time incremental learning. Simultaneously, visual target position is computed with respect to the animat's body. Estimates of target position activate a motor system to initiate approach movements toward the target. Motion cues from animat locomotion can elicit orienting head or camera movements to bring a never target into view. Approach and orienting movements arc alternately performed during animat navigation. Cumulative estimates of each movement, based on both visual and proprioceptive cues, arc stored within a motor working memory. Sensory cues are stored in a parallel sensory working memory. These working memories trigger learning of sensory and motor sequence chunks, which together control planned movements. Effective chunk combinations arc selectively enhanced via reinforcement learning when the animat is rewarded. The planning chunks effect a gradual transition from reactive to planned behavior. The model can read-out different motor sequences under different motivational states and learns more efficient paths to rewarded goals as exploration proceeds. Several volitional signals automatically gate the interactions between model subsystems at appropriate times. A 3-D visual simulation environment reproduces the animat's sensory experiences as it moves through a simplified spatial environment. The SOVEREIGN model exhibits robust goal-oriented learning of sequential motor behaviors. Its biomimctic structure explicates a number of brain processes which are involved in spatial navigation.

A Neural Network Model for the Spatial and Temporal Response of Retinal Ganglion Cells

This article introduces a quantitative model of early visual system function. The model is formulated to unify analyses of spatial and temporal information processing by the nervous system. Functional constraints of the model suggest mechanisms analogous to photoreceptors, bipolar cells, and retinal ganglion cells, which can be formally represented with first order differential equations. Preliminary numerical simulations and analytical results show that the same formal mechanisms can explain the behavior of both X (linear) and Y (nonlinear) retinal ganglion cell classes by simple changes in the relative width of the receptive field (RF) center and surround mechanisms. Specifically, an increase in the width of the RF center results in a change from X-like to Y-like response, in agreement with anatomical data on the relationship between α- and