All issues
- 2024 Vol. 16
- 2023 Vol. 15
- 2022 Vol. 14
- 2021 Vol. 13
- 2020 Vol. 12
- 2019 Vol. 11
- 2018 Vol. 10
- 2017 Vol. 9
- 2016 Vol. 8
- 2015 Vol. 7
- 2014 Vol. 6
- 2013 Vol. 5
- 2012 Vol. 4
- 2011 Vol. 3
- 2010 Vol. 2
- 2009 Vol. 1
-
Движение влекомых наносов над периодическим дном
Компьютерные исследования и моделирование, 2018, т. 10, № 1, с. 47-60Движение влекомых наносов по дну напорного канала может приводить к потере устойчивости донной поверхности, когда на дне канала возникают донные волны. Исследование процесса развития донных волн связано с возможностью определения характера движения влекомых наносов по дну периодической формы. Несмотря на большое внимание многих исследователей к данной проблеме, вопрос о развитии процесса донных волн остается открытым и в настоящее время. В значительной мере это связано с тем, что при анализе данного процесса многие исследователи используют в своих работах феноменологические формулы движения влекомых наносов. Полученные в таких моделях результаты позволяют лишь качественно оценить процесс развития донных волн. По этой причине представляет интерес проведение анализа развития донных волн с использованием аналитической модели движения влекомых наносов.
В работе предложена двумерная профильная математическая русловая модель, позволяющая исследовать движение влекомых наносов над периодическим дном. Особенностью математической модели является возможность расчета расхода влекомых наносов по аналитической модели с реологией Кулона–Прандтля, учитывающей влияние уклонов поверхности дна, придонных нормальных и касательных напряжений на процесс движения донного материала. Показано, что при движении донного материла по дну периодической формы диффузионные и напорные расходы влекомых наносов являются разнонаправленными и доминирующими по отношению к транзитному расходу. Изучались влияния параметра перекошенности донной волны на вклад транзитного, диффузионного и напорного расходов в полный расход влекомых наносов. Выполнено сравнение полученных результатов с численными решениями других авторов для донной поверхности косинусоидальной формы.
Ключевые слова: математическое моделирование, напорный канал, донные волны, влекомые наносы, расход влекомых наносов.
Movement of sediment over periodic bed
Computer Research and Modeling, 2018, v. 10, no. 1, pp. 47-60Views (last year): 9.The movement of bed load along the closed conduit can lead to a loss of stability of the bed surface, when bed waves arise at the bed of the channel. Investigation of the development of bed waves is associated with the possibility of determining of the bed load nature along the bed of the periodic form. Despite the great attention of many researchers to this problem, the question of the development of bed waves remains open at the present time. This is due to the fact that in the analysis of this process many researchers use phenomenological formulas for sediment transport in their work. The results obtained in such models allow only assess qualitatly the development of bed waves. For this reason, it is of interest to carry out an analysis of the development of bed waves using the analytical model for sediment transport.
The paper proposed two-dimensional profile mathematical riverbed model, which allows to investigate the movement of sediment over a periodic bed. A feature of the mathematical model is the possibility of calculating the bed load transport according to an analytical model with the Coulomb–Prandtl rheology, which takes into account the influence of bottom surface slopes, bed normal and tangential stresses on the movement of bed material. It is shown that when the bed material moves along the bed of periodic form, the diffusion and pressure transport of bed load are multidirectional and dominant with respect to the transit flow. Influence of the effects of changes in wave shape on the contribution of transit, diffusion and pressure transport to the total sediment transport has been studied. Comparison of the received results with numerical solutions of the other authors has shown their good qualitative initiation.
-
Развитие берегового откоса в русле трапециевидного канала
Компьютерные исследования и моделирование, 2022, т. 14, № 3, с. 581-592Сформулирована математическая модель эрозии берегового склона песчаного канала, происходящей под действием проходящей паводковой волны. Модель включает в себя уравнение движения квазиустановившегося гидродинамического потока в створе канала. Движение донной и береговой поверхности русла определяется из решения уравнения Экснера, которое замыкается оригинальной аналитической моделью движения влекомых наносов. Модель учитывает транзитные, гравитационные и напорные механизмы движения донного материала и не содержит в себе феноменологических параметров. Движение свободной поверхности гидродинамического потока определяется из решения дифференциальных уравнений баланса. Модель учитывает изменения средней по створу турбулентной вязкости при изменении створа канала.
На основе метода конечных элементов получен дискретный аналог сформулированной задачи и предложен алгоритм ее решения. Особенностью алгоритма является контроль влияния движения свободной поверхности потока и расхода потока на процесс определения турбулентной вязкости потока в процессе эрозии берегового склона. Проведены численные расчеты, демонстрирующие качественное и количественное влияние данных особенностей на процесс определения турбулентной вязкости потока и эрозию берегового склона русла.
Сравнение данных по береговым деформациям, полученных в результате численных расчетов, с известными лотковыми экспериментальными данными показали их согласование.
Bank slope evolution in trapezoidal channel riverbed
Computer Research and Modeling, 2022, v. 14, no. 3, pp. 581-592A mathematical model is formulated for the coastal slope erosion of sandy channel, which occurs under the action of a passing flood wave. The moving boundaries of the computational domain — the bottom surface and the free surface of the hydrodynamic flow — are determined from the solution of auxiliary differential equations. A change in the hydrodynamic flow section area for a given law of change in the flow rate requires a change in time of the turbulent viscosity averaged over the section. The bottom surface movement is determined from the Exner equation solution together with the equation of the bottom material avalanche movement. The Exner equation is closed by the original analytical model of traction loads movement. The model takes into account transit, gravitational and pressure mechanisms of bottom material movement and does not contain phenomenological parameters.
Based on the finite element method, a discrete analogue of the formulated problem is obtained and an algorithm for its solution is proposed. An algorithm feature is control of the free surface movement influence of the flow and the flow rate on the process of determining the flow turbulent viscosity. Numerical calculations have been carried out, demonstrating qualitative and quantitative influence of these features on the determining process of the flow turbulent viscosity and the channel bank slope erosion.
Data comparison on bank deformations obtained as a result of numerical calculations with known flume experimental data showed their agreement.
-
Устойчивость дна в напорных каналах
Компьютерные исследования и моделирование, 2015, т. 7, № 5, с. 1061-1068В работе на основе предложенной ранее русловой модели решена одномерная задача устойчивости песчаного дна напорного канала. Особенностью исследуемой задачи является используемое оригинальное уравнение русловых деформаций, учитывающее влияние физико-механических и гранулометрических характеристик донного материала и неровности донной поверхности при русловом анализе. Еще одной особенностью рассматриваемой задачи является учет влияния не только придонного касательного, но и нормального напряжения при изучении русловой неустойчивости. Из решения задачи устойчивости песчаного дна для напорного канала получена аналитическая зависимость, определяющая длину волны для быстрорастущих донных возмущений. Выполнен анализ полученной аналитической зависимости, показано, что она обобщает ряд известных эмпирических формул: Коулмана, Шуляка и Бэгнольда. Структура полученной аналитической зависимости указывает на существование двух гидродинамических режимов, характеризуемых числом Фруда, при которых рост донных возмущений может сильно или слабо зависеть от числа Фруда. Учитывая природную стохастичность процесса движения донных волн и наличие области определения решения со слабой зависимостью от чисел Фруда, можно сделать вывод о том, что экспериментальное наблюдение за процессом развития движения донных волн в данной области должно приводить к получению данных, имеющих существенную дисперсию, что и происходит в действительности.
Bottom stability in closed conduits
Computer Research and Modeling, 2015, v. 7, no. 5, pp. 1061-1068Views (last year): 1. Citations: 2 (RSCI).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.
-
Модель установившегося течения реки в поперечном сечении изогнутого русла
Компьютерные исследования и моделирование, 2024, т. 16, № 5, с. 1163-1178Моделирование русловых процессов при исследовании береговых деформаций русла требует вычисления параметров гидродинамического потока, учитывающих существование вторичных поперечных течений, формирующихся на закруглении русла. Трехмерное моделирование таких процессов на текущий момент возможно только для небольших модельных каналов, для реальных речных потоков необходимы модели пониженной размерности. При этом редукция задачи от трехмерной модели движения речного потока к двумерной модели потока в плоскости створа канала предполагает, что рассматриваемый гидродинамический поток является квазистационарным, и для него выполнены гипотезы об асимптотическом поведении потока по потоковой координате створа. С учетом данных ограничений в работе сформулирована математическая модель задачи о движении стационарного турбулентного спокойного речного потока в створе канала. Задача сформулирована в смешанной постановке скорости — «вихрь – функция тока». В качестве дополнительных условий для редукции задачи требуется задание граничных условий на свободной поверхности потока для поля скорости, определяемого в нормальном и касательном направлении к оси створа. Предполагается, что значения данных скоростей должны быть определены из решения вспомогательных задач или получены из данных натурных или экспериментальных измерений.
Для решения сформулированной задачи используется метод конечных элементов в формулировке Петрова – Галёркина. Получен дискретный аналог задачи и предложен алгоритм ее решения. Выполненные численные исследования показали в целом хорошую согласованность полученных решений при их сравнении с известными экспериментальными данными.
Полученные погрешности авторы связывают с необходимостью более точного определения циркуляционного поля скоростей в створе потока путем подбора и калибровки более подходящей модели вычисления турбулентной вязкости и граничных условий на свободной границе створа.
Model of steady river flow in the cross section of a curved channel
Computer Research and Modeling, 2024, v. 16, no. 5, pp. 1163-1178Modeling of channel processes in the study of coastal channel deformations requires the calculation of hydrodynamic flow parameters that take into account the existence of secondary transverse currents formed at channel curvature. Three-dimensional modeling of such processes is currently possible only for small model channels; for real river flows, reduced-dimensional models are needed. At the same time, the reduction of the problem from a three-dimensional model of the river flow movement to a two-dimensional flow model in the cross-section assumes that the hydrodynamic flow under consideration is quasi-stationary and the hypotheses about the asymptotic behavior of the flow along the flow coordinate of the cross-section are fulfilled for it. Taking into account these restrictions, a mathematical model of the problem of the a stationary turbulent calm river flow movement in a channel cross-section is formulated. The problem is formulated in a mixed formulation of velocity — “vortex – stream function”. As additional conditions for problem reducing, it is necessary to specify boundary conditions on the flow free surface for the velocity field, determined in the normal and tangential direction to the cross-section axis. It is assumed that the values of these velocities should be determined from the solution of auxiliary problems or obtained from field or experimental measurement data.
To solve the formulated problem, the finite element method in the Petrov – Galerkin formulation is used. Discrete analogue of the problem is obtained and an algorithm for solving it is proposed. Numerical studies have shown that, in general, the results obtained are in good agreement with known experimental data. The authors associate the obtained errors with the need to more accurately determine the circulation velocities field at crosssection of the flow by selecting and calibrating a more appropriate model for calculating turbulent viscosity and boundary conditions at the free boundary of the cross-section.
-
Моделирование эволюции песчано-гравийного дна канала в одномерном приближении
Компьютерные исследования и моделирование, 2015, т. 7, № 2, с. 315-328В работе предложена математическая модель для одномерного неравновесного руслового процесса. Модель учитывает движение наносов во взвешенном и влекомом состоянии. Транспорт влекомых наносов определен с помощью оригинальной формулы, аналитически полученной из уравнения движения тонкого придонного водогрунтового слоя. Данная формула не содержит новых феноменологических параметров и учитывает влияние уклона дна, физико-механических и гранулометрических параметров донного материала на процесс транспорта влекомых наносов. Для верификации предложенной модели был решен ряд классических тестовых задач. Выполнено сравнение результатов численных расчетов с известными экспериментальными данными и результатами других авторов. Показано, что, несмотря на относительную простоту предложенной математической модели, полученные численные решения хорошо согласуются с экспериментальными данными.
Modeling of sand-gravel bed evolution in one-dimension
Computer Research and Modeling, 2015, v. 7, no. 2, pp. 315-328In the paper the model for a one-dimensional non-equilibrium riverbed process is proposed. The model takes into account the suspended and bed-load sediment transport. The bed-load transport is determined by using the original formula. This formula was derived from the thin bottom layer motion equation. The formula doesn’t contain new phenomenological parameters and takes into account the influence of bed slope, granulometric and physical mechanical parameters on the bed-load transport. A number of the model test problems are solved for the verification of the proposed mathematical model. The comparison of the calculation results with the established experimental data and the results of other authors is made. It was shown, that the obtained results have a good agreement with the experimental data in spite of the relative simplicity of the proposed mathematical model.
Indexed in Scopus
Full-text version of the journal is also available on the web site of the scientific electronic library eLIBRARY.RU
The journal is included in the Russian Science Citation Index
The journal is included in the RSCI
International Interdisciplinary Conference "Mathematics. Computing. Education"