Attractors of the sine-Gordon equation in the field of a quasiperiodic external force

The well-known sine-Gordon equation, supplemented with small damping and small quasiperiodic external force, is studied under the zero Dirichlet boundary conditions at the endpoints of a finite interval. The main assumption is that all frequencies of the external force are in 1:1 resonance with certain eigenfrequencies of the unperturbed equation; i.e., the socalled fundamental multifrequency resonance is observed. It is shown that in this case, by an appropriate choice of the parameters of the external force, one can make it so that the boundary value problem has a stable invariant torus of any finite dimension that bifurcates from zero on any preassigned finite set of spatial modes. It is also shown (by numerical analysis) that in a number of cases the above-mentioned torus coexists with a chaotic attractor.     

Bifurcations of resonant solutions and chaos in physical pendulum equation with suspension axis vibrations

This paper is a continuation of ＂Complex Dynamics in Physical Pendulum Equation with Suspension Axis Vibrations＂[1].In this paper,we investigate the existence and the bifurcations of resonant solution for ω0：ω：Ω ≈ 1：1：n,1：2：n,1：3：n,2：1：n and 3：1：n by using second-order averaging method,give a criterion for the existence of resonant solution for ω0：ω：Ω ≈ 1：m：n by using Melnikov＇s method and verify the theoretical analysis by numerical simulations.By numerical simulation,we expose some other interesting dynamical behaviors including the entire invariant torus region,the cascade of invariant torus behaviors,the entire chaos region without periodic windows,chaotic region with complex periodic windows,the entire period-one orbits region;the jumping behaviors including invariant torus behaviors converting to period-one orbits,from chaos to invariant torus behaviors or from invariant torus behaviors to chaos,from period-one to chaos,from invariant torus behaviors to another invariant torus behaviors;the interior crisis;and the different nice invariant torus attractors and chaotic attractors.The numerical results show the difference of dynamical behaviors for the physical pendulum equation with suspension axis vibrations between the cases under the three frequencies resonant condition and under the periodic/quasi-periodic perturbations.It exhibits many invariant torus behaviors under the resonant conditions.We find a lot of chaotic behaviors which are different from those under the periodic/quasi-periodic perturbations.However,we did not find the cascades of period-doubling bifurcation.     

Solving elliptical equations with high-contrast coefficients in an unbounded region

An approximation model is proposed for an elliptical equation with complex rapidly varying coefficients. An efficient numerical method is developed and implemented. A problem of geoelectricity requiring solution of an equation in this setting is investigated. This research was partially supported by the Russian Foundation for Basic Research (grant No. 96-05-64340) and by the Interuniversity Scientific Program “Russian Universities: Basic Research.” Translated from Chislennye Metody v Matematicheskoi Fizike, Moscow State University, pp. 37–45, 1998.     

Applied aspects of the theory of continuation of solutions of wave-scattering boundary-value problems

The article compares two approaches to continuation of solutions of external boundary-value problems for the Helmholtz equation outside the original region in the solution space — analytical continuation and the generalized image method. The comparison focuses on naturalness and ease of use of the mathematical model for the construction of numerical algorithms. Translated from Chislennye Metody v Matematicheskoi Fizike, Published by Moscow University, Moscow, 1996, pp. 68–79.     

Similarity variables and reduction of the heat equation on torus

The problem of symmetry classification for the heat equation on torus is studied by means of classical Lie group theory. The Lie point symmetries are constructed and Lie algebra is formed for equation under consideration. Then these algebras are used to classify its subalgebras up to conjugacy classes. In general the heat equation on torus admits one-, two-, three- and four-dimensional algebras. For one-dimensional algebra £₁ and £₂ the heat equation on torus is reduced in independent variables whereas in two-dimensional algebras £₃ and £₄ the considered heat equation is investigated by quadrature. While three- and four-dimensional algebras lead to a trivial solution.     

Finite abelian group methods in iterative algorithms for solving boundary-value problems with a boundary surface of complex structure

We consider external boundary-value problems for the Laplace equation on surfaces of complex structure. Various iterative computational schemes are constructed for numerical solution of the corresponding integral equations using set-theoretical group methods. Simulation results of electron-optical images are reported. Translated from Chislennye Metody v Matematicheskoi Fizike, Published by Moscow University, Moscow, 1996, pp. 16–27.     

From Canards of Folded Singularities to Torus Canards in a Forced van der Pol Equation

In this article, we study canard solutions of the forced van der Pol equation in the relaxation limit for low-, intermediate-, and high-frequency periodic forcing. A central numerical observation made herein is that there are two branches of canards in parameter space which extend across all positive forcing frequencies. In the low-frequency forcing regime, we demonstrate the existence of primary maximal canards induced by folded saddle nodes of type I and establish explicit formulas for the parameter values at which the primary maximal canards and their folds exist. Then, we turn to the intermediate- and high-frequency forcing regimes and show that the forced van der Pol possesses torus canards instead. These torus canards consist of long segments near families of attracting and repelling limit cycles of the fast system, in alternation. We also derive explicit formulas for the parameter values at which the maximal torus canards and their folds exist. Primary maximal canards and maximal torus canards correspond geometrically to the situation in which the persistent manifolds near the family of attracting limit cycles coincide to all orders with the persistent manifolds that lie near the family of repelling limit cycles. The formulas derived for the folds of maximal canards in all three frequency regimes turn out to be representations of a single formula in the appropriate parameter regimes, and this unification confirms the central numerical observation that the folds of the maximal canards created in the low-frequency regime continue directly into the folds of the maximal torus canards that exist in the intermediate- and high-frequency regimes. In addition, we study the secondary canards induced by the folded singularities in the low-frequency regime and find that the fold curves of the secondary canards turn around in the intermediate-frequency regime, instead of continuing into the high-frequency regime. Also, we identify the mechanism responsible for this turning. Finally, we show that the forced van der Pol equation is a normal form-type equation for a class of single-frequency periodically driven slow/fast systems with two fast variables and one slow variable which possess a non-degenerate fold of limit cycles. The analytic techniques used herein rely on geometric desingularisation, invariant manifold theory, Melnikov theory, and normal form methods. The numerical methods used herein were developed in Desroches et al. (SIAM J Appl Dyn Syst 7:1131–1162, 2008, Nonlinearity 23:739–765 2010).     

Topological Bifurcations of Attracting 2-Tori of Quasiperiodically Driven Oscillators

We examine the solutions to a damped, quasiperiodic (QP) Mathieu equation with cubic nonlinearities. The system is suspended in a four-dimensional phase space ℝ² × T² in which there exist attracting, knotted 2-tori called torus braids. We develop a topological classification scheme in which a torus braid is characterized by closed braids that exist in two Poincare sections, ℝ² \times S¹ × {·} and ℝ² × {·} \times S¹. Based on the classification scheme, we develop numerical invariants that describe the linkedness of attractors and provide information about the global dynamics. Numerical simulations show that changes of a single parameter lead to a global bifurcation through which the attracting torus loses stability and locally "doubles," shedding a torus of different equivalence class. We call this a topological torus bifurcation of the doubling variety (TTBD). We provide a topological analysis of the doubling produced by TTBDs and we examine the qualitative dynamical changes that result. We also examine the effect of TTBDs on the spectrum of Lyapunov exponents and the time series power spectrum.     

Commuting Flows and Infinite-Dimensional Tori: Sine-Gordon

The present work concerns the periodic sine-Gordon equation. We explain why the complete set of conserved functionals for sine-Gordon is an infinite-dimensional torus; the periodic sine-Gordon solution is almost periodic in time on an infinite-dimensional torus.     

Inverse coefficient problem for a nonlinear differential equation and an iterative solution method

The article considers the determination of the solution-dependent coefficient of a nonautonomous ordinary differetial equation with a parameter. Reduction of the inverse problem to a nonlinear operator equation is used to prove existence and uniqueness theorems and to propose an iterative solution method. Translated from Obratnye Zadachi Estestvoznaniya, Published by Moscow University, Moscow, 1997, pp. 5–17.     

Normal form and long time analysis of splitting schemes for the linear Schrödinger equation with small potential

We consider the linear Schrödinger equation on a one dimensional torus and its time-discretization by splitting methods. Assuming a non-resonance condition on the stepsize and a small size of the potential, we show that the numerical dynamics can be reduced over exponentially long time to a collection of two dimensional symplectic systems for asymptotically large modes. For the numerical solution, this implies the long time conservation of the energies associated with the double eigenvalues of the free Schrödinger operator. The method is close to standard techniques used in finite dimensional perturbation theory, but extended here to infinite dimensional operators.     

Insulated antenna in a layered medium

The article presents a generalization of the integral equation of an insulated linear antenna immersed in a cylindrically layered lossy dielectric medium. The insulation is provided by a lossless dielectric layer. The kernel of the integral equation is represented as a superposition of the fundamental solutions of the wave equation with equivalent propagation constants for the given media. A generalization to a plane-layered medium is proposed. The problem of a vertical radiator above a layered half-space is considered. Translated from Chislennye Metody v Matematicheskoi Fiziki, Published by Moscow University, Moscow, 1996, pp. 80–88.     

Semidiscretisation Methods for Quasi‐periodic Solutions

The aim of our approach is a reliable numerical approximation of quasi‐periodic solutions of periodically forced ODEs without using a‐priori transformations into new coordinates [1]. The invariant torus is computed as a solution of a special invariance equation. In the case of two basic frequencies this system can be solved by semidiscretisation, which transforms the system into a higher dimensional autonomous ODE system with periodic solutions.     

Study of oscillations in a viscoelastic cylindrical shell in an elastic medium

An equation was obtained and solved for oscillations in a viscoelastic cylindrical shell located in the ground. A numerical example is examined.Moscow Electronic Engineering Institute. Translated from Mekhanika Polimerov, No. 1, pp. 178–181, January–February, 1974.     

Order Conditions for Sampling the Invariant Measure of Ergodic Stochastic Differential Equations on Manifolds

We derive a new methodology for the construction of high-order integrators for sampling the invariant measure of ergodic stochastic differential equations with dynamics constrained on a manifold. We obtain the order conditions for sampling the invariant measure for a class of Runge–Kutta methods applied to the constrained overdamped Langevin equation. The analysis is valid for arbitrarily high order and relies on an extension of the exotic aromatic Butcher-series formalism. To illustrate the methodology, a method of order two is introduced, and numerical experiments on the sphere, the torus and the special linear group confirm the theoretical findings.

     

Conformal pasting on a torus

We solve the problem concerning global conformal pasting on a torus given by the algebraic equation $$u^2 = (1 - z^2 ) (1 - k^z 2^2 ) (0< k< 1).$$ We obtain an algebraic equation for the new torus, and we find the function which accomplishes the conformal pasting.     

On the coexistence of periodic or quasi-periodic oscillations near a Hopf-pitchfork bifurcation in NFDE

Normal form method is first employed to study the Hopf-pitchfork bifurcation in neutral functional differential equation (NFDE), and is proved to be an efficient approach to show the rich dynamics (periodic and quasi-periodic oscillations) around the bifurcation point. We give an algorithm for calculating the third-order normal form in NFDE models, which naturally arise in the method of extended time delay autosynchronization (ETDAS). The existence of Hopf-pitchfork bifurcation in a van der Pol’s equation with extended delay feedback is given and the unfoldings near this critical point is obtained by applying our algorithm. Some interesting phenomena, such as the coexistence of several stable periodic oscillations (or quasi-periodic oscillations) and the existence of saddle connection bifurcation on a torus, are found by analyzing the bifurcation diagram and are illustrated by numerical method.     

A Five-Dimensional System of Navier-Stokes Equation for a Two-Dimensional Incompressible Fluid on a Torus

A five-mode truncation of Navier-Stokes equation for a two-dimensional incompressible fluid on a torus is studied. Its stationary solutions and stability are presented, the existence of attractor and the global stability of the system are discussed. The whole process, which shows a chaos behavior approached through an involved sequence of bifurcations with the changing of Reynolds number, is simulated numerically. Based on numerical simulation results of bifurcation diagram, Lyapunov exponent spectrum, Poincare section, power spectrum and return map of the system are revealed.     

A numerical method for solving the diffraction problem on a system of bodies of revolution

The article considers the diffraction of an acoustic wave on a system of bodies of revolution in a homogeneous space. The iteration method is applied to reduce the problem to diffraction on a single body of revolution, which is solved by the integral equation method with expansion of the unknown function in the harmonics of the azimuthal angle . A numerical implementation of the method is considered for a non-coaxial system consisting of a torus and two spheroids.Translated from Matematicheskoe Modelirovanie i Reshenie Obratnykh Zadach Matematicheskoi Fiziki, pp. 144–152, 1993.     

A canonical system of two differential equations with periodic coefficients and the Poincaré-Denjoy theory of differential equations on a torus

The passage from Cartesian to polar coordinates in a canonical system with periodic coefficients gives rise to a nonlinear differential equation whose right-hand side is periodic in time and the polar angle and thus this equation can be regarded as a differential equation on a torus. In accord with Poincaré-Denjoy theory, the behavior of a solution to a differential equation on a torus is characterized by the rotation number and some homeomorphic mapping of a circle onto itself. We study connections of strong stability (instability) of a canonical system, including the membership in the nth stability (instability) domain, with the rotation number and fixed points of this mapping.