All issues

This work is devoted to the numerical study of high-speed mixing layers of compressible flows. The problem under consideration has a wide range of applications in practical tasks and, despite its apparent simplicity, is quite complex in terms of modeling. Because in the mixing layer, as a result of the instability of the tangential discontinuity of velocities, the flow passes from laminar flow to turbulent mode. Therefore, the obtained numerical results of the considered problem strongly depend on the adequacy of the used turbulence models. In the presented work, this problem is studied based on the two-fluid approach to the problem of turbulence. This approach has arisen relatively recently and is developing quite rapidly. The main advantage of the two-fluid approach is that it leads to a closed system of equations, when, as is known, the long-standing Reynolds approach leads to an open system of equations. The paper presents the essence of the two-fluid approach for modeling a turbulent compressible medium and the methodology for numerical implementation of the proposed model. To obtain a stationary solution, the relaxation method and Prandtl boundary layer theory were applied, resulting in a simplified system of equations. In the considered problem, high-speed flows are mixed. Therefore, it is also necessary to model heat transfer, and the pressure cannot be considered constant, as is done for incompressible flows. In the numerical implementation, the convective terms in the hydrodynamic equations were approximated by the upwind scheme with the second order of accuracy in explicit form, and the diffusion terms in the right-hand sides of the equations were approximated by the central difference in implicit form. The sweep method was used to implement the obtained equations. The SIMPLE method was used to correct the velocity through the pressure. The paper investigates a two-liquid turbulence model with different initial flow turbulence intensities. The obtained numerical results showed that good agreement with the known experimental data is observed at the inlet turbulence intensity of $0.1 < I < 1 \%$. Data from known experiments, as well as the results of the $k − kL + J$ and LES models, are presented to demonstrate the effectiveness of the proposed turbulence model. It is demonstrated that the two-liquid model is as accurate as known modern models and more efficient in terms of computing resources.

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.

The work investigates the process of self-consistent relaxation of the region of disturbances created in a rarefied binary low-temperature plasma by a stationary charged ball or cylinder with an absorbing surface. A feature of such problems is their self-consistent kinetic nature, in which it is impossible to separate the processes of transfer in phase space and the formation of an electromagnetic field. A mathematical model is presented that makes it possible to describe and analyze the state of the gas, electric and thermal fields in the vicinity of the body. The multidimensionality of the kinetic formulation creates certain problems in the numerical solution, therefore a curvilinear system of nonholonomic coordinates was selected for the problem, which minimizes its phase space, which contributes to increasing the efficiency of numerical methods. For such coordinates, the form of the Vlasov kinetic equation has been justified and analyzed. To solve it, a variant of the large particle method with a constant form factor was used. The calculations used a moving grid that tracks the displacement of the distribution function carrier in the phase space, which further reduced the volume of the controlled region of the phase space. Key details of the model and numerical method are revealed. The model and the method are implemented as code in the Matlab language. Using the example of solving a problem for a ball, the presence of significant disequilibrium and anisotropy in the particle velocity distribution in the disturbed zone is shown. Based on the calculation results, pictures of the evolution of the structure of the particle distribution function, profiles of the main macroscopic characteristics of the gas — concentration, current, temperature and heat flow, and characteristics of the electric field in the disturbed region are presented. The mechanism of heating of attracted particles in the disturbed zone is established and some important features of the process of formation of heat flow are shown. The results obtained are well explainable from a physical point of view, which confirms the adequacy of the model and the correct operation of the software tool. The creation and testing of a basis for the development in the future of tools for solving more complex problems of modeling the behavior of ionized gases near charged bodies is noted.

The work will be useful to specialists in the field of mathematical modeling, heat and mass transfer processes, lowtemperature plasma physics, postgraduate students and senior students specializing in the indicated areas.

Safe carrying out of emergency recovery operations at emergency offshore gas condensate wells is possible when taking into account the hazardous factors that prevent anti-fontanning measures. One of such factors is the gassiness of the operation zone due to the release from the water column of a large amount of light, as compared to air, natural gas, as well as vapours of heavier components of gas condensate. To estimate the distribution of explosive concentration of petroleum product vapours in the near surface layer of the atmosphere, it is necessary to determine the characteristics of the source of the contamination. Based on the analysis of theoretical works concerning to the formation of the velocity field in the upper layer of the sea as a result of large amounts of gas coming to the surface, an analytical model is proposed to calculate the size of the area in which a significant amount of gas condensate coming to the surface is vaporised during accidents at shallow-water wells. The stationary regime of reservoir fluid flow during fountaining of offshore gas and oil wells with an underwater location of their mouths is considered. A low-parametric model of oil product evaporation from films of different thickness is constructed. It is shown that the size of the zone of intensive evaporation at shallow-water wells is determined by the volume flow of liquid fraction, its fractional composition and selected threshold for estimation of oil product vapour flow into the atmosphere. In the context of this work shallow water wells are wells with gas flow rate from 1 to 20 million cubic meters at sea depths of about 50–200 metres. In this case, the formation fluid jet from the wellhead on the seabed is transformed into a bubble plume, the stratification of the water column, typical for the summer-autumn period, does not limit the plume’s exit to the sea surface, and the velocity of bubble rise allows the gas dissolution process to be disregardded. The analysis was limited to almost calm hydrometeorological conditions. Such conditions are favourable for offshore operations, but unfavourable from the point of view of dispersion of high concentrations of oil product vapours in the near surface layer of the atmosphere. As a result of this work, an analytical dependence for an approximate assessment of the zone of intensive evaporation of gas condensate is proposed.

The article investigates the possibility of detecting traces of hazardous substances (explosives and narcotics) based on the detection of their vapors in the air. The relevance of the study stems from the need to counter terrorist threats and drug trafficking, where identifying even trace amounts of substances is critical. The focus is on mathematical modeling of the evaporation of a thin substance layer from a surface, based on molecular kinetic theory. A universal model is proposed, accounting for the physicochemical properties of substances, ambient temperature, adhesion to the surface, and the initial mass of the layer. Using the Hertz – Knudsen – Langmuir and Clausius – Clapeyron equations, analytical expressions are derived for the complete evaporation time, maximum vapor mass, and process dynamics. A dimensionless parameter, $\gamma$, is identified, determining the limiting conditions for evaporation. It is shown that substance adhesion (coefficient $\alpha$) affects the evaporation rate but not the final vapor mass. Calculations were performed for six model substances (TNT, RDX, PETN, amphetamine, cocaine, heroin) with a wide range of properties. At room temperature and a surface concentration of 100 ng/cm2, most substances evaporate completely, except for RDX, which remains on the surface at 84%. Evaporation times range from fractions of a second (amphetamine) to several hours (heroin). For low-volatility substances, the maximum mass capable of evaporating under given conditions is determined. The novelty of the work lies in the development of a universal model applicable to a broad class of hazardous substances and in identifying key parameters governing the evaporation process. The results enable the estimation of detection limits for trace substances using vapor-based methods and can be applied in the design of security systems.

The objective of this study is to determine the best formations for the joint movement of a group of small robots in an aquatic environment. Estimation of drag of the flow is a traditional and well-known area of research, but it is not always valid to extend the conclusions made for a single robot to a group of similar devices due to the physical effects that appear during joint movement, such as a wave shadow. For these reasons, it is necessary to study the hydrodynamic characteristics of certain robot formations as a stable structure. The hydrodynamic parameters of systems with two main types of propulsion were studied: locomotive (fishtails) and propellers. Formations similar in structure to schools of fish were mainly considered, and then their applicability for robots of different types was assessed. The relationship between the speed of movement of the group and the drag of each of its participants was also studied. Mathematical modeling of the flow around a group of robots was performed using the finite volume method using two software packages (FlowVision and OpenFoam). Robots with a screw propeller interfere with each other when packed into tight formations, and for the locomotive case, being in the disturbance zone, on the contrary, is preferable. Also, with poorly streamlined bodies, flows separating from the surface can turn into narrow turbulent jets that greatly interfere with the rear robots. It has been established that wake effect reduces energy costs only at low speeds of movement — about 5 cm/s; at high speeds, movement in columns becomes difficult for the rear robots. No large difference in frontal resistance was found between a single robot and a group for a fish-like tail. The studies made it possible to develop and substantiate recommendations for optimizing robot designs for group movement.

Currently, a significant development has been observed in the direction of distributed computing theory, where computational tasks are solved collectively by resource-constrained devices. In practice, this scenario is implemented when processing data in Internet of Things systems, with the aim of reducing system latency and network infrastructure load, as data is processed on edge network computing devices. However, the rapid growth and widespread adoption of IoT systems raise questions about the need to develop methods for reducing the resource intensity of computations. The resource constraints of computing devices pose the following issues regarding the distribution of computational resources: firstly, the necessity to account for the transit cost between different devices solving various tasks; secondly, the necessity to consider the resource cost associated directly with the process of distributing computational resources, which is particularly relevant for groups of autonomous devices such as drones or robots. An analysis of modern publications available in open access demonstrated the absence of proposed models or methods for distributing computational resources that would simultaneously take into account all these factors, making the creation of a new mathematical model for organizing distributed computing in IoT systems and its solution methods topical. This article proposes a novel mathematical model for distributing computational resources along with heuristic optimization methods, providing an integrated approach to implementing distributed computing in IoT systems. A scenario is considered where there exists a leader device within a group that makes decisions concerning the allocation of computational resources, including its own, for distributed task resolution involving information exchanges. It is also assumed that no prior knowledge exists regarding which device will assume the role of leader or the migration paths of computational tasks across devices. Experimental results have shown the effectiveness of using the proposed models and heuristics: achieving up to a 52% reduction in resource costs for solving computational problems while accounting for data transit costs, saving up to 73% of resources through supplementary criteria optimizing task distribution based on minimizing fragment migrations and distances, and decreasing the resource cost of resolving the computational resource distribution problem by up to 28 times with reductions in distribution quality up to 10%.

The effectiveness of communication and data transmission systems (CSiPS), which are an integral part of modern systems in almost any field of science and technology, largely depends on the stability of the frequency of the generated signals. The signals generated in the CSiPD can be considered as processes, the frequency of which changes under the influence of a combination of external influences. Changing the frequency of the signals leads to a decrease in the signal-tonoise ratio (SNR) and, consequently, a deterioration in the characteristics of the signal-to-noise ratio, such as the probability of a bit error and bandwidth. It is most convenient to consider the description of such changes in the frequency of signals as random processes, the apparatus of which is widely used in the construction of mathematical models describing the functioning of systems and devices in various fields of science and technology. Moreover, in many cases, the characteristics of a random process, such as the distribution law, mathematical expectation, and variance, may be unknown or known with errors that do not allow us to obtain estimates of the signal parameters that are acceptable in accuracy. The article proposes an algorithm for solving the problem of determining the characteristics of a random process (signal frequency) based on a set of samples of its frequency, allowing to determine the sample mean, sample variance and the distribution law of frequency deviations in the general population. The basis of this algorithm is the comparison of the values of the observed random process measured over a certain time interval with a set of the same number of random values formed on the basis of model distribution laws. Distribution laws based on mathematical models of these systems and devices or corresponding to similar systems and devices can be considered as model distribution laws. When forming a set of random values for the accepted model distribution law, the sample mean value and variance obtained from the measurement results of the observed random process are used as mathematical expectation and variance. The feature of the algorithm is to compare the measured values of the observed random process ordered in ascending or descending order and the generated sets of values in accordance with the accepted models of distribution laws. The results of mathematical modeling illustrating the application of this algorithm are presented.

This paper proposes an architectural solution for organizing a convolutional neural network (CNN) oriented towards hardware implementation on edge devices under limited resources. To this goal, an approach to compressing spectrograms to a given size (28 × 28) is proposed using discretization, monoconversion, windowed Fourier transform, and two-dimensional interpolation. A balanced convolution procedure is developed based on compact convolutional filters, the size of which provides the balance between computational complexity and accuracy required for edge devices. An algorithm that enables convolution operations and calculation of the error function gradient in the convolutional layer in a single cycle ensuring increased performance in both inference and training modes of the CNN is proposed. The tradeoff between network trainability and its resistance to overfitting is optimized by applying the Dropout regularization method with a dropout coefficient of 0.5 for the fully connected layer.

The effectiveness of the proposed solution was demonstrated using the example of recognizing audio spectrograms of car and airplane engine sounds. The CNN was trained on a balanced dataset consisting of 7160 audio recordings. The trained network demonstrated high recognition accuracy (95%), low loss values (< 0.2), and balanced precision/recall/F-metric, demonstrating the effectiveness of the developed CNN model.

For the first time, numerical modeling has determined the mechanism of the initiated-by-cavitation rise of the liquid level in tubes and capillaries, known as the laser-induced optocapillary effect, as well as its analogues — the acoustocapillary and plasmocapillary effects. It is shown that the key condition for the occurrence of the liquid rise is the asymmetric collapse of a single relatively large cavitation bubble inside a vertically oriented tube or capillary. The proximity of boundaries (the tube wall, the fiber optic tip, and others) disrupts the spherical symmetry of the bubble during its collapse, leading to the appearance of a liquid flow that rolls up into a long-lived toroidal vortex ring. Due to viscous entrainment of the surrounding medium, the vortex generates a directed liquid flow upward and also ensures the suction of a new portion of liquid through the open lower end of the tube. The simulation results show that the characteristic lifetime of the toroidal vortex significantly exceeds the duration of the growth and collapse stages of the cavitation bubble that generated it. It is demonstrated that the rise of the liquid level in the tube does not begin at the moment of bubble expansion, but after its complete disappearance, and continues over a relatively long period due to the inertia of the vortex motion. This result is in complete agreement with experimental data, confirming the validity of the proposed mechanism.

The study investigated the practically significant configuration of the laser-induced optocapillary effect using an optical fiber. This configuration opens broad prospects for technical and medical applications, particularly in laser surgery. The investigated mechanism can be used to create cavitation pumps — effective tools for cleaning technical surfaces and wound surfaces, where the process of removing damaged tissue and foreign bodies due to thermal exposure will be accompanied by the removal of debris through the tube, significantly increasing the efficiency and safety of the procedure.

The obtained results represent the first consistent explanation for a class of cavitation-induced capillary phenomena and create a foundation for their controlled application in biomedical and microfluidic technologies.