Integrability and non-integrability of periodic non-autonomous Lyness recurrences

This paper studies non-autonomous Lyness type recurrences of the form x_{n+2}=(a_n+x_n)/x_{n+1}, where a_n is a k-periodic sequence of positive numbers with prime period k. We show that for the cases k in {1,2,3,6} the behavior of the sequence x_n is simple(integrable) while for the remaining cases satisfying k not a multiple of 5 this behavior can be much more complicated(chaotic). The cases k multiple of 5 are studied separately.

Integrability and non-integrability of periodic non-autonomous Lyness recurrences (revised and enlarged version)

This paper studies non-autonomous Lyness type recurrences of the form xn+2 = (an+xn+1)=xn, where fang is a k-periodic sequence of positive numbers with primitive period k. We show that for the cases k 2 f1; 2; 3; 6g the behavior of the sequence fxng is simple (integrable) while for the remaining cases satisfying this behavior can be much more complicated (chaotic). We also show that the cases where k is a multiple of 5 present some di erent features.

Algebraic decay of velocity fluctuations near a wall

Computer simulations of the dynamics of a colloidal particle suspended in a fluid confined by an interface show that the asymptotic decay of the velocity correlation functions is algebraic. The exponents of the long-time tails depend on the direction of motion of the particle relative to the surface, as well as on the specific nature of the boundary conditions. In particular, we find that for the angular velocity correlation function, the decay in the presence of a slip surface is faster than the one corresponding to a stick one. An intuitive picture is introduced to explain the various long-time tails, and the simulations are compared with theoretical expressions where available.

Sierpinski curve Julia sets for quadratic rational maps

We investigate under which dynamical conditions the Julia set of a quadratic rational map is a Sierpiński curve.

Observation of a continuous interior crisis in the Hindmarsh-Rose neuron model.

Interior crises are understood as discontinuous changes of the size of a chaotic attractor that occur when an unstable periodic orbit collides with the chaotic attractor. We present here numerical evidence and theoretical reasoning which prove the existence of a chaos-chaos transition in which the change of the attractor size is sudden but continuous. This occurs in the Hindmarsh¿Rose model of a neuron, at the transition point between the bursting and spiking dynamics, which are two different dynamic behaviors that this system is able to present. Moreover, besides the change in attractor size, other significant properties of the system undergoing the transitions do change in a relevant qualitative way. The mechanism for such transition is understood in terms of a simple one-dimensional map whose dynamics undergoes a crossover between two different universal behaviors

Notas del primer seminario de integrabilidad de la Universitat Politècnica de Catalunya

Estas notas corresponden a las exposiciones presentadas en el \emph{Primer Seminario de Integrabilidad}, dentro de lo que se denomina \emph{Aula de Sistemas Din\'amicos}. Durante este evento se realizaron seis conferencias, todas presentadas por miembros del grupo de Sistemas Din\'amicos de la UPC. El programa desarrollado fue el siguiente:\\\begin{center}AULA DE SISTEMAS DIN\'AMICOS\end{center}\begin{center}\texttt{http://www.ma1.upc.es/recerca/seminaris/aulasd-cat.html}\end{center}\begin{center}SEMINARIO DE INTEGRABILIDAD\end{center}\begin{center}Martes 29 y Mi\'ercoles 30 de marzo de 2005\\Facultad de Matem\'aticas y Estad\'{\i}stica, UPC\\Aula: Seminario 1\end{center}\bigskip\begin{center}PROGRAMA Y RES\'UMENES\end{center}{\bf Martes 29 de marzo}\begin{itemize}\item15:30. Juan J. Morales-Ruiz. \emph{El problema de laintegrabilidad en Sistemas Din\'amicos}\medskip {\bf Resumen.} En esta presentaci\'on se pretende dar unaidea de conjunto, pero sin entrar en detalles, sobre las diversasnociones de integrabilidad, asociadas a nombres de matem\'aticostan ilustres como Liouville, Galois-Picard-Vessiot, Lie, Darboux,Kowalevskaya, Painlev\'e, Poincar\'e, Kolchin, Lax, etc. Adem\'astambi\'en mencionaremos la revoluci\'on que supuso en los a\~nossesenta del siglo pasado el descubrimiento de Gardner, Green,Kruskal y Miura sobre un nuevo m\'etodo para resolver en algunoscasos determinadas ecuaciones en derivadas parciales. \medskip\item16:00. David G\'omez-Ullate. \emph{Superintegrabilidad, pares deLax y modelos de $N-$cuerpos en el plano}\medskip{\bf Resumen.} Introduciremos algunas t\'ecnicas cl\'asicas paraconstruir modelos de N-cuerpos integrables, como los pares de Laxo la din\'amica de los ceros de un polinomio. Revisaremos lanoci\'on de integrabilidad Liouville y superintegrabilidad, ydiscutiremos un nuevo m\'etodo debido a F. Calogero para contruirmodelos de N-cuerpos en el plano con muchas \'orbitasperi\'odicas. La exposici\'on se acompa\~nar\'a de animaciones delmovimiento de los cuerpos, y se plantear\'an algunos problemasabiertos.\medskip\item17:00. Pausa\medskip\item17:30. Yuri Fedorov. \emph{An\'alisis de Kovalevskaya--Painlev\'ey Sistemas Algebraicamente Integrables}\medskip{\bf Resumen.} Muchos sistemas integrables poseen una propiedadremarcable: todas sus soluciones son funciones meromorfas deltiempo como una variable compleja. Tal comportamiento, que serefiere como propiedad de Kovalevskaya-Painleve (KP) y que se usafrecuentemente como una ensayo de integrabilidad, no es accidentaly tiene unas ra\'{\i}ces geom\'etricas profundas. En esta charladescribiremos una clase de tales sistemas (conocidos como lossistemas algebraicamente integrables) y subrayaremos suspropiedades geom\'etricas principales que permiten predecir laestructura de las soluciones complejas y adem\'as encontrarlasexpl\'{\i}citamente. Eso lo ilustraremos con algunos sistemas dela mec\'anica cl\'asica. Tambi\'en mencionaremos unasgeneralizaciones \'utiles de la noci\'on de integrabilidadalgebraica y de la propiedad KP.\end{itemize}\medskip{\bf Mi\'ercoles 30 de marzo}\begin{itemize}\item 15:30. Rafael Ram\'{\i}rez-Ros. \emph{El m\'etodo de Poincar\'e}\medskip{\bf Resumen.} Dado un sistema Hamiltoniano aut\'onomo cercano acompletamente integrable Poincar\'e prob\'o que, en general, noexiste ninguna integral primera adicional uniforme en elpar\'ametro de perturbaci\'on salvo el propio Hamiltoniano.Esbozaremos las ideas principales del m\'etodo de prueba ycomentaremos algunas extensiones y generalizaciones.\newpage\item16:30. Chara Pantazi. \emph{El M\'etodo de Darboux}\medskip{\bf Resumen.} Darboux, en 1878, present\'o su m\'etodo paraconstruir integrales primeras de campos vectoriales polinomialesutilizando sus curvas invariantes algebraicas. En estaexposici\'on presentaremos algunas extensiones del m\'etodocl\'asico de Darboux y tambi\'en algunas aplicaciones.\medskip\item17:30. Pausa\medskip\item18:00. Juan J. Morales-Ruiz. \emph{M\'etodos recientes paradetectar la no integrabilidad}\medskip{\bf Resumen.} En 1982 Ziglin utiliza la estructura de laecuaci\'on en variaciones de Poincar\'e (sobre una curva integralparticular) como una herramienta fundamental para detectar la nointegrabilidad de un sistema Hamiltoniano. En esta charla sepretende dar una idea de esta aproximaci\'on a la nointegrabilidad, junto con t\'ecnicas m\'as recientes queinvolucran la teor\'{\i}a de Galois de ecuaciones diferencialeslineales, haciendo \'enfasis en los ejemplos m\'as que en lateor\'{\i}a general. Ilustraremos estos m\'etodos con resultadossobre la no integrabilidad de algunos problemas de $N$ cuerpos enMec\'anica Celeste.\end{itemize}

A parameterization method for the computation of invariant tori and their whiskers in quasi-periodic maps: explorations and mechanisms for the breakdown of hyperbolicity

In two previous papers [J. Differential Equations, 228 (2006), pp. 530 579; Discrete Contin. Dyn. Syst. Ser. B, 6 (2006), pp. 1261 1300] we have developed fast algorithms for the computations of invariant tori in quasi‐periodic systems and developed theorems that assess their accuracy. In this paper, we study the results of implementing these algorithms and study their performance in actual implementations. More importantly, we note that, due to the speed of the algorithms and the theoretical developments about their reliability, we can compute with confidence invariant objects close to the breakdown of their hyperbolicity properties. This allows us to identify a mechanism of loss of hyperbolicity and measure some of its quantitative regularities. We find that some systems lose hyperbolicity because the stable and unstable bundles approach each other but the Lyapunov multipliers remain away from 1. We find empirically that, close to the breakdown, the distances between the invariant bundles and the Lyapunov multipliers which are natural measures of hyperbolicity depend on the parameters, with power laws with universal exponents. We also observe that, even if the rigorous justifications in [J. Differential Equations, 228 (2006), pp. 530-579] are developed only for hyperbolic tori, the algorithms work also for elliptic tori in Hamiltonian systems. We can continue these tori and also compute some bifurcations at resonance which may lead to the existence of hyperbolic tori with nonorientable bundles. We compute manifolds tangent to nonorientable bundles.

Sierpinski curve Julia sets for quadratic rational maps

We investigate under which dynamical conditions the Julia set of a quadratic rational map is a Sierpiński curve.

Persistent random walk model for transport trough thin slabs

We present a model for transport in multiply scattering media based on a three-dimensional generalization of the persistent random walk. The model assumes that photons move along directions that are parallel to the axes. Although this hypothesis is not realistic, it allows us to solve exactly the problem of multiple scattering propagation in a thin slab. Among other quantities, the transmission probability and the mean transmission time can be calculated exactly. Besides being completely solvable, the model could be used as a benchmark for approximation schemes to multiple light scattering.

SQualTrack: A Tool for Robust Fault Detection

One of the techniques used to detect faults in dynamic systems is analytical redundancy. An important difficulty in applying this technique to real systems is dealing with the uncertainties associated with the system itself and with the measurements. In this paper, this uncertainty is taken into account by the use of intervals for the parameters of the model and for the measurements. The method that is proposed in this paper checks the consistency between the system's behavior, obtained from the measurements, and the model's behavior; if they are inconsistent, then there is a fault. The problem of detecting faults is stated as a quantified real constraint satisfaction problem, which can be solved using the modal interval analysis (MIA). MIA is used because it provides powerful tools to extend the calculations over real functions to intervals. To improve the results of the detection of the faults, the simultaneous use of several sliding time windows is proposed. The result of implementing this method is semiqualitative tracking (SQualTrack), a fault-detection tool that is robust in the sense that it does not generate false alarms, i.e., if there are false alarms, they indicate either that the interval model does not represent the system adequately or that the interval measurements do not represent the true values of the variables adequately. SQualTrack is currently being used to detect faults in real processes. Some of these applications using real data have been developed within the European project advanced decision support system for chemical/petrochemical manufacturing processes and are also described in this paper

Chaotic, memory, and cooling rate effects in spin glasses: Evaluation of the Edwards-Anderson model

We investigate chaotic, memory, and cooling rate effects in the three-dimensional Edwards-Anderson model by doing thermoremanent (TRM) and ac susceptibility numerical experiments and making a detailed comparison with laboratory experiments on spin glasses. In contrast to the experiments, the Edwards-Anderson model does not show any trace of reinitialization processes in temperature change experiments (TRM or ac). A detailed comparison with ac relaxation experiments in the presence of dc magnetic field or coupling distribution perturbations reveals that the absence of chaotic effects in the Edwards-Anderson model is a consequence of the presence of strong cooling rate effects. We discuss possible solutions to this discrepancy, in particular the smallness of the time scales reached in numerical experiments, but we also question the validity of the Edwards-Anderson model to reproduce the experimental results.

Robust fault detection using consistency techniques for uncertainty handling

Often practical performance of analytical redundancy for fault detection and diagnosis is decreased by uncertainties prevailing not only in the system model, but also in the measurements. In this paper, the problem of fault detection is stated as a constraint satisfaction problem over continuous domains with a big number of variables and constraints. This problem can be solved using modal interval analysis and consistency techniques. Consistency techniques are then shown to be particularly efficient to check the consistency of the analytical redundancy relations (ARRs), dealing with uncertain measurements and parameters. Through the work presented in this paper, it can be observed that consistency techniques can be used to increase the performance of a robust fault detection tool, which is based on interval arithmetic. The proposed method is illustrated using a nonlinear dynamic model of a hydraulic system

SQualTrack: A Tool for Robust Fault Detection

One of the techniques used to detect faults in dynamic systems is analytical redundancy. An important difficulty in applying this technique to real systems is dealing with the uncertainties associated with the system itself and with the measurements. In this paper, this uncertainty is taken into account by the use of intervals for the parameters of the model and for the measurements. The method that is proposed in this paper checks the consistency between the system's behavior, obtained from the measurements, and the model's behavior; if they are inconsistent, then there is a fault. The problem of detecting faults is stated as a quantified real constraint satisfaction problem, which can be solved using the modal interval analysis (MIA). MIA is used because it provides powerful tools to extend the calculations over real functions to intervals. To improve the results of the detection of the faults, the simultaneous use of several sliding time windows is proposed. The result of implementing this method is semiqualitative tracking (SQualTrack), a fault-detection tool that is robust in the sense that it does not generate false alarms, i.e., if there are false alarms, they indicate either that the interval model does not represent the system adequately or that the interval measurements do not represent the true values of the variables adequately. SQualTrack is currently being used to detect faults in real processes. Some of these applications using real data have been developed within the European project advanced decision support system for chemical/petrochemical manufacturing processes and are also described in this paper

An Interval Intelligent-based Approach for Fault Detection and Modelling

Not considered in the analytical model of the plant, uncertainties always dramatically decrease the performance of the fault detection task in the practice. To cope better with this prevalent problem, in this paper we develop a methodology using Modal Interval Analysis which takes into account those uncertainties in the plant model. A fault detection method is developed based on this model which is quite robust to uncertainty and results in no false alarm. As soon as a fault is detected, an ANFIS model is trained in online to capture the major behavior of the occurred fault which can be used for fault accommodation. The simulation results understandably demonstrate the capability of the proposed method for accomplishing both tasks appropriately

Necessary and sufficient conditions for robust stability using modal intervals

In this paper, robustness of parametric systems is analyzed using a new approach to interval mathematics called Modal Interval Analysis. Modal Intervals are an interval extension that, instead of classic intervals, recovers some of the properties required by a numerical system. Modal Interval Analysis not only simplifies the computation of interval functions but allows semantic interpretation of their results. Necessary, sufficient and, in some cases, necessary and sufficient conditions for robust performance are presented