All issues

The article focuses on a nonlinear mathematical model that describes the interaction of the different age groups of the employed population. The interactions are treated by analogy with population relationship (competition, discrimination, assistance, oppression, etc). Under interaction of peoples we mean the generalized social and economic mechanisms that cause related changes in the number of employees of different age groups. Three age groups of the employed population are considered. It is young specialists (15–29 years), workers with experience (30–49 years), the employees of pre-retirement and retirement age (50 and older). The estimation of model’s parameters for the southern regions of the Far Eastern Federal District (FEFD) is executed by statistical data. Analysis of model scenarios allows us to conclude the observed number fluctuations of the different ages employees on the background of a stable total employed population may be a consequence of complex interactions between these groups of peoples. Computational experiments with the obtained values of the parameters allowed us to calculate the rate of decline and the aging of the working population and to determine the nature of the interaction between the age groups of employees that are not directly as reflected in the statistics. It was found that in FEFD the employed of 50 years and older are discriminated against by the young workers under 29, employed up to 29 and 30–49 years are in a partnership. It is shown in most developed regions (Primorsky and Khabarovsk Krai) there is “uniform” competition among different age groups of the employed population. For Primorsky Krai we were able to identify the mixing effect dynamics. It is a typical situation for systems in a state of structural adjustment. This effect is reflected in the fact the long cycles of employed population form with a significant decrease in migration inflows of employees 30–49 years. Besides, the change of migration is accompanied by a change of interaction type — from employment discrimination by the oldest of middle generation to discrimination by the middle of older generation. In less developed regions (Amur, Magadan and Jewish Autonomous Regions) there are lower values of migration balance of almost all age groups and discrimination by young workers up 29 years of other age groups and employment discrimination 30–49 years of the older generation.

In the paper the statistical relationships between the size and production characteristics of phytoplankton and zooplankton of the Vistula and Curonian lagoons, the Baltic Sea, were investigated. Research phytoplankton and zooplankton within the Russian part of the area of the Vistula and the Curonian lagoon was carried out on the monthly basis (from April to November) within the framework of long-term monitoring program on evaluating of ecological status of the lagoons. The size structure of plankton is the basis for understanding of the development of production processes, mechanisms of formation of the plankton species diversity and functioning of the lagoon ecosystems. As results of the work it was found that the maximum rate of photosynthesis and the integral value of the primary production with a change in cell volume of phytoplankton are changed according to a power law. The result shows that the smaller the size of algal cells in phytoplankton communities the more actively occur metabolism and the more effective they assimilate the solar energy. It is shown that the formation of plankton species diversity in ecosystems of lagoons is closely linked with the size structure of plankton communities and with features of development of the production processes. It is proposed the structure of a spatially homogenous mathematical model of the plankton food chain for the lagoon ecosystems taking into account the size spectrum and the characteristics of phytoplankton and zooplankton. The model parameters are the sizedependent indicators allometrically linked with average volumes of cells and organisms in different ranges of their sizes. In the model the algorithm for changes over time the coefficients of food preferences in the diet of zooplankton was proposed. Developed the size-dependent mathematical model of aquatic ecosystems allows to consider the impact of turbulent exchange on the size structure and temporal dynamics of the plankton food chain of the Vistula and Curonian lagoons. The model can be used to study the different regimes of dynamic behavior of plankton systems depending on the changes in the values of its parameters and external influences, as well as to quantify the redistribution of matter flows in ecosystems of the lagoons.

In the article is carried out the analysis of historical process with the use of methods of synergetics (science about the nonlinear developing systems in nature and the society), developed in the works of D. S. Chernavskii in connection with to economic and social systems. It is shown that social self-organizing depending on conditions leads to the formation of both the societies with the strong internal competition (Y-structures) and cooperative type societies (X-structures). Y-structures are characteristic for the countries of the West, X-structure are characteristic for the countries of the East. It is shown that in XIX and in XX centuries occurred accelerated shaping and strengthening of Y-structures. However, at present world system entered into the period of serious structural changes in the economic, political, ideological spheres: the domination of Y-structures concludes. Are examined the possible ways of further development of the world system, connected with change in the regimes of self-organizing and limitation of internal competition. This passage will be prolonged and complex. Under these conditions it will objectively grow the value of the civilizational experience of Russia, on basis of which was formed combined type social system. It is shown that ultimately inevitable the passage from the present do-mination of Y-structures to the absolutely new global system, whose stability will be based on the new ideology, the new spirituality (i.e., new “conditional information” according D. S. Chernavskii), which makes a turn from the principles of competition to the principles of collaboration.

We present results of a study aimed at identification of a controlled object’s channels based on postprocessing of measurements with development of a model of a multiple-input controlled object and subsequent active modelling experiment. The controlled object model is developed using approximation of its behavior by a neural network model using trends obtained during a passive experiment in the mode of normal operation. Recurrent neural network containing feedback elements allows to simulate behavior of dynamic objects; input and feedback time delays allow to simulate behavior of inertial objects with pure delay. The model was taught using examples of the object’s operation with a control system and is presented by a dynamic neural network and a model of a regulator with a known regulation function. The neural network model simulates the system’s behavior and is used to conduct active computing experiments. Neural network model allows to obtain the controlled object’s response to an exploratory stimulus, including a periodic one. The obtained complex frequency response is used to evaluate parameters of the object’s transfer system using the least squares method. We present an example of identification of a channel of the simulated control system. The simulated object has two input ports and one output port and varying transport delays in transfer channels. One of the input ports serves as a controlling stimulus, the second is a controlled perturbation. The controlled output value changes as a result of control stimulus produced by the regulator operating according to the proportional-integral regulation law based on deviation of the controlled value from the task. The obtained parameters of the object’s channels’ transfer functions are close to the parameters of the input simulated object. The obtained normalized error of the reaction for a single step-wise stimulus of the control system model developed based on identification of the simulated control system doesn’t exceed 0.08. The considered objects pertain to the class of technological processes with continuous production. Such objects are characteristic of chemical, metallurgic, mine-mill, pulp and paper, and other industries.

In various areas of science, technology, environment protection, construction, it is very important to study processes of porous materials interaction with different substances in different aggregation states. From the point of view of ecology and environmental protection it is particularly actual to investigate processes of porous materials interaction with water in liquid and gaseous phases. Since one mole of water contains 6.022140857 · 10²³ molecules of H₂O, macroscopic approaches considering the water vapor as continuum media in the framework of classical aerodynamics are mainly used to describe properties, for example properties of water vapor in the pore. In this paper we construct and use for simulation the macroscopic two-dimensional diffusion model [Bitsadze, Kalinichenko, 1980] describing the behavior of water vapor inside the isolated pore. Together with the macroscopic model it is proposed microscopic model of the behavior of water vapor inside the isolated pores. This microscopic model is built within the molecular dynamics approach [Gould et al., 2005]. In the microscopic model a description of each water molecule motion is based on Newton classical mechanics considering interactions with other molecules and pore walls. Time evolution of “water vapor – pore” system is explored. Depending on the external to the pore conditions the system evolves to various states of equilibrium, characterized by different values of the macroscopic characteristics such as temperature, density, pressure. Comparisons of results of molecular dynamic simulations with the results of calculations based on the macroscopic diffusion model and experimental data allow to conclude that the combination of macroscopic and microscopic approach could produce more adequate and more accurate description of processes of water vapor interaction with porous materials.

We study an appearance of gravitational convection in a porous medium saturated by the double-diffusive fluid. The rectangle heated from below is considered with anisotropy of media properties. We analyze Darcy – Boussinesq equations for a binary fluid with Soret effect.

Resulting system for the stream function, the deviation of temperature and concentration is cosymmetric under some additional conditions for the parameters of the problem. It means that the quiescent state (mechanical equilibrium) loses its stability and a continuous family of stationary regimes branches off. We derive explicit formulas for the critical values of the Rayleigh numbers both for temperature and concentration under these conditions of the cosymmetry. It allows to analyze monotonic instability of mechanical equilibrium, the results of corresponding computations are presented.

A finite-difference discretization of a second-order accuracy is developed with preserving of the cosymmetry of the underlying system. The derived numerical scheme is applied to analyze the stability of mechanical equilibrium.

The appearance of stationary and nonstationary convective regimes is studied. The neutral stability curves for the mechanical equilibrium are presented. The map for the plane of the Rayleigh numbers (temperature and concentration) are displayed. The impact of the parameters of thermal diffusion on the Rayleigh concentration number is established, at which the oscillating instability precedes the monotonic instability. In the general situation, when the conditions of cosymmetry are not satisfied, the derived formulas of the critical Rayleigh numbers can be used to estimate the thresholds for the convection onset.

We analyze variants of considering the inhomogeneity of the environment in computer modeling of the dynamics of a predator and prey based on a system of reaction-diffusion–advection equations. The local interaction of species (reaction terms) is described by the logistic law for the prey and the Beddington –DeAngelis functional response, special cases of which are the Holling type II functional response and the Arditi – Ginzburg model. We consider a one-dimensional problem in space for a heterogeneous resource (carrying capacity) and three types of taxis (the prey to resource and from the predator, the predator to the prey). An analytical approach is used to study the stability of stationary solutions in the case of local interaction (diffusionless approach). We employ the method of lines to study diffusion and advective processes. A comparison of the critical values of the mortality parameter of predators is given. Analysis showed that at constant coefficients in the Beddington –DeAngelis model, critical values are variable along the spatial coordinate, while we do not observe this effect for the Arditi –Ginzburg model. We propose a modification of the reaction terms, which makes it possible to take into account the heterogeneity of the resource. Numerical results on the dynamics of species for large and small migration coefficients are presented, demonstrating a decrease in the influence of the species of local members on the emerging spatio-temporal distributions of populations. Bifurcation transitions are analyzed when changing the parameters of diffusion–advection and reaction terms.

In a number of fundamental and practical problems, it is necessary to describe the dynamics of the motion of complexshaped particles in a high-speed gas flow. An example is the movement of coal particles behind the front of a strong shock wave during an explosion in a coal mine. The paper is devoted to numerical simulation of the dynamics of translational and rotational motion of a square-shaped body, as an example of a particle of a more complex shape than a round one, in a supersonic flow behind a passing shock wave. The formulation of the problem approximately corresponds to the experiments of Professor V. M. Boiko and Professor S. V. Poplavski (ITAM SB RAS).

Mathematical model is based on the two-dimensional Euler equations, which are solved in a region with varying boundaries. The defining system of equations is integrated using an explicit scheme and the Cartesian grid method which was developed and verified earlier. The computational algorithm at the time integration step includes: determining the step value, calculating the dynamics of the body movement (determining the force and moment acting on the body; determining the linear and angular velocities of the body; calculating the new coordinates of the body), calculating the gas parameters. To calculate numerical fluxes through the edges of the cell intersected by the boundaries of the body, we use a two-wave approximation for solving the Riemann problem and the Steger – Warming scheme.

The movement of a square with a side of 6 mm was initiated by the passage of a shock wave with a Mach number of 3,0 propagating in a flat channel 800 mm long and 60 mm wide. The channel was filled with air at low pressure. Different initial orientation of the square relative to the channel axis was considered. It is found that the initial position of the square with its side across the flow is less stable during its movement than the initial position with a diagonal across the flow. In this case, the calculated results qualitatively correspond to experimental observations. For the intermediate initial positions of a square, a typical mode of its motion is described, consisting of oscillations close to harmonic, turning into rotation with a constant average angular velocity. During the movement of the square, there is an average monotonous decrease in the distance between the center of mass and the center of pressure to zero.

An algorithm is proposed to identify parameters of a 2D vortex structure used on information about the flow velocity at a finite (small) set of reference points. The approach is based on using a set of point vortices as a model system and minimizing a functional that compares the model and known sets of velocity vectors in the space of model parameters. For numerical implementation, the method of gradient descent with step size control, approximation of derivatives by finite differences, and the analytical expression of the velocity field induced by the point vortex model are used. An experimental analysis of the operation of the algorithm on test flows is carried out: one and a system of several point vortices, a Rankine vortex, and a Lamb dipole. According to the velocity fields of test flows, the velocity vectors utilized for identification were arranged in a randomly distributed set of reference points (from 3 to 200 pieces). Using the computations, it was determined that: the algorithm converges to the minimum from a wide range of initial approximations; the algorithm converges in all cases when the reference points are located in areas where the streamlines of the test and model systems are topologically equivalent; if the streamlines of the systems are not topologically equivalent, then the percentage of successful calculations decreases, but convergence can also take place; when the method converges, the coordinates of the vortices of the model system are close to the centers of the vortices of the test configurations, and in many cases, the values of their circulations also; con-vergence depends more on location than on the number of vectors used for identification. The results of the study allow us to recommend the proposed algorithm for identifying 2D vortex structures whose streamlines are topologically close to systems of point vortices.

The paper simulates a mixture of gases in a multi-stage micro-pump and evaluates its effectiveness at separating the components of the mixture. A device in the form of a long channel with a series of transverse plates is considered. A temperature difference between the sides of the plates induces a radiometric gas flow within the device, and the differences in masses of the gases lead to differences in flow velocities and to the separation of the mixture. Modeling is based on the numerical solution of the Boltzmann kinetic equation, for which a splitting scheme is used, i. e., the advection equation and the relaxation problem are solved separately in alternation. The calculation of the collision integral is performed using the conservative projection method. This method ensures the strict fulfillment of the laws of conservation of mass, momentum, and energy, as well as the important asymptotic property of the equality of the integral of the Maxwell function to zero. Explicit first-order and second-order TVD-schemes are used to solve the advection equation. The calculations were performed for a neon-argon mixture using a model of solid spheres with real molecular diameters and masses. Software has been developed to allow calculations on personal computers and cluster systems. The use of parallelization leads to faster computation and constant time per iteration for devices of different sizes, enabling the modeling of large particle systems. It was found that the value of mixture separation, i. e. the ratio of densities at the ends of the device linearly depends on the number of cascades in the device, which makes it possible to estimate separation for multicascade systems, computer modeling of which is impossible. Flows and distributions of gas inside the device during its operation were analyzed. It was demonstrated that devices of this kind with a sufficiently large number of plates are suitable for the separation of gas mixtures, given that they have no moving parts and are quite simple in manufacture and less subject to wear.