All issues

This article is devoted to some problems of declaring and storing information for objects' visualization. The storage structure and resources control technology can be applied for real-time visualization of three-dimensional scenes. Such instruments as Sample Framework from DirectX SDK and Direct3D Extension Library (D3DX) were used in the implementation.

In this paper, we consider a quartic family of planar vector fields corresponding to a rational Holling system which models the dynamics of the populations of predators and their prey in a given ecological or biomedical system and which is a variation on the classical Lotka–Volterra system. For the latter system, the change of the prey density per unit of time per predator called the response function is proportional to the prey density. This means that there is no saturation of the predator when the amount of available prey is large. However, it is more realistic to consider a nonlinear and bounded response function, and in fact different response functions have been used in the literature to model the predator response. After algebraic transformations, the rational Holling system can be written in the form of a quartic dynamical system. To investigate the character and distribution of the singular points in the phase plane of the quartic system, we use our method the sense of which is to obtain the simplest (well-known) system by vanishing some parameters (usually field rotation parameters) of the original system and then to input these parameters successively one by one studying the dynamics of the singular points (both finite and infinite) in the phase plane. Using the obtained information on singular points and applying our geometric approach to the qualitative analysis, we study the limit cycle bifurcations of the quartic system. To control all of the limit cycle bifurcations, especially, bifurcations of multiple limit cycles, it is necessary to know the properties and combine the effects of all of the rotation parameters. It can be done by means of the Wintner–Perko termination principle stating that the maximal one-parameter family of multiple limit cycles terminates either at a singular point which is typically of the same multiplicity (cyclicity) or on a separatrix cycle which is also typically of the same multiplicity (cyclicity). Applying this principle, we prove that the quartic system (and the corresponding rational Holling system) can have at most two limit cycles surrounding one singular point.

This work is devoted to development of an algorithm for numerical integration of differential equations potentially-streaming method simulation of non-equilibrium processes. This method was developed by the author in his earlier published works. In this paper, consideration is limited to systems with lumped parameters. Also previously developed method for analyzing the correctness of the author of the approximate solution of the system potentially-streaming equations for systems in lumped parameters. The purpose of this article is to combine this technique with modern numerical methods for integrating systems of ordinary differential equations and the development of methods of numerical integration of systems of equations potentially-streaming method that allows to guarantee the correctness of the approximate solution.

The mathematical model of the magnetic memory cell MRAM with the in-plane anisotropy axis parallel to the edge of a free ferromagnetic layer (longitudinal anisotropy) has been constructed using approximation of uniform magnetization. The model is based on the Landau–Lifshits–Gilbert equation with the injection-current term in the Sloncžewski–Berger form. The set of ordinary differential equations for magnetization dynamics in a three-layered Co/Cu/Cu valve under the control of external magnetic field and spin-polarized current has been derived in the normal coordinate form. It was shown that the set of equations has two main stationary points on the anisotropy axis at any values of field and current. The stationary analysis of them has been performed. The algebraic equations for determination of additional stationary points have been derived. It has been shown that, depending on the field and current magnitude, the set of equations can have altogether two, four, or six stationary points symmetric in pairs relatively the anisotropy axis. The bifurcation diagrams for all the points have been constructed. The classification of the corresponding phase portraits has been performed. The typical trajectories were calculated numerically using Runge–Kutta method. The regions, where stable and unstable limit cycles exist, have been determined. It was found that the unstable limit cycles exist around the main stable equilibrium point on the axis that coincides with the anisotropy one, whereas the stable cycles surround the unstable additional points of equilibrium. The area of their existence was determined numerically. The new types of dynamics, such as accidental switching and non-complete switching, have been found. The threshold values of switching current and field have been obtained analytically. The estimations of switching times have been performed numerically.

In the last decades, universal scenarios of the transition to chaos in dynamic systems have been well studied. The scenario of the transition to chaos is defined as a sequence of bifurcations that occur in the system under the variation one of the governing parameters and lead to a qualitative change in dynamics, starting from the regular mode and ending with chaotic behavior. Typical scenarios include a cascade of period doubling bifurcations (Feigenbaum scenario), the breakup of a low-dimensional torus (Ruelle–Takens scenario), and the transition to chaos through the intermittency (Pomeau–Manneville scenario). In more complicated spatially distributed dynamic systems, the complexity of dynamic behavior growing with a parameter change is closely intertwined with the formation of spatial structures. However, the question of whether the spatial and temporal axes could completely exchange roles in some scenario still remains open. In this paper, for the first time, we propose a mathematical model of convection–diffusion–reaction, in which a spatial transition to chaos through the breakup of the quasi–periodic regime is realized in the framework of the Ruelle–Takens scenario. The physical system under consideration consists of two aqueous solutions of acid (A) and base (B), initially separated in space and placed in a vertically oriented Hele–Shaw cell subject to the gravity field. When the solutions are brought into contact, the frontal neutralization reaction of the second order A + B $\to$ C begins, which is accompanied by the production of salt (C). The process is characterized by a strong dependence of the diffusion coefficients of the reagents on their concentration, which leads to the appearance of two local zones of reduced density, in which chemoconvective fluid motions develop independently. Although the layers, in which convection develops, all the time remain separated by the interlayer of motionless fluid, they can influence each other via a diffusion of reagents through this interlayer. The emerging chemoconvective structure is the modulated standing wave that gradually breaks down over time, repeating the sequence of the bifurcation chain of the Ruelle–Takens scenario. We show that during the evolution of the system one of the spatial axes, directed along the reaction front, plays the role of time, and time itself starts to play the role of a control parameter.

In this paper, using our bifurcation-geometric approach, we study global dynamics and solve the problem of the maximum number and distribution of limit cycles (self-oscillating regimes corresponding to states of dynamical equilibrium) in a planar polynomial mechanical system of the Euler–Lagrange–Liйnard type. Such systems are also used to model electrical, ecological, biomedical and other systems, which greatly facilitates the study of the corresponding real processes and systems with complex internal dynamics. They are used, in particular, in mechanical systems with damping and stiffness. There are a number of examples of technical systems that are described using quadratic damping in second-order dynamical models. In robotics, for example, quadratic damping appears in direct-coupled control and in nonlinear devices, such as variable impedance (resistance) actuators. Variable impedance actuators are of particular interest to collaborative robotics. To study the character and location of singular points in the phase plane of the Euler–Lagrange–Liйnard polynomial system, we use our method the meaning of which is to obtain the simplest (well-known) system by vanishing some parameters (usually, field rotation parameters) of the original system and then to enter sequentially these parameters studying the dynamics of singular points in the phase plane. To study the singular points of the system, we use the classical Poincarй index theorems, as well as our original geometric approach based on the application of the Erugin twoisocline method which is especially effective in the study of infinite singularities. Using the obtained information on the singular points and applying canonical systems with field rotation parameters, as well as using the geometric properties of the spirals filling the internal and external regions of the limit cycles and applying our geometric approach to qualitative analysis, we study limit cycle bifurcations of the system under consideration.

This article presents the results of an analytical and computer study of the collective dynamic properties of a chain of self-oscillating systems (conditionally — oscillators). It is assumed that the couplings of individual elements of the chain are non-reciprocal, unidirectional. More precisely, it is assumed that each element of the chain is under the influence of the previous one, while the reverse reaction is absent (physically insignificant). This is the main feature of the chain. This system can be interpreted as an active discrete medium with unidirectional transfer, in particular, the transfer of a matter. Such chains can represent mathematical models of real systems having a lattice structure that occur in various fields of natural science and technology: physics, chemistry, biology, radio engineering, economics, etc. They can also represent models of technological and computational processes. Nonlinear self-oscillating systems (conditionally, oscillators) with a wide “spectrum” of potentially possible individual self-oscillations, from periodic to chaotic, were chosen as the “elements” of the lattice. This allows one to explore various dynamic modes of the chain from regular to chaotic, changing the parameters of the elements and not changing the nature of the elements themselves. The joint application of qualitative methods of the theory of dynamical systems and qualitative-numerical methods allows one to obtain a clear picture of all possible dynamic regimes of the chain. The conditions for the existence and stability of spatially-homogeneous dynamic regimes (deterministic and chaotic) of the chain are studied. The analytical results are illustrated by a numerical experiment. The dynamical regimes of the chain are studied under perturbations of parameters at its boundary. The possibility of controlling the dynamic regimes of the chain by turning on the necessary perturbation at the boundary is shown. Various cases of the dynamics of chains comprised of inhomogeneous (different in their parameters) elements are considered. The global chaotic synchronization (of all oscillators in the chain) is studied analytically and numerically.

We developed an algorithm for fast simulation of VLSI CMOS (Very Large Scale Integration with Complementary Metal-Oxide-Semiconductors) with an accuracy control. The algorithm provides an ability of parallel numerical experiments in multiprocessor computational environment. There is computation speed up by means of block-matrix and structural (DCCC) decompositions application. A feature of the approach is both in a choice of moments and ways of parameters synchronization and application of multi-rate integration methods. Due to this fact we have ability to estimate and control error of given characteristics.

The algorithm of fuzzy control system essentially nonlinear dynamic object is considered in this article. For solving nonlinear optimal control problem is proposed to use the method of linear quadratic regulation (LQR) with fuzzy Takagi–Sugeno model. The algorithm can be used for the design of deterministic optimal control of nonlinear objects. The algorithm of optimal control for controlling the rotational motion of a space vehicle is proposed.

The mathematical model of a three-layered Co/Cu/Co nanopillar for MRAM cell with one fixed and one free layer was investigated in the approximation of uniformly distributed magnetization. The anisotropy axis is perpendicular to the layers (so-called perpendicular anisotropy). Initially the magnetization of the free layer is oriented along the anisotropy axis in the position accepted to be “zero”. Simultaneous magnetic field and spinpolarized current engaging can reorient the magnetization to another position which in this context can be accepted as “one”. The mathematical description of the effect is based on the classical vector Landau–Lifshits equation with the dissipative term in the Gilbert form. In our model we took into account the interactions of the magnetization with an external magnetic field and such effective magnetic fields as an anisotropy and demagnetization ones. The influence of the spin-polarized injection current is taken into account in the form of Sloczewski–Berger term. The model was reduced to the set of three ordinary differential equations with the first integral. It was shown that at any current and field the dynamical system has two main equilibrium states on the axis coincident with anisotropy axis. It was ascertained that in contrast with the longitudinal-anisotropy model, in the model with perpendicular anisotropy there are no other equilibrium states. The stability analysis of the main equilibrium states was performed. The bifurcation diagrams characterizing the magnetization dynamics at different values of the control parameters were built. The classification of the phase portraits on the unit sphere was performed. The features of the dynamics at different values of the parameters were studied and the conditions of the magnetization reorientation were determined. The trajectories of magnetization switching were calculated numerically using the Runge–Kutta method. The parameter values at which limit cycles exist were determined. The threshold values for the switching current were found analytically. The threshold values for the structures with longitudinal and perpendicular anisotropy were compared. It was established that in the structure with the perpendicular anisotropy at zero field the switching current is an order lower than in the structure with the longitudinal one.