This thesis examines how novel invariants called the "stability index" as proposed by Podvigina and Ashwin can be used to characterize the local geometry of riddled basins of attraction for both skew and non-skew product systems. In particular, it would be interesting to understand how the stability index behaves on the basin boundary between multiple basins of attraction. Then we can ask this question: How can we identify when a basin is riddled? To answer this, we present three models with the presence of riddled basins. In the first model, we present a skew product system of a simple example of a piecewise linear map. We prove that the riddled basin occurs within a certain range of parameter and calculate the stability index analytically for this map. Our results for the stability index at a point show that for Lebesgue almost all points in the map, the index is positive and for some points the index may be negative. We verify these results with our numerical computation for this index. We also make a corollary claiming that the formula for the stability index at a point can be expressed in terms of the stability index for an attractor and Lyapunov exponents for this map. This suggests that this index could be useful as a diagnostic tool to study bifurcation of the riddled basins of attraction. In the second model, we refer to a skew product map studied by Keller. Previously, Keller computed the stability index for an attractor in his map whereas in this thesis, we use an alternative way to compute the index; that is on the basins of attraction for Keller's map, found by inverting his map. Using the same map, we also verify maximum and minimum measures as obtained in his paper by studying Birkhoff averages on periodic points of Markov map in his system. We also conjecture result by Keller and Otani on the dimension of zero sets of invariant graph (i.e. basin boundary) that appears in Keller's map to a complete range of a parameter in the map. The last model is a non-skew product map which is also has a riddled basin. For this map, we compute the stability index for an attractor on the baseline of the map. The result indicates that the index is positive for Lebesgue almost all points whenever the riddled basin occurs.

abstract: According to the CDC in 2010, there were 2.8 million emergency room visits costing $7.9 billion dollars for treatment of nonfatal falling injuries in emergency departments across the country. Falls are a recognized risk factor for unintentional injuries among older adults, accounting for a large proportion of fractures, emergency department visits, and urgent hospitalizations. The objective of this research was to identify and learn more about what factors affect balance using analysis techniques from nonlinear dynamics. Human balance and gait research traditionally uses linear or qualitative tests to assess and describe human motion; however, it is growing more apparent that human motion is neither a simple nor a linear task. In the 1990s Collins, first started applying stochastic processes to analyze human postural control system. Recently, Zakynthinaki et al. modeled human balance using the idea that humans will remain erect when perturbed until some boundary, or physical limit, is passed. This boundary is similar to the notion of basins of attraction in nonlinear dynamics and is referred to as the basin of stability. Human balance data was collected using dual force plates and Vicon marker position data for leans using only ankle movements and leans that were unrestricted. With this dataset, Zakynthinaki’s work was extended by comparing different algorithms used to create the critical curve (basin of stability boundary) that encloses the experimental data points as well as comparing the differences between the two leaning conditions. / Dissertation/Thesis / Masters Thesis Bioengineering 2016

When the states of a system can be described by the extrema of a high-dimensional function, the characterisation of its complexity, i.e. the enumeration of the accessible stable states, can be reduced to a sampling problem. In this thesis a robust numerical protocol is established, capable of producing numerical estimates of the total number of stable states for a broad class of systems, and of computing the a-priori probability of observing any given state. The approach is demonstrated within the context of the computation of the configurational entropy of two and three-dimensional jammed packings. By means of numerical simulation we show the extensivity of the granular entropy as proposed by S.F. Edwards for three-dimensional jammed soft-sphere packings and produce a direct test of the Edwards conjecture for the equivalent two dimensional systems. We find that Edwards’ hypothesis of equiprobability of all jammed states holds only at the (un)jamming density, that is precisely the point of practical significance for many granular systems. Furthermore, two new recipes for the computation of high-dimensional volumes are presented, that improve on the established approach by either providing more statistically robust estimates of the volume or by exploiting the trajectories of the paths of steepest descent. Both methods also produce as a natural by-product unprecedented details on the structures of high-dimensional basins of attraction. Finally, we present a novel Monte Carlo algorithm to tackle problems with fluctuating weight functions. The method is shown to improve accuracy in the computation of the ‘volume’ of high dimensional ‘fluctuating’ basins of attraction and to be able to identify transition states along known reaction coordinates. We argue that the approach can be extended to the optimisation of the experimental conditions for observing certain phenomena, for which individual measurements are stochastic and provide little guidance.

Central to tapping mode atomic force microscopy is an oscillating cantilever whose tip interacts with a sample surface. The tip-surface interactions are strongly nonlinear, rapidly changing, and hysteretic. We explore numerically a lumped-mass model that includes attractive, adhesive, and repulsive contributions as well as the interaction of the capillary fluid layers that cover both tip and sample in the ambient conditions common in experiment. To accomplish this, we have developed and used numerical techniques specifically tailored for discontinuous, nonlinear, and hysteretic dynamical systems. In particular, we use forward-time simulation with event handling and the numerical pseudo-arclength continuation of periodic solutions. We first use these numerical approaches to explore the nonlinear dynamics of the cantilever. We find the coexistence of three steady state oscillating solutions: (i) periodic with low-amplitude, (ii) periodic with high-amplitude, and (iii) high-periodic or irregular behavior. Furthermore, the branches of periodic solutions are found to end precisely where the cantilever comes into grazing contact with event surfaces in state space corresponding to the onset of capillary interactions and the onset of repulsive forces associated with surface contact. Also, the branches of periodic solutions are found to be separated by windows of irregular dynamics. These windows coexist with the periodic branches of solutions and exist beyond the termination of the periodic solution. We also explore the power dissipated through the interaction of the capillary fluid layers. The source of this dissipation is the hysteresis in the conservative capillary force interaction. We relate the power dissipation with the fraction of oscillations that break the fluid meniscus. Using forward-time simulation with event handling, this is done exactly and we explore the dissipated power over a range of experimentally relevant conditions. It is found that the dissipated power as a function of the equilibrium cantilever-surface separation has a characteristic shape that we directly relate to the cantilever dynamics. We also find that despite the highly irregular cantilever dynamics, the fraction of oscillations breaking the meniscus behaves in a fairly simple manner. We have also performed a large number of forward-time simulations over a wide range of initial conditions to approximate the basins of attraction of steady oscillating solutions. Overall, the simulations show a complex pattern of high and low amplitude periodic solutions over the range of initial conditions explored. We find that for large equilibrium separations, the basin of attraction is dominated by the low-amplitude periodic solution and for the small equilibrium separations by the high-amplitude periodic solution. / Ph. D.

Orientador: Michael John Brennan / Resumo: Nesta tese são investigadas, teórica e experimentalmente, a sincronização em fase e a sincronização em anti-fase de dois metrônomos oscilando sobre uma base móvel, a partir de um modelo aqui proposto. Uma descrição do funcionamento do mecanismo de escapamento dos metrônomos é feita, junto a um estudo da relação entre este e o oscilador de van der Pol. Também uma aproximação experimental do valor do amortecimento do metrônomo é fornecida. A frequência instantânea das respostas, numérica e experimental, do sistema é usada na analise. A diferença de outros trabalhos prévios, os dados experimentais têm sido adquiridos usando vídeos dos experimentos e extraídos com ajuda do software Tracker. Para investigar a relação entre as condições iniciais do sistema e seu estado final de sincronização, foram usados mapas bidimensionais chamados ‘basins of attraction’. A relação entre o modelo proposto e um modelo prévio também é mostrada. Encontrou-se que os parâmetros relevantes em relação a ambos os tipos de sincronização são a razão entre a massa do metrônomo e a massa da base, e o amortecimento do sistema. Tem-se encontrado, tanto experimental quanto teoricamente, que a frequência de oscilação dos metrônomos aumenta quando o sistema sincroniza-se em fase, e se mantém a mesma de um metrônomo isolado quando o sistema sincroniza-se em anti-fase. A partir de simulações numéricas encontrou-se que, em geral, incrementos no amortecimento do sistema levam ao sistema se sincronizar mais em fase d... (Resumo completo, clicar acesso eletrônico abaixo) / Doutor

Submitted by CAMILO CARRANZA (carranzacamilo@gmail.com) on 2017-07-25T19:58:22Z No. of bitstreams: 1 Camilo_PhD_Thesis.pdf: 11035322 bytes, checksum: efe400c07b13cabff41e927078789c59 (MD5) / Approved for entry into archive by Luiz Galeffi (luizgaleffi@gmail.com) on 2017-07-26T18:31:30Z (GMT) No. of bitstreams: 1 carranzalopez_jc_dr_ilha.pdf: 11035322 bytes, checksum: efe400c07b13cabff41e927078789c59 (MD5) / Made available in DSpace on 2017-07-26T18:31:30Z (GMT). No. of bitstreams: 1 carranzalopez_jc_dr_ilha.pdf: 11035322 bytes, checksum: efe400c07b13cabff41e927078789c59 (MD5) Previous issue date: 2017-06-05 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Nesta tese são investigadas, teórica e experimentalmente, a sincronização em fase e a sincronização em anti-fase de dois metrônomos oscilando sobre uma base móvel, a partir de um modelo aqui proposto. Uma descrição do funcionamento do mecanismo de escapamento dos metrônomos é feita, junto a um estudo da relação entre este e o oscilador de van der Pol. Também uma aproximação experimental do valor do amortecimento do metrônomo é fornecida. A frequência instantânea das respostas, numérica e experimental, do sistema é usada na analise. A diferença de outros trabalhos prévios, os dados experimentais têm sido adquiridos usando vídeos dos experimentos e extraídos com ajuda do software Tracker. Para investigar a relação entre as condições iniciais do sistema e seu estado final de sincronização, foram usados mapas bidimensionais chamados ‘basins of attraction’. A relação entre o modelo proposto e um modelo prévio também é mostrada. Encontrou-se que os parâmetros relevantes em relação a ambos os tipos de sincronização são a razão entre a massa do metrônomo e a massa da base, e o amortecimento do sistema. Tem-se encontrado, tanto experimental quanto teoricamente, que a frequência de oscilação dos metrônomos aumenta quando o sistema sincroniza-se em fase, e se mantém a mesma de um metrônomo isolado quando o sistema sincroniza-se em anti-fase. A partir de simulações numéricas encontrou-se que, em geral, incrementos no amortecimento do sistema levam ao sistema se sincronizar mais em fase do que em anti-fase. Adicionalmente se encontrou que, para dado valor de amortecimento, diminuir a massa da base leva a uma situação em que a sincronização em anti-fase é mais comum do que a sincronização em fase. / This thesis concerns a theoretical and experimental investigation into the synchronization of two coupled metronomes. A simplified model is proposed to study in-phase and anti-phase synchronization of two metronomes oscillating on a mobile base. A description of the escapement mechanism driving metronomes is given and its relationship with the van der Pol oscillator is discussed. Also an experimental value for the damping in the metronome is determined. The instantaneous frequency of the responses from both numerical and experimental data is used in the analysis. Unlike previous studies, measurements are made using videos and the time domain responses of the metronomes extracted by means of tracker software. Basins of attraction are used to investigate the relationship between initial conditions, parameters and both final synchronization states. The relationship between the model and a previous pendulum model is also shown. The key parameters concerning both kind of synchronization have been found to be the mass ratio between the metronome mass and the base mass, and the damping in the system. It has been shown, both theoretically and experimentally, that the frequency of oscillation of the metronomes increases when the system reaches in-phase synchronization, and is the same as an isolated metronome when the system synchronizes in anti-phase. From numerical simulations, it has been found that, in general, increasing damping leads the system to synchronize more in-phase than in anti-phase. It has also been found that, for a given damping value, decreasing the mass of the base results in the situation where anti-phase synchronization is more common than in-phase synchronization.

Bien que la dynamique des réseaux périodiques non-linéaires ait été investiguée dans les domainestemporel et fréquentiel, il existe un réel besoin d’identifier des relations pratiques avec lephénomène de la localisation d’énergie en termes d’interactions modales et topologies de bifurcation.L’objectif principal de cette thèse consiste à exploiter le phénomène de la localisation pourmodéliser la dynamique collective d’un réseau périodique de résonateurs non-linéaires faiblementcouplés.Un modèle analytico-numérique a été développé pour étudier la dynamique collective d’unréseau périodique d’oscillateurs non-linéaires couplés sous excitations simultanées primaire et paramétrique,où les interactions modales, les topologies de bifurcations et les bassins d’attraction ontété analysés. Des réseaux de pendules et de nano-poutres couplés électrostatiquement ont étéinvestigués sous excitation extérieure et paramétrique, respectivement. Il a été démontré qu’enaugmentant le nombre d’oscillateurs, le nombre de solutions multimodales et la distribution desbassins d’attraction des branches résonantes augmentent. Ce modèle a été étendu pour investiguerla dynamique collective des réseaux 2D de pendules couplés et de billes sphériques en compressionsous excitation à la base, où la dynamique collective est plus riche avec des amplitudes de vibrationplus importantes et des bandes passantes plus larges. Une deuxième investigation de cettethèse consiste à identifier les solitons associés à la dynamique collective d’un réseau périodique etd’étudier sa stabilité. / Although the dynamics of periodic nonlinear lattices was thoroughly investigated in the frequencyand time-space domains, there is a real need to perform profound analysis of the collectivedynamics of such systems in order to identify practical relations with the nonlinear energy localizationphenomenon in terms of modal interactions and bifurcation topologies. The principal goal ofthis thesis consists in exploring the localization phenomenon for modeling the collective dynamicsof periodic arrays of weakly coupled nonlinear resonators.An analytico-numerical model has been developed in order to study the collective dynamics ofa periodic coupled nonlinear oscillators array under simultaneous primary and parametric excitations,where the bifurcation topologies, the modal interactions and the basins of attraction havebeen analyzed. Arrays of coupled pendulums and electrostatically coupled nanobeams under externaland parametric excitations respectively were considered. It is shown that by increasing thenumber of coupled oscillators, the number of multimodal solutions and the distribution of the basinsof attraction of the resonant solutions increase. The model was extended to investigate the collectivedynamics of periodic nonlinear 2D arrays of coupled pendulums and spherical particles underbase excitation, leading to additional features, mainly larger bandwidth and important vibrationalamplitudes. A second investigation of this thesis consists in identifying the solitons associated tothe collective nonlinear dynamics of the considered arrays of periodic structures and the study oftheir stability.