In many applications from industry and technology computer simulations are performed using models which can be formulated by systems of differential equations. Often the equations underlie additional algebraic constraints. In this context we speak of descriptor systems. Very important characteristic values of such systems are the $L_\infty$-norms of the corresponding transfer functions. The main goal of this thesis is to extend a numerical method for the computation of the $L_\infty$-norm for standard state space systems to descriptor systems. For this purpose we develop a numerical method to check whether the transfer function of a given descriptor system is proper or improper and additionally use this method to reduce the order of the system to decrease the costs of the $L_\infty$-norm computation. When computing the $L_\infty$-norm it is necessary to compute the eigenvalues of certain skew-Hamiltonian/Hamiltonian matrix pencils composed by the system matrices. We show how we extend these matrix pencils to skew-Hamiltonian/Hamiltonian matrix pencils of larger dimension to get more reliable and accurate results. We also consider discrete-time systems, apply the extension strategy to the arising symplectic matrix pencils and transform these to more convenient structures in order to apply structure-exploiting eigenvalue solvers to them. We also investige a new structure-preserving method for the computation of the eigenvalues of skew-Hamiltonian/Hamiltonian matrix pencils and use this to increase the accuracy of the computed eigenvalues even more. In particular we ensure the reliability of the $L_\infty$-norm algorithm by this new eigenvalue solver. Finally we describe the implementation of the algorithms in Fortran and test them using two real-world examples.

The boom systems of mobile cranes and aerial platform vehicles are driven by hydraulic systems, to be specified, valve-controlled hydraulic cylinders. This hydraulic actuated boom system can accomplish the tasks such as lifting heavy loads or carrying personal to high position, by the design of a long boom structure. In practice, the boom structure is designed as light and slender as possible to control the structure self-weight. However, such structure is quite flexible and can be easily stimulated by the loads, including the driving force or torque from the hydraulic system. Our research focuses on trajectory planning for hydraulic actuated boom where both hydraulic driven system and boom structure deformation are considered. In this paper, the hydraulic actuated boom system is formulated as a port-Hamiltonian system which is a proper modelling method for multi-domain system. The problems of trajectory optimization and vibration control are formulated as optimal control problem based on port-Hamiltonian model and this procedure is tested on a model of hydraulic cylinder. A reasonable result is solved with the selected cost function and inputs.

[pt] O espaço de polígonos de um grupo de Lie é definido como a redução simplética em um produto de órbitas pela ação coadjunta. Neste trabalho comparamos alguns sistemas integráveis definidos em espaços de módulos de polígonos, matrizes e fibrados, tais como o sistema de Kapovich–Millson, o modelo de Gaudin e a aplicação de Hitchin. / [en] The Polygon Space of a Lie group is defined as the symplectic reduction of a product of orbits by the coadjoint action. In this work we compare integrable systems defined on different moduli spaces of polygons, matrices and bundles, such as Kapovich–Millson s system, Gaudin s model and the Hitchin s map.

There are many problems that lead to analysis of dynamical systems with phase variables of two types, slow and fast ones. Such systems are called slow-fast systems. The dynamics of such systems is usually described by means of different versions of perturbation theory. Many questions about accuracy of this description are still open. The difficulties are related to presence of resonances. The goal of the proposed thesis is to establish some estimates of the accuracy of the perturbation theory for slow-fast systems in the presence of resonances. We consider slow-fast Hamiltonian systems and study an accuracy of one of the methods of perturbation theory: the averaging method. In this thesis, we start with the case of slow-fast Hamiltonian systems with two degrees of freedom. One degree of freedom corresponds to fast variables, and the other degree of freedom corresponds to slow variables. Action variable of fast sub-system is an adiabatic invariant of the problem. Let this adiabatic invariant have limiting values along trajectories as time tends to plus and minus infinity. The difference of these two limits for a trajectory is known to be exponentially small in analytic systems. We obtain an exponent in this estimate. To this end, by means of iso-energetic reduction and canonical transformations in complexified phase space, we reduce the problem to the case of one and a half degrees of freedom, where the exponent is known. We then consider a quasi-linear Hamiltonian system with one and a half degrees of freedom. The Hamiltonian of this system differs by a small, ~ε, perturbing term from the Hamiltonian of a linear oscillatory system. We consider passage through a resonance: the frequency of the latter system slowly changes with time and passes through 0. The speed of this passage is of order of ε. We provide asymptotic formulas that describe effects of passage through a resonance with an improved accuracy O(ε3/2). A numerical verification is also provided. The problem under consideration is a model problem that describes passage through an isolated resonance in multi-frequency quasi-linear Hamiltonian systems. We also discuss a resonant phenomenon of scattering on resonances associated with discretisation arising in a numerical solving of systems with one rotating phase. Numerical integration of ODEs by standard numerical methods reduces continuous time problems to discrete time problems. For arbitrarily small time step of a numerical method, discrete time problems have intrinsic properties that are absent in continuous time problems. As a result, numerical solution of an ODE may demonstrate dynamical phenomena that are absent in the original ODE. We show that numerical integration of systems with one fast rotating phase leads to a situation of such kind: numerical solution demonstrates phenomenon of scattering on resonances, that is absent in the original system.

In this thesis we investigate the existence and stability of periodic solutions of Hamiltonian systems with a discrete symmetry. The global existence of periodic motions can be proven using the classical techniques of the calculus of variations; our particular interest is in how the stability type of the solutions thus obtained can be determined analytically using solely the variational problem and the symmetries of the system -- we make no use of numerical or perturbation techniques. Instead, we use a method introduced in [41] in the context of a special case of the three-body problem. Using techniques from symplectic geometry, and specifically the Maslov index for curves of Lagrangian subspaces along the minimizing trajectories, we verify conditions which preclude the existence of eigenvalues of the monodromy matrix on the unit circle. We study the applicability of this method in two specific cases. Firstly, we consider another special case from celestial mechanics: the hip-hop solutions of the 2N-body problem. This is a family of Z_2-symmetric, periodic orbits which arise as collision-free minimizers of the Lagrangian action on a space of symmetric loops [14, 53]. Following a symplectic reduction, it is shown that the hip-hop solutions are brake orbits which are generically hyperbolic on the reduced energy-momentum surface. Secondly we consider a class of natural Hamiltonian systems of two degrees of freedom with a homogeneous potential function. The associated action functional is unbounded above and below on the function space of symmetric curves, but saddle points can be located by minimization subject to a certain natural constraint of a type first considered by Nehari [37, 38]. Using the direct method of the calculus of variations, we prove the existence of symmetric solutions of both prescribed period and prescribed energy. In the latter case, we employ a variational principle of van Groesen [55] based upon a modification of the Jacobi functional, which has not been widely used in the literature. We then demonstrate that the (constrained) minimizers are again hyperbolic brake orbits; this is the first time the method has been applied to solutions which are not globally minimizing. / Thesis (Ph.D, Mathematics & Statistics) -- Queen's University, 2013-09-25 10:47:53.257

A pair of nonidentical two-level atoms, separated by a fixed distance R, interact through photon exchange. The system is described by a state vector which is assumed to be a superposition of four "essential states": (1) the first atom is excited, the second one is in the ground state, and no photon is present, (2) the first atom is in its ground state, the second one is excited, and no photon is present, (3) both atoms are in their ground states and a photon is present, and (4) both atoms are excited and a photon is also present. The system is initially in state (1). The probabilities of each atom being excited are calculated for both the minimally-coupled interaction and the multipolar interaction in the electric dipole approximation. For the minimally-coupled interaction Hamiltonian, the second atom has a probability of being instantaneously excited, so the interaction is not retarded. For the multipolar interaction Hamiltonian, the second atom is not excited before the retardation time, which agrees with special relativity. For the minimally-coupled interaction the nonphysical result occurs because the unperturbed Hamiltonian is not the energy operator in the Coulomb gauge. For the multipolar Hamiltonian in the electric dipole approximation the unperturbed Hamiltonian is the energy operator. An active view of unitary transformations in nonrelativistic quantum electrodynamics is used to derive transformation laws for the potentials of the electromagnetic field and the static Coulomb potential. For a specific choice of unitary transformation the transformation laws for the potentials are used in the minimally-coupled second-quantized Hamiltonian to obtain the multipolar Hamiltonian, which is expressed in terms of the quantized electric and magnetic fields.

Dans cette thèse, nous discutons l'influence d'un réseau qui possède une topologie non triviale sur les propriétés collectives d'un modèle hamiltonien pour spins,le modèle $XY$, défini sur ces réseaux.Nous nous concentrons d'abord sur la topologie des chaînes régulières et du réseau Petit Monde (Small World), créé avec le modèle Watt- Strogatz.Nous contrôlons ces réseaux par deux paramètres $\gamma$, pour le nombre d' interactions et $p$, la probabilité de ré-attacher un lien aléatoirement.On définit deux mesures, le chemin moyen $\ell$ et la connectivité $C$ et nous analysons leur dépendance de $(\gamma,p)$.Ensuite,nous considérons le comportement du modèle $XY$ sur la chaîne régulière et nous trouvons deux régimes: un pour $\gamma<1,5$,qui ne présente pas d'ordre longue portée et un pour $\gamma>1,5$ où une transition de phase du second ordre apparaît.Nous observons l'existence d'un état ​​métastable pour $\gamma_ {c} = 1,5$. Sur les réseaux Petit Monde,nous illustrons les conditions pour avoir une transition et comment son énergie critique $\varepsilon_{c}(\gamma,p)$ dépend des paramètres $(\gammap$).Enfin,nous proposons un modèle de réseau où les liens d'une chaîne régulière sont ré-attachés aléatoirement avec une probabilité $p$ dans un rayon spécifique $r$. Nous identifions la dimension du réseau $d(p,r)$ comme un paramètre crucial:en le variant,il nous est possible de passer de réseaux avec $d<2$ qui ne présentent pas de transition de phase à des configurations avec $d>2$ présentant une transition de phase du second ordre, en passant par des régimes de dimension $d=2$ qui présentent des états caractérisés par une susceptibilité infinie et une dynamique chaotique. / In this thesis we discuss the influence of a non trivial network topology on the collective properties of an Hamiltonian model defined on it, the $XY$ -rotors model. We first focus on networks topology analysis, considering the regular chain and a Small World network, created with the Watt-Strogatz model. We parametrize these topologies via $\gamma$, giving the vertex degree and $p$, the probability of rewiring. We then define two topological parameters, the average path length $\ell$and the connectivity $C$ and we analize their dependence on $\gamma$ and $p$. Secondly, we consider the behavior of the $XY$- model on the regular chain and we find two regimes: one for $\gamma<1.5$, which does not display any long-range order and one for $\gamma>1.5$ in which a second order phase transition of the magnetization arises. Moreover we observe the existence of a metastable state appearing for $\gamma_{c}=1.5$. Finally we illustrate in what conditions we retrieve the phase transition on Small World networks and how its critical energy $\varepsilon_{c}(\gamma,p)$ depends on the topological parameters $\gamma$ and $p$. In the last part, we propose a network model in which links of a regular chain are rewired according to a probability $p$ within a specific range $r$. We identify a quantity, the network dimension $d(p,r)$ as a crucial parameter. Varying this dimension we are able to cross over from topologies with $d<2$ exhibiting no phase transitions to ones with $d>2$ displaying a second order phase transition, passing by topologies with dimension $d=2$ which exhibit states characterized by infinite susceptibility and macroscopic chaotic dynamical behavior.

A presente tese analisa diferentes aspectos de sistemas de espalhamento clássico com caos. Espalhamento caótico é uma forma de caos transiente que ocorre em diversos sistemas físicos. Nestes sistemas o espaço de fase é aberto, mas o caos ocorre apenas em uma região restrita do espaço, chamada de região de espalhamento. Os efeitos desta dinâmica apresentam-se em qualquer relação de espalhamento pela presença de conjuntos fractais, que geram hiper-sensibilidade a condições iniciais. Em nosso primeiro trabalho, mostramos que as bifurcações que levam ao caos manifestam-se na Seção de Choque Diferencial (SCD) pela criação de infinitas singularidades arco-íris. Estas singularidades aparecem na forma de cascatas, registrando na SCD todas as transições sofridas pela sela caótica. O segundo trabalho mostra que a introdução de dissipação em sistemas de espalhamento pode limitar a autosimilaridade de conjuntos originalmente fractais. Uma partícula espalhada por potenciais repulsivos encontra regiões não acessíveis, que dependem do valor de sua energia. Estas regiões determinam a estrutura da sela caótica. Com a perda de energia, o cenário de órbitas presas é alterado e, dependendo do valor da dissipação, podem existir nas funções de espalhamento estruturas fractais truncadas. O terceiro estudo aborda a presença de advecção caótica em fluxos sanguíneos. Doenças circulatórias estão geralmente associadas a uma mudança de geometria de artérias ou veias. Essas deformações podem gerar espalhamento caótico das partículas sanguíneas carregadas pelo fluxo. Em nosso trabalho mostramos, a partir de simulações numéricas, que caos pode existir em fluxos sanguíneos e, assim, formar um ciclo no desenvolvimento de anomalias circulatórias. / In this thesis we study different scattering systems with chaos. Chaotic scattering, present in a large variety of physical systems, is a type of transient chaos. While the phase-space of such systems is unbounded, irregular motion occurs only in a bounded area, called the scattering region. Still, any (nontrivial) scattering function relating initial conditions to asymptotic variables contains fractal structures, resulting in a very sharp sensitivity to initial conditions. Our first work shows that bifurcations leading to chaos manifest themselves through an infinitely fine-scale structure of rainbow singularities in the cross section. These singularities appear as cascades, mirroring the bifurcation cascade undergone by the chaotic saddle. The second work shows that the presence of dissipation in scattering systems can limit the auto-similarity of originally fractal structures. Depending on the value of their energy, particles scattered by repulsive potentials find forbidden regions in the space-phase. These regions determinate the structure of the chaotic saddle. With friction, the scenario of trapped orbits changes and, depending on the ammount dissipation, scattering functions follow a truncated fractal structure. Our third study concerns the presence of chaotic advection in blood flows. Typically, circulatory diseases are due to sudden changes on the geometry of vessel walls. These deformations can generate chaotic scattering of blood particles carried by the flow. We show, with numerical simulations, that chaos can occur in blood flows and thus form a hazardous cycle in the further developing of circulatory anomalies.

Trente ans après la démonstration de la production d'harmoniques laser par interaction laser-gaz non linéaire, la génération d'harmoniques d’ordre élevées (HHG) est utilisée pour sonder la dynamique moléculaire et réalise son potentiel technologique comme source compacte d'impulsions attosecondes XUV à la gamme de rayons X. Malgré les progrès expérimentaux, le coût de calcul excessif des simulations fondées sur les premiers principes et la difficulté de dériver systématiquement des modèles réduits pour l'interaction non perturbatif et à échelles multiples d'une impulsion laser intense avec un gaz macroscopique d'atomes ont entravé les efforts théoriques. Dans cette thèse, nous étudions des modèles réduits de premier principe pour HHG utilisant la mécanique classique. En utilisant la dynamique non linéaire, nous élucidons le rôle indispensable joué par le potentiel ionique lors des recollisions dans la limite du champ fort. Ensuite, en empruntant une technique de la physique des plasmas, nous dérivons systématiquement une hiérarchie de modèles hamiltoniens réduits pour l’interaction cohérente entre le laser et les atomes lors de la propagation des impulsions. Les modèles réduits permettent une dynamique électronique soit classique, soit quantique. Nous construisons un modèle classique qui concorde quantitativement avec le modèle quantique pour la propagation des composantes dominantes du champ laser. Dans une géométrie simplifiée, nous montrons que le rayonnement à fréquence anormalement élevée observé dans les simulations résulte de l’interaction délicate entre le piégeage d’électrons et les recollisions de plus grande énergie provoqués par les effets de propagation. / Thirty years after the demonstration of the production of high laser harmonics through nonlinear laser-gas interaction, high harmonic generation (HHG) is being used to probe molecular dynamics in real time and is realizing its technological potential as a tabletop source of attosecond pulses in the XUV to soft X-ray range. Despite experimental progress, theoretical efforts have been stymied by the excessive computational cost of first-principles simulations and the difficulty of systematically deriving reduced models for the non-perturbative, multiscale interaction of an intense laser pulse with a macroscopic gas of atoms. In this thesis, we investigate first-principles reduced models for HHG using classical mechanics. Using nonlinear dynamics, we elucidate the indispensable role played by the ionic potential during recollisions in the strong-field limit. Then, borrowing a technique from plasma physics, we systematically derive a hierarchy of reduced Hamiltonian models for the self-consistent interaction between the laser and the atoms during pulse propagation. The reduced models can accommodate either classical or quantum electron dynamics. We build a classical model which agrees quantitatively with the quantum model for the propagation of the dominant components of the laser field. In a simplified geometry, we show that the anomalously high frequency radiation seen in simulations results from the delicate interplay between electron trapping and higher energy recollisions brought on by propagation effects.

In this thesis we will discuss connections between Hamiltonian systems with a periodic kick and planar diffeomorphisms. After a brief overview of Hamiltonian theory we will focus, as an example, on derivations of the Hâenon map that can be obtained by considering kicked Hamiltonian systems. We will conclude with examples of Hâenon maps of interest. / by Zalmond C. Barney. / Thesis (M.S.)--Florida Atlantic University, 2011. / Includes bibliography. / Electronic reproduction. Boca Raton, Fla., 2011. Mode of access: World Wide Web.