Publications

We suggest a universal map capable of recovering the behavior of a wide range of dynamical systems given by ODEs. The map is built as an artificial neural network whose weights encode a modeled system. We assume that ODEs are known and prepare training datasets using the equations directly without computing numerical time series. Parameter variations are taken into account in the course of training so that the network model captures bifurcation scenarios of the modeled system. The theoretical benefit from this approach is that the universal model admits applying common mathematical methods without needing to develop a unique theory for each particular dynamical equations. From the practical point of view the developed method can be considered as an alternative numerical method for solving dynamical ODEs suitable for running on contemporary neural network specific hardware. We consider the Lorenz system, the Rцssler system and also the Hindmarch – Rose model. For these three examples the network model is created and its dynamics is compared with ordinary numerical solutions. A high similarity is observed for visual images of attractors, power spectra, bifurcation diagrams and Lyapunov exponents.

We study the hyperchaos formation scenario in the modified Anishchenko–Astakhov generator. The scenario is connected with the existence of sequence of secondary torus bifurcations of resonant cycles preceding the hyperchaos emergence. This bifurcation cascade leads to the birth of the hierarchy of saddle-focus cycles with a two-dimensional unstable manifold as well as of saddle hyperchaotic sets resulting from the period-doubling cascades of unstable resonant cycles. Hyperchaos is born as a result of an inverse cascade of bifurcations of the emergence of discrete spiral Shilnikov attractors, accompanied by absorbing the cycles constituting this hierarchy.

Appearance of chaotic dynamics as a result of multi-frequency tori destruction is carried out on the example of a model of a multimode generator. Quasiperiodic bifurcations occurring with multi-frequency tori are discussed in the context of the Landau-Hopf scenario. Structure of the parameter space is studied, areas with various chaotic dynamics, including chaos and hyperchaos, are revealed. Scenarios of the development of chaotic dynamics are described, the features of chaotic signals of various types are revealed. 

In this paper, we consider a class of orientation-preserving Morse–Smale diffeomorphisms defined on an orientable surface. The papers by Bezdenezhnykh and Grines showed that such diffeomorphisms have a finite number of heteroclinic orbits. In addition, the classification problem for such diffeomorphisms is reduced to the problem of distinguishing orientable graphs with substitutions describing the geometry of a heteroclinic intersection. However, such graphs generally do not admit polynomial discriminating algorithms. This article proposes a new approach to the classification of these cascades. For this, each diffeomorphism under consideration is associated with a graph that allows the construction of an effective algorithm for determining whether graphs are isomorphic. We also identified a class of admissible graphs, each isomorphism class of which can be realized by a diffeomorphism of a surface with an orientable heteroclinic. The results obtained are directly related to the realization problem of homotopy classes of homeomorphisms on closed orientable surfaces. In particular, they give an approach to constructing a representative in each homotopy class of homeomorphisms of algebraically finite type according to the Nielsen classification, which is an open problem today.

In this paper, we consider regular topological flows on closed n-manifolds. Such flows have a hyperbolic (in the topological sense) chain recurrent set consisting of a finite number of fixed points and periodic orbits. The class of such flows includes, for example, Morse – Smale flows, which are closely related to the topology of the supporting manifold. This connection is provided by the existence of the Morse –Bott energy function for the Morse – Smale flows. It is well known that, starting from dimension 4, there exist nonsmoothing topological manifolds, on which dynamical systems can be considered only in a continuous category. The existence of continuous analogs of regular flows on any topological manifolds is an open question, as is the existence of energy functions for such flows. In this paper, we study the dynamics of regular topological flows, investigate the topology of the embedding and the asymptotic behavior of invariant manifolds of fixed points and periodic orbits. The main result is the construction of the Morse –Bott energy function for such flows, which ensures their close connection with the topology of the ambient manifold.

In this article we study the plasma motion in the transitional layer of a coronal loop randomly driven at one of its footpoints in the thin-tube and thin-boundary-layer (TTTB) approximation. We introduce the average of the square of a random function with respect to time. This average can be considered as the square of the oscillation amplitude of this quantity. Then we calculate the oscillation amplitudes of the radial and azimuthal plasma displacement as well as the perturbation of the magnetic pressure. We find that the amplitudes of the plasma radial displacement and the magnetic-pressure perturbation do not change across the transitional layer. The amplitude of the plasma radial displacement is of the same order as the driver amplitude. The amplitude of the magnetic-pressure perturbation is of the order of the driver amplitude times the ratio of the loop radius to the loop length squared. The amplitude of the plasma azimuthal displacement is of the order of the driver amplitude times Re1/6, where Re is the Reynolds number. It has a peak at the position in the transitional layer where the local Alfvén frequency coincides with the fundamental frequency of the loop kink oscillation. The ratio of the amplitude near this position and far from it is of the order of , where is the ratio of thickness of the transitional layer to the loop radius. We calculate the dependence of the plasma azimuthal displacement on the radial distance in the transitional layer in a particular case where the density profile in this layer is linear  

This paper is a continuation of research in the direction of energy function (a smooth Lyapunov function whose set of critical points coincides with the chain recurrent set of a system) construction for discrete dynamical systems. The authors established the existence of an energy function for any AA-diffeomorphism of a three-dimensional closed orientable manifold whose non-wandering set consists of a chaotic one-dimensional canonically embedded surface attractor and repeller.

Transition to chaos via the destruction of a two-dimensional torus is studied numerically using an example of the Hénon map and the Toda oscillator under quasiperiodic forcing and also experimentally using an example of a quasi-periodically excited RL–diode circuit. A feature of chaotic dynamics in these systems is the fact that the chaotic attractor in them has an additional zero Lyapunov exponent, which strictly follows from the structure of mathematical models. In the process of research, the influence of feedback is studied, in which the frequency of one of the harmonics of external forcing becomes dependent on a dynamic variable. Charts of dynamic regimes were constructed, examples of typical oscillation modes were given, and the spectrum of Lyapunov exponents was analyzed. Numerical simulations confirm that chaos resulting from the cascade of torus doubling has a close to the zero Lyapunov exponent, beside the trivial zero exponent.

The role of various long-wave approximations in the description of the wave field and bottom pressure caused by surface waves, and their relation to evolution equations are being considered. In the framework of the linear theory, these approximations are being tested on the well-known exact solution for the wave spectral amplitudes and pressure variations. The famous Whitham, Korteweg–de Vries (KdV) and Benjamin–Bona–Mahony (BBM) equations have been used as evolutionary equations. It has been shown that if the wave is long, though steep enough, the BBM approximation gives better results than the KdV approximation, and they are quite close to the exact results. The same applies to the description of rogue waves, though formed from smooth relatively long waves, are often short and steep, then they may be invisible in variations of the bottom pressure. Another advantage of the BBM approximation for calculating the bottom pressure is the ability to analyze noisy series without preliminary filtering, which is necessary when using the KdV approximation.

Meteotsunamis are long waves generated by displacement of a water body due to atmospheric pressure disturbances that have similar spatial and temporal characteristics to landslide tsunamis. NAMI DANCE that solves the nonlinear shallow water equations is a widely used numerical model to simulate tsunami waves generated by seismic origin. Several validation studies showed that it is highly capable of representing the generation, propagation and nearshore amplification processes of tsunami waves, including inundation at complex topography and basin resonance. The new module of NAMI DANCE that uses the atmospheric pressure and wind forcing as the other inputs to simulate meteotsunami events is developed. In this paper, the analytical solution for the generation of ocean waves due to the propagating atmospheric pressure disturbance is obtained. The new version of the code called NAMI DANCE SUITE is validated by comparing its results with those from analytical solutions on the flat bathymetry. It is also shown that the governing equations for long wave generation by atmospheric pressure disturbances in narrow bays and channels can be written similar to the 1D case studied for tsunami generation and how it is integrated into the numerical model. The analytical solution of the linear shallow water model is defined, and results are compared with numerical solutions. A rectangular shaped flat bathymetry is used as the test domain to model the generation and propagation of ocean waves and the development of Proudman resonance due to moving atmospheric pressure disturbances. The simulation results with different ratios of pressure speed to ocean wave speed (Froude numbers) considering sub-critical, critical and super-critical conditions are presented. Fairly well agreements between analytical solutions and numerical solutions are obtained. Additionally, basins with triangular (lateral) and stepwise shelf (longitudinal) cross sections on different slopes are tested. The amplitudes of generated waves at different time steps in each simulation are presented with discussions considering the channel characteristics. These simulations present the capability of NAMI DANCE SUITE to model the effects of bathymetric conditions such as shelf slope and local bathymetry on wave amplification due to moving atmospheric pressure disturbances.

On December 22, 2018, a destructive tsunami related to the phenomena caused by the volcanic eruption of Gunung Anak Krakatau (GAK) was generated following a partial collapse of the volcano that caused serious damage and killed more than 400 people. This recent event challenged the traditional understanding of tsunami hazard, warning and response mechanisms and raised the topic of volcanic tsunami hazard. The complex mechanism of this tsunamigenic volcano collapse still needs further investigation as Anak Krakatau is one of the potentially tsunamigenic volcanoes in the world. This study investigates the possible source mechanisms of this phenomenon and their contribution to explaining the observed sea level disturbances by considering the impacts of this destructive event. We configured a flank collapse scenario with a volume of 0.25 km3 as a combination of submarine and subaerial mass movement as the possible source scenarios to the December 22, 2018 Sunda Strait tsunami. A two-layer model is applied to simulate the tsunami generation by these landslides up to 420 s. The tsunami propagation and inundation are then simulated by NAMI DANCE model in GPU environment. The simulation results suggest that this scenario seems capable of generating such a tsunami observed along the coast of Sunda Strait. However, the contribution of the possible submarine mass movements in the close area between GAK and the surrounding islands either to this event or potential tsunami threat in the region is still questionable. We employed a bathymetric dataset through pre- and post-event analyses, which demonstrate submarine slope failures in the southwestern proximity of GAK. Hence, additional two scenarios of elliptical landslide sources on the slopes of bathymetry change area (could be triggered by seismic motion/volcanic eruption) are considered, searching for the possible effects of the tsunami that might be generated by these submarine landslides. The study may also provide some perspective for potential tsunami generation by combined sources and help to elucidate the extent of volcanic tsunami hazard in the region due to potential future eruptions of Gunung Anak Krakatau.

The present paper gives a partial answer to Smale’s question which diagrams can correspond to (A,B)-diffeomorphisms. Model diffeomorphisms of the two-dimensional torus derived by “Smale surgery” are considered, and necessary and sufficient conditions for their topological conjugacy are found. Also, a class G of (A,B)-diffeomorphisms on surfaces which are the connected sum of the model diffeomorphisms is introduced. Diffeomorphisms of the class G realize any connected Hasse diagrams (abstract Smale graph). Examples of diffeomorphisms from G with isomorphic labeled Smale diagrams which are not ambiently Ω-conjugated are constructed. Moreover, a subset G∗ ⊂ G of diffeomorphisms for which the isomorphism class of labeled Smale diagrams is a complete invariant of the ambient Ω-conjugacy is singled out.

The  existence problem for attractors of foliations with transverse linear connection is investigated. In general foliations with transverse linear connection do not admit attractors. A conditions that implies the existence of a global attractor which is a minimal set, is specified. An application to transversely similar pseudo-Riemannian foliations is obtained. The global structure of transversely similar Riemannian foliations is described. Different examples are constructed.

It is a well known fact that any smooth manifold admits a Morse function, whereas the problem of existence of a Morse function for a topological manifold stated by Marston Morse in 1959  is still open. In the present paper we prove that a topological manifold admits a continuous Morse function if it admits a topological flow with a finite hyperbolic chain recurrent set. We construct this function as a Lyapunov function whose set of the critical points coincides with the chain recurrent set of the flow.

We consider a class of gradient-like flows on three-dimensional closedmanifolds whose attractors and repellers belongs to a finite union of embedded surfaces and find conditions when the ambient manifold is Seifert.

The article is devoted to interrelations between an existence of trivial and nontrivial basic sets of A-diffeomorphisms of surfaces. We prove that if all trivial basic sets of a structurally stable diffeomorphism of surface $M^2$ are source periodic points $\alpha_1, …, \alpha_k$, then the non-wandering set of this diffeomorphism consists of points $\alpha_1, …, \alpha_k$ and exactly one one-dimensional attractor $\Lambda$. We give some sufficient conditions for attractor $\Lambda$ to be widely situated. Also, we prove that if a non-wandering set of a structurally stable diffeomorphism contains a nontrivial zero-dimensional basic set, then it also contains source and sink periodic points.

For a wide class of dynamical systems known as Pixton diffeomorphisms the topological conjugacy class is completely defined by the Hopf knot equivalence class, i.e. the knot whose equivalence class under homotopy of the loops is a generator of the fundamental group π1(S2×S1)π1(S2×S1). Moreover, any Hopf knot can be realized by a Pixton diffeomorphism. Nevertheless, the number of the classes of topological conjugacy of these diffeomorphisms is still unknown. This problem can be reduced to finding topological invariants of Hopf knots. In the present paper we describe a first order invariant for these knots. This result allows one to model countable families of pairwise non-equivalent Hopf knots and, therefore, infinite set of topologically non-conjugate Pixton diffeomorphisms.

The observation of a wave group persisting for more than 200 periods in the direct numerical simulation of nonlinear unidirectional irregular water waves in deep water is discussed. The simulation conditions are characterized by parameters realistic for broad-banded waves in the sea. Through solution of the associated scattering problem for the nonlinear Schr€odinger equation, the group is identified as the intense envelope soliton with remarkably stable parameters. Most of the extreme waves occur on top of this group, resulting in higher and longer rogue wave events.

For area-preserving H'enon-like maps and their compositions, we consider smooth perturbations that keep the reversibility of the initial maps but destroy their conservativity. For constructing such perturbations, we use two methods, a new method based on reversible properties of maps written in the so-called cross-form, and the classical Quispel-Roberts method based on a variation of involutions of the initial map. We study symmetry breaking bifurcations of symmetric periodic orbits in reversible families containing quadratic conservative orientable and nonorientable H'enon maps as well as a product of two H'enon maps whose Jacobians are mutually inverse.

An analytic reversible Hamiltonian system with two degrees of freedom is studied in a neighborhood of its symmetric heteroclinic connection made up of a symmetric saddle-center, a symmetric orientable saddle periodic orbit lying in the same level of a Hamiltonian, and two non-symmetric heteroclinic orbits permuted by the involution. This is a codimension one structure; therefore, it can be met generally in oneparameter families of reversible Hamiltonian systems. There exist two possible types of such connections depending on how the involution acts near the equilibrium. We prove a series of theorems that show a chaotic behavior of the system and those in its unfoldings, in particular, the existence of countable sets of transverse homoclinic orbits to the saddle periodic orbit in the critical level, transverse heteroclinic connections involving a pair of saddle periodic orbits, families of elliptic periodic orbits, homoclinic tangencies, families of homoclinic orbits to saddle-centers in the unfolding, etc. As a by-product, we get a criterion of the existence of homoclinic orbits to a saddle-center.  

A Shilnikov homoclinic attractor of a three-dimensional diffeomorphism contains a saddle-focus fixed point with a two-dimensional unstable invariant manifold and homoclinic orbits to this saddle-focus. The orientation-reversing property of the diffeomorphism implies a symmetry between two branches of the one-dimensional stable manifold. This symmetry leads to a significant difference between Shilnikov attractors in the orientation-reversing and orientation-preserving cases. We consider the three-dimensional Mirá map 𝑥¯=𝑦,𝑦¯=𝑧, and 𝑧¯=𝐵𝑥+𝐶𝑦+𝐴𝑧−𝑦2x¯=y,y¯=z, and z¯=Bx+Cy+Az−y2 with the negative Jacobian (𝐵<0B<0) as a basic model demonstrating various types of Shilnikov attractors. We show that depending on values of parameters 𝐴,𝐵A,B, and 𝐶C, such attractors can be of three possible types: hyperchaotic (with two positive and one negative Lyapunov exponent), flow-like (with one positive, one very close to zero, and one negative Lyapunov exponent), and strongly dissipative (with one positive and two negative Lyapunov exponents). We study scenarios of the formation of such attractors in one-parameter families.