All issues

Theoretical and experimental investigations of water vapor interaction with porous materials are carried out both at the macro level and at the micro level. At the macro level, the influence of the arrangement structure of individual pores on the processes of water vapor interaction with porous material as a continuous medium is studied. At the micro level, it is very interesting to investigate the dependence of the characteristics of the water vapor interaction with porous media on the geometry and dimensions of the individual pore.

In this paper, a study was carried out by means of mathematical modelling of the processes of water vapor interaction with suffering pore of the cylindrical type. The calculations were performed using a model of a hybrid type combining a molecular-dynamic and a macro-diffusion approach for describing water vapor interaction with an individual pore. The processes of evolution to the state of thermodynamic equilibrium of macroscopic characteristics of the system such as temperature, density, and pressure, depending on external conditions with respect to pore, were explored. The dependence of the evolution parameters on the distribution of the diffusion coefficient in the pore, obtained as a result of molecular dynamics modelling, is examined. The relevance of these studies is due to the fact that all methods and programs used for the modelling of the moisture and heat conductivity are based on the use of transport equations in a porous material as a continuous medium with known values of the transport coefficients, which are usually obtained experimentally.

In this paper, we address the mathematical modeling of robots based on tensegrity structures. The pivotal property of such structures is the forming elements working only for compression or tension, which allows the use of materials and structural solutions that minimize the weight of the structure while maintaining its strength.

Tensegrity structures hold several properties important for collaborative robotics, exploration and motion tasks in non-deterministic environments: natural compliance, compactness for transportation, low weight with significant impact resistance and rigidity. The control of such structures remains an open research problem, which is associated with the complexity of describing the dynamics of such structures.

We formulate an approach for describing the dynamics of such structures, based on second-order dynamics of the Cartesian coordinates of structure elements (rods), first-order dynamics for angular velocities of rods, and first-order dynamics for quaternions that are used to describe the orientation of rods. We propose a numerical method for solving these dynamic equations. The proposed methods are implemented in the form of a freely distributed mathematical package with open source code.

Further, we show how the provided software package can be used for modeling the dynamics and determining the operating modes of tensegrity structures. We present an example of a tensegrity structure moving in zero gravity with three rigid rods and nine elastic elements working in tension (cables), showing the features of the dynamics of the structure in reaching the equilibrium position. The range of initial conditions for which the structure operates in the normal mode is determined. The results can be directly used to analyze the nature of passive dynamic movements of the robots based on a three-link tensegrity structure, considered in the paper; the proposed modeling methods and the developed software are suitable for modeling a significant variety of tensegrity robots.

Changes in abundance for emerging populations can develop according to several dynamic scenarios. After rapid biological invasions, the time factor for the development of a reaction from the biotic environment will become important. There are two classic experiments known in history with different endings of the confrontation of biological species. In Gause’s experiments with ciliates, the infused predator, after brief oscillations, completely destroyed its resource, so its $r$-parameter became excessive for new conditions. Its own reproductive activity was not regulated by additional factors and, as a result, became critical for the invader. In the experiments of the entomologist Uchida with parasitic wasps and their prey beetles, all species coexisted. In a situation where a population with a high reproductive potential is regulated by several natural enemies, interesting dynamic effects can occur that have been observed in phytophages in an evergreen forest in Australia. The competing parasitic hymenoptera create a delayed regulation system for rapidly multiplying psyllid pests, where a rapid increase in the psyllid population is allowed until the pest reaches its maximum number. A short maximum is followed by a rapid decline in numbers, but minimization does not become critical for the population. The paper proposes a phenomenological model based on a differential equation with a delay, which describes a scenario of adaptive regulation for a population with a high reproductive potential with an active, but with a delayed reaction with a threshold regulation of exposure. It is shown that the complication of the regulation function of biotic resistance in the model leads to the stabilization of the dynamics after the passage of the minimum number by the rapidly breeding species. For a flexible system, transitional regimes of growth and crisis lead to the search for a new equilibrium in the evolutionary confrontation.

Based on the time-discrete model, we study the effect of selective proportional harvesting on the population dynamics with age and sex structure. When constructing the model, we assume that the population birth rate depends on the ratio of the sexes and the number of formed pairs. The regulation of population growth is carried out by limiting the juvenile’s survival when the survival of immature individuals decreases with an increase in the numbers of sex and age classes. We consider cases where the harvest is carried out only from a younger age class or from a group of mature females or males. We find that the harvesting of males or females at the optimal level is responsible for changing the ratio of females to males (taking into account the average size of the harem). We show that the maximum number of harvested males is achieved either at such a harvest rate when their excess number is withdrawn and the balance of sexes is established or at such an optimal catch quota at which the sex ratio is shifted towards breeding females. Optimal female harvesting, in which the highest number of them are taken, either maintains a preexisting shortage of adult males or leads to an excess of males or the fixing of a sex balance. We find that, depending on the population parameters for all considered harvesting strategies, the hydra effect can observe, i. e., the equilibrium size of the exploited sex and age-specific group (after reproduction) can increase with the growth of harvesting intensity. The selective harvesting, due to which the hydra effect occurs, simultaneously leads to an increase remaining population size and the number of harvested individuals. At the same time, the size of the exploited group after reproduction can become even more than without exploitation. Equilibrium harvesting with the optimal harvest rate that maximizes yield leads to a population size decrease. The effect of hydra is at lower values of the catch quota than the optimal harvest rate. At the same time, the consequence of the hydra effect may be a higher abundance of the age-sex group under optimal exploitation compared to the level observed in the absence of harvesting.

The model of market pricing in duopoly describing the prices dynamics as a two-dimensional map is presented. It is shown that the fixed point of the map coincides with the local Nash-equilibrium price in duopoly game. There have been numerically identified a bifurcation of the fixed point, shown the scheme of transition from periodic to chaotic mode through a doubling period. To ensure the sustainability of local Nashequilibrium price the controlling chaos mechanism has been proposed. This mechanism allows to harmonize the economic interests of the firms and to form the balanced pricing policy.

Transport of respiratory gases by respiratory and circulatory systems is one of the most important processes associated with living conditions of the human body. Significant and/or long-term deviations of oxygen and carbon dioxide concentrations from the normal values in blood can be a reason of significant pathological changes with irreversible consequences: lack of oxygen (hypoxia and ischemic events), the change in the acidbase balance of blood (acidosis or alkalosis), and others. In the context of a changing external environment and internal conditions of the body the action of its regulatory systems aimed at maintaining homeostasis. One of the major mechanisms for maintaining concentrations (partial pressures) of oxygen and carbon dioxide in the blood at a normal level is the regulation of minute ventilation, respiratory rate and depth of respiration, which is caused by the activity of the central and peripheral regulators.

In this paper we propose a mathematical model of the regulation of pulmonary ventilation parameter. The model is used to calculate the minute ventilation adaptation during hypoxia and hypercapnia. The model is developed using a single-component model of the lungs, and biochemical equilibrium conditions of oxygen and carbon dioxide in the blood and the alveolar lung volume. A comparison with laboratory data is performed during hypoxia and hypercapnia. Analysis of the results shows that the model reproduces the dynamics of minute ventilation during hypercapnia with sufficient accuracy. Another result is that more accurate model of regulation of minute ventilation during hypoxia should be developed. The factors preventing from satisfactory accuracy are analysed in the final section.

Respiratory function is one of the main limiting factors of the organism during intense physical activities. Thus, it is important characteristic of high performance sport and extreme physical activity conditions. Therefore, the results of this study have significant application value in the field of mathematical modeling in sport. The considered conditions of hypoxia and hypercapnia are partly reproduce training at high altitude and at hypoxia conditions. The purpose of these conditions is to increase the level of hemoglobin in the blood of highly qualified athletes. These conditions are the only admitted by sport committees.

We study the stochastic dynamics of the two-dimensional Hindmarsh–Rose model in the parametrical zone of coexisting stable equilibria and limit cycles. The phenomenon of noise-induced transitions between the attractors is investigated. Under the random disturbances, equilibrium and periodic regimes combine in bursting regime: the system demonstrates an alternation of small fluctuations near the equilibrium with high amplitude oscillations. This effect is analysed using the stochastic sensitivity function technique and a method of estimation of critical values for noise intensity is proposed.

Execution or violation of the principle of competitive exclusion in communities is the subject of many studies. The principle of competitive exclusion means that coexistence of species in community is impossible if the number of species exceeds the number of controlling mutually independent factors. At that time there are many examples displaying the violations of this principle in the natural systems. The explanations for this paradox vary from inexact identification of the set of factors to various types of spatial and temporal heterogeneities. One of the factors breaking the principle of competitive exclusion is intraspecific competition. This study holds the model of community with two species and one influencing factor with density-dependent mortality and spatial heterogeneity. For such models possibility of the existence of stable equilibrium is proved in case of spatial homogeneity and negative effect of the species on the factor. Our purpose is analysis of possible variants of dynamics of the system with spatial heterogeneity under the various directions of the species effect on the influencing factor. Numerical analysis showed that there is stable coexistence of the species agreed with homogenous spatial distributions of the species if the species effects on the influencing factor are negative. Density-dependent mortality and spatial heterogeneity lead to violation of the principle of competitive exclusion when equilibriums are Turing unstable. In this case stable spatial heterogeneous patterns can arise. It is shown that Turing instability is possible if at least one of the species effects is positive. Model nonlinearity and spatial heterogeneity cause violation of the principle of competitive exclusion in terms of both stable spatial homogenous states and quasistable spatial heterogeneous patterns.

We consider “predator – prey” model taking into account the competition of prey, predator for different from the prey resources, and their interaction described by the second type Holling trophic function. An analysis of the attractors is carried out depending on the coefficient of competition of predators. In the deterministic case, this model demonstrates the complex behavior associated with the local (Andronov –Hopf and saddlenode) and global (birth of a cycle from a separatrix loop) bifurcations. An important feature of this model is the disappearance of a stable cycle due to a saddle-node bifurcation. As a result of the presence of competition in both populations, parametric zones of mono- and bistability are observed. In parametric zones of bistability the system has either coexisting two equilibria or a cycle and equilibrium. Here, we investigate the geometrical arrangement of attractors and separatrices, which is the boundary of basins of attraction. Such a study is an important component in understanding of stochastic phenomena. In this model, the combination of the nonlinearity and random perturbations leads to the appearance of new phenomena with no analogues in the deterministic case, such as noise-induced transitions through the separatrix, stochastic excitability, and generation of mixed-mode oscillations. For the parametric study of these phenomena, we use the stochastic sensitivity function technique and the confidence domain method. In the bistability zones, we study the deformations of the equilibrium or oscillation regimes under stochastic perturbation. The geometric criterion for the occurrence of such qualitative changes is the intersection of confidence domains and the separatrix of the deterministic model. In the zone of monostability, we evolve the phenomena of explosive change in the size of population as well as extinction of one or both populations with minor changes in external conditions. With the help of the confidence domains method, we solve the problem of estimating the proximity of a stochastic population to dangerous boundaries, upon reaching which the coexistence of populations is destroyed and their extinction is observed.

The numerical methods developed by the author recently for calculating the molecular system based on the direct solution of the Schrodinger equation by the Monte Carlo method have shown a huge uncertainty in the choice of solutions. On the one hand, it turned out to be possible to build many new solutions; on the other hand, the problem of their connection with reality has become sharply aggravated. In ab initio quantum mechanical calculations, the problem of choosing solutions is not so acute after the transition to the classical format of describing a molecular system in terms of potential energy, the method of molecular dynamics, etc. In this paper, we investigate the problem of choosing solutions in the classical format of describing a molecular system without taking into account quantum mechanical prerequisites. As it turned out, the problem of choosing solutions in the classical format of describing a molecular system is reduced to a specific marking of the configuration space in the form of a set of stationary points and reconstruction of the corresponding potential energy function. In this formulation, the solution of the choice problem is reduced to two possible physical and mathematical problems: to find all its stationary points for a given potential energy function (the direct problem of the choice problem), to reconstruct the potential energy function for a given set of stationary points (the inverse problem of the choice problem). In this paper, using a computational experiment, the direct problem of the choice problem is discussed using the example of a description of a monoatomic cluster. The number and shape of the locally equilibrium (saddle) configurations of the binary potential are numerically estimated. An appropriate measure is introduced to distinguish configurations in space. The format of constructing the entire chain of multiparticle contributions to the potential energy function is proposed: binary, threeparticle, etc., multiparticle potential of maximum partiality. An infinite number of locally equilibrium (saddle) configurations for the maximum multiparticle potential is discussed and illustrated. A method of variation of the number of stationary points by combining multiparticle contributions to the potential energy function is proposed. The results of the work listed above are aimed at reducing the huge arbitrariness of the choice of the form of potential that is currently taking place. Reducing the arbitrariness of choice is expressed in the fact that the available knowledge about the set of a very specific set of stationary points is consistent with the corresponding form of the potential energy function.