Evans functions for integral neural field equations with Heaviside firing rate function

In this paper we show how to construct the Evans function for traveling wave solutions of integral neural field equations when the firing rate function is a Heaviside. This allows a discussion of wave stability and bifurcation as a function of system parameters, including the speed and strength of synaptic coupling and the speed of axonal signals. The theory is illustrated with the construction and stability analysis of front solutions to a scalar neural field model and a limiting case is shown to recover recent results of L. Zhang [On stability of traveling wave solutions in synaptically coupled neuronal networks, Differential and Integral Equations, 16, (2003), pp.513-536.]. Traveling fronts and pulses are considered in more general models possessing either a linear or piecewise constant recovery variable. We establish the stability of coexisting traveling fronts beyond a front bifurcation and consider parameter regimes that support two stable traveling fronts of different speed. Such fronts may be connected and depending on their relative speed the resulting region of activity can widen or contract. The conditions for the contracting case to lead to a pulse solution are established. The stability of pulses is obtained for a variety of examples, in each case confirming a previously conjectured stability result. Finally we show how this theory may be used to describe the dynamic instability of a standing pulse that arises in a model with slow recovery. Numerical simulations show that such an instability can lead to the shedding of a pair of traveling pulses.

Symmetric Interior Penalty DG Methods for the Compressible Navier-Stokes Equations I: Method Formulation

In this article we consider the development of discontinuous Galerkin finite element methods for the numerical approximation of the compressible Navier-Stokes equations. For the discretization of the leading order terms, we propose employing the generalization of the symmetric version of the interior penalty method, originally developed for the numerical approximation of linear self-adjoint second-order elliptic partial differential equations. In order to solve the resulting system of nonlinear equations, we exploit a (damped) Newton-GMRES algorithm. Numerical experiments demonstrating the practical performance of the proposed discontinuous Galerkin method with higher-order polynomials are presented.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular we are able to treat "patchy'" connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a "lattice-directed" traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs. Article published and (c) American Physical Society 2007

Finite time extinction for nonlinear fractional evolution equations and related properties

The finite time extinction phenomenon (the solution reaches an equilibrium after a finite time) is peculiar to certain nonlinear problems whose solutions exhibit an asymptotic behavior entirely different from the typical behavior of solutions associated to linear problems. The main goal of this work is twofold. Firstly, we extend some of the results known in the literature to the case in which the ordinary time derivative is considered jointly with a fractional time differentiation. Secondly, we consider the limit case when only the fractional derivative remains. The latter is the most extraordinary case, since we prove that the finite time extinction phenomenon still appears, even with a non-smooth profile near the extinction time. Some concrete examples of quasi-linear partial differential operators are proposed. Our results can also be applied in the framework of suitable nonlinear Volterra integro-differential equations.

Finite time extinction for nonlinear fractional evolution equations and related properties.

The finite time extinction phenomenon (the solution reaches an equilibrium after a finite time) is peculiar to certain nonlinear problems whose solutions exhibit an asymptotic behavior entirely different from the typical behavior of solutions associated to linear problems. The main goal of this work is twofold. Firstly, we extend some of the results known in the literature to the case in which the ordinary time derivative is considered jointly with a fractional time differentiation. Secondly, we consider the limit case when only the fractional derivative remains. The latter is the most extraordinary case, since we prove that the finite time extinction phenomenon still appears, even with a non-smooth profile near the extinction time. Some concrete examples of quasi-linear partial differential operators are proposed. Our results can also be applied in the framework of suitable nonlinear Volterra integro-differential equations.

Instabilities in threshold-diffusion equations with delay

The introduction of delays into ordinary or partial differential equation models is well known to facilitate the production of rich dynamics ranging from periodic solutions through to spatio-temporal chaos. In this paper we consider a class of scalar partial differential equations with a delayed threshold nonlinearity which admits exact solutions for equilibria, periodic orbits and travelling waves. Importantly we show how the spectra of periodic and travelling wave solutions can be determined in terms of the zeros of a complex analytic function. Using this as a computational tool to determine stability we show that delays can have very different effects on threshold systems with negative as opposed to positive feedback. Direct numerical simulations are used to confirm our bifurcation analysis, and to probe some of the rich behaviour possible for mixed feedback.

PRACTICAL RESOLUTION OF DIFFERENTIAL SYSTEMS

3. PRACTICAL RESOLUTION OF DIFFERENTIAL SYSTEMS by Marilia Pires, University of Évora, Portugal This practice presents the main features of a free software to solve mathematical equations derived from concrete problems: i.- Presentation of Scilab (or python) ii.- Basics (number, characters, function) iii.- Graphics iv.- Linear and nonlinear systems v.- Differential equations

Systemic Granulomatous Diseases Associated With Multiple Palpable Masses That May Involve The Breast: Case Presentation And An Approach To The Differential Diagnosis.

Palpable mass is a common complaint presented to the breast surgeon. It is very uncommon for patients to report breast mass associated with palpable masses in other superficial structures. When these masses are related to systemic granulomatous diseases, the diagnosis and initiation of specific therapy can be challenging. The purpose of this paper is to report a case initially assessed by the breast surgeon and ultimately diagnosed as granulomatous variant of T-cell lymphoma, and discuss the main systemic granulomatous diseases associated with palpable masses involving the breast.

Water Stress Reveals Differential Antioxidant Responses Of Tolerant And Non-tolerant Sugarcane Genotypes.

The biochemical responses of the enzymatic antioxidant system of a drought-tolerant cultivar (IACSP 94-2094) and a commercial cultivar in Brazil (IACSP 95-5000) grown under two levels of soil water restriction (70% and 30% Soil Available Water Content) were investigated. IACSP 94-2094 exhibited one additional active superoxide dismutase (Cu/Zn-SOD VI) isoenzyme in comparison to IACSP 95-5000, possibly contributing to the heightened response of IACSP 94-2094 to the induced stress. The total glutathione reductase (GR) activity increased substantially in IACSP 94-2094 under conditions of severe water stress; however, the appearance of a new GR isoenzyme and the disappearance of another isoenzyme were found not to be related to the stress response because the cultivars from both treatment groups (control and water restrictions) exhibited identical changes. Catalase (CAT) activity seems to have a more direct role in H2O2 detoxification under water stress condition and the shift in isoenzymes in the tolerant cultivar might have contributed to this response, which may be dependent upon the location where the excessive H2O2 is being produced under stress. The improved performance of IACSP 94-2094 under drought stress was associated with a more efficient antioxidant system response, particularly under conditions of mild stress.

Enhanced D1 And D2 Inhibitions Induced By Low-frequency Trains Of Conditioning Stimuli: Differential Effects On H- And T-reflexes And Possible Mechanisms.

Mechanically evoked reflexes have been postulated to be less sensitive to presynaptic inhibition (PSI) than the H-reflex. This has implications on investigations of spinal cord neurophysiology that are based on the T-reflex. Preceding studies have shown an enhanced effect of PSI on the H-reflex when a train of ~10 conditioning stimuli at 1 Hz was applied to the nerve of the antagonist muscle. The main questions to be addressed in the present study are if indeed T-reflexes are less sensitive to PSI and whether (and to what extent and by what possible mechanisms) the effect of low frequency conditioning, found previously for the H-reflex, can be reproduced on T-reflexes from the soleus muscle. We explored two different conditioning-to-test (C-T) intervals: 15 and 100 ms (corresponding to D1 and D2 inhibitions, respectively). Test stimuli consisted of either electrical pulses applied to the posterior tibial nerve to elicit H-reflexes or mechanical percussion to the Achilles tendon to elicit T-reflexes. The 1 Hz train of conditioning electrical stimuli delivered to the common peroneal nerve induced a stronger effect of PSI as compared to a single conditioning pulse, for both reflexes (T and H), regardless of C-T-intervals. Moreover, the conditioning train of pulses (with respect to a single conditioning pulse) was proportionally more effective for T-reflexes as compared to H-reflexes (irrespective of the C-T interval), which might be associated with the differential contingent of Ia afferents activated by mechanical and electrical test stimuli. A conceivable explanation for the enhanced PSI effect in response to a train of stimuli is the occurrence of homosynaptic depression at synapses on inhibitory interneurons interposed within the PSI pathway. The present results add to the discussion of the sensitivity of the stretch reflex pathway to PSI and its functional role.

Fundamentos geometricos e algebricos da calibração e reconstrução tridimensional : aplicação na analise cinematica de movimentos humanos

Universidade Estadual de Campinas . Faculdade de Educação Física

Wavelets And Adaptive Grids For The Discontinuous Galerkin Method

In this paper, space adaptivity is introduced to control the error in the numerical solution of hyperbolic systems of conservation laws. The reference numerical scheme is a new version of the discontinuous Galerkin method, which uses an implicit diffusive term in the direction of the streamlines, for stability purposes. The decision whether to refine or to unrefine the grid in a certain location is taken according to the magnitude of wavelet coefficients, which are indicators of local smoothness of the numerical solution. Numerical solutions of the nonlinear Euler equations illustrate the efficiency of the method. © Springer 2005.

Ressonância não linear de uma bússola em campos magnéticos

Fenômenos oscilatórios e ressonantes são explorados em vários cursos experimentais de física. Em geral os experimentos são interpretados no limite de pequenas oscilações e campos uniformes. Neste artigo descrevemos um experimento de baixo custo para o estudo da ressonância em campo magnético da agulha de uma bússola fora dos limites acima. Nesse caso, termos não lineares na equação diferencial são responsáveis por fenômenos interessantes de serem explorados em laboratórios didáticos.

Exact time-dependent solutions for a self-regulating gene

The exact time-dependent solution for the stochastic equations governing the behavior of a binary self-regulating gene is presented. Using the generating function technique to rephrase the master equations in terms of partial differential equations, we show that the model is totally integrable and the analytical solutions are the celebrated confluent Heun functions. Self-regulation plays a major role in the control of gene expression, and it is remarkable that such a microscopic model is completely integrable in terms of well-known complex functions.

Combining Legendre's Polynomials and Genetic Algorithm in the Solution of Nonlinear Initial-Value Problems

This work develops a method for solving ordinary differential equations, that is, initial-value problems, with solutions approximated by using Legendre's polynomials. An iterative procedure for the adjustment of the polynomial coefficients is developed, based on the genetic algorithm. This procedure is applied to several examples providing comparisons between its results and the best polynomial fitting when numerical solutions by the traditional Runge-Kutta or Adams methods are available. The resulting algorithm provides reliable solutions even if the numerical solutions are not available, that is, when the mass matrix is singular or the equation produces unstable running processes.