Все выпуски

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

Работа будет полезной специалистам в области математического моделирования, процессов тепло- и массообмена, физики низкотемпературной плазмы, аспирантам и студентам старших курсов, специализирующимся в указанных направлениях.

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.

Безопасное проведение аварийно-восстановительных работ на аварийных морских газоконденсатных скважинах возможно при учете опасных факторов, препятствующих проведению противофонтанных мероприятий. Одним из таких факторов является загазованность района работ вследствие выхода из водной толщи большого количества легкого, по сравнению с воздухом, природного газа, а также паров более тяжелых компонентов газового конденсата (ГК). Для оценки распределения взрывоопасных концентраций паров нефтепродукта в приводном слое атмосферы необходимо определить характеристики источника загазованности. На основании анализа теоретических работ, посвященных формированию поля скорости в верхнем слое моря вследствие выхода на поверхность большого количества газа, предложена аналитическая модель для расчета размеров области, в которой происходит испарение значительного количества поступающего на поверхность ГК при авариях на мелководных скважинах. Рассматривается стационарный режим истечения пластового продукта при открытом фонтанировании газонефтяных скважин морского базирования при подводном расположении их устья. Построена малопараметрическая модель испарения нефтепродуктов из пленок различной толщины. Показано, что размер зоны интенсивного испарения ГК при подводном выбросе на мелководных скважинах определяется объемным потоком жидкой фракции ГК, его фракционным составом и выбранным порогом для оценки потока паров нефтепродукта в атмосферу. В контексте данной работы мелководными называются скважины при дебите газа от 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.

Групповое движение малоразмерных подводных аппаратов — важная прикладная задача. В работе приводятся результаты исследования влияния формации группы на характер ее движения. Оценка лобового сопротивления подводных аппаратов и обтекания потоков вокруг них — традиционная и хорошо известная область исследований. Однако выводы, сделанные для единичного робота, не всегда применимы к группе однотипных устройств из-за появляющихся при совместном движении физических эффектов, например волновой тени. Исходя из этого были исследованы гидродинамические характеристики определенных формаций роботов, движущихся как единое целое. В ходе работы изучались гидродинамические параметры систем с двумя основными типами движителей: локомоторными (аналогами рыбьих хвостов) и гребными винтами. Из соображений природоподобия рассматривались формации, аналогичные по структуре рыбьим косякам, затем оценивалась их применимость для роботов разных видов. Была определена связь между скоростью движения группировки и лобовым сопротивлением каждого из ее участников. Математическое моделирование обтекания группировки роботов проводилось при помощи метода конечных объемов двумя программными комплексами (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.

Преобразование Фурье является основным инструментом для оценки спектральных характеристик турбулентного течения. Применение преобразования Фурье (как правило, дискретного) к первой степени пульсационной продольной или поперечной компоненты вектора скорости позволяет оценить спектральную плотность энергии (ESD) или мощности (PSD). Для оценки ESD и PSD турбулентного сигнала, полученного по результатам численного моделирования обтекания тела, создается массив значений сигнала с дискретизацией во временной или пространственной области. Референтное распределение спектральной характеристики ESD (закон масштабирования) в области волновых чисел инерциального поддиапазона следует из двух гипотез А.Н. Колмогорова и определяется законом $−\frac{5}{3}$. Закон $−\frac{5}{3}$ используется также в большинстве работ для оценки распределения ESD и PSD в частотной области. При этом распределение спектра мощности PSD получается из распределения спектра энергии ESD нормировкой к времени сканирования сигнала. Альтернативная спектральная характеристика энергии (ESS) пульсаций скорости может быть определена преобразованием Фурье для квадрата пульсаций скорости. Размерность ESS в пространственной области совпадает с размерностью закона $−\frac{5}{3}$ А.Н. Колмогорова в области волновых чисел. При дискретизации сигнала во временной области для ESS ранее получен закон масштабирования $−2$ в частотной области. В настоящей работе вводится альтернативная оценка спектра мощности сигнала (PSS), полученная с использованием преобразования Фурье для третьей степени пульсаций скорости. Из гипотез А.Н. Колмогорова следует, что в частотной области закон масштабирования спектра PSS определяется степенью $−\frac{5}{2}$. В качестве приложения рассматривается нестационарное обтекание отсека цилиндрической 3D-поверхности при числе Рейнольдса 3900. Численное моделирование выполнено с использованием пакета ANSYS Fluent на базе решения уравнений Навье – Стокса в несжимаемой постановке. Пространственно-временные характеристики вектора скорости турбулентного потока анализируются с применением правильной ортогональной декомпозиции (POD). Для оценки ESS и PSS используется преобразование Фурье с дискретизацией сигнала по времени.

The Fourier transformation is the basic tool for evaluating the spectral characteristics of a turbulent flow. The Fourier transform (usually discrete) of the first power of the longitudinal or transverse component of the velocity vector pulsations allows estimation of the energy spectral density (ESD) or power spectral density (PSD). To estimate the ESD and PSD of a turbulent signal obtained from numerical simulation, an array of signal values with discretization in the time or spatial domain is generated. The reference distribution of the ESD spectrum (scaling law) within the wave number domain of the inertial subrange is derived from two hypotheses proposed by A.N. Kolmogorov and is characterized by the $−\frac{5}{3}$ law. The $-\frac{5}{3}$ law is also used in most references to estimate the ESD distribution in the frequency domain. The distribution of the power spectrum PSD is derived from the distribution of the energy spectrum ESD by normalizing to the signal scanning time. An alternative energy spectral characteristic (ESS) of velocity fluctuations can be determined by the Fourier transform of the square of the velocity fluctuations. In the wave numbers domain, the dimension of ESS in the spatial domain coincides with the dimension of A.N. Kolmogorov's $−\frac{5}{3}$ law. When considering a signal sampled in the time domain, a scaling law of $−2$ for ESS was previously obtained in the frequency domain. An alternative estimate of the Power Signal Spectrum (PSS) is discussed in this paper based on the Fourier transform of the third-order velocity pulsations. Based on the hypotheses proposed by A.N.Kolmogorov, it can be inferred that in the frequency domain, the scaling law of the PSS spectrum is characterized by the power of $−\frac{5}{2}$. Unsteady incompressible flow around a 3D cylindrical surface section at the Reynolds number of 3900 is considered as an application. The numerical simulation is performed using ANSYS Fluent commercial code and based on the Navier – Stokes equations. The spatio-temporal characteristics of the turbulent flow velocity vector are analyzed using the Proper Orthogonal Decomposition (POD). The Fourier transform is used to estimate the ESS and PSS of a time-sampled signal.

В результате численного моделирования впервые определен механизм инициированного кавитацией подъема уровня жидкости в трубках и капиллярах, известного как лазерно-индуцированный светокапиллярный эффект, а также его аналоги — звукокапиллярный и плазмокапиллярный эффекты. Показано, что ключевым условием возникновения подъема жидкости является асимметричное схлопывание относительно крупного одиночного кавитационного пузырька внутри вертикально ориентированной трубки или капилляра. Близость границ (стенки трубки, торца оптоволокона и другие) нарушает сферическую симметрию пузырька при его коллапсе, что приводит к появлению потока жидкости, который замыкается в долгоживущее тороидальное вихревое кольцо (тороидальный вихрь). Вихрь за счет вязкого увлечения окружающей среды генерирует направленное течение жидкости вверх, а также обеспечивает всасывание новой порции жидкости через открытый нижний конец трубки. Результаты моделирования показывают, что характерное время жизни тороидального вихря значительно превышает время стадий роста и схлопывания кавитационного пузырька, который его породил. Показано, что подъем уровня жидкости в трубке начинается не в момент расширения пузырька, а после его полного исчезновения и за счет инерции вихревого движения продолжается на протяжении относительно длительного промежутка времени. Этот результат полностью согласуется с экспериментальными данными, что подтверждает достоверность предложенного механизма. В работе исследована практически значимая конфигурация лазерно-индуцированного светокапиллярного эффекта с использованием оптоволокна. Такая конфигурация открывает широкие перспективы для технических и медицинских приложений, в частности в лазерной хирургии. Исследованный механизм может быть использован для создания кавитационных насосов — эффективных инструментов очистки технических поверхностей и поверхностей ран, где процесс удаления поврежденных тканей и инородных тел в результате температурного воздействия будет сопровождаться удалением продуктов очистки через трубку, что существенно повышает эффективность и безопасность процедуры. Полученные результаты представляют собой первое непротиворечивое объяснение класса кавитационно-индуцированных капиллярных явлений и создают основу для их контролируемого применения в биомедицинских и микрофлюидных технологиях.

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.

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

Efficient method for numerical solving of the steady transport equation in x-y-z-geometry has been suggested. The equation is being solved on hexagonal mesh, reflecting real structure of the reactor active zone cross-section. Method of characteristics is used, that inherits all the outcomes from the two-dimensional r-z-geometry calculation. Two variants of the method of characteristics have been applied for solving the transport equation in a cell: method of short characteristics and its conservative modification. It has been confirmed that in three-dimensional geometry conservative method has advantage over pure characteristic and it produces highly accurate solution, especially for quasi-diffusion tensor components.

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

Based on contents of phytoplankton pigments, fluorescence samples and some physico-chemical characteristics of the Rybinsk reservoir waters, searching for connections between biological and physicalchemical characteristics is working out. The standard methods of statistical analysis (correlation, regression), methods of description of connection between qualitative classes of characteristics, based on deviation of the studied characteristics distribution from independent distribution, are studied. A method of searching for boundaries of quality classes by criterion of maximum connection coefficient is offered.

В рамках метода молекулярной динамики, с использованием многочастичного потенциаламежатомного взаимодействия, на модели свободной частицы Ni₆₀Ag₄₀ диаметром 40 Å, исследована структурная эволюция в процессе закалки. Проведен сравнительный анализ структурных перестроек для моделей частицы и массивного образца. Отмечено снижение температуры начала и конца превращений для частицы. Показано, что в образовании перколяционного кластера из взаимопроникающих и контактирующих между собой икосаэдров, для модели свободной частицы, задействовано на 10 % атомов больше, чем для модели массивного тела.

Molecular dynamics simulation using the embedded-atom method is applied to study the structural evolution of the particle diameter of 40 Å during the quenching process. Was carried comparative analysis of the structural reconstruction for the particle and the bulk models. Was a reduction in temperature of the beginning and end of the transformation of the particle. In formation of a percolation cluster from interpenetrating and contacting icosahedrons, for model of the particle, it is involved for 10 percent of atoms more, than for model of a bulk.

В работе на основе предложенной ранее русловой модели решена одномерная задача устойчивости песчаного дна напорного канала. Особенностью исследуемой задачи является используемое оригинальное уравнение русловых деформаций, учитывающее влияние физико-механических и гранулометрических характеристик донного материала и неровности донной поверхности при русловом анализе. Еще одной особенностью рассматриваемой задачи является учет влияния не только придонного касательного, но и нормального напряжения при изучении русловой неустойчивости. Из решения задачи устойчивости песчаного дна для напорного канала получена аналитическая зависимость, определяющая длину волны для быстрорастущих донных возмущений. Выполнен анализ полученной аналитической зависимости, показано, что она обобщает ряд известных эмпирических формул: Коулмана, Шуляка и Бэгнольда. Структура полученной аналитической зависимости указывает на существование двух гидродинамических режимов, характеризуемых числом Фруда, при которых рост донных возмущений может сильно или слабо зависеть от числа Фруда. Учитывая природную стохастичность процесса движения донных волн и наличие области определения решения со слабой зависимостью от чисел Фруда, можно сделать вывод о том, что экспериментальное наблюдение за процессом развития движения донных волн в данной области должно приводить к получению данных, имеющих существенную дисперсию, что и происходит в действительности.

In this paper on the basis of the riverbed model proposed earlier the one-dimensional stability problem of closed flow channel with sandy bed is solved. The feature of the investigated problem is used original equation of riverbed deformations, which takes into account the influence of mechanical and granulometric bed material characteristics and the bed slope when riverbed analyzing. Another feature of the discussed problem is the consideration together with shear stress influence normal stress influence when investigating the riverbed instability. The analytical dependence determined the wave length of fast-growing bed perturbations is obtained from the solution of the sandy bed stability problem for closed flow channel. The analysis of the obtained analytical dependence is performed. It is shown that the obtained dependence generalizes the row of well-known empirical formulas: Coleman, Shulyak and Bagnold. The structure of the obtained analytical dependence denotes the existence of two hydrodynamic regimes characterized by the Froude number, at which the bed perturbations growth can strongly or weakly depend on the Froude number. Considering a natural stochasticity of the waves movement process and the presence of a definition domain of the solution with a weak dependence on the Froude numbers it can be concluded that the experimental observation of the of the bed waves movement development should lead to the data acquisition with a significant dispersion and it occurs in reality.