Stochastic Hopf Bifurcation of MDOF Quasi-Integrable Hamiltonian Systems with Delayed Feedback Fractional-Order PD Controller

Hopf bifurcation, as the most representative dynamic bifurcation, is closely related to the stability of many engineering structures. In this work, the stochastic Hopf bifurcation (SHB) of a controlled quasi-integrable Hamiltonian system (H.S.) of multi-degree-of-freedom (MDOF) is investigated, where the system is subjected to wide-band noise and controlled by a Fractional-order Proportional-Derivative (FOPD) controller with time delay. By decoupling FOPD control force and simplifying it without time delay, the averaged Itô differential equations of the approximated system are derived with the technique of stochastic averaging. Then, the average bifurcation parameter expression of system is obtained, which can determine the criterion of the SHB deduced by the FOPD control force. Last, an illustration of coupled Rayleigh oscillators is given to demonstrate the validity of the procedure. The influences of time delay, noise intensities and fractional order on the system SHB are discussed.

Moduli Spaces of Generalized Hyperpolygons

We introduce the notion of generalized hyperpolygon, which arises as a representation, in the sense of Nakajima, of a comet-shaped quiver. We identify these representations with rigid geometric figures, namely pairs of polygons: one in the Lie algebra of a compact group and the other in its complexification. To such data, we associate an explicit meromorphic Higgs bundle on a genus-g Riemann surface, where g is the number of loops in the comet, thereby embedding the Nakajima quiver variety into a Hitchin system on a punctured genus-g Riemann surface (generally with positive codimension). We show that, under certain assumptions on flag types, the space of generalized hyperpolygons admits the structure of a completely integrable Hamiltonian system of Gelfand–Tsetlin type, inherited from the reduction of partial flag varieties. In the case where all flags are complete, we present the Hamiltonians explicitly. We also remark upon the discretization of the Hitchin equations given by hyperpolygons, the construction of triple branes (in the sense of Kapustin–Witten mirror symmetry), and dualities between tame and wild Hitchin systems (in the sense of Painlevé transcendents).

Numerical determination of the singularity order of a system of differential equations

We consider moving singular points of systems of ordinary differential equations. A review of Painlevé’s results on the algebraicity of these points and their relation to the Marchuk problem of determining the position and order of moving singularities by means of finite difference method is carried out. We present an implementation of a numerical method for solving this problem, proposed by N. N. Kalitkin and A. Al’shina (2005) based on the Rosenbrock complex scheme in the Sage computer algebra system, the package CROS for Sage. The main functions of this package are described and numerical examples of usage are presented for each of them. To verify the method, computer experiments are executed (1) with equations possessing the Painlevé property, for which the orders are expected to be integer; (2) dynamic Calogero system. This system, well-known as a nontrivial example of a completely integrable Hamiltonian system, in the present context is interesting due to the fact that coordinates and momenta are algebraic functions of time, and the orders of moving branching points can be calculated explicitly. Numerical experiments revealed that the applicability conditions of the method require additional stipulations related to the elimination of superconvergence points.

Numerical determination of the singularity order of a system of differential equations

We consider moving singular points of systems of ordinary differential equations. A review of Painlevé’s results on the algebraicity of these points and their relation to the Marchuk problem of determining the position and order of moving singularities by means of finite difference method is carried out. We present an implementation of a numerical method for solving this problem, proposed by N. N. Kalitkin and A. Al’shina (2005) based on the Rosenbrock complex scheme in the Sage computer algebra system, the package CROS for Sage. The main functions of this package are described and numerical examples of usage are presented for each of them. To verify the method, computer experiments are executed (1) with equations possessing the Painlevé property, for which the orders are expected to be integer; (2) dynamic Calogero system. This system, well-known as a nontrivial example of a completely integrable Hamiltonian system, in the present context is interesting due to the fact that coordinates and momenta are algebraic functions of time, and the orders of moving branching points can be calculated explicitly. Numerical experiments revealed that the applicability conditions of the method require additional stipulations related to the elimination of superconvergence points.

Numerical determination of the singularity order of a system of differential equations

We consider moving singular points of systems of ordinary differential equations. A review of Painlevé’s results on the algebraicity of these points and their relation to the Marchuk problem of determining the position and order of moving singularities by means of finite difference method is carried out. We present an implementation of a numerical method for solving this problem, proposed by N. N. Kalitkin and A. Al’shina (2005) based on the Rosenbrock complex scheme in the Sage computer algebra system, the package CROS for Sage. The main functions of this package are described and numerical examples of usage are presented for each of them. To verify the method, computer experiments are executed (1) with equations possessing the Painlevé property, for which the orders are expected to be integer; (2) dynamic Calogero system. This system, well-known as a nontrivial example of a completely integrable Hamiltonian system, in the present context is interesting due to the fact that coordinates and momenta are algebraic functions of time, and the orders of moving branching points can be calculated explicitly. Numerical experiments revealed that the applicability conditions of the method require additional stipulations related to the elimination of superconvergence points.

Nonlinear Stochastic Optimal Control Using Piezoelectric Stack Inertial Actuator

An optimal control strategy for the random vibration reduction of nonlinear structures using piezoelectric stack inertial actuator is proposed. First, the dynamic model of the nonlinear structure considering the dynamics of a piezoelectric stack inertial actuator is established, and the motion equation of the coupled system is described by a quasi-non-integrable-Hamiltonian system. Then, using the stochastic averaging method, this quasi-non-integrable-Hamiltonian system is reduced to a one-dimensional averaged system for total energy. The optimal control law is determined by establishing and solving the dynamic programming equation. The proposed control law is analytical and can be fully executed by a piezoelectric stack inertial actuator. The responses of optimally controlled and uncontrolled systems are obtained by solving the Fokker–Planck–Kolmogorov (FPK) equation to evaluate the control effectiveness of the proposed strategy. Numerical results show that our proposed control strategy is effective for random vibration reduction of the nonlinear structures using piezoelectric stack inertial actuator, and the theoretical method is verified by comparing with the simulation results.

Frequency analysis and representation of slowly diffusing planetary solutions

Context. Over short time-intervals, planetary ephemerides have traditionally been represented in analytical form as finite sums of periodic terms or sums of Poisson terms that are periodic terms with polynomial amplitudes. This representation is not well adapted for the evolution of planetary orbits in the solar system over million of years which present drifts in their main frequencies as a result of the chaotic nature of their dynamics. Aims. We aim to develop a numerical algorithm for slowly diffusing solutions of a perturbed integrable Hamiltonian system that will apply for the representation of chaotic planetary motions with varying frequencies. Methods. By simple analytical considerations, we first argue that it is possible to exactly recover a single varying frequency. Then, a function basis involving time-dependent fundamental frequencies is formulated in a semi-analytical way. Finally, starting from a numerical solution, a recursive algorithm is used to numerically decompose the solution into the significant elements of the function basis. Results. Simple examples show that this algorithm can be used to give compact representations of different types of slowly diffusing solutions. As a test example, we show that this algorithm can be successfully applied to obtain a very compact approximation of the La2004 solution of the orbital motion of the Earth over 40 Myr ([−35 Myr, 5 Myr]). This example was chosen because this solution is widely used in the reconstruction of the past climates.

Bifurcations of Liouville tori in a system of two vortices of positive intensity in a Bose-Einstein condensate

In this paper we consider a completely Liouville integrable Hamiltonian system with two degrees of freedom, which describes the dynamics of two vortex filaments in a Bose-Einstein condensate enclosed in a harmonic trap. For vortex pairs of positive intensity detected bifurcation of three Liouville tori into one. Such bifurcation was found in the integrable case of Goryachev-Chaplygin-Sretensky in the dynamics of a rigid body. For the integrable perturbation of the physical parameter of the intensity ratio, identified bifurcation proved to be unstable, which led to bifurcations of the type of two tori into one and vice versa.

Toric Geometry ofSL2(ℂ) Free Group Character Varieties from Outer Space

Culler and Vogtmann defined a simplicial spaceO(g), calledouter space, to study the outer automorphism group of the free groupFg. Using representation theoretic methods, we give an embedding ofO(g) into the analytification of X(Fg,SL2(ℂ)), theSL2(ℂ) character variety ofFg, reproving a result of Morgan and Shalen. Then we show that every pointvcontained in a maximal cell ofO(g) defines a flat degeneration of X(Fg,SL2(ℂ)) to a toric varietyX(PΓ). We relate X(Fg,SL2(ℂ)) andX(v) topologically by showing that there is a surjective, continuous, proper map Ξv:X(Fg,SL2(ℂ)) →X(v). We then show that this map is a symplectomorphism on a dense open subset of X(Fg, SL2(ℂ)) with respect to natural symplectic structures on X(Fg, SL2(ℂ)) andX(v). In this way, we construct an integrable Hamiltonian system in X(Fg, SL2(ℂ)) for each point in a maximal cell ofO(g), and we show that eachvdefines a topological decomposition of X(Fg, SL2(ℂ)) derived from the decomposition ofX(PΓ) by its torus orbits. Finally, we show that the valuations coming from the closure of a maximal cell inO(g) all arise as divisorial valuations built from an associated projective compactification of X(Fg, SL2(ℂ)).

Exploring the dynamics of ‘2P’ wakes with reflective symmetry using point vortices

Wakes formed behind bluff bodies frequently reveal complex patterns of coherent vortical structures, with emergence of streamwise spatial periodicity particularly in the mid-wake region. In some cases, the vortex positions also maintain symmetry about the wake centreline. For the case in which two pairs of vortices are generated per shedding cycle, thereby constituting the so-called ‘2P’ mode wake, assumptions of spatial periodicity and symmetry allow for development of a mathematically tractable model using the point-vortex approximation. Our previous work (Basu & Stremler, Phys. Fluids, vol. 27 (10), 2015, 103603) considered staggered 2P wake configurations with two glide-reflective pairs of vortices shed in each period. Here we investigate the dynamics of a spatially periodic point-vortex street consisting of two pairs of vortices arranged with reflective symmetry about the streamwise centreline. Because of the symmetry, it is possible to model the spatially periodic point-vortex dynamics as an integrable Hamiltonian system. For a particular choice of initial condition, the topological structure of the Hamiltonian level curves is determined by location in a circulation–impulse parameter space. These Hamiltonian level curves delineate multiple regimes of motion, with all vortex motions within one regime being qualitatively identical. This approach thus enables identification and a full classification of all possible vortex motions in this constrained system. There also exist a limited number of equilibrium configurations with no relative vortex motion; some of these relative equilibria are neutrally stable to (appropriate) perturbations. Only one such neutrally stable equilibrium configuration continues to preserve the distinct four-vortex array, and numerical experiments indicate that these configurations are also neutrally stable to small perturbations that break the spatial symmetry. We apply this analysis to identify the parameter values necessary for co-existence of two closely spaced, neutrally stable Kármán vortex streets that preserve the assumed symmetry. Finally, comparison of the model dynamics to a wake pattern reported in the literature suggests that the classification of exotic wakes should be based on more details than just the number of vortices periodically shed by the body.