Fluid DynamicsFluid DynamicsИзвестия Российской академии наук. Механика жидкости и газа1024-70843034-5340The Russian Academy of Sciences69196910.31857/S1024708425030094ArticlesСтатьиResearch ArticleSimulation of One-Sided Convection in a Porous Medium Using a Nonlinear Equation of StateМоделирование односторонней конвекции в пористой среде с использованием нелинейного уравнения состоянияSobolevaE. BСоболеваЕ. Бsoboleva@ipmnet.ruIshlinsky Institute for Problems in Mechanics of the Russian Academy of SciencesИнститут проблем механики им. А.Ю. Ишлинского РАН150720253NO3 (2025)№3 (2025)9410604102025Copyright ©; 2025, Russian Academy of SciencesCopyright ©; 2025, Российская академия наук2025Russian Academy of SciencesРоссийская академия наукhttps://rjsvd.com/1024-7084/article/view/691969

One-sided density-driven convection in a porous medium is simulated numerically with reference to hydrodynamic processes occurring during injection of carbon dioxide into underground porous formations. When carbon dioxide dissolves in water or oil, the density of solution increases. This leads to the growth of instability. A hydrodynamic model that includes the continuity equation, the equation of motion (in the form of Darcy equation), and the convection-diffusion equation has been used. The equation of state that relates the density of the fluid phase to the concentration of carbon dioxide is nonlinear. The density of solution reaches a maximum at a certain concentration, which varies. A new computational code based on the finite-difference method has been developed to solve the problem. The effect of the concentration that gives the maximum density on the parameters of convective motion and mass transfer is investigated. In particular, it is found that if the maximum density occurs at a higher concentration, the amount of carbon dioxide that is transported downward by the convective flow increases. This means that, in this case, convective dissolution is more effective in trapping of carbon dioxide at depth.

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

carbon dioxideporous mediumnonlinear equation of stateone-sided convectionconvective dissolutionnumerical simulationуглекислый газпористая среданелинейное уравнение состоянияодносторонняя конвекцияконвективное растворениечисленное моделированиеРабота выполнена при финансовой поддержке Российского Научного Фонда (код проекта 24-11-00222).Huppert H.E., Neufeld J.A. The fluid mechanics of carbon dioxide sequestration // Ann. Rev. Fluid Mech. 2014. V. 46. P. 255–272.Ennis-King J., Paterson L. Role of convective mixing in the long-term storage of carbon dioxide in deep saline formations // SPE Journal. 2005. V. 10. No. 3. P. 349–356.Paoli M.De. Convective mixing in porous media: a review of Darcy, pore-scale and Hele-Shaw studies // Eur. Phys. J. E. 2023. V. 46. P. 129.Soboleva E. Instability problems and density-driven convection in saturated porous media linking to hydrogeology: A review // Fluids. 2023. V. 8. No. 2. P. 36.Riaz A., Hesse M., Tchelepi H.A., Orr F.M. Onset of convection in a gravitationally unstable diffusive boundary layer in porous media // J. Fluid Mech. 2006. V. 548. P. 87–111.Pau G.S.H., Bell J.B., Pruess K., Almgren A.S., Lijewski M.J., Zhang K. High-resolution simulation and characterization of density-driven flow in CO2 storage in saline aquifers // Adv. Water Res. 2010. V. 33. P. 443–455.Soboleva E.B. Density-driven convection in an inhomogeneous geothermal reservoir // Int. J. Heat Mass Transf. 2018. V. 127 (C). P. 784–798.Neufeld J.A., Hesse M.A., Riaz A., Hallworth M.A., Tchelepi H.A., Huppert H.E. Convective dissolution of carbon dioxide in saline aquifers // Geophys. Res. Lett. 2010. V. 37. P. 22404.Daniel D., Riaz A. Effect of viscosity contrast on gravitationally unstable diffusive layers in porous media // Phys. Fluids. 2014. V. 26. P. 116601.Elenius M.T., Gasda S.E. Convective mixing driven by non-monotonic density // Transp. Porous Med. 2021. V. 138. P. 133–155.Leonard B.P. A stable and accurate convective modeling procedure based on quadratic upstream interpolation // Comput. Methods Appl. Mech. Eng. 1979. V. 19. P. 59–98.Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоиздат, 1984. 124 c.Nield D.A., Bejan A. Convection in Porous Media. New York: Springer, 2006. 640 p.Соболева Е.Б., Цыпкин Г.Г. Численное моделирование конвективных течений в грунте при испарении воды, содержащей растворенную примесь // Изв. РАН. МЖГ. 2014. №5. С. 88–99.Соболева Е.Б. Метод численного моделирования концентрационно-конвективных течений в пористых средах в приложении к задачам геологии // Ж. вычисл. матем. и матем. физ. 2019. Т. 59. №11. С. 162–173.Соболева Е.Б. Начало конвекции Рэлея–Тейлора в пористой среде // Изв. РАН. МЖГ. 2021. №2. С. 52–62.Соболева Е.Б. Влияние конечных возмущений плотности на развитие неустойчивости Рэлея–Тейлора в пористой среде // Теор. и матем. физ. 2022. Т. 211. №2. С. 333–346.Соболева Е.Б. Численное моделирование фильтрационных концентрационно-конвективных течений с контрастом вязкости // Ж. вычисл. матем. и матем. физ. 2022. Т. 62. №11. С. 1927–1939.Wang Zh., Hou J. Measurement of CO2 diffusion coefficients in both bulk liquids and carven filling porous media of fractured-vuggy carbonate reservoirs at 50 MPa and 393 K // RSC Adv. 2021. V. 11. P. 19712.Коллинз Р. Течения жидкостей через пористые материалы. М.: Мир, 1964. 352 с.