Все выпуски

В работе представлены результаты расчетного исследования локальной структуры восходящего газожидкостного потока в вертикальной трубе. Математическая модель основана на использовании двухжидкостного эйлерова подхода с учетом обратного влияния пузырьков на осредненные характеристики и турбулентность несущей фазы. Турбулентная кинетическая энергия жидкости рассчитывается с применением двухпараметрической изотропной модели турбулентности $k - \varepsilon$, модифицированной для двухфазных сред. Для описания динамики распределения пузырьков по размерам используются уравнения сохранения количества частиц для отдельных групп пузырьков с различными диаметрами для каждой фракции с учетом процессов дробления и коалесценции. Численно исследовано влияние изменения степени дисперсности газовой фазы, объемного расходного газосодержания, скорости дисперсной фазы на локальную структуру и поверхностное трение в двухфазном потоке. Сравнение результатов моделирования с экспериментальными данными показало, что разработанный подход позволяет адекватно описывать турбулентные газожидкостные течения в широком диапазоне изменения газосодержания и начальных размеров пузырьков.

The results of numerical simulation of the structure of a two-phase flow of a gas-liquid bubble mixture in a vertical ascending flow in a pipe are presented. The mathematical model is based on the use of the two-fluid Eulerian approach taking into account the inverse influence of bubbles on averaged characteristics and turbulence of the carrying phase. The turbulent kinetic energy of a liquid is calculated using equations for the transfer of Reynolds stresses. To describe the dynamics of bubble size distribution, the equations of particle number conservation for individual groups of bubbles with different constant diameters for each fraction are used taking into account the processes of breakup and coalescence. The influence of changes in the degree of dispersion of the gas phase, volume flow gas content and the velocity of the dispersed phase on the local structure and surface friction in the two-phase flow is numerically investigated. Comparison of simulation results with experimental data has shown that the developed approach allows an adequate description of turbulent gas-liquid flows in a wide range of changes in gas content and initial bubble sizes.

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

In the framework of the methods and equations of the mechanics of multiphase media, a mathematical model is constructed for the formation of gas hydrate during the injection of methane into a reservoir of finite length saturated with methane and ice. The problem under study is reduced to the problem of finding two moving boundaries of phase transitions. Based on the method of catching the front in the node of the spatial grid, numerical solutions of the problem are obtained. The distributions along the spatial coordinate of temperature, pressure, and hydrate saturation are found, and the time evolution of the coordinates of the phase transition boundaries is given. Analysis of the results of computational experiments has shown that the formation of methane gas hydrate can occur both at the frontal boundary and in the extended zone. It has been established that part of the gas hydrate formed in the extended region can be further decomposed into gas and water. In this case, the extended region of hydrate formation will degenerate into the frontal surface.