All issues

The article discusses the emergence of chaotic attractors in the system of three ordinary differential equations arising in the theory of «reaction-diffusion» systems. The dynamics of the corresponding one- and two-dimensional maps and Lyapunov exponents of such attractors are studied. It is shown that the transition to chaos is in accordance with a non-traditional scenario of repeated birth and disappearance of chaotic regimes, which had been previously studied for one-dimensional maps with a sharp apex and a quadratic minimum. Some characteristic features of the system — zones of bistability and hyperbolicity, the crisis of chaotic attractors — are studied by means of numerical analysis.

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.

Semiclassical approximation formalism is developed for the multidimensional Fokker–Planck–Kolmogorov equation with non-local and nonlinear drift vector with respect to a small diffusion coefficient D, D→0, in the class of trajectory concentrated functions. The Einstein−Ehrenfest system of (0, M)-type is obtained. A family of semiclassical solutions localized around a point driven by the Einstein−Ehrenfest system accurate to O(D^(M+1)/2) is found.

Investigation of occurrence of diffusion instability in a set of three reaction–diffusion equations is carried out. In the general case the condition for both Turing and wave instabilities are obtained. Qualitative properties of the system, in which the bifurcation of each of the two types can take place, are clarified. In numerical experiments it is shown that if the corresponding conditions are met in the nonlinear model, spatiotemporal patterns are formed, which are predicted by linear analysis.

The numerical solving of the system of high-temperature radiative gas dynamics (HTRGD) equations is a computationally laborious task, since the interaction of radiation with matter is nonlinear and non-local. The radiation absorption coefficients depend on temperature, and the temperature field is determined by both gas-dynamic processes and radiation transport. The method of splitting into physical processes is usually used to solve the HTRGD system, one of the blocks consists of a joint solving of the radiative transport equation and the energy balance equation of matter under known pressure and temperature fields. Usually difference schemes with orders of convergence no higher than the second are used to solve this block. Due to computer memory limitations it is necessary to use not too detailed grids to solve complex technical problems. This increases the requirements for the order of approximation of difference schemes. In this work, bicompact schemes of a high order of approximation for the algorithm for the joint solution of the radiative transport equation and the energy balance equation are implemented for the first time. The proposed method can be applied to solve a wide range of practical problems, as it has high accuracy and it is suitable for solving problems with coefficient discontinuities. The non-linearity of the problem and the use of an implicit scheme lead to an iterative process that may slowly converge. In this paper, we use a multiplicative HOLO algorithm named the quasi-diffusion method by V.Ya.Goldin. The key idea of HOLO algorithms is the joint solving of high order (HO) and low order (LO) equations. The high-order equation (HO) is the radiative transport equation solved in the energy multigroup approximation, the system of quasi-diffusion equations in the multigroup approximation (LO₁) is obtained by averaging HO equations over the angular variable. The next step is averaging over energy, resulting in an effective one-group system of quasi-diffusion equations (LO₂), which is solved jointly with the energy equation. The solutions obtained at each stage of the HOLO algorithm are closely related that ultimately leads to an acceleration of the convergence of the iterative process. Difference schemes constructed by the method of lines within one cell are proposed for each of the stages of the HOLO algorithm. The schemes have the fourth order of approximation in space and the third order of approximation in time. Schemes for the transport equation were developed by B.V. Rogov and his colleagues, the schemes for the LO₁ and LO₂ equations were developed by the authors. An analytical test is constructed to demonstrate the declared orders of convergence. Various options for setting boundary conditions are considered and their influence on the order of convergence in time and space is studied.

Numerical solutions for the two-dimensional reaction-diffusion equation with nonlocal nonlinearity are obtained. The solutions reveal formation of dissipative structures. Structures arising from initial distributions with one and several centers of localization are considered. Formation of extending circular structures is shown. Peculiarities of formation and interaction of extending circular structures depending on  nonlocal interaction are considered.

Optical systems with two-dimensional feedback demonstrate wide possibilities for studying the nucleation and development processes of dissipative structures. Feedback allows to influence the dynamics of the optical system by controlling the transformation of spatial variables performed by prisms, lenses, dynamic holograms and other devices. A nonlinear interferometer with a mirror image of a field in two-dimensional feedback is one of the simplest optical systems in which is realized the nonlocal nature of light fields.

A mathematical model of optical systems with two-dimensional feedback is a nonlinear parabolic equation with rotation transformation of a spatial variable and periodicity conditions on a circle. Such problems are investigated: bifurcation of the traveling wave type stationary structures, how the form of the solution changes as the diffusion coefficient decreases, dynamics of the solution’s stability when the bifurcation parameter leaves the critical value. For the first time as a parameter bifurcation was taken of diffusion coefficient.

The method of central manifolds and the Galerkin’s method are used in this paper. The method of central manifolds and the Galerkin’s method are used in this paper. The method of central manifolds allows to prove a theorem on the existence and form of the traveling wave type solution neighborhood of the bifurcation value. The first traveling wave born as a result of the Andronov –Hopf bifurcation in the transition of the bifurcation parameter through the сritical value. According to the central manifold theorem, the first traveling wave is born orbitally stable.

Since the above theorem gives the opportunity to explore solutions are born only in the vicinity of the critical values of the bifurcation parameter, the decision to study the dynamics of traveling waves of change during the withdrawal of the bifurcation parameter in the supercritical region, the formalism of the Galerkin method was used. In accordance with the method of the central manifold is made Galerkin’s approximation of the problem solution. As the bifurcation parameter decreases and its transition through the critical value, the zero solution of the problem loses stability in an oscillatory manner. As a result, a periodic solution of the traveling wave type branches off from the zero solution. This wave is born orbitally stable. With further reduction of the parameter and its passage through the next critical value from the zero solution, the second solution of the traveling wave type is produced as a result of the Andronov –Hopf bifurcation. This wave is born unstable with an instability index of two.

Numerical calculations have shown that the application of the Galerkin’s method leads to correct results. The results obtained are in good agreement with the results obtained by other authors and can be used to establish experiments on the study of phenomena in optical systems with feedback.

A classical for nonlinear dynamics model, Brusselator, is considered, being augmented by addition of a third variable, which plays the role of a fast-diffusing inhibitor. The model is investigated in one-dimensional case in the parametric domain, where two types of diffusive instabilities of system’s homogeneous stationary state are manifested: wave instability, which leads to spontaneous formation of autowaves, and Turing instability, which leads to spontaneous formation of stationary dissipative structures, or Turing structures. It is shown that, due to the subcritical nature of Turing bifurcation, the interaction of two instabilities in this system results in spontaneous formation of stationary dissipative structures already before the passage of Turing bifurcation. In response to different perturbations of spatially uniform stationary state, different stable regimes are manifested in the vicinity of the double bifurcation point in the parametric region under study: both pure regimes, which consist of either stationary or autowave dissipative structures; and mixed regimes, in which different modes dominate in different areas of the computational space. In the considered region of the parametric space, the system is multistable and exhibits high sensitivity to initial noise conditions, which leads to blurring of the boundaries between qualitatively different regimes in the parametric region. At that, even in the area of dominance of mixed modes with prevalence of Turing structures, the establishment of a pure autowave regime has significant probability. In the case of stable mixed regimes, a sufficiently strong local perturbation in the area of the computational space, where autowave mode is manifested, can initiate local formation of new stationary dissipative structures. Local perturbation of the stationary homogeneous state in the parametric region under investidation leads to a qualitatively similar map of established modes, the zone of dominance of pure autowave regimes being expanded with the increase of local perturbation amplitude. In two-dimensional case, mixed regimes turn out to be only transient — upon the appearance of localized Turing structures under the influence of wave regime, they eventually occupy all available space.

In the context of problems of hydrogen and thermonuclear power engineering intensive research of the hydrogen isotopes properties is being conducted. Mathematical models help to specify physical-chemical ideas about the interaction of hydrogen isotopes with structural materials, to discover the limiting factors. Classical diffusion models are often insufficient. The paper is devoted to the models and numerical solution of the boundary-value problems of hydrogen thermodesorption and permeability taking into account nonlinear sorption-desorption dynamics on the surface and reversible capture of hydrogen atoms in the bulk. Algorithms based on difference approximations. The results of computer simulation of the hydrogen flux from a structural material sample are presented.

We now review the main points of mean-field approximation (MFA) in its application to multicomponent stochastic reaction-diffusion systems.

We present the chemical reaction model under study — brusselator. We write the kinetic equations of reaction supplementing them with terms that describe the diffusion of the intermediate components and the fluctuations of the concentrations of the initial products. We simulate the fluctuations as random Gaussian homogeneous and spatially isotropic fields with zero means and spatial correlation functions with a non-trivial structure. The model parameter values correspond to a spatially-inhomogeneous ordered state in the deterministic case.

In the MFA we derive single-site two-dimensional nonlinear self-consistent Fokker–Planck equation in the Stratonovich's interpretation for spatially extended stochastic brusselator, which describes the dynamics of probability distribution density of component concentration values of the system under consideration. We find the noise intensity values appropriate to two types of Fokker–Planck equation solutions: solution with transient bimodality and solution with the multiple alternation of unimodal and bimodal types of probability density. We study numerically the probability density dynamics and time behavior of variances, expectations, and most probable values of component concentrations at various noise intensity values and the bifurcation parameter in the specified region of the problem parameters.

Beginning from some value of external noise intensity inside the ordered phase disorder originates existing for a finite time, and the higher the noise level, the longer this disorder “embryo” lives. The farther away from the bifurcation point, the lower the noise that generates it and the narrower the range of noise intensity values at which the system evolves to the ordered, but already a new statistically steady state. At some second noise intensity value the intermittency of the ordered and disordered phases occurs. The increasing noise intensity leads to the fact that the order and disorder alternate increasingly.

Thus, the scenario of the noise induced order–disorder transition in the system under study consists in the intermittency of the ordered and disordered phases.