All issues

Безопасное проведение аварийно-восстановительных работ на аварийных морских газоконденсатных скважинах возможно при учете опасных факторов, препятствующих проведению противофонтанных мероприятий. Одним из таких факторов является загазованность района работ вследствие выхода из водной толщи большого количества легкого, по сравнению с воздухом, природного газа, а также паров более тяжелых компонентов газового конденсата (ГК). Для оценки распределения взрывоопасных концентраций паров нефтепродукта в приводном слое атмосферы необходимо определить характеристики источника загазованности. На основании анализа теоретических работ, посвященных формированию поля скорости в верхнем слое моря вследствие выхода на поверхность большого количества газа, предложена аналитическая модель для расчета размеров области, в которой происходит испарение значительного количества поступающего на поверхность ГК при авариях на мелководных скважинах. Рассматривается стационарный режим истечения пластового продукта при открытом фонтанировании газонефтяных скважин морского базирования при подводном расположении их устья. Построена малопараметрическая модель испарения нефтепродуктов из пленок различной толщины. Показано, что размер зоны интенсивного испарения ГК при подводном выбросе на мелководных скважинах определяется объемным потоком жидкой фракции ГК, его фракционным составом и выбранным порогом для оценки потока паров нефтепродукта в атмосферу. В контексте данной работы мелководными называются скважины при дебите газа от 1 до 20 млн м³ на глубинах порядка 50–200 метров. В этом случае струя пластового флюида из устья скважины на морском дне трансформируется в пузырьковый шлейф, типичная для летне-осеннего периода стратификация водной толщи не ограничивает выход шлейфа на поверхность моря, а скорость подъема пузырьков позволяет не принимать во внимание процесс растворения газа. Проведенный анализ был ограничен условиями близкими к штилевым. Такие условия благоприятны для проведения морских операций, однако неблагоприятны с точки зрения рассеяния высоких концентраций паров нефтепродуктов в приводном слое атмосферы над морем. В результате проведенной работы предложено аналитическое соотношение для приближенной оценки зоны интенсивного испарения ГК.

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.

В статье исследуется возможность обнаружения следов опасных веществ (взрывчатых и наркотических) на основе детекции их паров в воздухе. Актуальность работы обусловлена задачами противодействия террористическим угрозам и наркотрафику, где критически важно определять даже следовые количества веществ. Основное внимание уделено математическому моделированию испарения тонкого слоя вещества с поверхности, основанному на молекулярно-кинетической теории. Предложена универсальная модель, учитывающая физико-химические свойства веществ, температуру окружающей среды, адгезию к поверхности и начальную массу слоя. На основе уравнений Герца – Кнудсена – Ленгмюра и Клаузиуса – Клапейрона получены аналитические выражения для времени полного испарения, предельной массы паров и динамики процесса. Выявлен безразмерный параметр $\gamma$, определяющий предельные условия испарения. Показано, что адгезия вещества (коэффициент $\alpha$) влияет на скорость испарения, но не на конечную массу паров. Проведены расчеты для шести модельных веществ (TNT, RDX, PETN, амфетамин, кокаин, героин) с широким диапазоном свойств. Установлено, что при комнатной температуре и поверхностной концентрации 100 нг/см2 большинство веществ испаряются полностью, за исключением RDX, который остается на поверхности на 84%. Время испарения варьируется от долей секунды (амфетамин) до нескольких часов (героин). Для веществ с низкой летучестью определена максимальная масса, способная испариться при заданных условиях. Новизна работы заключается в разработке универсальной модели, применимой для широкого класса опасных веществ, и в выявлении ключевых параметров, определяющих процесс испарения. Полученные результаты позволяют оценить пределы обнаружения следов веществ методами, основанными на регистрации паров, и могут быть использованы при проектировании систем безопасности.

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.

Групповое движение малоразмерных подводных аппаратов — важная прикладная задача. В работе приводятся результаты исследования влияния формации группы на характер ее движения. Оценка лобового сопротивления подводных аппаратов и обтекания потоков вокруг них — традиционная и хорошо известная область исследований. Однако выводы, сделанные для единичного робота, не всегда применимы к группе однотипных устройств из-за появляющихся при совместном движении физических эффектов, например волновой тени. Исходя из этого были исследованы гидродинамические характеристики определенных формаций роботов, движущихся как единое целое. В ходе работы изучались гидродинамические параметры систем с двумя основными типами движителей: локомоторными (аналогами рыбьих хвостов) и гребными винтами. Из соображений природоподобия рассматривались формации, аналогичные по структуре рыбьим косякам, затем оценивалась их применимость для роботов разных видов. Была определена связь между скоростью движения группировки и лобовым сопротивлением каждого из ее участников. Математическое моделирование обтекания группировки роботов проводилось при помощи метода конечных объемов двумя программными комплексами (FlowVision и OpenFoam). Показано, что роботы с винтовым движителем при размещении в тесных формациях мешают друг другу, а для локомоторного случая нахождение в зоне возмущения, наоборот, предпочтительно. Также при плохо обтекаемых корпусах отрывающиеся от поверхности потоки могут превращаться в узкие струи, сильно мешающие задним роботам. Установлено, что эффект водяной тени снижает затраты энергии только при малых скоростях движения — около 5 см/с; при больших скоростях движение в колоннах затрудняется для задних роботов. Кроме того, для рыбоподобного движителя не было выявлено большой разницы в лобовом сопротивлении между одиночным роботом и группой. Таким образом, программное моделирование позволило выработать и обосновать рекомендации по оптимизации построений роботов при групповом движении. Полученные результаты могут оказаться полезными для разработки подводных аппаратов, способных работать в группах, и средств управления ими.

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.

Эффективность систем связи и передачи данных (ССиПД), являющихся неотъемлемой составляющей современных систем практически в любой области науки и техники, во многом зависит от стабильности частоты формируемых сигналов. Формируемые в ССиПД сигналы могут рассматриваться как процессы, частота которых изменяется под действием совокупности внешних воздействий. Изменение частоты сигналов приводит к уменьшению отношения «сигнал/шум» (ОСШ) и, соответственно, ухудшению характеристик ССиПД, таких как вероятность битовой ошибки, пропускная способность. Описание таких изменений частоты сигналов наиболее удобно рассматривать как случайные процессы, аппарат которых находит широкое применение при построении математических моделей, описывающих функционирование систем и устройств в различных областях науки и техники. При этом во многих случаях характеристики случайного процесса, такие как закон распределения, математическое ожидание и дисперсия, могут являться неизвестными или известными с погрешностями, не позволяющими получить приемлемые по точности оценки параметров сигналов. В статье предлагается алгоритм решения задачи по определению характеристик случайного процесса (частоты сигнала) на основе набора отсчетов его частоты, позволяющих определить выборочное среднее, выборочную дисперсию и закон распределения отклонений частоты в генеральной совокупности. Основой данного алгоритма является сравнение измеренных на некотором временном интервале значений наблюдаемого случайного процесса с набором того же количества случайных значений, сформированных на основе модельных законов распределения. В качестве модельных законов распределения могут рассматриваться законы распределения, принятые на основе математических моделей этих систем и устройств или соответствующие аналогичным системам и устройствам. В качестве математического ожидания и дисперсии при формировании набора случайных значений для принятого модельного закона распределения принимаются выборочные среднее значение и дисперсия, полученные по результатам измерений наблюдаемого случайного процесса. Особенность алгоритма заключается в проведении сравнения упорядоченных по возрастанию или убыванию измеренных значений наблюдаемого случайного процесса и сформированных наборов значений в соответствии с принятыми моделями законов распределения. Приведены результаты математического моделирования, иллюстрирующие применение данного алгоритма.

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.

Для квантового описания поведения двуядерных систем на начальной стадии околобарьерного слияния тяжелых ядер использованы численные методы нахождения коллективных и одночастичных состояний. Коллективные возбужденные состояния в таких системах представляют собой согласованные колебания поверхностей сферических ядер. Одночастичные состояния внешних нейтронов аналогичны состояниям валентных электронов двухатомных молекул.

Numerical methods of obtaining collective and one-particle states were used for the quantum description of two-nuclear systems behavior at the initial stage of near-barrier heavy nuclei fusion. The collective exited states in such systems represent concordant oscillations of surfaces of spherical nuclei. The one-particle states of the external neutrons are similar to the states of valence electrons of diatomic molecules.

Построена математическая модель эволюции микробных популяций при длительном непрерывном культивировании на протоке. Модель представляет собой обобщение целого ряда известных математических моделей эволюции, в которых учитываются такие факторы генетической изменчивости как хромосомные мутации, мутации плазмидных генов, перенос плазмид между клетками микроорганизмов, потери плазмид при делении клеток и др. Для общей модели эволюции построена функция Ляпунова и на основании теоремы Ла-Салля доказано существование в пространстве состояний математической модели ограниченного, положительно инвариантного и глобально притягивающего множества. Дано аналитическое описание этого множества. Обсуждаются перспективы применения численных методов для оценки числа, местоположения и последующего исследования предельных множеств в математических моделях эволюции на протоке.

A mathematical model for the evolution of microbial populations during prolonged cultivation in a chemostat has been constructed. This model generalizes the sequence of the well-known mathematical models of the evolution, in which such factors of the genetic variability were taken into account as chromosomal mutations, mutations in plasmid genes, the horizontal gene transfer, the plasmid loss due to cellular division and others. Liapunov’s function for the generic model of evolution is constructed. The existence proof of bounded, positive invariant and globally attracting set in the state space of the generic mathematical model for the evolution is presented because of the application of La-Salle’s theorem. The analytic description of this set is given. Numerical methods for estimate of the number of limit sets, its location and following investigation in the mathematical models for evolution are discussed.

В работе предложен новый метод моделирования плотных материалов на основе алгоритма упаковки сферополиэдров, описана математическая модель сферополиэдра и обсуждены результаты вычислительных экспериментов на различных упаковках сферополиэдров. Результаты экспериментов показали сходимость метода. Проведенные эксперименты включают исследования упаковок сферополиэдров различной формы, полидисперсных и ориентированных структур. Метод может быть применен для виртуального проектирования плотных материалов, имеющих в составе несферические частицы.

The paper presents a new dense material modeling method based on spherepolyhedra packing algorithm, describes mathematical model of spherepolyhedra and discuss the results of computation experiments on different spherepolyhedra packs. The results of experiments show convergence of proposed method. Experiments include investigations of spherepolyhedra packs with different shapes, polydisperse and oriented structures. Presented method would be applied to virtual design of dense materials composed of non-spherical particles.

Данная работа рассматривает задачу оптимального управления гальваническим процессом в многоанодной ванне. Построена нестационарная математическая модель гальванического процесса, которая учитывает изменения концентрации компонентов электролита. Продемонстрировано обоснование выбора вида управляющих экстремалей на примере гальванического процесса хромирования в стандартном электролите.

This work considers the problem of optimal control galvanic process in multianode bath. The nonstationary mathematical model of galvanic process, which considers change concentrations of electrolyte components, is developed. Demonstrated rationale for the choice of the form to extremal control on example chrome galvanic process in the standard electrolyte.

В данной работе представлены результаты численного моделирования распределения магнитного поля некоторых модификаций спектрометрического магнита СП-40, используемого в экспериментальной установке НИС ОИЯИ. Основной особенностью такого магнита является прямоугольная апертура, следовательно, и область, в которой решается краевая задача, имеет всюду гладкую границу, за исключением угловой области ферромагнетика. В таких случаях решение задачи или производные решения могут иметь особенность. Изучено поведение магнитного поля в окрестности угловой области ферромагнетика и численным путем выбрана конфигурация магнита, для которой ширина области однородности (т. е. с $\Delta B/B < 1\,\%$) магнитного поля увеличилась с 0.5 м до 1.0 м, т. е. в два раза.

In this work, we propose a method of the numerical solution of the magnetostatic problem for domains with boundaries containing corners. With the help of this numerical method, the magnetic systems of rectangular configuration were simulated with high accuracy. In particular, the calculations of some modifications of the magnetic system SP-40 used in the NIS JINR experimental installation, are presented. The basic feature of such a magnet is a rectangular aperture, hence, the area in which the boundary-value problem is solved, has a smooth border everywhere, except for a finite number of angular points in the vicinity of which the border is formed by crossing two smooth curves. In such cases the solution to the problem or derivatives of the solution can have a special feature. A behavior of the magnetic field in the vicinity of an angular point is investigated, and the configuration of the magnet was chosen numerically. The width of the area of homogeneity of the magnetic field increased from 0.5 m up to 1.0 m, i. e. twice.

Проведен численный анализ влияния формы и размеров локального источника постоянной температуры на нестационарные режимы термогравитационной конвекции в квадратной полости с изотермическими вертикальными стенками. Рассматривался источник энергии прямоугольной, треугольной и трапециевидной формы. Краевая задача, сформулированная в безразмерных преобразованных переменных «функция тока – завихренность скорости – температура» в приближении Буссинеска, была реализована численно методом конечных разностей. Получены распределения изолиний функции тока и температуры, а также временные зависимости для среднего числа Нуссельта на поверхности источника энергии в широком диапазоне изменения определяющих параметров.

Numerical analysis of the effects of the local heat source shape on transient natural convection in a square enclosure has been carried out. The local heat source has rectangular, triangular and trapezoidal shape. The boundary value problem formulated in the dimensionless variables such as stream function, vorticity and temperature by using the Boussinesq approximation has been solved by means of finite difference method. Distributions of streamlines and isotherms and time dependences for the average Nusselt number along the heat source surface in a wide range of governing parameters have been obtained.