Title: A traffic flow with a bottelneck Author: Adam Kovařík Department: Department of Numerical Mathematics Supervisor: prof. RNDr. Vladimír Janovský, DrSc. Supervisor's e-mail address: janovsky@karlin.mff.cuni.cz Abstract: In this paper we study a microscopic follow-the-leader traffic model on a circu- lar road with a bottleneck. We assume that all drivers are identical and overtaking is not permitted. We sketch a small part of the rich dynamics of the model including Hopf and Neimark-Sacker bifurcations. We introduce so called POM and quasi-POM solutions and an algorithm how to search them. The main goal of this work is to investigate how the optimal velocity model with a bottleneck deals with so called aggressive behavior of dri- vers. The effect of variable reaction time and a combination of both named factors is also tested. Using numerical simulations we'll find out that aggressiveness and faster reactions have positive effect on traffic flow. In the end we discuss models with two bottlenecks and with one extraordinary driver. Keywords: dynamical systems, ODEs, traffic flow, bottleneck, aggressiveness. 1

L'objectif de cette thèse est d'étudier les comportements typiques d'une famille de transformations du cylindre discret Z×T (où T=R/Z est le cercle de longueur un). Appliquez une rotation à chaque cercle du cylindre puis coupez tous les cercles en deux et déplacez une moitié de chaque cercle d'un niveau vers le bas et une moitié d'un niveau vers le haut. Nous utilisons pour cela des résultats existants en théorie des échanges d'intervalles et en théorie des surfaces de translation compactes. Tout d'abord, nous avons prouvé que pour presque toute suite bi-infinie de rotations, le système obtenu est conservatif (c'est à dire il n'y a pas d'ensemble errant de mesure strictement positive). Ensuite, nous avons prouvé que pour un ensemble Gδ-dense de paramètres, le système est en même temps conservatif, minimal et ergodique. Ce système a un rapport heuristique avec une famille de billards planaires, ainsi qu'une traduction dans des flots sur des surfaces de translation de genre infini. Cela est expliqué dans la thèse. / The main objective of this thesis is the study of the typical dynamical behaviour of a family of transformations on the discrete cylinder Z×T (where T=R/Z is the length one circle). Apply a rotation to every single circle of the cylinder then cut every circle in two and move half of each circle one level down and the other half one level up. To achieve this goal, we use existing results about interval exchange transformations and about compact translation surfaces. First, we proved that for almost every bi-infinite sequence of rotations, the obtained system is conservative (i.e. there is not wandering set of positive measure). Next, we proved that for a Gδ-dense set of parameters, the described system is ergodic, minimal and conservative. This system is heuristically related to a family of planar billiards, it has also a translation into flows on infinite genus translation surfaces.

Investigaremos algumas propriedades dinâmicas e de transporte para um conjunto de partículas clássicas não interagentes em diversos sistemas físicos. Os sistemas descritos aqui, em sua maioria, apresentam estrutura mista no espaço de fase no sentido de que curvas invariantes do tipo spanning, mares de caos e ilhas periódicas estão presentes. A descrição de cada sistema será feita utilizando mapeamentos discretos não lineares. Detalharemos a forma de obter os mapeamentos assim como discutiremos algumas de suas propriedades dinâmicas. Expoentes de Lyapunov serão utilizados para caracterizar a região de caos nos sistemas. Hipóteses de escala são usadas para provar que certos observáveis, por exemplo a energia média ao longo de mares de caos, são invariantes de escala. Consideraremos também que quando uma partícula, ou de forma equivalente um conjunto delas atinge uma determinada altura no espaço de fases, ela pode escapar. Ao estudar o escape de partículas, vemos que o histograma do número de partículas que atingem uma certa altura (ou energia) h no espaço de fases em uma dada iterada n, ao qual observamos ser invariante de escala, cresce rapidamente até atingir um máximo e então tende à zero para n grande. Quando a altura h varia proporcionalmente a posição da primeira curva invariante spanning, podemos confirmar uma invariância de escala do histograma de frequências. O mesmo ocorre para a probabilidade de sobrevivência da partícula à dinâmica. Neste contexto, abordaremos os seguintes problemas: (1) Um guia de ondas senoidalmente corrugado; (2) Uma família de mapas Hamiltonianos bidimensionais que recupera diversos modelos; (3) Partículas confinadas em uma caixa com potenciais infinitos nas bordas e contendo em seu interior um poço de potencial dependente periodicamente do tempo; (4) Analisaremos um bilhar ovóide com dependência temporal introduzida através de giro, onde para certas condições observamos que este não apresenta um aparente crescimento ilimitado de energia (aceleração de Fermi), desta forma sendo um possível contra-exemplo da conjectura LRA. Esta tese é um resumo de 8 artigos que foram publicados em revistas internacionais. / We investigate some dynamical and transport properties for a set of non-interacting classical particles. The systems here described, for the most part, present mixed structure in the phase space in the sense that invariant spanning curves, chaotic seas and periodic islands are present. The dynamics of each model is described by using non-linear mappings. We show all the details to construct the mappings and discuss some of their dynamical properties including fixed points stability among others. Lyapunov exponents will be obtained to characterize the chaotic dynamics observed in the phase space. Moreover some scaling hypotheses are used to prove that certain observables, including the average energy, are scaling invariant. We consider also that when a particle or an ensemble of them reach a certain portion of the phase space, they can escape. When studying the escape, we see that the histogram for the number of particles that reach certain height (or energy) h in the phase space for the iteration n, for which we observe to be scaling invariant, grows quickly until reaching a maximum and then goes towards zero for large enough n. When changing the height h proportionally to the position of the first invariant spanning curve, we can confirm the scaling invariance. The same happens for the survival probability for a particle in the chaotic dynamics. In this way, we will discuss the following problems: (1) A corrugated waveguide; (2) A family of two-dimensional Hamiltonian mappings which can reproduce different scaling exponents; (3) Particles confined to bounce in the interior of a time-dependent potential well; (4) We will analyse a rotating oval billiard, where for certain conditions we observed that this system does not present the unbounded energy growth (Fermi acceleration), in this way it is a possible counterexample of the LRA conjecture. This thesis is as summary of eight papers already published.

A precisão de uma ação motora depende de vários fatores, como: 1) grau de variabilidade da força gerada por cada músculo envolvido, 2) velocidade de geração da força, 3) coordenação das ativações dos músculos. A geração e o controle da força muscular possuem mecanismos que ainda precisam ser mais bem estabelecidos, tanto para o aprimoramento das teorias de controle motor, quanto para o desenvolvimento de técnicas que permitam a prevenção ou a compensação de certas deficiências. A perda de desempenho motor pode ser decorrente de doenças que afetam o sistema neuromuscular ou de alterações associadas ao envelhecimento. Sabe-se, por exemplo, que idosos podem possuir maior variabilidade e menor velocidade de desenvolvimento da força, quando comparados com jovens. Uma das formas de se entender os mecanismos responsáveis pelos fenômenos observados em experimentos neurofisiológicos, em indivíduos saudáveis, em pacientes ou em idosos, é por meio de uma representação adequada de tais mecanismos em modelos matemáticos. Tais modelos podem, pela escolha de um conjunto de parâmetros e de sinais de entrada, ser simulados, explorando-se toda gama de cenários plausíveis para a geração de um determinado fenômeno, tendo como referência os dados obtidos experimentalmente. Resumidamente, o presente trabalho trata do estudo do sistema neuromuscular por modelagem matemática e simulação computacional, com particular interesse nos músculos do tríceps sural e no primeiro interósseo dorsal (um músculo intrínseco da mão), sendo estes músculos amplamente utilizados em estudos experimentais e de modelagem. Maior enfoque é dado em contrações isométricas (i.e., ângulo articular mantido fixo), avaliando-se a organização do núcleo motor, em termos anatômicos e fisiológicos, recebendo como entrada a atividade sináptica das vias pré-motoneuronais, e estudando como diferentes arranjos das propriedades neurais podem resultar em características encontradas experimentalmente para a força muscular. Inicialmente foi feita uma ampla expansão de um simulador existente (ReMoto), tanto em aspectos de modelagem quanto de interface. Em seguida, este modelo expandido foi empregado para um estudo da influência do grau de rigidez muscular nas respostas reflexas do tornozelo. Posteriormente, um novo modelo de pool de motoneurônios, com ampla representação de características biofísicas, foi desenvolvido e, por fim, este modelo foi utilizado, em conjunto com modelos de tratos descendentes e da geração de força muscular, para representar a geração de forças isométricas em jovens e idosos. / The precision of a motor action depends on several factors such as: 1) the level of force variability of each involved muscle, 2) the rate of force development, 3) the coordination of the activations of the muscles involved. Several mechanisms underlying the force generation in a muscle and its control by the nervous system remain to be fully comprehended. An appropriate description of these mechanisms would allow an improvement in motor control theories and could contribute to the development of techniques for the prevention or compensation of some disabilities. Losses in motor performance may be caused by diseases affecting the neuromuscular system or due to aging processes. For instance, old adults may exhibit higher force variability and lower velocity of force development than young adults. Proper representations of such mechanisms in mathematical models constitute a promising way to test hypotheses raised by neurophysiological experiments to explain an observed phenomenon. These models can be used to investigate aspects of health/disease or young/old subjects and, by varying their parameter sets, it is possible to explore a broad range of plausible scenarios under which the experimentally observed phenomena are replicated. This project deals with the study of the neuromuscular system by mathematical modeling and computer simulations, applied to the triceps surae and the first dorsal interosseus (two of the most experimentally and theoretically studied muscles). The principal focus is on isometric contractions (i.e., fixed joint angle) and the study of the organization of the motor nucleus (anatomical and physiological aspects) receiving inputs from premotoneuronal pathways. The study analyzes how different patterns of organization result in experimentally observed aspects of muscle force. Initially, an existing simulator of the neuromuscular system (ReMoto) was broadly extended to include new models and a friendly interface. The extended model was used to investigate the influence of muscle stiffness on the reflex responses in the ankle joint. Next, a new motoneuron pool mathematical model was developed based on known biophysics. Finally, this model was integrated with models of pre-motoneuronal neurons estabilishing synapses with motoneurons and of muscle force generation in order to represent isometric force generation in young and old adults.

The theory of impulsive dynamical systems is an important tool to describe the evolution of systems where the continuous development of a process is interrupted by abrupt changes of state. This phenomenon is called impulse. In many natural phenomena, the real deterministic models are often described by systems which involve impulses. The aim of this work is to investigate topological properties of impulsive semidynamical systems. We establish necessary and sufficient conditions to obtain uniform and orbital stability via Lyapunov functions. We solve a problem of Jake Hale for impulsive systems where we obtain the existence of a maximal compact invariant set. Also, we obtain results about almost periodic motions and asymptotically almost periodic motions in the context of impulsive systems. Some asymptotic properties for impulsive systems and for their associated discrete systems are investigated. The new results presented in this text are in the papers [11], [15] and [16]. / A teoria de sistemas dinâmicos com impulsos é apropriada para descrever processos de evolução que sofrem variações de estado de curta duração e que podem ser consideradas instantâneas. Este fenômeno é chamado impulso. Para muitos fenômenos naturais, os modelos determinísticos mais realistas são frequentemente descritos por sistemas que envolvem impulsos. O objetivo deste trabalho é estudar propriedades topológicas para sistemas semidinâmicos impulsivos. Estabelecemos condições necessárias e suficientes para obtermos estabilidade uniforme e estabilidade orbital utilizando funções do tipo Lyapunov. Resolvemos um problema de Jack Hale para os sistemas impulsivos, onde obtemos a existência de um conjunto invariante compacto maximal. Além disso, obtemos resultados de movimentos quase periódicos e movimentos assintoticamente quase periódicos para sistemas impulsivos. Algumas propriedades assintóticas são estabelecidas para um sistema impulsivo e para seu sistema discreto associado. Os resultados novos apresentados neste trabalho estão presentes nos artigos [11], [15] e [16].

A importância do planejamento de sistemas hidrotérmicos é largamente reconhecida. A questão constitui-se em um complicado problema de otimização que envolve aspectos não-lineares e estocásticos, e que se torna mais complexo quando várias usinas são consideradas em conjunto, em virtude da ocorrência de interações. Apesar de se conhecerem formalizações teóricas do problema e da solução, a complexidade acarreta sérios entraves computacionais e muitas vezes consideram-se apenas cenários simplificados, como o determinístico ou estocástico com um modelo de reservatório equivalente. Este trabalho tem como objetivo conciliar a facilidade numérica de problemas determinísticos com a consideração indireta da natureza estocástica do problema. Com esta finalidade, propõe-se o uso de um funcional de custo que reflete o risco de vertimento a cada período. Outra linha estudada neste trabalho foi a de rastreamento de alvo utilizando a técnica de controle Linear Quadrático (LQ). A eficácia das estratégias propostas é avaliada através de diversos casos de estudo, incluindo comparação com o planejamento determinístico usual / The scheduling of hydrothermal systems is an important and challenging problem. In fact, the problem is complex since it involves the optimization of a cost functional related to a nonlinear and stochastic system. Although mathematical formulations of the problem are available, there are computational issues (as the course of dimensionality in dynamical programming) that make difficult to obtain the solution. In this situation, some simplifyied equivalent models or deterministic formulations were considered in the literature, which are simpler to deal with. The aim of this work is to obtain a solution thet presents the ease of computaion of deterministic models and, simultaneously, it reflects the stochastic nature of the problem, at some extension. The strategy consists of taking into account a cost functional that is related to the risk of spillage at each period. Another strategy considered here is the linear quadratic control with a target. Several case estudies are presented, which employ different models for the stochastic inflow, allwing us to evaluate the effectiveness of the strategies and to compare them with the usual deterministic control

O filtro de Kalman é amplamente conhecido e utilizado em aplicações, em virtude de apresentar diversas propriedades interessantes. Este trabalho aborda uma das características mais importantes, a estabilidade do filtro de Kalman aplicado a sistemas lineares discretos com saltos Markovianos. Sistemas desta classe são muito empregados em problemas práticos. Neste trabalho mostramos que o conceito de controlabilidade fraca e detetabilidade estocástica são condições suficientes para estabilidade do filtro de Kalman com relação a condição inicial. No que se refere a estabilidade no sentido mais usual, apresentamos resultados parciais, dependentes de uma condição adicional sobre a cadeia de Markov, bem como uma conjectura. O estudo da estabilidade do filtro de Kalman é relevante, pois filtros instáveis oferecem estimativas de baixa qualidade. O tema tem interesse teórico inerente e é bastante relevante para aplicações.O filtro de Kalman é amplamente conhecido e utilizado em aplicações, em virtude de apresentar diversas propriedades interessantes. Este trabalho aborda uma das características mais importantes, a estabilidade do filtro de Kalman aplicado a sistemas lineares discretos com saltos Markovianos. Sistemas desta classe são muito empregados em problemas práticos. Neste trabalho mostramos que o conceito de controlabilidade fraca e detetabilidade estocástica são condições suficientes para estabilidade do filtro de Kalman com relação a condição inicial. No que se refere a estabilidade no sentido mais usual, apresentamos resultados parciais, dependentes de uma condição adicional sobre a cadeia de Markov, bem como uma conjectura. O estudo da estabilidade do filtro de Kalman é relevante, pois filtros instáveis oferecem estimativas de baixa qualidade. O tema tem interesse teórico inerente e é bastante relevante para aplicações / Kalman filters present several interesting features that make them relevant for many applications. In this work we study one of the main issues in Kalman filtering - stability. We deal with Kalman filters for Markov jump linear systems, a class of systems with applications in many different areas. We consider the concepts of weak controllability and stochastic detectability and we show that they ensure stability of the Kalman filter with respect to the initial condition. As for the stability, we present some results relying in a conjecture and an additional condition on the Markov chain. The study of the stability of the Kalman filter is important, since unstable filters may lead to poor estimates. The stability issue has inherent theoretical interest and is relevant for applications

Processos de sincronização são observados em uma imensa quantidade de sistemas físicos, biológicos, químicos, tecnológicos e sociais. Tais sistemas podem ser descritos e modelados utilizando a teoria das redes complexas, de forma que o completo entendimento da emergência do comportamento coletivo nestes sistemas complexos só é alcançado por teorias que englobam a interação entre seus elementos. Nesta dissertação, estudamos a emergência de transições de fase de primeira-ordem na sincronização de osciladores acoplados através de estruturas heterogêneas e não-triviais. Utilizando teorias de campo médio, obtemos a expressão analítica do acoplamento crítico necessário para a ocorrência de sincronização explosiva em redes livre-escala. Além disso, estudamos o comportamento de tais transições na presença de atrasos temporais e verificamos que é possível elevar o grau de sincronismo dos osciladores quando a interação se dá de forma não-instantânea. Os resultados obtidos contribuem para um melhor entendimento da relação entre topologia e dinâmica em redes. / Synchronization processes are observed in many physical, biological, chemical, technological and social systems. These systems can be described and modelled through the theory of complex networks, in a way that the full comprehension of the emergence of collective behavior in these complex systems will only be achieved by theories that encompass the interaction of its elements. In this thesis, we study the emergence of first-order phase transitions in the synchronization of oscillators coupled through heterogeneous and non-trivial structures. By using mean-field theories, we obtain an analytical expression for the critical coupling necessary for the occurrence of explosive synchronization in scale-free networks. Furthermore, we study the behavior of such transitions in the presence of time delays, verifying that is possible to enhance the synchronization level of the oscillators when the interaction is non-instantaneous. The obtained results contribute for the better understanding of the interplay between topology and dynamics in networks.

O estudo dos processos quebra espontânea de simetria na natureza têm atraído interesse em diversas áreas da física, como por exemplo em física quântica no estudo das energias de um átomo tal como em física de altas energias, no estudo das partículas elementares. Esses processos até então envolviam sistemas físicos microscópicos, em 1993 surge uma proposta de agregar as idéias de quebra de simetria à um sistema macroscópico, o código genético. A idéia básica é que os códons que formam código se diferenciam em um processo de quebra de simetria, preservando suas propriedades de degenerescência, nos dando uma \"picture\" de como se fez essa diferenciação que resultam nos 20 aminoácidos e do sinal de terminação que se conhece nos dias atuais. Esse modelo nos diz por exemplo, quantos eram os aminoácidos primordiais. Nosso interesse está na verificação dessa quebra de simetria e estudar as relações entre os códons do ponto de vista temporal para tanto, usamos aqui um sistema dinâmico. Esse sistema conserva no princípio de sua evolução a simetria proposta pelo modelo e através de um processo de quebra de simetria estudaremos se esse processo reproduz a cadeia de quebra de simetria proposta no modelo. Como primeiro passo estudamos a representação tridimensional do grupo Sp( 6), que serve como ponto de partida no processo de quebra de simetria no modelo, essa representação é conhecida como grupo de Weyl do Sp(6). É possível construir um sistema dinâmico ou mapa, que na verdade é uma função do R3, com as mesmas propriedades de simetria do grupo de Weyl do Sp(6). A construção desse sistema e seu estudo matemático acarreta no segundo passo deste trabalho. O mapa construído depende de parâmetros que variados de forma correta produzem uma cadeia de quebra de simetria. O estudo dessa quebras consiste no terceiro passo deste trabalho. Por fim determinamos a ação ou seja, como esse sistema muda a rotulação dos códons anteriormente proposta no modelo e mais ainda, que informação biológica poderá ser extraída desse sistema. Como resultado obtivemos em grande parte a ratificação do modelo proposto mostrando que a quebra proposta e a rotulação dos códons de acordo com a ordem evolutiva dada pela quebra de simetria segue também uma coerência dinâmica. / The study of natural symmetry breaking processes have attracted interest in many physics areas including energy atoms studies in quantum physics as well elementary particles in high energy physics. These processes were related with microscpics physic systems, in 1993 appears one propose to use the ideas of symmetry breaking in one macroscopic system the genetic code. The basic idea is that the differentiation of the codons wich are components of the code was done in a process of symmetry breaking preserving the degeneracy properties given to us one picture of how this process occour resulting in 20 aminoacids and termination sign known in the present days. With this model for example, we can predict how many aminoacids were primordial\'s. Our interest is in verify this symmetry breaking and study the codon temporal relations for this we use a dynamical system. The preservation of the starting symmetry proposed by the model is the main caracteristic of our system and through of the symmetry breaking we will study what relations between the symmetry breaking proposed by the model and the dynamical symmetry breaking. As a first step we will study the group Sp(6) in its tridimensional representation which is the starting point in the symmetry breaking process in the mode!. This representation in known as Weyl group of Sp(6). It is possible construct this dynamical system or map, which is one function in the R3 with the same symmetry properties of the Weyl group of Sp(6). The construction of this map and its mathematical study is our second step of this work. The map depends on parameter\'s which are changed in a correct way to produce some symmetry breaking chain. The symmetry breaking studies is our third step. At the end we look at the action of our map in the codons, in other words, how this action change the codons labelling proposed by the model Moreover, what kind of biological information can be extract from this action. As a result the symmetry breaking and the labelling of codons proposed by the model are isomophics, with little restrictions, when compared with the dynamical systems.

Sincronismo é um comportamento dinâmico que pode ser detectado na atividade dos neurônios que formam estruturas corticais, e essa atividade parece ser fundamental para a realização de processos relacionados a aprendizado, compreensão e reconhecimento. Estruturas corticais podem ser modeladas por redes de mapas acoplados. Nessas redes, a existência de solução síncrona é determinada pela equação que governa a atividade de cada mapa, pela topologia de acoplamento entre os mapas e pelos valores dos parâmetros. O enfoque deste trabalho é empregar redes de mapas acoplados voltados à detecção de simetria, e segmentação de imagens, via sincronismo. / Synchronism is a dynamical behavior that can be found in the neural activity of cortical structures, and this behavior seems to be fundamental in learning, understanding and recognition processes. Cortical structures can be modeled by using coupled map lattices. In such lattices, the existence of synchronous solution is determined by the equation that rules the activity of each map, by the coupling topology among the maps, and by the parameter values. The goal of this work is to employ coupled map lattices in symmetry detection and image segmentation, via synchronism.